快速哈达玛变换.docVIP

快速哈达玛变换基本理论2.1 沃尔什--哈达玛变换沃尔什--哈达玛变换因实现简单且性能比较好而得到广泛的应用；在数字水印中就有很好的应用，在数字水印技术中水印的生成基于m序列和快速沃尔什-哈达玛变换矩阵嵌入与检测时使用了一个私钥来产生伪随机序列作为原始水印，使得攻击者在没有密钥的情况下无法取得水印的信息，增强了保密性 [6] 。2.1.1 哈达玛变换哈达玛变换（FHT）是数字信号处理中的基本变换之一，在移动通信、多媒体编解码中得到了广泛的应用。设x(n),n=0,1, ...,N-1，为一实序列，其离散哈达玛变换（DHT）定义为???????? (2-1)将XH(k)分解为奇对称分量XHo(k)与偶对称分量XHe(k)之和，其次，可以把奇偶对称分量表示出来：?????????????????? (2-2)由DHT定义把XH（k）代入（2-2）式,可得到奇偶对称分量表示如下：?????????????????????? (2-3)仿照快速DFT的分解方法，可通过时域抽取或频域抽取的方式实现快速DHT。 x(n)的N=2M点DHT如下式：???? (2-4)2.1.2? 沃尔什--哈达玛变换???? 沃尔什函数是二值正交函数,它仅有可能的取值是+1和-1(或0和1),适合于数字信号的处理,可作为码分多址通信系统的地址码使用。沃尔什函数记作Ｗａｌ(ｎ,ｔ),其中ｎ称一般序数,ｔ是时间变数。按由小到大排列的前八个沃尔什函数分别是:Ｗａｌ(0,ｔ)、Ｗａｌ(1,ｔ)、Ｗａｌ(2,ｔ)、Ｗａｌ(3,ｔ)、Ｗａｌ(4,ｔ)、Ｗａｌ(5,ｔ)、Ｗａｌ(6,ｔ)、Ｗａｌ(7,ｔ)。可用一个8×8阶矩阵表示如下:?????????????????? (2-6)一般沃尔什函数矩阵记为Ｗ,它是一个Ｎ×Ｎ阶实数方阵。沃尔什函数也可用哈达玛矩阵Ｈ表示,Ｈ是一个正交方阵,它的每一行代表着一个沃尔什函数。二阶哈达玛矩阵为:?或???????????????????????????? (2-7)四阶、八阶哈达玛矩阵分别为:? 与?????????? (2-8)?从而可得 ???????????????????????????? (2-9)可见哈达玛矩阵Ｈ容易利用递推关系来构造,Ｈ中行的序数称为哈达玛序数,记作 。它不同于一般序数ｎ,而与所取矩阵的行数有关系。1．哈达玛到沃尔什的变换ｎ与 之间的关系哈达玛矩阵的行数Ｎ确定后,一般序数ｎ与哈达玛序数 的特定关系,可用格雷码推演出来。若Ｈ的行数为Ｎ,格雷码的位数ｍ随之确定(因Ｎ=2ｍ),则ｎ的格雷码为ｇ(ｎ)=(ｇｍ-1ｇｍ-2…ｇ1ｇ0 ),将ｇ(ｎ)中的二进制位上的数自右向左(或自左向右)翻转过来,得到格雷码的反序,其值就等于 的值,即:(ｇ0ｇ1…ｇｍ-2ｇｍ-1)2= ,Ｗａｌ(ｎ,ｔ)=? (ｔ)例如:Ｎ = 8 、ｍ = 3 、ｎ= 1时,ｇ(ｎ) = ( 001 ),则有? =( 100 )2 = 4 ,即Ｗａｌ(1,ｔ) = ｈ4(ｔ)，Ｈ8的第5行对应着沃尔什函数Ｗａｌ(1,ｔ)。而当Ｎ= 64 ｍ = 6 ｎ= 1时,ｇ(ｎ) = ( 000001 ) ,则? = ( 100000 )2 = 32 。即Ｗａｌ(1,ｔ) = ｈ32(ｔ),Ｈ64的第33行对应着沃尔什函数Ｗａｌ(1,ｔ)。当Ｈ的行数Ｎ = 4,8,16,32,64时,用此方法推算出的ｎ与 的对应关系如表2-1。从表中数据可知当0≤ｎＮ/2时,ｎ的值对应 的偶数;当Ｎ/2≤ｎＮ时,ｎ的值对应 的奇数。表2-1　Ｎ=4,8,16时ｎ与 的对应关系序数?对应数值n?0?? 1?? 2?? 3?? 4?? 5?? 6?? 7?? 8?? 9?? 10?? 11?? 12?? 13?? 14?? 15?0?? 2?? 3?? 1?0?? 4?? 6?? 2?? 3?? 7?? 5?? 1?0?? 8?? 12? 4?? 6?? 14? 10? 2?? 3?? 11?? 15?? 7??? 5??? 13?? 9??? 1 若将表2-1中的 先以从小到大的顺序排出,再得出ｎ的对应值如表2-2。表2-2　 =4,8,16时 与n的对应关系序数?对应数值?0?? 1?? 2?? 3?? 4?? 5?? 6?? 7?? 8?? 9?? 10?? 11?? 12?? 13?? 14?? 15?0?? 3?? 1?? 2?0?? 7?? 3?? 4?? 1?? 6?? 2?? 5?0?? 15?? 7?? 8?? 3?? 12? 4?? 11?? 1?? 14?? 6?? 9??? 2??? 13?? 5??? 10 2．沃尔什到哈达玛的变换Ｗ与Ｈ的相互转化方