24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1010  |  回复: 0

mlxqz

新虫 (初入文坛)

[交流] Matlab求解微分方程bvp4c

solinit = bvpinit(linspace(0,14,10),[1;0]);
sol=bvp4c(@ODEfun,@BCfun,solinit);
x = linspace(0,14,127);
y = deval(sol,x);
plot(x,y(2,)

function dydx=ODEfun(x,y)
a=pi*14.701^4/64;
b=pi*(24^4-14.701^4)/64;
c=3.54365802196686;
d=5.94058826800217;
e=2/14.701/pi;
f=14.701/2;
dydx(1,1)=y(2);
dydx(2,1)=(1/a+1/b)/(c+d)*f^2/e*(cos(x/(f*e)))^2*y(1)-(tan(x/(f*e)))/f/e*y(2);

function res=BCfun(ya,yb)
res=[ya(1); yb(1)-1];
哪位大神给看看我这个微分方程的数值解是多少。用bvp4c解出来的总感觉不对
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mlxqz 的主题更新
普通表情 高级回复(可上传附件)
信息提示
请填处理意见