| 查看: 1155 | 回复: 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 %参数设置 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; %物理常数及中间过程参数计算 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); %抽运光功率设置 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)]; 老是出错，清帮忙改一下，谢谢各位啦，初学者不大会 |
回复此楼» 猜你喜欢
【征稿】Springer Nature热电专题 Discover Materials / Discover Chemistry
已经有47人回复
-
25届硕士申请，材料计算方向
已经有2人回复
-
物理学I论文润色/翻译怎么收费?
已经有231人回复
-
课题组招1-2名硕士
已经有2人回复
-
机时越来越贵之后，我把一部分 DFT 流程搬到了 AI 计算平台上
已经有0人回复
-
斯坦福大学与Elsevier世界前2%学者新排名查询
已经有1人回复
-
27申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！
已经有19人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
apl投稿状态
已经有2人回复
-
华南师范大学+双一流学科物理学+广东省杰青团队2027年博士招生+机器学习驱动材料设计
已经有1人回复
-
日内瓦国际发明展申报完整流程（科研人员通用指南）
已经有11人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 各位大神，帮忙看一下这个谱
已经有4人回复
- 请大家帮忙看下是哪种菌类？
已经有4人回复
- 求大神帮忙看下我这篇文章投到哪里比较好
已经有3人回复
- 求大神帮检查下一个简单的UDF
已经有13人回复
- 求大神跑下1stOpt，一个sharemodel，谢谢！
已经有7人回复
- 求大神帮助小弟构建一下这个菌株序列的进化树
已经有5人回复
- 求大神帮忙解题啊！！！
已经有4人回复
- 大神们帮忙看看这是什么意思啊
已经有4人回复
- 求大神来帮助我,,,,关于C语言的.懂得快进来,,,
已经有7人回复
- 哪位大神帮帮忙
已经有4人回复
- 很简单的Fortran赋值出错啦，求大神们帮忙找bug...
已经有11人回复
- 请各位大神帮忙看看！！！急急
已经有7人回复
- 求大神帮忙解答！
已经有8人回复
- 求大神帮忙看一下断口
已经有6人回复
- 急急急，哪位大神帮忙分析一下透镜图~~看不懂~汗一个先
已经有20人回复
» 抢金币啦！回帖就可以得到:
-
海南大学徐月山老师招收2027年申请考核制博士研究生1名（国家青年人才团队）
+1/84
-
海南大学徐月山老师招收2027年申请考核博士研究生1名（国家青年人才团队）
+1/83
-
中国科学院过程工程研究所 科研助理（实验） 2026年招聘
+1/49
-
中国科学院过程工程研究所 全国重点实验室课题组 招收联合培养硕士研究生1-2名
+1/49
-
东北林业大学徐加廷教授招申请考核博士生
+3/45
-
宁波大学招聘科研助理1名
+1/38
-
昆明医科大学“柔性生物电子与智能传感课题组”招聘科研助理
+1/37
-
中国石油大学（北京）国家级大人才团队博士招生2名：化学、材料、石油工程：油田化学
+2/35
-
投了一篇PRL
+1/8
-
海南大学王志同&吴道雄&杨金霖老师团队招收计算+实验联合培养博士生
+1/7
-
加州大学宋吉奎教授团队WGBS多组学揭示DNMT1 CXXC结构域保护CpG岛低甲基化的双重机制
+1/7
-
海南大学海洋与极地地质团队2027年博士研究生招生启事
+1/6
-
Nature Communications｜Celcomen：空间因果解耦如何预测组织扰动
+1/5
-
中山大学深圳校区博导招生：转化生长因子TGF-β功能性抗体开发 （更新啦）
+1/5
-
招若干有分子生物，细胞培养背景的科研人员(中山大学)
+1/5
-
招若干有分子生物，细胞培养背景的科研人员(中山大学)
+1/4
-
以色列理工学院朱成招收联合培养博士生-催化或材料模拟
+1/4
-
华中科技大学金健课题组招收电催化CO2还原博士生
+1/3
-
光学/光子学方向博后招聘-新加坡科技与设计大学SUTD
+1/3
-
重症肺炎克雷伯菌肺炎的单核细胞炎症与适应性免疫功能紊乱
+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
