亚洲免费在线-亚洲免费在线播放-亚洲免费在线观看-亚洲免费在线观看视频-亚洲免费在线看-亚洲免费在线视频

使用python實現離散時間傅里葉變換的方法

系統 2121 0

我們經常使用傅里葉變換來計算數字信號的頻譜,進而分析數字信號,離散時間傅里葉變換的公式為:

可是自己動手實現一遍才是最好的學習。

在數字分析里面,傅里葉變換默認等時間間隔采樣,不需要時間序列,只需要信號數組即可分析。

分析過程如下:

  • 對于含有 n 個樣本值的數字信號序列,根據奈奎斯特采樣定律,包含的周期數最大為 n/2,周期數為 0 代表直流分量。所以,當周期數表示為離散的 0,1,2,3…n/2 ,總的數目為 n/2+1
  • 傅里葉變換之后的結果為復數, 下標為 k 的復數 a+b*j 表示時域信號中周期為 N/k 個取樣值的正弦波和余弦波的成分的多少, 其中 a 表示 cos 波形的成分, b 表示 sin 波形的成分
  • 首先產生一個長度為 n,一倍周期的 $e^{-jwn} $ (即為 $cos(wn)-jsin(wn) $ )波樣本序列.
  • 將數字信號序列中的每一個樣本與 1 倍周期的樣本波形序列相乘,得到 n 個乘積,將 n 個乘積相加,放入 f[1] 中。
  • 再產生一個長度為 n,兩倍周期的 $e^{-jwn} $ (即為 $cos(wn)-jsin(wn) $ )波樣本序列,再將數字信號序列中的每一個樣本與 2 倍周期的樣本波形序列相乘,得到 n 個乘積,將 n 個乘積相加,放入 f[2] 中。依次重復。
  • 對于 0 倍周期,即直流分量來說,可以認為產生的是 0 倍周期的樣本波形,重復操作,放入 f[0] 即可。
  • 這樣就得到了數字信號序列的傅里葉變換

使用方法:

從以上過程得到數字序列的傅里葉變換之后,如果想要得到真正頻譜,還需要做處理:

  • 計算出的每一個頻率下的幅值需要除以時間序列的長度,類似求平均的過程
  • 每一個頻率下的幅值是一個復數,需要對它求模,而且因為在負頻率處也有值,所以需要對于實信號需要乘 2
  • 頻率的序列為 0 到采樣率的一半,長度為 n/2+1

完整程序:

            
# 離散時間傅里葉變換的 python 實現
import numpy as np
import math
import pylab as pl
import scipy.signal as signal
import matplotlib.pyplot as plt

sampling_rate=1000
t1=np.arange(0, 10.0, 1.0/sampling_rate)
x1 =np.sin(15*np.pi*t1)

# 傅里葉變換
def fft1(xx):
#   t=np.arange(0, s)
  t=np.linspace(0, 1.0, len(xx))
  f = np.arange(len(xx)/2+1, dtype=complex)
  for index in range(len(f)):
    f[index]=complex(np.sum(np.cos(2*np.pi*index*t)*xx), -np.sum(np.sin(2*np.pi*index*t)*xx))
  return f

# len(x1)
          
            
xf=fft1(x1)/len(x1)
freqs = np.linspace(0, sampling_rate/2, len(x1)/2+1)
plt.figure(figsize=(16,4))
plt.plot(freqs,2*np.abs(xf),'r--')

plt.xlabel("Frequency(Hz)")
plt.ylabel("Amplitude($m$)")
plt.title("Amplitude-Frequency curve")

plt.show()
          

使用python實現離散時間傅里葉變換的方法_第1張圖片

            
plt.figure(figsize=(16,4))
plt.plot(freqs,2*np.abs(xf),'r--')

plt.xlabel("Frequency(Hz)")
plt.ylabel("Amplitude($m$)")
plt.title("Amplitude-Frequency curve")
plt.xlim(0,20)
plt.show()
          

使用python實現離散時間傅里葉變換的方法_第2張圖片

此處實現的是傳統的傅里葉變換,這種方法實際已經不用了,現在使用快速傅里葉變換,其實兩種是等價的,但是快速傅里葉變換時間復雜度要小很多。

以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。


更多文章、技術交流、商務合作、聯系博主

微信掃碼或搜索:z360901061

微信掃一掃加我為好友

QQ號聯系: 360901061

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描下面二維碼支持博主2元、5元、10元、20元等您想捐的金額吧,狠狠點擊下面給點支持吧,站長非常感激您!手機微信長按不能支付解決辦法:請將微信支付二維碼保存到相冊,切換到微信,然后點擊微信右上角掃一掃功能,選擇支付二維碼完成支付。

【本文對您有幫助就好】

您的支持是博主寫作最大的動力,如果您喜歡我的文章,感覺我的文章對您有幫助,請用微信掃描上面二維碼支持博主2元、5元、10元、自定義金額等您想捐的金額吧,站長會非常 感謝您的哦!!!

發表我的評論
最新評論 總共0條評論
主站蜘蛛池模板: 国产欧美日韩高清专区手机版 | 国产一区二区在线不卡 | 国产免费自拍视频 | 国产精品成人一区二区不卡 | 夜夜爽天天狠狠九月婷婷 | 黄色a级毛片 | 久久成人亚洲香蕉草草 | 精品久久天干天天天按摩 | free性欧美极度另类超级大 | 国产成人精品三级91在线影院 | a级毛片在线视频免费观看 a级免费网站 | 国产色产综合色产在线观看视频 | 久久国产精品只做精品 | 日韩啊啊啊 | 精品久久久久久久久久久久久久久 | 91最新免费地址入口 | 国产午夜精品福利视频 | 天天插天天干天天操 | 国产婷婷色综合成人精品 | 亚洲欧美另类在线观看 | 久久国产精品一国产精品金尊 | 日本一级毛片视频网站 | 在线观看国产一区二三区 | 亚洲另类第一页 | 亚洲高清国产一线久久 | 特级毛片免费观看视频 | 热久久影院 | 深夜激情网 | 欧美洲大黑香蕉在线视频 | 一本大道加勒比久久综合 | 久久精品这里精品 | 欧美一级毛片日本 | 男人天堂.com | 综合久久久久综合 | 亚洲久久草| 国产精品v欧美精品∨日韩 国产精品v一区二区三区 | a级毛片高清免费视频 | 四虎免费影院在线播放 | 欧美一级毛片aaa片 欧美一级毛片不卡免费观看 | 日韩最新中文字幕 | 精品国产一区二区二三区在线观看 |