怎样产生一组矩阵

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

怎样产生一组矩阵使得每个元素是\[0.5^{|i-j|}\]我的程序是：
  p=10
  rho=0.5
  a=NULL
  i=1
  for(j1 in 2:(p+1)){
    for (j2 in 2:(p+1)){
      a=rho^(abs(j1-j2))
      i=i+1
    }
  }
a


这个当p=1000时程序太慢了，请教大神！

扫码加我 拉你入群

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

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

分享1 收藏0 回帖

关键词：LaTeX atex Late null For 程序 元素

1. p=10
   
2. 0.5^abs(outer(1:p,1:p,"-"))

复制代码

首先，你先生成一个向量[LaTex]a=[1 0.5 0.5^2 ... 0.5^(p-1)] [/LaTex],则你矩阵的(i,j)元素就是a(abs(i-j)+1)，这样可以减少大量的乘法运算，进行赋值就行啦

ljchuan1988 发表于 2014-12-23 16:04
首先，你先生成一个向量,则你矩阵的(i,j)元素就是a(abs(i-j)+1)，这样可以减少大量的乘法运算，进行赋值就行 ...

这种方法似乎不行！谢谢！

ljchuan1988 发表于 2014-12-23 16:04
首先，你先生成一个向量,则你矩阵的(i,j)元素就是a(abs(i-j)+1)，这样可以减少大量的乘法运算，进行赋值就行 ...

这种方法似乎不行！谢谢！

a<-matrix(0,1000,1000)

for(i in 1:1000){
    for(j in 1:1000){
    a[i,j]<-0.5^abs(i-j)
   }
}

yywan0913 发表于 2014-12-23 15:56

这个很快啊，p=1000的时候，我已经试过了。

b<-matrix(rep(1:1000,1000),1000)
e<-t(b)
a<-0.5^(abs(b-e))

一诺9257 发表于 2014-12-23 16:24
这种方法似乎不行！谢谢！

其实是可以的，你要仔细想想，理解运算费时的原因，这种肯定非常快的

ljchuan1988 发表于 2015-1-9 10:22
其实是可以的，你要仔细想想，理解运算费时的原因，这种肯定非常快的

谢谢！哈哈