传输矩阵的Matlab简易编程o

《传输矩阵的Matlab简易编程o》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、n0=1;%n0为空气折射率a0=0;%初始入射角为0，即正入射n1=1.444;n2=1.7514;d1=3.4886;d2=4.6373;d=1.55;p1=2*pi*n1*d1/d;p2=2*pi*n2*d2/d;u0=4*pi*1e-7;%真空磁导率e0=1e-9/(36*pi);%真空介电常数b=sqrt(e0/u0);c1=b*n1;%波阻抗c2=b*n2;c0=b*n0;c00=b*n0;A1=cos(p1);B1=-i*sin(p1)/c1;C1=-i*c1*sin(p1);D1=cos(p1);M1=[A1B1;C1D1];%介质1的传输矩阵A2=cos(p2);B2=

2、-i*sin(p2)/c2;C2=-i*c2*sin(p2);D2=cos(p2);M2=[A2B2;C2D2];%介质2的传输矩阵M=M1*M2;%两个介质层的总传输矩阵A=M(1,1);B=M(1,2);C=M(2,1);D=M(2,2);r=(A*c0+B*c0*c00-C-D*c00)/(A*c0+B*c0*c00+C+D*c00)%反射系数t=(2*c0)/(A*c0+B*c0*c00+C+D*c00)%透射系数R=r*conj(r)T=t*conj(t)以上可作为TE波和TM波的一特列，即垂直入射n0=1;%n0为空气折射率a0=input('请输入入射角a0:');%初始入

3、射角n1=input('请输入介质1的折射率n1:');n2=input('请输入介质2的折射率n2:');d1=input('请输入介质1的厚度d1:');d2=input('请输入介质2的厚度d2:');d=input('请输入入射光波长d:');a1=asin(n0*sin(a0)/n1);%光在介质1中的传播角度a2=asin(n1*sin(a1)/n2);%光在介质2中的传播角度a3=asin(n2*sin(a2)/n0);%射出介质时的出射角p1=2*pi*n1*d1*cos(a1)/d;%相位厚度p2=2*pi*n2*d2*cos(a2)/d;u0=4*pi*1e-7;%真

4、空磁导率e0=1e-9/(36*pi);%真空介电常数b=sqrt(e0/u0);c1=b*n1*cos(a1);%波阻抗c2=b*n2*cos(a2);c0=b*n0*cos(a0);c00=b*n0*cos(a3);A1=cos(p1);B1=-i*sin(p1)/c1;C1=-i*c1*sin(p1);D1=cos(p1);M1=[A1B1;C1D1];%介质1的传输矩阵A2=cos(p2);B2=-i*sin(p2)/c2;C2=-i*c2*sin(p2);D2=cos(p2);M2=[A2B2;C2D2];%介质2的传输矩阵M=M1*M2;%两个介质层的总传输矩阵A=M(1,1

5、);B=M(1,2);C=M(2,1);D=M(2,2);r=(A*c0+B*c0*c00-C-D*c00)/(A*c0+B*c0*c00+C+D*c00)%反射系数t=(2*c0)/(A*c0+B*c0*c00+C+D*c00)%透射系数R=r*conj(r)T=t*conj(t)以上为TE波的例子。若为TM波则，只有波阻抗表达不同，即：n0=1;%n0为空气折射率a0=input('请输入入射角a0:');%初始入射角n1=input('请输入介质1的折射率n1:');n2=input('请输入介质2的折射率n2:');d1=input('请输入介质1的厚度d1:');d2=inpu

6、t('请输入介质2的厚度d2:');d=input('请输入入射光波长d:');a1=asin(n0*sin(a0)/n1);%光在介质1中的传播角度a2=asin(n1*sin(a1)/n2);%光在介质2中的传播角度a3=asin(n2*sin(a2)/n0);%射出介质时的出射角p1=2*pi*n1*d1*cos(a1)/d;%相位厚度p2=2*pi*n2*d2*cos(a2)/d;u0=4*pi*1e-7;%真空磁导率e0=1e-9/(36*pi);%真空介电常数b=sqrt(e0/u0);c1=b*n1/cos(a1);%波阻抗c2=b*n2/cos(a2);c0=b*n0/c

7、os(a0);c00=b*n0/cos(a3);A1=cos(p1);B1=-i*sin(p1)/c1;C1=-i*c1*sin(p1);D1=cos(p1);M1=[A1B1;C1D1];%介质1的传输矩阵A2=cos(p2);B2=-i*sin(p2)/c2;C2=-i*c2*sin(p2);D2=cos(p2);M2=[A2B2;C2D2];%介质2的传输矩阵M=M1*M2;%两个介质层的总传输矩阵A=M(1,1);B=M(1