white noise sequence.txt

1．产生[0,1]均匀分布的白噪声序列 (1) 打印出前50个数 (2) 分布检验 (3) 均值检验 (4) 方差检验 (5) 计算相关函数 Bx(i),i=0,±1,±2,…, ±10

for n=1:2000
xt(n)=unifrnd(0,1);	%产生2000个（0，1）均匀分布白序列
end
subplot(2,1,1);
plot(xt),title('2000个(0,1)均匀分布的白噪声');	
for i=1:5
    for j=1:10
        sc(j,i)=xt((i-1)*5+j);
    end;
end;	
disp([sc])			        %打印前50个数
mx=mean(xt)		        %求平均数并输出
dx=cov(xt)			        %求方差并输出
subplot(2,1,2);
p=hist(xt,10)				%将产生的2000个随机数分为10组
p=p/100; t=0.025:.1:.975;		%求概率密度
bar (t,p,1);
title('0-1均匀分布的白噪声直方图');
xlabel('x');
ylabel('f(x)');
[bx,i] = xcov(xt,10);			%τ取-10到10
Bx=bx/2000;				%求自相关函数Bx(τ)
figure 
subplot(2,1,1);
plot(i,Bx),title('自相关函数Bx分布图');
xlabel('τ');ylabel('Bx(τ)');
[tx,i] = xcorr(xt,10);			%τ取-10到10
Tx=tx/2000;
subplot(2,1,2);
plot(i, Tx),title('自相关函数Γx分布图');
xlabel('τ');ylabel('Γx(τ)');