| 查看: 357 | 回复: 1 | ||
wzycrfpxz新虫 (初入文坛)
|
[求助]
求纠错。。。
|
 function guiji clc;clear all; format long; global CL CD s m hs lo0 r0 mu E0 Ef rad CL=0.3515;CD=1.4595;s=15.9;m=2804;hs=9354;lo0=0.0145;r0=3395.76*10^3;mu=4.284e+013;E0=1.5600418*10^7;... Ef=1.0766379*10^5;rad=57.325; solinit=bvpinit(linspace(0,1,100),[1 1 1 1 1]); sol=bvp4c(@f,@bcfun,solinit); plot(sol.x(1,  ,sol.x(2,);grid; end function dx=f(E1,x) global CL CD s m hs lo0 r0 mu E0 Ef rad L=lo0*(2.7183^(-(x(3)-r0)/hs))*s*CL*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0))/m; D=lo0*(2.7183^(-(x(3)-r0)/hs))*s*CD*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0))/m; g=mu/x(3)^2; n=length(E1); if E1>=0&&E1<=0.24 sigma=-10/rad; elseif 0.24 elseif E1>0.4&&E1<=0.61 sigma=-80/rad; elseif E1>0.61&&E1<=0.8 sigma=(E1^2+840.69*E1-593.193)/rad; elseif E1>0.8&&E1<=0.96 sigma=80/rad; elseif E1>0.96&&E1<0.97 sigma=23.9733/rad; elseif E1>0.97&&E1<=1 sigma=-23.9733/rad; end %dx=zeros(5,10); dx=[-(Ef-E0)/x(3)*cos(x(5)/rad)*cos(x(4)/rad)/cos(x(2)/rad)/D; -(Ef-E0)/x(3)*cos(x(5)/rad)*sin(x(4)/rad)/D; -(Ef-E0)*sin(x(5)/rad)/D; (Ef-E0)*[L*sin(sigma)/D/cos(x(5)/rad)/(2*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0)))... +cos(x(5)/rad)*cos(x(4)/rad)*tan(x(2)/rad)/x(3)/D]; -(Ef-E0)/(2*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0)))*[L*cos(sigma)/D-(g... -(2*(E1*(Ef-E0)+E0+mu*(1/x(3)-1/r0)))/x(3))*cos(x(5)/rad)/D]]; end function res=bcfun(xa,xb) global CL CD s m hs lo0 r0 mu E0 Ef rad res=[xa(1)+90.072;xa(2)+43.898; xa(3)-3520.76*10^3;xa(4)-7.99; xa(5)+14.15;xb(1)+73.26; xb(2)+41.45;xb(3)-3403.76*10^3]; endSample Text |
回复此楼» 猜你喜欢
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有3人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有3人回复
-
基金中了
已经有14人回复
-
关于filecode
已经有7人回复
-
娱乐
已经有6人回复
-
化学口download_prp&fileCode的固定段好像这几天一直没变，有变的大神么？
已经有4人回复
-
关于代码变化问题，想知道的进来
已经有12人回复
-
固定端突然变了，今天
已经有10人回复
-
听说今天filecode变了
已经有47人回复
tian187725
银虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 363.5
- 散金: 20
- 红花: 3
- 帖子: 319
- 在线: 16小时
- 虫号: 2005035
- 注册: 2012-09-16
- 性别: MM
- 专业: 药物代谢与药物动力学
2楼2012-11-12 17:50:42
20