直波导与弯曲波导耦合功率计算.docxVIP

直波导与弯曲波导耦合功率计算 直波导与弯曲波导耦合功率计算 PAGE 直波导与弯曲波导耦合功率计算 设计工作三 clear all n1=;n2=;n3=n2;n4=n2;n5=n2;R=5000;a=2;k0=2*pi/; 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=(n1^2-n2^2)/(2*n1^2); V=k0*a*sqrt(n1^2-n2^2); ax=sqrt(k0^2*(n1^2-n2^2)-kx^2); d=::15; A0=1;B0=0; for z=-316::316 l=length(z); D=5000-sqrt(5000^2-z^2)+d+a; k=(sqrt(2*delta)*(kx*a)^2*(ax*a)^2*exp(-ax*(D-2*a)))/(a*(1+ax*a)*V^3); A=A0.*cos(k*-1i*B0.*sin(k*;%????·??￠?×è￥???? B=-1i*A0.*sin(k*+B0.*cos(k*; A0=A;B0=B; end Pb(:,1)=(abs(B0)).^2; figure plot(d,Pb) axis([ 10 0 1]); set(gca,XTick,(0::10));%éè??òa??ê?×?±ê?ì?è title(í??ú2¨μ???d±??ˉμ?1éò??ˉ1|?ê); xlabel(d(|ìm)); ylabel(P); 计算矩形波导与弯曲波导的耦合光功率,利用微元法将弯曲波导看成是很多段小的矩形波导相加起来,则由同相耦合器的公式 A B 因为两波导传播常数相同，所以δ=0，q=κ。两矩形波导的 κ= V= 而kx可由玛卡蒂尼近似方法求出 在近似时，由于当两波导相距超过10μm时，耦合功率可几乎为零，由设计工作二中得到波导芯长宽与模式个数的关系图可知在单模状况下取波导芯 a=2 μm比较合适，则最终可得图 由图可看出当d从0开始逐渐远离时，到，, μm时最高的B端口光功率最大，耦合效率最高，其后随着距离的增大逐渐减弱到0.