短時傅利葉變換

2021-10-22 09:55:27 字數 2415 閱讀 6993

首先用一張直觀的圖來看一下什麼是傅利葉變換

圖1:平穩訊號的傅利葉變換

做完fft(快速傅利葉變換)後,可以在頻譜上看到清晰的四條線,訊號包含四個頻率成分。

傅利葉變換是一種分析訊號的方法,它可分析訊號的成分,也可用這些成分合成訊號。在分析訊號時,主要應用於處理平穩訊號,通過傅利葉變換可以獲取一段訊號總體上包含哪些頻率的成分,但是對各成分出現的時刻無法得知。

圖2:非平穩訊號的傅利葉變換

最上方的圖是頻率始終不變的訊號,下方兩張**率隨時間變化,但是三張圖都包含同樣頻率的四個成分。這三個在時域上差異巨大的訊號,在頻譜上卻很相似。由此可見,傅利葉變換會忽略訊號的時間資訊,對非平穩訊號的處理有天生的缺陷。對於非平穩訊號,只知道包含哪些頻率成分是不夠的,還需要知道每個頻率出現的時間,即做時頻分析

由此出現了短時傅利葉變換。短時傅利葉變換的精髓就是加窗。「把整個時域過程分解成無數個等長的小過程,每個小過程近似平穩,再傅利葉變換,就知道在哪個時間點上出現了什麼頻率了。」

圖3:對語音頻號加窗

如圖,將訊號在時頻上分成一段一段,然後做fft,就可以知道訊號頻率隨時間變化的關係了。

但是這種方法也有缺陷,就是窗的寬窄難以確定。

圖4:窗過窄

窗太窄,每個窗內訊號太短,會導致頻率分辨不精準,頻率解析度差

圖5:窗過寬

窗太寬,時域上不夠精細,時間解析度低

【這個道理可以用海森堡不確定性原理來解釋。類似於我們不能同時獲取乙個粒子的動量和位置,我們也不能同時獲取訊號絕對精準的時刻和頻率。這也是 一對不可兼得的矛盾體。我們不知道在某個瞬間哪個頻率分量存在,我們知道的只能是在乙個時間段內某個頻帶的分量存在。 所以絕對意義的瞬時頻率是不存在的。)】

圖6:不同大小窗

同乙個訊號(包含四種頻率)採用不同寬度的窗做stft。用窄窗,食品土在時間軸上解析度很高,幾個峰基本成矩形。用寬窗,變成了綿延的矮山,但是頻率軸上,寬窗解析度更高。所以窄視窗時間解析度高,頻率解析度低,高頻適合小視窗,低頻適合大視窗。但短時傅利葉變換視窗大小確定,無法滿足非穩態訊號變化的頻率的需求

stft給訊號加窗,分段做fft。而小波直接把傅利葉變換的基給換了——將無限長的三角函式基換成了有限長的會衰減的小波基。這樣不僅能夠獲取頻率,還可以定位到時間了【暫未用到,以後用的的話補充】

對乙個時域的訊號做傅利葉變換,會得到

1. 一系列的頻率

2. 每個頻率下的振幅

3.每個頻率下的相位

再將這三個訊號組合起來反變換,可變成時域訊號

傅利葉級數的表示公式

任何乙個週期性函式都可以看成正余弦函式的疊加。

可得到每乙個訊號的頻率、振幅以及它的相位

標準正交基為1,sinnwt,cosnwt

橫軸單位為1,縱軸單位為i,對於任一點,橫座標是cosθ,縱座標是isinθ。點座標即cosθ+isinθ

傅利葉變換:

傅利葉逆變換:

1.聲音

大腦可以對聲音自動進行傅利葉變換,從而根據分出來的頻率大小確定說話人的性別。且人腦可以自動進行濾波,去噪

2.影象

橫座標是空間位置。將人的影象做傅利葉變換,得到低頻成分:表示人的輪廓。並接著得到高頻成分,表示人的細節。如果只要低頻,只要輪廓,可以直接把高頻成分濾掉。 

傅利葉-->短時傅利葉變換-->小波變換的方法

對音訊訊號作短時傅利葉變換(stft)/小波變換處理(python + matlab)

短時傅利葉變換

短時傅利葉變換 stft,short time fourier transform,或 short term fourier transform 是和傅利葉變換相關的一種數學變換,用以確定時變訊號其區域性區域正弦波的頻率與相位。它的思想是 選擇乙個時頻區域性化的窗函式,假定分析窗函式g t 在乙個短...

短時傅利葉變換

時間解析度和頻率解析度 時間解析度 訊號頻率隨時間變化,要將這種頻率變化分辨出來。自然,窗越短越好,以使得在窗內訊號頻率近似不變。例如,100hz的訊號和100.1hz的訊號疊加,一兩個週期恐怕看不出來,必須要足夠多的週期才能區別開。短時傅利葉變換可以看做移位訊號x n m 通過窗w m 的傅利葉變...

短時傅利葉變換原理解

前一段時間專案需要學習了短時傅利葉變換,今天我來總結一下現階段對短時傅利葉變換的理解。短時傅利葉變換是最常用的一種時頻分析方法,它通過時間窗內的一段訊號來表示某一時刻的訊號特徵。在短時傅利葉變換過程中,窗的長度決定頻譜圖的時間解析度和頻率解析度,窗長越長,擷取的訊號越長,訊號越長,傅利葉變換後頻率解...