应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 程序语言 » MATLAB/Mathematica » 求矩阵QR分解的源代码程序
北京石油化工学院2026年研究生招生接收调剂公告
51/1返回列表
查看: 4333  |  回复: 2
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

二驴

新虫 (小有名气)

[求助] 求矩阵QR分解的源代码程序 已有2人参与

不要matlab中自带的qr分解。需要源代码。那个大神帮帮忙。矩阵如图。
求矩阵QR分解的源代码程序
1.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-12-03 10:03:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

fish.yfyh

铜虫 (小有名气)

【答案】应助回帖

我提供一个Fortran版本的代码:
======================Fortran源码=============================
      subroutine qr_matrix_decomp(m,n,A, Q,R)
      ! 采用Gram-Schmit正交化方法进行QR分解
      implicit real*8 (a-h,o-z)
      integer ::m,n
      real*8  :: A(m,n),Q(m,m),R(m,n),vec_tmp(m)

      R(1,1) = dsqrt(dot_product(A(:,1),A(:,1)))
      Q(:,1) = A(:,1)/R(1,1)

      do k=2,n
          do j=1,k-1
              R(j,k) = dot_product(Q(:,j),A(:,k))
          end do
          vec_tmp = A(:,k)
          do j=1,k-1
              vec_tmp = vec_tmp - Q(:,j)*R(j,k)
          end do
          R(k,k) = dsqrt(dot_product(vec_tmp,vec_tmp))
          Q(:,k) = vec_tmp/R(k,k)
      end do
      return
      end subroutine


      program qr_main
      implicit real*8 (a-h,o-z)
      integer :: i
      real*8  :: A(4,4),Q(4,4),R(4,4)
      
      A = reshape( (/1.0, 2.0,-1.0, 1.0, &
                      2.0, 5.0, 0.0, 3.0, &
                      1.0, 7.0, 9.0, 2.0, &
                      8.0,-1.0,-2.0, 1.0 /),(/4,4/))
      
      call qr_matrix_decomp(4,4,A, Q,R)
      
      print*, 'Q='
      do i = 1,4
          print '(4(f14.8,x))',Q(i,
      end do
      print*, 'R='
      do i = 1,4
          print '(4(f14.8,x))',R(i,
      end do
      
      stop
      end program
============================================================
编译后运行，结果为
Q=
    0.37796447    -0.05902813     0.14256649     0.91287093
    0.75592895     0.29514067     0.45621276    -0.36514837
   -0.37796447     0.88542200     0.19959308     0.18257419
    0.37796447     0.35416880    -0.85539892     0.00000000
R=
    2.64575131     5.66946710     3.02371578     3.40168026
    0.00000000     2.42015348    10.68409218    -2.18404094
    0.00000000     0.00000000     3.42159569    -0.57026595
    0.00000000     0.00000000     0.00000000     7.30296743
Matlab的结果为：
>> A

A =

     1     2    -1     1
     2     5     0     3
     1     7     9     2
     8    -1    -2     1

>> [Q,R]=qr(A')

Q =

    0.3780   -0.0590    0.1426    0.9129
    0.7559    0.2951    0.4562   -0.3651
   -0.3780    0.8854    0.1996    0.1826
    0.3780    0.3542   -0.8554    0.0000


R =

    2.6458    5.6695    3.0237    3.4017
    0.0000    2.4202   10.6841   -2.1840
    0.0000    0.0000    3.4216   -0.5703
    0.0000    0.0000    0.0000    7.3030
============================================================
赞一下 回复此楼
3楼2013-12-13 14:13:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答

hppdyx

木虫 (知名作家)

【答案】应助回帖

CODE:
function [Q,R]=qrgv(A)
% 基于Givens变换，将方阵A分解为A=QR，其中Q为正交矩阵，R为上三角阵
%
% 参数说明
% A：需要进行QR分解的方阵
% Q：分解得到的正交矩阵
% R：分解得到的上三角阵
%
% 实例说明
% A=[-12 3 3;3 1 -2;3 -2 7];
% [Q,R]=qr(A) % 调用MATLAB自带的QR分解函数进行验证
% [q,r]=qrgv(A) % 调用本函数进行QR分解
% q*r-A % 验证 A=QR
% q'*q % 验证q的正交性
% norm(q) % 验证q的标准化，即二范数等于1
%
% 线性代数基础知识
% 1.B=P*A*inv(P)，称A与B相似，相似矩阵具有相同的特征值
% 2.Q*Q'=I，称Q为正交矩阵，正交矩阵的乘积仍为正交矩阵
%
% by dynamic of Matlab技术论坛
% see also http://www.matlabsky.com
% contact me matlabsky@gmail.com
% 2010-01-17 22:51:18
%
n=size(A,1);
R=A;
Q=eye(n);
for i=1:n-1
    for j=2:n-i+1
        x=R(i:n,i);
        rt=givens(x,1,j);
        r=blkdiag(eye(i-1),rt);
        Q=Q*r';
        R=r*R;
    end
end

function [R,y]=givens(x,i,j)
% 求解标准正交的Given变换矩阵R，使用Rx=y，其中y(j)=0，y(i)=sqrt(x(i)^2+x(j)^2)
%
% 参数说明
% x：需要进行Givens变换的列向量
% i：变为sqrt(x(i)^2+x(j)^2)的元素下标
% j：变为0的元素的下标
% R：Givens变换矩阵
% y：Givens变换结果
%
% 实例说明
% x=[1 3 5 9 6]'; % 将3等效到9上
% [R,y]=givens(x,4,2) % 注意3的下标为2，9的下标为4
% R*x-y % 验证Rx=y
% R'*R % 验证正交性
% norm(R) % 验证标准性，就是范数为1
%
% 关于Givens变换说明
% 1.Givens矩阵是标准正交矩阵，也叫平面旋转矩阵，它是通过坐标旋转的原理将元素j的数值等效到元素i上
% 2.Givens变换每次只能将一个元素变为0，而Householder变换则一次可以将任意个元素变为0
% 3.Givens变换常用于将矩阵A变为对角阵
%
xi=x(i);
xj=x(j);
r=sqrt(xi^2+xj^2);
cost=xi/r;
sint=xj/r;
R=eye(length(x));
R(i,i)=cost;
R(i,j)=sint;
R(j,i)=-sint;
R(j,j)=cost;
y=x(:);
y([i,j])=[r,0];
赞一下 回复此楼
不以风骚惊天下，但求淫荡动世人
2楼2013-12-13 13:10:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 3 个回答
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 323分（计算机视觉和大模型项目）能直接上手 +3 chaoxiicy 2026-04-01 3/150 2026-04-05 00:50 by chongya
[考研] 320求调剂 +3 一样圆 2026-04-04 3/150 2026-04-04 22:29 by 啵啵啵0119
[考研] 292分，材料与化工，申请调剂 +22 程晴之 2026-04-01 26/1300 2026-04-04 22:03 by hemengdong
[考研] 一志愿华北电力大学（北京），材料科学与工程学硕265，求调剂 +11 yelck 2026-04-03 12/600 2026-04-04 19:52 by dongzh2009
[论文投稿] 求文献 5+3 ys879651$ 2026-04-02 3/150 2026-04-04 17:22 by bobvan
[考研] 材料专硕(0856) 339分求调剂 +5 哈哈哈鹅哈哈哈 2026-04-04 5/250 2026-04-04 16:31 by dongzh2009
[考研] 085600，专业课化工原理，320分求调剂 +8 大馋小子 2026-03-29 8/400 2026-04-04 16:21 by dongzh2009
[考研] 复试调剂 +6 范根培 2026-04-04 6/300 2026-04-04 14:27 by 土木硕士招生
[考研] 材料专硕322分 +7 哈哈哈吼吼吼哈 2026-04-04 7/350 2026-04-04 12:09 by asdhh1991
[考研] 070300化学求调剂 +16 小黄鸭宝 2026-03-30 16/800 2026-04-04 11:49 by asdhh1991
[考研] 材料295 +13 小英11 2026-04-03 14/700 2026-04-04 09:02 by 来看流星雨10
[考研] 求调剂，一志愿北京中医药大学 +3 小小达不溜 2026-04-02 3/150 2026-04-03 22:55 by 冲矢昴星团
[考研] 320求调剂 +3 农业工程与信息� 2026-04-03 3/150 2026-04-03 11:40 by 土木硕士招生
[考研] 重庆大学材料与化工085600，初试370+，求求调剂建议 +8 shzhou_ 2026-04-01 9/450 2026-04-03 09:31 by 蓝云思雨
[考研] 260求调剂 +3 朱芷琳 2026-04-02 3/150 2026-04-03 08:44 by yulian1987
[考研] 一志愿大工学硕，求调剂 +4 yub0811 2026-04-02 4/200 2026-04-02 21:36 by 百灵童888
[考研] 求调剂求调剂 +7 121. 2026-04-02 7/350 2026-04-02 19:15 by dick_runner
[考研] 318求调剂 +3 笃行致远. 2026-03-31 4/200 2026-04-02 15:56 by Jaylen.
[考研] 材料工程322分 +8 哈哈哈吼吼吼哈 2026-04-01 8/400 2026-04-02 11:53 by 3041
[考研] 339求调剂 +5 zjjkt 2026-03-31 5/250 2026-04-01 09:18 by JourneyLucky
信息提示
关闭
请填处理意见
关闭