flat_metallic_mirror_reflectivity.m

用S-Matrix算法计算光栅衍射问题。

function flat_metallic_mirror_reflectivity
clc;clear all;close all

ni = 1.0;
nt = 0.135 + 1i*3.98;
theta_inc = [11:0.5:15]*pi/180;
sin_theta_inc = sin(theta_inc);
sin_theta_tran = ni/nt*sin_theta_inc;
cos_theta_inc = cos(theta_inc);
cos_theta_tran = sqrt(1 - sin_theta_tran.^2);

RTE = abs((ni*cos_theta_inc-nt*cos_theta_tran)./(ni*cos_theta_inc + nt*cos_theta_tran)).^2;
RTM = abs((nt*cos_theta_inc-ni*cos_theta_tran)./(ni*cos_theta_inc + nt*cos_theta_tran)).^2;

xx = theta_inc*180/pi;
plot(xx,RTE,'-*r',xx,RTM,'-ob')
axis([min(xx),max(xx),0.9,1.0])
xlabel('\theta/degree')
ylabel('Reflectivity')
legend('TE','TM')
end