专利摘要
专利摘要
本发明属于数字信号处理领域,公开了一种设计有限长脉冲响应数字低通滤波器的方法,包括:设置理想低通滤波器的频率响应函数;设置模拟低通滤波器的频率响应函数;设置粒子群算法的适应度函数;得到初始种群粒子的位置矩阵和初始种群粒子的速度矩阵;得到初始种群中每个粒子的最小适应度函数各自对应的粒子位置,以及初始种群中全局最小适应度函数对应的粒子位置;对速度矩阵、位置矩阵、每个粒子的最小适应度函数对应的粒子位置,以及全局最小适应度函数对应的粒子位置进行更新;根据全局最小适应度函数对应的粒子位置得到最终的模拟低通滤波器,使设计的滤波器有更为平坦的通带,而且拥有更强的阻带衰减。
权利要求
1.一种设计有限长脉冲响应数字低通滤波器的方法,其特征在于,所述方法包括如下步骤:
步骤1,设置理想低通滤波器的参数,根据所述理想低通滤波器的参数得到理想低通滤波器的频率响应函数;具体包括如下子步骤:
(1a)设置理想低通滤波器的参数,所述理想低通滤波器的参数至少包含通带内的最大允许衰减α
(1b)根据所述理想低通滤波器的参数,得到理想低通滤波器的频率响应函数H
其中,m表示采样点的索引,其范围为{1,2,…,M},
步骤2,设置模拟低通滤波器的参数,以及粒子群算法的算法参数;根据所述模拟低通滤波器的参数和所述粒子群算法的算法参数确定模拟低通滤波器的频率响应函数;
具体包括如下子步骤:
(2a)设置模拟低通滤波器的阶数为N-1阶,N>1;
(2b)设置粒子群算法的算法参数,所述粒子群算法的算法参数至少包含种群粒子个数N
其中,w_min为第k次迭代的惯性权重因子wei
(2c)计算得到第k次迭代中第i个粒子位置矢量
其中,
步骤3,设置粒子群算法的适应度函数,所述粒子群算法的适应度函数为所述模拟低通滤波器的采样点幅度与所述理想低通滤波器的采样点幅度的误差平方和;
步骤4,根据所述粒子群算法的算法参数和所述模拟低通滤波器的参数,得到初始种群粒子的位置矩阵X
初始化:令k∈{1,2,...,MaxDT},k的初值为1,MaxDT为最大迭代次数;
步骤5,在第k次迭代中,根据所述粒子群算法的算法参数、所述粒子群算法的适应度函数对速度矩阵进行更新,得到第k次迭代的速度矩阵V
步骤6,根据所述粒子群算法的适应度函数计算第k次迭代中每个粒子的适应度函数,从而得到第k次迭代中所有粒子的最小适应度函数对应的粒子位置pbest
步骤7,设置模拟低通滤波器的精度,若第k次迭代中全局最小适应度函数小于所述模拟低通滤波器的精度,则将第k次迭代中全局最小适应度函数对应的粒子位置gbest
否则,令k的值加1,依次重复执行步骤5和步骤6;
直到当k=MaxDT时,将第MaxDT次迭代中全局最小适应度函数对应的粒子位置gbest
根据所述最优解对应的粒子位置,计算得到最终的模拟低通滤波器,并将所述最终的模拟低通滤波器作为设计的有限长脉冲响应数字低通滤波器。
2.根据权利要求1所述的一种设计有限长脉冲响应数字低通滤波器的方法,其特征在于,所述步骤3具体为:
设置粒子群算法的适应度函数F,所述粒子群算法的适应度函数F为所述模拟低通滤波器的采样点幅度与所述理想低通滤波器的采样点幅度的误差平方和,记
其中,M为理想低通滤波器和模拟低通滤波器的采样点数,|H
3.根据权利要求1所述的一种设计有限长脉冲响应数字低通滤波器的方法,其特征在于,所述步骤4具体为:
根据种群粒子个数N
将初始种群中所有粒子的最小适应度函数各自对应的粒子位置记为pbest
到的最小适应度函数对应的粒子位置,其维数为1×N;初始种群中第i个粒子所找到的最小适应度函数为将初始种群中第i个粒子的N个位置值带入模拟低通滤波器的频率响应函数H(w)中,再根据此时模拟低通滤波器的频率响应函数计算得到的粒子群算法的适应度函数,
4.根据权利要求1所述的一种设计有限长脉冲响应数字低通滤波器的方法,其特征在于,所述步骤5具体为:
第k次迭代的速度矩阵V
V
其中,wei
第k次迭代的位置矩阵X
X
5.根据权利要求4所述的一种设计有限长脉冲响应数字低通滤波器的方法,其特征在于,得到第k次迭代的速度矩阵V
6.根据权利要求1所述的一种设计有限长脉冲响应数字低通滤波器的方法,其特征在于,步骤6具体包括:
第k次迭代中所有粒子的最小适应度函数对应的粒子位置为pbest
第k次迭代的第i个粒子所找到的最小适应度函数对应的粒子位置
其中,如果第k次迭代得到的第i个粒子的位置矢量
第k次迭代中全局最小适应度函数对应的粒子位置gbest
其中,
得粒子群算法的适应度函数F最小的元素φ;如果gbest′
说明书
技术领域
本发明属于数字信号处理技术领域,尤其涉及一种设计有限长脉冲响应数字低通滤波器的方法。
背景技术
有限长脉冲响应数字滤波器(FIR),是数字信号处理系统中最基本的元件,它可以在保证任意幅频特性的同时具有严格的线性相频特性,同时其单位抽样响应是有限长的,因而滤波器是稳定的系统。在通信、雷达、声纳、生物医学、地震勘探等方面有广泛的应用。传统的FIR数字滤波器设计方法有频率抽取法、窗函数法和一致逼近法。频率抽取法和窗函数法简单易行,但不易精确确定其通带和阻带的边界频率。一致逼近法能较好的获得好的通带和阻带性能,并能准确地指定通带和阻带的边缘,是一种有效的设计方法。近年来,有些学者提出了用优化算法来设计FIR数字滤波器,它是在一定的优化准则下,使设计的滤波器性能最优。基于混合整数规划技术的设计方法就是一种全局优化算法,是目前在最大误差最小化意义上能保证全局最优解的唯一解法,但是其运算量随滤波器阶数呈指数增加。基于局域邻域搜索技术的算法的主要特点是运算量小,但只能找到局部最优点,优化性能不理想。
发明内容
针对上述现有技术存在的不足,本发明的目的在于提供一种设计有限长脉冲响应数字低通滤波器的方法,通过对FIR滤波器参数不断跟踪两个极值(两个极值分别为通带最大衰减和阻带最小衰减),进而快速找到最优FIR滤波器参数设计,不但使设计的滤波器有更为平坦的通带,而且拥有更强的阻带衰减。
为达到上述目的,本发明采用如下技术方案予以实现。
一种设计有限长脉冲响应数字低通滤波器的方法,所述方法包括如下步骤:
步骤1,设置理想低通滤波器的参数,根据所述理想低通滤波器的参数得到理想低通滤波器的频率响应函数;
步骤2,设置模拟低通滤波器的参数,以及粒子群算法的算法参数;根据所述模拟低通滤波器的参数和所述粒子群算法的算法参数确定模拟低通滤波器的频率响应函数;
步骤3,设置粒子群算法的适应度函数,所述粒子群算法的适应度函数为所述模拟低通滤波器的采样点幅度与所述理想低通滤波器的采样点幅度的误差平方和;
步骤4,根据所述粒子群算法的算法参数和所述模拟低通滤波器的参数得到初始种群粒子的位置矩阵X
步骤5,在第k次迭代中,根据所述粒子群算法的算法参数、所述粒子群算法的适应度函数对速度矩阵进行更新,得到第k次迭代的速度矩阵V
步骤6,根据所述粒子群算法的适应度函数计算第k次迭代中每个粒子的适应度函数,从而得到第k次迭代中每个粒子的最小适应度函数对应的粒子位置pbest
步骤7,设置模拟低通滤波器的精度,若第k次迭代中全局最小适应度函数小于所述模拟低通滤波器的精度,则将第k次迭代中全局最小适应度函数对应的粒子位置gbest
否则,令k的值加1,依次重复执行步骤5和步骤6;
当k=MaxDT时,将第MaxDT次迭代中全局最小适应度函数对应的粒子位置gbest
根据所述最优解对应的粒子位置计算得到最终的模拟低通滤波器,并将所述最终的模拟低通滤波器作为设计的有限长脉冲响应数字低通滤波器。
与现有技术相比,本发明具有以下优点:
(1)本发明方法选用粒子群作为获得性能最好的FIR数字滤波器算法,其原理简单,可调控制参数比较少,因而简单的参数控制策略就能取得满意的优化效果,实施随机、并行、直接的全局搜索,易于求解和实现;(2)本发明采用更新惯性权重因子的方法来避免后期出现的粒子在全局最优解附近“震荡”,即让惯性权重因子从最大的w_max线性减小到w_min;(3)通过先验信息,对速度和位置矩阵进行边界控制,可以有效调整其搜索方向,提高收敛速度,同时又能有效防止算法陷入局部最优。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种设计有限长脉冲响应数字低通滤波器的方法的流程示意图;
图2为本发明实施例提供的方法设计的滤波器与等波纹一致逼近法设计的滤波器的对比示意图;
图3为本发明实施例提供的方法设计的滤波器的幅频特性与相频特性曲线示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种设计有限长脉冲响应数字低通滤波器的方法,如图1所示,所述方法包括如下步骤:
步骤1,设置理想低通滤波器的参数,根据所述理想低通滤波器的参数得到理想低通滤波器的频率响应函数。
所述步骤1具体包括如下子步骤:
(1a)设置理想低通滤波器的参数,所述理想低通滤波器的参数至少包含通带内的最大允许衰减αp,阻带内的最小允许衰减αs,通带截止频率Ωp,阻带截止频率Ωs以及理想低通滤波器的采样点数M;
(1b)从而根据所述理想低通滤波器的参数得到理想低通滤波器的频率响应函数Hd(w):
其中,m表示采样点的索引,其范围为{1,2,…,M}, 为理想低通滤波器在通带截止频率Ωp对应的幅度值,且 为理想低通滤波器在阻带截止频率Ωs对应的幅度值,且 m为理想低通滤波器的第m个采样点数,且m∈{1,2,…,M},m1为理想低通滤波器幅度值为1的截止采样点,且 月oor[·]表示向下取整。
步骤2,设置模拟低通滤波器的参数,以及粒子群算法的算法参数;根据所述模拟低通滤波器的参数和所述粒子群算法的算法参数确定模拟低通滤波器的频率响应函数。
所述步骤2具体包括如下子步骤:
(2a)设置模拟低通滤波器的阶数为N-1阶,N>1;
(2b)设置粒子群算法的算法参数,所述粒子群算法的算法参数至少包含种群粒子个数NP,最大迭代次数MaxDT,并记个体学习因子为c1、全局学习因子为c2,记第k次迭代的第i个粒子的位置为 的维数为1×N,第k次迭代的第i个粒子的速度为Vi
其中,w_min为第k次迭代的惯性权重因子wei
(2c)从而得到第k次迭代中第i个粒子位置矢量 所对应的模拟低通滤波器的频率响应函数H(w)为:
其中, 为第k次迭代的第i个粒子的第n个位置值,n∈{1,2,…,N},w为模拟低通滤波器的频率变量;对所述模拟低通滤波器的频率响应函数进行采样,得到M个采样点{H
步骤3,设置粒子群算法的适应度函数,所述粒子群算法的适应度函数为所述模拟低通滤波器的采样点幅度与所述理想低通滤波器的采样点幅度的误差平方和。
所述步骤3具体为:
设置粒子群算法的适应度函数F,所述粒子群算法的适应度函数F为所述模拟低通滤波器的采样点幅度与所述理想低通滤波器的采样点幅度的误差平方和,记 m∈(1,2,…,M);
其中,M为理想低通滤波器和模拟低通滤波器的采样点数,|H
步骤4,根据所述粒子群算法的算法参数和所述模拟低通滤波器的参数得到初始种群粒子的位置矩阵X
所述步骤4具体为:
根据种群粒子个数NP和模拟低通滤波器的阶数N-1,随机产生Np×N维的矩阵作为初始种群粒子的位置矩阵X
将初始种群中所有粒子的最小适应度函数各自对应的粒子位置记为 且pbest
步骤5,在第k次迭代中,根据所述粒子群算法的算法参数、所述粒子群算法的适应度函数对速度矩阵进行更新,得到第k次迭代的速度矩阵V
所述步骤5具体为:
第k次迭代的速度矩阵V
Vi
其中,Vi
第k次迭代的位置矩阵X
X
需要补充的是,得到第k次迭代的速度矩阵V
步骤6,根据所述粒子群算法的适应度函数计算第k次迭代中每个粒子的适应度函数,从而得到第k次迭代中每个粒子的最小适应度函数对应的粒子位置pbest
步骤6具体包括:
第k次迭代中所有粒子的最小适应度函数对应的粒子位置为pbest
第k次迭代的第i个粒子所找到的最小适应度函数对应的粒子位置 的更新公式为:
其中, 的更新依赖于第k-1次迭代的第i个粒子的结果 第k次迭代中第i个粒子的位置矢量 以及适应度函数F,如果第k次迭代得到的第i个粒子的位置矢量 所对应的适应度函数值小于第k-1次的迭代结果 对应的适应度函数值,则用 更新 反之,则用 更新
第k次迭代中全局最小适应度函数对应的粒子位置gbest
其中, 其中 表示从Φ集合中找出使得函数F最小的元素φ,gbest
步骤7,设置模拟低通滤波器的精度,若第k次迭代中全局最小适应度函数小于所述模拟低通滤波器的精度,则将第k次迭代中全局最小适应度函数对应的粒子位置gbest
否则,令k的值加1,依次重复执行步骤5和步骤6;
当k=MaxDT时,将第MaxDT次迭代中全局最小适应度函数对应的粒子位置gbest
根据所述最优解对应的粒子位置计算得到最终的模拟低通滤波器,并将所述最终的模拟低通滤波器作为设计的有限长脉冲响应数字低通滤波器。
具体的,设置模拟低通滤波器的精度eps为0.08。
本发明提出的设计有限长脉冲响应数字低通滤波器的方法,通过对FIR滤波器参数不断跟踪两个极值(两个极值分别为通带最大衰减和阻带最小衰减),进而快速找到最优FIR滤波器参数设计,不但使设计的滤波器有更为平坦的通带,而且拥有更强的阻带衰减。
下面通过仿真实验对本发明效果作进一步验证说明。
(一)仿真条件:
本发明的仿真实验在MATLAB 7.11软件下进行的,在本发明的仿真实验的详细的参数参照表1。
表1
(二)仿真内容
图2为基于粒子群算法设计FIR数字低通滤波器与等波纹一致逼近法设计的滤波器对比图,其中横轴表示频率,纵轴表示幅度。实线表示本发明频率响应曲线,虚线表示等波纹一致逼近法的频率响应曲线。从图2中可以看出基于粒子群算法设计FIR数字低通滤波器相比等波纹一致逼近法设计的滤波器,有更为平坦的通带同时拥有很强的阻带衰减。
图3为基于粒子群算法设计FIR数字低通滤波器的幅频特性与相频特性曲线。其中图3(a)表示其幅频特性曲线,横轴表示频率,纵轴表示幅度;图3(b)表示其相频特性曲线,横轴表示频率,纵轴表示相位。从图中可以看出本发明设计的滤波器通带平坦,阻带衰减>20dB,且在通带范围内其线性相位保持较好。
综上所述,仿真实验验证了本发明的正确性,有效性和可靠性。
以上所述,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以所述权利要求的保护范围为准。
一种设计有限长脉冲响应数字低通滤波器的方法专利购买费用说明
Q:办理专利转让的流程及所需资料
A:专利权人变更需要办理著录项目变更手续,有代理机构的,变更手续应当由代理机构办理。
1:专利变更应当使用专利局统一制作的“著录项目变更申报书”提出。
2:按规定缴纳著录项目变更手续费。
3:同时提交相关证明文件原件。
4:专利权转移的,变更后的专利权人委托新专利代理机构的,应当提交变更后的全体专利申请人签字或者盖章的委托书。
Q:专利著录项目变更费用如何缴交
A:(1)直接到国家知识产权局受理大厅收费窗口缴纳,(2)通过代办处缴纳,(3)通过邮局或者银行汇款,更多缴纳方式
Q:专利转让变更,多久能出结果
A:著录项目变更请求书递交后,一般1-2个月左右就会收到通知,国家知识产权局会下达《转让手续合格通知书》。
动态评分
0.0