作业帮matlab解非线性微分方程组小弟想求一个微分方程组的解析解,但是遇到了问题.模型：x,y,z为变量①dx/dt=ax②dy/dt=(bQx-S-T)*(1-z/Zm)*y③dz/dt=(cDx+P-M)*z+eyE我的程序：syms a b Q S T z Zm c A P M e E;SS=dsolve('Dx=a
来源：学生作业帮助网 编辑：作业帮 时间：2024/05/05 16:40:47
matlab解非线性微分方程组小弟想求一个微分方程组的解析解,但是遇到了问题.模型：x,y,z为变量①dx/dt=ax②dy/dt=(bQx-S-T)*(1-z/Zm)*y③dz/dt=(cDx+P-M)*z+eyE我的程序：syms a b Q S T z Zm c A P M e E;SS=dsolve('Dx=a
matlab解非线性微分方程组
小弟想求一个微分方程组的解析解,但是遇到了问题.
模型：x,y,z为变量
①dx/dt=ax
②dy/dt=(bQx-S-T)*(1-z/Zm)*y
③dz/dt=(cDx+P-M)*z+eyE
我的程序：
syms a b Q S T z Zm c A P M e E;
SS=dsolve('Dx=a*x','Dy=(b*Q*x-S-T)*y*(1-z/Zm)','Dz=(c*A*x+P-M)*z+e*E')；
结果：没有求到解析解,小弟是新手,不知道这个方程组是真的木有解,还是我的方法不对或有其他方法,
matlab解非线性微分方程组小弟想求一个微分方程组的解析解,但是遇到了问题.模型：x,y,z为变量①dx/dt=ax②dy/dt=(bQx-S-T)*(1-z/Zm)*y③dz/dt=(cDx+P-M)*z+eyE我的程序：syms a b Q S T z Zm c A P M e E;SS=dsolve('Dx=a
人家问的就是解析解,不知道上边几个在说什么
你这个方程高版本的matlab可以解,只是算出来之后带积分项目：
syms a b Q S T z Zm c A P M e E
SS=dsolve('Dx=a*x','Dy=(b*Q*x-S-T)*y*(1-z/Zm)','Dz=(c*A*x+P-M)*z+e*E')
SS =
y: [1x1 sym]
x: [1x1 sym]
z: [1x1 sym]
>> SS.x
ans =
C5*exp(a*t)
>> SS.y
ans =
C9*exp(int((C7*exp((P*a*t - M*a*t + A*C5*c*exp(a*t))/a)*(S + T - C5*Q*b*exp(a*t)))/Zm + (exp((P*a*t - M*a*t + A*C5*c*exp(a*t))/a)*int(e*exp(M*t - P*t - (A*C5*c*exp(a*t))/a + 1), t, IgnoreAnalyticConstraints)*(S + T - C5*Q*b*exp(a*t)))/Zm, t, IgnoreAnalyticConstraints) - T*t - S*t + (C5*Q*b*exp(a*t))/a)
>> SS.z
ans =
C7*exp((P*a*t - M*a*t + A*C5*c*exp(a*t))/a) + exp((P*a*t - M*a*t + A*C5*c*exp(a*t))/a)*int(e*exp(M*t - P*t - (A*C5*c*exp(a*t))/a + 1), t, IgnoreAnalyticConstraints)
IgnoreAnalyticConstraints是忽略分析限制,这个一般来说是可以忽略的