|
|
发表于 2009-9-28 15:40:28
|
显示全部楼层
x=randn(1,10000);
x=x/std(x);
x=x-mean(x);
a=0; %均值
b=1; %方差
y=a+sqrt(b)*x; %产生均值零,方差为1的高斯噪声
Y=fft(y);
subplot(121),plot(y); title('白噪声y') %此处可以看到高斯白噪声的频谱还是白噪声
subplot(122),plot(abs(Y)); title('y的幅值谱') %高斯函数的傅立叶变换还是高斯函数
S=mean(y.^3)./(mean(y.^2)).^1.5;%此处计算偏斜度
K=mean(y.^4)./(mean(y.^2)).^2-3;%此处计算峭度 |
|