matlab在线性代数中的应用金融数据挖掘.pdfVIP

附录二 Matlab 在线性代数中的应用 §1 向量组的线性相关性 求列向量组A 的一个最大线性无关组可用命令rref(A)将 A 化成 梯形的行最简形 式，其中单位向量对应的列向量即为最大线性无关组所含向量，其它列向量的坐标即 为其对应向量用最大线性无关组线性表示的系数。 例 1 求下列矩阵列向量组的一个最大无关组。  1 − 2 − 1 0 2    A = − 2 4 2 6 − 6   2 − 1 0 2 3     3 3 3 3 4  解 编写 M 文件 ex1.m 如下： format rat a=[1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4]; b=rref(a) 求得 b = 1 0 1/3 0 16/3 0 1 2/3 0 -1/9 0 0 0 1 -1/3 0 0 9 0 0 A α α α α α α α α 记矩阵 的五个列向量依次为 、 、 、 、 ，则 、 、 是列向量组的 1 2 3 4 5 1 2 4 一个最大无关组。且有 1 2 16 1 1 α3 = α1 + α2 ， α5 = α1 − α2 − α3 . 3 3 3 9 3 2 2 − 1 1 4         例 2 设A = [a ,a ,a ] = 2 − 1 2 ，B = [b ,b ] = 0 3 ， 1 2 3   1 2   − 1 2 2  − 4 2      验证a ,a ,a 是R 3 的一个基，并把b ,b 用这个基线性表示。 1 2 3 1 2 解 编写 M 文件 ex2.m 如下： format rat a=[2,2,-1;2,-1,2;-1,2,2];b=[1,4;0,3;-4,2]; c=rref([a,b]) 求得 c= 1 0 0 2/3 4/3 0 1 0 -2/3 1 0 0 1 -1 2/3 §2 线性方程组 Matlab 中解线性方程组可以使用“\” 。虽然表面上只是一个简简单单的符号，而它 的内部却包含许许多多的自适应算法，如对超定方程用最小二乘法，对欠定方程它将 给出范数最小的一个解，解三对角阵方程组时用追赶法等。 另外欠定方程组可以使用求矩阵A 的 梯形行最简形式命令 rref(A) ，求出所有的