不完全QR分解

最新推荐文章于 2022-12-06 20:02:47 发布

chder_白南

最新推荐文章于 2022-12-06 20:02:47 发布

阅读量265

点赞数

分类专栏：数值分析 fortran 文章标签：数值计算

本文链接：https://blog.csdn.net/chd_lkl/article/details/103598043

版权

数值分析同时被 2 个专栏收录

39 篇文章 11 订阅

订阅专栏

fortran

38 篇文章 16 订阅

订阅专栏

Program QR
  Implicit none
  Integer :: i, j
  Integer, parameter :: m = 4, n = 3
  Real(kind=8), parameter :: A(m,n) = reshape( [ 1.,2.,2.,-4.,3.,2.,2.,5.,-2.,6.,-4.,3. ],[m,n] )
  Real(kind=8) :: Q(m,n) = 0.d0, R(n,n) = 0.d0, y(m), A0(m,n)

  Do j = 1, n
    y = A(:,j)
    Do i = 1, j - 1
      R(i,j) = dot_product(Q(:,i),y)
      y = y - R(i,j)*Q(:,i)
    End do
    R(j,j) = sqrt( dot_product(y,y) )
    Q(:,j) = y / R(j,j)
  End do
  
  Write(*,'(1x,A)') 'The matrix Q is:'
  
  Do i = 1, m
    Write(*,'(*(f14.8))') Q(i,:)
  End do
  
  Write(*,'(1x,A)') 'The matrix R is:'

  Do i = 1, n
    Write(*,'(*(f14.8))') R(i,:)
  End do
  
  A0 = matmul(Q,R)
  Write(*,'(1x,A)') 'The matrix A is:'
  Do i = 1, m
    Write(*,'(*(f14.8))') A0(i,:)
  End do
  
End program QR