post_handle_corrugation.m

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

function post_handle_corrugation(ll,jj)
include_globals
include_flags

TE_ind_reflection = N+1;
TE_ind_transmission = N+1;
TE_reflection_angle = abs(acos(beta_b(P+1+1,TE_ind_reflection)/propagation_k(P+1+1))*180/pi);
TE_transmission_angle = abs(acos(beta_b(0+1,TE_ind_transmission)/propagation_k(0+1))*180/pi);
TE_ind_R = TE_ind_reflection - (N+1);
TE_ind_T = TE_ind_transmission - (N+1);
TE_reflection = [TE_ind_R ;TE_eta_reflection(TE_ind_reflection);sign(alpha_a(TE_ind_reflection)).*TE_reflection_angle]';
TE_transmission = [TE_ind_T;TE_eta_transmission(TE_ind_transmission);sign(alpha_a(TE_ind_transmission)).*TE_transmission_angle]';


TM_ind_reflection = N+1;
TM_ind_transmission = N+1;
TM_reflection_angle = abs(acos(beta_b(P+1+1,TM_ind_reflection)/propagation_k(P+1+1))*180/pi);
TM_transmission_angle = abs(acos(beta_b(0+1,TM_ind_transmission)/propagation_k(0+1))*180/pi);
TM_ind_R = TM_ind_reflection - (N+1);
TM_ind_T = TM_ind_transmission - (N+1);
TM_reflection = [TM_ind_R ;TM_eta_reflection(TM_ind_reflection);sign(alpha_a(TM_ind_reflection)).*TM_reflection_angle]';
TM_transmission = [TM_ind_T;TM_eta_transmission(TM_ind_transmission);sign(alpha_a(TM_ind_transmission)).*TM_transmission_angle]';


TE_reflection_eff_corrugation(ll,jj) = TE_eta_reflection(TE_ind_reflection);
TE_transmission_eff_corrugation(ll,jj) = TE_eta_transmission(TE_ind_transmission);

TM_reflection_eff_corrugation(ll,jj) = TM_eta_reflection(TM_ind_reflection);
TM_transmission_eff_corrugation(ll,jj) = TM_eta_transmission(TM_ind_transmission);
end