减反射膜matlab设计程序

%双层减反射膜程序,可用于太阳电池件反射膜的设计

%作者：ZHOU-JIANPENG

n1=1;%空气折射率

n2=1.38;%第一层膜折射率

n3=2.3;%第二层膜折射率

n4=3.58;%衬底折射率

a=100;%第一层膜厚度

b=60;%第二层膜厚度

c1=0;%入射角度

lambda1=300;%起始波长

lambda2=1200;%截止波长

f1='C:\';%数据存储位置

f2=datestr(now,30);%计算现有时间

f3=strcat(f1,f2);

f4='反射率数据.txt';

f5=strcat(f3,f4);

fid=fopen(f5,'a+');

fprintf(fid,'%s %s\n','LAMDA','R');%在文本第一行输入列标题'LAMDA''R' for d=lambda1:lambda2;

%计算各层折射角度

c2=asin(n1*sin(c1)/n2);

c3=asin(n2*sin(c2)/n3);

c4=asin(n3*sin(c3)/n4);

%相位因子

delta2=2*pi*n2*a*cos(c2)/d;

delta3=2*pi*n3*b*cos(c3)/d;

f=4*pi*1e-7;

e=1e-9/(36*pi);

m=sqrt(e/f);

z1=m*cos(c1)*n1;

z2=m*cos(c2)*n2;

z3=m*cos(c3)*n3;

z4=m*cos(c4)*n4;

m1=cos(delta2);

m2=-i*sin(delta2)/z2;

m3=-i*z2*sin(delta2);

m4=cos(delta2);

M=[m1 m2;m3 m4];

u1=cos(delta3);

u2=-i*sin(delta3)/z3;

u3=-i*z3*sin(delta3);

u4=cos(delta3);

U=[u1 u2;u3 u4];

W=M*U;

W11=W(1,1);

W12=W(1,2);

W13=W(2,1);

W14=W(2,2);

r=(z1*(W11+z4*W12)-W13-z4*W14)/(z1*(W11+z4*W12)+W13+z4*W14); r1=r*r';

s=d-(lambda1-1);

k(1,s)=r1*100;

fprintf(fid,'%f %f\n',d,r1);

end

d=lambda1:lambda2;

plot(d,k,'')%绘制反射率曲线

axis tight;

axis square;

xlabel('波长/nm');

ylabel('反射率/%');

title('减反射膜反射率图谱')

legend('反射率',1);

success='success'%计算完成显示success

运行结果如下:数据以txt格式保存在C盘