直波导与弯曲波导耦合功率计算(共3页).docx

设计工作三clear alln1=1.5150;n2=1.5000;n3=n2;n4=n2;n5=n2;R=5000;a=2;k0=2*pi/1.55;F=(x)(x(1)*a-atan(sqrt(k0*n1)2-(k0*n2)2-x(1)2)./x(1)-atan(sqrt(k0*n1)2-(k0*n4)2-x(1)2)./x(1); x(2)*a-atan(sqrt(n1/n3)2*sqrt(k0*n1)2-(k0*n3)2-x(2)2)/x(2)-atan(n1/n3)2*sqrt(k0*n1)2-(k0*n5)2-x(2)2)./x(2);tmp=fsolve(F,0 5);kx=tmp(1);delta=(n12-n22)/(2*n12);V=k0*a*sqrt(n12-n22);ax=sqrt(k02*(n12-n22)-kx2);d=0.1:0.1:15;A0=1