function gonglvfenbu
k0=4.051268e6;k1=4.6359e+005;
k2=6.038715e6;q=6.0799e+005;
r=[-0.8:pi/50:0.8];a=4.2e-6;
A=q*cos(k2*a*sin(r)/2);
B=k2*sin(r).*sin(k2*a*sin(r)/2);
C=A-B;
D=k1^2*q*cos(r);
F=k1^2-k2^2*sin(r).*sin(r);
G=q^2+k2^2*sin(r).*sin(r);
E=D.*C./F./G;
P=E.^2;
plot(r,P);title('平板波导内功率分布');
xlabel('衍射角θ (a=6μm)');ylabel('输入平板波导中的功率分布P');
a=4.1e-6;
A=q*cos(k2*a*sin(r)/2);
B=k2*sin(r).*sin(k2*a*sin(r)/2);
C=A-B;E=D.*C./F./G;
P=E.^2;
plot(r,P);title('平板波导内功率分布');
xlabel('衍射角θ (a=5μm)');ylabel('输入平板波导中的功率分布P');
a=4e-6;
A=q*cos(k2*a*sin(r)/2);
B=k2*sin(r).*sin(k2*a*sin(r)/2);
C=A-B;E=D.*C./F./G;
P=E.^2;
plot(r,P);title('平板波导内功率分布');
xlabel('衍射角θ (a=4μm)');ylabel('输入平板波导中的功率分布P');
评论0