% 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');
hold on;
%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');
hold on;
2023-03-16 14:59:28
53KB
薄膜干涉
1