y=zeros(10000,1);   %模拟10000次
for i=1:10000
  x=rand(0,1);  
  if  x<=p1
  y(i)=-1;
  else
  y(i)=1
end
end

不错，难得有实质性的交流帖子

二楼的方法太麻烦了，MATLAB最大的特点就是支持矩阵运算，所要要尽量少用循环，下面是我有一种比较简单的方法
*********************************************************************************************************
a=rand(1,10000)%产生一万个0-1均匀分布的随机数
a(a<=p1)=-1;
a(a>p1)=1;
a%就是楼主所说的随机数
*********************************************************************************************************
这种方法的运算速度要比二楼的快得多，命令也简单得多

4F code is neat, but a working code need pay attention to the initial condition.
For example insert 'rand('state', 0)' before random number generation. Or,
the computation result may Not be reproducible.