convergence_analysis.m

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

function convergence_analysis
clc;clear all;close all
data = load('convergence_data');
temp_data =data.TM_reflection_eff_convergence; % choose the wave you want to analyze
NN = data.NN;
[mm,nn,pp] = size(temp_data);

figure
hold on
jj = 2; % choose the row to analyze
yy = temp_data(jj,:,1);
plot(NN,yy,'-*')
xlabel('Truncation parameter N')
ylabel('Diffraction efficiency')
title('Convergence behavior of S-Matrix algorithm')
hold off
end