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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » Matlab 计算DBR折射率
61/1返回列表
查看: 1964  |  回复: 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返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 欢迎211本科同学，过A区国家线，A区非偏远一本，交叉学科课题组 +7 lisimayy 2026-03-04 11/550 2026-03-05 00:06 by 上岸上岸上岸_－
[考研] 求调剂院校 +6 云朵452 2026-03-02 12/600 2026-03-04 23:17 by 云朵452
[考研] 一志愿武汉理工大学-085602-总分296分-求调剂 +7 紫川葡柚 2026-03-04 7/350 2026-03-04 21:04 by kakakapanpan
[考研] 070300调剂一志愿南京大学化学专业，初试273分，有论文有专利有竞赛，求调剂B区211 +3 82206202 2026-02-28 8/400 2026-03-04 19:24 by 每天只摆一小会
[考研] 282求调剂 +5 2103240126 2026-03-02 8/400 2026-03-04 17:29 by 妄深于海
[考研] 0855机械工程324求调剂 +3 郑SK 2026-03-02 7/350 2026-03-04 17:24 by congmingdou
[考研] 085701环境工程 求调剂 +6 xiiiia 2026-03-04 6/300 2026-03-04 16:45 by yuningshan
[考研] 0703化学求调剂 +3 将未有期 2026-03-04 3/150 2026-03-04 16:19 by sslc1985
[考研] 一志愿314求调剂 +7 202111120625 2026-03-03 7/350 2026-03-04 15:56 by zhukairuo
[考研] 347分电子信息专硕求调剂 +3 咕嘟咕嘟1 2026-03-03 3/150 2026-03-04 12:32 by 热情沙漠
[考研] 0703化学调剂 +4 G212 2026-03-03 5/250 2026-03-04 09:34 by 每天只摆一小会
[考研] 复试调剂 +4 杰尼龟aaa 2026-03-03 4/200 2026-03-03 23:06 by zhukairuo
[考研] 材料学硕318求调剂 +15 February_Feb 2026-03-01 17/850 2026-03-03 14:47 by xiaomc_gzh
[考研] 调剂 +5 13853210211 2026-03-02 7/350 2026-03-03 11:42 by 13853210211
[考研] 化工335求调剂 +5 摸摸猫猫头 2026-03-02 5/250 2026-03-03 09:10 by houyaoxu
[考研] 0856材料调剂 +5 沿岸有贝壳OUC 2026-03-02 5/250 2026-03-02 20:31 by hypershenger
[考研] 085600材料工程一志愿中科大总分312求调剂 +9 吃宵夜1 2026-02-28 11/550 2026-03-02 20:14 by hypershenger
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师！
[考研] 304求调剂 +6 曼殊2266 2026-02-28 7/350 2026-03-01 15:14 by wjLi2017
[考研] 304求调剂 +3 52hz~~ 2026-02-28 5/250 2026-03-01 00:00 by 52hz~~
信息提示
关闭
请填处理意见
关闭