matlab noise level,Matlab生成白噪声

所谓白噪声，是指没有任何时迟相关时间序列，它只一个参数，就是方差，代表噪声的强度，也即迟后为0的滞后自协方差，而其余滞后自协方差统统为零。标准白噪声(standard

Gaussian white noise0，1)标准正态分布的白噪声。

我今天想要做的，就是生成100个150*1的列向量，每个列向量都是一个标准白噪声样本。在MATLAB中有两种生成方法：

noise=wgn(150,100,1);

noise1=randn(150,100);

从randn、wgn以及其后awgn的源码看，其实是wgn调用了randn，而awgn又调用了wgn。

虽然randn的help文件里面明确写着：The randn function generates arrays of

random numbers whose elements are normally distributed with mean 0,

variance σ2=1, and standard deviationσ=1.但实际情况并不理想。

考虑到150*1这个序列长度有些短，我实验了10000*1的序列，即：noise=wgnksdensity(10000,1,1); 及

noise1=randn(10000,1);

可以用hist或是ksdensity做出noise和noise1的概率分布图，附在了左面图区。从图中可以看出来，虽然分布形式很正态，但是mean(noise)及std(noise)都不是想像中的结果。wgn也就罢了，但是randn毕竟没有能做到它所承诺做到的事情。

这样之后就能得到：“mean(noise)= -3.8858e-017; 及 std(noise)=1.0000;”

这样比较漂亮的结果，用着也放心。

1、随机问题都有取样误差的问题，只有当样本趋于无穷大，样本才正确，而任何有限样本肯定存在误差。

2、randn是通过非线性叠代产生的伪随机数，并非真正的随机数。任何基于有限自动机理论上面的随机数生成函数都是伪随机数。