急求，用sas编程求解QR 分解，（不使用sas自带的函数）

矩阵a
a={1 0.2 0.3,
0.2 5 0.1,
0.3 0.1 6}

这个是计算机精度问题，(比如：你不可能让计算机显示根号3)吧实际上已经是上三角矩阵了，下三角部分无穷小了，实际上就是0了。如果你看不顺眼下三角部分的花，就把下三角部分变为0好了，方法如下：在print Q;语句上面加一段代码：

1.                 do i=1 to nrow(R);
   
2.                         do j=1 to i-1;
   
3.                                 R[i,j]=0;
   
4.                         end;
   
5.                 end;

复制代码

1. proc iml;
   
2.         A={1 0.2 0.3,0.2 5 0.1,0.3 0.1 6};
   
3.         Q=j(nrow(A),ncol(A),.);
   
4.         do i=1 to ncol(A);
   
5.                 T=j(ncol(A),1,0);
   
6.                 do j=1 to i-1;
   
7.                         T=T+(sum(A[,i]#Q[,j])/Q[##,j])#Q[,j];
   
8.                 end;
   
9.                 Q[,i]=A[,i]-T;
   
10.                 Q[,i]=Q[,i]/sqrt(Q[##,i]);
    
11.         end;
    
12.         R=Q`*A;
    
13. 
    
14.         print Q;
    
15.         print R;
    
16. quit;

复制代码

learsaas 发表于 2016-4-24 11:27

请问答主，用您的代码进行QR分解，得到的R矩阵不是具有正对角元的上三角矩阵啊？

这是根据代码得到的矩阵

learsaas 发表于 2016-4-22 18:31
这个是计算机精度问题，(比如：你不可能让计算机显示根号3)吧实际上已经是上三角矩阵了，下三角部分无穷小了 ...

太感谢了，答主666666