傅立叶变换家族中的关系

来自testwiki
跳转到导航 跳转到搜索

Template:NoteTA Template:傅里叶变换 在数学领域的谐波分析中,连续傅里叶变换(continuous Fourier transform, CFT)与傅里叶级数 (Fourier series, FS)有非常微妙的关系。而且连续傅里叶变换也与离散时间傅里叶变换(discrete time Fourier transform, DTFT)和离散傅里叶变换(discrete Fourier transform, DFT)有很近的关系。傅里叶变换家族通常就是指这四种变换。

通过利用Dirac delta函数 δ(t) ,CFT可以应用到时间离散 (time-discrete)或时间周期(time-periodic)信号。实际上,FS、 DTFT和DFT都可以由最广泛的CFT得到。从理论上看,它们也都是CFT的特殊情况。

在信号理论和数字信号处理(digital signal processing, DSP)中,DFT扩展用于近似计算连续信号的频谱,其变换的对象只是一个采样点的有限序列,而且可以由快速傅里叶变换(fast Fourier transform, FFT)实现。

家族中各个变换的定义

下表中左上、左下、右上和右下分别对应了傅里叶变换家族中CFT、FS、DTFT和DFT四个变换对的定义。

傅里叶变换家族中各种变换的定义
× 连续时间 离散时间
时间非周期 x(t)=X(f) ei2πftdf x[n]=Ts1/TsX¯(f) ei2πfnTs df
- X(f)=x(t) ei2πftdt X¯(f)=n=+x[n] ei2πfnTs
时间周期 x¯(t)=k=+X[k]ei2πkT0t xn=1Nk=0N1Xkei2πNkn,n=0,,N1.
- X[k]=1T0T0x¯(t)ei2πkT0tdt Xk=n=0N1xnei2πNkn,k=0,,N1.

显然,上表是从时域信号的角度来划分的:表的列区分了连续时间和离散时间的信号,而表的行则区分了时间上非周期的信号和时间上周期的信号。其中重要的参量符号解释为:

  • x[n]X[k] 都为无限序列,其采样间隔,即间隔时间和间隔频率分别为 Tsf0=1/T0
  • x¯(t)X¯(f) 都为周期函数,且时间周期和频率周期分别为 T0fs=1/Ts
  • xnXk 都为有限序列,且序列长度都为 N

关系推导所需的公式

前面表中的定义都可以通过Dirac delta函数 δ(t) 的扩展形式 ,即Dirac comb函数,由CFT引入或推导。为计算离散和/或周期信号的CFT,我们需要引入一些公式,并使用傅里叶变换的一些特性。以下集中给出:

1. Dirac comb函数的傅里叶变换

Dirac comb函数的定义为

ΔT(t)=defn=δ(tnT)

在电气工程中通常又称作冲击串(impulse train)或采样函数 (sampling function)。其重要的傅里叶变换为:

n=δ(tnT)=1Tk=ei2πkTt1Tk=δ(fkT)=n=ei2πnTf

这个变换在傅里叶变换家族中各个变换之间转换上起关键作用。

2. 傅里叶变换的卷积定理(convolution theorem)

这包括了傅里叶变换的时域卷积和频域卷积:

x1(t)x2(t)X1(f)X2(f)x1(t)x2(t)X1(f)X2(f)

3. 泊松求和公式(Poisson summation formula)

由Dirac comb函数的傅里叶变换和卷积定理,容易证明泊松求和公式:

1.n=x(tnT0)=1T0k=X(kT0)ei2πkT0t2.n=x(nT)ei2πnTf=1Tk=X(fkT)

若第1和第2公式中分别取 t=0f=0 ,得到相同等式:

n=x(nT)=1Tk=X(kT)

这表明,傅里叶变换时时域函数 x(t) 和频域函数 X(f) 分别以 T1/T 为间隔采样,则所有时域采样点的总和与所有频域采样点扩大 1/T 的总和相等。

各种变换之间的关系

图 1. 此“立方体”图形表示了连续傅里叶变换离散时间傅里叶变换傅里叶级数离散傅里叶变换之间的关系。


图中立方体包含了频域和时域两个平面上各种变换的关系,同时两平面相连的四个边则分别代表了CFT、FS、DTFT和DFT。其中参量符号与前面表中相同,另外增加:

  • X~k 为由FS和DTFT推导DFT得到的DFT'频域形式,与传统DFT的频域 Xk 有关系: Xk=T0X~k
  • 图中粗的双箭头()表示每个函数和其变换之间的联系;

总的说来,各种变换之间的转换是一个周期扩展或采样的过程:

  • 如果时域进行周期扩展,则频域为采样;如果时域进行采样,则频域为周期扩展;
  • 一个转换中,周期扩展的周期与采样的间隔有倒数关系;
  • 频域的周期扩展或者采样,都有一个周期或采样间隔作系数;

这里的周期扩展就是与Dirac comb函数相卷积,而采样则是与Dirac comb函数相乘。

从CFT分别到FS和DTFT的转换都容易推导,下面具体说明FS和DTFT到DFT/DFT'转换的推导,最后说明连续FT与DFT/DFT'的关系。

由DTFT推导DFT

设DTFT,及对应的CFT为:

x[n]𝒟𝒯𝒯X¯(f)n=x[n]δ(tnT)X¯(f)

在时域作周期为 T0 的扩展,有:

(n=x[n]δ(tnT))(n=δ(tnT0))=n=i=x(nT)δ(tnTiNT)=n=(i=x(nTiNT))δ(tnT)

其中代入了 T0=NT ,而由于 ni 的求和区间都为 ,可以用 niN 代替 n 得到最后一步推导。取:

xn=i=x(nTiNT)=i=x[niN]=i=x(nTiT0)

在频域作带系数 1/T0 且间隔也为 1/T0 的采样,有:

X¯(f)1T0k=δ(fkT0)=1T0T(k=X(fkT))(k=δ(fkT0))=1T0Tk=(i=X(fiT)δ(fkT0))=1T0Tk=(i=X(kT0iT))δ(fkT0)

取:

X~k=1T0Ti=X(kT0iT)=1T0X¯(kT0)


由FS推导DFT

设FS,及对应的CFT为:

x¯(t)𝒮X[k]x¯(t)k=X[k]δ(fkT0).

在时域作间隔为 T 采样,有:

x¯(t)n=δ(tnT)=(n=x(tnT0))(n=δ(tnT))=n=(i=x(tiT0)δ(tnT))=n=(i=x(nTiT0))δ(tnT)

取:

xn=i=x(nTiT0)=x¯(nT)

在频域作带系数 1/T 且周期也为 1/T 扩展,有:

(k=X[k]δ(fkT0))(1Tk=δ(fkT))=1T0Tk=i=X(kT0)δ(fkT0iNT0)=1T0Tk=(i=X(kiNT0))δ(fkT0)

其中也代入了 T0=NT ,而由于 ki 的求和区间都为 ,可用 kiN 替代 k 得到最后一步推导。 取:

X~k=1T0Ti=X(kiNT0)=1Ti=X[kiN]=1T0Ti=X(kT0iT)


CFT与DFT的关系

前面FS到DFT和DTFT到DFT的推导都得到相同的 xnX~k 。这里的 xnX~k 可看作一种DFT变换对,有关系:

X~k=1T0n=0N1 xne2πiNknxn=Tk=0N1 X~ke2πiNkn

记为:

xn𝒟𝒯X~k

对比传统DFT变换对的 xnXk ,显然有:

Xk=T0X~k.

这一对变换的等式右边系数的乘积为 T/T0=1/N ,符合我们在DFT中的说明,因而完全可以将这里的DFT'看作传统DFT的另一种变换形式 。

而由前面转换的推导过程可得到:

n=xnδ(tnT)k=X~kδ(fkT0)=1T0k=Xkδ(fkT0)

为一对CFT,其中要求 T0=NT 。加之如果 x(t)X(f),则有:

xn=i=x(nTiT0)𝒟𝒯X~k=1T0Ti=X(kT0iT)𝒟𝒯Xk=1Ti=X(kT0iT)

其中可以任选 T0T 。这样就建立了CFT和DFT之间的双向关系。但应注意到,此时我们已经将DFT'和DFT都做了周期拓展,即 n,k

参看

参考文献

  1. Oppenheim, Alan V.; Schafer, R. W.; and Buck, J. R., (1999). Discrete-time signal processing, Upper Saddle River, N.J. : Prentice Hall. ISBN 0137549202
  2. Sklar, B., (2001). Digital Communications: Foundamentals and Applicatons, 2nd Edition, Prentice Hall PTR. ISBN 0130847887

en:Fourier analysis#Variants of Fourier analysis