calc_the_reflection.m

FDTD TMz_with_different_PMLs

load cos_gaussian_12GHz_Bandwidth24GHz_ref18.mat
load cos_gaussian_12GHz_Bandwidth24GHz_spml18.mat
load cos_gaussian_12GHz_Bandwidth24GHz_upml18.mat
load cos_gaussian_12GHz_Bandwidth24GHz_cpml18.mat
load cos_gaussian_12GHz_Bandwidth24GHz_npml18.mat

reflection_t_A_spml = abs(uxdirect_spml - uxdirect_ref)/max(uxdirect_ref);
reflection_t_A_upml = abs(uxdirect_upml - uxdirect_ref)/max(uxdirect_ref);
reflection_t_A_cpml = abs(uxdirect_cpml - uxdirect_ref)/max(uxdirect_ref);
reflection_t_A_npml = abs(uxdirect_npml - uxdirect_ref)/max(uxdirect_ref);

reflection_t_B_spml = abs(uxycorner_spml - uxycorner_ref)/max(uxycorner_ref);
reflection_t_B_upml = abs(uxycorner_upml - uxycorner_ref)/max(uxycorner_ref);
reflection_t_B_cpml = abs(uxycorner_cpml - uxycorner_ref)/max(uxycorner_ref);
reflection_t_B_npml = abs(uxycorner_npml - uxycorner_ref)/max(uxycorner_ref);

reflection_tscale=(1e9)*stime_ref;

semilogy(reflection_tscale,reflection_t_A_spml,'r');
hold on;
semilogy(reflection_tscale,reflection_t_A_upml,'b.');
hold on;
semilogy(reflection_tscale,reflection_t_A_cpml,'g');
hold on;
semilogy(reflection_tscale,reflection_t_A_npml,'magenta.');
xlim([0 2]);
ylim([1e-8 1e-4]);

% semilogy(reflection_tscale,reflection_t_B_spml,'r');
% hold on;
% semilogy(reflection_tscale,reflection_t_B_upml,'b');
% hold on;
% semilogy(reflection_tscale,reflection_t_B_cpml,'g');
% hold on;
% semilogy(reflection_tscale,reflection_t_B_npml,'magenta');
% xlim([0 2]);
% ylim([1e-8 1e-4]);