没有合适的资源?快使用搜索试试~ 我知道了~
首页传输矩阵对应的MATLAB仿真程序.docx
传输矩阵对应的MATLAB仿真程序.docx
需积分: 47 81 下载量 152 浏览量
更新于2023-03-03
评论 11
收藏 19KB DOCX 举报
MATLAB程序对一维光子晶体模型的仿真在光子晶体研究方法上主要采用特性传输矩阵法,平面波展开法,球面波展开法、有限时域差分方程法和N阶(Order-N)法等。这些方法各有优缺点。 传输矩阵法简单实用,但是传统的传输矩阵法简化了非垂直入射以及多重散射的情况。
资源详情
资源评论
资源推荐

传输矩阵的 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)% 透射系数



















lanseyuguang
- 粉丝: 1
- 资源: 2
上传资源 快速赚钱
我的内容管理 收起
我的资源 快来上传第一个资源
我的收益
登录查看自己的收益我的积分 登录查看自己的积分
我的C币 登录后查看C币余额
我的收藏
我的下载
下载帮助

会员权益专享
安全验证
文档复制为VIP权益,开通VIP直接复制

评论0