如何向一个数据序列中添加已知频率范围的高频噪声信号？

caichengtao

如何向一个数据序列中添加已知频率范围的高频噪声信号？
和添加高斯噪声和随机噪声有什么区别啊?

没有人回答啊，版主给个建议，是不是根本就不可行啊？

[ 本帖最后由 ChaChing 于 2010-5-1 01:02 编辑 ]

eight

原帖由 caichengtao 于 2007-9-11 08:37 发表

                               
登录/注册后可看大图

没有人回答啊，版主给个建议，是不是根本就不可行啊？

注意请勿催贴。如果大伙懂，会给你回帖的，不懂的话，就不好乱发帖说：“这个我不懂”，是吧？论坛上面的讨论要给点耐性，而且在等待的过程中你也可以自己测试一下

zhlong

楼主可以试试直接用滤波器对randn产生的高斯白噪声进行滤波，得到你想要的频率段。
可以参考《全相位滤波白噪声产生有色噪声及其功率谱估计》，全相位滤波在信号处理板块有讨论。

caichengtao

感谢啊！感谢论坛，衷心接受批评！

zhlong

原帖由 花如月 于 2007-7-26 08:50 发表

                               
登录/注册后可看大图

clear;
x=randn(1,10000);
x=x/std(x);
x=x-mean(x);
a=0; %均值
b=1; %方差
y=a+sqrt(b)*x;  %产生均值零，方差为1的高斯噪声

>> K
K =
    0.0251
>> S
S =
  -0.0126

将白噪声通过下面的带通滤波器：

得到的“色噪声”信号及其频谱：

高斯色噪声峭度和偏度：
>> K
K =
    0.0676
>> S
S =
  9.2605e-005

下图表示白噪声、滤波得到的色噪声的幅值分布以及标准正态分布：

标准正态分布可以通过下式得到：
t=-5:0.01:5;
nnormal=1/sqrt(2*pi)*exp((-t.^2)/2);
figure;plot(t,nnormal)
下面的语句可以画出信号幅值分布：
[f,xi] = ksdensity(sig);
plot(xi,f); % 画概率密度曲线

花如月

讲的很详细，龙真是好样的

zhlong

呵呵，以前没有用到过有色噪声，所以找资料看了一下。

fuheshan

非常感谢
写的很详细

xiaoxiong702

看了一下zhlong君还有其他人整理相关噪声的的帖子。这儿有个问题想请教，你的白噪声频谱是怎么得到的？我用产生的白噪声后求自相关函数再作fft变换，画出来的图和你的有很大差异，不知道为什么。

[ 本帖最后由 xiaoxiong702 于 2009-5-15 22:30 编辑 ]

zengkeyeah

太详细了，佩服

200810100659

讲得不错。学习