post_handle_convergence.m

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

function post_handle_convergence(ll)
include_globals
include_flags

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


TE_reflection_eff_convergence(:,ll,1) = TE_reflection(:,2);
TE_reflection_eff_convergence(:,ll,2) = TE_reflection(:,3);


TE_transmission_eff_convergence(:,ll,1) = TE_transmission(:,2);
TE_transmission_eff_convergence(:,ll,2) = TE_transmission(:,3);

TM_reflection_eff_convergence(:,ll,1) = TM_reflection(:,2);
TM_reflection_eff_convergence(:,ll,2) = TM_reflection(:,3);

TM_transmission_eff_convergence(:,ll,1) = TM_transmission(:,2);
TM_transmission_eff_convergence(:,ll,2) = TM_transmission(:,3);
end