信号处理(6) 离散Fourier变换的进一步讨论
Nov 26, 2024
DFT和IDFT在实际应用中更可计算,因此FT计算数值解可以采取通过DFT作为跳板计算的方法
理论基础回顾
- 设连续时间信号f(t)的FT为F(ω)
- 对f(t)时域采样,得到数字域信号f(n)
- 于是,频域发生周期重复,采样信号f(t)n=−∞∑+∞δ(t−nTs)的FT变为((F(ω)/Ts))ωs,(⋅)T表示以T周期重复
- 数字域信号f(n)的DTFTX(ω)和((F(ω)/Ts))ωs值相等
- DFT谱X(k)由对X(ω)采样得到
讨论1:用DFT计算非周期信号的FT
设f(t)按Ts采样后的DFT为X(k),求其FT函数F(ω)的图像
根据上面的结论,X(k)是DTFT谱X(ω)的采样,因此
X(k)=X(ωk)
而
X(ω)=((F(ω)/Ts))ωs
于是
F(ωk)=TsX(k)
这里默认了周期重复不会对频谱值产生影响,即截止频率的存在
于是,只需要对TsX(k)连线,即得到F(ω)的近似图像
讨论2:用IDFT计算非周期信号频谱的IFT
已知F(ω)按采样频率ωs采样得到的IDFT(注:在计算IDFT时取的序列为F(ωk),其中ωk=kωs/N, k=0,1,⋯,N−1),求f(t)的近似图像
IDFT[F(ωk)]=IDTFT[((F(ω)))ωs]=Ts⋅IDTFT[((F(ω)/Ts))ωs]=Tsf(n)
注 上面三个等号的逻辑分别为:
- F(ωk)就是对F(ω)以ωs周期回绕后进行采样得到的,而因为DFT就是对DTFT进一步采样,因此由采样计算IDFT得到的信号和由原频谱计算IDTFT得到的信号不应该有所区别(事实上,由于二者都是从一个L长度的序列计算得到,因此虽然DTFT谱是连续的,但DTFT谱并不比DFT谱多出信息,只要DFT谱的长度N≥L)
- IDTFT是线性变换
- 根据DTFT[f(n)]=((F(ω)/Ts))ωs
这说明
f(n)=Ts1IDFT[F(ωk)]
即,仍然可以将IDFT看成是IFT的描点画图
讨论 讨论1和讨论2的两个公式有一个简单的记法
考虑到两个公式分别是用DFT表示FT和用IDFT表示IFT,而四者的量纲关系如下:
量纲 |
连续 |
离散 |
信号 |
IFT:Q/s |
IDFT:Q |
频谱 |
FT:Q/Hz=Q⋅s |
DFT:Q |
这里Q是一个表示强度的任意单位
根据上面的量纲分析很容易得出
FT=TsDFT
IFT=Ts1IDFT
讨论3:用DFT计算周期信号的FS
假设对周期信号进行频率为Nω1的采样,即一个周期内有N个抽样值,求N点DFT和原周期信号FS谱的关系
要求:周期信号满足采样定理,即Nω1要大于其频谱中的最大频率,即频谱只能出现在0, ±ω1, ⋯, ±(N−1)ω1位置
设周期信号FS系数为Fk
将周期信号取出一个周期,求其FT,设为F0(ω),则根据非周期信号FT与周期信号FS的关系可知
Fk=T11F0(kω1)
而由于
F0(kω1)=TsX(k)
于是
Fk=T1TsX(k)=N1X(k)