rayleigh_fourier_computation.m

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

function rayleigh_fourier_computation
include_globals
include_flags

Tuu = eye(DIM);
Rud = eye(DIM);
Rdu = eye(DIM);
Tdd = eye(DIM);

S_TE = get_S_cell(P,TE);
S_TM = get_S_cell(P,TM);

%cell2mat(S_TE)

Tuu_TE = cell2mat(S_TE(1,1));
Rud_TE = cell2mat(S_TE(1,2));
Rdu_TE = cell2mat(S_TE(2,1));
Tdd_TE = cell2mat(S_TE(2,2));

Tuu_TM = cell2mat(S_TM(1,1));
Rud_TM = cell2mat(S_TM(1,2));
Rdu_TM = cell2mat(S_TM(2,1));
Tdd_TM = cell2mat(S_TM(2,2));

Inc_down = [zeros(1,N),exp(-1i*beta_b(P+1+1,0+1)*bottom_y(P+1)),zeros(1,N)]';
E_transmission = Tdd_TE*Inc_down;   E_reflection = Rud_TE*Inc_down;
H_transmission = Tdd_TM*Inc_down;   H_reflection = Rud_TM*Inc_down;


TE_eta_reflection = (beta_b(P+1+1,:)/beta_b(P+1+1,N+1)).*[abs(E_reflection).^2]';
TE_eta_transmission = (-miu_u(P+1+1)/(-miu_u(0+1)*beta_b(P+1+1,N+1)))*beta_b(0+1,:).*[abs(E_transmission).^2]';
TM_eta_reflection = (beta_b(P+1+1,:)/beta_b(P+1+1,N+1)).*[abs(H_reflection).^2]';
TM_eta_transmission = (epsilon_e(P+1+1)/(epsilon_e(0+1)*beta_b(P+1+1,N+1)))*beta_b(0+1,:).*[abs(H_transmission).^2]';
end