fun.m

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

function a = fun(x,y,tx,ty,z)
% tx与ty以度表示
% lambda表示波长，wx0表示波导模场长半轴，wy0表示波导模场短半轴，wf0表示光纤模场半径
wavelength = 1.31;
w0x = 4.5;
w0y = 4.5;
wf0 = 4.6;

thetax=pi*tx/180;
thetay=pi*ty/180;

w0xz=sqrt(w0x^2*(1+(wavelength*z./(pi*w0x^2)).^2));
w0yz=sqrt(w0y^2*(1+(wavelength*z./(pi*w0y^2)).^2));

kx=4*w0x^2*wf0^2./((w0x^2+wf0^2)^2+wavelength^2*z.^2/pi^2);
ky=4*w0y^2*wf0^2./((w0y^2+wf0^2)^2+wavelength^2*z.^2/pi^2);

intax=kx.*exp(-kx.*(x.^2/2*(1/w0x^2+1/wf0^2)+pi^2*thetax.^2.*(w0xz.^2+wf0^2)./(2*wavelength^2)-x.*thetax.*z/w0x^2));
intay=ky.*exp(-ky.*(y.^2/2*(1/w0y^2+1/wf0^2)+pi^2*thetay.^2.*(w0yz.^2+wf0^2)./(2*wavelength^2)-y.*thetay.*z/w0y^2));

a=intax.*intay;
%Z=-10*log10(a);