1、线性代数在大二上学期,我们学习了线性代数这门课程。为了更好的理解这门课程,并且把线性代数作为一个工具来使用,我们现在初步使用MATLAB这款数学软件来做一些简单的实验。我们从以下六个方面来进行实验:1. MATLAB基础知识与入门2. 矩阵与行列式实验3. 线性方程组实验4. N元向量实验5. 矩阵特征值和特征向量6. 二次型向量MATLAB基础知识与入门MATLAB界面友好,易学易懂,操作方便简单,对于它的基本知识介绍这里就不与一一赘述,后面我们用到时相信读者也能看懂。矩阵与行列实验现在我们就来做一些简单的矩阵与行列实验。生成一个3*3的矩阵A生成一个3*2的矩阵B现在我们就来对这两个矩阵A
2、,B来做一些小实验1) 计算C=A*B2) 计算B+C3) 计算B-C4) 计算3*A5) 计算A36) 计算A的转置矩阵7) 计算方阵A的行列式8) 计算A的秩9) 我们看到A满秩,故A可逆。计算A的逆10) 计算A的行最简形式线性方程组实验我们在这个实验里面,通过MATLAB的高斯消元法(行初等变换)法来求解齐次线性方程组和非齐次线性方程组。1) 齐次线性方程组的解法教材上给出的解齐次线性方程组的方法即:a) 将系数矩阵实施初等行变换,变为行最简形式b) 取线性无关的自由向量,带入方程组中得到解向量,即方程组的基础解系c) 对基础解系进行线性组合,即得齐次线性方程组的解现在,我们使用MAT
3、LAB对齐次线性方程组进行同样的方法以解之。a) 将系数矩阵转换为符号矩阵,并解出其行最简形式b) 解出解空间的一个基础解系c) 对基础解系进行线性组合给出一个例子(教材95页例3.6)AX=0其中,X=(x1,x2,x3,x4,x5)我们先对系数矩阵进行第一步运算 转换为符号矩阵 变换为行最简形式 现在进行第二步运算 求其基础解系 由此我们得到其基础解系为 1=(-1 0 1 0 0)T2=(-1 0 0 1 0)T3=(-3 2 0 0 1)T 故该齐次线性方程组的通解为x=k1 *1+k2* 2+k3* 3(k1,k2,k3可取任意值)2) 非齐次线性方程组的解法(我们讨论有无穷多解的情
4、况) 教材上给出的解非齐次线性方程组的方法即:a) 对增广矩阵进行初等行变换,变为行最简形式b) 找出通解,由自由未知量解出方程组一个特解c) 将方程组变换为齐次线性方程组d) 解出通解e) 解出基础解系f) 对基础解系线性组合,再加上特解即可现在,我们使用MATLAB对齐次线性方程组进行同样的方法以解之。 a) 对增广矩阵符号化,并变为行最简形式b) 找出通解,由自由未知量解出方程组一个特解c) 对系数矩阵进行行最简变换d) 解出基础解系e) 对基础解系线性组合,再加上特解即可同样给出一个例子(教材99页,例3.11) AX=B先找出其增广矩阵 对其符号化 化为行最简形式 可知,x4,x5为
5、自由未知量,令x4=x5=0,带入解得一个特解 =(72/2 -11 -9/4 0 0 )T 系数矩阵 符号化 行最简变换 解出基础解系 1=(-19/2 -4 3/4 1 0)T2=(-4 1 0 0 1)T 对基础解系线性组合,再加上特解 即:x=+k1* 1+k2* 2(k1,k2为任意常数)N元向量实验下面我们对N元向量组进行一些实验。1) 计算向量组的秩先创建一个三元向量,并组合然后计算其秩2) 向量组的线性相关向量组的线性相关性我们同样可以用之前用过的”rref()”来计算创建5元向量组求秩所以,该向量组的秩为3,线性相关,极大无关组为R的前三列。矩阵特征值和特征向量 下面我们来做
6、一个特征值和和特征向量的实验。 创建一个矩阵 输出其多项式输出其特征值D和特征向量V 实对称矩阵的正交相似对角化 我们用”eig”函数就可以为实对称矩阵A找到正交矩阵Q.创建矩阵A 输出特征值D和正交矩阵Q二次型向量二次型的标准型 一个二次型可以等同于一个实对称矩阵,因此在上一个实验中有关实对称矩阵正交相似对角化的计算方法就可以来求二次型的标准型。因此在这里就不做过多介绍。 二次型和实对称矩阵的正定性 一个数值二次型的正定性可以用MATLAB来判定,方法如下:1) 用eig求其特征值,看是否全为正数2) 用poly求出其特征多项式,再用roots求出根,看是否全为正数3) 按照Sylvester定理,验证各阶顺序主子式是否都大与零、至此,我们就将MATLAB的实验介绍完了。
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1