网站首页 > 技术文章 正文
用FFT和iFFT实现滤波:
图1对应代码链接:https://gitee.com/anjiang2020_admin/fft/blob/master/fft.py
ys 为一个由频率为180,390,600的正弦波叠加而成的信号
ys_390是ys中频率为390的那个信号。用以下代码对ys进行了滤波,只剩下390频率的信号。图1中的ys filterd by H就是滤波后的信号,其与ys_390基本一样。
fft(ys)为ys的频谱。可以看到,频率在180,390,600时是峰值。
filtered signal's fft: 是 ys filtered by H的频谱,可以看到,因为被H滤波了,只剩下f=390的频率了。
the time signal of H:H的反傅立叶变换。这个信号要和ys做卷积,同样会产生滤波效果。
filter H:就是滤波器了,这是个带通滤波器,只剩下频率在390附近的信号,其它频率的信号过滤掉。
用卷积实现滤波:
图2中的convH就是用卷积实现滤波的效果,可以看起,其曲线与ys_390一致。
卷积的代码如下:
引入scipy.signal里的convolve函数
conv(ys,filter_h_org)实现了卷积滤波。这里的filter_h_org是由ifft(H)得到的。将ifft(H)的后半部分平移到前半部分即可。
总结:可以看到,用fft和ifft实现的滤波结果和卷积实现的滤波结果一致。
所以说:卷积神经网络里的卷积操作,其实是对输入信号的滤波。
关于卷积FFT与 DFT以及DFT的源码实现,请参考:
赵明明:[信号与处理]傅立叶变换中的DFT实现以及与scipy.fft的对比:https://www.toutiao.com/i7047374338911224357/
参考:理解Scipy卷积:https://link.zhihu.com/?target=https%3A//blog.csdn.net/weixin_39715997/article/details/110968470
猜你喜欢
- 2024-11-03 正点原子开拓者FPGA开发板资料连载第五十二章 低通滤波器实验
- 2024-11-03 图像频域及滤波处理 图像频域滤波变换的实验报告
- 2024-11-03 低通滤波和高通滤波的截止频率设定思路
- 2024-11-03 「科唛小课堂」如何解决录制的爆音、喷麦声?
- 2024-11-03 「正点原子NANO STM32开发板资料连载」第三十二章 DSP 测试实验
- 2024-11-03 最好的频谱分析仪基础知识 频谱分析仪的作用与性能
- 2024-11-03 OER的几何重构 几何分解重构图
- 2024-11-03 关于分频音箱问题的非专业理想解决方案
- 2024-11-03 STM32单片机从零开始使用教程(八) FIR滤波器
- 2024-11-03 深入剖析C++音频信号处理:ASL库的高级应用与实践
你 发表评论:
欢迎- 最近发表
-
- 在 Spring Boot 项目中使用 activiti
- 开箱即用-activiti流程引擎(active 流程引擎)
- 在springBoot项目中整合使用activiti
- activiti中的网关是干什么的?(activiti包含网关)
- SpringBoot集成工作流Activiti(完整源码和配套文档)
- Activiti工作流介绍及使用(activiti工作流会签)
- SpringBoot集成工作流Activiti(实际项目演示)
- activiti工作流引擎(activiti工作流引擎怎么用)
- 工作流Activiti初体验及在数据库中生成的表
- Activiti工作流浅析(activiti6.0工作流引擎深度解析)
- 标签列表
-
- oraclesql优化 (66)
- 类的加载机制 (75)
- feignclient (62)
- 一致性hash算法 (71)
- dockfile (66)
- 锁机制 (57)
- javaresponse (60)
- 查看hive版本 (59)
- phpworkerman (57)
- spark算子 (58)
- vue双向绑定的原理 (68)
- springbootget请求 (58)
- docker网络三种模式 (67)
- spring控制反转 (71)
- data:image/jpeg (69)
- base64 (69)
- java分页 (64)
- kibanadocker (60)
- qabstracttablemodel (62)
- java生成pdf文件 (69)
- deletelater (62)
- com.aspose.words (58)
- android.mk (62)
- qopengl (73)
- epoch_millis (61)
本文暂时没有评论,来添加一个吧(●'◡'●)