 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
小木蟲論壇-學術科研互動平臺 » 計算模擬區(qū) » 計算模擬 » 程序代碼 » Matlab 計算DBR折射率
| 5 | 1/1 | 返回列表 |
| 查看: 1970 | 回復: 5 | ||
| 當前只顯示滿足指定條件的回帖，點擊這里查看本話題的所有回帖 | ||
huishujie新蟲 (小有名氣)
|
[求助]
Matlab 計算DBR折射率 已有1人參與
|
|
|
根據(jù)傳輸矩陣法試著寫了一下程序，計算不出來，糾結找原因好幾天了，真的是找不到原因，因為我自己也是剛開始學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人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經有0人回復
-
物理學I論文潤色/翻譯怎么收費?
已經有102人回復
-
【新加坡】納米電子器件項目組有“聯(lián)合培養(yǎng)博士生”名額
已經有0人回復
-
幫我的英語口語老師找學生
已經有3人回復
-
計算三階力常數(shù)用thirdorder擴胞時參數(shù)如何選擇？臨近原子數(shù)怎么取值
已經有5人回復
-
2026年循環(huán)經濟功能材料國際會議（ICFMCE 2026）
已經有0人回復
-
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）
已經有0人回復
5樓2019-08-14 11:07:21
2樓2019-08-14 06:33:11
獨孤神宇
版主 (知名作家)
- 應助: 490 (碩士)
- 貴賓: 0.008
- 金幣: 31016.3
- 散金: 802
- 紅花: 122
- 沙發(fā): 1
- 帖子: 5600
- 在線: 856.5小時
- 蟲號: 3522474
- 注冊: 2014-11-06
- 性別: GG
- 專業(yè): 機械動力學
- 管轄: 計算模擬
【答案】應助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感謝參與，應助指數(shù) +1
月只藍: 金幣+5, 感謝指導！ 2019-08-14 08:33:19
huishujie: 金幣+10, ★有幫助 2019-08-14 09:54:34
感謝參與，應助指數(shù) +1
月只藍: 金幣+5, 感謝指導！ 2019-08-14 08:33:19
huishujie: 金幣+10, ★有幫助 2019-08-14 09:54:34
|
取大步長 降低矩陣次方 試了一下，可以運行 ，你可以改回去 慢慢算。。 另外，建議不要用 循環(huá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
新蟲 (小有名氣)
- 應助: 0 (幼兒園)
- 金幣: 340
- 帖子: 103
- 在線: 14.7小時
- 蟲號: 14243178
- 注冊: 2019-02-19
- 專業(yè): 半導體光電子器件
4樓2019-08-14 09:55:40
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考博] 售瀆SCI文章，我:8⊙ 55 1⊙ 54。備注【⊙=0】 +4 | ipys00zm58 2026-03-06 | 6/300 |
|
|---|---|---|---|---|
|
[考研] 一志愿中石油，化學工程，301求調劑 +4 | 豚饅枳枳 2026-03-04 | 4/200 |
|
|
|
[考研] 招調劑 +12 | 帆船哥 2026-03-04 | 14/700 |
|
|
|
[考研] 哈爾濱理工大學2026年研究生調劑，材料科學與化學工程學院研究生調劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
|
[考博] 2026年博士名額撿漏 +4 | 科研ya 2026-03-04 | 7/350 |
|
|
|
[考研] 085600材料與化工（高分子）290分求調劑 +9 | wengyujian 2026-03-04 | 9/450 |
|
|
|
[考研] 化學專業(yè)調劑 +5 | 好好好1233 2026-03-04 | 6/300 |
|
|
|
[考研] 材料學碩318求調劑 +13 | February_Feb 2026-03-01 | 13/650 |
|
|
|
[考研] 0856材料與化工求調劑！ +5 | 化工考生111 2026-03-04 | 11/550 |
|
|
|
[考研] 考研282分求調劑，接受跨專業(yè) +4 | 劉淄博 2026-03-04 | 7/350 |
|
|
|
[考研] 一志愿蘇大材料工程專碩293求調劑 +3 | 鋼鐵大炮 2026-03-04 | 3/150 |
|
|
|
[考研] 沒上岸的看過來 +3 | tangxiaotian 2026-03-01 | 5/250 |
|
|
|
[考研] 0703化學 學碩 理工科均可 不區(qū)分研究方向 總分279求調劑 +7 | 1一11 2026-03-03 | 7/350 |
|
|
|
[考研] 325求調劑 +5 | 學家科 2026-03-04 | 5/250 |
|
|
|
[考研] 求調劑 +11 | yunziaaaaa 2026-03-01 | 13/650 |
|
|
|
[考研] 0856求調劑285 +11 | 呂仔龍 2026-02-28 | 11/550 |
|
|
|
[考博] 博士自薦 +4 | kkluvs 2026-02-28 | 5/250 |
|
|
|
[考研] 303求調劑 +5 | 今夏不夏 2026-03-01 | 5/250 |
|
|
|
[考研] 265分求調劑不調專業(yè)和學校有行學上就 +6 | 禮堂丁真258 2026-02-28 | 9/450 |
|
|
|
[考研] 調劑 +3 | 簡木ChuFront 2026-02-28 | 3/150 |
|