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/06/27 22:36:42
![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](/uploads/image/z/1157040-0-0.jpg?t=matlab%E8%A7%A3%E9%9D%9E%E7%BA%BF%E6%80%A7%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E5%B0%8F%E5%BC%9F%E6%83%B3%E6%B1%82%E4%B8%80%E4%B8%AA%E5%BE%AE%E5%88%86%E6%96%B9%E7%A8%8B%E7%BB%84%E7%9A%84%E8%A7%A3%E6%9E%90%E8%A7%A3%2C%E4%BD%86%E6%98%AF%E9%81%87%E5%88%B0%E4%BA%86%E9%97%AE%E9%A2%98.%E6%A8%A1%E5%9E%8B%EF%BC%9Ax%2Cy%2Cz%E4%B8%BA%E5%8F%98%E9%87%8F%E2%91%A0dx%2Fdt%3Dax%E2%91%A1dy%2Fdt%3D%28bQx-S-T%29%2A%281-z%2FZm%29%2Ay%E2%91%A2dz%2Fdt%3D%28cDx%2BP-M%29%2Az%2BeyE%E6%88%91%E7%9A%84%E7%A8%8B%E5%BA%8F%EF%BC%9Asyms+a+b+Q+S+T+z+Zm+c+A+P+M+e+E%3BSS%3Ddsolve%28%27Dx%3Da)
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是忽略分析限制,这个一般来说是可以忽略的