如何在R中生成矩阵

Question

我尝试了两种不同的方法来计算矩阵的元素。 我认为这两种方法都是正确的，但是它们会产生不同的结果。 我很好奇，不知道哪一个是正确的。

例如，D是一个（t / 2 + 1）x（t / 2 + 1）矩阵，其元素读为

  D(k,m) = (2/t)*cos(((m-1)(k-1)*pi)/(t/2))

K = M = 1，...，（T / 2 + 1）。

我的第一个代码是：

   m<-seq(1,(t/2)+1,length=(t/2)+1);
   k<-m;

   D<-matrix(NA,nrow = length(k),ncol = length(k)); 

   for(i in 1:length(k)) # row
     for(j in 1:length(m)) # column
    {
     D[i,j]<-(2/t)*cos(((j-1)*(i-1)*pi)/(t*0.5))
    }

我的第二个是：

D<-(2/(t))*cos(as.matrix(seq(0,t/2,by=1))%*%
   t(as.matrix(seq(0,t/2,by=1)))*pi)/(0.5*t)

PI = 3.1415926

这两种方法应该产生相同的结果吗？谢谢您的建议！

Answer 1

将此更改为第二应该可以解决问题：

D<-(2/(t))*cos(as.matrix(seq(0,t/2,by=1))%*%
             t(as.matrix(seq(0,t/2,by=1)))*pi/(0.5*t))

请注意，您需要在cos()内部/(0.5*t) cos() 。

祝你好运！