funn.m

遗传算法与爬山法用于光纤与波导芯片对准仿真分析

function eta = fun(R0)

% tx与ty以度表示
% lambda表示波长，wx0表示波导模场长半轴，wy0表示波导模场短半轴，wf0表示光纤模场半径
lambda = 1.31;
wx0 = 4.5;
wy0 = 4.5;
wf0 = 4.6;
x = R0(1); y = R0(2); tx = R0(3)*pi/180; ty = R0(4)*pi/180; z = R0(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;