 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1975 | 回復(fù): 5 | |||
huishujie新蟲 (小有名氣)
|
[求助]
Matlab 計算DBR折射率 已有1人參與
|
|
根據(jù)傳輸矩陣法試著寫了一下程序，計算不出來，糾結(jié)找原因好幾天了，真的是找不到原因，因?yàn)槲易约阂彩莿傞_始學(xué)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); |
回復(fù)此樓» 猜你喜歡
【新加坡】納米電子器件項(xiàng)目組有“聯(lián)合培養(yǎng)博士生”名額
已經(jīng)有0人回復(fù)
-
幫我的英語口語老師找學(xué)生
已經(jīng)有4人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有55人回復(fù)
-
計算三階力常數(shù)用thirdorder擴(kuò)胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經(jīng)有7人回復(fù)
-
2026年循環(huán)經(jīng)濟(jì)功能材料國際會議（ICFMCE 2026）
已經(jīng)有0人回復(fù)
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經(jīng)有0人回復(fù)
-
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有2人回復(fù)
plasmaSim
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 424
- 帖子: 120
- 在線: 18.2小時
- 蟲號: 15510952
- 注冊: 2019-07-04
- 專業(yè): 等離子體物理
2樓2019-08-14 06:33:11
獨(dú)孤神宇
版主 (知名作家)
- 應(yīng)助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機(jī)械動力學(xué)
- 管轄: 計算模擬
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣+5, 感謝指導(dǎo)！ 2019-08-14 08:33:19
huishujie: 金幣+10, ★有幫助 2019-08-14 09:54:34
感謝參與，應(yīng)助指數(shù) +1
月只藍(lán): 金幣+5, 感謝指導(dǎo)！ 2019-08-14 08:33:19
huishujie: 金幣+10, ★有幫助 2019-08-14 09:54:34
|
取大步長 降低矩陣次方 試了一下，可以運(yùn)行 ，你可以改回去 慢慢算。。 另外，建議不要用 循環(huán)，直接 向量化 計算可以加快運(yùn)算效率 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
huishujie
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 340
- 帖子: 103
- 在線: 14.7小時
- 蟲號: 14243178
- 注冊: 2019-02-19
- 專業(yè): 半導(dǎo)體光電子器件
4樓2019-08-14 09:55:40
plasmaSim
鐵蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 424
- 帖子: 120
- 在線: 18.2小時
- 蟲號: 15510952
- 注冊: 2019-07-04
- 專業(yè): 等離子體物理
5樓2019-08-14 11:07:21
huishujie
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 340
- 帖子: 103
- 在線: 14.7小時
- 蟲號: 14243178
- 注冊: 2019-02-19
- 專業(yè): 半導(dǎo)體光電子器件
6樓2019-08-14 12:45:43
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 huishujie 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 歡迎211本科同學(xué)，過A區(qū)國家線，A區(qū)非偏遠(yuǎn)一本，交叉學(xué)科課題組 +22 | lisimayy 2026-03-04 | 35/1750 |
|
|---|---|---|---|---|
|
[考研] 086000生物與醫(yī)藥319分求調(diào)劑 +4 | Tolkien 2026-03-07 | 4/200 |
|
|
|
[考研] 材料與化工求調(diào)劑 +8 | 與冬清寧 2026-03-07 | 10/500 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實(shí)人 2026-03-09 | 4/200 |
|
|
|
[考研] 一志愿東北大學(xué)080500，英一數(shù)二343分求調(diào)劑 +7 | Seele_v5 2026-03-04 | 7/350 |
|
|
|
[考研] 293一志愿華東理工 0817化學(xué)工程與技術(shù) 調(diào)劑 +5 | fjj0912 2026-03-07 | 5/250 |
|
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
|
[考研] 290求調(diào)劑 +6 | Bananaiy 2026-03-04 | 6/300 |
|
|
|
[考研] 安徽農(nóng)業(yè)大學(xué)材料與化學(xué)學(xué)院0856材料招收調(diào)劑 +3 | akakk47 2026-03-05 | 9/450 |
|
|
|
[考研]
 化學(xué)工程與技術(shù)316求調(diào)劑
+7
|
Sixuan wang 2026-03-06 | 7/350 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)312分求調(diào)劑 +7 | T123 tt 2026-03-04 | 7/350 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | G212 2026-03-03 | 6/300 |
|
|
|
[考研] 一志愿南京大學(xué)資源與環(huán)境專業(yè)283分求調(diào)劑 +8 | 追云?? 2026-03-04 | 8/400 |
|
|
|
[考研] 287求調(diào)劑 +3 | 看看我. 2026-03-05 | 6/300 |
|
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財 2026-03-05 | 5/250 |
|
|
|
[考研] 080500材料科學(xué)與工程 +13 | 202114020319 2026-03-03 | 13/650 |
|
|
|
[考研] 求調(diào)劑，學(xué)校研究所都可以，材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
|
[考研] 一志愿武理085601專碩347分求調(diào)劑 +4 | 啊歐歐歐 2026-03-04 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +7 | 博斯特525 2026-03-04 | 7/350 |
|
|
|
[考研] 298求調(diào)劑 +3 | 人間唯你是清歡 2026-03-03 | 4/200 |
|