normalwhitenoise.txt

来自「产生 正态白噪声序列 (1) 打印出前50个数 (2) 分布检验 (3)」· 文本 代码 · 共 29 行

产生 正态白噪声序列 
(1) 打印出前50个数;       (2) 分布检验;
(3) 均值检验;              (4) 方差检验;
(5) 计算相关函数 Bx(i),i=0,±1,±2,…, ±10。
B(m)=1/1000 
MATLAB编程如下：
clear all;
for n=1:2000
xt(n)=normrnd(0,1);;		%产生2000个N(0,1)正态分布的白噪声
end
figure
for i=1:5
    for j=1:10
        sc(j,i)=xt((i-1)*5+j);
    end;
end;disp([sc]);
mx=mean(xt)					%求平均数并输出
dx=cov(xt)						%求方差并输出
subplot(2,1,1);
p=hist(xt,20)	;				%将产生的2000个随机数分为20组
p=p/100; t=-2.85:0.3:2.85;			%求概率密度
bar (t,p,1) ,title('N(0,1)正态分布的白噪声直方图');
xlabel('x');ylabel('f(x)');
[tx,i] = xcov(xt,10);				%τ取-10到10
Tx=tx/2000;						%求自相关函数Γx(τ)
subplot(2,1,2)
plot(i,Tx),title('自相关函数Γx(τ)分布图');
xlabel('τ');ylabel('Γx(τ)');