数学软件课程设计7100115.docVIP

东北大学秦皇岛分校 数学软件实习报告 设计题目 矩阵的应用 系 别 信息与计算科学 专 业 学号 姓名 指导教师 姜玉山 郭静梅 成 绩 教师评语： 指导教师签字： 2011年12月31日1 绪 论 Matlab软件对已知的矩阵对应的行（列）向量组、子矩阵及扩展矩阵，实施矩阵的初等变换及线性无关向量组的正交规范化，确定线性相关向量组的一个极大线性无关向量组，且将其余向量用极大线性无关向量组，并编辑文件来完成所有的实验题目。如调用矩阵A的部分行、BMatlab有关线性代数运算的指令，计算实验主要讨论了矩阵除法、线性方程组的通解和矩阵相似对角化问题；建模实验研究了投入产出分析、基因遗传等应用问题。 1.2【论文关键词】 矩阵 初等变换 向量 1.3【课题的背景 用矩阵的理论与方法来处理现代工程技术中的各种问题已越来越普遍。在工程技术中引进矩阵理论不仅使理论的表达极为简捷，而且对理论的实质刻画也更为深刻，这一点是不容置疑的，更由于计算机和计算方法的普及发展，不仅为矩阵理论的应用开辟了广阔的前景，也使工程技术的研究发生新的变化，开拓了崭新的研究途径，例如系统工程，优化方法，稳定性理论等，无不与矩阵理论发生紧密结合。因此矩阵的理论与方法已成为研究现代工程技术的数学基础。i 行作一个行向量：ai=A( i , : )；ai=A( : , j )； （3）选择A的某几行、某几列上的交叉元素作A的子矩阵： A(起始行：步长：终行，起始列：步长：终列)；当步长为1时可省略； （4）n阶单位阵：eye(n)；n阶零矩阵：zeros(n)； （5）作一个m维的索引向量L，然后取A( : , L )，将依次选取L中的元素对应的列形成子矩阵； （6）将非奇异矩阵A正交规范化：orth(A)；验证矩阵A是否为正交矩阵，只需做A*A ’ 看其是否得到单位矩阵E； （7）两个行向量的内积为：a1*a2 ’ ； （8）让A的第 i 行与第 j 行互换可用赋值语句：A( [i , j ], : )= A( [j , i ], : )； （9）让K乘A的第 i 行可用赋值语句：A( i , : ) =K*A( i , : )； （10）让A的第 i 行加上第 j 行的K倍可用赋值语句：A( i , : ) =A( i , : ) +K*A(j , : ) ； （11）求列向量组A的一个极大线性无关向量组可用命令：rref(A) 将A化为阶梯型的最简形式，其中单位向量对应的列向量即为极大线性无关向量组所含向量，其他列向量的坐标即为其对应向量用极大线性无关向量组线性表示的系数。 2 论文正文 2.1【程序及其运行结果】 A = [1 22 3 41 51 6 74 8 9 11 2 32 4 5 11 44 77 5 4 81 9 12 57 8 12 21 4 8 9 32 4 5 77 6 32 12] B = [2 45 21 2 3 22 4 56 45 11 2 1 51 4 12 11 5 4 4 55 66 44 1 2 1 11 22 55 2 3 7 88 55 4 6 11] % 1.做出A的行向量 A(1,:),A(2,:),A(3,:),A(4,:),A(5,:) A = 1 22 3 41 51 6 74 8 9 11 2 32 4 5 11 44 77 5 4 81 9 12 57 8 12 21 4 8 9 32 4 5 77 6 32 12 B = 2 45 21 2 3 22 4 56 45 11 2 1 51 4 12 11 5 4 4 55 66 44 1 2 1 11 22 55 2 3 7 88 55 4 6 11 ans = 1 22 3 41 51 6 ans = 74 8 9 11 2 32 ans = 4 5 11 44 77