Calculate when the slab waveguide connection interface in the AWG wavelength division multiplexer is coupled to the lower-level fiber

1. Problem description:

 Calculate when the slab waveguide connection interface in the AWG wavelength division multiplexer is coupled to the lower-level fiber

2. Part of the program:

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 when input the diffraction efficiency of the slab waveguide' );
xlabel('number of arrayed waveguides 2M+1'); ylabel('input normalized diffraction efficiency in slab waveguide');

3. Simulation conclusion:

D00003