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

登錄注冊

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進行實名認證將不得使用互聯(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)有57人回復(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)、準晶和納米晶
管轄: 計算模擬

【答案】應(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)、準晶和納米晶
管轄: 計算模擬

【答案】應(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)、準晶和納米晶
管轄: 計算模擬

【答案】應(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)、準晶和納米晶
管轄: 計算模擬

【答案】應(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é)，生物學(xué)，301分，求調(diào)劑 +4	1孫悟空 2026-03-17	4/200	2026-03-18 17:59 by fivewind
[考研] 一志愿武理材料305分求調(diào)劑 +5	想上岸的鯉魚 2026-03-18	6/300	2026-03-18 17:53 by 無際的草原
[考研] 材料與化工一志愿南昌大學(xué)327求調(diào)劑推薦 +8	Ncdx123456 2026-03-13	9/450	2026-03-18 14:40 by haxia
[考研] 312求調(diào)劑 +8	陌宸希 2026-03-16	9/450	2026-03-18 12:39 by Linda Hu
[考研] 0703化學(xué)調(diào)劑 +3	妮妮ninicgb 2026-03-17	3/150	2026-03-18 10:29 by macy2011
[考研] 265求調(diào)劑 +3	梁梁校校 2026-03-17	3/150	2026-03-18 09:12 by zhukairuo
[考研] 334求調(diào)劑 +3	志存高遠意在機?/a> 2026-03-16	3/150	2026-03-18 08:34 by lm4875102
[考研] 085601求調(diào)劑 +4	Du.11 2026-03-16	4/200	2026-03-17 17:08 by ruiyingmiao
[考研] 332求調(diào)劑 +6	Zz版 2026-03-13	6/300	2026-03-17 17:03 by ruiyingmiao
[考研] 283求調(diào)劑 +3	聽風(fēng)就是雨； 2026-03-16	3/150	2026-03-17 07:41 by 熱情沙漠
[考研] 0854控制工程 359求調(diào)劑可跨專業(yè) +3	626776879 2026-03-14	9/450	2026-03-16 17:42 by 626776879
[考研] 326求調(diào)劑 +3	mlpqaz03 2026-03-15	3/150	2026-03-16 07:33 by Iveryant
[考研] 266求調(diào)劑 +4	學(xué)員97LZgn 2026-03-13	4/200	2026-03-14 08:37 by zhukairuo
[考研] 26調(diào)劑/材料科學(xué)與工程/總分295/求收留 +9	2026調(diào)劑俠 2026-03-12	9/450	2026-03-13 20:46 by 18595523086
[考研] 工科278分求調(diào)劑 +5	周慢熱啊 2026-03-12	7/350	2026-03-13 15:49 by JourneyLucky
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3	T123 tt 2026-03-12	3/150	2026-03-13 10:49 by houyaoxu
[考研] 070303一志愿西北大學(xué)學(xué)碩310找調(diào)劑 +3	d如愿上岸 2026-03-13	3/150	2026-03-13 10:43 by houyaoxu
[考博] 福州大學(xué)楊黃浩課題組招收2026年專業(yè)學(xué)位博士研究生，2026.03.20截止 +3	Xiangyu_ou 2026-03-12	3/150	2026-03-13 09:36 by duanwu655
[考研] 283求調(diào)劑，材料、化工皆可 +8	蘇打水7777 2026-03-11	10/500	2026-03-13 09:06 by Linda Hu
[考研] 321求調(diào)劑（食品/專碩） +3	xc321 2026-03-12	6/300	2026-03-13 08:45 by xc321