离散与连续傅里叶变换之间的联系
在复习这学期的数字信号处理的时候,发现我的数字信号处理书对此的解释我完全看不懂啊,所以就花了几天时间自己推了一点公式,通过了解这些离散啊,连续啊什么的傅里叶变换的数学性质能比较直观的帮助我理解这科目到底在干什么,具体为什么要这么干
我用到的这本书上说,DTFT在信号与系统中有学,所以略讲,信号与系统里提都没提DTFT 而且我觉得DTFT还非常容易和DFT混,如果不是真去尝试证明真不清楚其中是什么含义
DTFT
想要理解DTFT得先从实际的采样开始,采样就是每隔采样周期T取一个值,就是把一个连续的信号在时间域上离散化,或者用数学一点的说法,用这种方法把一个定义在
所以这里的连续信号中的连续并不是说的数学上的函数连续,而是指定义域不是
)
其表达式为
很容易就会想知道对
我们把两个变换的定义的式子写出来
首先要说明一下我这边的符号,比如
可以看出他俩其实很相似,所以我们很容易就能联想到下面这个等式
Lemma 1:
其中
Proof:
那么输入的函数关系怎么样已经找出来了,输出的频域函数两个变换又有什么关系呢?
首先再来算一个东西
Lemma 2:
Proof:
首先对
将这个式子代入连续傅里叶变换的定义式,得到
可以看到,对原本的连续信号的采样进行连续傅里叶变换得到的结果是原本信号的傅里叶变换的结果在频率轴上移动
好了,由上面得到的两个小等式联立一下可以证明我们需要的东西了
Equation:
将两边的
- 用人话翻译一下这个句子,对一个离散信号进行DTFT,等价于对它采样前的连续信号进行周期延拓(连续傅里叶变换后的结果在频率轴上移动
的整数倍得到的所有的无穷多个函数的叠加的 倍)后再把整个频谱函数沿着频率轴放大T倍 - 或者说你也可以理解为采样前的连续信号先沿频率轴放大T倍在进行周期延拓
- 再或者你还能理解为原连续信号的周期延拓完了之后把他沿频率轴放缩到以
为周期(DTFT它是关于 成周期的,你可以从这个等式里一目了然的看见,也能从DTFT的定义式中证明这个)
那么你就可以理解DTFT诞生的原因了,本质上它也是一种近似获得连续信号的傅里叶变换的手段,只不过算出来的结果是原本连续傅里叶变换的移位之后又压缩了一下再叠加嘛。
从这里你就知道什么叫混叠了,在一些比较良好的情况下,比如去DTFT一个带限信号的采样,带限信号是在频率上截止频率为
所以你现在明白了所谓奈奎斯特采样定律就是从采样后离散信号的DTFT能约等于看作是原本连续信号的CTFT的”约等于”能成立的需要的条件
DFT/FFT
FFT本质就是一个DFT的快速算法,就是另一种计算DFT的方法,和DFT完全没有区别,唯一需要注意的是要会用就行了,简单来说就是背下来那个蝶形图是怎么个画法
而DFT就更有趣了,先看看它的定义式
这里没有写成
它和DTFT只有一点点差别,首先它对
这下看懂了,对一个有限长序列进行DFT就是在对它进行DTFT后
还没完,如果你给这个长度为N的有限长序列
好嘛,就是没添加0之前的结果采样的更加密集了一点,只要你M取特别特别大,你就可以用DFT完全逼近这个有限长序列的DTFT了,例子如下
上图的matlab代码如下
1 | pic_num = 1; |
所以,求有限长序列的DFT也是约等于求它的DTFT,这样做的优势就是,DFT还能让结果也变得离散化,方便存储的同时还有FFT这种神奇的算法能给计算加速
总结
- 当带限信号满足奈奎斯特采样定律时,从采样后离散信号的DTFT能约等于看作是原本连续信号的CTFT
- 对于有限长序列,计算DFT就是在逼近计算DTFT函数的一部分值
所以综上,我所使用的这本教材前三章的核心目的就是告诉我,可以使用某种手段(比如滤波器等等),利用FFT就能近似计算一个连续函数的CTFT,至于其他的都是具体的计算时的分支任务。比如设计数字滤波器,本质上就是在给采样后的信号去干连续信号一样的事情
Gitalk 加载中 ...