一类易于VLSI实现的对称双正交小波设计方法研究(完整版)实用资料.doc
《一类易于VLSI实现的对称双正交小波设计方法研究(完整版)实用资料.doc》由会员分享,可在线阅读,更多相关《一类易于VLSI实现的对称双正交小波设计方法研究(完整版)实用资料.doc(88页珍藏版)》请在淘文阁 - 分享文档赚钱的网站上搜索。
1、一类易于VLSI实现的对称双正交小波设计方法研究(完整版)实用资料(可以直接使用,可编辑 完整版实用资料,欢迎下载)第32卷第4期电子与信息学报Vol.32No.4 2021年4月 Journal of Electronics & Information Technology Apr. 2021一类易于VLSI实现的对称双正交小波设计方法研究张池军彭安金王厚军李本亮(电子科技大学自动化工程学院成都 610054(西南民族大学电气与信息工程学院成都 610041(中国人民解放军95989部队北京 100076摘要:该文提出了一类对称双正交小波的设计方法。该类双正交小波的小波滤波器组具有格形结构,
2、实现该小波变换的分析滤波器组和综合滤波器组满足双正交条件和正则性条件,且设计的各滤波器均为实数二进制系数,因而该小波变换易于高速VLSI实现。文中的理论推导和设计实例,均验证了该设计方法的有效性。关键词:对称双正交小波;格形结构;双正交滤波器组;线性相位中图分类号:TN713 文献标识码:A 文章编号:1009-5896(202104-0973-05The Design Method of a Class of Symmetric BiorthogonalWavelets Suitable for VLSI ImplementationZhang Chi-junPeng An-jin Wang
3、 Hou-junLi Ben-liang(School of Automation Engineering, University of Electronic Science and Technology of China, Chengdu 610054,China(School of Electrical and Information Engineering, Southwest University for Nationalities, Chengdu 610041, China(The 95989 Unit of PLA, Beijing 100076,ChinaAbstract: T
4、he design method of a class of symmetric biorthogonal wavelets is proposed in this paper. The filterbanks of the wavelets possess lattice structure, the analysis and synthesis filter banks for wavelets meet biorthogonality and regularly conditions, and the filters are all real binary coefficients. T
5、herefore, the wavelet transform is suitable for high-speed VLSI implementation. Both the mathematical derivations and the design examples in the paper verify the effectiveness of proposed method.Key words: Symmetric biorthogonal wavelets; Lattice structure; Biorthogonal filter banks; Linear phase1引言
6、小波变换在图像处理,信号与信息处理,地震勘探,电力系统,生物医学以及军事领域具有广阔的应用前景。而完成小波变换的前提是找到合适的小波基,因此小波基的构造方法一直是小波理论研究的热点。Daubechies 构造了紧支正交小波基1 , 尽管该小波基在信号处理领域得到了广泛应用,但紧支集的正交小波, 除Haar小波外都是非对称的。要获得对称小波,就必须使小波滤波器具有线性相位,故必须放弃小波的正交性,选择双正交性1,2。Vetterli 等研究了线性相位的理想重建滤波器组(Linear Phase Perfect Reconstruction Filter Banks(LPPRFB 和双正交小波基的
7、关系2,并提出了用LPPRFB构造双正交小波基的方法,尤其对双正交小波的实系数FIR滤波器的对称性与长度选2021-04-13收到,2021-10-16改回国家部委基金(A1420061264和国家自然科学基金(60673011资助课题通信作者:张池军 zcj1120163 取方法进行了讨论。文献3设计了一种小波滤波器组为梯形结构的小波,其小波滤波器具有线性相位和实数二进制系数,小波变换易于VLSI实现,但其设计方法不易控制小波滤波器的阶数。文献4提出了一种小波滤波器组具有格形结构的小波,但其小波滤波器不具有线性相位。深入分析发现:这不是由于采用了格形结构的缘故,而是格形结构中滤波器的设计方法
8、原因。国内学者对双正交小波的工程应用进行了探索5,也有学者对双正交小波滤波器设计开展了研究6,但未见有文献报道在构成小波的滤波器组结构方面出现突破。为此,本文在格形滤波器组的基础上,找到了格形结构中滤波器的另外一种设计方法,进而提出了一类小波滤波器组为格形结构的对称双正交小波的设计方法。该设计方法在算法上保证了格形小波滤波器组的线性相位和实数二进制系数(即k/2n,其中k,n均为整数。故该类小波的变换仅用移位寄存器和累加器就可完成,因而易于高速VLSI实现。当然,在实时性要求不苛刻时,该类小波变换更适合微控制器的软件实现。974 电 子 与 信 息 学 报 第32卷2 对称双正交小波与两通道的
9、LPPRFB双正交小波有两个互为对偶的尺度函数(t 和(t ,有两个互为对偶的小波函数(t 和(t ,其尺度函数方程和小波函数方程如下:0(100(2(2N n t h n t n =(11100(2(2N n t f n t n = (21(110(2(2N n t h n t n =(30110(2(2N n t f n t n = (4 其中0(h n 和1(h n 分别为0(h n 和1(h n 的时间反转;0(h n 和0(f n 是低通滤波器,长度分别为0N 和1N ;1(h n 和1(f n 为高通滤波器,长度分别为1N 和0N 。如果00(f k h k =,11(f k h
10、k =,则(t = (t ,(t t =,显然,双正交小波包含正交小波,但正交小波除Harr 小波外,不具有线性相位,故本文不做讨论。就多速率滤波器组而言,0(h n 和1(h n 可构成分析滤波器组,而0(f n 和1(f n 则构成综合滤波器组,二者可构成两通道的多速率滤波器组(Filter banks ,FB,本文称之为小波滤波器组。文献2,7,8对小波滤波器组与两通道LPPRFB 进行了研究,结果表明:对称双正交小波的小波滤波器组是一类特殊两通道LPPRFB ,即除了满足线性相位和PR(perfect reconstruction特性外,还必须满足正则性(regularity。下面先分
11、析两通道LPPRFB 的PR 特性,然后讨论对称双正交小波的正则性。 2.1 两通道的LPPRFB 的PR 特性分析两通道的LPPRFB 由分析滤波器组(k H z ,抽取器,0值内插器和综合滤波器组(k F z 构成,其结构如图1所示。 图1 两通道的LPPRFB 的结构两通道的LPPRFB 的PR 特性可表示为010011(20(c n H z H z F z z H z H z F z = (5 其中c n 为奇数。显然式(5包含了滤波器系数0(h n ,1(h n ,0(f n ,1(f n ,但这些滤波器系数间有何种联系?为了寻求滤波器系数0(h n ,1(h n ,0(f n ,1
12、(f n 间的关系, 根据式(5可选取0(F z 、1(F z 为01(F z H z = (6 10(F z H z = (7由式(5-式(7得0110(2c n H z H z H z H z z = (8 两通道LPPRFB 的设计,可归结为根据式(8,寻求满足设计要求的0(h n 和1(h n 。 2.2 对称双正交小波的正则性分析小波变换将信号分解成小波系数,而小波的正则性则有利于小波系数更好地表示信号。小波的正则性是对尺度函数(t 的“光滑性”要求,它可用(t 的连续性及其导数的连续性定义。该正则性与两通道LPPRFB 的PR 特性无关,但由低通滤波器(c n 迭代产生的尺度函数1
13、0(2n t c n t n = (9对(c n 的正则性要求可表示为两个约束条件9,10:(1归一化条件 将式(9两端积分后即可得归一化条件为10(1N n c n = (10(2(C z 在1z =处有零点 (c n 的Z 变换(C z 在1z =处有零点,对尺度函数(t 的连续性及其导数的连续性很有必要。但对正则性而言,并不是(C z 在1z =的零点越多,正则性越好(正则性阶数越高。但(C z 最好在1z =处至少有3个零点3。上述分析表明:对称双正交小波的设计,从多速率滤波器组的角度而言,就是设计具有正则性的两通道LPPRFB 。3 格形双正交滤波器组文献4中提出了一类格形双正交滤波
14、器组,其结构如图2所示。图中0(H z ,1(H z ,0(F z ,1(F z 为小波滤波器;0(T z ,1(T z ,0(R z ,1(R z 为传输滤波器;0(L z ,1(L z ,0(W z ,1(W z 为交叉格形滤波器。文献4的思想是采用格形网络来改善小波滤波器的选频特性。本文的目的是研究小波滤波器组具有格形网络结构的一类对称双正交小波设计方法。选择0(1H z =,11(H z z = 10(F z z =,1(1F z =,并将图2等效成图1结构,则对应的等效滤波器为l 212000(H z T z z L z =+ (11 l 212111(H z L z z T z =
15、+ (12 l 21201(Fz W z z R z =+ (13l 212110(F z R z z W z =+ (14第4期 张池军等:一类易于VLSI 实现的对称双正交小波设计方法研究 975 图2 格形双正交滤波器组的结构由于等效滤波器也应满足式(6和式(7,由式(11-式(14,可得2202121022002211(R z T z R z T z W z L z W z L z =(15 式(15表明,只要求得20(T z ,21(T z ,20(L z ,21(L z ,则综合部分(即图2中(b的格形结构便可确定。但关键是如何求得20(T z ,21(T z ,20(L z ,2
16、1(L z ?根据式(8,式(11,式(12可得0101(k T z T z L z L z z = (16 其中k 为正整数。显然式(16就是格形双正交滤波器组的理想重构条件。文献4设计的小波滤波器不具有线性相位,其重要原因是将式(16中20(T z ,21(T z ,20(L z ,21(L z 进一步约束为 01(L z T z =和10(L z T z =所致。故本文对此做出改进,放弃该约束条件,根据式(16用代数方法直接求解。4 格形对称双正交小波的设计4.1 格形对称双正交小波的设计方法小波滤波器组具有格形结构的双正交小波,本文称之为格形双正交小波。格形双正交小波的设计,其实可以归
17、结为设计具有正则性的格形双正交滤波器组,即设计的格形双正交滤波器组要同时满足正则性条件和PR 条件。等效低通滤波器0(hn 和0(f n 要满足的正则性条件分别为(10(hn 的正则性条件: l 0100(1(0,0,1,1Nn m n n hn m p = (17其中l 0N为0(h n 的长度,p 为整数,它表示l 0(H z 在z =-1处的零点个数。(20(f n的正则性条件: l 0110(0,0,1,1Nm n n hn m p = (18显然l 0(Fz 在z =-1处的1p 阶零点与l 1(H z 在z =1处的1p 阶零点等价。有了上述正则性条件式(17、式(18和PR 条件
18、式(16,就可以设计格形对称双正交小波。具体设计步骤如下:(1根据104N N k =+(k 为整数,确定0(hn 和1(hn 的长度; (2根据0(h n 和1(h n 的长度确定0(t n ,1(t n ,0(l n ,1(l n 的长度,根据0(hn 和1(h n 的对称性(奇对称或偶对称,找到0(t n ,1(t n ,0(l n ,1(l n 系数之间的关系;(3根据式(16得到01(/4N N +个方程;(4根据l l 01(Hz i H z i =得到1个方程; (5根据0(2Nn h n =得到1个方程; (6根据式(17,式(18得到若干方程;(7用MATLAB 的符号变量表
19、示所有方程,然后用solve( 函数求解。4.2 格形对称双正交小波的设计实例设计实例中,“SE AE 01N N ”表示0(hn 和1(hn 的类型分别为SE(偶对称,偶长度和AE(奇对称,偶长度,长度分别为0N 和1N 。设计实例1 设计 SE AE 6 10 类型小波。 (10(h n 和1(h n 的长度分别为6,10,类型分别为SE ,AE ;(2根据0(h n 和1(h n 的长度确定0(t n ,1(t n ,0(l n ,1(l n 的长度均分别为3,5,3,5,并用符号设定各滤波器的系数向量为= (19 根据0(h,a a a c c c d d d d d b b b b
20、b =(20式(20表明向量a 由向量c 表示,向量d 由向量b 表示,显然要求解的未知量减少到01(/2N N +个。由式(11,式(12,式(19,式(20可得+=+=(22976 电 子 与 信 息 学 报 第32卷(4根据l l 01(Hz i H z i =得到1个方程: 432102100b b b b b c c c += (23(5根据0(2Nn h n =,由式(21得 2101c c c += (24(6根据式(17,式(18得到2个方程;+=(25 (7用MATLAB 的符号变量表示所有方程,然后用solve(函数求解得0(t n =-1/8, 1, 1/8;1(t n
21、=-1/64, -1/8, -31/32, 1/8, -1/64; 0(l n =1/8, 1, -1/8;1(l n =1/64, -1/8, 31/32, 1/8, 1/64; 0(hn =-1/8, 1/8, 1, 1, 1/8, -1/8; 1(hn =1/64, -1/64, -1/8, -1/8, 31/32, -31/32, 1/8, 1/8, 1/64, -1/64;0(f n = 1/64, 1/64, -1/8, 1/8, 31/32, 31/32, 1/8, -1/8, 1/64, 1/64;1(f n = 1/8, 1/8, -1, 1, -1/8, -1/8。 用MA
22、TLAB 小波工具箱中的 bswfun(函数得到尺度函数(t 和(t ,小波函数(t 和i (t 如图3,图4所示。从图3和图4中可以看出,设计的对称双正交小波的分析尺度函数和小波函数其光滑性较好,l 0(Hz 在z =-1处有3个零点。综合尺度函数和小波函数光滑性也较好,l 0(Fz 在z =-1处也有3个零点。 设计实例2 设计 SE AE 8 12 类型小波 求解过程不再详述,得到解为0(t n =1/16, 1/16, 15/16, -1/16;1(t n =-1/128, 23/128, -67/64, -11/64, 7/128, -1/128;0(l n =-1/16, 15/1
23、6, 1/16, 1/16;1(l n =1/128, -7/128, 11/64, 67/64, -23/128,1/128。用MATLAB 小波工具箱中的 bswfun(函数得到尺度函数(t 和(t ,小波函数(t 和(t 如图5,图6所示。图5和图6表明,无论综合尺度函数和小波函数,还是分析尺度函数和小波函数,光滑性都好,而且l 0(Fz 和l 0(H z 在z =-1处,零点数目均超过3。 5 结论本文从多速率滤波器组的角度,提出了一类易于VLSI 实现的格形对称双正交小波的设计方法。该方法选择小波滤波器组具有图2所示的格形结构,设计的双正交滤波器组同时满足正则性条件和PR 条件。在求
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 一类 易于 VLSI 实现 对称 正交 设计 方法 研究 完整版 实用 资料
限制150内