matlab解微分方程组

答案:2 悬赏:50 手机版

解决时间 2021-03-03 02:27

提问者网友：嗝是迷路的屁
2021-03-02 22:22

[a,b]=dsolve('Dy=-2*2.825678*10^-7*x*Dx','22050000*sin(atan(Dy/Dx))/(7056000-7500*t)-1.624=-5.652*10^-7*Dx^2-6.652*10^-1442537.6*cos(atan(Dy/Dx))/(7056000-7500*t)','t')
无解，求大神给求

最佳答案

五星知识达人网友：woshuo
2021-03-02 23:06

>> s=dsolve('Dy=-2*2.825678*10^-7*x*Dx','22050000*sin(atan(Dy/Dx))/(7056000-7500*t)-1.624=-5.652*10^-7*Dx^2-6.652*10^-1442537.6*cos(atan(Dy/Dx))/(7056000-7500*t)','t')
Warning: Explicit solution could not be found.
> In dsolve at 197

s =

[ empty sym ]

>> s=dsolve('Dy=-2*2.825678*10^-7*x*Dx','22050000*sin(atan(Dy/Dx))/(7056000-7500*t)-1.624=-5.652*10^-7*Dx^2-6.652*10^-1442537.6*cos(atan(Dy/Dx))/(7056000-7500*t)')
Warning: Explicit solution could not be found.
> In dsolve at 197

s =

[ empty sym ]
是不是就是无解啊！

全部回答

1楼网友：北城痞子
2021-03-03 00:46

这个好像是常微分方程组把，先建立一个调用式子的函数： function f=hanshu(y,x) f=[a(x(3)-x(2)); b(x(1)-0.5*sqrt(x(2)-c)); d(2.5-x(3))]; 然后： [x,y]=ode45('hanshu',[t_start,t_final],x0); plot(t,x); 把初值和区间写上。感觉你的x1应该是x(1)把。我也正好做长微分方程，有空加我qq724918113

我要举报

如以上回答内容为低俗、色情、不良、暴力、侵权、涉及违法等信息，可以点下面链接进行举报！

点此我要举报以上问答信息