灵鹊做题网matlab常微分方程组的解法
来源：学生作业帮 编辑：灵鹊做题网作业帮 分类：综合作业 时间：2024/05/01 03:25:12
matlab常微分方程组的解法
matlab怎么求一阶常微分方程组（solve做不出来的那种）
matlab怎么求一阶常微分方程组（solve做不出来的那种）
s=dsolve(‘方程1’, ‘方程2’,…,’初始条件1’,’初始条件2’ …,’自变量’)
用字符串方程表示,自变量缺省值为t.导数用D表示,2阶导数用D2表示,以此类推.S返回解析解.在方程组情形,s为一个符号结构.
[tout,yout]=ode45(‘yprime’,[t0,tf],y0) 采用变步长四阶Runge-Kutta法和五阶Runge-Kutta-Felhberg法求数值解,yprime是用以表示f(t,y)的M文件名,t0表示自变量的初始值,tf表示自变量的终值,y0表示初始向量值.输出向量tout表示节点(t0,t1, …,tn)T,输出矩阵yout表示数值解,每一列对应y的一个分量.若无输出参数,则自动作出图形
http://jpkc.yrcti.edu.cn/2006/gdsx/link/sxsj/sxsy6.htm
再问： 还是不会，Dy1=y2 Dy2=400*(1-y1^2)*y2-y1 (0《x《1000) y1(0)=0 y2(0)=1 怎么求出来并且画图
再答： odefun=@(t,x)[x(2);400*(1-x(1)^2)*x(2)-x(1)]; [t,x]=ode23(odefun,[0 1000],[0 1]); plot(t,x(:,1),'r',t,,x(:,2)'b') % x(:,1)----y1 x(:,2)--y2 t----x
再问： 我会了，非常感谢，能不能再问一下，这样求出来的是不是满足 3（是个符号对称过来）
用字符串方程表示,自变量缺省值为t.导数用D表示,2阶导数用D2表示,以此类推.S返回解析解.在方程组情形,s为一个符号结构.
[tout,yout]=ode45(‘yprime’,[t0,tf],y0) 采用变步长四阶Runge-Kutta法和五阶Runge-Kutta-Felhberg法求数值解,yprime是用以表示f(t,y)的M文件名,t0表示自变量的初始值,tf表示自变量的终值,y0表示初始向量值.输出向量tout表示节点(t0,t1, …,tn)T,输出矩阵yout表示数值解,每一列对应y的一个分量.若无输出参数,则自动作出图形
http://jpkc.yrcti.edu.cn/2006/gdsx/link/sxsj/sxsy6.htm
再问： 还是不会，Dy1=y2 Dy2=400*(1-y1^2)*y2-y1 (0《x《1000) y1(0)=0 y2(0)=1 怎么求出来并且画图
再答： odefun=@(t,x)[x(2);400*(1-x(1)^2)*x(2)-x(1)]; [t,x]=ode23(odefun,[0 1000],[0 1]); plot(t,x(:,1),'r',t,,x(:,2)'b') % x(:,1)----y1 x(:,2)--y2 t----x
再问： 我会了，非常感谢，能不能再问一下，这样求出来的是不是满足 3（是个符号对称过来）