首页 > 其他分享 >03. 离散傅里叶变换

03. 离散傅里叶变换

时间:2023-11-30 16:14:43浏览次数:35  
标签:03 frac limits sum 离散 tag tilde pi 傅里叶

离散傅里叶级数

对一个周期为N的序列,其离散傅里叶级数有:

\[\tilde{x}(n) = \frac{1}{N}\sum\limits_{k=0}^{N-1}\tilde{X}[k]e^{j\frac{2\pi}{N}kn} \tag{1.1} \]

两边同时乘以\(e^{-j \frac{2\pi}{N}rn}\),并累加有:

\[\sum\limits_{n=0}^{N-1}\tilde{x}(n)e^{-j \frac{2\pi}{N}rn} = \sum\limits_{k=0}^{N-1}\sum\limits_{n=0}^{N-1} \frac{1}{N}\tilde{X}(k)e^{j \frac{2\pi}{N}n(k-r)} \tag{1.2} \]

有:

\[e^{j \frac{2\pi}{N}n(k-r)} = \left \{ \begin{matrix} N &k-r=m2\pi \\ 0 \end{matrix} \right. \tag{1.3} \]

所以有:

\[\sum\limits_{k=0}^{N}\tilde{X}[r] = \sum\limits_{n=0}^{N-1}\tilde{x}(n)e^{-j \frac{2\pi}{N}rn} \tag{1.4} \]

即:

\[\tilde{X}[k] = \sum\limits_{n=0}^{N-1}\tilde{x}(n)e^{-j \frac{2\pi}{N}kn} \tag{1.5} \]

旋转因子\(W_{N}^{kn} = e^{-j \frac{2\pi}{N}kn}\)
所以有周期函数的离散傅里叶级数:

\[\tilde{x}(n) = \frac{1}{N}\sum\limits_{n=0}^{N-1}\tilde{X}[k]W_{N}^{kn} \tag{1.6} \]

\[\tilde{X}[k] = \sum\limits_{n=0}^{N-1}\tilde{x}(n)W_{N}^{-kn} \tag{1.7} \]

周期函数的离散时间傅里叶变换:
有结论:

\[\tilde{X}(e^{j\omega}) = DTFT[\tilde{x}(n)] = \sum\limits_{k=-\infty}^{\infty} \frac{2\pi}{N}\tilde{X}[k]\delta\left(\omega - \frac{2\pi}{N}k\right) \tag{1.8} \]

证明:只需要将两边同时进行IDTFT即可得到相等。

假设去上述周期函数的一个主值函数\(x(n)\),有:

\[\tilde{x}(n) = x(n)*\tilde{p}(n) \tag{1.9} \]

其中\(\tilde{p}(n) = \sum\limits_{r=-\infty}^{\infty}\delta(n-rN)\)
并且:

\[\tilde{P}[k] = \sum\limits_{k=0}^{N-1}\sum\limits_{r=-\infty }^{\infty}\delta(n-rN)W_{N}^{kn} = 1 \tag{1.10} \]

根据上面(1.8)结论有:

\[\tilde{P}(e^{j\omega}) = \sum\limits_{k=-\infty}^{\infty} \frac{2\pi}{N}\delta(\omega- \frac{2\pi}{N}k) \tag{1.11} \]

由时域卷积定理有:

\[\tilde{X}(e^{j\omega}) = \sum\limits_{k=-\infty}^{\infty} \frac{2\pi}{N}X(e^{j \frac{2\pi}{N}k})\delta(\omega - \frac{2\pi}{N}k) \tag{1.12} \]

对比上面两个式子(1.8),(1.12)有:

\[\tilde{X}[k] = X\left(e^{j \frac{2\pi}{N}k}\right)= X(e^{j\omega})|_{\omega= \frac{2\pi}{N}k} = X(z)|_{z= e^{j \frac{2\pi}{N}k}} \tag{1.13} \]

从上面式子得出离散傅里叶级数就是对离散时间傅里叶变换的采样。

离散傅里叶变换

对一个有限长序列\(x(n)\),我们可以将其进行时域周期延拓,并且其在频域有对应的周期离散傅里叶级数,因此定义离散傅里叶变换如下:

\[X[k] = DFT[x(n)] = \sum\limits_{k=0}^{N-1}x(n)W_{N}^{kn} \tag{1.14} \]

性质

和其他傅里叶变换性质类似,只不过注意旋转因子中已经包含有一个负号。

意义

  我们已经知道离散时间傅里叶变换的频域是连续的,可能(我目前也不知道有啥好处,好像后面FFT有用,汗。。)会带来些影响。所以我们引入离散傅里叶变换,通过离散傅里叶变换,我们将原本的数字序列变成同样离散的频域信息,根据这些信息,我们仍能复原出我们的离散数字序列.

标签:03,frac,limits,sum,离散,tag,tilde,pi,傅里叶
From: https://www.cnblogs.com/C-qian/p/17867605.html

相关文章

  • 查看一个 '函数' 到底是【函数】还是【方法】
    需要用到FunctionType和MethodTypefromtypesimportFunctionType,MethodType 准备测试用代码写一个函数:defadd():pass写一个类:classPerson:defrun(self):pass@classmethoddefxx(cls):pass@staticmethod......
  • 重装vs2022 nuget添加包报错: Unexpected character encountered while parsing value:
    工具--》选项--》Nuget包管理器,点击清除所有Nuget存储 参考文献:关于VSNuGet包无法更新,设置包源映射无效的问题-CSDN博客         微软官方文献 ......
  • 达梦数据库,使用存储过程切分','分隔字段,生成多条数据
      CREATEORREPLACEPROCEDUREsplit_listASv_nameVARCHAR2(255);v_listVARCHAR2(255);v_valueVARCHAR2(255);v_posNUMBER;v_lengthNUMBER;BEGINFORrecIN(SELECTZLY_NAME,LIST_IDFROMyour_table)LOOPv_name:=re......
  • linux awk 多分隔符 -F'[ :,=]'
    摘自:https://blog.csdn.net/whatday/article/details/90415997awk的-F参数可以指定新的记录分隔符,有些时候可能需求指定多个分隔符,比如下面的内容width:720height:360如果需要取出width和height后面的值的话,一般大家会这样做,即做两次awk操作#echo"width:720height:360"......
  • day03
    day031.今日内容介绍1.编程语言介绍编程语言分类,总结2.python介绍3.解释器多版本共存4.运行python程序的两种方式5.一个python应用程序的运行的三个步骤(*****)对比文本编辑器运行的三个步骤6.注释7.IDE集成开发环境pycharm8.虚拟环境的简单使用9.注......
  • (查找)03-寻找峰值
    1importjava.util.*;23publicclassSolution{4/**5*@paramnumsint整型一维数组6*@returnint整型7*/8publicintfindPeakElement(int[]nums){9//申请左指针10intleft=0;11//申......
  • 麻烦问一下Python采集到的文本列表中有大量的 ', ' 符号 想这种符号怎么删除
    大家好,我是皮皮。一、前言前几天在Python铂金流群【泅渡】问了一个Python字符处理的问题,一起来看看吧。问题描述:麻烦问一下Python采集到的文本列表中有大量的  ','  符号 想这种符号怎么删除?二、实现过程这里【不上班能干啥!】和【瑜亮老师】分别给了一个指导,如下......
  • # yyds干货盘点 # 麻烦问一下Python采集到的文本列表中有大量的 ', ' 符号 想这种符号
    大家好,我是皮皮。一、前言前几天在Python铂金流群【泅渡】问了一个Python字符处理的问题,一起来看看吧。问题描述:麻烦问一下Python采集到的文本列表中有大量的  ','  符号 想这种符号怎么删除?二、实现过程这里【不上班能干啥!】和【瑜亮老师】分别给了一个指导,如下图所示:......
  • #P1031. 三连击
    next_permutation全排列公式数据不大,直接暴力全排列模拟出数据即可usingnamespacestd;intp[9]={1,2,3,4,5,6,7,8,9};intmain(){ inta,b,c; cin>>a>>b>>c; intcnt=0; do{ intx=p[0]*100+p[1]*10+p[2]; inty=p[3]*100+p[4]*10+p[5]; intz=p[6]*100+p[7]*10+......
  • FileNotFoundError: [WinError 2] 系统找不到指定的文件。: '0054243eb93327df4b59023
    importos#指定目录directory='E:\\pythonProject\\a'#获取当前目录下所有图片文件image_files=[fforfinos.listdir(directory)iff.endswith('.jpg')orf.endswith('.png')orf.endswith('.jpeg')]#重命名图片文件fori,fileinenumer......