数学建模 matlab在线性代数中的应用.docxVIP

数学建模 matlab在线性代数中的应用 数学建模-matlab在线性代数中的应用 精品文档你我共享 -310- 附录二matlab在线性代数中的应用 §1向量组的线性相关性 求列向量组a的一个最大线性无关组可用命令rref(a)将a化成阶梯形的行最简形 式，其中单位向量对应的列向量即为为最小线性毫无关系胸八色含向量，其它列于向量的座标即为其对应向量用最小线性毫无关系组线性则表示的系数。基准1谋以下矩阵列于向量组的一个最小毫无关系组。 ?????? ???????????= 33334210232426612102a 求解撰写m文件ex1.m如下：formatrat 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=101/3016/3012/30-1/90001-1/300000 记矩阵a的五个列于向量依次为1 α、2 α、3α、4α、5α，则1α、2α、4 腹有诗书气自华 精品文档你我共享资源 α是列向 量组的一个最小毫无关系组。且存有 312 3 23 α=1α+α，512431913 α=16α?α?α. 基准2设 ?????????????==122212221 [,,]123aaaa，????? ??????==420314 [,]12bbb， 验证123a,a,a是r3的一个基，并把12b,b用这个基线性表示。 求解撰写m文件ex2.m如下：formatrat 腹有诗书气自华 精品文档你我共享资源 a=[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解以下方程组 ??????? +?+=?+=???=+?+=4760225369258 12342341241234 xxxxxxxxxxxxxx 解编写m文件ex3.m如下：formatrat a=[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求超定方程组 ?? 腹有诗书气自华 精品文档你我共享资源 ?????+=+=?=+=27263532411 xxxxxxxx 求解撰写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求解方程组 ??????? ??+=??+?=??+=.22131,0, 123412341234 xxxx 深通诗书气自华 精品文档你我共享 xxxxxxxx 解编写m文件ex5.m如下：formatrat a=[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 故方程组存有求解，并存有 ???????=+=++22121 34124 xxxxx -312- §3相似矩阵及二次型 有时我们须要准确的特征值和特征向量，就须利用matlab的符号运算功能。在matlab中建立符号矩阵和建立数值矩阵的形式很相近，只不过必须使用符号定义函数sym。下面了解采用此函数建立符号函数的几种形式。3.1采用sym函数轻易建立符号矩阵 此方法和直接创建数值矩阵的方法几乎完全相同。矩阵元素可以是符号表达式，各符号表达式的长度可以不同，矩阵元素之间可用空格或逗号分隔。例如