薄膜干涉的仿真程序.docx
- 文档编号:10963201
- 上传时间:2023-02-24
- 格式:DOCX
- 页数:14
- 大小:15.03KB
薄膜干涉的仿真程序.docx
《薄膜干涉的仿真程序.docx》由会员分享,可在线阅读,更多相关《薄膜干涉的仿真程序.docx(14页珍藏版)》请在冰豆网上搜索。
薄膜干涉的仿真程序
1:
%n1=1.2
n0=1;
n2=1.5;
n1=1.2;
langmuda=400*10^(-9);
h=0:
0.0001*langmuda/1.2:
langmuda/1.2;
x=n1*h;
j0=0;
j1=asin(n0*sin(j0)/n1);
f=4*pi*n1*h*cos(j1)/langmuda;
a=n0-n2;
b=n0+n2;
c=((n0.*n2)./n1)-n1;
d=((n0.*n2)./n1)+n1;
A=a^2*(cos(f/2).^2)+c^2*(sin(f/2).^2);
B=b^2*(cos(f/2).^2)+d^2*(sin(f/2).^2);
R=A./B;
plot(x,R,'-r');
holdon;
%n1=1
n0=1;
n2=1.5;
n1=1;
langmuda=400*10^(-9);
h=0:
0.0001*langmuda:
langmuda;
x=n1*h;
j0=0;
j1=asin(n0*sin(j0)/n1);
f=4*pi*n1*h*cos(j1)./langmuda;
a=n0-n2;
b=n0+n2;
c=((n0.*n2)./n1)-n1;
d=((n0.*n2)./n1)+n1;
A=a.^2.*(cos(f./2).^2)+c.^2.*(sin(f./2).^2);
B=b.^2.*(cos(f./2).^2)+d.^2.*(sin(f./2).^2);
R=A./B;
plot(x,R,'-k');
holdon;
%n1=1.5
n0=1;
n2=1.5;
n1=1.5;
langmuda=400*10^(-9);
h=0:
0.0001*langmuda/1.5:
langmuda/1.5;
x=n1*h;
j0=0;
j1=asin(n0*sin(j0)/n1);
f=4*pi*n1*h*cos(j1)/langmuda;
a=n0-n2;
b=n0+n2;
c=((n0.*n2)./n1)-n1;
d=((n0.*n2)./n1)+n1;
A=a^2*(cos(f/2).^2)+c^2*(sin(f/2).^2);
B=b^2*(cos(f/2).^2)+d^2*(sin(f/2).^2);
R=A./B;
plot(x,R,'-.g');
holdon;
%n=1.4
n0=1;
n2=1.5;
n1=1.4;
langmuda=400*10^(-9);
h=0:
0.0001*langmuda/1.4:
langmuda/1.4;
x=n1*h;
j0=0;
j1=asin(n0*sin(j0)/n1);
f=4*pi*n1*h*cos(j1)/langmuda;
a=n0-n2;
b=n0+n2;
c=((n0.*n2)./n1)-n1;
d=((n0.*n2)./n1)+n1;
A=a^2*(cos(f/2).^2)+c^2*(sin(f/2).^2);
B=b^2*(cos(f/2).^2)+d^2*(sin(f/2).^2);
R=A./B;
plot(x,R,'-b');
holdon;
%n=1.7
n0=1;
n1=1.5;
n1=1.7;
langmuda=400*10^(-9);
h=0:
0.0001*langmuda/1.7:
langmuda/1.7;
x=n1*h;
j0=0;
j1=asin(n0*sin(j0)/n1);
f=4*pi*n1*h*cos(j1)/langmuda;
a=n0-n2;
b=n0+n2;
c=((n0.*n2)./n1)-n1;
d=((n0.*n2)./n1)+n1;
A=a^2*(cos(f/2).^2)+c^2*(sin(f/2).^2);
B=b^2*(cos(f/2).^2)+d^2*(sin(f/2).^2);
R=A./B;
plot(x,R,'-.b');
holdon;
%n=2.0
n0=1;
n2=1.5;
n1=2;
langmuda=400*10^(-9);
h=0:
0.0001*langmuda/2:
langmuda/2;
x=n1*h;
j0=0;
j1=asin(n0*sin(j0)/n1);
f=4*pi*n1*h*cos(j1)/langmuda;
a=n0-n2;
b=n0+n2;
c=((n0.*n2)./n1)-n1;
d=((n0.*n2)./n1)+n1;
A=a^2*(cos(f/2).^2)+c^2*(sin(f/2).^2);
B=b^2*(cos(f/2).^2)+d^2*(sin(f/2).^2);
R=A./B;
plot(x,R,'-.g');
holdon;
%n=3.0
n0=1;
n2=1.5;
n1=3;
langmuda=400*10^(-9);
h=0:
0.0001*langmuda/3:
langmuda/3;
x=n1*h;
j0=0;
j1=asin(n0*sin(j0)/n1);
f=4*pi*n1*h*cos(j1)/langmuda;
a=n0-n2;
b=n0+n2;
c=((n0*n2)./n1)-n1;
d=((n0*n2)./n1)+n1;
A=a^2*(cos(f/2).^2)+c^2*(sin(f/2).^2);
B=b^2*(cos(f/2).^2)+d^2*(sin(f/2).^2);
R=A./B;
plot(x,R,'-k');
holdon;
legend('n1=1.2','n1=1','n1=1.5','n1=1.4','n1=1.7','n1=2','n1=3');
2
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
%i0=0:
pi/10000:
pi/2;%入射角为i0
%w=4*exp(-7);%给定的为可见光的一个波长400nm
w=370:
1/10:
1000;
R=0:
1/9:
8;
lamd=550;
q=lamd/4;
w0=q;%w0为光学厚度
i0=0;
n0=1;
n1=1.38;
i1=asin((n0*sin(i0))/n1);%折射角为i1
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
%f是相邻两个出射光束间的相位差
r1s=-(n1*cos(i1)-n0*cos(i0))/(n1*cos(i1)+n0*cos(i0));
r1p=(n1/cos(i1)-n0/cos(i0))/(n1/cos(i1)+n0/cos(i0));
R1s=r1s.^2;
R1p=r1p.^2;
R1R=(R1s+R1p)/2;
r1=sqrt(R1R);
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
i2=asin((n1.*sin(i1))./n2);%折射角为i1
r2s=(n1*cos(i1)-n2*cos(i2))/(n1*cos(i1)+n2*cos(i2));
r2p=(n2/cos(i1)-n1/cos(i2))/(n2/cos(i1)+n1/cos(i2));
R2s=r2s.^2;
R2p=r2p.^2;
R2R=(R2s+R2p)/2;
r2=sqrt(R2R);
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
A=r1.^2+r2.^2;
B=1+(r1.^2)*(r2.^2);
C=2*r1*r2.*cos(f);
R1=(A+C)./(B+C);
R=R1/100;
plot(w,R,'r')
holdon
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
%i0=0:
pi/10000:
pi/2;%入射角为i0
%w=4*exp(-7);%给定的为可见光的一个波长400nm
w=370:
1/10:
1000;
R=0:
1:
8;
lamd=550;
q=lamd/4;
w0=q;%w0为光学厚度
i0=pi/6;
n0=1;
n1=1.38;
i1=asin((n0*sin(i0))/n1);%折射角为i1
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
%f是相邻两个出射光束间的相位差
r1s=-(n1*cos(i1)-n0*cos(i0))/(n1*cos(i1)+n0*cos(i0));
r1p=(n1/cos(i1)-n0/cos(i0))/(n1/cos(i1)+n0/cos(i0));
R1s=r1s.^2;
R1p=r1p.^2;
R1R=(R1s+R1p)/2;
r1=sqrt(R1R);
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
i2=asin((n1.*sin(i1))./n2);%折射角为i1
r2s=(n1*cos(i1)-n2*cos(i2))/(n1*cos(i1)+n2*cos(i2));
r2p=(n2/cos(i1)-n1/cos(i2))/(n2/cos(i1)+n1/cos(i2));
R2s=r2s.^2;
R2p=r2p.^2;
R2R=(R2s+R2p)/2;
r2=sqrt(R2R);
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
A=r1.^2+r2.^2;
B=1+(r1.^2)*(r2.^2);
C=2*r1*r2.*cos(f);
R1=(A+C)./(B+C);
R=R1/100;
plot(w,R,'g')
holdon
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
%i0=0:
pi/10000:
pi/2;%入射角为i0
%w=4*exp(-7);%给定的为可见光的一个波长400nm
w=370:
1/10:
1000;
R=0:
1:
8;
lamd=550;
q=lamd/4;
w0=q;%w0为光学厚度
i0=40*pi/180;
n0=1;
n1=1.38;
i1=asin((n0*sin(i0))/n1);%折射角为i1
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
%f是相邻两个出射光束间的相位差
r1s=-(n1*cos(i1)-n0*cos(i0))/(n1*cos(i1)+n0*cos(i0));
r1p=(n1/cos(i1)-n0/cos(i0))/(n1/cos(i1)+n0/cos(i0));
R1s=r1s.^2;
R1p=r1p.^2;
R1R=(R1s+R1p)/2;
r1=sqrt(R1R);
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
i2=asin((n1.*sin(i1))./n2);%折射角为i1
r2s=(n1*cos(i1)-n2*cos(i2))/(n1*cos(i1)+n2*cos(i2));
r2p=(n2/cos(i1)-n1/cos(i2))/(n2/cos(i1)+n1/cos(i2));
R2s=r2s.^2;
R2p=r2p.^2;
R2R=(R2s+R2p)/2;
r2=sqrt(R2R);
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
A=r1.^2+r2.^2;
B=1+(r1.^2)*(r2.^2);
C=2*r1*r2.*cos(f);
R1=(A+C)./(B+C);
R=R1/100;
plot(w,R,'b')
holdon
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
%i0=0:
pi/10000:
pi/2;%入射角为i0
%w=4*exp(-7);%给定的为可见光的一个波长400nm
w=370:
1/10:
1000;
R=0:
1:
8;
lamd=550;
q=lamd/4;
w0=q;%w0为光学厚度
i0=50*pi/180;
n0=1;
n1=1.38;
i1=asin((n0*sin(i0))/n1);%折射角为i1
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
%f是相邻两个出射光束间的相位差
r1s=-(n1*cos(i1)-n0*cos(i0))/(n1*cos(i1)+n0*cos(i0));
r1p=(n1/cos(i1)-n0/cos(i0))/(n1/cos(i1)+n0/cos(i0));
R1s=r1s.^2;
R1p=r1p.^2;
R1R=(R1s+R1p)/2;
r1=sqrt(R1R);
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
i2=asin((n1.*sin(i1))./n2);%折射角为i1
r2s=(n1*cos(i1)-n2*cos(i2))/(n1*cos(i1)+n2*cos(i2));
r2p=(n2/cos(i1)-n1/cos(i2))/(n2/cos(i1)+n1/cos(i2));
R2s=r2s.^2;
R2p=r2p.^2;
R2R=(R2s+R2p)/2;
r2=sqrt(R2R);
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
A=r1.^2+r2.^2;
B=1+(r1.^2)*(r2.^2);
C=2*r1*r2.*cos(f);
R1=(A+C)./(B+C);
R=R1/100;
plot(w,R,'k')
holdon
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
%i0=0:
pi/10000:
pi/2;%入射角为i0
%w=4*exp(-7);%给定的为可见光的一个波长400nm
w=300:
1/10:
1000;
R=0:
1:
8;
lamd=550;
q=lamd/4;
w0=q;%w0为光学厚度
i0=pi/3;
n0=1;
n1=1.38;
i1=asin((n0*sin(i0))/n1);%折射角为i1
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
%f是相邻两个出射光束间的相位差
r1s=-(n1*cos(i1)-n0*cos(i0))/(n1*cos(i1)+n0*cos(i0));
r1p=(n1/cos(i1)-n0/cos(i0))/(n1/cos(i1)+n0/cos(i0));
R1s=r1s.^2;
R1p=r1p.^2;
R1R=(R1s+R1p)/2;
r1=sqrt(R1R);
n0=1;%已知n0和n2,得到R随着f即R随n1h的变化规律
n2=1.5;
n1=1.38;
i2=asin((n1.*sin(i1))./n2);%折射角为i1
r2s=(n1*cos(i1)-n2*cos(i2))/(n1*cos(i1)+n2*cos(i2));
r2p=(n2/cos(i1)-n1/cos(i2))/(n2/cos(i1)+n1/cos(i2));
R2s=r2s.^2;
R2p=r2p.^2;
R2R=(R2s+R2p)/2;
r2=sqrt(R2R);
f=(4*pi*w0*cos(i1))./w;%w为wavelength,即波长的简称
A=r1.^2+r2.^2;
B=1+(r1.^2)*(r2.^2);
C=2*r1*r2.*cos(f);
R1=(A+C)./(B+C);
R=R1/100;
plot(w,R,'m')
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 薄膜 干涉 仿真 程序
![提示](https://static.bdocx.com/images/bang_tan.gif)