您所在位置:网站首页 > 海量文档
 > 计算机 > matlab
MATLAB论文_矩阵的QR分解及其MATLAB实现.doc8页
本文档一共被下载:次,您可全文免费在线阅读后下载本文档。
下载提示
1.本站不保证该用户上传的文档完整性,不预览、不比对内容而直接下载产生的反悔问题本站不予受理。
2.该文档所得收入(下载+内容+预览三)归上传者、原创者。
3.登录后可充值,立即自动返金币,充值渠道很便利
《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阶,且m>n,则运行后,输出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阶,且m>n,则运行后,输出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阶,且m>n,则运行后,输出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
发表评论
请自觉遵守互联网相关的政策法规,严禁发布色情、暴力、反动的言论。
用户名:
验证码:
匿名?
发表评论