 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 972 | 回復(fù): 1 | |||
[交流]
請大神幫忙看下這個程序
|
|
function Pout = fiberlaser_twoend global R1 R2 Ppl Ppr sigma_ap sigma_ep sigma_as sigma_es gamma_s ... gamma_p N alpha_p alpha_s Pssat Ppsat Ppsp Ppsv mu kappa elta %參數(shù)設(shè)置 lambda_s = 1090 * 1e-9; lambda_p = 975 * 1e-9; tau = 1e-3; sigma_ap = 25e-21*1e-4; sigma_ep = 25e-21*1e-4; sigma_as = 1.4e-23*1e-4; sigma_es = 2e-21*1e-4; A_c = 5e-11; N = 4e+025; alpha_p = 3e-5*1e2; alpha_s = 5e-5*1e2; gamma_s = 0.82; gamma_p = 0.0012; R1 =.98; R2 =.04; L = 50; %物理常數(shù)及中間過程參數(shù)計算 c = 3e8; h = 6.626e-34; nu_s = c/lambda_s; nu_p = c/lambda_p; Pssat = h * nu_s * A_c/( gamma_s * (sigma_es+sigma_as) * tau); Ppsat = h * nu_p * A_c/( gamma_p * (sigma_ep+sigma_ap) * tau); %抽運(yùn)光功率設(shè)置 Ppl =1:1:25; Ppr =0; Pout=zeros(1,25); for i=1:25 OPTION = bvpset('Stats','ON'); solinit = bvpinit(linspace(0,L,10),[Ppl(i) Ppr 10 Ppr]); sol = bvp4c(@f,@fsbc,solinit); x = [sol.x]; y = [sol.y]; Pout(i) = y(3,end)*(1-R2); end plot(Ppl,Pout,'r'); grid on; title('Pump and laser powers'); legend('Pp+(z)','Pp-(z)','Ps+(z)','Ps-(z)'); xlabel('Position z (m)'); ylabel('Power (W)'); %端面抽運(yùn)的光纖激光器速率方程組 function dy = f(x,y) global sigma_ap sigma_ep sigma_as sigma_es gamma_s gamma_p N ... alpha_p alpha_s Pssat Ppsat Ppsp Ppsv dy = zeros(4,1); N21=N*(sigma_ap/(sigma_ap+sigma_ep)*(y(1)+y(2))/Ppsat+sigma_as/... (sigma_as+sigma_es)*(y(3)+y(4))/Pssat)/((y(1)+y(2))/Ppsat+... 1+(y(4)+y(3))/Pssat); dy(1)=(-gamma_p*(sigma_ap*N-(sigma_ap+sigma_ep)*N21)-alpha_p)*y(1); dy(2)=-(-gamma_p*(sigma_ap*N-(sigma_ap+sigma_ep)*N21)-alpha_p)*y(2); dy(3)=(gamma_s*((sigma_as+sigma_es)*N21-sigma_as*N)-alpha_s)*y(3); dy(4)=-(gamma_s*((sigma_as+sigma_es)*N21-sigma_as*N)-alpha_s)*y(4); %端面抽運(yùn)的光纖激光器邊界條件 function res = fsbc(y0,yL) global R1 R2 Ppl Ppr res = [y0(1)-Ppl yL(2)-Ppr y0(3)-R1*y0(4) yL(4)-R2*yL(3)]; 老是出錯，清幫忙改一下，謝謝各位啦，初學(xué)者不大會 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學(xué)物理系招收物理學(xué)考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學(xué)I論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有291人回復(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)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運(yùn)，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學(xué)芯片設(shè)計徐瓏真實(shí)水平咋樣
已經(jīng)有1人回復(fù)
-
調(diào)劑
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 各位大神，幫忙看一下這個譜
已經(jīng)有4人回復(fù)
- 請大家?guī)兔聪率悄姆N菌類？
已經(jīng)有4人回復(fù)
- 求大神幫忙看下我這篇文章投到哪里比較好
已經(jīng)有3人回復(fù)
- 求大神幫檢查下一個簡單的UDF
已經(jīng)有13人回復(fù)
- 求大神跑下1stOpt，一個sharemodel，謝謝！
已經(jīng)有7人回復(fù)
- 求大神幫助小弟構(gòu)建一下這個菌株序列的進(jìn)化樹
已經(jīng)有5人回復(fù)
- 求大神幫忙解題�。。�！
已經(jīng)有4人回復(fù)
- 大神們幫忙看看這是什么意思啊
已經(jīng)有4人回復(fù)
- 求大神來幫助我,,,,關(guān)于C語言的.懂得快進(jìn)來,,,
已經(jīng)有7人回復(fù)
- 哪位大神幫幫忙
已經(jīng)有4人回復(fù)
- 很簡單的Fortran賦值出錯啦，求大神們幫忙找bug...
已經(jīng)有11人回復(fù)
- 請各位大神幫忙看看！��！急急
已經(jīng)有7人回復(fù)
- 求大神幫忙解答！
已經(jīng)有8人回復(fù)
- 求大神幫忙看一下斷口
已經(jīng)有6人回復(fù)
- 急急急，哪位大神幫忙分析一下透鏡圖~~看不懂~汗一個先
已經(jīng)有20人回復(fù)
» 搶金幣啦！回帖就可以得到:
-
上海真誠脫單，這是一篇很認(rèn)真的相親貼，找一個靈魂伴侶
+1/180
-
電力全國重點(diǎn)實(shí)驗(yàn)室雙一流A類長江學(xué)者團(tuán)隊(duì)招2026年全日制博士1名（4月3日截止）
+2/96
-
中飛院航空鋰電池團(tuán)隊(duì)接受碩士生調(diào)劑
+1/87
-
材料、化學(xué)類碩士研究生調(diào)劑機(jī)會，學(xué)碩專碩均可，有意請聯(lián)系，非誠勿擾！
+1/84
-
飄在美國，有點(diǎn)想家
+1/51
-
Alicat 層流壓差質(zhì)量流量控制器在新型硅碳負(fù)極流化床的應(yīng)用- 艾里卡特（Alicat）
+3/43
-
招收材料、化學(xué)類調(diào)劑生，專碩學(xué)碩均可
+1/40
-
哈爾濱工業(yè)大學(xué)（深圳）-何自開教授團(tuán)隊(duì)誠招化學(xué)工程與技術(shù)專業(yè)博士研究生
+1/39
-
湘潭大學(xué)“過程強(qiáng)化與綠色化工”創(chuàng)新團(tuán)隊(duì)補(bǔ)招2026年秋入學(xué)博士生
+1/35
-
中南大學(xué)喻桂朋教授團(tuán)隊(duì)新增2026年博士研究生
+1/16
-
招收專業(yè)代碼08的學(xué)碩！
+1/15
-
【深圳理工大學(xué)】院士工作站招聘博士后-電催化和有機(jī)電化學(xué)
+1/12
-
【2026考研調(diào)劑】南方科技大學(xué)-深圳理工大學(xué)聯(lián)培 086000 免疫方向(26.9入學(xué))
+1/9
-
中國民航飛行學(xué)院先進(jìn)連接團(tuán)隊(duì)招收航空宇航、航空機(jī)械專業(yè)碩士研究生（學(xué)碩專碩均可）
+1/8
-
江西贛州-贛南師范大學(xué)-學(xué)校推薦
+1/7
-
我寫了一篇關(guān)于“健康永生”的論文，被VITA秒拒過程公開了，歡迎討論
+1/4
-
招收化工與材料學(xué)科點(diǎn)研究生
+1/2
-
“機(jī)械、材料”招生
+1/1
-
中科院上海微系統(tǒng)所材料生長理論與仿真課題組-博士后招聘（材料計算與AI4S）
+1/1
-
溫州醫(yī)科大學(xué)招收特種醫(yī)學(xué)博士研究生1名（專業(yè)不限，有實(shí)驗(yàn)基礎(chǔ)）
+1/1
|
??? Error using ==> vertcat CAT arguments dimensions are not consistent. Error in ==> p300_fiberlaser_twoend>fsbc at 67 res = [y0(1)-Ppl Error in ==> bvparguments at 117 testBC = bc(ya,yb,bcExtras{:}); Error in ==> bvp4c at 130 [n,npar,nregions,atol,rtol,Nmax,xyVectorized,printstats] = ... Error in ==> p300_fiberlaser_twoend at 38 sol = bvp4c(@f,@fsbc,solinit); 老是顯示上面的錯誤 |
2樓2015-03-18 10:21:38
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 LCP1212 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
 一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 一志愿華中農(nóng)業(yè)071010，總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 269專碩求調(diào)劑 +6 | 金恩貝 2026-03-21 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)070300 學(xué)碩336求調(diào)劑 +5 | vv迷 2026-03-21 | 8/400 |
|
|
|
[考研] 286求調(diào)劑 +10 | Faune 2026-03-21 | 10/500 |
|
|
|
[考研] 初試 317 +7 | 半拉月丙 2026-03-20 | 7/350 |
|
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 0805材料320求調(diào)劑 +3 | 深海物語 2026-03-20 | 3/150 |
|
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
|
[考研] 一志愿南昌大學(xué)，327分，材料與化工085600 +9 | Ncdx123456 2026-03-19 | 9/450 |
|
|
|
[考研] 考研調(diào)劑求學(xué)校推薦 +3 | 伯樂29 2026-03-18 | 5/250 |
|
|
|
[考研] 材料學(xué)碩297已過四六級求調(diào)劑推薦 +11 | adaie 2026-03-19 | 11/550 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
|
[考研] 281求調(diào)劑（0805） +14 | 煙汐憶海 2026-03-16 | 25/1250 |
|
|
|
[考研] 0856調(diào)劑，是學(xué)校就去 +8 | sllhht 2026-03-19 | 9/450 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 考研求調(diào)劑 +3 | 橘頌. 2026-03-17 | 4/200 |
|
|
|
[論文投稿] 有沒有大佬發(fā)小論文能帶我個二作 +3 | 增銳漏人 2026-03-17 | 4/200 |
|