小波变换第1讲基本概念
64
1999年1月
继电器 第27卷 第1期 RE LATY
专家讲座
小波变换
第1讲 基本概念
董新洲 贺家李 天津大学电力系 (300072)
葛耀中 西安交通大学电力系 (710049)
【摘要】 从信号处理及工程应用的角度出发, 介绍了小波及小波变换的常用数学术语、换的时频局部化特性, 说明了小波变换是分析非平稳变化信号的理想工具, 程应用奠定了基础。
【关键词】 小波 小波变换 时频局部化
1 引言
, 一种表示, 、变换后的形式下变得简洁和特征明确。
信号有两类:一类是稳定变化的信号; 一类是具有突变性质的、非稳定变化的信号。
对于稳定变化的信号, 工程上最常使用的一种变换就是F ourier 变换。F ourier 变换把一个周期变化的信号表示成一族具有不同频率的正弦波的线形叠加。从数学上讲,F ourier 变换是通过一个被称为基函数的函数w (x ) =e ix 的整数膨胀而生成任意一个周期平方可积函数f (x ) ∈L 2(0, 2π) , 其中2
) 称为平方可积函数空间。通过F ourier 变L (0, 2π
换, 在时域中连续变化的信号转化为频域中的信号, 因此F ourier 变换是一种纯频域分析方法。
对于具有突变性质的、非稳定变化的信号, 人们不只感兴趣该信号的频率, 而且尤其关心该信号在不同时刻的频率, 换句话说, 需要时间和频率两个指标来刻划信号。显然, F ourier 变换是无能为力的。这是因为:Fourier 变换在频域上是完全局部化了的(能把信号分解到每个频率细节) , 但在时域上却没有任何分辨能力。因此需要时频分析方法来分析这种信号。
时频分析方法的典型例子是窗口F ourier 变换。一个具有有限能量的模拟信号f (t ) 的窗口Fourier 变换被定义为
) =G (ω, τ
) , 则上式就是Fourier 1g (t -τ
。因此窗口Fourier 变换可描述为:对于待分析的信号f (t ) 先开窗再做Fourier 变换, 随着窗的移动, f (t ) 被一部分一部分地分解。其中的时限函数g (t ) 因此被称为窗函数[1]。
) 既是频率ω的函数, 又由式(1) 可见, G (ω, τ
是时间τ的函数, 因此窗口Fourier 变换提供了信号f (t ) 在时间τ的频率信息, 它是一种时频分析方法。
实际信号是由多种频率分量组成的, 当信号尖锐变化时, 需要有一个短的时间窗为其提供更多的频率信息; 当信号变化平缓时, 需要一个长的时间窗用于描述信号的整体行为。换句话说, 希望能有一个灵活可变的时间窗。而窗口Fourier 变换无法做到这一点。这是因为窗口Fourier 变换的窗函数g (t ) 的大小和形状是固定不变的, 不能适应不同频率分量信号的变化。这就导致了小波变换的出现。
2 连续小波变换
受Fourier 变换和窗口Fourier 变换的启发, 可以寻找另外一个单一函数的膨胀和平移来表示信号f (t ) , 这样的函数被称为基小波Ψ(t ) , 它必须满足容许性条件
C Ψ=
∫
∞
-∞
|ω|
-1
2
|Ψ^(ω) |d ω
(2)
或者等价地
Ψ(t ) d t =0
∫
R
(3)
∫
R
) e 2i ωt d t f (t ) g (t -τ
(1)
由基小波的伸缩和平移所生成的函数族
Ψa , b (t ) 被称为连续小波
Ψa , b (t ) =|a |
-
其中, g (t ) 是具有紧支集的时限函数。显然, 窗口Fourier 变换和Fourier 变换的区别就是前者多了一个时限函数g (t ) 。
2Ψ(
) , a
(4)
a , b ∈R , a ≠0
其中:a 称为尺度因子, b 称为平移因子。
小波变换65
信号f (t ) ∈L 2(R ) (L 2(R ) 又被称为能量有限信号空间[2]) 关于小波Ψa , b (t ) 的连续小波变换被定义为
(W Ψf ) (a , b ) =|a |
-
a
Ψ^,
3
a
+
a
Ψ^)
3
a
(8)
) , 窗的高度(频窗) 和宽度
2
∞
其中心在(b +at 3,
(时窗) 分别为2
-∞
f (t ) Ψ(
) d t a
(5)
或者写成内积形式[2]
(W Ψf ) (a , b ) =
(6)
a
ΔΨΨ。^, 2a
信号f (t ) 可以由它的小波变换重构, 重构公式为
f (t ) =
窗函数决定的窗口是对信号f (t ) 局部性的一
次刻化, 小波窗函数提供了信号f (t ) 在时段(b +
333
ΔΔat -a +a -Ψ, b +at Ψ) 和频带(Ψ^, ω
3
a a
C Ψ
∫
-
∞∞
-∞-∞
(W Ψf ) (a , b ) Ψa , b
a
2d b
a
(7)
a
, 小波变换具有时频) a
外, , 同样的函数|a |
2Ψ(
个对偶。
, 可以看出它和Fourier 变换的异同:
①小波变换和Fourier 变换都是使用一个被称为基函数的单一函数Ψ(t ) 和e it 来表示原信号;
②小波变换是用尺度因子a 对基小波Ψ(t ) 进行伸缩, Fourier 变换是用膨胀因子ω对基函数e it 进行伸缩;
③小波变换还用平移参数b 对基小波进行平移, 而Fourier 变换只有伸缩没有平移, 因此小波变换具有时频局部化性能, 而Fourier 变换是一种纯粹的频域分析法;
作为时频分析方法, 小波变换和窗口Fourier 变换都能给出信号在某一时刻的频率信息, 但两者有本质的差别, 可以通过下述看出。
() 知, 小波窗函数的窗口形状是变。对于高频信号, 时窗变窄, 频窗变宽, 有利于描述信号的细节; 对于低频信号, 时窗变宽, 频窗变窄, 有利于描述信号的整体行为。正是由于小波函数的这种变窗特性, 使它能够表示各种不同频率分量的信号, 特别是具有突变性质的信号。
窗口Fourier 变换不同于小波变换。
设窗口Fourier 变换的时限函数g (t ) 和它的Fourier 变换^g (ω) 都是窗函数, 其中心与半径分别为t 3, ω3, Δg , Δ^若令g 。
ωi t
) w ω, τ=e g (t -τ
则w ω, τ和它的Fourier 变换w ^ω, τ也是窗函数, 它
们也定义了一个时频窗
(t 3+τ-Δg , t 3+τ+Δg ) ×(ω3+ω-Δ^g , ω3+ω+Δ^(9) g )
) 和频率范由上式可见, 除了时间上的移动(τ
围(ω) 的变化外, 窗口的大小和形状是不变的。因此, 窗口Fourier 变换不能适应不同频率信号的变化。但在实际中, 为了检测高频信号, 必须选择足够窄的时间窗, 而在检测低频信号时, 必须选择足够宽的时间窗, 这个矛盾在窗口Fourier 变换中是无法解决的。而这正是小波变换的优点。
3 小波变换的时频局部化性能
为了说明小波变换的时频局部化性能, 首先给
出窗函数的定义。
定义1 非平凡函数w (t ) 被称为是一个窗函数:如果w (t ) ∈L 2(R ) , 且tw (t ) ∈L 2(R ) 。表征窗函数的两个参数是窗函数的中心与半径[3].
设基小波Ψ(t ) 及其Fourier 变换Ψ^(w ) 都是窗函数, 其中心与半径分别为t 3, ω3, ΔΨ, ΔΨ^, 则小波函数Ψa , b 和它的Fourier 变换Ψ^a , b (ω) 也是窗函数, 它们一起在时间—频率平面上定义了一个矩形窗(时频窗)
(b +at
3
4 两类重要的小波变换
连续小波变换(W Ψf ) (a , b ) 是信号f (t ) 的一种表示, 在这里, 参数a , b 取遍整个实轴。若对参数a , b 的取值作一些限制, 则有不同类的小波变换[3]。常用的小波变换有两类:离散小波变换[3]和二进小波变换[4]。
(1) 离散小波变换
取a =j , b =j ; j , k ∈Z 。即尺度参数a 使用
22
Δ-a Ψ, b +at
3
Δ+a Ψ) ×(
3
a
-
66
1999年1月
继电器 第27卷 第1期 RE LATY
其中, Ψ 是Ψ的对偶, 它是前述共轭概念的推
广。
(16) , 信号f (t ) 可以表示为不同频 根据(15) 、
率的小波分量的和, 即
f (t ) =∑g j =…+g -1(t ) +g 0(t ) +
j
2的幂把频率轴剖分为二进的、相互毗邻的频带, 同
时, 平移参数b 只在时间轴上的二进位值取值。此
时, 式(5) 的连续小波变换转换为离散的小波变换(W Ψf ) (j , j ) =
22
∫
∞
-∞
f (t ) [2
j/2
Ψ(2j t -k ) ]d t (10)
(11)
Ψj , k 就是小波函数, 它被写成Ψj , k =2j/2Ψ(2j t -k ) , j , k ∈Z
(2) 二进小波变换
g 1(t ) +…(17)
取a =
, j ∈Z ; b ∈R 。即只对尺度参数a 进2j
行二进离散, 而平移参数b 保持连续变化。此时, 式(5) 的连续小波变换转换为半离散的小波变换或者称为二进小波变换
(W Ψf ) (j , b ) =
2
∫
∞
-∞
f (t ) {2
j/2
j t -b (12)
小波函数Ψj , b 被写成
Ψj , b =2j/2Ψ[2j (t -b ) ], j ∈Z , b ∈R
(13)
由于二进小波变换具有一个重要的特性—平移不变性[4], 因此被广泛应用于模式识别和信号的奇异性检测中。
小波分析中还有一类变换, 就是所谓的小波包变换[5]。前述的小波变换把信号按照二进频带分解成小波分量, 小波包变换则是对于小波分量的再分解。以下仅给出正交小波包的定义。
定义 函数Ψn , n =2l , 2l +1; l =0, 1, …, 为正交小波包, 若
Ψ2l (x ) =
P Ψ0(2x -∑
k
k
当信号f (t ) 被分解为小波后, 对信号的研究就转化为对其小波分量或在某一尺度(不同的j ) 下的小波变换的研究。
(2) 、不同的研究目的、采用不同。譬如, 要求小波函数具有正交性、一定的对称性和光滑性等, 这些要求经常矛盾, 需要在应用中合理予以取舍。在两类重要的小波变换中, 因为只有部分连续小波变换的值用于重构原信号, 因而对小波函数提出了更高的要求, 与之对应的小波有两类:R 小波和二进小波。
定义2 一个R 函数[3]Ψ被称为是一个R 小波, 如果Ψ的对偶Ψ 存在。
定义3 一个函数Ψ被称为是一个二进小波[4], 如果存在常数0
∞
A ≤
-∞
∑|
j
Ψ) |≤B (2-ω
(18)
小波和小波变换虽然种类繁多, 但本文将局限
于介绍离散小波更换和二进小波变换、R 小波和二进小波。
k )
(14)
k )
3
6 小波分析的发展史及应用概况
如前所述, 小波变换的本质是用一个小波基函数来表示一个能量有限的信号。小波分析的思想来源于伸缩与平移。
小波分析方法的提出可以追溯到1910年Harr 提出的第一个小波规范正交基;1984年法国地质学家M orlet 和理论物理学家G rossman 提出了连续小波变换的概念、1986年法国数学家Meyer 创造性地构造出了具有一定衰减性的光滑化函数—正交小波函数, 标志着小波热潮的开始。
1987年, 法国人Mallat [6]提出了多分辨分析的概念, 为统一地构造小波函数奠定了基础, 同时给出了以他的名字命名的小波分解与重构算法—Mallat 算法。
Ψ2l +1(x ) =
q Ψ0(2x -∑
k
k
而系数q k =(-1) k P -
k +1
。
5 信号的小波表示
(1) 信号表示为小波分量的叠加
研究小波变换的目的在于用小波表示信号。对于离散小波变换和二进小波变换, 这种表示可由它们的逆变换直观看出
∞
f (t ) =f (t ) =
j , k =-∞
∑(W Ψf ) (2
[2∑∫∞
∞
j
,
2j
) Ψ j , k (15)
∞
j/2
j , k =-∞-
(W Ψf ) (2-j , b ) ]×
(16)
-j
{2Ψ [2j (t -b ) ]}d b
小波变换67
学出版社,19941
1988年,Daubechies [7]构造了具有有限支集的正
交小波基, 至此小波分析的系统理论初步得到建立。
1990年, 崔锦泰和王建中构造了基于样条的半正交小波函数[8], 使得小波分析的系统理论得到完善。
小波分析虽然是一个数学领域, 但自始至终与应用科学交织在一起。1987年Mallat 首先把小波变换应用于计算机图象的分解与重构[6], 之后, 他和H wang Wenliang 、Zhong Sifen 等人系统地建立起了一套信号奇异性检测的理论和方法[4]; 同时, 国内外众多学者也开始把小波变换应用于自己的研究领域, 像分形[2]、医学成象与诊断、地震勘探与处理、机械振动与噪声[9]、行波信号分析[10]、, 得了大批成果。
2 龚怀云, 寿纪麟, 王绵森1应用泛函分析1西安:西安交
通大学出版社,19851
3 崔锦泰[美]著, 程正兴译1小波分析导论1西安:西安交
通大学出版社,19961
4 Mallat ,S. and H wang ,W. L. S ingularity detection and processing with wavelets. IEEE T rans. on In formation Theory ,1992,38. 5 程正兴1小波分析算法与应用1西安:西安交通大学出
版社,19981
6 Mallat ,S. A theory of multi -res signal decom position :thewavelet Pattern Anal. Machine In 2117, I. of com pactly supported and Appl. Math. 41,1988.
C K and Wang J Z 1A cardinal spline approach to wavelets. Proc. Amer. Math. S oc. ,1991.
9 赵纪元, 何正嘉等1小波包—自回归谱分析及在振动诊
断中的应用1振动工程学报,1995, (8) 1
10 董新洲, 耿中行, 葛耀中, 张伏生, 徐丙垠1小波变换应
7 结论
(1) 小波变换是在F ourier 变换基础上发展起来
用于电力系统故障信号分析初探1中国电机工程学报,
1997, (6) .
收稿日期:1998—10—22
董新洲 副教授, 博士后, 天津大学电力系, 主要从事行波保护、行波测距以及小波变换在电力系统中的应用等研究工作。
贺家李 教授, 博士生导师, 长期从事电力系统故障分析、继电保护领域的科学研究工作。
葛耀中 教授, 博士生导师, 长期从事电力系统继电保护基础理论与应用的研究。
的一个新的数学分支;
(2) 小波变换的思想来自于伸缩与平移; (3) 小波变换与小波函数具有多样性;
(4) 小波变换具有时频局部化性能, 是分析具有突变性质的、非平稳变化信号的理想工具;
(5) 小波变换在电力系统、特别是在故障分析与继电保护中有着广阔的应用前景。
参考文献
1 秦前清, 杨宗凯1实用小波分析1西安:西安电子科技大
WAVE LET TRANSFOR M
P art One B asic Concept
Dong XinZhou , He Jia Li (P ower Department of T ianjin University ,T ian Jin ,300072,China )
G e Yaozhong (X i ’an Jiaotong University ,710049,China )
Abstract As viewed from signal processing and engineering application ,the often -used mathematical terms and basical concept of wavelet and wavelet trans form and the time -frequence localization feature of the wavelet trans former are introduced in this paper. I t shows the wavelet trans form is a desired tool to analyze non -steady variable signal and gives a basis to further understand the wavelet trans form and to apply it in engineering.
K eyw ords Wavelet Wavelet trans form T ime -frequence localization
(上接63页)
DEVE LOPMENT OF CABINET STRUCTURE AN D OUTG OING LINE T YPE FOR PR OTECTIVE RE LAY
Zhang Jianjie (Xuchang Relay Research Institute ,Xuchang ,461000,China )
Abstract The development of cabinet structure and outg oing line type for power system is systematically analyzed and com pared. A com plete scheme is presented and has been applied in the new generation of hardware of X J E lectric C orporation. K eyw ords Cabinet structure Outg oing line type