布鲁塞尔振子控制问题（公式推导）请指教

notemple

布鲁塞尔振子：
dx/dt=A-(B+1)x+x^2*y+alpha*cos(w*t)
dy/dt=Bx-x^2*y

对该系统做比例微分控制，控制量反馈量:
deltay=k*dy/dt+y

则系统为：
dx/dt=A-(B+1)x+x^2*deltay+alpha*cos(w*t)
dy/dt=Bx-x^2*y
deltay=k*dy/dt+y

问：该式怎么化为dX＝f(X,t)的形式？若不能化为该形式，怎么用matlab中的ode45等函数画相图，poincare图等？

notemple

有人回答吗？

就是当一组方程组dX=f(X,t)，引入一个控制量deltay=k*dy/dt+y后如何写它的运动方程？必须要代入吗？
同理，若dX=f(X,t)中有很多的重复项，能否提取出来以减小计算量？

请指教，谢谢