R软件中有没有把一组向量关于施密特正交化的软件包？谢谢大神赐教。

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

R软件中有没有把一组向量关于施密特正交化的软件包？谢谢大神赐教。

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：施密特正交化 有没有 施密特 软件包 r软件 施密特 软件包

直接用qr就可以了

ltx5151 发表于 2014-3-5 10:12
直接用qr就可以了

QR分解吗？

kfywt 发表于 2014-3-5 10:38
QR分解吗？

是的，直接输入?qr查看帮助文档

写了个R函数，加了个例子，见下：

1. Schmid_orthogonalization<-function(alpha){
   
2.   n<-dim(alpha)[1]
   
3.   s<-dim(alpha)[2]
   
4.   beta1<-alpha[,1]
   
5.   result<-beta1
   
6.   if(s>1){
   
7.     for(i in 2:s){
   
8.       vec<-rep(0,n)
   
9.       for(j in 1:(i-1)){
   
10.         beta0<-get(paste0("beta",j))
    
11.         vec<-vec-(sum(alpha[,i]*beta0)/sum(beta0^2))*beta0
    
12.       }
    
13.       assign(paste0("beta",i),alpha[,i]+vec)
    
14.       result<-cbind(result,get(paste0("beta",i)))
    
15.     }
    
16.   }
    
17.   return(result)
    
18. }
    
19. alpha<-matrix(c(1,1,1,1,
    
20.                 1,-2,-3,-4,
    
21.                 1,2,2,3),byrow=F,ncol=3)#向量组为3个四维列向量排成的矩阵alpha
    
22. library(MASS)
    
23. fractions(Schmid_orthogonalization(alpha))#化为分数形式

复制代码

学而好思 发表于 2020-4-23 12:13
写了个R函数，加了个例子，见下：

怎么发上来是这个格式的？

wuhan0130 发表于 2020-4-27 07:20
怎么发上来是这个格式的？

可以编辑，发帖编辑栏里点<>就行