一个数字信号处理的简单问题

[jun lin]

dear all forks:
最近发了很多的邮件，在询问如何做数字信号处理的，但苦于没有实际的数据。
现在终于被我整到一个实际的数据啦！附件是图片和数据（已经整理成一个data的tuple。）
现在我可以给它叠加，fft，截取一个部分，玩数据还是满有意思的。
但是最后还是要做事。图片上可以看出来，这个东西有好几个波峰，需要给识别出来。
现在我的方案：
1.numpy.fft后，看到有几个频率的幅度很大，考虑滤波。但是怎么滤呢？
现在在看:http://pyscin.appspot.com/html/frequency_process.html中。。。
2.现在暂时用叠加的方法来整出了一个平滑的曲线了（可惜第一个特征峰被抵消掉了。。）
然后就是曲线求峰。现在用的方法是求numpy.diff，然后判断0的位置。。
好吧，差分曲线波动还是挺大的，看到波峰位置的幅度很大，但是不知道怎么编程取得。。

有哪位是专家？给个学习链接？

[jun lin]

多么漂亮的FFT。。。

2009/12/16 jun lin <linjun...@gmail.com>

[jun lin]

我可怜的数学知识愣是没有求到最大的那个峰。。。

2009/12/16 jun lin <linjun...@gmail.com>

[abracadabra]

你来错地方了吧？
http://groups.google.com/group/comp.dsp/topics

你抽样频率，希望得到的频段等参数都不清楚。可以用Matlab，玩不同的滤波器。信号处理的方法多得很。

> ssss.PNG
> 40KViewDownload
>
> sss.py
> 37KViewDownload

[Shuo Chen]

你对信号处理了解多少? 比如咒语"时域卷积等于频域相乘"听着顺耳吗？
要滤波,最简单的是FIR, 就是找一个信号, 让它的频域是个低通 (截至频率取决于你想滤掉几次谐波), 然后把这个信号和你的数据卷积, 就得到
了低通滤波后的信号.
FIR 的设计是个优化问题, 最简单的就是窗函数法, 你先用教科书上的办法试一试, 看看效果如何.

> ssss.PNG
> 40KViewDownload
>
> sss.py
> 37KViewDownload

[Tinyfool]

同意，感觉楼主需要了解的是信号处理的基础知识，先去找本书学学，然后再写这个程序吧。

2009/12/17 Shuo Chen <gian...@gmail.com>

--
Tinyfool的开发日记 http://www.tinydust.net/dev/
代码中国网 http://www.codechina.org
myTwitter: http://twitter.com/tinyfool

[Shuo Chen]

或者先试试提取信号包络？ signal envelope

On Dec 17, 12:34 pm, Tinyfool <tinyf...@gmail.com> wrote:
> 同意，感觉楼主需要了解的是信号处理的基础知识，先去找本书学学，然后再写这个程序吧。
>

> 2009/12/17 Shuo Chen <giantc...@gmail.com>

[jun lin]

OK,还是老实回去恳信号与系统，以及数字信号处理。
不知道2个月能否看出点眉头来。。。

2009/12/17 Shuo Chen <gian...@gmail.com>

[ZHOU Guibin]

推荐你看看Richard G. Lyons 的 “Understanding Digital Signal Processing (2nd Edition)”  ,最好的DSP入门书。Amazon五星推荐

2009/12/17 jun lin <linjun...@gmail.com>