 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 973 | 回復(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); %抽運光功率設(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)'); %端面抽運的光纖激光器速率方程組 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); %端面抽運的光纖激光器邊界條件 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)]; 老是出錯，清幫忙改一下，謝謝各位啦，初學者不大會 |
回復(fù)此樓» 猜你喜歡
哈爾濱理工大學物理系招收物理學考研調(diào)劑
已經(jīng)有7人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑
已經(jīng)有5人回復(fù)
-
物理學I論文潤色/翻譯怎么收費?
已經(jīng)有286人回復(fù)
-
求調(diào)劑
已經(jīng)有0人回復(fù)
-
0702一志愿吉大B區(qū)求調(diào)劑有論文
已經(jīng)有0人回復(fù)
-
請問還有沒有用Latex寫文章的小伙伴們？
已經(jīng)有0人回復(fù)
-
光學工程學碩調(diào)劑信息
已經(jīng)有31人回復(fù)
-
歡迎加入課題組
已經(jīng)有0人回復(fù)
-
散金幣，求好運，祝面上順利！
已經(jīng)有72人回復(fù)
-
華東師范大學芯片設(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īng)有5人回復(fù)
- 求大神幫忙解題�。。�！
已經(jīng)有4人回復(fù)
- 大神們幫忙看看這是什么意思啊
已經(jīng)有4人回復(fù)
- 求大神來幫助我,,,,關(guān)于C語言的.懂得快進來,,,
已經(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ù)
» 搶金幣啦！回帖就可以得到:
-
福建師范大學葉應(yīng)祥課題組招收2026級博士研究生
+2/140
-
中山大學環(huán)境材料團隊招聘
+1/85
-
材料、化學類碩士研究生調(diào)劑機會，學碩專碩均可，有意請聯(lián)系，非誠勿擾！
+1/84
-
雙一流南京醫(yī)科大學招計算機、AI、統(tǒng)計、生物信息等方向26年9月入學博士
+1/80
-
有沒有在香港的蟲友啊
+1/72
-
光學工程學碩調(diào)劑信息
+2/42
-
招收2026年材料與化工專碩調(diào)劑，山東省屬重點大學，過國家線即可調(diào)劑，名額充足
+2/38
-
天津城建大學招收物理、化學或者生物調(diào)劑生一名，前往天津大學聯(lián)合培養(yǎng)
+1/35
-
福建師范大學閩江學者陳巍教授團隊招收2026年化學、化工或材料專業(yè)的學碩6名
+1/31
-
下周見一個合作者，應(yīng)該問些什么
+1/19
-
電子信息求調(diào)劑
+1/17
-
湖南大學微生物結(jié)構(gòu)與功能實驗室2026年計劃招收博士研究生
+1/14
-
【2026考研調(diào)劑】南方科技大學-深圳理工大學聯(lián)培 086000 免疫方向(26.9入學)
+1/9
-
寧波東方理工大學余鵬課題組2026年上海交大聯(lián)培博士生、博士后、科研助理招聘
+1/9
-
北京航空航天大學招生微流控與智能打印技術(shù)方向博士研究生
+2/8
-
多倫多大學 & Synetic Labs 招博士后/研究員
+1/8
-
天津工業(yè)大學雙一流高校招博士、碩士、博士后
+1/6
-
211/雙一流石河子大學化學化工學院--電催化，CO2轉(zhuǎn)化
+1/3
-
2026考研的同學，分數(shù)不夠？想調(diào)劑上岸的同學，來看
+1/2
-
招收調(diào)劑生 國家重點實驗室（廣西某師范大學、B區(qū)）
+1/2
|
??? 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ā)表 | |
|---|---|---|---|---|
|
[考研] 287求調(diào)劑 +8 | 晨昏線與星海 2026-03-19 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 289求調(diào)劑 +7 | 懷瑾握瑜l 2026-03-20 | 7/350 |
|
|
|
[考研] 305分求調(diào)劑（食品工程） +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 求調(diào)劑院校信息 +6 | CX 330 2026-03-21 | 6/300 |
|
|
|
[考研] 資源與環(huán)境 調(diào)劑申請(333分) +5 | holy J 2026-03-21 | 5/250 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | 妮妮ninicgb 2026-03-21 | 4/200 |
|
|
|
[考研] 工科0856求調(diào)劑 +3 | 沐析汀汀 2026-03-21 | 3/150 |
|
|
|
[考研] 336求調(diào)劑 +5 | rmc8866 2026-03-21 | 5/250 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
|
[考研] 一志愿山大07化學 332分 四六級已過 本科山東雙非 求調(diào)劑！ +3 | 不想理你 2026-03-16 | 3/150 |
|
|
|
[考研] 303求調(diào)劑 +5 | 睿08 2026-03-17 | 7/350 |
|
|
|
[考研] 304求調(diào)劑 +7 | 司空. 2026-03-18 | 7/350 |
|
|
|
[考研] 350求調(diào)劑 +5 | weudhdk 2026-03-19 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +3 | @taotao 2026-03-20 | 3/150 |
|
|
|
[考研] 環(huán)境工程調(diào)劑 +9 | 大可digkids 2026-03-16 | 9/450 |
|
|
|
[考研] 工科材料085601 279求調(diào)劑 +7 | 困于星晨 2026-03-17 | 9/450 |
|
|
|
[考研] 一志愿中國海洋大學，生物學，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 0703化學調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
|
[考研] 有沒有道鐵/土木的想調(diào)劑南林，給自己招師弟中～ +3 | TqlXswl 2026-03-16 | 7/350 |
|