matlab---可运行--色散图和透射图

透过率图•clc;clear;•na=2.10;•nb=1.46;•n1=1;•n2=1;•ford=400:1600•c3=0;c1=asin(n1*sin(c3)/na);c2=asin(na*sin(c1)/nb);c4=asin(nb*sin(c2)/n2);•d1=1064;•a=d1/(4*na);b=d1/(4*nb);•Bb=2*pi*nb*b*cos(c2)/d;•f=4*pi*1e-7;•e=1e-9/(36*pi);•m=sqrt(e/f);•za=m*cos(c1)*na;zb=m*cos(c2)*nb;z1=f*cos(c3)*n1;z2=f*cos(c4)*n2;•p1=cos(Bb);p2=-i*sin(Bb)/zb;p3=-i*zb*sin(Bb);p4=cos(Bb);•P=[p1p2;p3p4];•q1=cos(Ba);q2=-i*sin(Ba)/za;q3=-i*za*sin(Ba);q4=cos(Ba);•Q=[q1q2;q3q4];•O=Q*P;O1=O^27;•O11=O1(1,1);O12=O1(1,2);O13=O1(2,1);O14=O1(2,2);•z1=sqrt(e/f)*n1*cos(c3);z2=sqrt(e/f)*n2*cos(c4);•t=2*z1/(z1*(O11+z2*O12)+O13+z2*O14);•r=(z1*(O11+z2*O12)-O13-z2*O14)/(z1*(O11+z2*O12)+O13+z2*O14);•t1=abs(t);•r1=abs(r);•s=d-399;•k(1,s)=t1•end•d=400:1600;•plot(d,k);色散图•%不同波导芯厚度h对应的TE模式也TM模式色散图程序•k0=2*pi/1.55;•form=0:1:3•n1=1.537;%TE模式下的折射率分布•n2=1.444;•n3=1.51;•N1=1.510:0.0001:1.537;%TE模式下有效折射率的范围•%TE模式下波导芯厚度•ifn1.^2-N1.^2•h_te=(m*pi.*ones()+atan(((N1.^2-n2.^2)./(n1.^2-N1.^2)).^0.5)+atan(((N1.^2-n3.^2)./(n1.^2-N1.^2)).^0.5))./(k0.*(n1.^2-N1.^2).^0.5);•n1=1.530;%TM模式下的折射率分布•n2=1.444;•n3=1.5095;•N2=1.5095:0.0001:1.530;%TM模式下有效折射率的范围•%TM模式下波导芯厚度•h_tm=(m*pi.*ones()+atan((((N2.^2-n2.^2)./(n1.^2-N2.^2)).^0.5).*n1.^2./n2.^2)+atan((((N2.^2-n3.^2)./(n1.^2-N2.^2)).^0.5).*n1.^2./n3.^2))./(k0.*(n1.^2-N2.^2).^0.5);•title('色散图')•plot(h_te,N1,'-b',h_tm,N2,'-r')•axis([0,15,1.509,1.537])•legend('TE','TM','Location','Northwest')•xlabel('h(um)','fontsize',10)•ylabel('N','fontsize',10)•holdon•gridon•gridminor•end