post_handle.m

来自「用S-Matrix算法计算光栅衍射问题。」· M 代码 · 共 23 行

function post_handle
include_globals
include_flags

TE_ind_reflection = find(real(TE_eta_reflection) ~= 0);
TE_ind_transmission = find(real(TE_eta_transmission) ~= 0);
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 = find(real(TM_eta_reflection) ~=0);
TM_ind_transmission = find(real(TM_eta_transmission) ~=0);
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]';
end