应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » Matlab 计算DBR折射率
61/1返回列表
查看: 1963  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

huishujie

新虫 (小有名气)

[求助] 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);
回复此楼

» 猜你喜欢
1楼 2019-08-13 14:19:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plasmaSim

铁虫 (小有名气)
把原理公式也写出来

发自小木虫IOS客户端
赞一下 回复此楼
等离子体放电、电磁波与等离子体相互作用等仿真
2楼2019-08-14 06:33:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

独孤神宇

版主 (知名作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +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);
赞一下(2人) 回复此楼
数值计算
3楼2019-08-14 08:22:01
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huishujie

新虫 (小有名气)
感谢感谢，可以算出来了，后面有啥问题 再请教老兄！
赞一下 回复此楼
4楼2019-08-14 09:55:40
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

plasmaSim

铁虫 (小有名气)
引用回帖:
4楼: Originally posted by huishujie at 2019-08-14 09:55:40
感谢感谢，可以算出来了，后面有啥问题 再请教老兄！

看起来问题是这个R(i)=real(R);重复使用了R？
赞一下 回复此楼
等离子体放电、电磁波与等离子体相互作用等仿真
5楼2019-08-14 11:07:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

huishujie

新虫 (小有名气)
引用回帖:
5楼: Originally posted by plasmaSim at 2019-08-14 11:07:21
看起来问题是这个R(i)=real(R);重复使用了R？...

对的呢，就是这个问题。
赞一下 回复此楼
6楼2019-08-14 12:45:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 huishujie 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 0856材料与化工，270求调剂 +16 YXCT 2026-03-01 18/900 2026-03-04 22:05 by 学zh
[考研] 化学工程求调剂 +10 化工人999 2026-03-04 10/500 2026-03-04 21:21 by zhukairuo
[考研] 267化工调剂求助 +5 聪少OZ 2026-03-04 5/250 2026-03-04 20:33 by kakakapanpan
[考研] 一志愿中石油，化学工程，301求调剂 +3 豚馒枳枳 2026-03-04 3/150 2026-03-04 18:32 by houyaoxu
[考研] 282求调剂 +5 2103240126 2026-03-02 8/400 2026-03-04 17:29 by 妄深于海
[考研] 化学 0703求调剂 总分293 一志愿211 +8 土土小虫 2026-03-03 8/400 2026-03-04 16:45 by xwxstudy
[考研] 本科太原理工采矿工程，求调剂 +3 onlx 2026-03-01 3/150 2026-03-04 15:57 by Stephen_ym
[考研] 347分电子信息专硕求调剂 +3 咕嘟咕嘟1 2026-03-03 3/150 2026-03-04 12:32 by 热情沙漠
[考研] 289求调剂 +8 yang婷 2026-03-02 10/500 2026-03-04 09:10 by 无际的草原
[考研] 0856材料工程，初试313调剂 +7 卖个关子吧 2026-03-03 7/350 2026-03-03 21:33 by L135790
[考研] 理学，工学，农学调剂，少走弯路，这里欢迎您！ +8 likeihood 2026-03-02 11/550 2026-03-03 19:39 by 一叶凡尘?
[考研] 085600材料与化工调剂 280分 +10 yyqqhh 2026-03-03 10/500 2026-03-03 18:37 by 386661696
[考研] 283求调剂 +9 鹿沫笙 2026-03-02 9/450 2026-03-03 15:58 by tgxtgxtgx9
[考研] 计算机学硕 分数285求调剂 +4 glwshine 2026-03-02 5/250 2026-03-03 14:27 by king呀
[考研] 307求调剂 +6 wyyyqx 2026-03-01 6/300 2026-03-03 09:24 by 2235787770
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 16/800 2026-03-03 00:36 by 好好好1233
[考研] 一志愿山东大学材料与化工325求调剂 +5 半截的诗0927 2026-03-02 5/250 2026-03-02 18:37 by 明亮9527
[考研] 295求调剂 +8 19171856320 2026-02-28 8/400 2026-03-02 11:19 by yuchj
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 0856材料求调剂 +4 麻辣鱿鱼 2026-02-28 4/200 2026-03-01 16:51 by caszguilin
信息提示
关闭
请填处理意见
关闭