| 查看: 1826 | 回复: 2 | ||
[求助]
利用Matlab计算DBR反射率的程序
|
| 有没有大神用Matlab 计算过垂直腔面发射激光器（VCSEL）DBR反射率啊？老师让我自己编写程序计算不同波长下的折射率，（最后会出现一个图 波长---折射率），怎么弄都出不来图啊，我现在只能算出来特定波长、特定层数的DBR反射率，不同波长下的怎么编写程序都出不来图啊，关键就是要出现那个波长---折射率的图，有没有大神弄过这些？帮帮我啊，，，，，真的是心累，真的是快被代码程序整疯了。。。。。谢谢谢谢！！！ |
回复此楼» 猜你喜欢
浙大物院表面物理课题组长期招收普博及直博学生
已经有259人回复
-
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
物理学I论文润色/翻译怎么收费?
已经有84人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
课题组招1-2名硕士
已经有1人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有0人回复
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31030.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5603
- 在线: 856.8小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
2楼2019-08-11 19:04:15
|
>> n0=1; a0=0*pi/180; %入射角 n1=3.1; %介质1的折射率 n2=3.5; %介质2的折射率 d0=850; %波长 d1=d0/4/n1; %介质1的厚度 d2=d0/4/n2; %介质1的厚度 N=28; %层数 a1=asin(n0*sin(a0)/n1); a2=asin(n1*sin(a1)/n2); a3=asin(n2*sin(a2)/n0); p1=2*pi*n1*d1*cos(a1)/d0; p2=2*pi*n2*d2*cos(a2)/d0; 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]; A2=cos(p2);B2=-i*sin(p2)/c2;C2=-i*c2*sin(p2);D2=cos(p2); M2=[A2 B2;C2 D2]; M=(M1*M2)^N; 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); 这就是传输矩阵模型，只能算出来特定波长的，有什么建议吗？不同波长的不能 算出来，谢谢啦! |
3楼2019-08-12 08:42:27
