【MATLAB】线性代数方程组及矩阵特征值

      矩阵是由英国数学家凯莱于1855年作为一个独立的概念引入数学中的，在经济工作中应用广泛． 矩阵是线性代数的一个重要研究对象，在数学理论和社会生产中都有着广泛的应用，尤其是计算机技术的发展为矩阵应用开辟了广阔的前景。因此，学习和掌握矩阵的基本理论和计算方法是必不可少的。

一、对角阵与三角阵

1. 对角阵：求矩阵的对角元素用diag(A)函数，其还有更进一步的形式diag(A,k)，其功能是提取第k条对角线的元素。
(1)提取矩阵的对角线元素
设A为m×n矩阵，diag(A)函数用于提取矩阵A主对角线元素产生一个具有min(m,n)个元素的列向量。
(2)构造对角矩阵
设V为具有m个元素的向量，diag(V)将产生一个m×m对角矩阵，其主对角线元素即为向量V的元素。

    diag(V)函数也有更进一步的形式diag(V,k)，其功能是产生一个n×n(n=m+k)对角阵，其第k条对角线的元素即为向量V的元素。
例1.1 先建立5×5矩阵A，然后将A的第1行元素乘以1，第2行乘以2，…，第5行乘以5。
命令如下：
A=[17,0,1,0,15;23,5,7,14,16;4,0,13,0,22;10,12,19,21,3;11,18,25,2,19];
D=diag([1,2,3,4,5]);

二. 矩阵的三角阵

 (1)下三角矩阵
    求矩阵A的下三角阵的MATLAB函数是tril(A)。
    tril(A)函数也有更进一步的一种形式，即tril(A,k)，其功能是求矩阵A的第k条对角线以下的元素。
 (2)上三角矩阵
    在MATLAB中，提取矩阵A的上三角矩阵的函数是triu(A)和triu(A,k)，其用法与提取下三角矩阵的函数tril(A)和tril(A,k)完全相同。

三、 矩阵结构变换

1. 矩阵的转置
    转置运算符是单撇号(')。
2. 矩阵的旋转
    矩阵的旋转利用函数rot90(A,k)，功能是将矩阵A逆时针旋转90º的k倍，当k为1时可省略。
3. 矩阵的左右翻转
    对矩阵A实施左右翻转的函数是fliplr(A)。
4. 矩阵的上下翻转
   对矩阵A实施上下翻转的函数是flipud(A)。
 

四、 矩阵的逆与伪逆

1. 矩阵的逆
   求一个矩阵的逆非常容易。求方阵A的逆可调用函数inv(A)。
例1.4 用求逆矩阵的方法解线性方程组。
    命令如下：
     A=[1,2,3;1,4,9;1,8,27]; b=[5,–2,6]';  
     x=inv(A)*b
     一般情况下，用左除比求矩阵的逆的方法更有效，即x=A\b。

五、 方阵的行列式

    求方阵A所对应的行列式的值的函数是det(A)。
例1.7用克莱姆(Cramer)方法求解线性方程组。
程序如下：
D=[2,2,-1,1;4,3,-1,2;8,5,-3,4;3,3,-2,2];  %定义 系数矩阵
b=[4;6;12;6];      %定义常数项向量
D1=[b,D(:,2:4)]; %用方程组的右端向量置换D的第1列
 D2=[D(:,1),b,D(:,3:4)];        %用方程组的右端向量置换D的第2列
D3=[D(:,1:2),b,D(:,4:4)];%用方程组的右端向量置换D的第3列
D4=[D(:,1:3),b];%用方程组的右端向量置换D的第4列
 DD=det(D);
 x1=det(D1)/DD;
 x2=det(D2)/DD;
 x3=det(D3)/DD;
 x4=det(D4)/DD;
 [x1,x2,x3,x4]

六、 矩阵的秩

    MATLAB中，求矩阵秩的函数是rank(A)。例如：求例1.7中方程组系数矩阵D的秩，命令是：
       D=[2,2,-1,1;4,3,-1,2;8,5,-3,4;3,3,-2,2];
       r=rank(D)
       r =
            4
       说明D是一个满秩矩阵。

七、 向量和矩阵的范数

1. 计算向量3种常用范数的函数
   (1)norm(V)或norm(V,2)  计算向量V的2—范数
   (2)norm(V,1)  计算向量V的1—范数
   (3)norm(V,inf)  计算向量V的∞—范数
   例1.8 已知V，求V的3种范数。
   命令如下：
   V=[-1,1/2,1];
   v1=norm(V,1)              %求V的1—范数
   v2=norm(V)               %求V的2—范数
   vinf=norm(V,inf)           %求∞—范数
 

---

MATLAB中，计算矩阵A的特征值和特征向量的函数是eig(A)，常用的调用格式有3种：
(1)E=eig(A)  求矩阵A的全部特征值，构成向量E。
(2)[V,D]=eig(A)  求矩阵A的全部特征值，构成对角阵D，并求A的特征向量构成V的列向量。
(3)[V,D]=eig(A,'nobalance')  与第2种格式类似，但第2种格式中先对A作相似变换后求矩阵A的特征值和特征向量，而格式3直接求矩阵A的特征值和特征向量。