请教大侠怎么求解这个耦合方程?

julian

请教大侠怎么求解这个耦合方程?
a*y..+b*y.=c*q..
q..+d*q=e*y..

(注:..为y和q上的两点,即对时间的导数)

sogooda

做个变量代换，x(1)=y, x(2)=y., x(3)=q, x(4)=q.
之后用x的元素就可以表示出x.,即可转化为ode45或者类似的ode求解器可解的形式。