相关文章推荐
发财的楼房  ·  CTFT、DTFT、DFT、FFT - ·  14 小时前    · 
内向的鞭炮  ·  “汉语拼音之父”周有光在112岁生日后第二天 ...·  2 年前    · 
自信的楼梯  ·  吕鹏教授受聘北京通用人工智能研究院研究员-公 ...·  2 年前    · 
考研的海豚  ·  传说还是史实:有关“五谷”的考古发现_植物·  2 年前    · 
忐忑的眼镜  ·  路亚竿硬度、调性、感度是什么? - 知乎·  2 年前    · 
刚失恋的香蕉  ·  积极防范化解重大经济金融风险---党建网·  2 年前    · 
小百科  ›  CTFT、DTFT、DFT、FFT -
DFT 关系运算 信号频率 dft fft
发财的楼房
14 小时前
打开 App

CTFT、DTFT、DFT、FFT

  • 学习 ·
  • 2022-07-14 02:42 ·
  • 672浏览
毕导RAM
关注
  • 欧拉公式

    sinx%3D%5Cfrac%7B(e%5E%7Bjx%7D-e%5E%7B-%20jx%7D)%7D%7B2j%7D
    cosx%3D%5Cfrac%7B(e%5E%7Bjx%7D%2Be%5E%7B-jx%7D)%7D%7B2%7D

    e%5E%7B%C2%B1ix%7D%3Dcosx%C2%B1jsinx

【 1. CTFT 】

  • 连续时间傅里叶变换 (Continuous Time Fourier transform)

  • 处理对象: 连续时间信号

    CTFT%EF%BC%9AX_a(j%CE%A9)%3D%5Cint_%7B-%E2%88%9E%7D%5E%7B%2B%E2%88%9E%7Dx_a(t)e%5E%7B-j%CE%A9t%7Ddt

    ICTFT%EF%BC%9Ax_a(t)%3D%5Cfrac%7B1%7D%7B2%CF%80%7D%5Cint_%7B-%E2%88%9E%7D%5E%7B%2B%E2%88%9E%7DX_a(j%CE%A9)e%5E%7Bj%CE%A9t%7Dd%CE%A9

【 2. DTFT 】

  • 离散时间傅里叶变换 (Discrete Time Fourier transform)

  • 处理对象: 离散时间信号

    DTFT%EF%BC%9AX(e%5E%7Bjw%7D)%EF%BC%9D%5Csum_%7Bn%EF%BC%9D-%E2%88%9E%7D%5E%7B%2B%E2%88%9E%7Dx%5Bn%5De%5E%7B-jwn%7D

    IDTFT%EF%BC%9Ax%5Bn%5D%EF%BC%9D%5Cfrac%7B1%7D%7B2%CF%80%7D%5Cint_%7B-%CF%80%7D%5E%7B%CF%80%7DX(e%5E%7Bjw%7D)e%5E%7Bjwn%7Ddw

【 3. DFT 】

  • 离散傅里叶变换 (Discrete Fourier transform)

  • 处理对象: 离散时间信号

  • DFT 和 DTFT 的关系
    在 DTFT 基础上,w轴上的 [0,2π] 范围内,对 X(e%5E%7Bjw%7D) 以 w_k%3D%5Cfrac%7B2%5Cpi%20k%7D%7BN%7D 的频率均匀抽样N次,即序列X[k]就是 X(e%5E%7Bjw%7D) 在N个等间隔频率 w_k%3D%5Cfrac%7B2%5Cpi%20k%7D%7BN%7D 上的一组频率样本。

    DFT%EF%BC%9A%CE%A7%5Bk%5D%3D%5Csum_%7Bn%3D0%7D%5E%7BN-1%7Dx%5Bn%5De%5E%7B-%5Cfrac%7Bj2%5Cpi%20kn%7D%7BN%7D%7D%3D%5Csum_%7Bn%3D0%7D%5E%7BN-1%7Dx%5Bn%5DW_%7BN%7D%5E%7Bkn%7D%EF%BC%8C0%E2%89%A4k%E2%89%A4N-1

    IDFT%EF%BC%9Ax%5Bk%5D%3D%5Cfrac%7B1%7D%7BN%7D%5Csum_%7Bk%3D0%7D%5E%7BN-1%7DX%5Bk%5DW_%7BN%7D%5E%7B-kn%7D%EF%BC%8C0%E2%89%A4n%E2%89%A4N-1

【 4. CTFT、DTFT、DFT 之间的关系 】

【 5. FFT 】

  • 快速傅里叶变换 (Fast Fourier transform)

  • FFT 和 DFT 之间的关系 :FFT是一种 DFT 的高效算法

5.1 Matlab 对 fft 函数的阐述

  • 采样点数N(xn的长度) > FFT运算点数M

    fft() 函数对信号 xn 进行截断操作即只取N个采样点中的前M个点进行分析。

    幅度轴:

    X%3D%5Cfrac%7BX*2%7D%7BM%7D%EF%BC%8C%E7%9B%B4%E6%B5%81%E5%88%86%E9%87%8F%E5%8D%B3X(0)%20%E7%9A%84%E5%B9%85%E5%BA%A6%E5%85%B3%E7%B3%BB%E5%BC%8F%E4%B8%BA%20X(0)%3D%5Cfrac%7BX(0)%7D%7BM%7D%20

    频率轴:

    %5Cfrac%7Bfs*(0%3A1%3AM-1)%7D%7BM-1%7D

  • 采样点数N(xn的长度) < FFT运算点数M

    fft() 函数对信号 xn 进行尾补零操作即在该信号尾部添加多个值为0的数据点以使信号总点数N增至FFT运算所需点数M。

    幅度轴:

    X%3D%5Cfrac%7BX*2%7D%7BN%7D%EF%BC%8C%E7%9B%B4%E6%B5%81%E5%88%86%E9%87%8F%E5%8D%B3X(0)%20%E7%9A%84%E5%B9%85%E5%BA%A6%E5%85%B3%E7%B3%BB%E5%BC%8F%E4%B8%BA%20X(0)%3D%5Cfrac%7BX(0)%7D%7BN%7D%20

    频率轴:

    %5Cfrac%7Bfs*(0%3A1%3AM-1)%7D%7BM-1%7D

5.2 FFT运算点数对结果的影响

为了更好的分析单一变量对结果的影响,本例不考虑频谱泄漏。

采样频率、采样点数固定时,随 FFT运算点数增大,精度不断提高 。

  • 图1、图2和图3 :FFT运算点数M < 采样点数N,此时由由 fft() 函数对信号 xn 进行截断操作即只取N个采样点中的前M个点进行分析。随着FFT运算点数的增大,精度逐渐提高。幅度的关系式为 X=X*2/M。

  • 图4、图5和图6 :FFT运算点数M > 采样点数N,此时由 fft() 函数对信号 xn 进行尾补零操作即在该信号尾部添加多个值为0的数据点以使信号总点数N增至FFT运算所需点数M。随着FFT运算点数的增大,精度逐渐提高。幅度的关系式为 X=X*2/N。

5.3 采样点数对结果的影响

采样频率、FFT运算点数 固定时,随着 采样点数 的提高,精度逐渐提高直至不变。

  • 图1、图2和图3 :采样点数N < FFT运算点数M,此时由 fft() 函数对信号 xn 进行尾补零操作即在该信号尾部添加多个值为0的数据点以使信号总点数N增至FFT运算所需点数M。而随着采样点数的提高,原信号的点数增多,分析时所包含的信息也相应增多,精度逐渐提高。幅度的关系式为 X=X*2/N。

  • 图4、图5和图6 :采样点数N > FFT运算点数M,此时由 fft() 函数对信号 xn 进行截断操作即只取N个采样点中的前M个点进行分析。因此,随着采样点数的提高,精度不变。幅度的关系式为 X=X*2/M。

【 6. 分辨率 】

6.1 时间分辨率

  • 时间分辨率是指在同一区域进行的相邻两次观测的最小时间间隔。时间间隔越小,时间分辨率越高。

  • 反映在FFT中,就相当于对一段序列的进行采样的采样周期。

6.2 频率分辨率

  • 频率分辨率是指将两个相邻谱峰分开的能力,在实际应用中是指分辨两个不同频率信号的最小间隔。

  • 反映在FFT中,频率分辨率可以理解为在频率轴上的所能得到的最小频率间隔。

    %E5%86%B3%E5%AE%9A%E5%BC%8F%EF%BC%9Af_0%3D%5Cfrac%7B1%7D%7BL%7D%2C%20L%E4%B8%BA%E8%A2%AB%E9%87%87%E6%A0%B7%E4%BF%A1%E5%8F%B7%E7%9A%84%E6%97%B6%E9%97%B4%E9%95%BF%E5%BA%A6

    %E5%AE%9A%E4%B9%89%E5%BC%8F%EF%BC%9Af_0%3D%5Cfrac%7Bf_s%7D%7BN%7D%2C%E5%85%B6%E4%B8%ADN%E4%B8%BAFFT%E8%BF%90%E7%AE%97%E7%82%B9%E6%95%B0%2C%20f_s%E4%B8%BA%E9%87%87%E6%A0%B7%E9%A2%91%E7%8E%87

  • 所以,当被采样信号长度一确定,信号的最大周期就是该长度,信号的最小频率(即频率分辨率)就是原始信号时间长度的倒数。

6.3 频率分辨率和准确度的关系

  • 只有当原信号频率 f0 正好等于频率分辨率的整数倍时才能获得该信号频率及其幅值的准确值。

【 7. 信号补零 】

  • FFT运算点数M > 采样点数N 时,fft(xn,M) 函数对信号 xn 进行尾补零操作即在该信号尾部添加多个值为0的数据点以使信号总点数N增至FFT运算所需点数M。

  • 补1次零相当于在原始频谱图中每两个频率之间插入1个频率值,补2次零相当于在原始频谱图中每两个频率之间插入2个频率值,并且原始频率值的位置及其幅值保持不变。因此, 补零会使频谱图中的频率点的数量增加,从而使得频谱图更加的光滑连续,但是补零不能对频谱图中的 频率分辨率、频率值以及幅值有所改善 。

  • 在补零操作中,所补的零值点的数量所占分析信号比例越大,分析结果中的幅度越小。

  • 补零的好处:
    ① 使数据N为2的整次幂,便于使用 FFT。
    ② 补零后,其实是对DFT结果做了插值,克服“栅栏"效应,使谱外观平滑化。我把“栅栏"效应形象理解为,就像站在栅栏旁边透过栅栏看外面风景,栅栏会挡住比较多风景,此时就可能漏掉较大频域分量,但是补零以后,相当于你站远了,风景就看的越来越清楚了。
    ③ 由于对时域数据的截短必然造成频谱泄露,因此在频谱中可能出现难以辨认的谱峰,补零在一定程度上能消除这种现象。

【 8. 频谱泄漏 】

  • 对频率为 f0 的原信号 xt 以 Fs 的采样频率进行采样,采样点数为 N ,得到离散序列值 xn ,FFT 运算点数为M,由于 M≠N 时,fft() 函数会进行尾补零或截断操作,因此 令 fft() 函数经过上述操作后所得分析信息中原信号 f0 的周期数为 k,若 k 为非整数,则相当于对信号进行非整周期截断,此时会造成严重的频谱泄露现象,造成频率值和幅度值的偏差。

  • 通过增加 FFT 所分析信息中原信号 f0 的周期数 k 可以改善频谱泄露的现象,减小偏差 。当周期数增多时,非整周期部分占整体信号的比例逐渐减少,对整个信号的影响也逐渐变小,频谱泄露逐渐减少,频率值逐渐靠近真实频率值;当周期数足够多时能够得到较准确的频率值,幅值会收敛到一个常数。

【 9. 信号的周期延拓 】

  • 有时我们采集到的信号比较短,这样得到的频率分辨率比较低,这时可能会想如果将信号进行周期延拓可以增加信号的长度。

  • 整周期延拓时两段信号之间是连续的。 非整周期截断会导致频谱中的频率及其幅值不准确,对信号进行延拓后相邻两段信号之间会出现间断点,使得信号之间的连接不连续。

  • 综合整周期截断和非整周期截断后的延拓次数对频谱的影响情况来看, 对信号进行延拓相当于增加了采样时间,能够提高频率分辨率,但是不会改善频谱中频率及其幅值的准确性。


本文禁止转载或摘编

本文为我原创

  • FFT
  • DFT
  • DTFT
  • CTFT
  • 信号补零
  • FFT运算点数
  • 抽样点数
  • 频谱泄漏
  • --
  • --
  • --
展开阅读全文
打开App,看更多精彩内容
 
推荐文章
发财的楼房  ·  CTFT、DTFT、DFT、FFT -
14 小时前
内向的鞭炮  ·  “汉语拼音之父”周有光在112岁生日后第二天去世,一个名副其实的 ...
2 年前
自信的楼梯  ·  吕鹏教授受聘北京通用人工智能研究院研究员-公共管理学院门户网站
2 年前
考研的海豚  ·  传说还是史实:有关“五谷”的考古发现_植物
2 年前
忐忑的眼镜  ·  路亚竿硬度、调性、感度是什么? - 知乎
2 年前
刚失恋的香蕉  ·  积极防范化解重大经济金融风险---党建网
2 年前
今天看啥   ·   Py中国   ·   codingpro   ·   小百科   ·   link之家   ·   卧龙AI搜索
删除内容请联系邮箱 2879853325@qq.com
小百科 - 百科知识指南
© 2024 ~ 沪ICP备11025650号