信号时频分析-讲义.

- -
从Fourier 分析到小波分析
1 Fourier 分析
所有客观存在的事物都包含着大量标志其本身所存的时间空间特征的数据，这就是该事物的信息。

当人们要了解事物某方面的情况时，通常要以各种手段把所需的信息表达出来，供人们观测和分析，这种对信息的表达形式称之为“信号”，所以信号是信息的载体。

信号是无处不在的。

如我们随时可听到的语音信号，随时可看到的视频图像信号，发电机组运行时的温度信号和振动信号等。

对一个给定的信号或过程，如)(t x ，我们可以用众多的方法来描述它，
如)(t x 的函数表达式，通过Fourier 变换所得到的)(t x 的频谱，即)(ˆωx
，再如)(t x 的相关函数，其能量谱或功率谱等。

在这些众多的描述方法中，有两个最基本的物理量，即时间和频率。

Fourier 变换和反Fourier 变换作为
桥梁建立了信号)(t x 与其频谱)(ˆωx
之间的一对一映射关系，从时域到频域的映射关系为Fourier 变换：
⎰∞
∞--=dt e t x x t j ωω)()( (1-1) 反过来，从频域到时域的映射关系为反Fourier 变换：
⎰∞
∞-=ωωπωd e x t x t j )(21
)( (1-2) Fourier 变换的本质思想是用一些简单的基本函数的加权和来近似和表示一个复杂的函数，这样的近似和表示有很多优点，它给我们分析和认识复杂现象提供了一种有效的途径，一些在时域内难以观察的现象和规律，在频域内往往能十分清楚地显示出来。

Fourier 变换和反Fourier 变换属于整体或全局变换，即只能从整体信号的时域表示得到其频谱，或者只能从整体信号的频域表示得到信号的时
域表示。

也就是说频谱)(ˆωx
的任一频点值都是由时间过程)(t x 在整个时域(-∞,∞)上的贡献所决定；反之，过程)(t x 在某一时刻的状态也是由其频谱)(ˆωx
在整个频域(-∞,∞)上的贡献所决定。

也就是说，)(t x 在任何时刻的微
小变化都会牵动整个频谱，而任何有限频段上的信息都不足确定任意小时间范围内的过程)
(t
x。

因此，Fourier变换建立的只是一个域到另一个域的桥梁，并没有把时域和频域组合在一起，所以频谱)
(ˆω
x只是显示了信号)
(t
x
中各频率分量的振幅和相位，而无法表现信号各频率分量随时间变换的关系。

t / s t / s
(a) x1(t) (b) x2(t)
图2.1 信号x1(t)和x2(t)
图2.1中的两个信号x1(t)、x2(t)可很好地说明Fourier变换的局限性，它们的时域表示如下：
)
18
sin(
)
12
sin(
)
6
sin(
)(
1
t
t
t
t
xπ
π
π+
+
=0≤t≤4s (1-3)
⎩
⎨
⎧
≤
≤
+
<
≤
+
=
s4
t
2
)
(18
sin
2
)
(12
sin
s2
t
)
(12
sin
)
(6
sin
2
)(
2t
t
t
t
t
x
π
π
π
π
(1-4) 这两个信号都是由三种频率分量组成，但它们的持续过程是不一样的，在x1(t)中，三种分量一直存在；而在x2(t)中，只有一个分量一直存在，另两个只是分别占信号整个过程的前一半和后一半。

f / Hz f / Hz
(a)2
1
|)
(ˆ|ω
x(b)2
2
|)
(
ˆ|ω
x
图2.2 信号x1(t)和x2(t)的频谱
图2.2是这两个信号的频谱2
1
|)
(ˆ|ω
x、2
2
|)
(
ˆ|ω
x，显然这两个不同的信号有相同的频谱，这说明Fourier分析不能将这两个信号区分开。

--
- -
2 短时Fourier 分析
2.1 基本定义
为了克服Fourier 变换不能同时进行时间——频率局域性分析的缺点，因发明全息照相技术而获诺贝尔奖的Gabor 于1946年提出了短时Fourier 变换(STFT)。

短时Fourier 变换的思想是把非平稳过程看成是一系列短时平稳信号的叠加，而短时性则是通过时间域加窗来实现，所以也称为加窗Fourier 变换，定义如下：
⎰∞
∞---=dt e t g t x t j x ωτωτ)()(),(ST FT (2-6) 式中)(t g 是分析窗函数，它在时域是紧支的，一般选用能量集中在低频处的实偶函数。

随着τ的不断变化，由g 所确定的窗口在时间轴上移动，使分析信号)(t x 逐步进入被分析的状态，因此该变换反映了信号)(t x 在时刻为τ、频率为ω的分量的相对含量。

Gabor 采用Gauss 函数)(t g a (式2-7)作为分析窗函数，因此用Gauss 函数作为窗函数的短时Fourier 变换也称Gabor 变换),(τωx G 。

Gauss 函数是
紧支的，它的Fourier 变换也是Gauss 函数)(ˆωa g
(式2.2-8)，从而保证了Gabor 变换在时域和频域都具有局域化功能。

a t a e a
t g 4/221)(-=π (2-7) 2
)(ˆωωa a e g -=
(2-8)
t / s rad/s (a))(t g a (b))(ˆωa g
图2.3 Gauss 函数)(t g a 及其Fourier 变换)(ˆωa g
可以证明，对于Gabor 变换，式2-9是成立的
⎰∞
∞-=)(ˆ),(ωττωx
d G x (2-9)
- -
这说明信号x (t )的Gabor 变换按窗口宽度精确地分解了x (t )的频谱)(ˆωx
，提取了它的局部频谱信息，当τ在整个时间轴上平移时，就给出了x (t )的完整的Fourier 变换，因此没有损失x (t )在频域上的任何信息。

短时Fourier 变换是能量守恒变换，对于任何窗函数下式都成立
⎰⎰⎰∞∞-∞
∞-∞∞-=τωωτπd d dt t x x 22|),(STFT |21|)(| (2-10) 在归一化条件下，即
dt t x ⎰∞
∞-2|)(|＝1，短时Fourier 变换是可逆的，其逆变换公式如下 ⎰⎰∞∞-∞∞--=
τωτωτπωd d e t g t x t j x )(),(STFT 21)( (2-11) 这里用短时Fourier 变换对上节中的信号x 1(t )和x 2(t )进行了分析，图2.4和图2.5给出了变换结果和相应的3D 显示效果。

t / s 图2.4 信号x 1(t )的Gabor 变换(窗口长度为1/8的信号长度)
t / s 图2.5 信号x 2(t )的Gabor 变换(窗口长度为1/8的信号长度)
通过图2.4和图2.5可以很容易地辨识出信号x 1(t )和x 2(t )，它们各个频率分量的持续时间也可轻易地知道，如对x 1(t )，它的三个频率成份就一直
f / H z f / H z
- -
存在，而x 2(t )中，只有一个频率成份一直存在，而其它两个频率成份只是占据了信号整个过程的前一半和后一半。

2.2.2 短时Fourier 变换的时频分辨率
短时Fourier 变换是一种时频变换，从上面的例子可以看出，它可以方便地分析非平稳信号，现在很自然会产生这么一个问题，是不是窗口越小越好呢？先看两个极端的例子。

当窗口函数选择为)(τδ时，这时
ωττωτj x e x -=)(),(S TFT (2-12)
信号的STFT 变成了信号x (t )，它保持了信号的所有时间特征，有完美的时域分辨率，却无任何频域分辨率。

另外，当取无限宽的窗函数时，即g (t )=1时，此时的短时Fourier 变换退化成一般的Fourier 变换，这时
)(ˆ),(S TFT ωωτx
x = (2-13) 信号的STFT 变成了信号)(ˆωx
，它有极好的频率分辨率，但没有任何时间分辨率。

为了分析Fourier 变换的时频局部化特性，引入相空间的概念。

所谓一个相空间是指以“时间”为横坐标，以“频率”为纵坐标的欧氏空间，而相空间中的有限区域被称为窗口。

相空间的作用是用来刻画一定的物理状态，因此它具有很强的工程背景。

从数学上来说，如果函数)()(2R L t g ∈，且)()(2R L t tg ∈，则)(t g 被称为窗口函数，相空间的点),(00ωt
⎪⎪⎩
⎪⎪⎨⎧==⎰⎰∞∞-∞∞-ωωωωωd g g dt t g t t g t 220220|)(ˆ|)(ˆ1 |)(|)(1 (2-14) 被称为窗函数)(t g 的中心。

式中，)(ˆωg
为窗函数)(t g 的Fourier 变换(下同)。

定义： ⎪⎪⎪⎩⎪⎪⎪⎨⎧⎪⎪⎭⎫ ⎝⎛-=∆⎪⎪⎭⎫ ⎝⎛-=∆⎰⎰∞∞-∞∞-2/122022/12202|)(ˆ|)()(ˆ1ˆ)()()(1ωωωωωd g g g dt t g t t t g g (2-15)
- -
为窗函数)(t g 的时宽和频宽。

相空间中以),(00ωt 为中心，以长为g ∆2，宽
为g
ˆ2∆的平行于坐标轴的矩形称为由g (t )所确定的时频窗口。

若g ∆越小，则说明)(t g 在时域上的局部化程度越高，当用这么一个窗进行短时Fourier
变换时，将取得较好的时域分辨率；同样，若g
ˆ∆越小，则说明)(t g 在频域上的局部化程度越高，用于短时Fourier 变换时，将取得较好的频域分辨率。

可以证明
g
g ˆ∆∆≥1/2 (2-16) 这就是Heisenberg 测不准原理在时频变换中的表现，它表明g ∆和g
ˆ∆之间存在一定的制约关系，两者不可能同时都任意小。

当且仅当g (t )取Gauss 函数时，式(2.2-16)中等号才成立。

从物理的直观意义上讲，信号的频率必须至少在一个周期内(t ∆≥0/1ω)进行测量，精确测量并认定某一时刻的频率是多少是没有意义的。

图2.6 短时Fourier 变换的相空间表示
图2.6给出了短时Fourier 变换的相空间表示。

很明显窗口函数g (t )一
旦选定，g ∆和g ˆ∆也随之确定。

因此，对于任意给定的0t 和0ω，短时Fourier
变换的时频分辨率可由尺度固定的分辨基元)]ˆ()[(0g
g t o ∆±⨯∆±ω来表示，也就是说，短时Fourier 变换在相空间中任何一点(00,ωx )给出的关于信号
x (t )的信息，都是由g ∆和g
ˆ∆这两个不确定量限定的。

由于短时Fourier 变换的时频窗口有相同的时宽和频宽，也就是窗口的大小形状是固定不变的，它在时域和频域的分辨率是固定不变的，即在高频段和低频段有同样的分辨率，这在图2.4和图2.5上有很好的表现，可看出信号中的三个不同的频率成份在图上表现出了同样的带宽。

为了更好
地说明g ∆和g
ˆ∆之间的相互制约性，这里用不同宽度的窗函数对前面的信号x 1和信号x 2进行了分析，如图2.7所示。

- -
t / s
t / s 图2.7 信号x 1(t )和x 2(t )的Gabor 变换(窗口长度为1/4的信号长度) 和图2.4和图2.5相比，那里的窗口宽度为1/8的信号长度，而现在的窗口用的是1/4的信号长度。

可明显看出，图2.7中的频率分辨率提高了，每个频率分量都在上面表现为比图2.4和图2.5中更窄的带。

而它在时域的分辨精度下降了，这通过信号x 2(t )的分析可看出。

在图2.5中，两个各占信号过程一半的分量基本上以2s 点为分界，而在图2.7中，这两个分量在时间轴上却出现了交叠。

通过上面的分析可知，当用短时Fourier 变换分析信号时，如果想对高频分量分析取得很好的时域分辨率，就必须选择宽带的短时窗；如果想对低频分量分析取得很好的频域分辨率，就必须选择窄带的宽时窗，但无法同时达到这两个目的。

而实际的信号过程是很复杂的，无论是单一的还是多分量的信号，为了提取高频分量或速变成份的信息，时域窗口g ∆应尽量
窄，而同时容许频域窗g
ˆ∆适当放宽，因为更高频率分量即使有较大的绝对频率误差，仍可以使相应的相对误差保持不变；对于慢变信号或低频成份，
频域窗口g
ˆ∆就应当尽量缩小，保证有较高的频率分辨率，以保证频率的相对误差满足提取信息的基本需要。

简而言之，实际的信号的分析需要时频
窗口具有自适应性，它可按上面的情形自动改变g ∆和g
ˆ∆的大小，高频时频窗宽，时窗窄；低频时则频窗窄，时窗宽。

这么一个自适应窗口的相空间特性可用图2.8表示。

f / H z
- -
图2.8 自适应窗的相空间表示
短时Fourier 变换还有一个缺点就是它的离散形式没有正交展开，难于实现高效算法，这也大大限制了其应用范围。

3 连续小波变换
从Fourier 变换到短时Fourier 变换再到对自适应窗口的需求反映了信号分析处理过程中一个共同的基本要求，这就是具有自适应窗口特性和平移功能，为了实现高效算法，要求对信号x (t )进行变换处理的积分核应具有正交基。

归结起来，变窗口、平移和正交性是作为信号分析最有效的数学工具的主要条件。

小波变换(Wavelet Transform)正是为了满足这个需求而发展起来的。

3.1基本定义
设x (t )是一有限能量函数，即)()(2R L t x ∈，则该函数的小波变换定义为以函数族)(,t b a ψ为积分核的积分变换，如下式所示
);,(ψb a W x ,()()a b x t t dt ψ∞
-∞=⎰ a >0 (3-17) 函数族)(,t b a ψ由基本小波函数)(t ψ通过伸缩和平移产生，如下所示
)()(2/1,a
b t a t b a -=-ψψ (3-18) 式中a 是尺度参数，b 是定位参数，2/1-a 因子是归一化常数，用来保证变
- -
换的能量守恒，即
⎰⎰∞
∞
-∞∞-==dt t dt t t b a b a 22,2,|)(||)(|||)(||ψψψ (3-19) 小波函数的频域表示如下 )(ˆ)(ˆ,ωψωψ
ωa e a j b a -= (3-20) 可以看出，当a 减小时，小波函数的时宽减小，频宽增大，且)(,t b a ψ的窗口中心向|ω|增大方向移动；当a 增大时，小波函数的频宽减小，时宽增大，且)(,t b a ψ的窗口中心向|ω|减小方向移动。

这说明连续小波变换的局部化是变化的，在高频处时域分辨率高，频域分辨率低；而在低频处时域分辨率低，频域分辨率高，即具有“变焦”的性质，这也正是我们追求的自适应窗的性质。

下面以常用的Morlet 小波函数来说明。

rad/s rad/s rad/s
a = 0.5 a = 1 a = 2 图2.9 不同尺度下的Morlet 小波函数及其频谱(50=ω)
Morlet 小波是最常用的复值小波，由下式给出
2/2/4/12200)()(t t j e e e t -----=ωωπ
ψ (3-21)
其Fourier 变换为 ][)(ˆ2/2/2/)(4/122020ωωωωπωψ
------=e e
e (3-22) 当0ω≥5时，020≈-ωe ，则Morlet 小波函数可简化如下 2/4/120)(t
t j e e t ---=ωπψ (3-22)
- - 其Fourier 变换相应地变为
2/)(4/120)(ˆωωπωψ
---=e (3-23) 图2.9 给出了不同尺度参数a 下的Morlet 小波函数及其频谱，可以看出，当a 增大时，小波函数被延展，对应的Fourier 变换缩小，频域窗口中心减小，反之亦然。

另外，可以看出Morlet 小波函数的频谱具有带通的特征，事实上，所有满足下面条件的小波函数的频谱都具有这样的性质。

小波)(t ψ的选择不是唯一的，很多函数可用来作为小波基函数，但也不是任意的，它的选择应该满足以下条件：
1. 定义域是紧支撑的，即在一个很小的区间之外，函数值为零，它保证了
函数的速降特性，以便获得时域局域化。

2. 平均值为零，也就是⎰∞
∞-=0)(dt t ψ，甚至)(t ψ的高阶矩也应该为零，即
⎰∞∞-=0)(dt t t k
ψ k = 0,1,…N -1 (3-24) 我们称满足这项要求的小波函数具有k 阶消失矩，它可以消除信号x (t )的多项式展开中k t (k <N )各项在小波变换中的贡献，以便突出信号的高阶起伏和高阶导数中可能存在的奇点，因此小波变换能突出信号的高阶变化。

均值为零的条件也称小波容许条件，即
∞<=⎰∞
∞-ωψωωψd C |||)(ˆ|2 (3-25) 式⎰∞
∞--=dt e t t j ωψωψ)()(ˆ，这个条件使函数)(t ψ的波形必定具有振荡性，并且随着k 的增大，)(t ψ的振荡性会越来越强。

从小波变换的定义可看出，小波变换是一线性变换，它的物理图案就是用一族频率不同的振荡函数作为窗口函数)(,t b a ψ对信号x (t )进行扫描和平移，其中a 为改变振荡频率的伸缩参数，b 为平移参数。

这时小波变换在某种意义上类似于短时Fourier 变换，但不同的是，小波变换的时域和频域分辨率与频率有关。

在高频段，小波变换能达到高时域分辨率，而频域分辨率比较差，对低频段则刚好相反。

而短时Fourier 变换在所有频段的时域和频域分辨率都是不变的。

对于所有的x (t )，)()(2R L t ∈ψ，x (t )的连续小波变换的逆变换可由下式
- -
给出
⎰⎰∞∞-∞∞--=dadb t b a W a C t x b a x )();,(1)(,2ψψψ (3-26)
利用Parseval 公式和Fourier 变换的相似性，很容易证明上述逆变换公式。

这说明信号x (t )的小波变换并没有损失任何信息，变换是守恒的，因而下式成立
⎰⎰⎰∞∞-∞
∞--∞∞-=db b a W da a C dt t x x 222|);,(|1|)(|ψψ (3-27) 2.3 小波变换的分辨率
我们同样可以通过相空间来分析小波变换的分辨率，定义相空间的点),(0ˆ0,,b a b
a t ψψω为小波函数)(,t
b a ψ的中心。

⎪⎪⎪⎩
⎪⎪⎪⎨⎧==⎰⎰⎰⎰∞∞∞∞-∞∞
-02,02,02,2,0|)(ˆ||)(ˆ||)(||)(|,,ωωψωωψωωψψψψd d dt t dt t t t b a b a b a b a b a b a (3-28) 定义
⎪⎪⎩
⎪⎪⎨⎧⎪⎭⎫ ⎝⎛-=∆⎪⎭⎫ ⎝⎛-=∆⎰⎰∞∞-∞∞-2/12,20,2/12,20,|)(ˆ|)(ˆ|)(|)(,,ωωψωωψψψψψd dt t t t b a b a b a b a b a b a (3-29) 为小波函数)(,t b a ψ的时宽和频宽。

不难证明下面各式成立
b at t b a +=000
,1,ψψ (3-30) 000,1,1ψψωωa b a =
(3-31) 0,1,ψψ∆=∆a b a (3-32)
0,1,ˆ1ˆψψ
∆=∆a b a (3-33) 由(3-32)和(3-33)还可推出
=∆∆=∆∆0,10,1,,ˆˆψψψ
ψb a b a 常数 (3-34)
- -
对小波函数)(,t b a ψ而言，它的频谱)(ˆ,ωψb a 具有带通特性，而0ˆ,b a ψω就是它
的通频带中心，带宽BW 则为2b a ,ˆψ
∆。

由上述公式可看出，随着a 的增大，0ˆ,b
a ψω减小，这表明带通的中心向低频分量偏移，这时小波变换分析的是信号的低频分量，而此时
b a ,ˆψ
∆相应减小，b a ,ψ∆相应增大，因此在低频段，小波变换可达到较高的频率分辨率和较低的时域分辨率，反之亦然。

这正是我们所需要的图2.8所示的自适应窗的特点。

另外，由式(3-34)可知，在小波函数的相空间中，即使每个分辨基元的宽度(确定了时域分辨率)和高度(确定了频域分辨率)在各处不一样，但它的面积是一常数，这也正是Heisenberg 测不准原理在小波变换中的体现。

再看看反映小波函数滤波特性的品质因子Q ，由(2.3-31)和(2.3-33)还可推出
常数＝带宽中心频率=∆=∆=0,10,0ˆ2ˆ20,1,ψωψ
ωψψb a b a Q (3-35) Q 为常数说明小波变换相当于一个恒Q 的带通滤波器，可用图2.10来表示小波变换的滤波特性。

与此相对应，短时Fourier 变换的滤波特性可用2.11来表示。

ω0 ω2 ω4 ω8 ω 图2.10 小波变换的滤波特性
ω0 ω1 ω2 ω3 ω4 ω5 ω6 ω7 ω8 ω 图2.11 短时Fourier 变换的滤波特性
图2.12则给出了小波变换的相空间表示。

- -
图2.12 小波变换的相空间表示
这里用小波变换对上节中的信号x 1(t )和x 2(t )进行了分析，如图2.13所示的变换结果。

t / s t / s 图2.13 信号x 1(t )和x 2(t )的小波变换(Morlet 小波 60=ω)
为了对比，这里进行了尺度和中心频率之间的转换(以后的所有有关小波变换的图如未经标注，都默认经过此转换)，对于任一小波，它的频谱都具有带通性质，因此这样的变换是可行的，变换公式为
a b a /0ˆ0ˆ0,1,ψψωω= (3-35)
可以看出，小波变换和短时Fourier 变换的结果非常相似，但不同的是，小波变换对高频分量的频域分辨率没有在低频的高，很明显频率越高的分量在小波变换图上对应着越宽的频带。

另外，从信号x 2(t )的小波变换也可看出，小波变换对低频分量的时域分辨率没有对高频分量的高，很明显，本应只是占据0～2s 的那个分量在小波变换图上直到近3s 处还出现，而占据2～4s 的分量在小波变换图上只是延伸到了将近1.8s 处。

值得注意的是小波变换存在边界扭曲现象，上面的小波变换虽然采取
- -
了一些措施（对信号两端进行了对称拓展），但图中还是可以看出明显的不光滑的边缘效应。

之所以发生边界扭曲，是因为进行小波变换需要作卷积，而通常对有限长度的信号进行卷积时，都会产生边界扭曲。

这种边界扭曲带来的误差是与实际的信号有关的。

一般信号处理都是对实际信号加矩形窗进行截断，所以在卷积时矩形窗外的部分没有值(一般取零值卷积)，这与实际情况不符合，因此就带来了误差。

为了解决边界扭曲问题，在进行小波变换之前，需要将信号扩展。

对于不同的信号，可能会用到不同的方法来解决边界问题。

2.3.3 确定最大和最小尺度参数
用连续小波变换分析信号x (t )时，取小的尺度参数时分析的是高频率成份，但为了分析更高的频率，尺度参数a 的取值是否要取得非常小呢？相应地，是否可用非常大的尺度参数来分析信号中的超低频率成份呢？事实上尺度参数的选择范围不是任意的，为了节约计算量，也为了消除小波变换结果中的无用信息，尺度参数的选择应该遵守一些原则。

以Morlet 小波为例，它的频谱具有带通特性，当尺度参数为a 时，它的中心频率为a a /0ωω=，通过公式(2.3-23)不难得出它的半功率带宽为
a
BW a 2ln 2= (3-36) 因此，在尺度a 下，小波函数频带的上下截止频率分别为
a
a H 2ln /0+=ωω (3-37) a
a L 2ln /0-=ωω (3-38) 如果信号的采样间隔是T s ，则它的采用频率为f s =1/ T s ，根据Shannon 采样定理，采样频率f s 必须高于信号x (t )的最高频率的两倍，因此
s
s T f ππ22=≥)2ln /(220a a H +=ωω (3-39) 即
a ≥π
ωs T )2ln (0
+ (3-40) 这就确定了尺度参数的最小取值m i n a ，计算比min a 还小的尺度参数下的小波
- -
变换将没有任何意义。

在实际计算中也可以根据信号的截止频率c f (2c f ≤f s )来确定最小尺度参数m i n a ，这时
c f a πω22ln 0min +≈
(3-41) 最大的尺度参数可以通过小波函数)(,t b a ψ的时宽来确定，为了保证小波变换具有一定的时域分辨率，考虑到|)(,t b a ψ|的幅值在3b a ,ψ∆之处将降到本身的99.9%，因此一般要求
3b a ,ψ∆≤信号长度/2
对Morlet 小波来说，就是
a ≤信号长度/6 (3-42)
这就确定了尺度参数的最大取值m a x a 。

从上面的分析可知，信号的采样频率或截止频率决定了最小尺度参数m i n a ，即最大分析频率；而信号的采样长度决定了最大尺度参数m a x a ，即最小分析频率；采样频率和采样长度不会影响小波变换的分辨率，小波变换的频域分辨率和时域分辨率只和所选择的小波函数及尺度参数有关。

4 各种变换之间的比较[12]
下面用表格的形式给出Fourier 变换、短时Fourier 变换和小波变换之间的比较。

表2.1 Fourier 变换
变换类型
频率 分析函数
三角函数（是时域支撑区无限的振荡函数） 变量
频率 提供的信息
信号包含的频率成份 适用场合
平稳信号（其频率成份不随时间变化） 备注 FFT 的计算量为n n log
表2.2 短时Fourier变换
变换类型时间－频率
分析函数有限长度的窗函数和三角振荡函数的复合函数(每次分析时，窗函数的长度是固定的，但窗函数所包络的三角函数的频率是变化的)
变量频率；窗口位置
提供的信息窗函数越小，分析得到的时域信息越好，但可能损失在低频的频域信息；反之，窗函数越大，可以得到很好的频域信息，但时域
的分析精度将降低
适用场合准平稳信号（在窗口尺度下信号是平稳的）备注也称加窗Fourier变换，当取Gauss函数为窗函数时，也称Gabor 变换。

Fourier变换存在正交基，但短时Fourier变换不存在正交基，
难于实现高速算法
表2.3 小波变换
变换类型时间－尺度
分析函数具有固定振荡次数，时域有限支撑的波；小波函数通过伸缩来改变“窗口”大小，同时也改变尺度，因此可以在不同的尺度下观测信
号；由于振荡次数是不变的，因此小波函数的频率只是随着尺度的
变换而变化
变量频率；小波位置
提供的信息大尺度对应窄的小波函数，它能提供好的时域局部化信息，但频域局部化信息较差；小尺度对应宽的小波函数，它能提供好的频
域局部化信息，但时域局部化信息较差
适用场合非平稳信号，如一些简单的信号(如线性调频信号)和一些在不同尺度下具有不同频率的信号(如分形信号)
备注小波变换具有连续和离散两种形式；正交小波变换和双正交小波变换是离散小波变换的两种特殊形式；对于正交小波变换，当用
快速小波变换时，对长度为n的信号，其计算量是2kn，其中k
是小波支集的长度
--
5 小结
1)Fourier变换是一个全局变换，缺乏局域分析性能，不能提供任何时域
信息，只适合于平稳信号分析。

2)短时Fourier变换通过加窗，使之可以在时域和频域内同时对信号进行
分析，因此可以分析非平稳信号；但其窗口大小是固定的，缺乏自适应性，因此它在任何频段的频域分辨率和时域分辨率是固定的；它的离散形式缺乏正交基，难于实现高速算法。

3)小波变换通过伸缩和平移对信号进行分析，从而可以实现多分辨率分
析；它的窗口具有自适应性，在高频段可以获得好的时域分辨率和较低的频域分辨率，而在低频段则刚好相反；它的离散形式具有正交基，可以实现快速计算。

4)小波变换存在边界扭曲现象，可以通过边界拓展来减轻其影响。

5)小波变换的最小尺度参数由采样频率决定，最大尺度参数则可由信号的
采样长度决定。

--。