| 查看: 1988 | 回复: 5 | ||
[求助]
Matlab 计算DBR折射率 已有1人参与
|
|
根据传输矩阵法试着写了一下程序，计算不出来，纠结找原因好几天了，真的是找不到原因，因为我自己也是刚开始学Matlab, 所以有么有大神给看看到底哪里有问题啊？真心挺着急的。。。。。 n0=1; a0=0*pi/180; n1=3.1; n2=3.6; N=28; Lamda=(700:1:900); for i=1:length(Lamda) lamda=Lamda(i); d1=68.5; d2=59.02; 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)/lamda; p2=2*pi*n2*d2*cos(a2)/lamda; 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); R(i)=real(R); end figure; plot(Lamda,R); |
回复此楼» 猜你喜欢
欢迎加入课题组
已经有0人回复
-
散金币，求好运，祝面上顺利！
已经有80人回复
-
物理学I论文润色/翻译怎么收费?
已经有130人回复
-
华东师范大学芯片设计徐珑真实水平咋样
已经有1人回复
-
调剂
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有0人回复
-
中德博士后2026项目5月30日截止- 国内博士后入站、德国两年全职
已经有0人回复
-
面向全球申请人的中德博士后2026项目5月30日截止- 国内博后入站德国两年全职
已经有0人回复
-
【新加坡】纳米电子器件项目组有“联合培养博士生”名额
已经有2人回复
-
重庆交大材料学院26年硕士招生复试实施细则已出。欢迎加入光子学微结构与器件课题组。
已经有0人回复
-
北工大宋晓艳教授团队招收2026年9月入学博士【学博和专博名额均有】！
已经有1人回复
2楼2019-08-14 06:33:11
独孤神宇
版主 (知名作家)
- 应助: 490 (硕士)
- 贵宾: 0.008
- 金币: 31016.3
- 散金: 802
- 红花: 122
- 沙发: 1
- 帖子: 5600
- 在线: 856.5小时
- 虫号: 3522474
- 注册: 2014-11-06
- 性别: GG
- 专业: 机械动力学
- 管辖: 计算模拟
【答案】应助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
月只蓝: 金币+5, 感谢指导！ 2019-08-14 08:33:19
huishujie: 金币+10, ★有帮助 2019-08-14 09:54:34
感谢参与，应助指数 +1
月只蓝: 金币+5, 感谢指导！ 2019-08-14 08:33:19
huishujie: 金币+10, ★有帮助 2019-08-14 09:54:34
|
取大步长 降低矩阵次方 试了一下，可以运行 ，你可以改回去 慢慢算。。 另外，建议不要用 循环，直接 向量化 计算可以加快运算效率 n0=1; a0=0*pi/180; n1=3.1; n2=3.6; N=2; Lamda=700:5:900; RR=zeros(1,length(Lamda)); for i=1:length(Lamda) lamda=Lamda(i); d1=68.5; d2=59.02; 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)/lamda; p2=2*pi*n2*d2*cos(a2)/lamda; 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); RR(i)=real(R); end figure; plot(Lamda,RR); |
3楼2019-08-14 08:22:01
4楼2019-08-14 09:55:40
5楼2019-08-14 11:07:21
6楼2019-08-14 12:45:43
