[sqqssqq]

下面是程序,急啊~~~~~
T=20:20:200;
a1=5.35583;
a2=0.100473;
a3=0.20692;
a4=100;
a5=11.34929;
a6=1.5334e-2;
b1=4.629e-7;
b2=3.862e-8;
b3=-0.89e-8;
b4=2.657e-5;
p=0.532;%波长532nm
n=30;
f=(T-24.5).*(T+570.832);
mp=(a1+b1*f+(a2+b2*f)./(p^2-(a3+b3*f).^2)+(a4+b4*f)./(p^2-a5^2)-a6*p^2).^(1/2);
mi=(a1+b1*f+(a2+b2*f)./(i.^2-(a3+b3*f).^2)+(a4+b4*f)./(i.^2-a5^2)-a6*i.^2).^(1/2);
s=(1/p-1./i).^(-1);
ms=(a1+b1*f+(a2+b2*f)./(s.^2-(a3+b3*f).^2)+(a4+b4*f)./(s.^2-a5^2)-a6*s.^2).^(1/2);
i=mi.*(1/n+ms./s-mp/p);
plot(T,s)
xlabel('T(℃)')
ylabel('λs(μm)')