matlab 可运行 色散图和透射图

透过率图

• clc;clear;

• na=2.10;

• nb=1.46;

• n1=1 ;

• n2=1 ;

• for d=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=[p1 p2;p3 p4];

• q1=cos(Ba);q2=-i*sin(Ba)/za;q3=-i*za*sin(Ba);q4=cos(Ba);

• Q=[q1 q2;q3 q4];

• 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;

• for m=0:1:3

• n1=1.537; %TE模式下的折射率分布

• n2=1.444;

• n3=1.51;

• N1=1.510:0.0001:1.537; %TE模式下有效折射率的范围

• %TE模式下波导芯厚度

• if n1.^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)

• hold on

• grid on

• grid minor

• end