2024年5月19日发(作者：袭阑)

①

：

% 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;

%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');

hold on;

%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);

2024年5月19日发(作者：袭阑)

①

：

% 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;

%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');

hold on;

%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);