CN108011615B - 一种信号处理的方法和装置 - Google Patents

一种信号处理的方法和装置 Download PDF

Info

Publication number
CN108011615B
CN108011615B CN201711420773.5A CN201711420773A CN108011615B CN 108011615 B CN108011615 B CN 108011615B CN 201711420773 A CN201711420773 A CN 201711420773A CN 108011615 B CN108011615 B CN 108011615B
Authority
CN
China
Prior art keywords
frequency response
amplitude
sequence
cic
response sequence
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201711420773.5A
Other languages
English (en)
Other versions
CN108011615A (zh
Inventor
汪啸尘
庄志
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BMC Medical Co Ltd
Original Assignee
BMC Medical Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by BMC Medical Co Ltd filed Critical BMC Medical Co Ltd
Priority to CN201711420773.5A priority Critical patent/CN108011615B/zh
Publication of CN108011615A publication Critical patent/CN108011615A/zh
Application granted granted Critical
Publication of CN108011615B publication Critical patent/CN108011615B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0219Compensation of undesirable effects, e.g. quantisation noise, overflow

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Digital Transmission Methods That Use Modulated Carrier Waves (AREA)

Abstract

本发明实施例提供了一种信号处理的方法和装置,所述方法包括:确定CIC滤波器降采样后的幅频响应;依据所述降采样后的幅频响应生成初始幅频响应序列;确定期望幅频响应序列;依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列;依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数;采用所述滤波系数进行信号补偿,以输出目标幅频响应。本实施例可以根据需求将最终输出的信号的幅频响应配置为理想状态,实现对幅频响应的灵活调节,并能够保证通带和阻带的性能。

Description

一种信号处理的方法和装置
技术领域
本发明涉及信号处理技术领域,特别是涉及一种信号处理的方法、一种信号处理的装置和一种存储介质。
背景技术
在生理信号的采集中,往往会涉及到对采样信号进行降采样的过程。为了防止在降采样过程中产生信号混叠,就需要对信号进行滤波处理。通常,可以选择CIC滤波器对信号进行降采样。
由于CIC滤波器的通带衰减较快,幅频响应不理想,-3dB带宽只有最终采样频率的四分之一,在高频段会造成一定的波形失真,因而限制了这一滤波方式在高精度信号采集中的应用。为了克服上述问题,通常需要对其进行滤波补偿。
目前,最常见的CIC补偿滤波器是ISOP补偿滤波器。ISOP补偿滤波器是一种为进行CIC补偿滤波而设计的、具有固定形式系统函数的滤波器,通过所需的通带特性调节其中的参数即可完成滤波器设计。但是,由于ISOP补偿滤波器具有固定形式的系统函数,可调节参数仅有两个,也就无法灵活地调节出最为理想的幅频响应,进而也会对阻带信号产生不必要的放大效果,削弱了阻带性能。
发明内容
鉴于上述问题,提出了本发明实施例以便提供一种克服上述问题或者至少部分地解决上述问题的一种信号处理的方法、一种信号处理的装置和相应的一种存储介质。
为了解决上述问题,本发明实施例公开了一种信号处理的方法,包括:
确定CIC滤波器降采样后的幅频响应;
依据所述降采样后的幅频响应生成初始幅频响应序列;
确定期望幅频响应序列;
依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列;
依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数;
采用所述滤波系数进行信号补偿,以输出目标幅频响应。
可选地,所述确定CIC滤波器降采样后的幅频响应的步骤包括:
获取CIC滤波器的降采样倍数、级联环节数量,以及,各个环节的延迟信息;
采用所述降采样倍数、级联环节数量,以及,各个环节的延迟信息,计算所述CIC滤波器降采样后的归一化幅频响应。
可选地,所述确定CIC滤波器降采样后的幅频响应的步骤包括:
确定所述CIC滤波器的归一化幅频响应和降采样倍数;
根据所述CIC滤波器的归一化幅频响应和降采样倍数,计算所述CIC滤波器降采样后的归一化幅频响应。
可选地,所述依据所述降采样后的幅频响应生成初始幅频响应序列的步骤包括:
确定抽样间隔;
按照所述抽样间隔,对所述CIC滤波器降采样后的归一化幅频响应进行离散化抽样,以获得初始幅频响应序列,所述初始幅频响应序列包括多个序列点。
可选地,所述确定期望幅频响应序列的步骤包括:
分别确定在各个序列点上的期望幅频响应;
根据所述各个序列点上的期望幅频响应生成期望幅频响应序列。
可选地,所述依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列的步骤包括:
分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,以获得所述CIC补偿滤波器的目标幅频响应序列。
可选地,所述依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数的步骤包括:
对所述目标幅频响应序列进行离散傅里叶逆变换,以获得时域序列,以所述时域序列作为所述CIC补偿滤波器的滤波系数。
为了解决上述问题,本发明实施例公开了一种信号处理的装置,包括:
幅频响应确定模块,用于确定CIC滤波器降采样后的幅频响应;
初始幅频响应序列生成模块,用于依据所述降采样后的幅频响应生成初始幅频响应序列;
期望幅频响应序列确定模块,用于确定期望幅频响应序列;
目标幅频响应序列计算模块,用于依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列;
滤波系数确定模块,用于依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数;
信号补偿模块,用于采用所述滤波系数进行信号补偿,以输出目标幅频响应。
可选地,所述幅频响应确定模块包括:
信息获取子模块,用于获取CIC滤波器的降采样倍数、级联环节数量,以及,各个环节的延迟信息;
第一归一化幅频响应计算子模块,用于采用所述降采样倍数、级联环节数量,以及,各个环节的延迟信息,计算所述CIC滤波器降采样后的归一化幅频响应。
可选地,所述幅频响应确定模块包括:
初始信息确定子模块,用于确定所述CIC滤波器的归一化幅频响应和降采样倍数;
第二归一化幅频响应计算子模块,用于根据所述CIC滤波器的归一化幅频响应和降采样倍数,计算所述CIC滤波器降采样后的归一化幅频响应。
可选地,所述初始幅频响应序列生成模块包括:
抽样间隔确定子模块,用于确定抽样间隔;
离散化抽样子模块,用于按照所述抽样间隔,对所述CIC滤波器降采样后的归一化幅频响应进行离散化抽样,以获得初始幅频响应序列,所述初始幅频响应序列包括多个序列点。
可选地,所述期望幅频响应序列确定模块包括:
期望幅频响应确定子模块,用于分别确定在各个序列点上的期望幅频响应;
期望幅频响应序列生成子模块,用于根据所述各个序列点上的期望幅频响应生成期望幅频响应序列。
可选地,所述目标幅频响应序列计算模块包括:
目标幅频响应序列计算子模块,用于分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,以获得CIC补偿滤波器的目标幅频响应序列。
可选地,所述滤波系数确定模块包括:
离散傅里叶逆变换子模块,用于对所述目标幅频响应序列进行离散傅里叶逆变换,以获得时域序列,以所述时域序列作为所述CIC补偿滤波器的滤波系数。
为了解决上述问题,本发明实施例公开了一种信号处理的装置,包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述计算机程序被所述处理器执行时实现上述的信号处理的方法的步骤。
为了解决上述问题,本发明实施例公开了一种计算机可读存储介质,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现上述的信号处理的方法的步骤。
与背景技术相比,本发明实施例包括以下优点:
本发明实施例,通过确定CIC滤波器在降采样后的幅频响应,从而可以依据降采样后的幅频响应生成初始幅频响应序列,然后,在确定出期望幅频响应序列后,可以进一步依据上述初始幅频响应序列和期望幅频响应序列,计算得到CIC补偿滤波器的目标幅频响应序列,通过对目标幅频响应序列进行处理,可以确定出CIC补偿滤波器的滤波系数,在采用上述滤波系数进行信号补偿后,可以输出所需的目标幅频响应。本实施例可以根据需求将最终输出的信号的幅频响应配置为理想状态,实现对幅频响应的灵活调节,并能够保证通带和阻带的性能。
其次,本发明实施例可以修正CIC滤波造成的信号幅频响应畸变,恢复真实的采集信号;与现有的CIC补偿滤波器相比本实施例还可以对补偿后信号的幅频响应曲线进行任意配置,在恢复真实信号的同时获得理想的幅频响应曲线,便于后续对信号的分析和处理。
附图说明
图1是本发明一个实施例的一种信号处理的方法的步骤流程示意图;
图2是本发明一个实施例的另一种信号处理的方法的步骤流程示意图;
图3是本发明一个实施例的一种CIC滤波器的归一化幅频响应曲线的示意图;
图4是本发明一个实施例的一种CIC滤波器降采样后的归一化幅频响应曲线的示意图;
图5是本发明一个实施例的一种初始幅频响应序列的示意图;
图6是本发明一个实施例的一种期望幅频响应序列的示意图;
图7是本发明一个实施例的一种目标幅频响应序列的示意图;
图8是本发明一个实施例的一种时域序列的示意图;
图9是本发明一个实施例的一种目标幅频响应的示意图;
图10是本发明一个实施例的又一种信号处理的方法的步骤流程示意图;
图11是本发明一个实施例的另一种期望幅频响应序列的示意图;
图12是本发明一个实施例的另一种目标幅频响应序列的示意图;
图13是本发明一个实施例的另一种时域序列的示意图;
图14是本发明一个实施例的另一种目标幅频响应的示意图;
图15是本发明一个实施例的一种信号处理的装置实施例的示意性结构框图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
参照图1,示出了本发明一个实施例的一种信号处理的方法的步骤流程示意图,具体可以包括如下步骤:
步骤101,确定CIC滤波器降采样后的幅频响应;
需要说明的是,本实施例中的滤波器可以是CIC滤波器(Cascaded integrator-comb filter,级联积分梳状滤波器),可以通过CIC滤波器对信号作降采样处理。由于CIC滤波器具有线性相位,且仅需延迟和加法运算,具有计算量少、易于实现等优点,是降采样过程中最常见的一种滤波方式。
在本发明实施例中,为了实现对CIC滤波器降采样过程中的滤波补偿,可以首先确定CIC滤波器在降采样后的幅频响应。
在本发明实施例中,可以首先获取CIC滤波器的降采样倍数、级联环节数量,以及,各个环节的延迟信息,然后采用上述降采样倍数、级联环节数量,以及,各个环节的延迟信息等,计算该CIC滤波器降采样后的归一化幅频响应。
或者,也可以通过首先确定CIC滤波器的归一化幅频响应和降采样倍数,然后再根据上述CIC滤波器的归一化幅频响应和降采样倍数,计算该CIC滤波器降采样后的归一化幅频响应。
例如,对于某个CIC滤波器,其降采样倍数为R,各个环节的延迟为M,级联环节的数量为N,则可以确定该CIC滤波器的归一化幅频响应为:
Figure BDA0001522999390000061
由于信号经过CIC滤波器降采样后,采样频率发生了变化,则在上式(1)的基础上,可以得到该CIC滤波器降采样后的归一化幅频响应为:
Figure BDA0001522999390000071
步骤102,依据所述降采样后的幅频响应生成初始幅频响应序列;
在本发明实施例中,可以通过对降采样后的归一化幅频响应进行离散化抽样的方式来获得初始幅频响应序列,初始幅频响应序列中可以包括多个序列点。
在具体实现中,可以首先确定抽样间隔,然后按照上述抽样间隔,对CIC滤波器降采样后的归一化幅频响应进行离散化抽样,从而获得初始幅频响应序列。
例如,可以以1/K为抽样间隔,对H(fd)进行离散化抽样,得到总数为K+1个的初始幅频响应序列H(k)。其中,K的数值大小可以根据实际需要,由本领域技术人员具体确定,本实施例对此不作限定。
需要说明的是,一旦K的数值确定,则在本实施例中设计的CIC补偿滤波器的补偿阶数也就确定了,即为2K+1。
步骤103,确定期望幅频响应序列;
在本发明实施例中,期望幅频响应序列是对于初始幅频响应序列中的各个序列点最终期望能够达到的幅频响应。
例如,在按照1/K的抽样间隔对H(fd)进行离散化抽样获得总数为K+1个的初始幅频响应序列H(k)后,可以分别确定在各个序列点上的期望幅频响应,然后根据上述各个序列点上的期望幅频响应生成期望幅频响应序列。期望幅频响应序列可以表示为HE(k)。
步骤104,依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列;
在本发明实施例中,目标幅频响应序列可以通过在K+1个序列点上分别对对应的期望幅频响应序列HE(k)和初始幅频响应序列H(k)求商得到。
具体地,可以分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,从而获得CIC补偿滤波器的目标幅频响应序列HC(k)。
步骤105,依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数;
在本发明实施例中,CIC补偿滤波器的滤波系数可以是一个序列,包括多个序列点上的滤波系数值。
在本发明实施例中,可以对目标幅频响应序列HC(k)进行离散傅里叶逆变换,获得时域序列hC(k),从而可以以该时域序列hC(k)作为上述CIC补偿滤波器的滤波系数。
步骤106,采用所述滤波系数进行信号补偿,以输出目标幅频响应。
在具体实现中,获得滤波系数后,可以直接利用该CIC补偿滤波器对输出信号进行滤波即可实现信号补偿,获得所需的目标幅频响应。
在本发明实施例中,通过确定CIC滤波器在降采样后的幅频响应,从而可以依据降采样后的幅频响应生成初始幅频响应序列,然后,在确定出期望幅频响应序列后,可以进一步依据上述初始幅频响应序列和期望幅频响应序列,计算得到CIC补偿滤波器的目标幅频响应序列,通过对目标幅频响应序列进行处理,可以确定出滤波器的滤波系数,在采用上述滤波系数进行信号补偿后,可以输出所需的目标幅频响应。本实施例可以根据需求将最终输出的信号的幅频响应配置为理想状态,实现对幅频响应的灵活调节,并能够保证通带和阻带的性能。
为了便于理解,下面以一个完整的示例对本发明的信号处理的方法作一介绍。
参照图2,示出了本发明一个实施例的另一种信号处理的方法的步骤流程示意图,具体可以包括如下步骤:
步骤201,确定CIC滤波器的归一化幅频响应和降采样倍数;
作为本发明的一种示例,该CIC滤波器的降采样倍数可以为64,各个环节的延迟为1,级联环节的数量为3,则参考公式(1)可知,该CIC滤波器的归一化幅频响应可以表示为:
Figure BDA0001522999390000091
如图3所示,是本发明一个实施例的一种CIC滤波器的归一化幅频响应曲线的示意图,图3中的曲线即是对应于上述公式(3)的CIC滤波器的归一化幅频响应曲线。
步骤202,根据所述CIC滤波器的归一化幅频响应和降采样倍数,计算所述CIC滤波器降采样后的归一化幅频响应;
由于信号经过CIC滤波器按照64被进行降采样后,采样频率发生了变化,则在上式(3)的基础上,可以得到CIC滤波器降采样后的归一化幅频响应为:
Figure BDA0001522999390000092
如图4所示,是本发明一个实施例的一种CIC滤波器降采样后的归一化幅频响应曲线的示意图,图4中的曲线即是对应于上述公式(4)的CIC滤波器降采样后的归一化幅频响应曲线。
步骤203,按照所述抽样间隔,对所述CIC滤波器降采样后的归一化幅频响应进行离散化抽样,以获得初始幅频响应序列,所述初始幅频响应序列包括多个序列点;
在本发明实施例中,可以取K=32,然后按照1/32的抽样间隔对上述H(fd)即图4中所示的曲线进行离散化抽样,得到如图5所示的总数为33的初始幅频响应序列H(k)。
步骤204,分别确定在各个序列点上的期望幅频响应;
步骤205,根据所述各个序列点上的期望幅频响应生成期望幅频响应序列;
在本发明实施例中,可以配置在各个序列点上的期望幅频响应为HE(k)=1,即期望将全部的幅频响应补偿至无衰减状态。
如图6所示,是本发明一个实施例的一种期望幅频响应序列的示意图。在图6中,各个序列点上的期望幅频响应序列数值均被补偿至无衰减状态。
步骤206,分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,以获得CIC补偿滤波器的目标幅频响应序列;
如图7所示,是本发明一个实施例的一种目标幅频响应序列的示意图。图7中的目标幅频响应序列即是通过对33个序列点上的期望幅频响应序列HE(k)和初始幅频响应序列H(k)求商而得到的。
步骤207,对所述目标幅频响应序列进行离散傅里叶逆变换,以获得时域序列,以所述时域序列作为所述CIC补偿滤波器的滤波系数;
在本发明实施例中,可以对图7中的目标幅频响应序列HC(k)进行离散傅里叶逆变换得到如图8所示的时域序列hC(k)的示意图,图8中的时域序列hC(k)即为上述CIC补偿滤波器的滤波系数。
步骤208,采用所述滤波系数进行信号补偿,以输出目标幅频响应。
在本发明实施例中,利用该CIC补偿滤波器对输出信号进行滤波即可实现信号补偿,获得所需的目标幅频响应。
如图9所示,是本发明一个实施例的一种目标幅频响应的示意图。在图9中,上扬的曲线为补偿的幅频响应曲线,下沉的曲线为原CIC补偿滤波器的幅频响应曲线,而中间的曲线即为滤波补偿后所需的目标幅频响应。
参照图10,示出了本发明一个实施例的又一种信号处理的方法的步骤流程示意图,具体可以包括如下步骤:
步骤1001,确定CIC滤波器的归一化幅频响应和降采样倍数;
步骤1002,根据所述CIC滤波器的归一化幅频响应和降采样倍数,计算所述CIC滤波器降采样后的归一化幅频响应;
步骤1003,按照所述抽样间隔,对所述CIC滤波器降采样后的归一化幅频响应进行离散化抽样,以获得初始幅频响应序列,所述初始幅频响应序列包括多个序列点;
在本发明实施例中,还可以实现对幅频响应曲线的任意配置。在任意配置幅频响应曲线时,其步骤1001-步骤1003与上述实施例中步骤201-步骤203类似,可以相互参阅,本实施例对此不再赘述。
步骤1004,确定期望幅频响应序列;
在本发明实施例中,在经过离散化抽样获得的33个序列点上分别给出每个序列点最终期望达到的幅频响应,从而可以构成期望幅频响应序列HE(k)。
在本实施例中,可以构造一个截止频率为0.8的理想低通幅频特性曲线。如图11所示,是本发明一个实施例的另一种期望幅频响应序列的示意图,图11中所示的曲线即为截止频率为0.8的理想低通幅频特性曲线。
步骤1005,分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,以获得CIC补偿滤波器的目标幅频响应序列;
类似地,可以对图11中各个序列点对应的期望幅频响应序列数值与初始幅频响应序列数值求商,从而得到如图12所示的目标幅频响应序列HC(k)。
步骤1006,对所述目标幅频响应序列进行离散傅里叶逆变换,以获得时域序列,以所述时域序列作为所述CIC补偿滤波器的滤波系数;
在对图12中的目标幅频响应序列HC(k)进行离散傅里叶逆变换后,可以得到如图13所示的时域序列hC(k)的示意图。
步骤1007,采用所述滤波系数进行信号补偿,以输出目标幅频响应。
然后,按照上述滤波系数对输出信号进行滤波即可实现相应的信号补偿,获得如图14所示的需要的目标幅频响应。
本发明实施例可以修正CIC滤波造成的信号幅频响应畸变,恢复真实的采集信号;与现有的CIC补偿滤波器相比本实施例还可以对补偿后信号的幅频响应曲线进行任意配置,在恢复真实信号的同时获得理想的幅频响应曲线,便于后续对信号的分析和处理。
需要说明的是,对于方法实施例,为了简单描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明实施例并不受所描述的动作顺序的限制,因为依据本发明实施例,某些步骤可以采用其他顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作并不一定是本发明实施例所必须的。
参照图15,示出了本发明一个实施例的一种信号处理的装置实施例的示意性结构框图,具体可以包括如下模块:
幅频响应确定模块1501,用于确定CIC滤波器降采样后的幅频响应;
初始幅频响应序列生成模块1502,用于依据所述降采样后的幅频响应生成初始幅频响应序列;
期望幅频响应序列确定模块1503,用于确定期望幅频响应序列;
目标幅频响应序列计算模块1504,用于依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列;
滤波系数确定模块1505,用于依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数;
信号补偿模块1506,用于采用所述滤波系数进行信号补偿,以输出目标幅频响应。
在本发明实施例中,所述幅频响应确定模块1501具体可以包括如下子模块:
信息获取子模块,用于获取CIC滤波器的降采样倍数、级联环节数量,以及,各个环节的延迟信息;
第一归一化幅频响应计算子模块,用于采用所述降采样倍数、级联环节数量,以及,各个环节的延迟信息,计算所述CIC滤波器降采样后的归一化幅频响应。
在本发明实施例中,所述幅频响应确定模块1501还可以包括如下子模块:
初始信息确定子模块,用于确定所述CIC滤波器的归一化幅频响应和降采样倍数;
第二归一化幅频响应计算子模块,用于根据所述CIC滤波器的归一化幅频响应和降采样倍数,计算所述CIC滤波器降采样后的归一化幅频响应。
在本发明实施例中,所述初始幅频响应序列生成模块1502具体可以包括如下子模块:
抽样间隔确定子模块,用于确定抽样间隔;
离散化抽样子模块,用于按照所述抽样间隔,对所述CIC滤波器降采样后的归一化幅频响应进行离散化抽样,以获得初始幅频响应序列,所述初始幅频响应序列包括多个序列点。
在本发明实施例中,所述期望幅频响应序列确定模块1503具体可以包括如下子模块:
期望幅频响应确定子模块,用于分别确定在各个序列点上的期望幅频响应;
期望幅频响应序列生成子模块,用于根据所述各个序列点上的期望幅频响应生成期望幅频响应序列。
在本发明实施例中,所述目标幅频响应序列计算模块1504具体可以包括如下子模块:
目标幅频响应序列计算子模块,用于分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,以获得CIC补偿滤波器的目标幅频响应序列。
在本发明实施例中,所述滤波系数确定模块1505具体可以包括如下子模块:
离散傅里叶逆变换子模块,用于对所述目标幅频响应序列进行离散傅里叶逆变换,以获得时域序列,以所述时域序列作为所述CIC补偿滤波器的滤波系数。
本发明实施例还公开了一种信号处理的装置,包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述计算机程序被所述处理器执行时实现如上述的信号处理的方法的步骤。
本发明实施例还公开了一种计算机可读存储介质,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现如上述的信号处理的方法的步骤。
对于装置实施例而言,由于其与方法实施例基本相似,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。
本领域内的技术人员应明白,本发明实施例的实施例可提供为方法、装置、或计算机程序产品。因此,本发明实施例可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明实施例可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明实施例是参照根据本发明实施例的方法、终端设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理终端设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理终端设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理终端设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理终端设备上,使得在计算机或其他可编程终端设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程终端设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
尽管已描述了本发明实施例的优选实施例,但本领域内的技术人员一旦得知了基本创造性概念,则可对这些实施例做出另外的变更和修改。所以,所附权利要求意欲解释为包括优选实施例以及落入本发明实施例范围的所有变更和修改。
最后,还需要说明的是,在本文中,诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来,而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。而且,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者终端设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、物品或者终端设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括所述要素的过程、方法、物品或者终端设备中还存在另外的相同要素。
以上对本发明所提供的一种信号处理的方法、一种信号处理的装置和一种存储介质,进行了详细介绍,本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种信号处理的方法,其特征在于,包括:
确定CIC滤波器降采样后的幅频响应;
依据所述降采样后的幅频响应生成初始幅频响应序列;
确定期望幅频响应序列;所述期望幅频响应序列为初始幅频响应序列中的各个序列点期望达到的幅频响应;
依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列;
依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数;
采用所述滤波系数进行信号补偿,以输出目标幅频响应;
所述依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列的步骤包括:分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,以获得所述CIC补偿滤波器的目标幅频响应序列;
其中,所述依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数的步骤包括:
对所述目标幅频响应序列进行离散傅里叶逆变换,以获得时域序列,以所述时域序列作为所述CIC补偿滤波器的滤波系数。
2.根据权利要求1所述的方法,其特征在于,所述确定CIC滤波器降采样后的幅频响应的步骤包括:
获取CIC滤波器的降采样倍数、级联环节数量,以及,各个环节的延迟信息;
采用所述降采样倍数、级联环节数量,以及,各个环节的延迟信息,计算所述CIC滤波器降采样后的归一化幅频响应。
3.根据权利要求1所述的方法,其特征在于,所述确定CIC滤波器降采样后的幅频响应的步骤包括:
确定所述CIC滤波器的归一化幅频响应和降采样倍数;
根据所述CIC滤波器的归一化幅频响应和降采样倍数,计算所述CIC滤波器降采样后的归一化幅频响应。
4.根据权利要求2或3所述的方法,其特征在于,所述依据所述降采样后的幅频响应生成初始幅频响应序列的步骤包括:
确定抽样间隔;
按照所述抽样间隔,对所述CIC滤波器降采样后的归一化幅频响应进行离散化抽样,以获得初始幅频响应序列,所述初始幅频响应序列包括多个序列点。
5.根据权利要求4所述的方法,其特征在于,所述确定期望幅频响应序列的步骤包括:
分别确定在各个序列点上的期望幅频响应;
根据所述各个序列点上的期望幅频响应生成期望幅频响应序列。
6.一种信号处理的装置,其特征在于,包括:
幅频响应确定模块,用于确定CIC滤波器降采样后的幅频响应;
初始幅频响应序列生成模块,用于依据所述降采样后的幅频响应生成初始幅频响应序列;
期望幅频响应序列确定模块,用于确定期望幅频响应序列;所述期望幅频响应序列为初始幅频响应序列中的各个序列点期望达到的幅频响应;
目标幅频响应序列计算模块,用于依据所述初始幅频响应序列和期望幅频响应序列,计算CIC补偿滤波器的目标幅频响应序列;所述目标幅频响应序列计算模块包括:计算模块,用于分别计算各个序列点上对应的期望幅频响应序列数值与初始幅频响应序列数值之间的商值,以获得所述CIC补偿滤波器的目标幅频响应序列;
滤波系数确定模块,用于依据所述目标幅频响应序列,确定所述CIC补偿滤波器的滤波系数;
信号补偿模块,用于采用所述滤波系数进行信号补偿,以输出目标幅频响应;
所述滤波系数确定模块包括:
离散傅里叶逆变换子模块,用于对所述目标幅频响应序列进行离散傅里叶逆变换,以获得时域序列,以所述时域序列作为所述CIC补偿滤波器的滤波系数。
7.一种信号处理的装置,其特征在于,包括处理器、存储器及存储在所述存储器上并可在所述处理器上运行的计算机程序,所述计算机程序被所述处理器执行时实现如权利要求1至5中任一项所述的信号处理的方法的步骤。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质上存储计算机程序,所述计算机程序被处理器执行时实现如权利要求1至5中任一项所述的信号处理的方法的步骤。
CN201711420773.5A 2017-12-25 2017-12-25 一种信号处理的方法和装置 Active CN108011615B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711420773.5A CN108011615B (zh) 2017-12-25 2017-12-25 一种信号处理的方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711420773.5A CN108011615B (zh) 2017-12-25 2017-12-25 一种信号处理的方法和装置

Publications (2)

Publication Number Publication Date
CN108011615A CN108011615A (zh) 2018-05-08
CN108011615B true CN108011615B (zh) 2020-06-02

Family

ID=62061083

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711420773.5A Active CN108011615B (zh) 2017-12-25 2017-12-25 一种信号处理的方法和装置

Country Status (1)

Country Link
CN (1) CN108011615B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111258264B (zh) * 2020-02-24 2021-06-15 北京龙鼎源科技股份有限公司 现场噪声的滤波方法及装置、存储介质和处理器
CN112769410B (zh) * 2020-12-25 2024-06-11 西安讯飞超脑信息科技有限公司 滤波器构建方法、音频处理方法及电子设备、存储装置
CN113325215B (zh) * 2021-04-28 2022-04-26 哈尔滨工业大学 一种基于多级梳状滤波的光纤电流互感器输出控制器的设计方法
CN116318463A (zh) * 2023-03-27 2023-06-23 苏州联讯仪器股份有限公司 一种矫正失真信道响应的方法、系统及装置

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102055435A (zh) * 2010-12-23 2011-05-11 中国科学院武汉物理与数学研究所 一种窄带数字滤波器
CN104836547A (zh) * 2015-06-05 2015-08-12 中国科学院武汉物理与数学研究所 一种短群延时数字滤波方法
CN105066990A (zh) * 2015-07-07 2015-11-18 哈尔滨工程大学 一种适用于捷联惯性导航的高精度数字滤波器
CN106452388A (zh) * 2016-09-29 2017-02-22 电子科技大学 一种基于并行运算的cic滤波器设计方法
CN106533394A (zh) * 2016-11-11 2017-03-22 江西师范大学 一种基于自适应滤波器幅频响应的高精度频率估计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101523738B1 (ko) * 2010-12-23 2015-06-04 한국전자통신연구원 다중 대역폭을 지원하는 디지털 프론트 엔드를 갖는 신호 수신장치 및 이를 이용한 신호 수신방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102055435A (zh) * 2010-12-23 2011-05-11 中国科学院武汉物理与数学研究所 一种窄带数字滤波器
CN104836547A (zh) * 2015-06-05 2015-08-12 中国科学院武汉物理与数学研究所 一种短群延时数字滤波方法
CN105066990A (zh) * 2015-07-07 2015-11-18 哈尔滨工程大学 一种适用于捷联惯性导航的高精度数字滤波器
CN106452388A (zh) * 2016-09-29 2017-02-22 电子科技大学 一种基于并行运算的cic滤波器设计方法
CN106533394A (zh) * 2016-11-11 2017-03-22 江西师范大学 一种基于自适应滤波器幅频响应的高精度频率估计方法

Also Published As

Publication number Publication date
CN108011615A (zh) 2018-05-08

Similar Documents

Publication Publication Date Title
CN108011615B (zh) 一种信号处理的方法和装置
Novak et al. Synchronized swept-sine: Theory, application, and implementation
JP6573869B2 (ja) 仮想サンプルレートを増大させた音声フィルタリング
KR101687658B1 (ko) 처프-지 역변환 방법 및 시스템
Kopparapu et al. Identifying optimal Gaussian filter for Gaussian noise removal
RU2015140752A (ru) Модуль определения коэффициентов частотной корректирующей фильтрации, устройство, процессор коэффицентов частотной корректирующей фильтрации, система и способы
RU2008134898A (ru) Эффективная фильтрация банком комплексно-модулированных фильтров
DE102017203804A1 (de) Umwandlung digitaler Abtastrate
CN109741762B (zh) 声音活动检测方法及装置和计算机可读存储介质
JPS6051017A (ja) アナログ信号の解析及び検索方法及び装置
CN107785026B (zh) 一种用于机顶盒室内回声消除的时延估计方法
Pei et al. Design of variable comb filter using FIR variable fractional delay element
US9036752B2 (en) Low-delay filtering
Okoniewski et al. A concept of IIR filters with time-varying coefficients and equalised group delay response
EP2495874A2 (en) Apparatus and method for adaptive signal processing
Ferreira et al. Enhancement of the comb filtering selectivity using iterative moving average for periodic waveform and harmonic elimination
CN110896303B (zh) 滤波方法和滤波装置以及存储介质
CN109951174B (zh) 一种快速傅里叶变换滤波方法
CN105099396B (zh) 一种滤波器切换方法、装置以及医疗设备
CN114047499A (zh) 传感器阵列信号幅度不一致时的信号归一化方法及设备
JP6154777B2 (ja) 高速畳込近似装置、高速畳込近似方法、プログラム
JP2022045086A (ja) 残響を求めるためのシステム
Ahmed et al. Resource efficient pre-processor for drift removal in neurochemical signals
Borkowski et al. Multi-Rate Signal Processing Issues in Active Safety Algorithms
Stanciu et al. B-Spline Functions for the Design of Crossover Digital Filters

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant