数值分析方阵的QR分解

function[]=iqr()
% 实验名称：方阵的QR分解
% 实验描述：先将方阵化为上海申博格阵，再用QR分解法求上海申博格阵的特征值,则所得到的特征值也是方阵的特征值
% 作者:newplan
% 实验完成日期：6月10号
%下面的A为测试三阶的方阵
A=[5,-3,2;6,-4,4;4,-4,5]
%下面的A为测试四阶的方阵
%A = [1 2 1 2;2 2 -1 1;1 -1 1 1;2 1 1 1]
%通过调用malab的自带的函数求得A的所有特征值和特征向量
%特征值保存在v中，特征向量保存的在d中,将其打印出来和我们的算法算出来的特征值进行对比
[v,d]=eig(A)
%求出行和列的大小
msize=size(A);
%取得矩阵的列数，其实行数和列数都为n
n=msize(1);
%生成n阶单位阵
Q=eye(n);
%用household的方法求矩阵A的上海森伯格阵
for i=1:n-2%从第一列开始到倒数第三列 %求出每一列的最大值d=max(abs(A(i+1:n,i)));%规范化U(i+1:n,i)=A(i+1:n,i)/d;delta=U(i+1,i)*norm(U(i+1:n,i))/abs(U(i+1,i));U(i+1,i)=U(i+1,i)+delta;beta = delta*U(i+1,i);%求出R矩阵根据课本316P例题三 R = eye(n-i,n-i)-inv(beta)*U(i+1:n,i)*U(i+1:n,i)';u=eye(n,n);for j =i+1:nfor k =i+1:nu(j,k)=R(j-i,k-i);endendA=u*A*u;%生成新的A=u×A×u
end
%error为我们设定的误差限制
error = 0.0000001;
%flag为判断QR法是否继续进行的标志位
flag =1;
while flag==1
flag =0 ;
R =A;
Q = eye(n,n);
%按照QR分解法求出cos，sin 然后计算V，最终得到R和Q
for i=1:n-1r = norm(R(i:i+1,i));icos=R(i,i)/r;isin=R(i+1,i)/r;v=eye(n,n);v(i,i)=icos;v(i+1,i+1)=icos;v(i,i+1)=isin;v(i+1,i)=-isin;R=v*R;Q=Q*v';
end
%用R*Q的结果去替换A
A =R*Q;
%下面这个循环检测A的精度时候足够，去看A的次对角线各个元素的绝对值是否小于误差限制
for w =2:nif abs(A(w,w-1))>errorflag = 1 ;break;%若有其中一个元素的绝对值还是大于误差限制则还要继续进行QR分解end
end
%判断的过程完毕
end
%把A打印出来
A

转载出处：http://www.cppblog.com/newplan/archive/2008/06/24/54439.html

http://emuch.net/html/201004/1981106.html

数值分析方阵的QR分解相关推荐

数值分析中的QR分解及其代码实现
QR分解若 A∈Cm×k A ∈ C m × k A\in C^{m×k} 是一个列满秩的矩阵,rank(A) = k,则矩阵A 可以分解为 A=QR A = Q R A = QR , Q∈Cm×k ...
矩阵分解——QR分解
文章目录满秩方阵的QR分解矩阵QR分解例题列满秩矩阵的QR分解满秩方阵的QR分解可以看到,该证明过程是构造性的,即通过构造出了 Q Q Q, R R R的方式,证明了QR分解的存在性,不仅证 ...
数值分析（五）：C++实现一般实矩阵的QR分解
先简单介绍一下一般实矩阵的QR分解是什么: 再是源码,不过这次的C++源码是教材上的,提前说明,才五十几行代码就实现了一般实矩阵的QR分解,精炼且巧妙,若是我来写肯定是办不到这么简洁的,正如<C ...
双步位移求解特征值matlab,数值分析——带双步位移的QR分解求特征值算法
C语言实现数值分析中带双步位移的QR分解求特征值算法. 数值分析(B) 大作业(二) 1.算法设计: ①矩阵的拟上三角化: 对实矩阵A进行相似变换化为拟上三角矩阵A(n 1),其变换矩阵采用 ...
基于QR分解迭代求解方阵特征值和特征向量
基于QR分解迭代求解方阵特征值和特征向量一.特征值与特征向量求解的难点线性代数的知识告诉我们如果要求一个方阵的特征值,只需要求解如下的特征方程的根即可: f ( λ ) = ( λ − λ 1 ) ...
Julia 矩阵QR分解和特征值
Julia 矩阵QR分解和特征值前言 1. 施密特正交 (1) 利用施密特正交求出正交矩阵Q (2) 求出上三角矩阵R (3) 改进的消减QR分解 2. 完全QR分解 3. 矩阵QR分解的作用 (1 ...
5.9 QR分解--Gram-Schmidt 分解
5.9 QR分解–Gram-Schmidt 分解最小二乘法需要解方程 ATAx=ATbA^TA\mathbf{x} = A^T\mathbf{b}ATAx=ATb ,需要计算矩阵乘法 ATAA^TA ...
机器学习（十一）——机器学习中的矩阵方法（1）LU分解、QR分解
http://antkillerfarm.github.io/ 因子分析的EM估计(续) 去掉和各参数无关的部分后,可得: ∑i=1mE[logp(x(i)|z(i);μ,Λ,Ψ)]=∑i=1mE[1 ...
crout分解计算例题_矩阵与数值计算（2）——矩阵三角分解LU、PALU、Cholesky三角分解、QR分解...
前言矩阵分解是设计算法的主要技巧,通过分解可以将复杂问题转化为几个简单问题求解,通常完成这一转化任务的主要技巧就是矩阵分解.例如,我们知道上三角矩阵和下三角矩阵是容易求解的,或者对角矩阵是最理想的求 ...

数值分析方阵的QR分解

数值分析方阵的QR分解相关推荐

最新文章

热门文章