MATLAB求解常系数时变系统微分方程组dxm(t)=A*xm(t)+B*um(t)dxp(t)=A*xp(t)+B*[up(t)+r(t)]dr(t)=C*[um(t)-up(t)]这个方程组怎么求解啊?

来源:学生作业帮助网 编辑:作业帮 时间:2024/04/29 15:50:17
MATLAB求解常系数时变系统微分方程组dxm(t)=A*xm(t)+B*um(t)dxp(t)=A*xp(t)+B*[up(t)+r(t)]dr(t)=C*[um(t)-up(t)]这个方程组怎么求解啊?

MATLAB求解常系数时变系统微分方程组dxm(t)=A*xm(t)+B*um(t)dxp(t)=A*xp(t)+B*[up(t)+r(t)]dr(t)=C*[um(t)-up(t)]这个方程组怎么求解啊?
MATLAB求解常系数时变系统微分方程组
dxm(t)=A*xm(t)+B*um(t)
dxp(t)=A*xp(t)+B*[up(t)+r(t)]
dr(t)=C*[um(t)-up(t)]
这个方程组怎么求解啊?

MATLAB求解常系数时变系统微分方程组dxm(t)=A*xm(t)+B*um(t)dxp(t)=A*xp(t)+B*[up(t)+r(t)]dr(t)=C*[um(t)-up(t)]这个方程组怎么求解啊?
这个题目并不难,只是你的题目写太晦涩难懂.
我也是看了好长时间才猜出来的,还不知道对不对呢.
题意:xm(t),xp(t),r(t)是要求解的未知函数,而um(t),up(t)是已知的关于t的函数.程序如下:
clear;clc;
syms A B C t;
[xm,xp,r]=dsolve('Dxm=A*xm+B*um','Dxp=A*xp+B*(up+r)','Dr=C*(um-up)','t')
结果:
xm =
C*t*um-C*t*up+C2
xp =
(-B*um+exp(A*t)*C3*A)/A
r =
-(B*up*A+B*C*um*t*A+C*B*um-B*C*up*t*A-C*B*up+C2*B*A-exp(A*t)*C1*A^2)/A^2
>>