x=100*ones(100,1);
max=1;
while max<=100
i=randint(1,1,[1,100]);
x(i)=x(i)-1;
i=randint(1,1,[1,100]);
x(i)=x(i)+1;
max=max+1;
end
subplot(3,1,1)
plot(x)
subplot(3,1,2)
hist(x)
subplot(3,1,3)
histfit(x)
V=var(x);%样本方差
sigma=sqrt(V);%=std(x),std(x)=std(x,1)*sqrt(n/n-1)