24小时热门版块排行榜

返回列表

mlxqz

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 24
帖子: 6
在线: 17.9小时
虫号: 7600233
注册: 2017-12-10
专业: 机械结构强度学

[交流] 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解出来的总感觉不对

回复此楼

» 猜你喜欢

1楼2018-12-24 22:34:00

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

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

返回列表