0
  • 聊天消息
  • 系统消息
  • 评论与回复
登录后你可以
  • 下载海量资料
  • 学习在线课程
  • 观看技术视频
  • 写文章/发帖/加入社区
会员中心
创作中心

完善资料让更多小伙伴认识你,还能领取20积分哦,立即完善>

3天内不再提示

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

0BFC_eet_china 2017-11-29 15:03 次阅读

这里做一下记录,关于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()

声明:本文内容及配图由入驻作者撰写或者入驻合作网站授权转载。文章观点仅代表作者本人,不代表电子发烧友网立场。文章及其配图仅供工程师学习之用,如有内容侵权或者其他违规问题,请联系本站处理。 举报投诉
  • FFT
    FFT
    +关注

    关注

    15

    文章

    438

    浏览量

    59769
  • python
    +关注

    关注

    56

    文章

    4812

    浏览量

    85236

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

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

收藏 人收藏

    相关推荐

    变换

    LABVIEW是怎样进行非周期波形的变换的,各路大神飘过的,求。在线急等。
    发表于 10-21 14:59

    快速傅立叶变换FFT)的Nios II实现

    快速傅立叶变换FFT)的Nios II实现 随着数字电子技术的发展,数字信号处理的理论和技术广泛地应用于通讯、语音处理、计算机和多媒体等领域。快速
    发表于 02-09 09:38 81次下载

    分数(p>1)级变换的无透镜光学实现

    摘要用球面波照明物体的自由空间菲涅耳衍射3实现分数(p1)级变换,提供了分数
    发表于 02-23 09:35 29次下载
    分数(p>1)级<b class='flag-5'>傅</b><b class='flag-5'>里</b><b class='flag-5'>叶</b>逆<b class='flag-5'>变换</b>的无透镜光学实现

    小波与分析基础

    向读者展示分析和小波的许多基础知识以及在信号分析方面的应用。全书分为8章和3个附录,第0章是学习第1章至第7章的准备知识,即内积空间;
    发表于 07-14 11:25 0次下载
    小波与<b class='flag-5'>傅</b><b class='flag-5'>里</b><b class='flag-5'>叶</b><b class='flag-5'>分析</b>基础

    小波与分析基础_中文版_

    电子发烧友网站提供《小波与分析基础_中文版_.txt》资料免费下载
    发表于 07-07 13:32 0次下载

    小波与分析基础

    小波与分析基础 有需要的朋友下来看看
    发表于 12-30 15:33 0次下载

    级数和傅里叶变换的关系

    级数对周期性现象做数学上的分析傅里叶变换可以看作
    发表于 11-24 14:32 4w次阅读
    <b class='flag-5'>傅</b><b class='flag-5'>里</b><b class='flag-5'>叶</b>级数和傅里叶<b class='flag-5'>变换</b>的关系

    周期信号的频谱分析——级数

    周期信号的频谱分析——级数
    发表于 12-06 14:27 0次下载

    级数电路分析级数表示简介

    了解级数在电路分析级数方程中的重要性,
    的头像 发表于 01-27 14:11 1271次阅读
    <b class='flag-5'>傅</b><b class='flag-5'>里</b><b class='flag-5'>叶</b>级数电路<b class='flag-5'>分析</b> — <b class='flag-5'>傅</b><b class='flag-5'>里</b><b class='flag-5'>叶</b>级数表示简介

    傅里叶变换级数的关系

    傅里叶变换级数的关系  傅里叶变换
    的头像 发表于 09-07 16:39 4591次阅读

    傅里叶变换变换的关系

    傅里叶变换变换的关系  傅里叶变换
    的头像 发表于 09-07 16:43 8204次阅读

    如何由傅里叶变换推出变换

    如何由傅里叶变换推出变换  傅里叶变换
    的头像 发表于 09-07 17:04 2491次阅读

    电子硬件中的变化

    是法国数学家,物理学家。曾推导出著名的热传导方程,并在求解的过程中发现了解函数可以由三角函数构成的级数形式表示。 从而提出了任一函数都可以变成三角函数的无穷级数。傅里叶
    的头像 发表于 11-06 11:03 794次阅读
    电子硬件中的<b class='flag-5'>傅</b><b class='flag-5'>里</b><b class='flag-5'>叶</b>变化

    连续时间信号的分析

    电子发烧友网站提供《连续时间信号的分析.pdf》资料免费下载
    发表于 11-18 15:25 0次下载

    快速傅里叶变换-FFT分析仪基础知识

    FFT频谱分析仪的概念是围绕快速傅里叶变换建立的,该变换基于约瑟夫·
    发表于 01-16 14:26 1461次阅读