syms s t
A = [5 0 3;0 5 2;9 1 3];
Fs = inv(seye(3)-A); %求预解矩阵FS=〖(sI-A)〗^(-1),eye(3)为33矩阵
eAt = ilaplace(Fs,s,t); %求ⅇAt=L(-1) [〖(sI-A)〗^(-1) ]
eAt = simplify(eAt); %化简ⅇ^At的表达式
解得
eAt =
[ exp(4t), 0, 0]
[ 0, (exp(2t)(exp(2t) + 1))/2, (exp(2t)(exp(2t) - 1))/2]
[ 0, (exp(2t)(exp(2t) - 1))/2, (exp(2t)(exp(2*t) + 1))/2]
应用matlab计算线性定常系统的矩阵指数
最新推荐文章于 2023-12-11 17:29:05 发布