main_basic.m

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

% Content: Code for diffraction grating computation
% Method: Rayleigh-Fourier Method
% Algorithm: S-Matrix
% By: Barak DEP
% Date: 28/04/2009
function main_basic
clc;clear all;close all
include_globals
format long
%   Parameters of Incident Wave
%lambda0 = 1.0;
%theta = 30*pi/180; % incident angle, [rad]
lambda0 = 0.6328;
theta = 12*pi/180;
%   Parameters of Grating
P = 1;
%thick_e = [0,.25*lambda0]; % take care: e0 = y0 - y(-1) = y0 = 0
corrugation_h = [0.0,0.0];
thick_e = [0,0.0-corrugation_h(2)/2]; % take care: e0 = y0 - y(-1) = y0 = 0
bottom_y = [0, thick_e(1)];  
%period_d = 1.3*lambda0;
period_d = 0.37;
period_k = 2*pi/period_d;
%epsilon_e = [(1.45)^2,(1.6)^2,(1.0)^2];
epsilon_e = [(1.0)^2,(1.5)^2,(1.0)^2];
miu_u = [1.0,1.0,1.0];
%corrugation_h = [0.1*period_d,0.1*period_d];

%   N: Truncation Parameter
N = 5;     
DIM = 2*N+1;
para_setting
rayleigh_fourier_computation    
post_handle

%   Display the results
TE_reflection
TE_transmission
TM_reflection
TM_transmission
end