matlab求全相应

求全相应

%求全相应
clear
syms s t Fs
A = [0 1; -2 -3];
B = [2; 0];
X0 = [0;1];
%输入
u = exp(-t);

%预解矩阵Fs
Fs = inv(s*eye(2)-A);
%预解矩阵拉普拉斯反变换得到矩阵指数函数
eAt = ilaplace(Fs,s,t);
%计算零输入响应
X0u = eAt * X0;

%求零状态
syms tao F
F = simple(subs(eAt,t,t-tao) * B * subs(u,t,tao));
X0s = int(F, tao,0,t);
%求全响应
Xt = X0u + X0s;
Xt = expand(Xt)

计算结果如下:

因为解题过程中用到化简公式,下面整理了一些常用的