para_setting.m

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

function para_setting
include_globals

%   Parameters of Diffraction Orders
  k0 = 2*pi/lambda0;
  propagation_k = sqrt(epsilon_e)*k0;
  temp_alpha0 = k0*sin(theta);
  temp_m = [-N:1:N];
  alpha_a = temp_alpha0 + temp_m*period_k;
  
  beta_b = zeros((P+2),(2*N+1));
  for jj = 0:(P+1 )
      beta_b(jj+1,:) = sqrt(propagation_k(jj+1)^2 - alpha_a.^2);
  end
  
end