-310-附录二Matlab在线性代数中的应用§1向量组的线性相关性求列向量组A的一个最大线性无关组可用命令rref(A)将A化成阶梯形的行最简形式，其中单位向量对应的列向量即为最大线性无关组所含向量，其它列向量的坐标即为其对应向量用最大线性无关组线性表示的系数。例1求下列矩阵列向量组的一个最大无关组。⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡−−−−−=43333320126624220121A解编写M文件ex1.m如下：formatrata=[1,-2,-1,0,2;-2,4,2,6,-6;2,-1,0,2,3;3,3,3,3,4];b=rref(a)求得b=101/3016/3012/30-1/90001-1/300000记矩阵A的五个列向量依次为1α、2α、3α、4α、5α，则1α、2α、4α是列向量组的一个最大无关组。且有2133231ααα+=，42153191316αααα−−=.例2设⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−−−==221212122],,[321aaaA，⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡−==243041],[21bbB，验证321,,aaa是3R的一个基，并把21,bb用这个基线性表示。解编写M文件ex2.m如下：formatrata=[2,2,-1;2,-1,2;-1,2,2];b=[1,4;0,3;-4,2];c=rref([a,b])求得c=1002/34/3010-2/31001-12/3§2线性方程组Matlab中解线性方程组可以使用“\”。虽然表面上只是一个简简单单的符号，而它的内部却包含许许多多的自适应算法，如对超定方程用最小二乘法，对欠定方程它将给出范数最小的一个解，解三对角阵方程组时用追赶法等。另外欠定方程组可以使用求矩阵A的阶梯形行最简形式命令rref(A)，求出所有的-311-基础解系。例3求解下列方程组⎪⎪⎩⎪⎪⎨⎧=+−+−=+−=−−=+−+067452296385243214324214321xxxxxxxxxxxxxx解编写M文件ex3.m如下：formatrata=[2,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6];b=[8;9;-5;0];solution=a\b求得solution=[3-4-11]'。例4求超定方程组⎪⎪⎩⎪⎪⎨⎧=+=+=−=+7262353114221212121xxxxxxxx解编写M文件ex4.m如下：a=[2,4;3,-5;1,2;2,1];b=[11;3;6;7];solution=a\b求得solution=[3.04031.2418]'。例5求解方程组⎪⎪⎩⎪⎪⎨⎧−=+−−=−+−=+−−.212,13,0432143214321xxxxxxxxxxxx解编写M文件ex5.m如下：formatrata=[1,-1,-1,1,0;1,-1,1,-3,1;1,-1,-2,3,-1/2];b=rref(a)求得：b=1-10-11/2001-21/200000故方程组有解，并有⎪⎪⎩⎪⎪⎨⎧+=++=2122143421xxxxx-312-§3相似矩阵及二次型有时我们需要精确的特征值和特征向量，就须利用Matlab的符号运算功能。在Matlab中创建符号矩阵和创建数值矩阵的形式很相似，只不过要用到符号定义函数sym。下面介绍使用此函数创建符号函数的几种形...
