python傅里叶逆变换_Python傅里叶逆变换
问题
我试图从一组从音频文件中获取的数据中删除一个频率。在
为了简化我的问题,我创建了下面的代码,它创建了一组wave并将它们合并到一个复杂的wave中。然后求出复波的傅里叶变换并求逆。在
我希望看到原始波形的结果,因为应该没有数据丢失,但是我收到了一个非常不同的波形。在
代码:import numpy as np
import matplotlib.pyplot as plt
import random
#Get plots
fig, c1 = plt.subplots()
c2 = c1.twinx()
fs = 100 # sample rate
f_list = [5,10,15,20,100] # the frequency of the signal
x = np.arange(fs) # the points on the x axis for plotting
# compute the value (amplitude) of the sin wave for each sample
wave = []
for f in f_list:
wave.append(list(np.sin(2*np.pi*f * (x/fs))))
#Adds the sine waves together into a single complex wave
wave4 = []
for i in range(len(wave[0])):
data = 0
for ii in range(len(wave)):
data += wave[ii][i]
wave4.append(data)
#Get frequencies from complex wave
fft = np.fft.rfft(wave4)
fft = np.abs(fft)
#Note: Here I will add some code to remove specific frequencies
#Get complex wave from frequencies
waveV2 = np.fft.irfft(fft)
#Plot the complex waves, should be the same
c1.plot(wave4, color="orange")
c1.plot(waveV2)
plt.show()
结果:(橙色为产生波,蓝色为原始波)
预期结果:
蓝色和橙色线条(原始和新创建的波浪)的值应该完全相同
python傅里叶逆变换_Python傅里叶逆变换相关推荐
- Python学习教程(Python学习视频_Python学些路线):Day05 总结和练习
Python学习教程(Python学习视频_Python学些路线):总结和练习 练习清单 寻找"水仙花数". 寻找"完美数". "百钱百鸡" ...
- 浅谈傅里叶——5. 短时傅里叶的缺点与卷积的基本概念
在<浅谈傅里叶 4>里,我引入了卷积的概念来简要的解释了傅里叶的改进方法--短时傅里叶.它类似一种滑动的滤波器,只不过与我们熟知的与对空间滤波.频域滤波不同,这是一个类似时域滤波的滤波器, ...
- 视频教程-快速入门Python基础教程_Python基础知识大全-Python
快速入门Python基础教程_Python基础知识大全 十余年计算机技术领域从业经验,在中国电信.盛大游戏等多家五百强企业任职技术开发指导顾问,国内IT技术发展奠基人之一. 杨千锋 ¥99.00 立即 ...
- Python学习教程(Python学习视频_Python学习路线):Day04循环结构
Python学习教程(Python学习视频_Python学习路线):循环结构 循环结构的应用场景 如果在程序中我们需要重复的执行某条或某些指令,例如用程序控制机器人踢足球,如果机器人持球而且还没有进入 ...
- 视频教程-快速入门Python基础教程_Python基础进阶视频-Python
快速入门Python基础教程_Python基础进阶视频 十余年计算机技术领域从业经验,在中国电信.盛大游戏等多家五百强企业任职技术开发指导顾问,国内IT技术发展奠基人之一. 杨千锋 ¥199.00 立 ...
- Python学习教程(Python学习路线_Python基础学习教程_Python视频教程):初学者新手怎样快速入门Python
Python学习教程(Python学习路线_Python基础学习教程_Python视频教程):初学者新手怎样快速入门Python? 人生苦短,我用Python!!!短短几个字,现在在各大学习类平台随处 ...
- python傅里叶逆变换_C# 傅里叶变换 逆变换 调用MathNet包|简明python教程|python入门|python教程...
https://www.xin3721.com/eschool/pythonxin3721/ 上司给我讲采样率,我觉得自己有点钻牛角尖了,想来想去还是不明白. 因为说的是日语,サンプリング周波数,我就 ...
- 用python画校徽,傅里叶不会用,点动成线也可以
前几天想着用py画一个校徽,找了好多方法,都说要用傅里叶,但是那个看着不简单,就没有去研究.就自己整了整,用了 点动成线这个原理画了一下,没想到还成了,只不过绕着半圆打圈的字没有往上整,那个有点内容. ...
- python时间序列小波分析_python时间序列分析
什么是时间序列 时间序列简单的说就是各时间点上形成的数值序列,时间序列分析就是通过观察历史数据预测未来的值.在这里需要强调一点的是,时间序列分析并不是关于时间的回归,它主要是研究自身的变化规律的(这里 ...
- python傅里叶变换库_python的numpy库和cv2库实现图像傅里叶变换
码字不易,如果对您有所帮助,记着点赞哦! 一. 图像傅里叶变换原理: 对二维图像进行傅里叶变换用如下公式进行: 图像长M,高N.F(u,v)表示频域图像,f(x,y)表示时域图像.u的范围为[0,M- ...
最新文章
- javax消息队列_Java面试—消息队列
- python3 issubset 判断列表子集
- 对应生成树的基本回路_数据结构与算法——最小生成树
- 信息通信建设工程预算定额_通信建设工程概预算
- 当root.sh与ORA-15031相遇
- ScriptManager和ClientScriptManager
- sqlserver2014内存数据库特性介绍
- mysql 时间_MySQL 日期时间
- 怎么解锁blockinput锁定的键盘_“键盘云”--锐角云母公司(Acute angle)又一款矿渣神器...
- 微信公众号维护用服务器吗,微信公众号商城开发必须用服务器吗?
- 苹果app充值限制解除_2020还在充值退款?正规苹果app手游充值折扣来了!
- 鼠标右键没有新建选项怎么办
- android pc投屏,安卓投屏助手pc版
- this command with --force, or --legacy-peer-deps
- 什么软件可以测试网络的稳定性,手机怎么测试网络稳定性
- 粘包问题和struck解决
- 10.sprd背光流程
- 分子内电荷转移有哪些最新发表的毕业论文呢?
- 【html】屏蔽浏览器自带的输入提示
- JSP-java服务器端页面【学习笔记】