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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 计算模拟 » 程序代码 » Matlab 计算DBR折射率
51/1返回列表
查看: 1883  |  回复: 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的回帖

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的回帖
查看全部 6 个回答

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的回帖
查看全部 6 个回答
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭