DFT和IDFT在实际应用中更可计算,因此FT计算数值解可以采取通过DFT作为跳板计算的方法

理论基础回顾

  • 设连续时间信号f(t)f(t)的FT为F(ω)F(\omega)
  • f(t)f(t)时域采样,得到数字域信号f(n)f(n)
  • 于是,频域发生周期重复,采样信号f(t)n=+δ(tnTs)f(t)\sum\limits_{n=-\infty}^{+\infty}\delta(t-nT_s)的FT变为((F(ω)/Ts))ωs((F(\omega)/T_s))_{\omega_s}()T(\cdot)_{T}表示以TT周期重复
  • 数字域信号f(n)f(n)的DTFTX(ω)X(\omega)((F(ω)/Ts))ωs((F(\omega)/T_s))_{\omega_s}值相等
  • DFT谱X(k)X(k)由对X(ω)X(\omega)采样得到

讨论1:用DFT计算非周期信号的FT

f(t)f(t)TsT_s采样后的DFT为X(k)X(k),求其FT函数F(ω)F(\omega)的图像

根据上面的结论,X(k)X(k)是DTFT谱X(ω)X(\omega)的采样,因此

X(k)=X(ωk)X(k)=X(\omega_k)

X(ω)=((F(ω)/Ts))ωsX(\omega)=((F(\omega)/T_s))_{\omega_s}

于是

F(ωk)=TsX(k)F(\omega_k)=T_sX(k)

这里默认了周期重复不会对频谱值产生影响,即截止频率的存在

于是,只需要对TsX(k)T_sX(k)连线,即得到F(ω)F(\omega)的近似图像

讨论2:用IDFT计算非周期信号频谱的IFT

已知F(ω)F(\omega)按采样频率ωs\omega_s采样得到的IDFT(注:在计算IDFT时取的序列为F(ωk)F(\omega_k),其中ωk=kωs/N, k=0,1,,N1\omega_k=k\omega_s/N,\ k=0,1,\cdots,N-1),求f(t)f(t)的近似图像

IDFT[F(ωk)]=IDTFT[((F(ω)))ωs]=TsIDTFT[((F(ω)/Ts))ωs]=Tsf(n){\rm IDFT}[F(\omega_k)]={\rm IDTFT}[((F(\omega)))_{\omega_s}]=T_s\cdot{\rm IDTFT}[((F(\omega)/T_s))_{\omega_s}]=T_sf(n)

上面三个等号的逻辑分别为:

  1. F(ωk)F(\omega_k)就是对F(ω)F(\omega)ωs\omega_s周期回绕后进行采样得到的,而因为DFT就是对DTFT进一步采样,因此由采样计算IDFT得到的信号和由原频谱计算IDTFT得到的信号不应该有所区别(事实上,由于二者都是从一个LL长度的序列计算得到,因此虽然DTFT谱是连续的,但DTFT谱并不比DFT谱多出信息,只要DFT谱的长度NLN\geq L
  2. IDTFT是线性变换
  3. 根据DTFT[f(n)]=((F(ω)/Ts))ωs{\rm DTFT}[f(n)]=((F(\omega)/T_s))_{\omega_s}

这说明

f(n)=1TsIDFT[F(ωk)]f(n)=\frac{1}{T_s}{\rm IDFT}[F(\omega_k)]

即,仍然可以将IDFT看成是IFT的描点画图


讨论 讨论1和讨论2的两个公式有一个简单的记法

考虑到两个公式分别是用DFT表示FT和用IDFT表示IFT,而四者的量纲关系如下:

量纲 连续 离散
信号 IFT:Q/s\rm Q/s IDFT:Q\rm Q
频谱 FT:Q/Hz=Qs\rm Q/Hz=Q\cdot s DFT:Q\rm Q

这里Q\rm Q是一个表示强度的任意单位

根据上面的量纲分析很容易得出

FT=TsDFT{\rm FT}=T_s{\rm DFT}

IFT=1TsIDFT{\rm IFT}=\frac{1}{T_s}{\rm IDFT}

讨论3:用DFT计算周期信号的FS

假设对周期信号进行频率为Nω1N\omega_1的采样,即一个周期内有NN个抽样值,求NN点DFT和原周期信号FS谱的关系

要求:周期信号满足采样定理,即Nω1N\omega_1要大于其频谱中的最大频率,即频谱只能出现在0, ±ω1, , ±(N1)ω10,\ \pm\omega_1,\ \cdots,\ \pm(N-1)\omega_1位置

设周期信号FS系数为FkF_k

将周期信号取出一个周期,求其FT,设为F0(ω)F_0(\omega),则根据非周期信号FT与周期信号FS的关系可知

Fk=1T1F0(kω1)F_k=\frac{1}{T_1}F_0(k\omega_1)

而由于

F0(kω1)=TsX(k)F_0(k\omega_1)=T_sX(k)

于是

Fk=TsT1X(k)=1NX(k)F_k=\frac{T_s}{T_1}X(k)=\frac{1}{N}X(k)