(3)反射系数r，反射率R与入射角x和折射率之比n的关系

clear;

clc;

clf;

[n,zeta1]=meshgrid(0.6:0.05:1.5,0:pi/50:pi/2);

zeta2=real(asin(n.*sin(zeta1))); %折射角

x=zeta1*180/pi; %转换成角度

rpz=-n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分子

rpm=n.*cos(zeta2)+cos(zeta1); %平行分量反射部分分母

rp=rpz./rpm; %平行分量反射系数

rsz=n.*cos(zeta1)-cos(zeta2); %垂直分量反射部分分子

rsm=n.*cos(zeta1)+cos(zeta2); %垂直分量反射部分分母

rs=rsz./rsm; %垂直分量反射系数

Rp=rp.^2; %平行分量反射率

Rs=rs.^2; %垂直分量反射率

%brewster角所在的位置

Brewster=acot(n)*180/pi; %求布鲁斯特角

zeta3=real(asin(n.*sin(zeta1))); %以布鲁斯特角入射的折射角rpz0=-n.*cos(zeta3)+cos(Brewster);

rpm0=n.*cos(zeta3)+cos(Brewster);

rp0=rpz0./rpm0;

rsz0=n.*cos(Brewster)-cos(zeta3);

rsm0=n.*cos(Brewster)+cos(zeta3);

rs0=rsz0./rsm0;

Rp0=rp0.^2;

Rs0=rs0.^2;