码农知识堂 - 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}

  • 相关阅读:
    电商数据的采集标准
    代码随想录算法训练营第三天| LeetCode ● 203.移除链表元素 ● 707.设计链表 ● 206.反转链表
    URL endoce 和 decode
    手把手教你实现buffer(四)——webrtc中的CopyOnWriteBuffer
    MySQL5.7安装详细过程--window系统
    关于c++中成员函数做友元的详细解释
    一篇文章教会你什么是高度平衡二叉搜索(AVL)树
    前端开发规范
    从零实现带RLHF的类ChatGPT:逐行解析微软DeepSpeed Chat
    目标检测 YOLOv5网络v6 0版本总结
  • 原文地址:https://blog.csdn.net/qq_40206371/article/details/126772853
  • 最新文章
  • 【FHE】我们如何实现同态加密推理(八):SiLU 的密文化——两条路径,和一个 8 字节的开关
    SQLcl 中 c/old/new 替换变量为何失效?先定位再替换
    不让模型写作文,直接从它脑子里读答案:Jev 决策在 .NET 的两条路线
    Halo插件:支持GitHub的图床附件策略
    电视投屏界的“一股清流”!SimpleTVCast 6.0.0-beta 抢先体验:无需TV端App,界面大换血!✨
    华三VXLAN集中式网关实战指南:Spine多活网关 + 双活接入(M-LAG)
    软件更好写之后口碑会更重要
    大 O、大 Omega、大 Theta 的区别,O(n) 复杂度是什么意思
    Java标签跳转:一个鲜为人知但却很实用的小技巧
    多智能体系统论文速读
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号