24小时热门版块排行榜

返回列表

LCP1212

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 104.2
帖子: 42
在线: 18.5小时
虫号: 3620048

[交流] 请大神帮忙看下这个程序

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)];

老是出错，清帮忙改一下，谢谢各位啦，初学者不大会

回复此楼

» 猜你喜欢

26申博推荐：南京航空航天大学国际前沿院光学方向招收博士生！已经有1人回复
如何从铁电相到顺电相。已经有1人回复
物理学I论文润色/翻译怎么收费? 已经有296人回复
各位大佬，求一份最新Wien2k版本作为学习用途已经有0人回复
南科大活性流体和软物质课题组诚招2027级博士、硕士研究生和博后已经有31人回复
瑞典林雪平大学博士后招聘｜PEC Water Splitting 方向已经有24人回复
PRB投稿，Acknowledgment sent to author状态十天了已经有2人回复
PVA溶解已经有1人回复
基金申请已经有40人回复

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到:

查看全部散金贴

1楼 2015-03-18 10:17:38

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

LCP1212

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 104.2
帖子: 42
在线: 18.5小时
虫号: 3620048

??? 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

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

相关版块跳转我要订阅楼主 LCP1212 的主题更新

返回列表

24小时热门版块排行榜

LCP1212

[交流] 请大神帮忙看下这个程序

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:

» 抢金币啦！回帖就可以得到: 查看全部散金贴

LCP1212

» 抢金币啦！回帖就可以得到:

查看全部散金贴