corrugation_analysis.m

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

function corrugation_analysis
clc;clear all;close all

data = load('corrugation_data');
hh = data.hh;
angang = data.angang*180/pi;
temp_data1 =data.TE_reflection_eff_corrugation;
temp_data2 =data.TM_reflection_eff_corrugation;

[XX,YY] = meshgrid(hh,angang);
figure
hold on
mesh(XX,YY,temp_data1')
title(' 0^t^h Reflection Order of TE Wave ')
xlabel('h/\mum')
ylabel('\theta/degree')
zlabel('\eta')
hold off

figure
hold on
mesh(XX,YY,temp_data2')
title(' 0^t^h Reflection Order of TM Wave ')
xlabel('h/\mum')
ylabel('\theta/degree')
zlabel('\eta')
hold off

% search for minimum of TM efficiency
min_efficiency = min(min(temp_data2));
[ll,jj] = find(temp_data2==min_efficiency);
min_angle = angang(jj);
min_hh = hh(ll);
min_hh
min_angle
min_efficiency
end