利索能及
我要发布
收藏
专利号: 2019105420948
申请人: 杭州电子科技大学
专利类型:发明专利
专利状态:已下证
更新日期:2026-06-16
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.两通道正交镜像滤波器组的多相结构,其特征在于:包括H0,M(z)模块、q个 模块、模块连接部分、1个总输入端口Input、总输出端口Output0和总输出端口Output1,i=1,

2,…q;

模块包括1个乘加复用单元、1个Xi偶数单元和1个Xi奇数单元;Xi偶数单元包括乘法器xi,2、乘法器xi,4、……、乘法器xi,d、d-2个延时为t的延时器和d-2个加法器,t=2T;乘法器xi,2、乘法器xi,4、……、乘法器xi,d的输入端均与总输入端口Input相连接;乘法器xi,2、乘法器xi,4、……、乘法器xi,d的输出信号依次通过d/2-1个加法器叠加输出至端口Ki0;乘法器xi,d、乘法器xi,d-2、……、乘法器xi,2的输出信号依次通过d/2-1个加法器叠加输出至延时为(Rd+2M)*T的延时器1i,延时器1i输出至端口Ii1;Xi偶数单元内各加法器的输入端处均设置有延时为t的延时器;

Xi奇数单元包括乘法器xi,1、乘法器xi,3、……、乘法器xi,d-1、d-2个延时为t的延时器和d-2个加法器;乘法器xi,1、乘法器xi,3、……、乘法器xi,d-1的输入端均与总输入端口Input相连接;乘法器xi,1、乘法器xi,3、……、乘法器xi,d-1的输出信号依次通过d/2-1个加法器叠加输出至端口Ki1;乘法器xi,d-1、乘法器xi,d-3、……、乘法器xi,1的输出信号依次通过d/2-1个加法器叠加输出至延时为(Rd+2M)*T的延时器2i,延时器2i输出至端口Ii0;Xi奇数单元内各加法器的输入端处均设置有延时为t的延时器;

乘加复用单元包括乘法器yi,R、乘法器yi,R-1、……、乘法器yi,1、2R-2个延时为T1的延时器(T1=d T)、2R-1个加法器、输入端口Ki、输入端口Ii以及输出端口Fi;

输入端口Ki的输入信号经过j个延时为T1的延时器,输入端口Ii输入信号经过R-1-j个延时为T1的延时器,共同输出至同一个加法器,该加法器输出至乘法器yi,R-j,j=0,1,2,…,R-1;乘法器yi,R、乘法器yi,R-1、……、乘法器yi,1的输出信号依次通过R-1个加法器叠加输出至端口Fi;

H0,M(z)模块分为M为偶数和M为奇数两种情况;

M为偶数的情况下,H0,M(z)模块包括1个HM偶数单元、1个HM奇数单元、输出端口H0,M0(z2)和输出端口H0,M1(z2);

HM偶数单元包括乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-2)、M-2个延时为t的延时器和M-2个加法器;乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-2)的输入端连接在一起,并通过延时为Rd*T的延时器与总输入端口Input相连接;乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-2)的输出信号依次通过M/2-1个加法器和延时器叠加输出至端口S1;乘法器h'(M-2)、乘法器h'(M-4)、……、乘法器h'(0)的输出信号依次通过M/2-1个加法器和延时器叠加输出至第一个延时为M T的延时器,第一个延时为M T的延时器输出至端口S0;

HM奇数单元包括乘法器h'(1)、乘法器h'(3)、……、乘法器h'(M-1)、M-2个延时为t的延时器和M-2个加法器;乘法器h'(1)、乘法器h'(3)、……、乘法器h'(M-1)的输入端均与乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-2)的输入端连接在一起;乘法器h'(1)、乘法器h'(3)、……、乘法器h'(M-1)的输出信号依次通过M/2-1个加法器和延时器叠加输出至端口P0;乘法器h'(M-1)、乘法器h'(M-3)、……、乘法器h'(1)的输出信号依次通过M/2-1个加法器和延时器叠加输出至第二个延时为M T的延时器;第二个延时为M T的延时器输出至端口P1;

端口P0和端口S0通过加法器,叠加输出至H0,M(z)模块的输出端口H0,M0(z2);端口P1和端口S1通过加法器,叠加输出至H0,M(z)模块的输出端口H0,M1(z2);

M为奇数的情况下,H0,M(z)模块包括1个HM偶数单元、1个HM奇数单元、输出端口H0,M0(z2)和输出端口H0,M1(z2);

HM偶数单元包括乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-1)、M-1个延时为t的延时器和M-1个加法器;乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-1)的输入端连接在一起,并通过延时为Rd*T的延时器与总输入端口Input相连接;乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-1)的输出信号依次通过(M-1)/2个加法器和延时器叠加输出至端口S1;乘法器h'(M-1)、乘法器h'(M-3)、……、乘法器h'(0)的输出信号依次通过(M-1)/2个加法器和延时器叠加输出至延时为(M-1)T的延时器,延时为(M-1)T的延时器输出至端口S0;

HM奇数单元包括乘法器h'(1)、乘法器h'(3)、……、乘法器h'(M-2)、M-3个延时为t的延时器和M-3个加法器;乘法器h'(1)、乘法器h'(3)、……、乘法器h'(M-2)的输入端与乘法器h'(0)、乘法器h'(2)、……、乘法器h'(M-1)的输入端连接;乘法器h'(1)、乘法器h'(3)、……、乘法器h'(M-2)的输出信号依次通过(M-3)/2个加法器和延时器叠加输出至端口P0;乘法器h'(M-2)、乘法器h'(M-4)、……、乘法器h'(1)的输出信号依次通过(M-3)/2个加法器和延时器叠加输出至延时为(M+1)T的延时器;延时为(M+1)T的延时器输出至端口P1;

端口P0和端口S0通过加法器,叠加输出至H0,M(z)模块的输出端口H0,M1(z2);端口P1和端口S1通过加法器,叠加输出至H0,M(z)模块的输出端口H0,M0(z2);

所述的模块连接部分包括2个降2采样器、q+1个加法器和2个延时为T的延时器;将q个模块中乘加复用单元的输出端口Fi的输出信号依次通过q-1个加法器叠加输出后,分成相同的两路信号,一路通过第一个降2采样器,输出至端口H0,H0(z2);另一路通过延时为T的延时器和第二个降2采样器,输出至端口H0,H1(z2);端口H0,H0(z2)与端口H0,M0(z2)经过加法器叠加输出至总输出端口Output0;端口H0,H1(z2)与端口H0,M1(z2)经过加法器叠加输出至总输出端口Output1,端口H0,M1(z2)与对应的加法器输入端之间连接有延时为T的延时器;

M、d、R均为非负整数,d为偶数,且满足M+Rd=N/2,N为原型滤波器阶数;T为系统单位时延。

2.根据权利要求1所述的两通道正交镜像滤波器组的多相结构,其特征在于:M等于12、d等于4,R等于7。

3.根据权利要求1所述的两通道正交镜像滤波器组的多相结构,其特征在于:q等于2。

4.如权利要求1所述的两通道正交镜像滤波器组的多相结构的系数设计方法,其特征在于:步骤一、确定原型滤波器阶数N,阻带截止频率fs,频率点数LF=8N;

步骤二、设定原型滤波器初始系数h0=[h0(0),h0(1),...,h0(N-1)],简化系数h'(n)=[h'(0),h'(1),...,h'(N/2-1)]=[h0(N/2),h0(N/2+1),...,h0(N-1)],n=0,1,...,N/2-1;

步骤三、将h'(n)表示成一个d×R的系数矩阵H和一个行向量HM的形式如式(1a)、(1b)所示;

HM=[h'(0),h'(1),...,h'(M-1)]      (1b)

步骤四、将系数矩阵H进行LU分解,系数矩阵H等价于矩阵L和矩阵U的乘积,其中矩阵L、矩阵U分别是下三角、上三角矩阵,H=LU;矩阵L、矩阵U的表达式如式(2a)、(2b)所示;

式(2a)、(2b)中,r表示系数矩阵H的秩r(H)的大小;

步骤五、取矩阵L的前q列和矩阵U的前q行计算新系数矩阵H'如式(3)所示;

H'=[X1 X2 ... Xq]·[Y1 Y2 ... Yq]T=X·YT                   (3)步骤六、将X、Y和HM三部分系数作为初始值,利用基于线性化的迭代梯度搜索技术进行优化,具体优化方法如下:

6-1.定义一个列向量h(z)如式(4)所示;

h(z)=[HM,H”]T                              (4)式(4)中,H”是矩阵H'中的元素按列排成行的行向量;

原型滤波器H0(z)的频率响应H0(w h(z))如式(5)所示

式(5)中, 表

示[0,π]内的一系列等间隔的离散频率点,即 i=1,2,…,LF;

6-2.根据式(5)定义全部频带点上的重构误差er(h(z))如式(6)所示

式(6)中,H0(w+πh(z))=e-j(w+π)(N-1)/2·cT(w+π)·h(z);

定义全部阻带频率点上的阻带衰减Hs(h(z))如式(7)所示;

6-3.设定迭代限制步长ε,迭代终止系数ε1,迭代次数k,k的初始值为1;将式(4)得到的h(z)作为优化设计迭代的系数初始值h(z0);

6-4.第k-1次迭代时的重构误差er(h(zk-1))如式(8a)所示,第k-1次迭代时的阻带衰减Hs(h(zk-1))如式(8b)所示式(8)中,h(zk-1)表示第k-1次迭代的系数zk-1的映射向量,hin(zk-1)表示第k-1次迭代的辅助系数hin(z),每次迭代中的辅助系数hin(z)由式(9)得到;

式(9)中,i是线性平滑系数;

6-5.定义

由一阶泰勒近似,得出第k次迭代时的重构误差er(h(zk))如式(10)所示

式(10)中, 表示第k次迭代时z的增量

定义

由一阶泰勒近似,得出第k次迭代时的阻带衰减Hs(h(z))如式(11)所示

6-6.建立关于系数增量 的凸优化问题,如式(12a)至式(12e)所示;

Subject to:

zH(z|0)==0    (12c)

zH(z|1)==1    (12d)

式(12a)至式(12e)中,||·||∞表示无穷范数运算,||·||2表示2范数运算,α是一个加权常数值;H(z|0)表示系数z中等于0的元素的位置集合,H(z|1)表示系数z中等于1的元素的位置集合;

求出式(12a)至式(12e)中第k次迭代时的系数增量 则第k次迭代的系数

同时更新hin(zk+1)=(1-i)·h(zk)+i·hin(zk);进入步骤6-7;

6-7.设置一个迭代终止条件如式(13)所示

式(13)中,δk表示式(12a)至式(12e)中第k次迭代得出的目标函数最小值δ;

如果(13)式不成立,则将k增大1,并执行6-4到6-6;若式(13)成立,则终止迭代,进入步骤七;

步骤七、用步骤六所获得的HM、X和Y,将HM中的具体系数h'(0)、h'(1)、……、h'(M-3)、h'(M-1)分别作为H0,M(z)模块中的乘法器h'(0)、乘法器h'(1)、……、乘法器h'(M-3)、乘法器h'(M-1)的系数;将X中的具体系数xi,1、xi,2、……、xi,d-1、xi,d分别作为 模块中的乘法器xi,1、乘法器xi,2、……、乘法器xi,d-1、乘法器xi,d的系数,将Y中的具体系数yi,R、yi,R-1、……、yi,1分别作为 模块中的乘法器yi,R、乘法器yi,R-1、……、乘法器yi,1的系数。

5.根据权利要求4所述的两通道正交镜像滤波器组的多相结构的系数设计方法,其特-3

征在于:步骤6-3中迭代限制步长ε等于1,迭代终止系数ε1等于10 。

6.根据权利要求4所述的两通道正交镜像滤波器组的多相结构的系数设计方法,其特征在于:步骤6-4中线性平滑系数i等于0.4。

7.根据权利要求4所述的两通道正交镜像滤波器组的多相结构的系数设计方法,其特征在于:步骤6-6中α等于0.102。

8.根据权利要求4所述的两通道正交镜像滤波器组的多相结构的系数设计方法,其特征在于:步骤6-6中,第k次迭代时的系数增量 借助matlab软件的CVX工具箱求出。