solution.m

来自「关于病态系统的识别以及求解的方法。其中包括几种解法的求解精度和速度的比较。」· M 代码 · 共 21 行

%-----------------------------------
%求解解析解
%-----------------------------------
syms s ;            %定义字符变量
A=[-21,19,-20;
    19,-21,20;
    40,-40,-40];
root=eig(A)         %求矩阵A的特征根，并判断是否是病态系统
B=[s,0,0;
    0,s,0;
    0,0,s];
C=B-A;              %（SI-A）
D=inv(C);           %（SI-A）的逆
eAT=ilaplace(D);    %exp(AT)=L((SI-A)-1)
X=[0;0;0];
X= eAT * [1;0;-1]   %解析解
t=0:0.01:2;         %绘图采样点间隔和区间长度
ezplot(X(1),t)      
%ezplot(X(2),t)
%ezplot(X(3),t)