版塊導(dǎo)航: 正在加載中...

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號能夠正常使用，請盡快對帳號進行手機號驗證，感謝您的理解與支持！

24小時熱門版塊排行榜

返回列表

惜夢尋草

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時
蟲號: 2118505
注冊: 2012-11-10
性別: MM
專業(yè): 光學(xué)

[求助] 利用matlab畫出級聯(lián)OPO波長調(diào)諧曲線

求高手指點！如題，模擬PPLN晶體級聯(lián)OPO+OPA溫度調(diào)諧曲線，但是不知道該怎么模擬！要得到最后面的那個圖。
下面是我的程序，希望哪位高人能幫忙調(diào)一下！十分感激！

clear
clc
n=1;
m=1;
a1=5.756;
a2=0.0983;
a3=0.2020;
a4=189.32;
a5=12.52;
a6=1.32e-2;
b1=2.860e-6;
b2=4.700e-8;
b3=6.113e-8;
b4=1.516e-4;
A=31.5;
lamdap=1.064;
T = [25:0.1:250];
lamdas = [1.065:0.01:3.2];
lamdass = [2.5:0.01:4.5];
for i_T = 1:length(T)
f=(T(i_T)-24.5)*(T(i_T)+570.82);
ni = zeros(size(lamdas));
np = ni;
ns = ni;
nss = ni;
nii = ni;
for i_lamdas = 1:length(lamdas)
      lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));
      ni(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdai(i_lamdas).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdai(i_lamdas).^2-a5^2)-a6.*lamdai(i_lamdas).^2).^(1/2);
      np(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdap.^2-(a3+b3*f)^2)+(a4+b4.*f)./(lamdap.^2-a5^2)-a6.*lamdap.^2).^(1/2);
      ns(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdas(i_lamdas).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdas(i_lamdas).^2-a5^2)-a6.*lamdas(i_lamdas).^2).^(1/2);
end
for i_lamdass = 1:length(lamdass)
      lamdaii(i_lamdass)=1/(1/lamdas-1/lamdass(i_lamdass));
      nii(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdai2(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdai2(i_lamdass).^2-a5^2)-a6.*lamdai2(i_lamdass).^2).^(1/2);
      ns(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdas(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdas(i_lamdass).^2-a5^2)-a6.*lamdas(i_lamdass).^2).^(1/2);
      nss(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdass(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdass(i_lamdass).^2-a5^2)-a6.*lamdass(i_lamdass).^2).^(1/2);
end
   temp1 = abs(np./lamdap-ns./lamdas-ni./lamdai-1/A);
   [x1,i_x1] = min(temp1);
   lamdai_opt(i_T) = lamdai(i_x1);
   lamdas_opt(i_T) = lamdas(i_x1);
   temp2 = abs(ns./lamdas-nss./lamdass-nii./lamdaii-1/A);
   [x2,i_x2] = min(temp2);
   lamdaii_opt(i_T) = lamdaii(i_x2);
   lamdass_opt(i_T) = lamdass(i_x2);
end
plot(T,lamdai_opt,'b.')
hold on
plot(T,lamdas_opt,'b.')
hold on
plot(T,lamdaii_opt,'b.')
hold on
plot(T,lamdass_opt,'b.')

級聯(lián)OPO+OPA溫度調(diào)諧曲線.jpg

回復(fù)此樓

» 猜你喜歡

2026年循環(huán)經(jīng)濟功能材料國際會議（ICFMCE 2026）已經(jīng)有0人回復(fù)
2026年第五屆電氣、電子與信息工程國際會議（ISEEIE 2026）已經(jīng)有0人回復(fù)
物理學(xué)I論文潤色/翻譯怎么收費? 已經(jīng)有269人回復(fù)
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑已經(jīng)有5人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑已經(jīng)有5人回復(fù)
求調(diào)劑已經(jīng)有0人回復(fù)
0702一志愿吉大B區(qū)求調(diào)劑有論文已經(jīng)有0人回復(fù)
請問還有沒有用Latex寫文章的小伙伴們？已經(jīng)有0人回復(fù)
光學(xué)工程學(xué)碩調(diào)劑信息已經(jīng)有29人回復(fù)
歡迎加入課題組已經(jīng)有0人回復(fù)
散金幣，求好運，祝面上順利！已經(jīng)有52人回復(fù)

1樓 2013-01-04 20:25:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★
感謝參與，應(yīng)助指數(shù) +1
惜夢尋草: 金幣+3, ★有幫助 2013-01-05 15:56:02

你這問題好多，比如
lamdai2是多少
lamdas和lamdass長度不一樣怎么一對一的計算
lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));到底是除元素還是除矩陣？我估計除元素

贊一下

回復(fù)此樓

showmethemoney

2樓2013-01-05 10:44:09

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

惜夢尋草

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時
蟲號: 2118505
注冊: 2012-11-10
性別: MM
專業(yè): 光學(xué)

引用回帖:

2樓: Originally posted by csgt0 at 2013-01-05 10:44:09
你這問題好多，比如
lamdai2是多少
lamdas和lamdass長度不一樣怎么一對一的計算
lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));到底是除元素還是除矩陣？我估計除元素

lamdai2就是lamdaii，我在編程序的時候疏忽了
lamdas和lamdass需要長度一樣才能算嗎？這個我不太清楚
lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));正如你想的，是要除元素。
我知道這個程序的錯誤很多，我不知道應(yīng)該怎么算這個，你能幫我寫一下能得到下面那張圖的程序嗎？謝謝啦！

贊一下

回復(fù)此樓

3樓2013-01-05 15:55:28

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

引用回帖:

3樓: Originally posted by 惜夢尋草 at 2013-01-05 15:55:28
lamdai2就是lamdaii，我在編程序的時候疏忽了
lamdas和lamdass需要長度一樣才能算嗎？這個我不太清楚
lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));正如你想的，是要除元素。
我知道這個程序的錯誤很多 ...

因為你的temp2 = abs(ns./lamdas-nss./lamdass-nii./lamdaii-1/A);
在一個公式里面必須保證lamdas,lamdass,lamdaii都是一樣長得

贊一下

回復(fù)此樓

showmethemoney

4樓2013-01-05 16:49:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

惜夢尋草

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時
蟲號: 2118505
注冊: 2012-11-10
性別: MM
專業(yè): 光學(xué)

引用回帖:

4樓: Originally posted by csgt0 at 2013-01-05 16:49:51
因為你的temp2 = abs(ns./lamdas-nss./lamdass-nii./lamdaii-1/A);
在一個公式里面必須保證lamdas,lamdass,lamdaii都是一樣長得...

我將lamdas范圍改成[1.2:0.01:3.2],lamdass范圍改成[2.5:0.01:4.5]了，按你的說法，公式temp2=abs(ns./lamdas-nss./lamdass-nii./lamdaii-1/A)長度應(yīng)該是一樣的了，可是還是出不來圖啊。
仍然是Error using ==> mrdivide
Matrix dimensions must agree.
你要是能幫我畫出來，我可以再加金幣的。拜托了！

贊一下

回復(fù)此樓

5樓2013-01-06 08:43:35

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

★ ★ ★ ★ ★ ★ ★
惜夢尋草: 金幣+7, ★★★很有幫助 2013-01-06 14:40:14

能算，就是不知道圖對不對，不對就看你的算法錯沒錯

CODE:

clear

 clc

 n=1;

 m=1;

 a1=5.756;

 a2=0.0983;

 a3=0.2020;

 a4=189.32;

 a5=12.52;

 a6=1.32e-2;

 b1=2.860e-6;

 b2=4.700e-8;

 b3=6.113e-8;

 b4=1.516e-4;

 A=31.5;

 lamdap=1.064;

 T = [25:0.1:250];

 lamdas = [1.2:0.01:3.2];

 lamdass = [2.5:0.01:4.5];

 for i_T = 1:length(T)

     f=(T(i_T)-24.5)*(T(i_T)+570.82);

     ni = zeros(size(lamdas));

     np = ni;

     ns = ni;

     nss = ni;

     nii = ni;

     for i_lamdas = 1:length(lamdas)

         lamdai(i_lamdas)=1./(1./lamdap-1./lamdas(i_lamdas));

         ni(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdai(i_lamdas).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdai(i_lamdas).^2-a5^2)-a6.*lamdai(i_lamdas).^2).^(1/2);

         np(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdap.^2-(a3+b3*f)^2)+(a4+b4.*f)./(lamdap.^2-a5^2)-a6.*lamdap.^2).^(1/2);

         ns(i_lamdas)=(a1+b1.*f+(a2+b2.*f)./(lamdas(i_lamdas).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdas(i_lamdas).^2-a5^2)-a6.*lamdas(i_lamdas).^2).^(1/2);

     end

     for i_lamdass = 1:length(lamdass)

         lamdaii(i_lamdass)=1/(1/lamdas(i_lamdass)-1/lamdass(i_lamdass));

         nii(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdaii(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdaii(i_lamdass).^2-a5^2)-a6.*lamdaii(i_lamdass).^2).^(1/2);

         ns(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdas(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdas(i_lamdass).^2-a5^2)-a6.*lamdas(i_lamdass).^2).^(1/2);

         nss(i_lamdass)=(a1+b1.*f+(a2+b2.*f)./(lamdass(i_lamdass).^2-(a3+b3.*f)^2)+(a4+b4.*f)./(lamdass(i_lamdass).^2-a5^2)-a6.*lamdass(i_lamdass).^2).^(1/2);

     end

        temp1 = abs(np./lamdap-ns./lamdas-ni./lamdai-1/A); 

       [x1,i_x1] = min(temp1);

        lamdai_opt(i_T) = lamdai(i_x1);

        lamdas_opt(i_T) = lamdas(i_x1);

        temp2 = abs(ns./lamdas-nss./lamdass-nii./lamdaii-1/A);   

        [x2,i_x2] = min(temp2);

        lamdaii_opt(i_T) = lamdaii(i_x2);

        lamdass_opt(i_T) = lamdass(i_x2);     

 end

 plot(T,lamdai_opt,'b',T,lamdas_opt,'r',T,lamdaii_opt,'g',T,lamdass_opt,'b')

grid

贊一下

回復(fù)此樓

showmethemoney

6樓2013-01-06 11:12:53

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

惜夢尋草

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時
蟲號: 2118505
注冊: 2012-11-10
性別: MM
專業(yè): 光學(xué)

引用回帖:

6樓: Originally posted by csgt0 at 2013-01-06 11:12:53
能算，就是不知道圖對不對，不對就看你的算法錯沒錯

clear
clc
n=1;
m=1;
a1=5.756;
a2=0.0983;
a3=0.2020;
a4=189.32;
a5=12.52;
a6=1.32e-2;
b1=2.860e-6;
b2=4.700e-8;
b3=6.113e-8; ...

非常感謝！
不過我還有兩個問題:
第一，為什么我的程序里lamdas = [1.2:0.01:3.2];lamdass = [2.5:0.01:4.5];修改以后還是出不了圖呢？
第二，為什么你寫的程序理論上沒有問題，但是輸出的圖形和我要的不一樣呢？

贊一下

回復(fù)此樓

7樓2013-01-06 14:39:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

csgt0

榮譽版主 (著名寫手)

彩色掛圖

專家經(jīng)驗: +2
計算強帖: 2
應(yīng)助: 367 (碩士)
貴賓: 0.457
金幣: 7807
散金: 8897
紅花: 20
沙發(fā): 7
帖子: 2809
在線: 515.8小時
蟲號: 449046
注冊: 2007-11-02
專業(yè): 金屬非晶態(tài)、準(zhǔn)晶和納米晶
管轄: 計算模擬

【答案】應(yīng)助回帖

引用回帖:

7樓: Originally posted by 惜夢尋草 at 2013-01-06 14:39:20
非常感謝！
不過我還有兩個問題:
第一，為什么我的程序里lamdas = ;lamdass = ;修改以后還是出不了圖呢？
第二，為什么你寫的程序理論上沒有問題，但是輸出的圖形和我要的不一樣呢？...

可能有的地方你沒有改正。
至于圖形不一樣那就看你的算法對不對，你可以檢查計算方法和過程。

贊一下

回復(fù)此樓

showmethemoney

8樓2013-01-06 14:41:30

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

惜夢尋草

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時
蟲號: 2118505
注冊: 2012-11-10
性別: MM
專業(yè): 光學(xué)

引用回帖:

8樓: Originally posted by csgt0 at 2013-01-06 14:41:30
可能有的地方你沒有改正。
至于圖形不一樣那就看你的算法對不對，你可以檢查計算方法和過程。...

過程應(yīng)該不會有錯的，還是沒找出來原因。

贊一下

回復(fù)此樓

9樓2013-01-07 08:19:34

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

惜夢尋草

新蟲 (初入文壇)

應(yīng)助: 0 (幼兒園)
金幣: 875.5
帖子: 25
在線: 5.7小時
蟲號: 2118505
注冊: 2012-11-10
性別: MM
專業(yè): 光學(xué)

引用回帖:

8樓: Originally posted by csgt0 at 2013-01-06 14:41:30
可能有的地方你沒有改正。
至于圖形不一樣那就看你的算法對不對，你可以檢查計算方法和過程。...

我的問題已經(jīng)解決了，不用費心了~謝謝啦！

贊一下

回復(fù)此樓

10樓2013-01-07 10:50:44

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主惜夢尋草的主題更新

返回列表

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 化工學(xué)碩306求調(diào)劑 +10	42838695 2026-03-12	10/500	2026-03-18 14:42 by haxia
[考研] 298-一志愿中國農(nóng)業(yè)大學(xué)-求調(diào)劑 +7	手機用戶 2026-03-17	7/350	2026-03-18 14:34 by vgtyfty
[考研] 一志愿西南交大，求調(diào)劑 +4	材化逐夢人 2026-03-18	4/200	2026-03-18 14:22 by 007_lilei
[考研] 288求調(diào)劑，一志愿華南理工大學(xué)071005 +4	ioodiiij 2026-03-17	4/200	2026-03-18 12:36 by Linda Hu
[考研] 0703化學(xué)求調(diào)劑總分331 +3	ZY-05 2026-03-13	3/150	2026-03-18 10:58 by macy2011
[考研] 考研求調(diào)劑 +3	橘頌. 2026-03-17	4/200	2026-03-17 21:43 by 有只貍奴
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級食品/生物方向碩士 +3	1喜春8 2026-03-17	5/250	2026-03-17 17:18 by ber川cool子
[考研] 290求調(diào)劑 +3	p asserby. 2026-03-15	4/200	2026-03-17 16:35 by wangkm
[考研] 274求調(diào)劑 +5	時間點 2026-03-13	5/250	2026-03-17 07:34 by 熱情沙漠
[考研] 318求調(diào)劑 +3	Yanyali 2026-03-15	3/150	2026-03-16 16:41 by houyaoxu
[考研] 0703一志愿211 285分求調(diào)劑 +5	ly3471z 2026-03-13	5/250	2026-03-16 16:16 by 哦哦123
[考研] 283求調(diào)劑 +10	小樓。 2026-03-12	14/700	2026-03-16 16:08 by 13811244083
[考研] 070303 總分349求調(diào)劑 +3	LJY9966 2026-03-15	5/250	2026-03-16 14:24 by xwxstudy
[考研] 求老師收留調(diào)劑 +4	jiang姜66 2026-03-14	5/250	2026-03-15 20:11 by Winj1e
[考研] 26考研一志愿中國石油大學(xué)(華東)305分求調(diào)劑 +3	嘉年新程 2026-03-15	3/150	2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3	麥田耕者 2026-03-14	3/150	2026-03-14 20:04 by 外星文明
[考研] 招收0805（材料）調(diào)劑 +3	18595523086 2026-03-13	3/150	2026-03-14 00:33 by 123%、
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3	土木小天才 2026-03-12	3/150	2026-03-13 15:01 by JourneyLucky
[考研] 328化工專碩求調(diào)劑 +4	。，。，。，。i 2026-03-12	4/200	2026-03-13 14:44 by JourneyLucky
[考研] 081200-11408-276學(xué)碩求調(diào)劑 +3	崔wj 2026-03-12	4/200	2026-03-12 19:33 by 求調(diào)劑zz