资源描述：
《附录二matlab在线性代数中的应用》由会员上传分享，免费在线阅读，更多相关内容在工程资料-天天文库。
1、附录二Matlab在线性代数中的应用§1向量组的线性相关性求列向量组A的一个最大线性无关组可用命令rref(A)将A化成阶梯形的行最简形式，其中单位向量对应的列向量即为最大线性无关组所含向量，其它列向量的坐标即为其对应向量用最大线性无关组线性表示的系数。例1求下列矩阵列向量组的一个最大无关组。⎡1−2−102⎤⎢⎥−2426−6A=⎢⎥⎢2−1023⎥⎢⎥⎣33334⎦解编写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=rre
2、f(a)求得b=101/3016/3012/30-1/90001-1/300000记矩阵A的五个列向量依次为α、α、α、α、α，则α、α、α是列向12345124量组的一个最大无关组。且有121611α=α+α，α=α−α−α.312512433393⎡22−1⎤⎡14⎤⎢⎥⎢⎥例2设A=[a,a,a]=2−12，B=[b,b]=03，123⎢⎥12⎢⎥⎢⎣−122⎥⎦⎢⎣−42⎥⎦3验证a,a,a是R的一个基，并把b,b用这个基线性表示。12312解编写M文件ex2.m如下：formatrata=[2,2,-1;2,-
3、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)，求出所有的-310-基础解系。例3求解下列方程组⎧2x1+x2−5x3+x4=8⎪⎪x1−3x2
4、−6x4=9⎨2x−x+2x=−5⎪234⎪x+4x−7x+6x=0⎩1234解编写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=ab求得solution=[3-4-11]'。例4求超定方程组⎧2x1+4x2=11⎪⎪3x1−5x2=3⎨x+2x=6⎪12⎪⎩2x1+x2=7解编写M文件ex4.m如下：a=[2,4;3,-5;1,2;2,1];b=[11;3;6;7];solution=ab求得so
5、lution=[3.04031.2418]'。例5求解方程组⎧⎪x−x−x+x=0,1234⎪⎨x1−x2+x3−3x4=1,⎪1⎪x−x−2x+x=−.1234⎩2解编写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故方程组有解，并有⎧1x=x+x+124⎪⎪2⎨1⎪x=2x+⎪⎩342-311-§3相似矩阵及二次型有时我们需要精确的特征值和特征向量，就须利用Matla
6、b的符号运算功能。在Matlab中创建符号矩阵和创建数值矩阵的形式很相似，只不过要用到符号定义函数sym。下面介绍使用此函数创建符号函数的几种形式。3.1使用sym函数直接创建符号矩阵此方法和直接创建数值矩阵的方法几乎完全相同。矩阵元素可以是符号表达式，各符号表达式的长度可以不同，矩阵元素之间可用空格或逗号分隔。例如：x=sym('[a+sin(d),b;1/c,d]');y=det(x)求得y=(d*c*a+d*c*sin(d)-b)/c3.2将数值矩阵转化为符号矩阵在Matlab中，数值矩阵不能直接参与符号运算，必须
7、先转化为符号矩阵。例如：a=[2/3,sqrt(2);3,1]a=0.66671.41423.00001.0000b=sym(a)b=[2/3,sqrt(2)][3,1]3.3符号矩阵的索引和修改Matlab的符号矩阵索引和修改同数值矩阵的索引和修改完全相同。例如：对上例中的矩阵b进行修改b(2,2)='log(9)'3.4举例例6求一个正交变换x=Py，把二次型f=2xx+2xx−2xx−2xx+2xx+2xx121314232434化为标准形。解：二次型的矩阵为⎡011−1⎤⎢⎥10−11A=⎢⎥⎢1−101⎥⎢⎥⎣
8、−1110⎦由A=[0,1,1,-1;1,0,-1,1;1,-1,0,1;-1,1,1,0];[P,D]=eig(A)求得P=0.78870.21130.5000-0.28870.21130.7887-0.50000.28870.5774-0.5774-0.50000.2887000.50000.8660D=1.00
