计算AWG波分复用器中平板波导连接接口在耦合到下级光纤时

1.问题描述：

 计算AWG波分复用器中平板波导连接接口在耦合到下级光纤时

2.部分程序：

function yanshexiaolv
k0=4.051268e6;k1=4.6359e+005;
k2=6.038715e6;q=6.0799e+005;
syms r K;a1=6e-6;a2=11e-6;a3=15.5e-6;
M=[0:160];f=7948.707992e-6;d=15.5e-6;
E1=cos(k2*a1*r/2)-k2/q*(cos(19.1161*r)+cos(17.1161*r));
E3=cos(k2*a3*r/2)-k2/q*(cos(47.8*r)+cos(45.8*r));
P1=E1^2;P3=E3^2;
X1=int(P1,r,(K*d/f-a1/2/f),(K*d/f+a1/2/f));
X3=int(P3,r,(K*d/f-a3/2/f),(K*d/f+a3/2/f));
Q=double(symsum(X3,K,0,160));
Y1=[];
for n=M
    Y1=[Y1,double(symsum(X1,K,0,n))];
end
G1=double(Y1/Q);
plot(M,G1);grid on;title('w=6μm时输入平板波导的衍射效率');
xlabel('阵列波导数2M+1');ylabel('输入平板波导中的归一化衍射效率');

3.仿真结论：

D00003