一，FFT解釋

FFT（Fast Fourier Transformation）是離散傅氏變換（DFT）的快速算法。即為快速傅氏變換。它是根據(jù)離散傅氏變換的奇、偶、虛、實(shí)等特性，對(duì)離散傅立葉變換的算法進(jìn)行改進(jìn)獲得的。

二、基于python畫出時(shí)域頻域波形

1.python代碼

依據(jù)快速傅里葉算法得到信號(hào)的頻域，將其各個(gè)頻率分量的幅值繪制成圖。為便于計(jì)算，我們將采樣頻率8000次近似設(shè)為8192次
正弦波表達(dá)式為： s(t) = 0.6 sin( 2π 50t ) 和s(t) = 0.6 sin( 2π 500t )

import numpy as np#導(dǎo)入一個(gè)數(shù)據(jù)處理模塊import matplotlib.pyplot as plt#導(dǎo)入一個(gè)繪圖模塊# 依據(jù)快速傅里葉算法得到信號(hào)的頻域 def test_fft():sampling_rate = 8192 # 采樣率fft_size = 8192 # FFT取樣長度t = np.arange(0, 8.192, 1.0 / sampling_rate)#np.arange(起點(diǎn)，終點(diǎn)，間隔)產(chǎn)生8.192s長的取樣時(shí)間x=0.6*np.sin(2*np.pi*500*t)+0.6*np.sin(2*np.pi*50*t)# 兩個(gè)正弦波疊加，500HZ和50HZ# N點(diǎn)FFT進(jìn)行精確頻譜分析的要求是N個(gè)取樣點(diǎn)包含整數(shù)個(gè)取樣對(duì)象的波形。# 因此N點(diǎn)FFT能夠完美計(jì)算頻譜對(duì)取樣對(duì)象的要求是n*Fs/N（n*采樣頻率/FFT長度），# 因此對(duì)8KHZ和512點(diǎn)而言，完美采樣對(duì)象的周期最小要求是8000/512=15.625HZ,# 所以156.25的n為10,234.375的n為15。xs = x[:fft_size]# 從波形數(shù)據(jù)中取樣fft_size個(gè)點(diǎn)進(jìn)行運(yùn)算xf = np.fft.rfft(xs) / fft_size # 返回fft_size/2+1 個(gè)頻率#利用np.fft.rfft()進(jìn)行FFT計(jì)算，rfft()是為了更方便對(duì)實(shí)數(shù)信號(hào)進(jìn)行變換，# 由公式可知 / fft_size為了正確顯示波形能量# rfft函數(shù)的返回值是N/2+1個(gè)復(fù)數(shù)，分別表示從0(Hz)到sampling_rate/2(Hz)的分。# 于是可以通過下面的np.linspace計(jì)算出返回值中每個(gè)下標(biāo)對(duì)應(yīng)的真正的頻率：freqs = np.linspace(0, sampling_rate*10, fft_size/2+1 ) # 表示頻率#freqs = np.linspace(0, sampling_rate/2 , fft_size/2 + 1) # 表示頻率xfp = 20 * np.log10(np.clip(np.abs(xf), 1e-20, 1e100))#xfp = np.abs(xf) * 2 # 代表信號(hào)的幅值，即振幅# 最后我們計(jì)算每個(gè)頻率分量的幅值，并通過 20*np.log10()將其轉(zhuǎn)換為以db單位的值。# 為了防止0幅值的成分造成log10無法計(jì)算，我們調(diào)用np.clip對(duì)xf的幅值進(jìn)行上下限處理plt.figure(figsize=(8, 4))plt.subplot(211)plt.plot(t[:fft_size], xs)plt.xlabel(u"時(shí)間(秒)", fontproperties='FangSong')plt.title(u"500Hz和50Hz的波形和頻譜", fontproperties='FangSong')plt.subplot(212)plt.plot(freqs, xfp)plt.xlabel(u"頻率(Hz)", fontproperties='FangSong')#字體FangSongplt.ylabel(u'幅值', fontproperties='FangSong')plt.subplots_adjust(hspace=0.4)'''subplots_adjust(left=None, bottom=None, right=None, top=None, wspace=None, hspace=None)有六個(gè)可選參數(shù)來控制子圖布局。值均為0~1之間。其中l(wèi)eft、bottom、right、top圍成的區(qū)域就是子圖的區(qū)域。wspace、hspace分別表示子圖之間左右、上下的間距。實(shí)際的默認(rèn)值由matplotlibrc文件控制的。'''plt.show()test_fft()

2.結(jié)果顯示

總結(jié)

以上是生活随笔為你收集整理的频谱分析：基于python画出时域频域波形的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。