main_corrugation_h.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_corrugation_h
clc;clear all;close all
format long

include_globals

lambda0 = 0.6328; % micrometer
%   Parameters of Grating
P = 0;
thick_e = [0];
bottom_y = [0];  
epsilon_e = [(0.135+1i*3.98)^2,(1.0)^2];
miu_u = [1.0,1.0];
period_d = 0.5;
period_k = 2*pi/period_d;
N = 15;     %   N: Truncation Parameter
DIM = 2*N+1;
%******************************************************************
hh = [0.01:0.005:0.05];
angang = [11:0.2:15]*pi/180;
%******************************************************************
NUM1 = length(hh);
NUM2 = length(angang);
for ll = 1:NUM1   
    corrugation_h = hh(ll);
    for jj = 1:NUM2
        theta = angang(jj);
        para_setting
        rayleigh_fourier_computation    
        post_handle_corrugation(ll,jj)
    end
end

savefile = 'corrugation_data.mat';
save(savefile, 'hh','angang','TM_reflection_eff_corrugation','TM_transmission_eff_corrugation','TE_reflection_eff_corrugation','TE_transmission_eff_corrugation')
% Modify and run "corrugation_analysis.mat" to analyze the result
end