det(A)trace(A)rank(A)norm(A)poly(A)这是数值法求解,解析法可以用
charppoly,新版本方法可能有改变A=[16,2,3,13;5,11,10,8;9,7,6,12;4,14;15,1]; p1=poly(A); p2=charpoly(sym(A));
- 1
- 2
- 3
poly(a,A),a是多项式系数的行向量
[
a
1
,
a
2
,
⋯
,
a
n
+
1
]
[a_1,a_2,\cdots,a_{n+1}]
[a1,a2,⋯,an+1][V,D]=eig(A),D是特征值对角矩阵,V是D对应的特征向量矩阵。expm(A)矩阵指数 e A t e^{\mathbf{A}t} eAt
expm(A*t)
- 1
矩阵的其他函数如 cos ( A ) \cos(\mathbf{A}) cos(A)
funm(A,@cos)
- 1
注意
funm使用了特征值和特征向量的求解方法,如果矩阵有重根,特征向量矩阵可能是奇异矩阵,求解会失效,这时候应该用泰勒幂级数展开求解。
[L,U]=lu(A),L是上三角矩阵,U是下三角矩阵。[D,P]=chol(A),
A
=
D
T
D
\mathbf{A}=\mathbf{D}^T\mathbf{D}
A=DTDQ=orth(A)[L,B,M]=svd(A),
A
=
L
B
M
T
\mathbf{A}=\mathbf{LBM}^T
A=LBMTc=cond(A)C=inv(A)B=pinv(A)X=A\B,对于线性方程组
A
X
=
B
\mathbf{AX}=\mathbf{B}
AX=BX=lyap(A,C),对于
A
X
+
X
A
T
=
−
C
\mathbf{AX+XA}^T=\mathbf{-C}
AX+XAT=−C,
C
\mathbf{C}
C是对称矩阵X=dlyap(A,C),对于
X
A
X
T
−
X
+
C
=
0
\mathbf{XAX}^T-\mathbf{X}+\mathbf{C}=\mathbf{0}
XAXT−X+C=0X=lyap(A,B,C),对于
A
X
+
X
B
=
−
C
\mathbf{AX+XB}=\mathbf{-C}
AX+XB=−C,解析解可以用lyapsym函数.X=are(A,B,C),对于
A
T
X
+
X
A
−
X
B
X
+
C
=
0
\mathbf{A}^T\mathbf{X+XA-XBX+C}=\mathbf{0}
ATX+XA−XBX+C=0,离散系统用dare函数。