clear;
clc;
X=rand(1,2000);%产生在[0,1]均匀分布的白色序列
X50=X(1:50);%提取序列的前50个数
b=reshape(X50,5,10);
fprintf('%s\n','白色序列的前50个数:');
b1=b';
disp(b1);
EX=mean(X);%求取平均值
fprintf('EX=%g\n',EX);
DX=var(X);%求取方差
fprintf('DX=%g\n',DX);
subplot(2,1,1);
hist(X,10);%画出分布检验的直方图
title('均匀分布的白色序列的分布检验图');
y=zeros(1,21);
%以下求取白色序列的相关函数过程
for m=-10:10
absm=abs(m);%取数m的绝对值
sum=0;
for n=1:2000-absm
sum=sum+(X(n+absm)-EX)*(X(n)-EX);
end
y(m+11)=sum/2000;
end
m=(-10:10);
subplot(2,1,2);
plot(m,y);%画出白色序列的相关检验图
title('均匀分布的白色序列的相关检验图');