funn.m

来自「遗传算法与单纯形算法复合用于光纤与波导芯片对准仿真技术分析」· M 代码 · 共 16 行

function eta = funn(x0)

% tx与ty以度表示
% lambda表示波长，wx0表示波导模场长半轴，wy0表示波导模场短半轴，wf0表示光纤模场半径
lambda = 1.31;
wx0 = 4.5;
wy0 = 4.5;
wf0 = 4.6;
x = x0(1); y = x0(2); tx = x0(3)*pi/180; ty = x0(4)*pi/180; z = x0(5);
wx_z = wx0*sqrt(1+(lambda*z/(pi*wx0^2))^2);
wy_z = wy0*sqrt(1+(lambda*z/(pi*wy0^2))^2);
kx = 4*wx0^2*wf0^2/((wx0^2+wf0^2)^2+lambda^2*z^2/pi^2);
ky = 4*wy0^2*wf0^2/((wy0^2+wf0^2)^2+lambda^2*z^2/pi^2);
eta_x = kx*exp(-kx*(x^2*(1/wx0^2+1/wf0^2)/2+pi^2*tx^2*(wx_z^2+wf0^2)/(2*lambda^2)-x*tx*z/wx0^2));
eta_y = ky*exp(-ky*(y^2*(1/wy0^2+1/wf0^2)/2+pi^2*ty^2*(wy_z^2+wf0^2)/(2*lambda^2)-y*ty*z/wy0^2));
eta = eta_x*eta_y;