传输矩阵的Matlab简易编程o.doc

下载文档

156
0
约4.59千字
约 5页
2016-08-22 发布于浙江
举报
版权申诉
保障服务

传输矩阵的Matlab简易编程o.doc

1、本文档共5页，可阅读全部内容。
2、有哪些信誉好的足球投注网站（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击这里二次下载。
4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。

传输矩阵的Matlab简易编程o

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=[A1 B1;C1 D1];%介质1的传输矩阵 A2=cos(p2);B2=-i*sin(p2)/c2;C2=-i*c2*sin(p2);D2=cos(p2); M2=[A2 B2;C2 D2];%介质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:);%初始入射角 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;%真空磁导率 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=[A1 B1;C1 D1];%介质1的传输矩阵 A2=cos(p2);B2=-i*sin(p2)/c2;C2=-i*c2*sin(p2);D2=cos(p2); M2=[A2 B2;C2 D2];%介质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:);%初始入射角 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;%真空磁导率 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=[A1 B1;C1 D1];%介质1的传输矩阵 A2=cos(p2);B2=