求助：用R怎么产生Hadamard矩阵，有c语言程序，改编

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

Hadamard Matrix在matlab里面可以用函数hadamard(n)实现，可是我的程序是在R编的，希望可以用R来实现。
本人水平太低，搞了很久都没有弄出来，希望各位大虾帮忙看看。
我有一个c语言的程序是这么编的：
/* 产生Hadamard矩阵(hmat)*/
  for (u = 0;u < N; u++)
    for (v = 0;v < N;v++)
      {        
cnt = 0;
for (uv=u&v; uv ; uv >>= 1 ) {
   cnt += uv & 1;
}
        *(hmat+(N*u)+v) = ((cnt & 1) ? -1 : 1);
      }
/***** 将Hadamard矩阵转换成顺序形式sequency order ******/
/* ncount - vector of number of zero crossings in each row of hmat*/
  for (i = 0;i < (N-1);i++)
    for (j = 0;j < N;j++)
      if (*(hmat+i+(j*N)) != (*(hmat+(i+1)+(j*N)))) ncount[j] += 1;

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：c语言程序 Adam Had dam C语言 矩阵 Matrix matlab c语言

survey包里有个hadamard函数，不知道是不是你想要的。

qoiqpwqr 发表于 2012-4-12 09:51
survey包里有个hadamard函数，不知道是不是你想要的。

太感谢了，有希望！！本来想去下载那个包，可是r-project一直打不开，也不知道是怎么了。。。

qoiqpwqr 发表于 2012-4-12 09:51
survey包里有个hadamard函数，不知道是不是你想要的。

请问您是在哪里找到这个函数的，我安装了程序包servey但是调用这个函数说没有哦。

mengqinqing 发表于 2012-4-12 14:53
请问您是在哪里找到这个函数的，我安装了程序包servey但是调用这个函数说没有哦。

install.packages("survey")
library(survey)

help(hadamard)



Hadamard matrices
Description
Returns a Hadamard matrix of dimension larger than the argument.

ywh19860616 发表于 2012-4-12 15:05
install.packages("survey")
library(survey)

找到了，的确是我要的。非常感谢！

ywh19860616 发表于 2012-4-12 15:05
install.packages("survey")
library(survey)

R的hadamard矩阵好奇怪，我认为hadamard(8)就会输出一个8阶矩阵，结果输出了一个12阶矩阵。不知道您知不知这是怎么回事？

不知道这个矩阵的具体构造，按照命令提示是说：
Returns a Hadamard matrix of dimension larger than the argument.



你看下这个pdf文档：http://designtheory.org/library/encyc/topics/had.pdf

ywh19860616 发表于 2012-4-16 11:55
不知道这个矩阵的具体构造，按照命令提示是说：
Returns a Hadamard matrix of dimension larger than the ...

好的，谢谢！

mengqinqing 发表于 2012-4-16 11:43
R的hadamard矩阵好奇怪，我认为hadamard(8)就会输出一个8阶矩阵，结果输出了一个12阶矩阵。不知道您知不知 ...

%H=hadamard(n),  H'*H = N*EYE(N).

%%%%in Matlab
n=8;
hadamard(n)

ans =

     1     1     1     1     1     1     1     1
     1    -1     1    -1     1    -1     1    -1
     1     1    -1    -1     1     1    -1    -1
     1    -1    -1     1     1    -1    -1     1
     1     1     1     1    -1    -1    -1    -1
     1    -1     1    -1    -1     1    -1     1
     1     1    -1    -1    -1    -1     1     1
     1    -1    -1     1    -1     1     1    -1


#####in R
library(survey)
n=8
2*hadamard(n-4)-1

     [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8]
[1,]    1    1    1    1    1    1    1    1
[2,]    1   -1    1   -1    1   -1    1   -1
[3,]    1    1   -1   -1    1    1   -1   -1
[4,]    1   -1   -1    1    1   -1   -1    1
[5,]    1    1    1    1   -1   -1   -1   -1
[6,]    1   -1    1   -1   -1    1   -1    1
[7,]    1    1   -1   -1   -1   -1    1    1
[8,]    1   -1   -1    1   -1    1    1   -1