1、傅里叶变换和傅里叶级数的收敛问题
由于傅里叶级数是一个无穷级数,因而存在收敛问题。这包含两方面的意思:是否任何周期信号都可以表示为傅里叶级数;如果一个信号能够表示为傅里叶级数,是否对任何t值级数都收敛于原来的信号。关于傅里叶级数的收敛,有两组稍有不同的条件。
第一组条件:如果周期信号x?t?在一个周期内平方可积,即
?2T0x?t?dt??
则其傅里叶级数表达式一定存在。
第二组条件,与第一组条件稍有不同,就是狄里赫利条件,它包括以下三点: (1)在任何周期内,x必须绝对可积,即
?x?t?dt??
T0(2)在任何周期内,x?t?只有有限个极值点,且在极值点处的级值为有限值。 (3)在任何有限区间内,x?t?只有有限个间断点,且在这些不连续点处,x?t?为有限值。
傅里叶变换的收敛问题也有两组类似的条件: 第一组条件:如果x?t?平方可积,即
?x?t????2dt??
则x?t?的傅里叶变换存在。满足上式可以保证X???为有限值。 第二组条件也称为狄里赫利条件,这就是: (1)x?t?绝对可积,即
?x?t?dt??
???(2)在任何有限区间内,x?t?只有有限个极值点,且在这些极值点处的极值是有限值。
(3)在任何有限区间内,x?t?只能有有限个间断点,而且这些间断点都必须是有限值。
吉布斯现象:
当简单地把信号频谱截断时,相当于给信号频谱加上了一个矩形窗口函数,正是由于矩形窗口函数的时域特性导致了在间断点处的吉布斯现象的产生。
2、周期序列的傅里叶级数展开和傅里叶变换之间的问题
假定x?t?是一个长度为N的有限长序列,将x?t?以N为周期延拓而成的周期序列为x?n?,则有
x?n??~~~~r????x?n?rN?
?或表示为x?n??x??n??N。于是x?n?与x?n?的关系表示为:
x?n??x??n??N x?n??x?n?RN?n?
~~将x?n?表示为离散时间傅里叶级数有:
1N?1~?kn x?n???X?k??WNNn?0~kn X?k???x?n??WNn?0~N?1~~其中X?k?是傅里叶级数的系数,这样做的目的是使其表达形式与离散时间傅里叶变换的形式相类似。如果将X?k?的主值周期记为X?k?,0?k?N?1,由于以上两式中的求和范围均取为区间0~N-1,在次区间内x?n?=x?n?,因此可以得到:
kn, 0?k?N?1 X?k???x?n?WNn?0N?1~~~~1x?n??N?X?k?Wn?0N?1~?knN, 0?n?N?1
表明时域N点有限长序列x?n?可以变换成频域N点有限长序列X?k?。显然,DFT与DFS之间存在以下关系:
X?k??X??k??N X?k??X?k?RN?k?
~~
3、频率分辨率的问题
若信号最高频率为fh,按抽样定理,抽样频率应满足
fs?2fh
也就是抽样间隔为T满足
T?11 ?fs2fh一般取
fs??2.5~3.0?fh
如果不满足fs?2fh的要求,就会产生频率响应的周期延拓分量相互重叠的现象,也就是产生频率响应的混叠失真。
对于DFT来说,频率函数也要抽样,变成离散的序列,其抽样间隔为F0,这就是我们能得到得的频率分辨力,有它可引出时间函数的周期,也就是所取的记录长度T0为
T0?1 F0从以上T和T0两个公式来看,信号的最高频率分量fh与频率分辨力F0之间有着矛盾关系,要想fh增加,则时域抽样间隔T就一定减小??T?f?2f??,而fssh??(抽样频率)就增加,由于抽样点数满足
fsT0??N F0T?11?则此时fs增加,若是N固定的情况下,必然要F0增加,即分辨率下降。 反之,要提高频率分辨力(减少F0),就要增加T0,当N给定时,必然导致T的增加(fs减小)。要不产生混叠失真,则必然会减小高频容量(信号的最高频率分量)fh。
要想兼顾高频容量fh与频率分辨力F0,即一个性能提高而另一个性能不变(或也得以提高)的惟一办法就是增加记录长度的点数N,即要满足
N?fs2fh ?F0F0这个公式是未采用任何特殊数据处理(例如加窗处理)的情况下,为实现基
本DFT算法所必须满足的最低条件。如果加窗处理,相当于时域相乘,则频域卷积,必然加宽频谱分量,频率分辨力就可能变坏,为了保证频率分辨力不变,则须增加记录长度,也就是增加数据长度T0。
4、MATLAB的图示说明:有效观察时间与补零后的DFT之间的关系,以及与DTFT之间的关系
对8点正弦离散序列求8点、32点和64点DFT,观察频域变化(分别用绿、黄、红色表示)。结果:
矩形窗序列后补零的时、频域示意图
从图中可以看出:序列后补零可以降低栅栏效应;信号频谱的形状只取决于时域信号,与补零个数无关。补零并不能提高频谱分辨率,因为频谱分辨率只与时域数据的有效长度有关。
DTFT与DFT(或DFS)的关系:DFT时域序列为周期序列,周期为N;频域序列也是周期序列,周期也是N点。当N不断增大时,频域包络不变,但谱线变密;显然,N??时,时域序列变为非周期序列,频域为连续的频谱,即变化为DTFT。
5、教材
《信号与线性系统》,阎鸿森、王新凤、田惠生编,西安交通大学出版社 《数字信号处理教程》,程佩青编,清华大学出版社
(后附连续信号傅里叶变换的DFT近似计算)
傅里叶变换的DFT近似计算
连续时间非周期信号x?t?的傅里叶变换对为
X?j????x?t?e?j?tdt (1)
???1x?t??2?????X?j??ej?td? (2)
用DFT方法计算这一对变换的方法如下:
(1)将x?t?在t轴上等间隔(宽度为T)分段,每一段用一个矩形脉冲代替,脉冲的幅度为其起始点的抽样值x?t?t?nT?x?nT??x?n?,然后把所有矩形脉冲的面积相加。由于
t?nT
dt?T ?dt??n?1?T?nT?
?则得频谱密度X?j??????dt?n????T
??????x?t?e?j?tdt的近似值为
?X?j????x?nT??e?j?nT?T (3)
??(2)将序列x?n??x?nT? 截断成从t?0开始长度为t0的有限长序列,包含有N个抽样(即时域取N个样点),则上式成为
X?j???T?x?nT??e?j?nT (4)
n?0N?1由于时域抽样,抽样频率为fs?1T,则频域产生以fs为周期的周期延拓,如果频域是限带信号,则可能不产生混叠,成为连续周期频谱序列,频域周期为fs?1T(即时域的抽样频率)。
(3)为了数值计算,再频域上也要离散化(抽样)即在频域的一个周期(fs)中也分成N段,即取N个样点fs?NF0,每个样点间的间隔为F0。频域抽样,那么频域的积分式(2)式就变成求和式,而时域就得到原已截断的离散时间序列的周期延拓序列,其时域周期为T0?1F0。这时??k?0
d???k?1??0?k?0??0