make_cookie.m

利用电磁场的源激发方法来计算光子晶体波导例如光子晶体光纤

clear
oGd.nSources = 75;
oGd.nTestingPoints = round(oGd.nSources*2);
oGd.alphaIn = 1.5;
oGd.alphaOut = 1.5;
R = 2000e-9;
mod_amp = 0.3;
n = 4;
curveArray = primitive('cookie', R, mod_amp, n);

oGd.curveArray = curveArray;
oGd.sdArray = compute_sd_array(curveArray);
oGd.nSd = length(oGd.sdArray);
% set material parameters
oGd.sdArray(1).Er = 1;
oGd.sdArray(2).Er = 'sellmeier';
oGd.offsetX = 0;
oGd.offsetY = 0;
save cookie.gd oGd