MATLAB论文12010245292岳安健.docVIP

《MATLAB语言》课程论文 矩阵的QR分解及其MATLAB实现 姓 名：岳安健 学 号：12010245292 专 业：电子信息工程 班 级：2010级电子班 指导老师：汤全武 学 院：物电学院 完成日期：2011.12.11 矩阵的QR分解及其MATLAB实现 （岳安健 12010245292 2010级电子班） [摘 要] 矩阵的QR分解是计算中小型矩阵的全部特征值问题最有效，最简便的方法之一，收敛速度快，算法稳定。在MATLAB系统中，矩阵的QR分解也是比较方便，实用的，主要应用qr函数来进行分解。 [关键字] MATLAB程序 矩阵 QR分解 一、问题的提出 MATLAB是英文MATrix LABoratory矩阵实验室）的缩写。MATLAB的核心采用C语言程序编写，功能强大，拥有可靠的数值计算和符号计算、功能强大的绘图功能、简单易学的语言体系以及为数众多的应用工具箱。 在MATLAB系统中的数值计算以矩阵为数据操作的基本单位，MATLAB提供了丰富的计算函数，使矩阵运算变得非常简捷、方便、高效。 QR分解是将n阶方阵A分解为A=QR的形式，其中Q为n阶正交矩阵，R为n阶上三角矩阵，故QR分解又称正交三角分解。如果我们可以利用MATLAB求解QR矩阵，那么，就可以简洁、高效的处理QR矩阵问题。所以，我们要知道矩阵QR分解的输入命令及命令功能，以及使用的实例。 二、QR分解 矩阵QR分解的输入命令及命令功能。 1.命令：[Q,R]=qr(A) 功能： a.如果输入的n阶矩阵A是非奇异的，运行后输出的n阶非奇异的上三角形矩阵R和n阶正交矩阵Q使得A=QR; b.如果输入的n阶矩阵A是奇异的, 运行后输出的n阶奇异的上三角形矩阵R和n阶正交矩阵Q使得A=QR; c.如果输入的矩阵A是mn阶，且mn，则运行后，输出mn阶上三角形矩阵R（其中R的m-n行和这行以后的全为零）和m阶正交矩阵Q，使得A=QR，其中实际上仅仅是Q的前n列与R计算。 2.命令：[Q,R,E]=qr(A) 功能： a.如果输入的n阶矩阵A是非奇异的，运行后输出的n阶非奇异上三角形矩阵R和n阶正交矩阵Q使得AE=QR; b.如果输入的n阶矩阵A是奇异的, 运行后输出的n阶奇异的上三角形矩阵R和n阶正交矩阵Q使得AE=QR; c.如果输入的矩阵A是mn阶，且mn，则运行后，输出mn阶上三角形矩阵R（其中R的m-n行和这行以后的全为零）和m阶正交矩阵Q，使得AE=QR，其中实际上仅仅是Q的前n列与R计算。 3.命令：[Q,R]=qr(A,0) 功能：[Q,R] =qr(A,0)产生“econormy size”分解。 a.如果输入的是n阶方阵A，则运行后输出结果与[Q,R]=qr(A)相同。 b.如果输入的矩阵A是mn阶，且mn，则运行后，输出mn阶上三角形矩阵R和mn阶矩阵Q（其中Q的列向量两两正交），使得Q的前n列与R的积等于A。 4.命令：[Q,R,E]=qr(A,0) 功能：[Q,R,E]=qr(A,0)产生“econormy size”分解，使得QR=A（：,E），其中E是置换向量。 在了解QR分解在MATIAB程序中的命令及其功能后，来应用MATLAB程序解决实际学习、生活中的应用问题。 下面先通过一个简单的例子来了解下QR分解： 问题1： 对矩阵 A= 进行QR分解。 输入命令如下： A=[1,-1,1;5,-4,3;2,7,10]; %输入矩阵 [Q,R]=qr(A) %对矩阵A进行QR分解 屏幕显示结果如下： Q = -0.1826 -0.0956 -0.9785 -0.9129 -0.3532 0.2048 -0.3651 0.9307 -0.0228 R = -5.4772 1.2780 -6.5727 0 8.0229 8.1517