[求助]一个发酵模型的拟合与优化

[dur08]

需要根据测得的数据建模进行优化，找到了一个例子，但是里面的模型那块怎么改也不对，求达人帮助，感激感激！
function xdot=moxing(t,x)
global A B C E

if (A*x(2)*x(1)/(B+x(2)))>0
xdot=[A*x(2)/(B+x(2)),0,0
-C*A*x(2)/(B+x(2)),0,0
E*x(1),0,0]*x;
else xdot=[A*x(2)/(B+x(2)),0,0
0,0,0
E*x(1),0,0]*x;
end
例子中的模型是：
dx(1)/dt=A*x(2)*x(1)/(B+x(2))
dx(2)/dt=C*A*x(2)*x(1)/(B+x(2))
dx(3)/dt=E*x(1)^2
我的模型是：
dx(1)/dt=A*x(2)*x(1)/(B+x(2))
dx(2)/dt=C*A*x(2)*x(1)/(B+x(2))
dx(3)/dt=E*x(1)+D*A*x(2)*x(1)/(B+x(2))
急盼哪位高人指点一下，或者解释一下怎样把模型中的函数转换成xdot=[],万分感谢！！！