Python快速变换傅里叶(FFT)过程实例分析

您所在的位置:网站首页 交流电与直流电的根本区别是什么 Python快速变换傅里叶(FFT)过程实例分析

Python快速变换傅里叶(FFT)过程实例分析

#Python快速变换傅里叶(FFT)过程实例分析| 来源: 网络整理| 查看: 265

这里做一下记录,关于FFT就不做介绍了,直接贴上代码,有详细注释的了:import numpy as npfrom scipy.fftpack import fft,ifftimport matplotlib.pyplot as pltimport seaborn#采样点选择1400个,因为设置的信号频率分量最高为600赫兹,根据采样定理知采样频率要大于信号频率2倍,所以这里设置采样频率为1400赫兹(即一秒内有1400个采样点,一样意思的)x=np.linspace(0,1,1400) #设置需要采样的信号,频率分量有180,390和600y=7*np.sin(2*np.pi*180*x) + 2.8*np.sin(2*np.pi*390*x)+5.1*np.sin(2*np.pi*600*x)yy=fft(y) #快速傅里叶变换yreal = yy.real # 获取实数部分yimag = yy.imag # 获取虚数部分yf=abs(fft(y)) # 取绝对值yf1=abs(fft(y))/len(x) #归一化处理yf2 = yf1[range(int(len(x)/2))] #由于对称性,只取一半区间xf = np.arange(len(y)) # 频率xf1 = xfxf2 = xf[range(int(len(x)/2))] #取一半区间plt.subplot(221)plt.plot(x[0:50],y[0:50]) plt.title('Original wave')plt.subplot(222)plt.plot(xf,yf,'r')plt.title('FFT of Mixed wave(two sides frequency range)',fontsize=7,color='#7A378B') #注意这里的颜色可以查询颜色代码表plt.subplot(223)plt.plot(xf1,yf1,'g')plt.title('FFT of Mixed wave(normalization)',fontsize=9,color='r')plt.subplot(224)plt.plot(xf2,yf2,'b')plt.title('FFT of Mixed wave)',fontsize=10,color='#F08080')plt.show() 结果:

再添加一个简单的例子

# -*- coding: utf-8 -*-import matplotlib.pyplot as pltimport numpy as npimport seabornFs = 150.0; # sampling rate采样率Ts = 1.0/Fs; # sampling interval 采样区间t = np.arange(0,1,Ts) # time vector,这里Ts也是步长ff = 25; # frequency of the signaly = np.sin(2*np.pi*ff*t)n = len(y) # length of the signalk = np.arange(n)T = n/Fsfrq = k/T # two sides frequency rangefrq1 = frq[range(int(n/2))] # one side frequency rangeYY = np.fft.fft(y) # 未归一化Y = np.fft.fft(y)/n # fft computing and normalization 归一化Y1 = Y[range(int(n/2))]fig, ax = plt.subplots(4, 1)ax[0].plot(t,y)ax[0].set_xlabel('Time')ax[0].set_ylabel('Amplitude')ax[1].plot(frq,abs(YY),'r') # plotting the spectrumax[1].set_xlabel('Freq (Hz)')ax[1].set_ylabel('|Y(freq)|')ax[2].plot(frq,abs(Y),'G') # plotting the spectrumax[2].set_xlabel('Freq (Hz)')ax[2].set_ylabel('|Y(freq)|')ax[3].plot(frq1,abs(Y1),'B') # plotting the spectrumax[3].set_xlabel('Freq (Hz)')ax[3].set_ylabel('|Y(freq)|')plt.show()

原文标题:Python实现快速傅里叶变换(FFT)

文章出处:【微信号:eet-china,微信公众号:电子工程专辑】欢迎添加关注!文章转载请注明出处。

收藏 人收藏

扫一扫,分享给好友

复制链接分享 评论

发布评论请先 登录

相关推荐

pip 的超全使用指南 1. 查询软件包查询当前环境安装的所有软件包 [code]$ pip list[/code]查询 pypi 上含有某名字的包 [code]$ pip search pkg[/co... 发表于 03-11 16:03 • 3659次 阅读 一文了解pip的超全使用指南 由于默认情况下,wheel 包的平台是运行 pip download 命令 的平台,所以可能出现平台.... 的头像 python爬虫知识分享 发表于 03-11 16:03 • 174次 阅读 python常规包与命名空间包 python常规包与命名空间包1. 常规包在 Python 3.3 之前或者说 Python 2 中,一个包想要被导入使用,那么该包内必须要有 __init... 发表于 03-11 15:46 • 2060次 阅读 详解python常规包与命名空间包 python常规包与命名空间包 1. 常规包 在 Python 3.3 之前或者说 Python 2.... 的头像 python爬虫知识分享 发表于 03-11 15:46 • 181次 阅读 python花式导包的八种方法 python花式导包的八种方法1. 直接 import人尽皆知的方法,直接导入即可 [code]>>> import os >>> os.getcwd() '/home/xxx'[/c... 发表于 03-10 16:51 • 756次 阅读 python花式导包的八种方法 python花式导包的八种方法 1. 直接 import 人尽皆知的方法,直接导入即可 import.... 的头像 python爬虫知识分享 发表于 03-10 16:48 • 165次 阅读 python安装第三方包的八种方法 python安装第三方包的八种方法1. 使用 easy_installeasy_install 这应该是最古老的包安装方式了,目前基本没有人使用了。下面是 ... 发表于 03-10 16:29 • 1027次 阅读 python安装第三方包的八种方法 python安装第三方包的八种方法 1. 使用 easy_install easy_install .... 的头像 python爬虫知识分享 发表于 03-10 16:27 • 138次 阅读 今天咱们用Python整一个超级玛丽游戏 | 附带源码 开发环境Python版本:3.7.8相关模块:requests模块;tqdm模块;pyfreeproxy模块;pyecharts模块;以及一些python自带的模块。... 发表于 03-10 14:26 • 743次 阅读 今天咱们用Python整一个超级玛丽游戏 | 附带源码 如何在RK3288上安装Opencv? 如何在RK3288上安装Opencv? 发表于 03-10 06:49 • 181次 阅读 如何才能在RV1109平台上编写基于USB Camera的demo程序? 如何才能在RV1109平台上编写基于USB Camera的demo程序?... 发表于 03-10 06:15 • 82次 阅读 什么是python包、模块和库? 1. 模块以 .py 为后缀的文件,我们称之为 模块,英文名 Module。 模块让你能够有逻辑地组织你的 Python 代码段,把相关的代码... 发表于 03-09 16:48 • 461次 阅读 python包、模块和库是什么 1. 模块 以 .py 为后缀的文件,我们称之为 模块,英文名 Module。 模块让你能够有逻辑地.... 的头像 python爬虫知识分享 发表于 03-09 16:47 • 199次 阅读 python类的多态和类的property属性 python类的多态多态,是指在同一类型下的不同形态。 比如下面这段代码 [code]class People:     def speak(self):... 发表于 03-09 16:38 • 742次 阅读 python类的多态和类的property属性 python类的多态 多态,是指在同一类型下的不同形态。 比如下面这段代码 class People.... 的头像 python爬虫知识分享 发表于 03-09 16:37 • 166次 阅读 用Python学习科学编程 用Python学习科学编程,Python经典教材。 发表于 03-09 15:00 • 26次 阅读 怎样在debian下安装python3+pymssql呢 怎样在debian下安装python3+pymssql呢?有哪些安装步骤?... 发表于 03-09 06:32 • 79次 阅读 python类的继承详解 python类的继承 类的继承,跟人类繁衍的关系相似。 被继承的类称为基类(也叫做父类),继承而得的.... 的头像 python爬虫知识分享 发表于 03-08 16:40 • 333次 阅读 python私有变量和私有方法 python私有变量和私有方法 1. 下划线妙用 在 Python 中,下划线可是非常推荐使用的符号.... 的头像 python爬虫知识分享 发表于 03-08 16:30 • 384次 阅读 python静态方法与类方法 python静态方法与类方法 1. 写法上的差异 类的方法可以分为: 静态方法:有 staticme.... 的头像 python爬虫知识分享 发表于 03-07 16:56 • 405次 阅读 python类的理解与使用 python类的理解与使用 1. 通俗理解类 类(英文名 class),是具有相同特性(属性)和行为.... 的头像 python爬虫知识分享 发表于 03-07 16:51 • 386次 阅读 python如何捕获异常和主动抛出异常 python如何主动抛出异常和捕获异常 1. 如何抛出异常? 异常的产生有两种来源: 一种是程序自动.... 的头像 python爬虫知识分享 发表于 03-04 17:09 • 660次 阅读 Python中有哪些常见的错误和异常 python常见异常类型 在程序运行过程中,总会遇到各种各样的问题和错误。 有些错误是我们编写代码时.... 的头像 python爬虫知识分享 发表于 03-04 16:58 • 665次 阅读 python变量的作用域 python变量的作用域 1. 作用域 Python的作用域可以分为四种: L (Local) 局部.... 的头像 python爬虫知识分享 发表于 03-03 16:50 • 400次 阅读 python偏函数和泛型函数详解 python偏函数 假如一个函数定义了多个位置参数,那你每次调用时,都需要把这些个参数一个一个地传递.... 的头像 python爬虫知识分享 发表于 03-03 16:43 • 413次 阅读 python高阶函数详解 python高阶函数 1. map 函数 map 函数,它接收两个参数,第一个参数是一个函数对象(当.... 的头像 python爬虫知识分享 发表于 03-02 16:47 • 175次 阅读 python高阶函数详解 python匿名函数的使用 python匿名函数的使用 匿名函数(英语:anonymous function)是指一类无需定义标.... 的头像 python爬虫知识分享 发表于 03-02 16:42 • 170次 阅读 11个案例讲解python函数参数 函数,在定义的时候,可以有参数的,也可以没有参数。 的头像 python爬虫知识分享 发表于 03-01 16:39 • 241次 阅读 详解python普通函数创建与调用 函数是一种仅在调用时运行的代码块。您可以将数据(称为参数)传递到函数中,然后由函数可以把数据作为结果.... 的头像 python爬虫知识分享 发表于 03-01 16:32 • 238次 阅读 python推导式是什么 python推导式 推导式(英文名:comprehensions),也叫解析式,是Python的一种.... 的头像 python爬虫知识分享 发表于 02-28 17:13 • 193次 阅读 python while循环详解 python while循环 while 语句用于循环执行程序,即在某条件下,循环执行某段程序,以处.... 的头像 python爬虫知识分享 发表于 02-28 16:39 • 171次 阅读 《深入Python3》中文版pdf 《深入Python3》中文版pdf 发表于 02-28 09:45 • 35次 阅读 《Dive Into Python》中文翻译版.pdf 《Dive Into Python》中文翻译版.pdf 发表于 02-25 17:22 • 29次 阅读 python for循环的案例说明 python for循环 for 循环可以遍历任何序列的项目,如一个列表或者一个字符串。 它的基本语.... 的头像 python爬虫知识分享 发表于 02-25 16:27 • 215次 阅读 python判断语句的详细说明 python判断语句:if 1. 简单小例子 如果满足条件 A,则执行代码块 a,否则执行代码块 b.... 的头像 python爬虫知识分享 发表于 02-25 16:22 • 221次 阅读 python生成器是什么 python生成器 1. 什么是生成器? 生成器(英文名 Generator ),是一个可以像迭代器.... 的头像 python爬虫知识分享 发表于 02-24 15:53 • 235次 阅读 python迭代器详解 python迭代器 1. 可迭代对象 可以利用 for 循环的对象,都叫可迭代对象。 列表、元组、字.... 的头像 python爬虫知识分享 发表于 02-24 15:42 • 203次 阅读 初学者学Python必看的几个练手小项目 Python是一种面向对象的解释型编程语言,源代码与解释器CPython遵守GPL协议,Python.... 的头像 叶枫架构师 发表于 02-23 17:06 • 291次 阅读 4个Python经典项目实战,练手必备! Python是一种极具可读性和通用性的编程语言。Python这个名字的灵感来自于英国喜剧团体Mont.... 的头像 叶枫架构师 发表于 02-23 17:06 • 273次 阅读 python集合是什么 python集合 集合(英文名 set),它是一个无序的不重复元素序列。 这里面有两个重点: 无序,.... 的头像 python爬虫知识分享 发表于 02-23 17:01 • 258次 阅读 python字典是什么 python字典 字典(英文名 dict),它是由一系列的键值(key-value)对组合而成的数据.... 的头像 python爬虫知识分享 发表于 02-23 16:54 • 237次 阅读 python元组的详细介绍 python元组 元组(英文名 tuple),和列表非常的相似,它也是由一系列元素按顺序进行排列而成.... 的头像 python爬虫知识分享 发表于 02-22 15:19 • 159次 阅读 python列表的详细介绍 python列表 列表(英文名 list),是由一系列元素按顺序进行排列而成的容器。 这里面有两个重.... 的头像 python爬虫知识分享 发表于 02-22 15:14 • 159次 阅读 python cookbook(第三版)免费下载 python cookbook(第三版)免费下载。 发表于 02-22 13:57 • 47次 阅读 python运算符是什么 python运算符 0. 什么是运算符? 本章节主要说明Python的运算符。举个简单的例子 4 +.... 的头像 python爬虫知识分享 发表于 02-21 16:44 • 203次 阅读 python字符串格式化 python字符串格式化 格式化输出,主要有三种方式 使用 % 进行格式化 使用 format 函数.... 的头像 python爬虫知识分享 发表于 02-21 16:28 • 178次 阅读 python字符串格式化 python输入与输出详解 python输入与输出 无论是从我们一开始的“hello world”,还是前面章节的里各种例子,基.... 的头像 python爬虫知识分享 发表于 02-21 16:18 • 224次 阅读 python布尔值是什么 python布尔值 1. 什么是布尔值 但在Python语言中,布尔类型只有两个值 True :表示.... 的头像 python爬虫知识分享 发表于 02-21 16:10 • 190次 阅读 python整数与浮点数分解 2.3 python整数与浮点数 Python 支持的数字类型有三种:整数、浮点数和复数。 1. 整.... 的头像 python爬虫知识分享 发表于 02-18 09:09 • 293次 阅读 2.2 python字符串类型 2.2 python字符串类型 1. 如何定义字符串? 字符串是Python中最常用的数据类型之一。.... 的头像 python爬虫知识分享 发表于 02-17 17:08 • 183次 阅读 深入了解python常量与变量   Python 的变量和常量不需要事先声明类型,这是根据Python的动态语言特性而来。 的头像 python爬虫知识分享 发表于 02-16 18:22 • 340次 阅读 深入了解python常量与变量 解析两种运行Python程序方法 首先打开你的终端,直接输入 python3 回车,然后输入 print("Hello, World".... 的头像 python爬虫知识分享 发表于 02-16 18:17 • 459次 阅读 解析两种运行Python程序方法 1.2 Python 开发环境的搭建 1.2 Python 开发环境的搭建 理论上只需要你安装了 CPython 解释器后,就可以开始写 .... 的头像 python爬虫知识分享 发表于 02-15 17:16 • 316次 阅读 1.2 Python 开发环境的搭建 快速安装 Python 解释器 1.1【环境】快速安装 Python 解释器 Python 是一门解释性脚本语言,因此要想让你编写的.... 的头像 python爬虫知识分享 发表于 02-15 16:57 • 542次 阅读 快速安装 Python 解释器 《Python编程入门》.pdf 《Python编程入门》.pdf 发表于 02-11 16:03 • 124次 阅读 python获取历届春晚数据 不知道今年的春晚大家看了吗?小编表示没有看,而且是已经很多年没有完整的看过春晚了,关于看春晚的热情都.... 的头像 yiniuyun 发表于 02-07 17:52 • 314次 阅读 Python语言主要有哪些特点及Python图形界面框架推荐 Python由荷兰数学和计算机科学研究学会的吉多·范罗苏姆 于1990 年代初设计,作为一门叫做AB.... 的头像 Les 发表于 01-30 16:44 • 259次 阅读 PikaScript入选2021年度Gitee最有价值开源项目 PikaScript是一个完全重写的超轻量级python引擎,具有完整的解释器,字节码和虚拟机架构,.... 的头像 RTThread物联网操作系统 发表于 01-26 15:59 • 591次 阅读 入坑ESP32,首次MicroPython项目移植 原先一直在ARM架构下编写Cpython,先来说说ARM架构优缺点:1、系统资源足够,虽然pytho.... 发表于 01-18 11:24 • 68次 阅读 入坑ESP32,首次MicroPython项目移植 树莓派3B+使用STM32扩展板实现按键驱动 前几天入了个树莓派3B+,之前用过树莓派zero W觉得很方便,但是外围接口少了点。准备用3B+做个.... 发表于 01-18 11:05 • 61次 阅读 树莓派3B+使用STM32扩展板实现按键驱动


【本文地址】


今日新闻


推荐新闻


    CopyRight 2018-2019 办公设备维修网 版权所有 豫ICP备15022753号-3