码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 数学笔记;离散傅里叶变化 DFT


    1 原理介绍 

    • 离散傅里叶变化是连续傅里叶变化在如下信号下的等价形式:
      • N个样本
      • 每个样本的采样间隔是T
    • 令f(t)是数据源的持续信号,将N个样本表示成f[0],f[1],f[2],.....f[k],....f[N-1]
      • ——>原始数据源信号f(t)的傅里叶变换,应该是

       

    • 因为信号脉冲仅在样本点采集,所以可以如下近似:

     

    • 这里,DFT假设波形是周期性的(周期正好是T个interval的跨度)
      • 也就是这里我们是从0到N-1采样的f(0)到f(N-1),和从N到2N-1采样得到的f(N)到f(2N-1)是一样的
      • 比如我们采样了10个点(0~10),那么DFT会隐式地人为周期就是10

       

    • 所以上式中的ω被表示为:

    • 所以离散傅里叶级数可以表示为:

       

    • 上式可以被重写成: 

     1.1 欧拉公式 (复习)

     

     

     2 举例

    令连续信号为:f[k]=5+2cos(\frac{\pi}{2}k-90^\cdot)+3cos(4\pi t)

    采样k=0~k=3 四个点,得到f[0]=8,f[1]=4,f[2]=8,f[3]=0

    于是 

    可以写成:

    F[n]=\sum_{k=0}^3 f[k] e^{-j \frac{2\pi}{4}kn}=\sum_{k=0}^3 f[k] (-j)^{kn}

    用矩阵的形式,可以写成

    \begin{bmatrix} F[0]]\\ F[1]\\ F[2]\\ F[3] \end{bmatrix} =\begin{bmatrix} 1 & 1 & 1 & 1\\ (-j)^0 &(-j)^1 &(-j)^2 &(-j)^3 \\ [(-j)^2]^0 & [(-j)^2]^1 & [(-j)^2]^2 & [(-j)^2]^3 \\ [(-j)^3]^0 & [(-j)^3]^1 & [(-j)^3]^2 & [(-j)^3]^3 \end{bmatrix} \begin{bmatrix} f[0]\\ f[1]\\ f[2]\\ f[3] \end{bmatrix}

                  =\begin{bmatrix} 1 & 1 & 1 & 1\\ 1 & -j & -1 & j \\ 1 & -1 & 1 & -1 \\ 1 & j & -1 & -j \end{bmatrix} \begin{bmatrix} 8\\ 4\\ 8\\ 0\end{bmatrix} =\begin{bmatrix} 20\\ -4j\\ 12\\ 4j \end{bmatrix}

                  

    4 逆傅里叶变化 

    的逆傅里叶变化是:

     

     如果这里令W=exp(j2\pi/N)的话,那么逆傅里叶变化的矩阵形式为

     4.1 举例

    还是用前面一个例子

    已知F[0]=20,F[1]=-4j,F[2]=12,F[3]=4j

    希望得到:f[0]=8,f[1]=4,f[2]=8,f[3]=0

    \begin{bmatrix} f[0]]\\ f[1]\\ f[2]\\ f[3] \end{bmatrix} =\frac{1}{4}\begin{bmatrix} 1 & 1 & 1 & 1\\ (j)^0 &(j)^1 &(j)^2 &(j)^3 \\ [(j)^2]^0 & [(j)^2]^1 & [(j)^2]^2 & [(j)^2]^3 \\ [(j)^3]^0 & [(j)^3]^1 & [(j)^3]^2 & [(j)^3]^3 \end{bmatrix} \begin{bmatrix} F[0]\\ F[1]\\ F[2]\\ F[3] \end{bmatrix}

    =\frac{1}{4}\begin{bmatrix} 1 & 1 & 1 & 1\\ 1 & j & -1 & -j \\ 1 & -1 & 1 & -1 \\ 1 & -j & -1 & j \end{bmatrix} \begin{bmatrix} 20\\ -4j\\ 12\\ 4j\end{bmatrix} =\begin{bmatrix} 8\\ 4\\ 8\\ 0 \end{bmatrix}

  • 相关阅读:
    制作OpenCV相机标定板棋盘格图像
    异常解决!前后端交互时,CORS跨域问题的解决过程
    Python_数据容器_集合set
    Python学习记录 异常处理
    原生js做打地鼠游戏
    Ngrok实现内网穿透
    商品分类,汉码批发进销存管理软件
    【科技素养】蓝桥杯STEMA 科技素养组模拟练习试卷F
    科技云报道:汽车云,云计算换挡后的下一个“赛点”?
    机器学习评价指标(分类、目标检测)
  • 原文地址:https://blog.csdn.net/qq_40206371/article/details/126772853
  • 最新文章
  • 沪漂五周年了:我越来越迷茫了
    Agentic Skill Routing 实战:别再把所有 Skill 塞进 AI Agent 上下文
    MySQL-Seconds_behind_master的精度误差
    [MAF预定义ChatClient中间件-03]CachingChatClient——利用缓存省钱省时间
    AI的至暗历史:从万众期待到被政府撤资,AI的两次死亡徘徊
    Agent OS :五种驯服不确定性的范式
    PortSwigger SQL注入LAB11
    数据库即时编译JIT
    [Begin]AI Learn Data Day 0
    深度学习进阶(二十七)现代 LLM 的核心架构设计其二:SwiGLU
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号