族谱网 头条 人物百科

离散时间傅里叶变换

2020-10-16
出处:族谱网
作者:阿族小谱
浏览:796
转发:0
评论:0
定义一组离散的实数或复数:x[n](n为所有整数)的离散时间傅里叶变换是产生以频率为变量的周期函数的一个傅里叶级数。当频率变量ω的单位是归一化的弧度/样本时,周期为2π,而傅里叶级数为:此频率域函数的性质源于泊松求和公式(英语:Poissonsummationformula)。令X(f)为任意函数x(t)的傅里叶变换,采样间隔为T(秒),等价于序列x[n](或与之成正比),即T⋅⋅-->x(nT)=x[n]{\displaystyleT\cdotx(nT)=x[n]}。则以傅里叶级数表示的周期函数是X(f)的周期求和。赫兹以赫兹(周期/秒)为单位的频率f{\displaystyle\textstylef}的话就会是:图一.傅立叶变换(左上)和左下的其周期求和(DTFT)的图示。右下角显示了用离散傅里叶变换(DFT)计算DTFT的采样。整数k的单位为转/样本,采样频率是1/T,fs(样本/秒...

定义

一组离散的实数或复数: x [ n ]( n 为所有整数)的离散时间傅里叶变换是产生以频率为变量的周期函数的一个傅里叶级数。当频率变量 ω 的单位是归一化的 弧度/样本 时,周期为 2π,而傅里叶级数为:

此频率域函数的性质源于 泊松求和公式 ( 英语 : Poisson summation formula ) 。令 X(f) 为任意函数 x(t) 的傅里叶变换,采样间隔为 T( 秒 ),等价于序列 x[n](或与之成正比),即 T ⋅ ⋅ --> x ( n T ) = x [ n ] {\displaystyle T\cdot x(nT)=x[n]} 。则以傅里叶级数表示的周期函数是 X(f) 的周期求和。赫兹以赫兹( 周期/秒 )为单位的频率 f {\displaystyle \textstyle f} 的话就会是:

离散时间傅里叶变换

图一. 傅立叶变换(左上)和左下的其周期求和(DTFT)的图示。右下角显示了用离散傅里叶变换(DFT)计算DTFT的采样。

整数 k 的单位为 转/样本 ,采样频率是 1/ T , f s ( 样本/秒 )。因而 X 1/ T ( f ) 含有移位 f s 倍数赫兹了的 X ( f ) 的精确副本,并加和在一起。对于足够大的 f s ,可以在区间 [− f s /2, f s /2] 有很少或没有失真(混叠)地观察到 k = 0 项。在图1中,该左上角分布的末端在左下图中被周期求和的混叠遮盖住了。

我们还注意到 e − − --> i 2 π π --> f T n {\displaystyle e^{-i2\pi fTn}} 是 δ δ --> ( t − − --> n T ) {\displaystyle \textstyle \delta (t-nT)} 的傅里叶变换。因此,DTFT的另一个定义为:

调制的 狄拉克梳状函数 ( 英语 : Dirac comb ) 是一个数学抽象,有时被称为 脉冲采样 。

频谱的周期性与混叠

频谱周期性

F D T F T ( e i ω ω --> T ) {\displaystyle F_{DTFT}(e^{i\omega T})} 具有周期性:

显然有: F D T F T ( e i ω ω --> T ) = F D T F T ( e i ( ω ω --> + 2 π π --> / T ) T ) {\displaystyle F_{DTFT}(e^{i\omega T})\,\!=F_{DTFT}(e^{i(\omega +2\pi /T)T})}

频谱混叠

根据DTFT的定义,有

即, f(nT) 的DTFT是 f(t) 的傅里叶变换以Ω为周期的延拓,这也从另一个角度证明了DTFT的周期性。很显然,如果 f(t) 的频谱带不限于Nyquist间隔([-Ω/2, Ω/2]), f(nT) 的DTFT必然发生混叠(aliasing),如右图所示。混叠使得信号的低频部分被高频部分“污染”,造成信号的失真。为避免这种情况,通常在进行进一步的数字信号处理之前要对采样序列进行抗混叠滤波(anti-aliasing filtering),这一处理通常是由低通滤波器除去高频分量实现的。

DTFT与DFT

DFT(离散傅里叶变换)是对离散周期信号的一种傅里叶变换,对于有限长信号,则相当于对其周期延拓进行变换。在频域上,DFT的离散谱是对DTFT连续谱的等间隔采样。

DTFT与DFT频率分辨率

离散时间傅里叶变换

图2 DTFT与DFT。上图为10点DFT,下图为补零到12点的DFT。

N 点序列 f(n) (n=0, ... ,N-1)的DFT离散谱对应于对 f(nT) 连续谱(即DTFT)的N点采样,因此DFT的频率分辨率 Δ Δ --> ω ω --> = 2 π π --> / N {\displaystyle \Delta \omega =2\pi /N} 。为了提高频率分辨率,可以考虑增加在DTFT频域上的采样点数,对偶在时域就是增加对时域信号 f(n) 的采样数。对于有限长信号 f(n) ,在时刻 0 至 N-1 以外的值实际上是已知的——都为 0 。因此,只要在序列 f(n) 前后 补零 就增加了在时域的采样,假设在f(n)前后补上 M-N (其中 M>N )个零,则补零之后序列的DFT的频率分辨率就相应提高到 Δ Δ --> ω ω --> = 2 π π --> / M {\displaystyle \Delta \omega =2\pi /M} 。相关证明如下:

File:DFT 20 points.jpg 图3 20点采样的DFT能够分辨 f 1 与 f 2

通过上式可以清楚的看到, f(n) 补零之后的DFT增加了在 f(nT) 连续频域谱上的采样。采样点数从N增加到M,从而提高了DFT频谱的分辨率。另一方面,补零之后在频域采样的位置发生了变化,因此可以观察到其他的频点。

如图2所示,对信号 f ( t ) = 0.5 cos ⁡ ⁡ --> ( 2 π π --> f 1 t ) + cos ⁡ ⁡ --> ( 2 π π --> f 2 t ) {\displaystyle f(t)=0.5\cos(2\pi f_{1}t)+\cos(2\pi f_{2}t)} (其中 f 1 = 300 H z {\displaystyle f_{1}=300Hz} , f 2 = 250 H z {\displaystyle f_{2}=250Hz} )按照 f s = 1000 H z {\displaystyle f_{s}=1000Hz} 采样。 f 1 {\displaystyle f_{1}} 与 f 2 {\displaystyle f_{2}} 对应的数字角频率分别为 ω ω --> 1 = 0.6 π π --> {\displaystyle \omega _{1}=0.6\pi } , ω ω --> 2 = 0.5 π π --> {\displaystyle \omega _{2}=0.5\pi } 。上图为取10个采样点作DFT,可在0.6π处看到对应的频率分量,然而由于采样点少,看不到0.5π处的分量。下图为补2个零之后的DFT离散谱,可以见到离散谱的分辨率提高到了π/6,而且能够观察到10点DFT无法看到的0.5π频率分量。另外,虚线为DTFT连续谱,可见,DFT确实是在频域对DTFT的采样。

但是我们可以看到,即使是10点DTFT的连续谱也不能分辨 f 1 {\displaystyle f_{1}} 与 f 2 {\displaystyle f_{2}} (只有一个峰)。这是因为10点DTFT的分辨率为 f s /10 = 100Hz ,大于 f 1 - f 2 = 50 Hz 。所以,只有采样的点数超过20(即分辨率小于50Hz),才能分辨出 f 1 {\displaystyle f_{1}} 与 f 2 {\displaystyle f_{2}} 这两个频率分量(如图3所示)。而前面提到的对有限长信号补零作DFT以提高频谱分辨率的说法,也只是针对在DTFT连续谱上采样而言,只有增加采样点数和提高采样频率才能真正提高离散谱的分辨率。

以DFT近似DTFT

前面提到,在时间序列前后补零之后作DFT可以增加在DTFT上的采样点数。可以想见,如果补上无穷多个零,则可以得到无穷多个DTFT连续谱上的采样点,从而以DFT逼近DTFT。即,使得离散谱的分辨率足够小,即为连续谱。

这种想法显然是错误的,首先应该注意到,DFT的周期型,虽然进行DFT的序列x(k)是有限长序列,但是DFT隐含对x(k)的周期延拓,x(k)是信号序列,周期延拓后还是原来的信号,如果对x(k)进行补零后,1.补零后序列长度不是原信号序列长度的整数倍,那么补零后的序列就不再是原信号的序列,已经破坏了原信号的信息。2.补零后序列长度是原信号序列的整数倍,原信号信息没有被破坏,仍能通过DFT后获得原信号信息,但是补零后序列的DFT虽然频率步长小了,可以观察的频谱信息更丰富了,但是这些额外的频谱并不是原信号的频谱,而是窗函数的频谱信息,也就是包络线是窗函数的形状。所以补零并没有给出关于原信号频谱的更多信息。补零并不能使补零后的DFT结果的包络线和原信号频谱的包络线一致。

DTFT与Z变换

离散时间傅里叶变换可以被看作Z变换的特例。Z变换被定义为:

如果在z平面的单位圆( z = e i ω ω --> {\displaystyle z=e^{i\omega }} )上对信号 f ( n ) {\displaystyle f(n)} 做Z变换:

此即为 f ( n ) {\displaystyle f(n)} 的离散时间傅里叶变换。因此通常用 F ( e i ω ω --> ) {\displaystyle F(e^{i\omega })} ,而不是 F ( ω ω --> ) {\displaystyle F(\omega )} 表示DTFT。

参看

离散傅里叶变换

Z变换

傅里叶级数

采样定理

参考文献

Alan V. Oppenheim and Ronald W. Schafer : Discrete-Time Signal Processing , Prentice Hall Signal Processing Series, ISBN 978-0-13-754920-7

Boaz Porat : A Course in Digital Signal Processing , pp. 27-29, pp. 104-105, John Wiley and Sons, ISBN 978-0-471-14961-3

Sophocles J. Orfanidis : Introduction to Signal Processing , Prentice Hall International, Inc., ISBN 7-32-03059-6


免责声明:以上内容版权归原作者所有,如有侵犯您的原创版权请告知,我们将尽快删除相关内容。感谢每一位辛勤著写的作者,感谢每一位的分享。

——— 没有了 ———
编辑:阿族小谱

更多文章

更多精彩文章
评论 {{commentTotal}} 文明上网理性发言,请遵守《新闻评论服务协议》
游客
发表评论
  • {{item.userName}} 举报

    {{item.content}}

    {{item.time}} {{item.replyListShow ? '收起' : '展开'}}评论 {{curReplyId == item.id ? '取消回复' : '回复'}}

    回复评论
加载更多评论
打赏作者
“感谢您的打赏,我会更努力的创作”
— 请选择您要打赏的金额 —
{{item.label}}
{{item.label}}
打赏成功!
“感谢您的打赏,我会更努力的创作”
返回
打赏
私信

推荐阅读

· 离散傅里叶变换
定义对于N点序列{x[n]}0≤≤-->n<N{\displaystyle\left\{x[n]\right\}_{0\leqn傅里叶的离散傅里叶变换(DFT)为其中e{\displaystylee}是自然对数的底数,i{\displaystylei}是虚数单位。通常以符号F{\displaystyle{\mathcal{F}}}表示这一变换,即离散傅里叶变换的逆变换(IDFT)为:可以记为:实际上,DFT和IDFT变换式中和式前面的归一化系数并不重要。在上面的定义中,DFT和IDFT前的系数分别为1和1/N。有时会将这两个系数都改成1/N{\displaystyle1/{\sqrt{N}}}。从连续到离散连续时间信号x(t)以及对应的连续傅里叶变换x^^-->(ωω-->){\displaystyle{\hat{x}}(\omega连续函数都是连续函数。由于数字系...
· 离散时间信号
基本信号单位脉冲信号单位阶跃信号
· 离散傅里叶变换矩阵
定义N点的离散傅里叶变换可以用一个n××-->m{\displaystylen\timesm}的矩阵乘法来表示,即X=Wx{\displaystyleX=Wx},其中x{\displaystylex}是原始的输入信号,X{\displaystyleX}是经过离散傅里叶变换得到的输出信号。一个n××-->n{\displaystylen\timesn}的变换矩阵W{\displaystyleW}可以定义成W=(ωω-->ij)i,j=0,……-->,N−−-->1/N{\displaystyleW=(\omega^{ij})_{i,j=0,\ldots,N-1}/{\sqrt{N}}},或等效如下:其中ωω-->{\displaystyle\omega}是1{\displaystyle1}的n{\displa...
· 离散盘
形成阋神星,已知最大的黄道离散天体,和它的卫星阋卫一(迪丝诺美亚,中央偏左的小光点)。我们对离散盘的所知非常有限,虽然天文学的主流观点认为它是太阳系形成的早期过程中,因为海王星向外迁徙造成的引力扰动才被从柯伊伯带散射入高倾斜和高离心率的轨道内。相比之下,柯伊伯带像是一个相对“圆”和“平坦”的甜甜圈,以平和的圆轨道(QB1天体)和略为椭圆的共振轨道(类冥天体),将天体约束在30至44天文单位的圆环内;离散盘内的黄道离散天体轨道环境就比较怪异了。黄道离散天体,就以矮行星阋神星为例,在垂直黄道方向上的距离几乎和平行方向上与太阳的距离一样远;轨道模拟也显示黄道离散天体的轨道是怪异且不稳定的,并且最终会从太阳的核心区域抛至奥尔特云甚至更遥远的地方。有些迹象显示半人马群只是单纯的从柯伊伯带被向内抛射,而不是向外抛射的天体,可以称为“内海王星天体”(cis-Neptunianobject)。事实上,有些...
· 离散信号
参见采样采样定理信号连续信号离散时间信号

关于我们

关注族谱网 微信公众号,每日及时查看相关推荐,订阅互动等。

APP下载

下载族谱APP 微信公众号,每日及时查看
扫一扫添加客服微信