产生高斯分布的随机变量

for j=1:1:1000
    aa(j) = sqrt(-2*log(rand(1,1)))*cos(2*pi*rand(1,1));
end
bb=randn(1,1000) ;

figure(1)
histfit(aa);
normplot(aa);
figure(2)
 y=aa;
x=linspace(min(y),max(y),100) ;
yy=hist(y,x) ;
yy=yy/length(y); %计算各个区间的个数
bar(x,yy) %画出概率密度分布图

figure(3)
histfit(bb);

figure(4)
normplot(bb);

 [muhat,sigmahat,muci,sigmaci]=normfit(aa);
 muhat
 sigmahat
 

参考：

[1]. http://blog.sina.com.cn/s/blog_70a14458010155b8.html