矩阵的分析与处理
1 矩阵的行列式、秩与迹及特征值分析
- 矩阵的行列式
运算符:det( )
,用于求方阵的行列式值 - 矩阵的秩与迹
运算符:rank( )
,矩阵的秩是矩阵的列向量组(或行向量组)的任一极大线性无关组所含向量的个数
运算符:trace( )
,矩阵的迹等于矩阵主对角线元素的总和。也等于矩阵特征值的总和。(要求矩阵是方阵) - 矩阵的特征值分析
运算符:E=eig(A)
,求矩阵A的全部特征值,并构成向量E
运算符:[V,D]=eig(A)
,求矩阵A的全部特征值,构成对角矩阵D;求A的特征向量构成列向量V。
2 矩阵的逆与线性方程组的求解
- 矩阵的逆
inv( )
:用于求满秩方阵的逆
pinv( )
:用于求不是方阵或非满秩方阵的逆——伪逆(如果ABA=A,BAB=B,则称B为A的伪逆,或广义逆矩阵)。 - 线性方程组求解
x1- 2x2+3x3=1
2x1+3x2+x3=2
3x1- x2- x3=4
A=[1,-2, 3;2,3,1;3,-1 ,-1];
B=[1;2;4];
X=inv(A)*B
X1 = A\B %%矩阵左除
3 矩阵的分解与变换
- 矩阵的分解
- 三角分解(方阵):
[l,u]=lu(a)
- 正交分解:
[q,r]=qr(a)
。
a(n,m)
q:n阶正交方阵
r:与a同阶的上三角矩阵 - 奇异值分解:
[u,s,v]=svd(a)
u: n阶正交方阵
s: n×m阶的对角阵,对角线元素为a的奇异值,长度为n、m的较小者
v:m阶正交方阵
- 矩阵的变换
- 矩阵的共轭转置:
'
- 矩阵的共轭:
conj
- 矩阵的转置:
conj'
- 复数矩阵的赋值
对元素逐个赋值:z=[1+2i,3+4i;5+6i,7+8i]
对实部和虚部矩阵分别赋值:z=[1,3;5,7]+[2,4;6,8]*i
注意:只有数字和i的乘积可省略乘号。 - 矩阵的行、列扩展
c a=[1,2,3; 4,5,6;7,8,9]
行扩展:a(4,3)=6.5
;a(5,:)=[5,4,3]
列扩展:a(:,4)=[5;4;3;2;1]