24小时热门版块排行榜

>论坛更新日志 (4617)
>虫友互识 (537)
>文献求助 (439)
>休闲灌水 (401)
>导师招生 (399)
>考博 (136)
>硕博家园 (133)
>博后之家 (104)
>基金申请 (101)
>招聘信息布告栏 (100)
>论文投稿 (95)
>考研 (88)
>教师之家 (65)
>找工作 (50)
>公派出国 (48)
>绿色求助(高悬赏) (46)

返回列表

wzycrfpxz

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1255.8
红花: 1
帖子: 15
在线: 21.6小时
虫号: 1947568
注册: 2012-08-20
专业: 系统科学与系统工程

[求助] 求纠错。。。

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 sigma=(E1^2-483.14*E1+95.096)/rad;
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

回复此楼