仿真

clc;clear
lax=1.2 : 0.02 : 2;
s1=(1+(2.8939*lax.^2)./(lax.^2-0.13967^2)).^1/2;
% plot(t,Wt,'r-o','markersize',6,'markeredgecolor','r','markerfacecolor','r')
plot(lax,s1,'r')
xlabel('wavelength(μm)')
ylabel('Refractive index')
box off
grid on

在这里插入代码片clc; 
clear all; 
close all; 
clf;

s1=0
lax=1.2 : 0.01 : 2;
b=[0.6961663,0.4079426,0.8974794];
la=[0.0684046,0.1162414,9.896161];
for i=1:3
    s1=s1+(b(i)*lax.^2)./(lax.^2-la(i).^2);
end
s1=(1+s1).^1/2;
figure(1);
plot(lax,s1,'b');
xlabel('wavelength(μm)')
ylabel('Refractive index')

clc;clear
lax=1.2 : 0.02 : 2;
la=11.686;
a=0.9398;
b=8.1046*10^-3;
lax1=0.0684046;
laxg=1.1071;
s1=(la+a./lax.^2+(b*lax1^2)./(lax.^2-laxg^2)).^1/2;
% plot(t,Wt,'r-o','markersize',6,'markeredgecolor','r','markerfacecolor','r')
plot(lax,s1,'r')
xlabel('wavelength(μm)')
ylabel('Refractive index')
box off
grid on