II. Short-time Fourier Transform

Slides:



Advertisements
Similar presentations
1 時 頻 分 析 近 年 來 的 發 展時 頻 分 析 近 年 來 的 發 展 丁 建 均 國立台灣大學電信工程學研究所 Recent Development of Time-Frequency Analysis.
Advertisements

663 Chapter 14 Integral Transform Method Integral transform 可以表示成如下的積分式的 transform  kernel Laplace transform is one of the integral transform 本章討論的 integral.
第 1 章 信號與系統簡介 by 胡興民老師 連續時間信號與離散時間信號 連續時間信號 (continuous-time signal) :連續時間 信號以函數 x(t) 表示之,其中 t 是連續時間變數 。 離散時間信號 (discrete-time signal) :離散時間信 號只定義在離散的時間點上,一般以離散時間變數.
仪 容. 一、化妆的技巧 眼部的化妆 唇部化妆 眉部化妆 鼻部化妆 根据脸型化妆 根据脸型选发型.
Final Review Chapter 1 Discrete-time signal and system 1. 模拟信号数字化过程的原理框图 使用 ADC 变换器对连续信号进行采样的过程 使用 ADC 变换器对连续信号进行采样的过程 x(t) Analog.
精神疾病与社区处理.
第四章 平稳过程.
司法体制改革与律师执业前景瞻望 黄太云
加油添醋話擴寫 日新國小 鄒彩完.
创新大赛经验浅谈 高二(18)班 黄佳淇.
數位訊號處理 第4章 離散時間訊號與LTI系統之傅利葉分析
乳腺中心实验室 2012级硕士 李满秀.
師資培育中心外埠教育參觀.
劳动统计专业年报培训 社会科 洪惠娟 2009年11月.
数字图像处理 Digital Image Processing.

他們,與眾不同…….
XI. Hilbert Huang Transform (HHT)
Time and frequency domain
III. Gabor Transform III-A Definition Standard Definition:
3-3 Modeling with Systems of DEs
AN INTRODUCTION TO OFDM
IV. Implementation IV-A Method 1: Direct Implementation 以 STFT 為例
Applications of Digital Signal Processing
期末考的範圍遠遠多於期中考,要了解的定理和觀念也非常多
V. Homomorphic Signal Processing
Differential Equations (DE)
第五章 数字滤波器设计 Filtering Beijing Institute of Technology 数字信号处理.
期末考的範圍遠遠多於期中考,要了解的定理和觀念也非常多
GoldWave 音樂編輯軟體 -演算法與技巧- 學生: 冼達 指導教授: 丁建均 DISP, NTU 2018/11/20.
§5.6 Hole-Burning and The Lamb Dip in Doppler- Broadened Gas Laser
X. Other Applications of Time-Frequency Analysis
第二單元 面積與黎曼和.
第二十九單元 方向導數與梯度.
信号与图像处理基础 An Introduction to Signal and Image Processing 中国科学技术大学 自动化系
無線通訊系統概論 行動通訊與網路 Chapter 7 多重分工技術.
加油添醋話擴寫 鄒彩完.
第2章 短时傅立叶变换 2.1 连续信号的短时傅立叶变换 2.2 短时傅立叶反变换 2.3 离散信号的短时傅立叶变换
語音處理簡介.
本章大綱 2.1 The Limit of a Function函數的極限 2.2 Limit Laws極限的性質
VI. Brief Introduction for Acoustics
Principle and Application of Digital Television
消費者偏好與效用概念.
一般論文的格式 註:這裡指的是一般 journal papers 和 conference papers 的格式。
Semantic Navigation Liang Zheng.
Advanced Digital Signal Processing 高等數位訊號處理
Chp.4 The Discount Factor
第三章 付里叶分析 离散付氏级数的数学解释(The Mathematical Explanation of DFS)
1 离散信号 2019/4/10.
XIV. Orthogonal Transform and Multiplexing
Chp.4 The Discount Factor
虚 拟 仪 器 virtual instrument
通信工程专业英语 Lesson 13 Phase-Locked Loops 第13课 锁相环
Chp.4 The Discount Factor
第4章 连续时间傅立叶变换 The Continuous-Time Fourier Transform
XI. Hilbert Huang Transform (HHT)
(二)盲信号分离.
IV. Implementation IV-A Method 1: Direct Implementation 以 STFT 為例
第 8 章 計量與質性預測變數之迴歸模型.
96學年度第二學期電機系教學助理課後輔導進度表(一)(查堂重點)
5. Combinational Logic Analysis
本講義為使用「訊號與系統,王小川編寫,全華圖書公司出版」之輔助教材
II. Short-time Fourier Transform
Surface wave dispersion measurements using Hilbert-Huang Transform
定语从句(4).
Principle and application of optical information technology
4.理財規劃者適格性分析與實作 理財規劃重點 生涯階段 「就業前準備階段」(學習階段) 「初入社會階段」 「確定職涯階段」 「維持職涯階段」
Voronoi Diagram and Delaunay Triangulation
Gaussian Process Ruohua Shi Meeting
Hybrid fractal zerotree wavelet image coding
Presentation transcript:

II. Short-time Fourier Transform II-A Definition Short-time Fourier transform (STFT) Alternative definition 參考資料 [1] S. Qian and D. Chen, Section 3-1 in Joint Time-Frequency Analysis: Methods and Applications, Prentice-Hall, 1996. [2] S. H. Nawab and T. F. Quatieri, “Short time Fourier transform,” in Advanced Topics in Signal Processing, pp. 289-337, Prentice Hall, 1987.

STFT Inverse of the STFT: To recover x(t), where w(t1 – t)  0. For the alternative definition,

The mask function w(t) always has the property of (a) even: w(t) = w(t), (通常要求這個條件要滿足) (b) max(w(t)) = w(0), w(t1)  w(t2) if |t2| > |t1| (c) w(t)  0 when |t| is large w(t) = (t) (triangular function) w(t) = exp(a|t|b) (hyper-Laplacian function) t =  1 t = 1 Max[(t)] = 1

II-B Rec-STFT Rectangular mask STFT (rec-STFT) Inverse of the rec-STFT where t – B < t1 < t + B The simplest form of the STFT Other types of the STFT may require more computation time than the rec-STFT.

II-C Properties of the Rec-STFT (1) Integration (recovery): (a) (b) when v B < t < v + B, = 0 otherwise

(2) Shifting property (橫的方向移動) (3) Modulation property (縱的方向移動)

(4) Special inputs: (1) When x(t) = (t), when –B < t < B, otherwise (2) When x(t) = 1 思考: B 值的大小,對解析度的影響是什麼?

(5) Linearity property If h(t) =  x(t) +  y(t) and H(t, f ), X(t, f ) and Y(t, f ) are their rec-STFTs, then H(t, f ) =  X(t, f ) +  Y(t, f ). (6) Power integration property (7) Energy sum property (Parseval’s theorem)

思考: (1) 哪些性質 Fourier transform 也有? (2) 其他型態的 STFT 是否有類似的性質? Shifting Modulation

Example: x(t) = cos(2 t) when t < 10,

II-D Advantage and Disadvantage  Compared with the Fourier transform: All the time-frequency analysis methods has the advantage of: The instantaneous frequency can be observed. All the time-frequency analysis methods has the disadvantage of: Higher complexity for computation

 Compared with other types of time-frequency analysis: The rec-STFT has an advantage of the least computation time for digital implementation but its performance is worse than other types of time-frequency analysis.

II-E STFT with Other Windows (1) Rectangle (2) Triangle (3) Hanning (4) Hamming (5) Gaussian -B B -B B

(6) Asymmetric window -B1 t = 0 B2 t-axis B1  B2 應用: seismic wave analysis, collision detection (The applications that require real-time processing) onset detection

動腦思考: Are there other ways to choose the mask of the STFT? Which mask is better? 沒有一定的答案

II-F Spectrogram STFT 的絕對值平方,被稱作 Spectrogram

附錄二:使用 Matlab 將時頻分析結果 Show 出來 可採行兩種方式: (1) 使用 mesh 指令畫出立體圖 (但結果不一定清楚,且執行時間較久) (2) 將 amplitude 變為 gray-level,用顯示灰階圖的方法將結果表現出來 假設 y 是時頻分析計算的結果 image(abs(y)/max(max(abs(y)))*C) % C 是一個常數,我習慣選 C=400 colormap(gray(256)) % 變成 gray-level 的圖 set(gca,‘Ydir’,‘normal’) % 若沒這一行, y-axis 的方向是倒過來的

set(gca,‘Fontsize’,12) % 改變橫縱軸數值的 font sizes xlabel('Time (Sec)','Fontsize',12) % x-axis ylabel('Frequency (Hz)','Fontsize',12) % y-axis title(‘STFT of x(t)','Fontsize',12) % title 計算程式執行時間的指令: tic (這指令如同按下碼錶) toc (show 出碼錶按下後已經執行了多少時間) 註:通常程式執行第一次時,由於要做程式的編譯,所得出的執行時間會比較長 程式執行第二次以後所得出的執行時間,是較為正確的結果