%產(chǎn)生n個N(a,b)正態(tài)分布隨機數(shù)
%其中a為均值，b為方差
%function x=normal(a,b,n)
function x=normal(a,b,n)
m=48;%
for i=1:n
    r=rand(1,m);
    x(i)=a+sqrt(b)*(sum(r)-m/2)/sqrt(m/12);
end