傅里叶变换和傅里叶级数的收敛问题
1、傅里叶变换和傅里叶级数的收敛问题
由于傅里叶级数是一个无穷级数,因而存在收敛问题。这包含两方面的意思:是否任何周期信号都可以表示为傅里叶级数;如果一个信号能够表示为傅里叶级数,是否对任何t值级数都收敛于原来的信号。关于傅里叶级数的收敛,有两组稍有不同的条件。
第一组条件:如果周期信号在一个周期内平方可积,即 ,,xt
2
x,,tdt,, ,T0
则其傅里叶级数表达式一定存在。
第二组条件,与第一组条件稍有不同,就是狄里赫利条件,它包括以下三点:
(1)在任何周期内,x必须绝对可积,即
,,xtdt,, ,T0
(2)在任何周期内,只有有限个极值点,且在极值点处的级值为有限值。 ,,xt
(3)在任何有限区间内,只有有限个间断点,且在这些不连续点处,,,,,xtxt为有限值。
傅里叶变换的收敛问题也有两组类似的条件:
第一组条件:如果平方可积,即 ,,xt
,2,,xtdt,, ,,,
则的傅里叶变换存在。满足上式可以保证为有限值。 ,,,,X,xt
第二组条件也称为狄里赫利条件,这就是:
(1)绝对可积,即 ,,xt
,
,,xtdt,, ,,,
(2)在任何有限区间内,只有有限个极值点,且在这些极值点处的极值,,xt
是有限值。
(3)在任何有限区间内,只能有有限个间断点,而且这些间断点都必须,,xt
是有限值。
吉布斯现象:
当简单地把信号频谱截断时,相当于给信号频谱加上了一个矩形窗口函数,正是由于矩形窗口函数的时域特性导致了在间断点处的吉布斯现象的产生。
- 1 -
2、周期序列的傅里叶级数展开和傅里叶变换之间的问题
假定是一个长度为N的有限长序列,将以N为周期延拓而成的周期,,,,xtxt
~
序列为,则有 ,,xn
,~
,,,,xn,xn,rN ,r,,,
~~
或表示为。于是与的关系表示为: ,,,,,,xn,xn,,,,xnxnN
~
,,,,,,xn,xnN
~
,,,,,,xn,xnRnN~
将表示为离散时间傅里叶级数有: ,,xn
N,1~~1,kn,,,,xn,Xk,W ,NNn,0
,1N~~kn,,,,Xk,xn,W ,Nn,0
~
其中是傅里叶级数的系数,这样做的目的是使其表达形式与离散时间傅里,,Xk
~
0,k,N,1叶变换的形式相类似。如果将的主值周期记为,,由于以,,,,XkXk
~上两式中的求和范围均取为区间0~N-1,在次区间内=,因此可以得到: ,,,,xnxn
N,1~kn,,,,Xk,xnW0,k,N,1, ,Nn,0
1N,~1,kn,,,,xn,XkW0,n,N,1, ,NNn,0
表明时域N点有限长序列可以变换成频域N点有限长序列。显然,DFT,,,,xnXk
与DFS之间存在以下关系:
~
,,,,,, Xk,XkN
~
,,,,,, Xk,XkRkN
- 2 -
3、频率分辨率的问题
若信号最高频率为,按抽样定理,抽样频率应满足 fh
f,2fsh
也就是抽样间隔为T满足
11T,, f2fsh
一般取
,,f,2.5~3.0fsh
如果不满足的要求,就会产生频率响应的周期延拓分量相互重叠的现f,2fsh
象,也就是产生频率响应的混叠失真。
对于DFT来说,频率函数也要抽样,变成离散的序列,其抽样间隔为,这F0就是我们能得到得的频率分辨力,有它可引出时间函数的周期,也就是所取的
记录
混凝土 养护记录下载土方回填监理旁站记录免费下载集备记录下载集备记录下载集备记录下载
长度为 T0
1T, 0F0
从以上和两个公式来看,信号的最高频率分量与频率分辨力之间有TTfF0h0
,,11,,T着矛盾关系,要想增加,则时域抽样间隔就一定减小,,,而Tffhs,,f2fsh,,(抽样频率)就增加,由于抽样点数满足
fTs0,,N FT0
则此时增加,若是N固定的情况下,必然要增加,即分辨率下降。 fFs0
反之,要提高频率分辨力(减少),就要增加,当N给定时,必然导致TFT00
的增加(减小)。要不产生混叠失真,则必然会减小高频容量(信号的最高频fs
率分量)。 fh
要想兼顾高频容量与频率分辨力,即一个性能提高而另一个性能不变fFh0
(或也得以提高)的惟一办法就是增加记录长度的点数N,即要满足
- 3 -
f2fsh,, NFF00
这个公式是未采用任何特殊数据处理(例如加窗处理)的情况下,为实现基本DFT算法所必须满足的最低条件。如果加窗处理,相当于时域相乘,则频域卷积,必然加宽频谱分量,频率分辨力就可能变坏,为了保证频率分辨力不变,则须增加记录长度,也就是增加数据长度。 T0
4、MATLAB的图示说明:有效观察时间与补零后的DFT之间的关
系,以及与DTFT之间的关系
对8点正弦离散序列求8点、32点和64点DFT,观察频域变化(分别用绿、黄、红色表示)。结果:
矩形窗序列后补零的时、频域示意图
从图中可以看出:序列后补零可以降低栅栏效应;信号频谱的形状只取决于时域信号,与补零个数无关。补零并不能提高频谱分辨率,因为频谱分辨率只与时域数据的有效长度有关。
DTFT与DFT(或DFS)的关系:DFT时域序列为周期序列,周期为N;频域序列也是周期序列,周期也是N点。当N不断增大时,频域包络不变,但谱
N,,线变密;显然,时,时域序列变为非周期序列,频域为连续的频谱,即变化为DTFT。
5、教材
《信号与线性系统》,阎鸿森、王新凤、田惠生编,西安交通大学出版社 《数字信号处理教程》,程佩青编,清华大学出版社
(后附连续信号傅里叶变换的DFT近似计算)
- 4 -
傅里叶变换的DFT近似计算
连续时间非周期信号的傅里叶变换对为 ,,xt
,jt,,,,,,Xj,,xtedt (1) ,,,
,1jt, (2) ,,,,xt,Xj,ed,,,,2,
用DFT方法计算这一对变换的方法如下:
(1)将在t轴上等间隔(宽度为T)分段,每一段用一个矩形脉冲代替,脉冲的,,xt
,,,,,,xt,xnT,xn幅度为其起始点的抽样值,然后把所有矩形脉冲的面积相加。由于 t,nT
t,nT
dt,T ,,,,dt,n,1T,nT
,,dt,T ,,,,n,,,
,,,j,t,,,,Xj,,xtedt则得频谱密度的近似值为 ,,,
,jnT,,,,,,Xj,,xnT,e,T (3) ,,,
t,0(2)将序列 截断成从开始长度为的有限长序列,包含有N个抽t,,,,xn,xnT0样(即时域取N个样点),则上式成为
N,1,j,nT,,,,Xj,,TxnT,e (4) ,n,0
由于时域抽样,抽样频率为,则频域产生以为周期的周期延拓,如果频域f,1Tfss是限带信号,则可能不产生混叠,成为连续周期频谱序列,频域周期为(即时域的f,1Ts抽样频率)。
(3)为了数值计算,再频域上也要离散化(抽样)即在频域的一个周期()中也分fs成N段,即取N个样点,每个样点间的间隔为。频域抽样,那么频域的积分f,NFFs00式(2)式就变成求和式,而时域就得到原已截断的离散时间序列的周期延拓序列,其时域
周期为。这时 ,,k,T,1F000
,,d,,k,1,,k,,,000
- 5 -
N,1,d,,, ,0,,,k,0各参量的关系为
1NT,,,NT 0FF00又
,,2,F00则
FT,1,22,00T,,,,,,,,,,,,,,222 ,,,000ffTN,,ssss0
这样,经过上面三个步骤后,时域、频域都是离散周期的序列,推导如下:
第1,2两步:时域抽样、截断
N,1,j,nT,,,,Xj,,xnT,e,T (5) ,n,0
,1sjnT, (6) ,,,,xnT,Xj,,ed,,02,
第3步:频域抽样,得到
2,N,1N,1,jnk,jk,nT0N ,,,,,,,,,,Xjk,,TxnT,e,Txn,e,T,DFTxn,,0n,0n,0
2,N,1N,1jnk,jk,nT00N ,,,,,,xnT,Xjk,,e,FXjk,,e,,0002,k,0k,0
2,N,1jnk1N ,,,F,N,Xjk,,e,00Nk,0
2,N,1jnk1N ,,,f,Xjk,,e,s0Nk,0
,,,,,f,IDFTXjk,s0
,,,,,,Xjk,,Xj,,T,DFT,,xn (7) 0,,k,0
1 (8) ,,,,,,,,xnxtIDFTXjk,,,,t,nT0T
这就是从离散傅里叶变换法求连续非周期信号的傅里叶变换的抽样值的方法。
由,,及的上两个近似式求连续的及的方法,则可分别用频域抽样Xjk,,,,,xnXj,,,xt0
定理的插值公式和时域抽样定理的插值公式求得。
- 6 -