CN114448390A - 一种Biquad数字滤波器装置及实现方法 - Google Patents

一种Biquad数字滤波器装置及实现方法 Download PDF

Info

Publication number
CN114448390A
CN114448390A CN202210340180.2A CN202210340180A CN114448390A CN 114448390 A CN114448390 A CN 114448390A CN 202210340180 A CN202210340180 A CN 202210340180A CN 114448390 A CN114448390 A CN 114448390A
Authority
CN
China
Prior art keywords
biquad
data
module
digital filter
signal data
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.)
Pending
Application number
CN202210340180.2A
Other languages
English (en)
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.)
Zhejiang Xinsheng Electronic Technology Co Ltd
Original Assignee
Zhejiang Xinsheng Electronic Technology 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 Zhejiang Xinsheng Electronic Technology Co Ltd filed Critical Zhejiang Xinsheng Electronic Technology Co Ltd
Priority to CN202210340180.2A priority Critical patent/CN114448390A/zh
Publication of CN114448390A publication Critical patent/CN114448390A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0202Two or more dimensional filters; Filters for complex signals
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0283Filters characterised by the filter structure
    • H03H17/0292Time multiplexed filters; Time sharing filters
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H2017/0072Theoretical filter design
    • H03H2017/009Theoretical filter design of IIR filters

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Complex Calculations (AREA)

Abstract

本发明提供一种Biquad数字滤波器装置及实现方法,该Biquad数字滤波器装置包括数据输入存储模块、数据选择模块、Biquad运算模块和控制模块,数据输入存储模块用于接收和存储信号数据和滤波器系数;数据选择模块用于选择参与运算的信号数据和滤波器系数;Biquad运算模块用于进行滤波运算并将最终的滤波结果进行输出;Biquad运算模块包括一个乘法器和一个加法器;控制模块用于控制数据输入存储模块、数据选择模块和Biquad运算模块。本发明中的Biquad数字滤波器装置通过时分复用同一个乘法器减少了乘法器资源,且其控制逻辑简单,可以用于形成多级、多通道Biquad数字滤波器,该Biqaud数字滤波器装置更接近于DSP的应用特性,可以实现灵活的调用。

Description

一种Biquad数字滤波器装置及实现方法
技术领域
本发明涉及数字信号处理领域,尤其涉及一种Biquad数字滤波器装置及实现方法。
背景技术
双二阶滤波器(即Biquad滤波器)是一种常见的滤波器,被广泛地应用于数字音频信号处理,以及其他工程控制领域,如飞行控制技术中的陀螺仪滤波,以及手机光学防抖等。
Biquad数字滤波器是一种特殊的二阶IIR滤波器(具有两个极点和两个零点),Biquad数字滤波器有多种形式,其中最为常用的是二阶差分方程的方式:
Figure DEST_PATH_IMAGE002
,也被称之为直接I型。在上述二阶差分方程中,b0、b1、b2、a1和a2为IIR滤波器的系数,这些系数是提前由算法仿真确定;
Figure DEST_PATH_IMAGE004
Figure DEST_PATH_IMAGE006
Figure DEST_PATH_IMAGE008
为输入信号,其中,
Figure 955273DEST_PATH_IMAGE004
为当前时刻的输入信号,
Figure 934730DEST_PATH_IMAGE006
为前一个时刻的输入信号,
Figure 776784DEST_PATH_IMAGE008
为前两个时刻的输入信号;
Figure DEST_PATH_IMAGE010
Figure DEST_PATH_IMAGE012
Figure DEST_PATH_IMAGE014
为输出信号,
Figure 732846DEST_PATH_IMAGE010
为当前时刻的输出信号,
Figure 609535DEST_PATH_IMAGE012
为前一个时刻的输出信号,
Figure 431997DEST_PATH_IMAGE014
为前两个时刻的输出信号。
如图1中所示,为一个典型的Biquad数字滤波器示意图,其表示直接I型的逻辑运算结构,该结构包含五个乘法器、四个延迟单元和一个加法器,该结构具有单个求和点的特性,因此,可易于实现在定点计算时的中间过程值的溢出控制。
此外,若需要设计一个截止频率更高的滤波器,可以采用两个或者多个Biquad滤波器级联的方式,相比于单个多阶的IIR滤波器,采用了级联方式的Biquad滤波器结构将可以显著地降低滤波器对抽头系数的敏感性。例如,一个双二阶低通滤波器有12dB/octave斜度的截止频点,可以用于音调控制;如果需要设计一个24dB/octave的斜度,可以级联两个双二阶滤波器,这个结构比一个单4阶滤波器对参数更不敏感。因此,Biquad滤波器还常被用作为更复杂滤波器的基本构建单元,如图2中所示,为一个四级Biquad滤波器级联的滤波器结构。从图2中的四级Biquad滤波器级联的滤波器结构可以看出,这样的结构中将存在大量的乘法器、延时单元和加法器,不利于硬件结构的小型化和集成化。目前,也有多级Biquad滤波器复用一个Biquad单元块的结构,如图3中所示,为两级单通道Biquad运算结构,采用时分复用的方式,增加多路选择控制逻辑,选择当前通路的信号源、中间运算缓存以及滤波器的系数,图3中同时给出了相对于系统时钟CLOCK,多路选择器输入信号Q的状态信号时序图,多路选择器输入信号Q的状态在系统时钟的每个周期之间在0和1之间交替。在图3中所示的结构中,包含有多个多路选择器及其控制逻辑,占用了较多的运算资源,且这样的结构中,其时序是固定的,仅适用于两级Biquad滤波器的情况,并不适用于其他多级结构。
在很多应用场景中,通常需要多通道的并行滤波处理,例如,音频滤波时需要处理左右声道(L/R-channel),陀螺仪(Gyroscope)需要处理X、Y以及Z轴的信号。各通道的输入信号各不相同,各通道的滤波器系数也可能存在差异。以双通道滤波为例:通常的实现方式有两种,一种是采用并行的方式,即单独使用两个传统结构的Biquad滤波器,如图4中所示,为典型的并行Biquad双通道结构,这样并行的滤波器结构中,每增加一个通道,将需要相应地增加一个Biquad滤波器,因此,这样的结构会随着通道数的增加,其硬件元件数量急剧增加,不利于硬件结构的小型化和集成化。另一种方式是采用时分复用的方式,增加多路选择控制逻辑,选择当前通路的信号源、中间运算缓存以及滤波器的滤波系数,如图5中所示,为一种常见的双通道复用Biquad结构,多路选择器根据时钟顺序向乘法器提供两个不同系数中的其中之一,多路选择器具有两种功能,一是输入信号的数据选择,二是输出信号的数据保持。图5中同时给出了相对于系统时钟CLOCK,多路选择器输入信号Q的状态信号时序图,多路选择器输入信号Q的状态在系统时钟的每个周期之间在0和1之间交替。在图5中所示的结构中,包含有多个多路选择器及其控制逻辑,占用了较多的运算资源,且这样的结构中,其时序是固定的,图5中所示的结构仅适用于x0通道和x1通道执行交替运算的形式,对于其他的顺序并不适用。
通过上述介绍可知,目前的多级Biquad滤波器以及多通道Biquad滤波器存在以下问题:(1)采用多个Biquad单元进行级联或者并行工作,这样的结构中,每增加一级或者每增加一个通道将会增加相应的硬件和软件资源,所消耗的资源较多,不利于硬件结构的小型化和集成化;(2)采用分时复用结构的多级Biquad滤波器或者多通道Biquad滤波器,其中包含多个多路选择器以及多路选择器的控制逻辑,占用了较多的运算资源;此外,这样的结构只适用于特定级数的滤波器以及特定通道数量的滤波器,灵活性较差。
因此,亟需要设计一种Biquad滤波器,其能减小对硬件资源的消耗,且具有较好的灵活性,能适用于不同级数以及不同通道数量的滤波器。
发明内容
为了解决上述技术问题,本发明提供一种Biquad数字滤波器装置,其特征在于,所述Biquad数字滤波器装置包括数据输入存储模块、数据选择模块、Biquad运算模块和控制模块,其中:所述数据输入存储模块用于接收和存储信号数据和滤波器系数,所述信号数据包括当前的输入信号数据、之前的输入信号数据以及之前Biquad数字滤波器装置的输出信号数据;所述数据选择模块用于从数据输入存储模块中选择参与运算的信号数据和滤波器系数,并将该信号数据和滤波器系数传输到Biquad运算模块;所述Biquad运算模块用于进行滤波运算并将最终的滤波结果进行输出;所述Biquad运算模块包括一个乘法器和一个加法器,滤波器系数与信号数据的乘法运算时分复用同一个乘法器,并通过加法器将乘法计算结果进行累加;所述控制模块用于控制所述数据输入存储模块、所述数据选择模块和所述Biquad运算模块;所述控制模块控制所述数据输入存储模块进行信号数据和滤波器系数的读取与存储;所述控制模块控制所述数据选择模块选择每次参与运算的信号数据和滤波器系数;所述控制模块控制Biquad运算模块进行乘累加运算,并将滤波结果输出。
采用本发明中提供的Biquad数字滤波器装置,具有如下优势:(1)该Biquad数字滤波器装置结构较为简单,其中的Biquad运算模块只包括一个乘法器和一个加法器,通过时分复用同一个乘法器,即可实现Biquad滤波运算的功能,减少了乘法器的资源;(2)可以通过该Biqaud数字滤波器装置实现多级Biquad级联滤波器结构、实现多通道Biquad滤波器结构甚至实现多通道多级级联的Biquad数字滤波器结构,在实现多级Biquad级联滤波器结构、实现多通道Biquad滤波器结构甚至实现多通道多级级联的Biquad数字滤波器结构时,只需要增加数据输入存储模块中的滤波器系数存储通路和信号数据的配置通路,不再需要增加额外的计算单元,硬件成本较低;(3)控制逻辑简单,实现多级Biquad级联滤波器结构时,其控制逻辑与单级Biquad滤波器基本相同;当实现多通道Biquad滤波器结构时,通过分时复用同一个Biquad数字滤波器装置,只是改变其中的参数,控制逻辑完全相同。因此,本发明中的Biquad数字滤波器装置更接近于数字信号处理(DSP)的应用特性,其将Biquad滤波模块例化成函数体的形式,可以实现自由灵活的调用实现不同的设计目的;(4)本发明中提供的Biquad数字滤波器装置尤其适用于音频和传感器等低速率数据信号的处理,因为系统的运行频率将远高于一般音频和传感器等信号输入频率,采用本发明中提供的Biquad数字滤波器装置,可以将对输入信号的乘累加运算分拆到几个系统时钟内完成,其对信号的处理速率与先前技术中提供的运用多个乘法器的结构相当。
优选地,所述控制模块在控制所述数据输入存储模块时,在每次滤波之前,控制所述数据输入存储模块获取信号数据和滤波器系数,在单次的Biquad滤波运算完成时,所述控制模块控制所述数据输入存储模块更新信号数据。
优选地,所述数据选择模块包括若干个多路选择器,用于分别对信号数据和滤波器系数进行选择。
优选地,所述数据输入存储模块从外部存储器获取滤波器系数和缓存的信号数据,所述缓存的信号数据包括之前的输入信号数据和之前Biquad数字滤波器装置的输出信号数据,所述数据输入存储模块通过输入接口获取当前的输入信号数据。
优选地,所述数据输入存储模块将滤波器的计算结果存储至外部存储器。
优选地,所述数据输入存储模块与外部存储器之间的数据传输由数据控制接口完成。
优选地,所述控制模块通过系统时钟和计时器对数据输入存储模块、数据选择模块和Biquad运算模块进行控制。
优选地,所述控制模块控制所述数据选择模块在每个时钟周期内,选取一组对应的滤波器系数和信号数据。
优选地,所述控制模块控制所述Biquad运算模块在每个时钟周期内对所述数据选择模块选取的滤波器系数和信号数据进行乘法运算,并在下一个时钟周期内将该乘法运算的结果与之前加法器的计算结果进行累加。
优选地,所述Biquad数字滤波器装置用于构成多级Biquad数字滤波器、多通道Biquad数字滤波器以及多级多通道数字滤波器。
优选地,所述Biquad数字滤波器装置用于处理以下信号其中之一:音频信号、传感器信号。
本发明的另一方面,提供一种Biquad数字滤波器的实现方法,其特征在于,所述Biquad数字滤波器的实现方法包含以下步骤:步骤S01:数据输入存储模块获取或更新滤波器系数以及信号数据,所述信号数据包括当前的输入信号数据、之前的输入信号数据以及之前Biquad数字滤波器装置的输出信号数据;步骤S02:数据选择模块选择参与运算的滤波器系数和信号数据;步骤S03:Biquad运算模块进行乘法运算和累加运算;步骤S04:判断该单次滤波是否已经完成,若没有完成,则返回至步骤S02中;若已经完成,则执行步骤S05;步骤S05:判断所有的运算是否已经完成,若没有完成,则返回至步骤S01中更新数据输入存储模块中的信号数据;若已经完成,则执行步骤S06;步骤S06:更新数据输入存储模块中的信号数据,并输出滤波结果。
采用本发明中提供的Biquad数字滤波器的实现方法,具有如下优势:(1)该方法具有很好的通用性,该实现方法可以用于实现多级级联的Biquad数字滤波器结构、多通道的Biquad数字滤波器结构甚至是多级级联的多通道Biquad数字滤波器结构,采用该实现方法的Biquad数字滤波器可以将Biquad滤波模块例化成函数体,可以实现自由灵活的调用;(2)采用该实现方法的Biquad数字滤波器硬件资源占用较少,其中的Biquad运算模块只包含一个乘法器和一个加法器,通过时分复用的运算结构有效地减少了乘法器资源,减少了硬件资源的占用。
优选地,控制模块执行所述步骤S04和所述步骤S05中的判断操作。
附图说明
图1为一个典型的Biquad数字滤波器示意图。
图2为一个四级Biquad滤波器级联的滤波器结构。
图3为两级单通道Biquad运算结构。
图4为典型的并行Biquad双通道结构。
图5为一种常见的双通道复用Biquad结构。
图6为本发明中提供的一种Biquad数字滤波器装置示意图。
图7为本发明中单级Biquad数字滤波器装置的运算流程图。
图8为本发明中的Biquad数字滤波器装置实现2级Biquad级联结构的运算控制示意图。
图9为本发明中2级Biquad级联结构的运算硬件电路图。
图10为双通道、2级级联的Biquad滤波器结构的数据存储通路示意图。
图11为本发明中一种Biquad数字滤波器实现方法的流程图。
具体实施方式
以下配合图式及本发明的较佳实施例,进一步阐述本发明为达成预定发明目的所采取的技术手段。
如图6中所示,为本发明中提供的一种Biquad数字滤波器装置,该装置包括数据输入存储模块、数据选择模块、Biquad运算模块和控制模块。其中,数据输入存储模块用于接收和存储信号数据和滤波器系数,并将信号数据和滤波器系数通过数据选择模块输入至Biquad运算模块。数据输入存储模块包括若干个输入接口和数据存储单元,若干个输入接口用于接收不同的信号数据和滤波器系数,并将这些输入信号和滤波器系数存储至数据存储单元。数据选择模块用于从数据输入存储单元中选择参与运算的信号数据和滤波器系数,并将其传输到Biquad运算模块,数据选择模块包括若干个多路选择器,该若干个多路选择器与数据存储单元和控制模块相连接。Biquad运算模块用于进行滤波运算并将最终的滤波结果进行输出。在本实施例中,Biquad运算模块包括一个加法器和一个乘法器,滤波器系数与信号数据的乘法运算时分复用同一个乘法器,在每一个时钟周期内,乘法器将输入信号和滤波器系数进行乘法运算,并在下一个时钟周期内将该乘法计算结果与之前的加法器的计算结果进行累加,将新的累加结果进行存储,当所有的结果计算完成时,Biquad运算模块将累加结果进行输出。控制模块用于控制数据输入存储模块、数据选择模块和Biquad运算模块。具体而言,控制模块根据Biquad滤波器的结构,控制数据输入存储模块进行输入信号和滤波器系数的输入和存储,控制模块控制数据存储单元在单次Biquad滤波运算完成时,开始进行信号数据的更新,此处所指的单次Biquad滤波运算表示进行了一次单级Biquad滤波运算;控制模块控制数据选择模块在每个时钟周期选择计算所需的滤波器系数和信号数据;控制模块用于控制Biquad运算模块在每一个时钟周期内将多路选择器选择的滤波器系数和信号数据进行乘法运算,在下一个时钟周期内将该乘法计算结果与之前加法器的计算结果进行累加,并在所有的结果计算完成时,控制Biquad运算模块将滤波结果进行输出。
如图7中所示,图7为本发明中提供的单级Biquad数字滤波器装置的运算控制示意图。以一个单级Biquad数字滤波器为例,说明单级Biquad数字滤波器装置的运算过程,其中,滤波器系数与信号数据的对应关系与所有采用直接I型的Biquad结构相同。首先,在开始进行计算之前,数据输入存储模块存储计算所需的参数,这些参数包括滤波器系数b0、b1、b2、a1和a2,以及信号数据X0、X0_D1、X0_D2、BIQUAD_OUT_D1和BIQUAD_OUT_D2,其中,X0表示输入信号,X0_D1表示上一次的输入信号,X0_D2表示上两次的输入信号,BIQUAD_OUT_D1表示上一次的Biquad输出信号,BIQUAD_OUT_D2表示上两次的Biquad输出信号。
在开始计算时,控制模块通过系统时钟和计时器对数据输入存储模块、数据选择模块和Biquad运算模块进行控制,计时器用于在Biquad数字滤波器装置开始工作后,对系统时钟进行计数。数据选择模块在每个时钟周期选择计算所需的参数输入到乘法器进行乘运算,乘运算在每个时钟周期计算所需的参数包括滤波器系数和信号数据,滤波器系数和信号数据两者一一对应,其对应关系如下:滤波器系数b0对应输入信号X0,滤波器系数b1对应上一次的输入信号X0_D1,滤波器系数b2对应上两次的输入信号X0_D2,滤波器系数a1对应上一次的Biquad输出信号BIQUAD_OUT_D1以及滤波器系数a2对应上两次的Biquad输出信号BIQUAD_OUT_D2。乘运算的结果进行自累加,即在每一个时钟周期内将数据选择模块选择的滤波器参数和信号数据进行乘法运算,并在下一个时钟周期内将该乘法运算结果与之前加法器的计算结果进行累加,该自累加运算受控于控制模块中的计时器calc_cnt。具体的计算过程如下:在第一个时钟周期内,即计时器calc_cnt=1时,滤波器系数b0与对应输入信号X0进行乘运算,得到乘运算结果MUL_RESULT。在第二个时钟周期内,即计时器calc_cnt=2时,将该乘运算结果MUL_RESULT与加法器的计算结果MUL_SUM进行累加,并将结果保存在MUL_SUM中,其中加法器的计算结果MUL_SUM的初始值为0;并同时进行第二组乘运算,即滤波器系数b1与对应的上一次输入信号X0_D1进行乘运算,得到乘运算结果MUL_RESULT。在第三个时钟周期内,即计时器calc_cnt=3时,将该乘运算结果MUL_RESULT与加法器的计算结果MUL_SUM进行累加,并将结果保存在MUL_SUM中;并同时进行第三组乘运算,即滤波器系数b2与对应的上两次输入信号X0_D2进行乘运算,得到乘运算结果MUL_RESULT。在第四个时钟周期内,即计时器calc_cnt=4时,将该乘运算结果MUL_RESULT与加法器的计算结果MUL_SUM进行累加,并将结果保存在MUL_SUM中;并同时进行第四组乘运算,即滤波器系数a1与对应的上一次Biquad输出信号BIQUAD_OUT_D1进行乘运算,得到乘运算结果MUL_RESULT。在第五个时钟周期内,即计时器calc_cnt=5时,将该乘运算结果MUL_RESULT与加法器的计算结果MUL_SUM进行累加,并将结果保存在MUL_SUM中;并同时进行第五组乘运算,即滤波器系数a2与对应的上两次Biquad输出信号BIQUAD_OUT_D2进行乘运算,得到乘运算结果MUL_RESULT。在第六个时钟周期内,即计时器calc_cnt=6时,将该乘运算结果MUL_RESULT与加法器的计算结果MUL_SUM进行累加,并将结果保存在MUL_SUM中,此时,MUL_SUM的值即为该单级Biquad数字滤波器的输出值,控制模块控制Biquad运算模块将该输出值进行输出,并控制数据输入存储模块中的值进行更新,该运算周期结束。
采用本发明中提供的Biquad数字滤波器装置具有如下优势:(1)本发明中提供的Biquad数字滤波器装置结构较为简单,其中的Biquad运算模块只包括一个乘法器和一个加法器,通过时分复用同一个乘法器,即可实现Biquad滤波运算的功能,减少了乘法器的资源;(2)本发明中提供的Biquad数字滤波器装置尤其适用于音频和传感器等信号的处理,因为系统的运行频率将远高于一般音频和传感器等信号输入频率,采用本发明中提供的Biquad数字滤波器装置,可以将对输入信号的乘累加运算分拆到几个系统时钟内完成,其对信号的处理速率与先前技术中提供的运用多个乘法器的结构相当。
本发明中提供的Biquad数字滤波器装置不仅适用于单级Biquad数字滤波器,更可以使用该Biquad数字滤波器装置实现多级Biquad滤波器和多通道Biquad滤波器。
如图8中所示,为采用本发明中的Biquad数字滤波器装置实现2级Biquad级联结构的运算控制示意图,其中,系数与数据的对应关系与所有采用直接I型的Biquad结构相同。首先,在开始进行计算之前,数据输入存储模块存储计算所需的参数,这些参数包括第一级滤波器系数b0,0、b1,0、b2,0、a1,0和a2,0,第二级滤波器系数b3,0、b4,0、b5,0、a3,0和a4,0以及信号数据X0、X0_D1、X0_D2、BIQUAD_1ST _OUT_D1、BIQUAD_1ST_OUT_D2, BIQUAD_1ST_OUT,BIQUAD_2ND_OUT_D1和BIQUAD_2ND_OUT_D2,其中,X0表示输入信号,X0_D1表示上一次的输入信号,X0_D2表示上两次的输入信号,BIQUAD_1ST_OUT_D1表示上一次的第一级Biquad输出信号,BIQUAD_1ST_OUT_D2表示上两次的第一级Biquad输出信号,BIQUAD_1ST_OUT表示当前第一级Biquad输出信号,BIQUAD_2ND_OUT_D1表示上一次的第二级Biquad输出信号,以及BIQUAD_2ND_OUT_D2表示上两次的第二级Biquad输出信号。
在开始计算时,数据选择模块在每个时钟周期选择计算所需的参数输入到乘法器进行乘运算,乘运算在每个时钟周期计算所需的参数包括滤波器系数和信号数据,滤波器系数和信号数据两者一一对应,其对应关系如下:滤波器系数b0,0对应输入信号X0;滤波器系数b1,0对应上一次的输入信号X0_D1;滤波器系数b2,0对应上两次的输入信号X0_D2;滤波器系数a1,0对应上一次的第一级Biquad输出信号BIQUAD_1ST_OUT_D1;滤波器系数a2,0对应上两次的第一级Biquad输出信号BIQUAD_1ST_OUT_D2;滤波器系数b3,0对应当前第一级Biquad输出信号BIQUAD_1ST_OUT;滤波器系数b4,0对应上一次的第一级Biquad输出信号BIQUAD_1ST_OUT_D1;滤波器系数b5,0对应上两次的第一级Biquad输出信号BIQUAD_1ST_OUT_D2;滤波器系数a3,0对应上一次的第二级Biquad输出信号BIQUAD_2ND_OUT_D1以及滤波器系数a4,0对应上两次的第二级Biquad输出信号BIQUAD_2ND_OUT_D2。乘运算的结果进行自累加,即在每一个时钟周期内将数据选择模块选择的滤波器参数和信号数据进行乘法运算,并在下一个时钟周期内将该乘法运算结果与之前加法器的计算结果进行累加,该自累加运算受控于控制模块中的计时器calc_cnt。具体的计算过程如下:在第一个时钟周期内,即计时器calc_cnt=1时,滤波器系数b0,0与对应输入信号X0进行乘运算,得到乘运算结果MUL_RESULT。在第二个时钟周期内,即计时器calc_cnt=2时,将该乘运算结果MUL_RESULT与加法器的计算结果MUL_SUM进行累加,并将结果保存在MUL_SUM中,其中加法器的计算结果MUL_SUM的初始值为0。此后的9组乘法运算结果MUL_RESULT都累加到MUL_SUM上,即在每个时钟周期内,均进行计算MUL_SUM=MUL_SUM+MUL_RESULT。在第五个时钟周期内,即计时器calc_cnt=5时,得到滤波器系数a2,0与上两次的第一级Biquad输出信号BIQUAD_1ST_OUT_D2的相乘结果MUL_RESULT,并在下一个时钟周期内,即计时器calc_cnt=6时,叠加到MUL_SUM,此时MUL_SUM的值即为第一级Biquad滤波器的输出值。该输出值一方面是作为第二级滤波器的输入,与滤波器系数b3,0相乘,另一方面用于刷新数据输入存储模块中BIQUAD_1ST_OUT_D1的值。同理,在第十一个时钟周期内,即计时器calc_cnt=11时,当系数a4,0与上两次的第二级Biquad输出信号BIQUAD_2ND_OUT_D2的相乘结果MUL_RESULT叠加到MUL_SUM时,此时MUL_SUM的值即为第二级Biquad滤波器的输出值,该值也用于刷新数据存储模块中BIQUAD_2ND_OUT_D1的值。
图9为采用本发明中2级Biquad级联结构的运算硬件电路图。如图9中所示,在该运算硬件电路图中,数据输入存储模块分为两部分,一部分用于存储滤波器系数,另一部分用于存储信号数据,因为信号数据中包含了一些滤波器运算的中间结果,因此,用于存储信号数据的硬件电路也在图9中表示出来。数据选择模块包括两个多路选择器,两个多路选择器分别用于选择滤波器系数和信号数据,在每个时钟周期内,多路选择器从数据存储单元中分别选择信号数据和与之对应的滤波器系数,并将这些数据传输到Biquad运算模块。Biquad运算模块中包括一个乘法器,一个加法器以及一个寄存器,乘法器用于将信号数据和与之对应的滤波器系数进行乘法运算,加法器和寄存器用于实现对上述乘法器的运算结果实现累加。在第六个时钟周期内,此时MUL_SUM的值即为第一级Biquad滤波器的输出值。该输出值一方面是作为第二级滤波器的输入,与滤波器系数b3,0相乘,另一方面用于刷新数据输入存储模块中BIQUAD_1ST_OUT_D1的值。同理,在第十一个时钟周期内,即计时器calc_cnt=11时,当系数a4,0与上两次的第二级Biquad输出信号BIQUAD_2ND_OUT_D2的相乘结果MUL_RESULT叠加到MUL_SUM时,此时MUL_SUM的值即为第二级Biquad滤波器的输出值,该值也用于刷新数据存储模块中BIQUAD_2ND_OUT_D1的值。
从图9中2级Biquad级联结构的运算硬件电路图可以看出,本发明中的Biquad数字滤波器装置在实现多级Biquad滤波器时,其只需要通过时分复用同一个乘法器即可以实现多级Biquad滤波器的功能。本发明中提供的Biquad数字滤波器装置不仅适用于2级Biquad级联滤波器结构,也适用于多级Biquad级联滤波器结构,只要相应地增加数据输入存储模块中的数据,并在得到每一级滤波器结果的输出时,刷新数据输入存储模块中的信号数据即可。
因此,本发明中提供的Biquad数字滤波器装置可以在控制逻辑相同的条件下,增加数据输入存储模块中的滤波器系数存储通路和信号数据的配置通路,即可以实现多级Biquad级联滤波器的功能。采用本发明中提供的Biquad数字滤波器装置实现多级Biquad级联滤波器结构时,具有如下优点:(1)硬件成本较小,在实现多级Biquad级联滤波器结构时,只需要增加数据输入存储模块中的滤波器系数存储通路和信号数据的配置通路,不再需要增加额外的计算单元;(2)控制逻辑简单,实现多级Biquad级联滤波器结构时,其控制逻辑与单级Biquad滤波器基本相同,因此,本发明中的Biquad数字滤波器装置更接近于数字信号处理(DSP)的应用特性,其将Biquad滤波模块例化成函数体的形式,可以实现自由灵活的调用。
进一步地,本发明中提供的Biquad数字滤波器装置通过增加数据输入存储模块中的滤波器系数存储通路和信号数据的配置通路,即可以灵活地实现Biquad滤波器的多通道、多级级联的运算结构。如图10中所示,为双通道、2级级联的Biquad滤波器结构的数据存储通路示意图。在该双通道、2级级联的Biquad滤波器结构的数据存储通路示意图中,Biquad数字滤波器装置中的数据输入存储模块分别与输入信号X0和X1以及外部存储器相连接,外部存储器用于存储滤波器系数以及缓存的信号数据,缓存的信号数据包括之前的输入信号以及之前滤波器的输出信号。该双通道、2级级联的Biquad滤波器的运算过程如下:(1)将输入信号数据X0更新到Biquad数字滤波器装置中的输入信号X;(2)Biquad数字滤波器装置中的数据输入存储模块向外部存储器获取通路0的滤波器系数和信号数据,通路0的滤波器系数包括:第一级滤波器系数B0,0、B1,0、B2,0、A1,0和A2,0;第二级滤波器系数B3,0、B4,0、B5,0、A3,0和A4,0;通路0的信号数据上一次的输入信号X0_D1、上两次的输入信号X0_D2、上一次通路0的第一级Biquad输出信号X0_BIQUAD_1ST_OUT_D1、上两次通路0的第一级Biquad输出信号X0_BIQUAD_1ST_OUT_D2、上一次通路0的第二级Biquad输出信号X0_BIQUAD_2ND_OUT_D1以及上两次通路0的第二级Biquad输出信号X0_BIQUAD_2ND_OUT_D2;(3)开始进行通路0中的时分复用累加运算,其运算过程如图8中所示,在此不进行赘述;(4)在通路0计算完成时,Biquad数字滤波器装置中产生的新的缓存数据保存至外部存储器中,Biquad数字滤波器装置中的数据输入存储模块对外部存储器中的以下数据进行更新:Biquad数字滤波器装置中的第一级Biquad的输出信号更新到外部存储器中的X0_BIQUAD_1ST_OUT_D1;Biquad数字滤波器装置中的X0_BIQUAD_1ST_OUT_D1更新到外部存储器中的X0_BIQUAD_1ST_OUT_D2;Biquad数字滤波器装置中的第二级Biquad的输出信号Y更新到外部存储器中的X0_BIQUAD_2ND_OUT_D1以及Biquad数字滤波器装置中的X0_BIQUAD_2ND_OUT_D1更新到外部存储器中的X0_BIQUAD_2ND_OUT_D2。此时,对于通道0的数据X0的2级Biquad滤波过程已经完成,在以上滤波过程完成之后,再以同样的方式执行通道1中的数据X1的滤波,其运算过程如下:(5)将输入信号数据X1更新到Biquad数字滤波器装置中的输入信号X;(2)Biquad数字滤波器装置中的数据输入存储模块向外部存储器获取通路1的滤波器系数和信号数据,通路1的滤波器系数包括:第一级滤波器系数B0,1、B1,1、B2,1、A1,1和A2,1;第二级滤波器系数B3,1、B4,1、B5,1、A3,1和A4,1;通路1的信号数据上一次的输入信号X1_D1、上两次的输入信号X1_D2、上一次通路1的第一级Biquad输出信号X1_BIQUAD_1ST_OUT_D1、上两次通路1的第一级Biquad输出信号X1_BIQUAD_1ST_OUT_D2、上一次通路1的第二级Biquad输出信号X1_BIQUAD_2ND_OUT_D1以及上两次通路1的第二级Biquad输出信号X1_BIQUAD_2ND_OUT_D2;(3)开始进行通路1中的时分复用累加运算,其运算过程如图8中所示,在此不进行赘述;(4)在通路1计算完成时,Biquad数字滤波器装置中产生的新的缓存数据保存至外部存储器中,Biquad数字滤波器装置中的数据输入存储模块对外部存储器中的以下数据进行更新:Biquad数字滤波器装置中的第一级Biquad的输出信号更新到外部存储器中的X1_BIQUAD_1ST_OUT_D1;Biquad数字滤波器装置中的X1_BIQUAD_1ST_OUT_D1更新到外部存储器中的X1_BIQUAD_1ST_OUT_D2;Biquad数字滤波器装置中的第二级Biquad的输出信号Y更新到外部存储器中的X1_BIQUAD_2ND_OUT_D1以及Biquad数字滤波器装置中的X1_BIQUAD_2ND_OUT_D1更新到外部存储器中的X1_BIQUAD_2ND_OUT_D2。此时,对于通道1中的数据X的2级Biquad滤波过程已经完成。再重复上述对于输入X0和输入X1的运算,即可以通过本发明中的Biquad数字滤波器装置实现双通道2级级联的Biquad滤波运算。
在图10中,外部存储器的存储介质可以是FLASH存储器,静态随机读写存储器(SRAM)、动态随机存取存储器(DRAM)、同步动态随机存取存储器(SDRAM)、一级缓存(L1Cache)、二级缓存(L2 Cache)、三级缓存(L3 Cache)等中的一种或者多种,在此不进行限制,只要其能达到存储目的即可。外部存储器与Biquad数字滤波器装置中的数据输入存储模块之间的数据传输由数据控制接口完成,该接口可以是AHB标准总线,也可以是自定义的传输协议。Biquad数字滤波器装置中的数据输入存储模块向外部存储器读取数据,可以是由Biquad数字滤波器装置自动发起,如Biquad数字滤波器装置中的控制模块检测到输入信号中的X已经被更新,则Biquad数字滤波器装置向外部存储器读取数据;也可以是响应软件指令。同样,Biquad数字滤波器装置中的数据输入存储模块向外部存储器存储数据,可以是由Biquad数字滤波器装置自动发起,如Biquad数字滤波器装置中的控制模块已经控制Biquad运算模块完成了所有的乘累加计算,则Biquad数字滤波器装置向外部存储器存储数据;也可以是响应软件指令。
综上所述,本发明中提供的Biquad数字滤波器装置具有如下优势:(1)该Biquad数字滤波器装置结构较为简单,其中的Biquad运算模块只包括一个乘法器和一个加法器,通过时分复用同一个乘法器,即可实现Biquad滤波运算的功能,减少了乘法器的资源;(2)在实现多级Biquad级联滤波器结构、实现多通道Biquad滤波器结构甚至实现多通道多级级联的Biquad数字滤波器结构时,只需要增加数据输入存储模块中的滤波器系数存储通路和信号数据的配置通路,不再需要增加额外的计算单元,硬件成本较低;(3)控制逻辑简单,实现多级Biquad级联滤波器结构时,其控制逻辑与单级Biquad滤波器基本相同;当实现多通道Biquad滤波器结构时,通过分时复用同一个Biquad数字滤波器装置,只是改变其中的参数,控制逻辑完全相同。因此,本发明中的Biquad数字滤波器装置更接近于数字信号处理(DSP)的应用特性,其将Biquad滤波模块例化成函数体的形式,可以实现自由灵活的调用实现不同的设计目的;(4)本发明中提供的Biquad数字滤波器装置尤其适用于音频和传感器等低速率数据信号的处理,因为系统的运行频率将远高于一般音频和传感器等信号输入频率,采用本发明中提供的Biquad数字滤波器装置,可以将对输入信号的乘累加运算分拆到几个系统时钟内完成,其对信号的处理速率与先前技术中提供的运用多个乘法器的结构相当。
本发明的另一方面,提供一种Biquad数字滤波器的实现方法,如图11中所示,图11为本发明中一种Biquad数字滤波器实现方法的流程图,该方法包含以下步骤:步骤S01:数据输入存储模块获取或更新滤波器系数以及信号数据,其中,信号数据包括输入信号数据、之前的输入信号以及之前滤波器的输出信号;步骤S02:数据选择模块选择参与运算的滤波器系数和信号数据,在每个时钟周期内,滤波器系数与信号数据一一对应,滤波器系数与信号数据的对应关系与所有采用直接I型的Biquad结构相同;步骤S03:Biquad运算模块进行乘法运算和累加运算,在该步骤中,将步骤S02中选择的滤波器系数和信号数据进行乘法运算,并在下一个时钟周期内将该乘法运算结果与之前加法器的计算结果进行累加运算;步骤S04:判断该单次滤波是否已经完成,若没有完成,则返回至步骤S02中,继续选择滤波器系数和信号数据进行运算,若已经完成,则进行步骤S05,此处所指的单次滤波表示进行了一次单级Biquad滤波运算;步骤S05:判断所有运算是否已经完成,若没有完成,则返回至步骤S01中更新滤波器系数及信号数据,若已经完成,则执行步骤S06,更新数据输入存储模块中的信号数据,并输出滤波结果。其中,控制模块执行步骤S04和步骤S05中的判断操作。
采用本发明中的Biquad数字滤波器的实现方法,具有如下优势:(1)该方法具有很好的通用性,该实现方法可以用于实现多级级联的Biquad数字滤波器结构、多通道的Biquad数字滤波器结构甚至是多级级联的多通道Biquad数字滤波器结构,采用该实现方法的Biquad数字滤波器可以将Biquad滤波模块例化成函数体,可以实现自由灵活的调用;(2)采用该实现方法的Biquad数字滤波器硬件资源占用较少,其中的Biquad运算模块只包含一个乘法器和一个加法器,通过时分复用的运算结构有效地减少了乘法器资源,减少了硬件资源的占用。
以上所述仅是本发明的优选实施例而已,并非对本发明做任何形式上的限制,虽然本发明已以优选实施例揭露如上,然而并非用以限定本发明,任何熟悉本专业的技术人员,在不脱离本发明技术方案的范围内,当可利用上述揭示的技术内容作出些许更动或修饰为等同变化的等效实施例,但凡是未脱离本发明技术方案的内容,依据本实用发明的技术实质对以上实施例所作的任何简单修改、等同变化与修饰,均仍属于本发明技术方案的范围内。

Claims (13)

1.一种Biquad数字滤波器装置,其特征在于,所述Biquad数字滤波器装置包括数据输入存储模块、数据选择模块、Biquad运算模块和控制模块,其中:
所述数据输入存储模块用于接收和存储信号数据和滤波器系数,所述信号数据包括当前的输入信号数据、之前的输入信号数据以及之前Biquad数字滤波器装置的输出信号数据;
所述数据选择模块用于从数据输入存储模块中选择参与运算的信号数据和滤波器系数,并将该信号数据和滤波器系数传输到Biquad运算模块;
所述Biquad运算模块用于进行滤波运算并将最终的滤波结果进行输出;所述Biquad运算模块包括一个乘法器和一个加法器,滤波器系数与信号数据的乘法运算时分复用同一个乘法器,并通过加法器将乘法计算结果进行累加;
所述控制模块用于控制所述数据输入存储模块、所述数据选择模块和所述Biquad运算模块;所述控制模块控制所述数据输入存储模块进行信号数据和滤波器系数的读取与存储;所述控制模块控制所述数据选择模块选择每次参与运算的信号数据和滤波器系数;所述控制模块控制Biquad运算模块进行乘累加运算,并将滤波结果输出。
2.根据权利要求1中所述的Biquad数字滤波器装置,其特征在于,所述控制模块在控制所述数据输入存储模块时,在每次滤波之前,控制所述数据输入存储模块获取信号数据和滤波器系数,在单次的Biquad滤波运算完成时,所述控制模块控制所述数据输入存储模块更新信号数据。
3.根据权利要求1中所述的Biquad数字滤波器装置,其特征在于,所述数据选择模块包括若干个多路选择器,用于分别对信号数据和滤波器系数进行选择。
4.根据权利要求1中所述的Biquad数字滤波器装置,其特征在于,所述数据输入存储模块从外部存储器获取滤波器系数和缓存的信号数据,所述缓存的信号数据包括之前的输入信号数据和之前Biquad数字滤波器装置的输出信号数据,所述数据输入存储模块通过输入接口获取当前的输入信号数据。
5.根据权利要求4中所述的Biquad数字滤波器装置,其特征在于,当滤波计算完成时,所述数据输入存储模块将滤波器的计算结果存储至外部存储器。
6.根据权利要求4中所述的Biquad数字滤波器装置,其特征在于,所述数据输入存储模块与外部存储器之间的数据传输由数据控制接口完成。
7.根据权利要求1中所述的Biquad数字滤波器装置,其特征在于,所述控制模块通过系统时钟和计时器对数据输入存储模块、数据选择模块和Biquad运算模块进行控制。
8.根据权利要求7中所述的Biquad数字滤波器装置,其特征在于,所述控制模块控制所述数据选择模块在每个时钟周期内,选取一组对应的滤波器系数和信号数据。
9.根据权利要求8中所述的Biquad数字滤波器装置,其特征在于,所述控制模块控制所述Biquad运算模块在每个时钟周期内对所述数据选择模块选取的滤波器系数和信号数据进行乘法运算,并在下一个时钟周期内将该乘法运算的结果与之前加法器的计算结果进行累加。
10.根据权利要求1-8任一项所述的Biquad数字滤波器装置,其特征在于,所述Biquad数字滤波器装置用于构成多级Biquad数字滤波器、多通道Biquad数字滤波器以及多级多通道数字滤波器。
11.根据权利要求1-8任一项所述的Biquad数字滤波器装置,其特征在于,所述Biquad数字滤波器装置用于处理以下信号其中之一:音频信号、传感器信号。
12.一种Biquad数字滤波器的实现方法,其特征在于,所述Biquad数字滤波器的实现方法包含以下步骤:
步骤S01:数据输入存储模块获取或更新滤波器系数以及信号数据,所述信号数据包括当前的输入信号数据、之前的输入信号数据以及之前Biquad数字滤波器装置的输出信号数据;
步骤S02:数据选择模块选择参与运算的滤波器系数和信号数据;
步骤S03:Biquad运算模块进行乘法运算和累加运算;
步骤S04:判断该单次滤波是否已经完成,若没有完成,则返回至步骤S02中;若已经完成,则执行步骤S05;
步骤S05:判断所有的运算是否已经完成,若没有完成,则返回至步骤S01中更新数据输入存储模块中的信号数据;若已经完成,则执行步骤S06;
步骤S06:更新数据输入存储模块中的信号数据,并输出滤波结果。
13.根据权利要求12中所述的Biquad数字滤波器的实现方法,其特征在于,控制模块执行所述步骤S04和所述步骤S05中的判断操作。
CN202210340180.2A 2022-04-02 2022-04-02 一种Biquad数字滤波器装置及实现方法 Pending CN114448390A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210340180.2A CN114448390A (zh) 2022-04-02 2022-04-02 一种Biquad数字滤波器装置及实现方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210340180.2A CN114448390A (zh) 2022-04-02 2022-04-02 一种Biquad数字滤波器装置及实现方法

Publications (1)

Publication Number Publication Date
CN114448390A true CN114448390A (zh) 2022-05-06

Family

ID=81359835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210340180.2A Pending CN114448390A (zh) 2022-04-02 2022-04-02 一种Biquad数字滤波器装置及实现方法

Country Status (1)

Country Link
CN (1) CN114448390A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115102524A (zh) * 2022-07-07 2022-09-23 武汉市聚芯微电子有限责任公司 一种滤波器

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001160736A (ja) * 1999-12-01 2001-06-12 Nec Ic Microcomput Syst Ltd デジタルフィルタ回路
US20010032226A1 (en) * 1998-01-29 2001-10-18 Koji Takano Digital filters
JP2002217684A (ja) * 2001-01-18 2002-08-02 Mitsubishi Electric Corp ディジタルフィルタ
US20050076073A1 (en) * 2003-08-29 2005-04-07 Texas Instruments Incorporated Biquad digital filter operating at maximum efficiency
JP2005094378A (ja) * 2003-09-17 2005-04-07 Ricoh Co Ltd Iirフィルタ装置
CN101047367A (zh) * 2006-03-29 2007-10-03 三洋电机株式会社 滤波装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20010032226A1 (en) * 1998-01-29 2001-10-18 Koji Takano Digital filters
JP2001160736A (ja) * 1999-12-01 2001-06-12 Nec Ic Microcomput Syst Ltd デジタルフィルタ回路
JP2002217684A (ja) * 2001-01-18 2002-08-02 Mitsubishi Electric Corp ディジタルフィルタ
US20050076073A1 (en) * 2003-08-29 2005-04-07 Texas Instruments Incorporated Biquad digital filter operating at maximum efficiency
JP2005094378A (ja) * 2003-09-17 2005-04-07 Ricoh Co Ltd Iirフィルタ装置
CN101047367A (zh) * 2006-03-29 2007-10-03 三洋电机株式会社 滤波装置

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115102524A (zh) * 2022-07-07 2022-09-23 武汉市聚芯微电子有限责任公司 一种滤波器
CN115102524B (zh) * 2022-07-07 2023-08-08 武汉市聚芯微电子有限责任公司 一种滤波器

Similar Documents

Publication Publication Date Title
US5179531A (en) Accelerated digital signal processor
CN109194307B (zh) 数据处理方法及系统
JPH0370411B2 (zh)
US20040252829A1 (en) Montgomery modular multiplier and method thereof using carry save addition
CN114448390A (zh) 一种Biquad数字滤波器装置及实现方法
CN113556101B (zh) Iir滤波器及其数据处理方法
JPH082014B2 (ja) 多段デジタル・フィルタ
KR0147758B1 (ko) Mpeg-2 오디오 복호화기의 합성 필터
US6401106B1 (en) Methods and apparatus for performing correlation operations
JP4630056B2 (ja) 畳み込み演算回路
JP2513218B2 (ja) Firデイジタルフイルタ
KR100378192B1 (ko) 연산 프로세서를 이용한 디지털 베이스 부스터
US7603400B2 (en) Method and system for filter loop with saturation
US6944217B1 (en) Interleaved finite impulse response filter
JPH0998069A (ja) Fir型ディジタルフィルタ
JP2000293357A (ja) マイクロプロセッサ
JPH0773022A (ja) ディジタル信号処理方法及びその装置
US20020010727A1 (en) Area efficient fir filter with programmable coefficients
JP2004128858A (ja) Firデジタルフィルタ
KR0176468B1 (ko) 대역 통합 필터
US6993551B2 (en) Sparse-coefficient functions for reducing computational requirements
CN115603707A (zh) 通用滤波器、协处理器及滤波方法
JPH06309164A (ja) ディジタル信号処理回路
JPH03211910A (ja) ディジタルフィルタ
CN115865044A (zh) 支持浮点系数的定点iir滤波器和电子设备

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20220506

RJ01 Rejection of invention patent application after publication