matlab求解微分方程组
发布网友
发布时间:2022-03-25 19:59
我来回答
共1个回答
热心网友
时间:2022-03-25 21:28
程序本身没什么问题,就是算的很慢而已。你想快一点可以用ode15s、ode23s等刚性求解器。
--------------
A=dsolve('Dy1/dt=P0-r6*y1/I2-y2/c7','Dy2/dt=y1/I2-y2/(r4*c7)+y3/(r4*c8)','Dy3/dt=y2/(r4*c7)-y3/(r4*c8)-A0*y4/Im','Dy4/dt=A0*y3/c8-Sf-r5*Pm/Im')
然后A.y1到A.y4就是解。解里含有rootof的形式,指的是里面的根,因为都是字母,所以写不出显式的根,只能这样表示。