没有合适的资源?快使用搜索试试~ 我知道了~
首页离散沃尔什变换及其快速算法
资源详情
资源评论
资源推荐

离散正交变换导论
第十一章 离散沃尔什变换及其快速算法
第一节 沃尔什级数表示
全体沃尔什函数 ",2,1,0 )},,({ =kkWal
θ
以及它的子集 都形
成群,并且是可交换群。这说明它们满足封闭性的条件。它们还满足完备性条件。
12,,2,1,0 )},,({ −=
p
kkWal "
θ
一个绝对可积信号 x(t),即 ∞<
∫
∞
0
)( dttx 可以展开成完备的正交函数系的级数:
∑
∞
=
=
0
),()()(
n
tnnatx
ϕ
∫
=
b
a
t
t
dttntxna ),()()(
ϕ
即广义傅里叶级数,其中
{}
),( tn
ϕ
为任意正交函数系。
因为沃尔什函数也是完备的正交函数系,那么满足条件的连续函数 x(t)也能展开成沃尔什
级数。为了方便,我们把 x(t)定义在归一化的半开区间 上,以
)1,0[ )(
θ
x
表示,则
∑
∞
=
=
0
),()(
k
Wk
kWaldx
θθ
(11.1)
",2,1,0 ),()(
1
0
==
∫
kdkWalxd
Wk
θθθ
(11.2)
使用沃尔什排列的 ),(
θ
kWal
W
,有下列关系
()
(
)
() ()
⎭
⎬
⎫
−==
==
12 ,,
2 ,,
ikiSalkWal
ikiCalkWal
W
W
θθ
θθ
其中 i 为 ),(
θ
kWal
W
的归一化列率
⎪
⎩
⎪
⎨
⎧
+
=
=
为奇数
为偶数
kk
kk
k
i
2)1(
2
0 0
所以(11.1)式变成
∑
∞
=
++=
1
0
)],(),([),0()(
i
ii
iSalbiCalaCalax
θθθθ
(11.3)
式中
⎪
⎩
⎪
⎨
⎧
−==
==
==
12
2
0
00
ikdb
ikda
kda
ki
ki
有限项级数展开式:
- 1 -

离散正交变换导论
∑
−
=
=
1
0
),()(
N
i
Wk
kWaldx
θθ
(11.4)
),2()],(),([),0()(
2
12
1
0
θθθθθ
NSalbiSalbiCalaCaldx
N
N
i
ii
+++=
∑
−
=
(11.5)
若 )(
θ
x 在 内连续,则(11.1)和(11.3)式一致收敛于)1,0[ )(
θ
x 的值,即
0lim
=
∞→
k
k
d (11.6)
因此,对给定的任意小的 0>
ε
,都能找到一个 N,使 k>N 时,
ε
<
k
d 。也就是 N 以后的各
项可以忽略,即(11.4)和(11.5)式是足够准确的。
第二节 离散沃尔什变换的定义
一、按沃尔什排列的离散沃尔什变换(DWT)
W
(DWT)
W
可直接由前 N 项沃尔什函数表示式(11.4)导出,其系数由(11.2)式的前 N 个决定。
在这里,把 ,则由(11.2)式 )(kWd
xk
换写成
1,,1,0 ),()()(
1
0
−==
∫
NkdkWalxkW
Wx
"
θθθ
(11.7)
由(11.4)式
∑
−
=
=
1
0
),()()(
N
k
Wx
kWalkWx
θθ
(11.8)
上两式中的最高列率为 N/2。根据列率域的取样定理,把
θ
的定义区间 N 等分,每一
等份的时隙为
)1,0[
N1=Δ
θ
,且令
() ( ) ( )
j
j
n
xn xn x
θ
θ
θθ
=
⋅Δ
=⋅Δ=
(,) (, ) (, )
j
j
WW W
n
Walkn Walkn Walk
θ
θ
θθ
=
⋅Δ
=⋅Δ=
式中,
1,,1,0
−
== Njn "
则(11.7)变成
∑
=
Δ=
1-N
0n
),()()(
θ
nkWalnxkW
Wx
代入 N1=Δ
θ
,得
- 2 -

离散正交变换导论
1,,1,0 ),()(
1
)(
1-N
0n
−=⋅=
∑
=
NknkWalnx
N
kW
Wx
" (11.9a)
1,,1,0 ),()()(
1
0
−==
∑
−
=
NnnkWalkWnx
N
k
Wx
" (11.9b)
(11.9a)
和(11.9b)两式构成一对沃尔什正、反变换。
说明:离散沃尔什变换 DWT 与离散傅里叶变换 DFT 非常相似,其区别在于用列率域代
替了频域:(1)以权函数 代替了 ;(2)以变换系数 代替了 X(k);(3)因子 1/N
的位置不同,前者在逆变换中,后者在正变换中。
),( nkWal
kn
N
W )(kW
x
)(
kW
x
实质上是离散沃尔什级数前 N 项的系数。以 表示 阶矩阵,其中 p
为正整数,则(11.9)式可以写成矩阵形式
)]([
pWal
W
p
N 2=
[]
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
⋅=
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
)1(
)1(
)0(
)(
1
)1(
)1(
)0(
Nx
x
x
pWal
N
NW
W
W
W
x
x
x
#
#
(11.11a)
[]
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
⋅=
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
)1(
)1(
)0(
)(
)1(
)1(
)0(
NW
W
W
pWal
Nx
x
x
x
x
x
W
#
#
(11.11b)
以(10.28)式代入(11.9)式可以得到(DWT)
W
的指数形式
1,,1,0 )1)((
1
)1)((
1
)(
1-N
0n
)(
1-N
0n
)(
1
1
1
1
−=
∑
−⋅=
∑
−⋅=
∑∑
=
⊕⋅
=
⋅
=
−−
=
−−
Nknx
N
nx
N
kW
p
m
mmmp
p
m
mmp
kknkgn
x
" (11.12a)
1,,1,0 )1)(()1)(()(
1
0
)(
1
0
)(
1
1
1
1
−=
∑
−=
∑
−=
∑∑
−
=
⊕⋅
−
=
⋅
=
−−
=
−−
NnkWkWnx
N
k
kkn
x
N
k
kgn
x
p
m
mmmp
p
m
mmp
" (11.12b)
式中, ,p 为正整数
p
N 2=
位位二进制表示式的第的是序号 )1(
−
−
mpnn
mp
,见(10.26)式;
位位格雷码的第的是序号 1)(
1
−
−
mpkkg
m
01
1
=
=
−
− pmm
kpmmmpkkk 时,位。当、位二进制表示式的第的是序号、
二、按阿达玛排列的离散沃尔什变换(DWT)
H
由于按阿达玛排列的沃尔什函数 实际上是阿达玛矩阵从上而下按行编号(k)所
得,而阿达玛矩阵有简单的递推关系,即从低阶阿达玛矩阵很容易得到高阶阿达玛,因此按
阿达玛排列的沃尔什变换应用得更广泛一些。
),( nkWal
H
- 3 -

离散正交变换导论
- 4 -
按阿达玛排列的沃尔什变换
) 与按沃尔什排列的沃尔什变换 )n 只是排
列次序不同,本质上是一样的,所以只要以
)n 代替(DWT)
W
定义式中的 )n 就
可以得到(DWT)
H
的定义式。
,( nkWal
H
,(kWal
W
(Wal
W
,(kWal
H
,k
1,,1,0 ),()(
1
)(
1-N
0n
−=⋅=
∑
=
NknkWalnx
N
kB
Hx
" (11.13a)
1,,1,0 ),()()(
1
0
−==
∑
−
=
NnnkWalkBnx
N
k
Hx
" (11.13b)
同样,以阿达玛矩阵 )]([][)]([ ppWal
NH
HH
=
= , ,p 为正整数,代替(11.11)式中
的 ,可以得到(DWT)
H
的矩阵形式定义式
p
N 2=
)]([ pWal
W
[]
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
⋅=
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
)1(
)1(
)0(
1
)1(
)1(
)0(
Nx
x
x
N
NB
B
B
N
x
x
x
#
#
H
(11.14a)
[]
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
⋅=
⎥
⎥
⎥
⎥
⎦
⎤
⎢
⎢
⎢
⎢
⎣
⎡
−
)1(
)1(
)0(
)1(
)1(
)0(
NB
B
B
Nx
x
x
x
x
x
N
#
#
H
(11.14b)
以(10.30)式代入(11.13)式可以得到(DWT)
H
的指数形式定义式
1,,1,0 )1)((
1
)(
1-N
0n
1
11
−=
∑
−⋅=
∑
=
⋅
=
−−
Nknx
N
kB
p
m
mm
kn
x
" (11.15a)
1,,1,0 )1)(()(
1
0
1
11
−=
∑
−=
∑
−
=
⋅
=
−−
NnkBnx
N
k
kn
x
p
m
mm
" (11.15b)
式中, ,p 为正整数,
p
N 2=
位位二进制表示式的第的、是序号、 1
11
−
−−
mpnkkn
mm
)(kB
x
实质上是离散沃尔什级数前 N 项的系数,常称为变换系数。
)(kB
x
所对应的沃尔什分量的列率简称 的列率,以符号 i
B
表示。下面从)(kB
x
),(
θ
kWal
H
与
),(
θ
k
W
Wal 的关系来说明 i
B
。由(10.9)式
若 ,则 ) ,(),((
WH
kgk >=< )
θ
θ
WWHH
kWalkWal
=
其中, 是序号 的 p 位二进制表示式的码位倒置后的表示式, 是序号 的 p 位
格雷码,其中 ,p 为正整数。
><
H
k
H
k
p
)(
W
kg
W
k
N 2=

离散正交变换导论
令 b( )为( )内的格雷码至二进码的转换,则
WWH
kkgbkb
=
=
>
<
))(()( ,则
)),((),(
θ
θ
>
<
=
kbWalkWal
WH
(11.16)
根据 )),((
θ
>< kbWal
W
与序号 的关系可得 )( >< kb
⎪
⎪
⎩
⎪
⎪
⎨
⎧
><
+><
><><
=
==
为奇数
为偶数的列率
)kb(
2
1)kb(
)kb( )/2kb(
0 0
),(
k
kWali
HB
θ
(11.17)
以 为例,依照(11.17)式计算出的 的列率 i 如表 11.1 所示。 82 ,3
3
=== Np )(kB
x
表 11.1 (DWT)
H
的变换系数 的列率 i )(kB
x
(DWT)
H
系数
k k
2
k
1
k
0
<k> b(<k>) b(<k>)
十
i
B
)0(
x
B 0 0 0 0 000 000 0 0
)1(
x
B 1 0 0 1 100 111 7 4
)2(
x
B 2 0 1 0 010 011 3 2
)3(
x
B 3 0 1 1 110 100 4 2
)4(
x
B 4 1 0 0 001 001 1 1
)5(
x
B 5 1 0 1 101 110 6 3
)6(
x
B 6 1 1 0 011 010 2 1
)7(
x
B 7 1 1 1 111 101 5 3
第三节 阿达玛排列的快速沃尔什变换
沃尔什排列的快速沃尔什变换是从阿达玛排列的快速沃尔什变换修改而来的,所以首
先介绍阿达玛排列的快速沃尔什变换。
阿达玛排列的快速沃尔什变换以符号(FWT)
H
表示。
我们从(11.14a)式导出(FWT)
H
的两种算法流图。
一、第一种算法
取 N=8 为例,则根据(11.14a)式
7,,2,1,0 )](][[
8
1
)]([
8
"== k,nnxkB
x
H (11.18)
分解 8 阶阿达玛矩阵 ][
8
H
- 5 -
剩余21页未读,继续阅读





安全验证
文档复制为VIP权益,开通VIP直接复制

评论0