首页 > 其他分享 >有关DFT的条理不清晰理解

有关DFT的条理不清晰理解

时间:2023-11-09 10:12:00浏览次数:38  
标签:采样 frac 条理 DFT 频域 序列 pi 清晰

有关DFT的条理不清晰理解

本来今年暑假前就写好了,但写的md文件直接上传博客园一直不能正常显示公式,今天刚好解决就上传了一下,解决方法准备附在另外一篇博客

寄发现点进随笔还得刷新一下才能正常显示

以下内容是在对DFT以及其他变换有基础理解的基础上阐述,因此基本的符号解释和原理不再赘述

1. DFT、DTFT和DFS

首先是公式

\(DFS: \widetilde{X}[k]=\sum_{n=0}^{N-1} \widetilde{x}[n]e^{-j\frac{2\pi}{N}kn}\)

\(DTFT: X(e^{jw})=\sum_{n=0}^{N-1} x[n]e^{-jwn}\)

\(DFT: X[k]=\sum_{n=0}^{N-1}{x}[n]e^{-j\frac{2\pi}{N}kn}= X(e^{jw})\mid_{w=\frac{2\pi}{N}k}\)

其次是关系
  • DFT和DTFT

    假设对离散时间序列\(x[n]\)做DTFT后的结果是\(X(e^{jw})\),DFT相当于对\(X(e^{jw})\)在\([0,2\pi]\)的范围上以\(\frac{2\pi}{N}\)的间隔采样\(N\)个点得到\(X[k]\),即DFT是对离散信号进行傅里叶变换后频域采样的结果

    继续对上述DFT的结果\(X[k]​\)进行IDFT得到\(y[n]​\),其和原序列\(x[n]​\)的关系是,\(y[n]​\)为\(x[n]​\)以采样点数\(N​\)为周期进行周期延拓后截取主值(\(0\le n\le N-1​\))的结果,因此如果希望IDFT能够不失真的得到原序列,需要DFT点数不小于原序列长度(假设长度\(M​\)),即\(N \ge M​\)

  • DFT和DFS

    假设离散时间序列\(x[n]\)的长度为\(N\),\(\widetilde{x}[n]\)为\(x[n]\)以\(N\)为周期进行周期延拓,对\(\widetilde{x}[n]\)进行DFS后得到\(\widetilde X[k]\),截取主值(\(0\le k\le N-1\))后即为\(x[n]\)进行DFT的结果

2. 连续信号谱和连续信号采样谱的关系(为什么有T or 1/T)

先说结论
  1. 连续信号的频谱可以通过对连续信号采样,并进行 DFT 再对幅度乘以 \(T\)(采样周期) 的近似方法得到。
  2. 连续信号的时域采样信号可以通过对其频谱函数进行采样,并进行 IDFT 再对幅度乘以\(\frac{1}{T}\)的近 似方法得到。
推导过程和理解
  1. 假设一段连续信号(非周期)\(x_a(t)\), 持续时间为\(T_p\), 其最高频率为\(f_c\), \(\Omega\)为模拟角频率

    \(x_a(t) FT: X_a(j\Omega)=\int_{-\infty}^{+\infty} x_a(t)e^{-j\Omega t}dt=\int_{-\infty}^{+\infty} x_a(t)e^{-j2\pi ft}dt\)

  2. )对\(x_a(t)\)以\(f_s=\frac{1}{T}\ge 2f_c\)的频率进行采样,\(\omega=\Omega T\)为数字角频率,得到\(x[n]=x_a(nT)\),在持续时间\(T_p\)内共采\(N\)点,此时积分间隔\(dt=T\),且\(t=nT\) (这部分不确定能不能这么理解)

    \(x_a(nT) FT: X[j\Omega]=\sum_{n=0}^{N-1} x_a(nT)e^{-j\Omega nT}T=\sum_{n=0}^{N-1} x[n]e^{-j2\pi fnT}T\)

  3. 对\(X[j\Omega]\)在区间\([0,f_s]\)上等间隔采样\(N\)点,\(F=\frac{f_s}{N}=\frac{1}{NT}=\frac{1}{T_p}\)

    \(X[jkF]=\sum_{n=0}^{N-1} x[n]e^{-j2\pi kFnT}T=\sum_{n=0}^{N-1} x[n]e^{-j\frac{2\pi}{N} kn}T=T\bullet DFT(x[n])\)

\(T\)倍的形成可以理解成,时域采样相当于原连续信号和冲击序列相乘得到,相当于时域卷积,而冲击序列的傅里叶变换结果乘上了\(\frac{1}{T}\)的系数,因此需要对DFT结果再乘上\(T\)来抵消

3. 频谱泄露

主要现象

对连续信号进行采样后截断后进行DFT后,频谱中出现了原本信号中不存在的频率分量

比如用相同的频率采样,但截断不同的点数\(N\)后做FFT的结果,上图是非整数周期,下图是整数周期
image

对此我的理解是,原离散序列和对离散序列进行DFT后再进行IDFT的关系是,后者是原序列以频域采样点数\(N\)为周期周期延拓后截取主值(此处假定离散序列的长度和DFT点数相同)

那么如果离散序列包含非整数周期,进行周期延拓后就会导致拼接不上(最简单的正弦信号采样成离散序列,只取1.5个周期,之后周期延拓显然拼接不上),从而产生时域上的突变,频域带宽增大,最终表现为频域上有效频点向周围“泄露”

4.栅栏效应

DFT是在频域以点的形式进行采样,就行从N个栏栅缝里去观察频域信号,其余看不见的地方显然会丢失信息,导致有点频率无法分辨,常见的方法是补0

首先是两种分辨率定义

\(\Delta R_w=\frac{1}{T_p} \\\Delta R_{fft}=\frac{F_S}{N_{fft}}\)

前一个是波形分辨率,\(T_p=\frac{N}{F_S}\)表示时域采样的时长,\(N\)为采样长度,后者是FFT分辨率

  • 时域采样信号补零相当于频域插值,能让FFT结果更加光滑、采样频点更加密集,但无法改变波形分辨率,该丢失的频点还是丢失

  • 增加时域采样时长\(T_p\),也就是增加原序列长度,可以分辨率增大,即\(\Delta R_w\)减小,从而分辨出两个相近的频点

  • 增加频域采样点数\(N_{fft}\),增大FFT分辨率,可以增大波形分辨率

标签:采样,frac,条理,DFT,频域,序列,pi,清晰
From: https://www.cnblogs.com/Dmmmy/p/17819057.html

相关文章

  • WOA-CNN基于鲸鱼算法优化卷积神经网络的多变量回归预测 可直接运行 注释清晰适合新手
     ✅作者简介:热爱科研的Matlab仿真开发者,修心和技术同步精进,代码获取、论文复现及科研仿真合作可私信。......
  • 属性很多时如何很清晰的通过充血模型创建数据合法的对象
    原因:在构建充血模型时,为构建合法的对象,对象赋值都在私有的。属性少的时候可以直接写构造函数,属性多时就不太合适了。如何解决这个问题呢?直接上代码Book类:1publicclassBook2{3publiclongId{get;privateset;}4publiclongAuthorId{g......
  • DDD能够帮我们设计出清晰的领域和应用边界
    DDD能够帮我们设计出清晰的领域和应用边界DDD包括战略设计和战术设计两部分。战略设计主要从业务视角出发,建立业务领域模型,划分领域边界,建立通用语言的限界上下文,这些限界上下文可以作为微服务设计的参考边界。而战术设计则从技术视角出发,着重于领域模型的技术实现,包括聚合根、......
  • 半导体测试DFT
    如下图所示,一颗芯片最终做到终端产品上,一般需要经过芯片设计、晶圆制造、晶圆测试、封装、成品测试、板级封装等这些环节。 在整个价值链中,芯片公司需要主导的环节主要是芯片设计和测试,其余的环节都可以由相应的partner来主导或者完成。 半导体测试定义与基本工作机制:半导体......
  • Unity完美像素Sprite:怎么让图片变得清晰(转载) Unity Pixel Perfect Sprite: How To Ach
    https://gamedevelopertips.com/unity-pixel-perfect-sprite/SoIwasmakingalittleprototypeformynewgamewhenIjustcameacrossalittleproblem.ThespritethatIwasloadingintoUnitywasjustlookingtoopixelatedandnotsharpatall.SoIspen......
  • VISIO画的论文图如何清晰地放在LATEX论文中?
    1、将VISIO图Microsoftprinttopdf2、利用AdobeAcrobat进行裁剪3、(optional)如果visio图中有图片,步骤1的pdf中图片可能存在不清晰,利用AdobeAcrobat编辑替换掉不清晰的图片4、(optional)如果进行了第3步,还需要再进行Microsoftprinttopdf一次,否则步骤5中的图片仍然不清晰5、......
  • 针对直播痛点的关键技术解析——首帧秒开、清晰度、流畅度
    ​一、直播体验升级的一把利刃——「畅直播」良好的直播体验一直是直播产品留住用户的关键流量密码。国内直播行业已呈现一片红海,各直播产品的用户增长速度在逐渐放缓,行业内竞争进入白热化阶段,急需通过在直播体验上的精细化运营来保障平台营收。如针对重大活动、头部主播进行......
  • 数字信号处理-DFT
    1.线性卷积和循环卷积线性卷积可由循环卷积计算而来,循环卷积可由DFT计算,DFT可由FFT进行计算。当量序列长度相当的时候,可以通过循环卷积来减少线性卷积的计算量。如果序列长度相差较大,直接计算线性卷积的计算量会更少。举例求x[k]和h[k]的卷积:......
  • winform 解决 窗体清晰度问题
    在C#中让Winform应用程序窗口自动适应高DPI缩放不使用其他特殊组件库的情况下添加清单文件在末尾添加:<applicationxmlns="urn:schemas-microsoft-com:asm.v3"><windowsSettings><dpiAwarexmlns="http://schemas.microsoft.com/SMI/2005/WindowsSettings">true&l......
  • DFT 在信号频谱分析中的应用
    DFT在信号频谱分析中的应用实验目的熟悉DFT的性质。DFT是离散傅里叶变换的缩写,是一种将时域信号转换为频域信号的数学工具。下面是DFT的一些基本性质:线性性:DFT是线性的,即它满足叠加原理。如果x1(n)和x2(n)是两个长度为N的离散时间信号,那么它们的DFT可以表示为:X(k)=D......