%生成兩組互相獨立的高斯隨機數 期望為m,方差為n^2,個數為num
function res=spRand(m,n,num)
   rand( 'state',sum(100*clock) );
   d1=rand(1,num);
   rand( 'state',7*sum(100*clock)+3);
   d2=rand(1,num);
   xi=n*sqrt(-2*log(d1)).*cos(2*pi*d2)+m;%%正交獨立的高斯序列 N(m,n^2)
   xq=n*sqrt(-2*log(d1)).*sin(2*pi*d2)+m;%%
   res(1,:)=xi;
   res(2,:)=xq;