CN1706100A - 数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器 - Google Patents

数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器 Download PDF

Info

Publication number
CN1706100A
CN1706100A CNA2003801018286A CN200380101828A CN1706100A CN 1706100 A CN1706100 A CN 1706100A CN A2003801018286 A CNA2003801018286 A CN A2003801018286A CN 200380101828 A CN200380101828 A CN 200380101828A CN 1706100 A CN1706100 A CN 1706100A
Authority
CN
China
Prior art keywords
numerical string
filter
filtering device
mentioned
tap
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
CNA2003801018286A
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.)
Neuro Solution Corp
Original Assignee
Neuro Solution Corp
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 Neuro Solution Corp filed Critical Neuro Solution Corp
Publication of CN1706100A publication Critical patent/CN1706100A/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/06Non-recursive filters
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0223Computation saving measures; Accelerating measures
    • 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
    • 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/0286Combinations of filter structures

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Filters That Use Time-Delay Elements (AREA)
  • Picture Signal Circuits (AREA)

Abstract

通过借助任意组合将给定的非对称型的数值串作为滤波器系数H1~H3的单元滤波器(L10’、H10’)来设计数字滤波器,能够自动地获得仅通过单元滤波器的组合求得的数字滤波器的系数。另外,通过令作为非对称型的滤波器系数H1~H3,使用与将对称型的数值串{-1,0,9,16,9,0,-1}/32从其中央分成两半后的其中一半,从而减少所设计的数字滤波器所需的抽头数,同时不须使用窗函数,并且所获得的滤波器特性中不会产生截尾误差。

Description

数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器
技术领域
本发明是关于数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器,特别是关于具备有由多个迟延器构成的带抽头的延迟线,并在将各抽头信号分别变为数倍之后,相加后输出的FIR滤波器及其设计方法。
背景技术
在提供给通信、测量、声音、影像信号处理、医疗、地震学等各种领域的各种电子机器中,经常要在其内部实施某种数字信号处理。数字信号处理的最重要基本操作中,有从混合着各种信号或噪声的输入信号之中,仅取出必要的某个频段信号的滤波处理。所以,在进行数字信号处理的电子机器中,经常采用数字滤波器。
作为数字滤波器,经常采用IIR(Infinire Impulse Response,无限脉冲响应)滤波器、或FIR(Finite Impulse Response,有限脉冲响应)滤波器。其中,FIR滤波器具有下述优点。第1,由于FIR滤波器的传递函数的极点仅在z平面的原点处才有,因此电路总处于稳定。第2,若滤波器系数为对称型,便可实现完全正确的直线相位特性。
若将滤波器根据通过区域与阻止区域的配置进行分类的话,主要可分为:低通滤波器、高通滤波器、带通滤波器、带阻滤波器这四种。基本用IIR滤波器与FIR滤波器构成的是低通滤波器,而其他的高通滤波器、带通滤波器、带阻滤波器,则通过进行频率变换等处理从低通滤波器导出。
但是,FIR滤波器,由有限时间长度表示的脉冲响应,直接当作滤波器的系数。所以,所谓设计FIR滤波器,是指决定滤波器系数以获得期望的频率特性。
以往,在设计构成基本的低通滤波器时,通过根据采样频率与截止频率的比率,进行利用窗函数与切比雪夫近似法等的卷积运算,求取与FIR滤波器的各抽头对应的滤波器系数。然后,通过用所求得滤波器系数进行模拟,一边确认频率特性,一边适当修正系数值,获得所需特性的低通滤波器。
另外,在设计高通滤波器、带通滤波器、带阻滤波器等其他滤波器时,首先按照上述次序设计多个基本的低通滤波器。然后,通过将其组合后进行频率变换等操作,设计出具有期望的频率特性的FIR滤波器。
然而,由于依现有设计方法所获得滤波器的频率特性,依赖于窗函数或切比雪夫近似式,因此若未将其良好设定,无法获得良好的频率特性。可是,一般来说难以适当设定窗函数与近似式。即,上述现有的滤波器设计方法,熟练的技术人员需要大量时间和精力进行设计,造成设计所需特性的FIR滤波器较为不便的问题。
再者,即便假设可设计出接近所需特性的滤波器,所设计的滤波器的抽头数也过于庞大,且其系数值为非常复杂随机的数值。因此问题就是,若要实现此抽头数与系数值,需要大规模的电路结构(加法器、乘法器)。此外,实际使用所设计的FIR滤波器时,运算量非常巨大,导致处理负荷过重的问题。
本发明正是为解决上述问题,其目的在于,能较为方便地设计具所需频率特性的FIR数字滤波器。
再者,本发明的目的在于,提供一种可用较小的电路规模高精度实现的FIR滤波器。
发明内容
本发明下的数字滤波器的设计方法,其特征在于:
通过使用:具有将值设定为数值串为非对称型、上述数值串的总和非零、并且上述数值串的隔1相加之和为同号且相等的基本的滤波器系数的第1单元滤波器,和具有将值设定为数值串为非对称型、该数值串的总和为零、并且该数值串的隔1相加之和为反号且绝对值相等的基本的滤波器系数的第2单元滤波器之中的至少一种,并借助任意级联连接0个以上的上述第1单元滤波器和0个以上的上述第2单元滤波器,来进行滤波器设计。
本发明的另一个方式中,通过使用:将具有将值设定为数值串的总和非零、该数值串的隔1相加之和为同号且相等的非对称型的滤波器系数的多个基本滤波器,级联连接为整体上的数值串构成对称型而形成的第1单元滤波器,和将具有将值设定为数值串的总和为零、该数值串的隔1相加之和为异号且相等的非对称型的滤波器系数的多个基本滤波器,级联连接为整体上的数值串构成对称型而形成的第2单元滤波器之中的至少一种,并借助任意级联连接0个以上的上述第1单元滤波器和0个以上的上述第2单元滤波器,来进行滤波器设计。
在本发明的另一个形式中,其特征为:通过在与构成上述第1及第2单元滤波器的非对称型的滤波器系数相对应的各抽头间,插入n时钟量的延迟,来调整滤波器的通过频带。
另外,本发明下的数字滤波器,其特征为:具备由多个延迟器所形成的带抽头的延迟线,通过用上述之一的滤波器设计方法求得的滤波器系数、分别将各抽头的信号变为数倍后,相加输出。
本发明的另一个形式中,其特征为:用于在与上述非对称型的滤波器系数相对应的各抽头间插入n时钟量的延迟延迟模块,如下构成。即具备按照周期为基准时钟的1/n倍的第2时钟、将输入的数据保持输出的延迟电路
如上述说明,依据本发明,由于通过任意组合1个以上的以给定的非对称型的数值串为滤波器系数的单元滤波器,来设计数字滤波器,因此能够仅通过单元滤波器的组合,就能自动地获得具有期望的频率特性的数字滤波器的系数,即不是熟练的技术人员,也可以简单地进行滤波器的设计。
另外,由于依据本发明,上述给定的非对称型的数值串,对应于将给定的对称型的数值串从其中央分出的一半,因此设计出的数字滤波器所需要的抽头数非常少,而且,对于各抽头输出来说必要的滤波器系数的种类也非常少。而且,不须使用窗函数,所获得的滤波器系数也不会产生截尾误差。因此,可实现大幅缩减电路元件数(特别是乘法器)减小电路规模、降低耗电量、减轻运算负荷等,同时,能高精度地实现期望的数字滤波器的频率特性。此外,由于设计的数字滤波器,为通过反复基本滤波器这个相同样式来构成具有极单纯的结构,因此在集成化时可减少工序数,便于IC化。
附图说明
图1,是表示5抽头单元滤波器L10、H10的电路结构及滤波器系数的数值串的图。
图2是表示5抽头单元滤波器L10、H10的滤波器系数的生成算法的图。
图3是用于说明5抽头低通单元滤波器L10的滤波器系数的意义的图。
图4是用于说明5抽头高通单元滤波器H10的滤波器系数的意义的图。
图5是表示5抽头低通单元滤波器L11的滤波器系数的生成算法的图。
图6是表示5抽头低通单元滤波器L10、L11的频率-增益特性的图。
图7是表示5抽头高通单元滤波器H11的滤波器系数的生成算法的图。
图8是表示5抽头高通单元滤波器H10、H11的频率-增益特性的图。
图9是表示5抽头低通单元滤波器(L10)m的频率-增益特性的图。
图10是表示5抽头高通单元滤波器(H10)m的频率-增益特性的图。
图11是用于说明频段的抽取的图。
图12是表示另一个频段的抽取示例的图。
图13是表示用第1实施方式下的滤波器设计方法、构成最基本的2种3抽头单元滤波器的电路结构及滤波器系数数值串的图。
图14是表示将简单数值串{8、-9、0、1}/16作为滤波器系数的情况下、3抽头高通滤波器H10’的频率特性的图。
图15是表示当仅调整了滤波器系数H3时、3抽头高通单元滤波器H10’的频率特性的图。
图16是表示当调整了滤波器系数H2、H3时、3抽头高通单元滤波器H10’的频率特性的图。
图17是表示第1实施方式下的带通滤波器的电路示例的图。
图18是表示图17所示的带通滤波器的频率特性的图。
图19是表示第2实施方式下的单元滤波器L10”、H10”的电路结构及滤波器系数的数值串的图。
图20是表示第2实施方式下的低通单元滤波器L10”的频率特性的图。
图21是表示第2实施方式下的高通单元滤波器H10”的频率特性的图。
图22是表示5抽头单元2次滤波器L20、H20的滤波器系数的生成算法的图。
具体实施方式
以下,参照附图对本发明一个实施方式进行说明。
本实施方式的数字滤波器,是具备由多个延迟器所构成的带抽头的延迟线,将各抽头的输出信号通过被赋给的滤波器系数分别变为数倍之后,相加输出类型的FIR滤波器。
(第1实施方式)
第1实施方式下的滤波器设计方法,制作下述所说明的2种3抽头单元滤波器L1n’、H1n’,并能仅通过它们的组合设计具期望频率特性的FIR滤波器。另外,在表示单元滤波器的符号后所附加的“n”字,是表示插入于各抽头间的延迟的时钟数,即插入于各滤波器系数间的“0”的数量(容后详述)。
首先,对作为用于理解上述2种3抽头单元滤波器L1n’、H1n’的参考的2种5抽头单元滤波器L1n、H1n进行说明。图1是表示5抽头单元滤波器L10、H10图,(a)表示其电路结构,(b)表示滤波器系数的数值串。
如图1(a)所示,5抽头单元滤波器L10、H10,是利用级联连接的6个D型触发器1-1~1-6,将输入信号逐个时钟CK依次延迟。然后,对于从各D型触发器1-1~1-6的给定抽头中抽取出的信号,将图1(b)所示的滤波器系数h1~h5分别与5个系数器2-1~2-5相乘,并将其乘法结果全部利用4个加法器3-1~3-4相加后输出。
上述2种5抽头单元滤波器L10、H10的电路结构,均如图1(a)所示,仅滤波器系数(系数器2-1~2-5的乘数值h1~h5)不同,如图1(b)所示。
由图1(b)得知,5抽头低通单元滤波器L10的的滤波器系数,是由极单纯的数值串{-1,0,9,16,9,0,-1}/32所构成(其中,数值为“0”的部分如图1(a)所示不是抽头输出,没有用作滤波器系数)。此种滤波器系数,是其数值串为对称型,具有数值串的总和非零、且数值串隔1相加的总和同号且相等的性质(-1+9+9-1=16,0+16+0=16)。
另外,5抽头高通单元滤波器H10的滤波器系数,由极单纯的数值串{1,0,-9,16,-9,0,1}/32构成(其中,数值“0”的部分不是抽头输出,没有用作滤波器系数)。此种滤波器系数,其数值串为对称型,具有数值串的总和为零、且数值串隔1相加之和反号且绝对值相等的性质(1-9-9+1=-16,0+16+0=16)。
在此,对构成这些滤波器系数的数值串的涵义,用图2~4进行说明。
图2为,用于对构成5抽头单元滤波器L10、H10的滤波器系数的数值串的生成法进行说明的图。如图2(a)所示,构成5抽头低通单元滤波器L10的滤波器系数的数值串,通过对在每1时钟CK数据值变化为{-1,1,8,8,1,-1}/16的给定数字基本函数,进行1次移动平均运算求得。
另外,如图2(b)所示,构成5抽头高通单元滤波器H10的滤波器系数的数值串,通过对在每1时钟CK数据值变化为{1,-1,-8,8,-1,1}/16、和每1时钟CK数据值变化为{1,-1,8,-8,-1,1}/16的数字基本函数,进行移动平均运算求得。
图3是表示对构成5抽头低通单元滤波器L10的滤波器系数的数值串,实施4倍过采样与卷积(convolution)运算的结果图。另外,在此为了令说明便于理解,仅表示的是对将原数值串变为32倍的整数数值串{-1,0,9,16,9,0,-1}实施过采样与卷积运算的示例。
在图3(a)中,最左列所示的一串数值串,是对原数值串{-1,0,9,16,9,0,-1}实施4倍过采样后的数值。此外,从最左边起向右侧4列的数值串,将最左列所示数值串向下逐一移位。图3(a)的列方向表示时间轴,所谓将数值串向下移位,与将最左列所示数值串逐渐延迟相对应。
换句话说,左起第2列数值串,表示将最左列表示的数值串偏移4倍频率的时钟4CK的1/4相位的数值串。此外,左起第3列的数值串,表示将左起第2列所示的数值串偏移4倍频率的时钟4CK的1/4相位的数值串,左起第4列数值串,表示将左起第3列所示的数值串再偏移4倍频率的时钟4CK的1/4相位的数值串。
另外,左起第5列数值串,是将对应第1~4列的各数值串的各行相加后的数值。通过到左起第5列为止的处理,数字式地执行伴随4相卷积运算的4倍过采样。
从第5列起向右4列的数值串,将第5列所示数值串向下逐一移位。此外,左起第9列数值串,是将与第5~8列各数值串对应的各行相加后的数值。通过至左起第9列为止的处理,数字式地2次执行伴随4相卷积运算的4倍过采样。
另外,左起第10列的数值串,将第9列所示的数值串向下移1位。此外,左起第11列(最右列)数值串,是用对应第9列数值串与第10列数值串的各行相加后的数值。
将该图3(a)的最右列所示的最后获得的数值串图形化,则如图3(b)所示。具有此图3(b)所示波形的函数,是仅当沿横坐标轴的采样位置处于t1至t4间时,有“0”以外的有限数值,处此之外的区域的数值全部为“0”的函数,即、是在采样位置t1、t4中数值收敛于“0”的函数。将上述函数的数值在局部区域为“0”以外的有限数值、而在此外的区域中为“0”的情况,称之为“有限基(finitebase)”。
另外,图3(b)所示函数,是具有如下特征,即仅在中央的采样位置t5取得极大值,而在t1、t2、t3、t4这4个采样位置中数值为“0”的采样化函数,为了获得平滑的波形数据,通过全部必要的采样。
然后,图4是,表示对构成5抽头高通单元滤波器H10的滤波器系数的数值串,实施4倍过采样与卷积运算的结果的图。另外,这里也为了令说明便于理解,仅表示对将原数值串变为32倍的整数数值串{1,0,-9,16,-9,0,1}实施过采样与卷积运算的示例。
图4(a)表示与上述图3(a)相同的运算过程。将图4(a)的最右列所示的最终获得的数值串图形化,如图4(b)所示。此图4(b)所示的函数,也是仅在中央的采样位置t7’处取得极大值的采样化函数,是可在整个区域中进行1次微分,并且在采样位置t1’、t6’中收敛于0的有限基的函数。
下面,对将插入各抽头间的延迟时钟数n设定为n≥1的情况进行说明。图5是,表示5抽头低通单元滤波器L11(n=1的情况)的滤波器系数的图。如图5所示,5抽头低通单元滤波器L11的滤波器系数,是通过在上述5抽头低通单元滤波器L10的各滤波器系数之间,各插入1个“0”生成。
同样的,5抽头低通单元滤波器L1n(n=2,3,…)的滤波器系数,通过在上述5抽头低通单元滤波器L10的各滤波器系数之间,各插入n个“0”生成。
图6是,表示将5抽头低通单元滤波器L10、L11的数值串,FFT变换(Fast Fourier Transfer:高速傅立叶变换)后的结果的频率-增益特性的图。这里,将增益与频率用“1”标准化。
由此图6可知,5抽头低通单元滤波器L10、L11可获得良好的低通滤波器特性,即中心频率处增益为0.5,且不存在低频区域下的超调量(overshoot)或高频区域下的振铃(ringing)。此外,若令插入各滤波器系数之间的“0”的数量为n,其频率-增益特性的频率坐标轴(频率方向对应的周期)为1/n。
作为实现此种低通滤波器特性的基础的、上述图2(a)的数值串,构成图3(b)所示的有限基的采样化函数的基础。与现在一般所采用的采样化函数在t=±∞的采样位置收敛于“0”不同,图3(b)所示的采样化函数,是在t=t1、t4的有限采样位置处收敛于“0”。
因此,当对上述图2(a)的数值串实施FFT变换时,仅相当于t=t1~t4范围内的数据有意义。由于对于相当于t=t1~t4范围外的数据来说,本来应当考虑并不是要忽略,但因为理论上并无考虑的必要,因此不会发生截尾误差。所以,若将上述图2(a)所示的数值串使用作为滤波器系数,不必用窗函数进行系数的截尾,能够获得良好的低通滤波器特性。
图7是,表示5抽头高通单元滤波器H11的滤波器系数的图。如图7所示,5抽头高通单元滤波器H11的滤波器系数,通过在上述5抽头高通单元滤波器H10的各滤波器系数之间,各插入1个“0”生成。
同样,5抽头高通单元滤波器H1n(n=2,3,…)的滤波器系数,通过在上述5抽头高通单元滤波器H10的各滤波器系数之间,各插入n个“0”生成。
图8是,表示5抽头高通单元滤波器H10、H11的频率-增益特性的图。这里,也将增益与频率利用“1”标准化。由此图8可知,5抽头高通单元滤波器H10、H11可获得良好的高通滤波器特性,即在中心频率处增益为0.5,且不存在高频区域下的超调量或低频区域下的振铃。此外,若令插入各滤波器系数间的“0”的数量为n,其频率-增益特性的频率坐标轴(频率方向对应的周期)为1/n。
作为实现此种高通滤波器特性的基础、上述图2(b)的数值串,也构成图4(b)所示的有限基采样化函数的基础。所以,通过将此数值串使用作为滤波器系数,无需用窗函数进行系数截尾,并可获得良好的高通滤波器特性。
下面,对5抽头单元滤波器的级联连接进行说明。通过将5抽头单元滤波器级联连接,在各单元滤波器的系数间进行乘法·加法生成新的滤波器系数。以下,例如将5抽头低通单元滤波器L10的级联数设定为m时,并将其记为(L10)m
图9是,表示5抽头低通单元滤波器L10、(L10)2、(L10)4、(L10)8的频率-增益特性的图。此图9也是将增益及频率用“1”标准化。当只有1个5抽头低通单元滤波器L10时,振幅为0.5的位置的时钟为0.25。与此相对,若级联数m变多,滤波器的带通宽度变窄。例如,当m=8时,振幅为0.5的位置的时钟为0.125。
图10表示的是,5抽头高通单元滤波器H10、(H10)2、(H10)4、(H10)8的频率-增益特性。此图10也将增益及频率用“1”标准化。当只有1个5抽头高通单元滤波器H10时,振幅为0.5的位置的时钟为0.25。与此相对,若级联数m变多,滤波器的带通宽度变窄。例如当m=8时,振幅为0.5的位置的时钟为0.375。
下面,对期望频带的抽取进行说明。图11是用于对频带抽取进行说明的图。频带的抽取,用级联连接4个以上上述单元滤波器的产物来进行。图11(a)是,将5抽头低通单元滤波器(L10)4、(L11)4、(L13)4、(L17)4的频率-增益特性,整合在1个图形上表示的图。此图11(a)也将增益及频率用“1”标准化。
将这多种5抽头低通单元滤波器(L10)4、(L11)4、(L13)4、(L17)4组合后,各特性值之间相互抵消结合来进行频带抽取。此外,通过将这些波形作为基本来组合、并根据需要向反转频率坐标轴方向移动,从而能够作成仅期望频带为通过频段的滤波器。
图11(b)是,将这4种5抽头低通单元滤波器(L10)4、(L11)4、(L13)4、(L17)4级联连接时所获得的频率特性。据此,获得低频的极窄区域为通过频段、具有基本良好的衰减特性的低通滤波器。虽然发生些许振铃,但是因为此振铃发生在落入-106dB以上的部分,因此几乎可忽视。
图12是,表示频带的抽取的另一个示例的图。如图12(a)所示,将5抽头高通单元滤波器(H10)8、与5抽头低通单元滤波器(L11)8、(L13)8、(L17)8组合并将其级联连接后,可获得给定频带为通过频段的高通滤波器。
再者,如图12(b)所示,将5抽头低通单元滤波器(L10)8、(L11)8、(L13)8组合并将其级联连接后,可获得给定频带为通过频段的低通滤波器。此外,如图12(c)所示,将5抽头高通单元滤波器(H11)8、与5抽头低通单元滤波器(L13)4、(L17)4、(L115)8组合并将其级联连接后,可获得给定频带为通过频段的带通滤波器。
如上述,通过用将给定的基本数值串设定为滤波器系数的5单元滤波器并将其任意组合,仅由单元滤波器的组合,便可自动生成具有期望的频率特性的FIR滤波器的滤波器系数。所以,简化了滤波器设计方法,即便不是熟练的技术者,也能很容易地进行滤波器设计。
再者,因为在采用上述手法所设计的滤波器电路中,所需要的抽头数量非常少,且对各抽头输出所需的滤波器系数种类也极少,因此能够简化滤波器电路的运算部的结构。所以,能够大幅削减电路元件数(特别是乘法器)、缩小滤波器电路的规模,而且可实现降低消耗功率、减轻运算负荷等。
再者,因为采用上述手法所设计的滤波器电路,通过反复相同的样式形成、具有极简单的结构,因此在集成化时可缩短工序数量,且具有便于IC化的优点。此外,在特性面上,截止特性可获得极大的改善,相位特性也为直线、能够获得良好的滤波器特性。
对于上述所说明的内容,本案申请人已然提出专利申请案(日本专利特愿2001-321321号)。本实施方式,是对此已提出申请的内容进行进一步改良,如上所述,能只通过2种3抽头单元滤波器L1n’、H1n’的组合,设计具有期望的频率特性的FIR滤波器。上述3抽头单元滤波器L1n’、H1n’的滤波器系数,通过将上述的5抽头单元滤波器L1n、H1n的滤波器系数的数值串从其中央分为两半、并对其中一半进行进一步调整后得到。
图13是,表示3抽头单元滤波器L10’、H10’的图,(a)表示其电路结构,(b)表示滤波器系数数值串。如图13(a)所示,3抽头单元滤波器L10’、H10’,通过级联连接的3个D型触发器11-1~11-3,将输入信号依次逐个时钟CK延迟。然后,对从各D型触发器11-1~11-3的给定抽头中抽取出的信号,将图13(b)所示的滤波器系数H1~H3,用3个系数器12-1~12-3分别相乘,并将其全部乘法结果用2个加法器13-1~13-2相加后再输出。
上述2种3抽头单元滤波器L10’、H10’的电路结构,均如图13(a)所示,仅滤波器系数(系数器12-1~12-3的乘数值H1~H3)不同,如图13(b)所示。
3抽头低通单元滤波器L10,的滤波器系数,其数值串为非对称型,该数值串的总和非零,且具有该数值串隔1相加之和同号相等的性质。此外,3抽头高通单元滤波器H10的滤波器系数,其数值串为非对称型,数值串的总和为零,且具有该数值串隔1相加之和异号且绝对值相等的性质。
如上所述,这3抽头单元滤波器L10’、H10’的滤波器系数,除了非对称型这点之外,其它性质均与上述5抽头单元滤波器L10、H10相同。但是,与5抽头单元滤波器L10、H10的滤波器系数相比,数值变得有些复杂。以下,说明此理由。
图14,是表示将5抽头高通单元滤波器H10的滤波器系数{1,0,-9,16,-9,0,1}/32,从其中央处分为两半,并将相当于其单侧的简单数值串{8,-9,0,1}/16作为滤波器系数时(其中,数值为“0”部分不是抽头输出,不作为滤波器系数使用)的频率特性的图。在此,将增益及频率用“1”标准化。
如该图14所示,若将5抽头高通单元滤波器H10的数值串仅简单分出一半的数值串作为滤波器系数使用,在其频率特性中构成通过频段的部分的峰值将波动并产生多个极大值,且其极大值会超过“1”。具有此种频率特性的单元滤波器,不适用于将其多个级联连接来对期望的FIR滤波器进行设计的方式。
因此,调整上述简单的数值串{8,-9,0,1}/16。首先,缩小决定高频区域侧的频率特性的滤波器系数H3的绝对值。换句话说,通过将其系数值从“1”变更为“1-N/8”(N=1,2,…,8之一)来减少高频区域成分,从而仅令通过频段的正中央为增益最大值。
图15是,表示在滤波器系数H3设定为N=3时的频率特性的图。由此图15可知,通过频段的波动现象获得改善,且仅中央部分为增益最大值。但是,最大值仍然超过“1”。所以,接着进一步调整系数值,以令此增益的最大值恰好为“1”。
最大值的调整,可用与高频区域成分的调整所用的滤波器系数H3符号相反的滤波器系数H2来进行。这里,将滤波器系数H2的绝对值减小,从“-9”变为“-(9-N/8)”(N=1,2,…,8之一)。
此时,通过将增益调整用的N值(滤波器系数H2的N值),设定为与上述的高频区域调整用的N值(滤波器系数H3的N值)相同,令调整后的数值串的总和与调整前的数值串的总和相同(调整前:8-9+0+1=0,调整后:8-(9-N/8)+0+(1-N/8)=0)。此外,令数值串的隔1相加之和也在调整前后不发生改变。
图16,是表示对于滤波器系数H2、H3设定为N=3时的频率特性的图,(a)将增益用直线刻度表示,(b)将增益用对数刻度表示。此图16中,也将增亦及频率用“1”标准化。
由此图16中可知,调整滤波器系数H2、H3后,通过频段中的波动消失且仅有中央部分为增益最大值,同时其最大值恰好为“1”。此外,还获得约-55dB的良好的衰减量。具此种频率特性的3抽头高通单元滤波器H10’,适用于将其级联连接来设计期望的FIR滤波器的方式。
再者,与上述5抽头高通单元滤波器H1n(n=1,2,…)相同,通过在3抽头高通单元滤波器H10’的各滤波器系数之间插入n个“0”,可生成3抽头高通单元滤波器H1n’。
与上述3抽头高通单元滤波器H10’相同,对3抽头低通单元滤波器L10’的滤波器系数也进行适当调整。换句话说,对于将5抽头低通单元滤波器L10的数值串仅简单分出一半的数值串{8,9,0,-1}/16,减小绝对值,以令决定高频区域侧的频率特性的滤波器系数H3的数值从“-1”变为“-(1-N/8)”。
另外,用符号与此滤波器系数H3相反的滤波器系数H2,调整增益最大值。换句话说,将滤波器系数H2值从“9”减小为“9-N/8”。此时,通过将高频区域调整用的N值设定为与增益调整用的N值相同,令调整后的数值串的总和与调整前的数值串的总和为相同数值(调整前:8+9+0-1=16,调整后:8+(9-N/8)+0-(1-N/8)=16)。另外,令数值串的隔1相加之和也在调整前后不发生改变。
如上所述,当对3抽头低通单元滤波器L10’的滤波器系数进行调整后,消除了通过频段中的波动、且仅中央部分为增益最大值,并且,能够获得其最大值恰好为“1”的低通滤波器特性。具有此种频率特性的3抽头低通单元滤波器L10’,也适用于将其级联连接来对期望的FIR滤波器进行设计的方式。
另外,与上述5抽头低通单元滤波器L1n(n=1,2,…)相同,通过在3抽头低通单元滤波器L10’的各滤波器系数之间插入n个“0”,生成3抽头低通单元滤波器L1n’。
如上所述,通过将如图13所示构成的2种3抽头单元滤波器L10’、H10’适当组合后级联连接,可与将5抽头单元滤波器L10、H10组合时相同,简单地设计具有期望频率特性的FIR滤波器。
另外,3抽头单元滤波器L10’、H10’的滤波器系数,也可以将5抽头单元滤波器L10、H10的滤波器系数的数值串从其中央分为两半后,通过调整另一半数值串{-1,0,9,8}/16、{1,0,-9,8}/16来生成。
下面,表示使用了上述3抽头单元滤波器L10’、H10’的带通滤波器的设计例。设计的带通滤波器的目标规格如下,即带通滤波器的中心频率Fc为675KHz,信号的采样频率Fs为2.7MHz(Fs=4*Fc),-3dB的带宽为100KHz,-80dB的带宽为200KHz,带宽外衰减量为-80dB以上。
图17是表示实现此目标规格的带通滤波器的电路示例的图。图17所示的电路,是通过将6个3抽头高通单元滤波器H11’、4个3抽头低通单元滤波器L13’、3个3抽头高通单元滤波器H15’、3个3抽头低通单元滤波器L17’、2个5抽头高通单元滤波器H19,依此顺序级联连接构成。
图18是表示图17所示的带通滤波器的频率特性的图,(a)以直线刻度表示增益,(b)以对数刻度表示增益。由图18可知,图17所示的带通滤波器满足上述目标规格。另外,由于3抽头单元滤波器L1n’、H1n’,其滤波器系数的数值串为非对称型,因此其单独无法保证如图16所示相位的直线性。但是,将其聚集多个并如图17所示级联连接后,相位的偏差消除,如图18所示相位特性获得一些改善。
通过图17及图18可知,通过本实施方式的滤波器设计法,以抽头数更少的3抽头单元滤波器L1n’、H1n’为基础,并通过将其级联连接得到的极简单的结构,获得期望特性的带通滤波器。此时,整体上所必要的抽头数,每一位仅为51个抽头。
即,此带通滤波器的滤波器系数,依据上述有限基的函数而产生,在级联连接的两侧附近,系数值小至可忽视的程度。由于根据有限基的性质,即使忽视此部分,也不会产生截止误差,因此将该部分去除。去除后剩余的部分,构成了关于满足目标规格的带通滤波器求出的最终的滤波器系数。
因此,实际上,可将此滤波器系数构成作为硬件,且必要的抽头数仅为51。从而,与只用5抽头单元滤波器L1n、H1n时相比,整体上所需的抽头数减少,能显著简化滤波器电路的结构。
另外,此处虽就利用3抽头单元滤波器L1n’、H1n’以设计带通滤波器的示例进行了说明,但是也可以通过同样的方法设计高通滤波器、低通滤波器、带阻滤波器等。
(第2实施形态)
接着,对本发明的第2实施方式进行说明。在上述第1实施方式中,例如作为3抽头低通单元滤波器L10’的滤波器系数,使用的是将构成5抽头低通单元滤波器L10的滤波器系数的数值串{-1,0,9,16,9,0,-1}/32分成两半后、进一步调整其中一半的数值串(8,9,0,-1)/16得到的数值串。
另外,作为3抽头高通单元滤波器H10’的滤波器系数,使用的是将构成5抽头高通单元滤波器H10的滤波器系数的数值串{1,0,-9,16,-9,0,1}/32分成两半后、进一步调整其中一半的数值串{8,-9,0,1}/16得到的数值串。由于这3抽头单元滤波器L10’、H10’的滤波器系数均是非对称型,因此如上所述无法确保相位的直线性。
与此相对,第2实施方式,使用抽头数减少的3抽头单元滤波器,也能实现直线相位特性。
图19是,表示第2实施方式下的单元滤波器L10”、H10”的图,(a)表示其电路结构,(b)表示滤波器系数的数值串。
如图19(a)所示,本实施方式的单元滤波器L10”、H10”均具有相同结构。换句话说,低通单元滤波器L10”,级联连接2个3抽头低通单元滤波器1L10、2L10构成。另外,高通单元滤波器H10”,级联连接2个3抽头高通单元滤波器1H10、2H10构成。
如图19(b)所示,构成低通单元滤波器L10”的一个3抽头低通单元滤波器2L10,将上述5抽头低通单元滤波器L10的滤波器系数数值串{-1,0,9,16,9,0,-1}/32分为两半,再进一步调整其中一半的数值串{8,9,0,-1}/16后,将其设定为滤波器系数H4~H6。此与第1实施方式所说明的3抽头低通单元滤波器L10’相同。此外,另一个3抽头低通单元滤波器1L10,再进一步调整如上述分为两半的另一半的数值串{-1,0,9,8}/16后,将其设定为滤波器系数H1~H3。
另外,构成高通单元滤波器H10”的一个3抽头高通单元滤波器2H10,将上述5抽头高通单元滤波器H10的滤波器系数的数值串{1,0,-9,16,-9,0,1}/32分为两半,再进一步调整其中一半的数值串{8,-9,0,1}/16后,将其设定为滤波器系数H4~H6。此与第1实施方式所说明的3抽头高通单元滤波器H10’相同。另外,另一个3抽头高通单元滤波器1H10,再进一步调整如上述分为两半的另一半的数值串{1,0,-9,8}/16后,将其设定为滤波器系数H1~H3。
再者,上述两个3抽头低通单元滤波器1L10、2L10的连接关系、和上述两个3抽头高通单元滤波器1H10、2H10的连接关系,可与图19(a)所示的形态左右相反。
由于若如上述构成单元滤波器L10”、H10”,滤波器系数变为对称型,因此相位特性呈直线。图20是表示低通单元滤波器L10”的频率特性的图,图21是表示高通单元滤波器H10”的频率特性图,并将增益用对数刻度表示。此图20与图21,也将增益及频率用“1”标准化。
由这些图可知,频率-增益特性中、通过频段下没有波动,且最大值恰好为“1”。此外,获得约-55dB的良好的衰减量。还有,频率-相位特性获得非常好的直线相位特性。所以,通过将具有此种频率特性的单元滤波器L1n”、H1n”如第1实施方式所说明那样级联连接,能够用非常少的抽头数构成系数对称的直线相位滤波器。
但是,图19所示单元滤波器L10”、H10”,总共具备6个抽头,抽头数比5抽头单元滤波器L10、H10多。可是,相对于5抽头单元滤波器L10、H10由1段构成,单元滤波器L10”、H10”由两个3抽头单元滤波器级联连接构成。所以,带宽已然比5抽头单元滤波器L10、H10窄。
因此第2实施方式,特别适用于要设计带宽较窄的FIR滤波器的情况。换句话说,为了实现期望的窄带宽,整体所需要的单元滤波器级联数,比采用5抽头单元滤波器L10、H10的情况大幅减少。因此,从整体来看能够缩小电路规模。
用以实现上述说明的第1及第2实施方式下的数字滤波器的设计方法的装置,可通过硬件结构、DSP、软件之一实现。例如当用软件实现时,滤波器设计装置由CPU或MPU、RAM、ROM等构成,并通过运行存储于RAM、ROM或硬盘等中的程序来实现。
所以,可通过将会被运行以令计算机具有上述本实施方式的功能的程序,储存于例如CD-ROM之类的记录媒体中,并将其读入计算机来实现。作为记录上述程序的记录媒体,除CD-ROM之外,还可采用例如:软盘、硬盘、磁带、、光盘、光磁盘、DVD、非易失性存储卡等。此外,还可通过将上述程序通过互联网等网路,下载到计算机中实现。
即,将与各种单元滤波器的滤波器系数作为信息,储存于RAM或ROM等存储器中,用户指示与单元滤波器相关的任意组合后,CPU用上述存储器中储存的滤波器系数的信息,运算指示的组合所对应的滤波器系数,求得FIR滤波器。
例如,可预先将各种单元滤波器图标化(对应各图标将滤波器系数作为信息储存),用户通过在显示器画面上任意组合配置这些图标,CPU自动运算并求取对应该排列的滤波器系数。此外,若将所求得的滤波器系数自动地进行FFT变换,并将其结果显示为频率-增益特性图,能够确认所设计的滤波器的特性,更方便地进行滤波器设计。
再者,除了通过计算机执行被提供的程序,实现上述实施方式的功能,还可用该程序与电脑中运行的OS(操作系统)或其他应用程序软件等共同实现上述实施方式的功能,或者将被提供的程序的处理的全部或一部分,用计算机的功能扩展口或功能扩展单元实施,来实现上述实施方式的功能,此时相关的程序也包含于本发明的实施方式。
(其他实施方式)
虽然在上述第1及第2实施方式中,将图2(a)及(b)所示的对称型数值串分出一半作为非对称型3抽头单元滤波器的滤波器系数,但是基础对称型数值串并不仅限于此。也可使用如图22(a)所示生成的5抽头低通单元2次滤波器L20、或如图22(b)所示生成的5抽头高通单元2次滤波器H20的对称型数值串。
这些个5抽头单元2次滤波器L20、H20的频率-增益特性中,截止特性虽然陡峭,但是在肩部的部分(通过频带)产生一些隆起。因此,这些个5抽头单元2次滤波器L20、H20,在级联连接的滤波器块的最终段中,优选作为波形调整用。
可将构成具有此种特性的5抽头单元2次滤波器L2n、H2n的滤波器系数的数值串分成两半、并进一步对其中一半进行调整,作为3抽头单元2次滤波器L2n’、H2n’的滤波器系数。另外,通过将3抽头单元2次滤波器L2n’、H2n’镜连接为对称型,也可以构成实现直线相位特性的单元2次滤波器L2n”、H2n”。
另外,除上述图2及图24所示之外,还可用绝对值为“1”和“8”的数值,形成不同于上述数值串的数值串作为5抽头单元滤波器的滤波器系数,并将此数值串分出一半作为3抽头单元滤波器的滤波器系数。
另外,在上述实施形态中,是就通过组合多数种单元滤波器,借助各特性值彼此相抵消来进行频带的抽取的示例进行了说明。相对于此,当随意决定带通滤波器的中心频率Fc或者信号的采样频率Fs的其中一种时,通过将频率抽取的条件最佳化,可以更进一步简化滤波器的结构。
现在,将带通滤波器的中心频率Fc和信号的采样频率Fs的关系设为:
Fs=4*Fc。当Fc=450KHz时,Fs=1.8MHz。在如此设定的情况下,不用通过组合低通单元滤波器和高通单元滤波器的特性值彼此的相抵消来抽取必要的频带,而例如只通过高通单元滤波器的级联连接便可以直接设计带通滤波器。
例如,只通过5抽头高通单元滤波器H11的级联连接,能够构成FM用带通滤波器,可通过其连接数m来调整带宽(通过带宽)。可将此应于3抽头高通单元滤波器H11’或者高通单元滤波器H11”等中。
再有,通过在如上所述构成的FM用带通滤波器的后段、级联连接5抽头低通单元滤波器(L1n)m,也可以构成AM用带通滤波器。对此来说也相同,可以应用于3抽头低通单元滤波器L11’或者低通单元滤波器L11”等中。
另外,也可以通过用上述的单元滤波器设计比较简单的低通滤波器,并将其通过频带在频率坐标轴方向上移位,来设计将期望的频带作为通过频带的高通滤波器、带通滤波器、带阻滤波器等。
一般来说,低通滤波器的系数H0、H1、H-1、H2、H-2、…(系数H0处于中央构成对称型)乘上2cos(mω0)后,可获得带通滤波器的系数。此处,ω0=2πFc/Fs(Fc是滤波器的中心频率,Fs是信号的采样频率),m=0、1、-1、2、-2、…。
此外,上述实施方式,不过是表示实施本发明的具体化的一例,不可用来对本发明的技术范围作限定解释。即,本发明在不脱离其精神或主要特征的前提下,可用各种形式实施。
本发明适用于:具备由多个延迟器构成的带抽头的延迟线,并将各抽头的信号分别变为数倍后,相加输出类型的FIR数字滤波器。

Claims (30)

1.一种数字滤波器的设计方法,是对将由多个延迟器构成的带抽头的延迟线中的各抽头信号、用获取的滤波器系数分别变为数倍后、相加输出的数字滤波器进行设计的方法,其特征在于:
通过使用:具有将值设定为数值串为非对称型、上述数值串的总和非零、并且上述数值串的隔1相加之和为同号且相等的基本的滤波器系数的第1单元滤波器,和
具有将值设定为数值串为非对称型、该数值串的总和为零、并且该数值串的隔1相加之和为反号且绝对值相等的基本的滤波器系数的第2单元滤波器之中的至少一种,
并借助任意级联连接0个以上的上述第1单元滤波器和0个以上的上述第2单元滤波器,来进行滤波器设计。
2.一种数字滤波器的设计方法,是对将由多个延迟器构成的带抽头的延迟线中的各抽头信号、用获取的滤波器系数分别变为数倍后、相加输出的数字滤波器进行设计的方法,其特征在于:
通过使用:将具有将值设定为数值串的总和非零、该数值串的隔1相加之和为同号且相等的非对称型的滤波器系数的多个基本滤波器,级联连接为整体上的数值串构成对称型而形成的第1单元滤波器,和
将具有将值设定为数值串的总和为零、该数值串的隔1相加之和为异号且相等的非对称型的滤波器系数的多个基本滤波器,级联连接为整体上的数值串构成对称型而形成的第2单元滤波器之中的至少一种,
并借助任意级联连接0个以上的上述第1单元滤波器和0个以上的上述第2单元滤波器,来进行滤波器设计。
3.根据权利要求1或2的任一项所述的数字滤波器的设计方法,其特征在于:
构成上述第1单元滤波器的非对称型的滤波器系数,是将给定的对称型的数值串从其中央分成两半后、对其中一半的数值串再进行调整后的数值串构成,
上述给定的对称型的数值串,其值设定为该数值串的总和非零、该数值串的隔1相加之和为同号且相等。
4.根据权利要求3所述的数字滤波器的设计方法,其特征在于:
上述给定的对称型的数值串,是由-1,0,9,16,9,0,-1的比率构成,
上述非对称型的滤波器系数,其数值串是由-(1-N/8),0,(9-N/8),8、或者8,(9-N/8),0,-(1-N,8)的比率构成,其中N为0≤N≤8。
5.根据权利要求1或2的任一项所述的数字滤波器的设计方法,其特征在于:
构成上述第2单元滤波器的非对称型的滤波器系数,是将给定的对称型的数值串从其中央分成两半后、对其中一半的数值串再进行调整后的数值串构成,
上述给定的对称型的数值串,其值设定为该数值串的总和为零、该数值串的隔1相加之和为异号且绝对值相等。
6.根据权利要求5所述的数字滤波器的设计方法,其特征在于:
上述给定的对称型的数值串,是由1,0,-9,16,-9,0,1的比率构成,
上述非对称型的滤波器系数,其数值串是由(1-N/8),0,-(9-N/8),8、或者8,-(9-N/8),0,(1-N/8)的比率构成,其中N为0≤N≤8。
7.根据权利要求1~6的任一项所述的数字滤波器的设计方法,其特征在于:
通过在与构成上述第1及第2单元滤波器的非对称型的滤波器系数相对应的各抽头间,插入n时钟量的延迟,来调整滤波器的通过频带。
8.一种数字滤波器的设计装置,其特征在于,具备:
第1单元滤波器保持模块,保持与第1单元滤波器相关的信息,该第1单元滤波器,具有将值设定为数值串的总和非零、该数值串的隔1相加之和为同号且相等的滤波器系数,并通过上述滤波器系数、分别将由多个延迟器所形成的带抽头的延迟线中的各抽头的信号变为数倍后,将其相乘结果相加输出;
第2单元滤波器保持模块,保持与第2单元滤波器相关的信息,该第2单元滤波器,具有将值设定为数值串的总和为零、该数值串的隔1相加之和为异号且绝对值相等的滤波器系数,并通过上述滤波器系数、分别将由多个延迟器所形成的带抽头的延迟线中的各抽头的信号变为数倍后,将其相乘结果相加输出;
组合模块,指示0个以上的上述第1单元滤波器和0个以上的上述第2单元滤波器的任意组合;以及,
滤波器系数运算模块,使用由上述第1以及第2单元滤波器保持模块所保持的信息,求取与由上述组合模块所指示的组合相对应的滤波器系数。
9.一种数字滤波器的设计装置,其特征在于,具备:
第1单元滤波器保持模块,保持与第1单元滤波器相关的信息,该第1单元滤波器,通过将具有将值设定为数值串的总和非零、该数值串的隔1相加之和为同号且相等的非对称型的滤波器系数的多个基本滤波器级联连接为整体上的数值串为对称型而形成;
第2单元滤波器保持模块,保持与第2单元滤波器相关的信息,该第2单元滤波器,通过将具有将值设定为数值串的总和为零、该数值串的隔1相加之和为异号且绝对值相等的非对称型的滤波器系数的多个基本滤波器级联连接为整体上的数值串为对称型而形成;
组合模块,指示0个以上的上述第1单元滤波器和0个以上的上述第2单元滤波器的任意组合;以及,
滤波器系数运算模块,使用由上述第1以及第2单元滤波器保持模块所保持的信息,求取与由上述组合模块所指示的组合相对应的滤波器系数。
10.根据权利要求8或9所述的数字滤波器的设计装置,其特征在于:
具备延迟模块,通过在与构成上述第1及第2单元滤波器的非对称型的滤波器系数相对应的各抽头间,插入n时钟量的延迟,来调整滤波器的通过频带。
11.一种数字滤波器,其特征在于:
具备由多个延迟器所形成的带抽头的延迟线,通过用权利要求1~7的任一项所述的滤波器设计方法求得的滤波器系数、分别将各抽头的信号变为数倍后,相加输出。
12.一种数字滤波器,具备由多个延迟器形成的带抽头的延迟线,通过获得的滤波器系数分别将各抽头的信号变为数倍后,相加输出,其特征为:
上述滤波器系数,将值设定为其数值串为非对称型、上述数值串的总和非零、且上述数值串的隔1相加之和为同号且相等。
13.一种数字滤波器,具备由多个延迟器形成的带抽头的延迟线,通过获得的滤波器系数分别将各抽头的信号变为数倍后,相加输出,其特征为:
将具有将值设定为数值串的总和非零、该数值串的隔1相加之和为同号且相等的非对称型的滤波器系数的多个基本滤波器,级联连接为整体上的数值串构成对称型来形成。
14.根据权利要求11或12所述的数字滤波器,其特征在于:
上述非对称型的滤波器系数,是将给定的对称型的数值串从其中央分成两半后、对其中一半的数值串再进行调整后的数值串构成,
上述给定的对称型的数值串,其值设定为该数值串的总和非零、该数值串的隔1相加之和为同号且相等。
15.根据权利要求14所述的数字滤波器,其特征在于,上述给定的对称型的数值串,是由:-1,0,9,16,9,0,-1的比率构成。
16.根据权利要求15所述的数字滤波器,其特征在于,上述非对称型的滤波器系数,其数值串是由-(1-N/8),0,(9-N/8),8、或者8,(9-N/8),0,-(1-N,8)的比率构成,其中N为0≤N≤8。
17.一种数字滤波器,具备由多个延迟器形成的带抽头的延迟线,通过获得的滤波器系数分别将各抽头的信号变为数倍后,相加输出,其特征为:
上述滤波器系数,将值设定为其数值串为非对称型、上述数值串的总和为零、且上述数值串的隔1相加之和为异号且绝对值相等。
18.一种数字滤波器,具备由多个延迟器形成的带抽头的延迟线,通过获得的滤波器系数分别将各抽头的信号变为数倍后,相加输出,其特征为:
将具有将值设定为数值串的总和为零、该数值串的隔1相加之和为异号且绝对值相等的非对称型的滤波器系数的多个基本滤波器,级联连接为整体上的数值串构成对称型来形成。
19.根据权利要求17或18所述的数字滤波器,其特征在于:
上述非对称型的滤波器系数,是将给定的对称型的数值串从其中央分成两半后、对其中一半的数值串再进行调整后的数值串构成,
上述给定的对称型的数值串,其值设定为该数值串的总和为零、该数值串的隔1相加之和为异号且绝对值相等。
20.根据权利要求19所述的数字滤波器,其特征在于:上述给定的对称型的数值串,是由1,0,-9,16,-9,0,1的比率构成。
21.根据权利要求20所述的数字滤波器,其特征在于:上述非对称型的滤波器系数,其数值串是由(1-N/8),0,-(9-N/8),8、或者8,-(9-N/8),0,(1-N/8)的比率构成,其中N为0≤N≤8。
22.根据权利要求12~21的任一项所述的数字滤波器,其特征在于,具备延迟模块,其在与上述非对称型的滤波器系数对应的各抽头间,插入n时钟量的延迟。
23.一种数字滤波器,其特征在于:
将权利要求12或13所述的数字滤波器,或者在与上述非对称型的滤波器系数相对应的各抽头间插入了n时钟量的延迟的数字滤波器,作为基本的单元滤波器,将多段该单元滤波器级联连接而构成。
24.一种数字滤波器,其特征在于:
将权利要求17或18所述的数字滤波器,或者在与上述非对称型的滤波器系数相对应的各抽头间插入了n时钟量的延迟的数字滤波器,作为基本的单元滤波器,将多段该单元滤波器级联连接而构成。
25.一种数字滤波器,其特征在于:
将权利要求12或13所述有数字滤波器,或者在与上述非对称型的滤波器系数相对应的各抽头间插入了n时钟量的延迟的数字滤波器作为基本的低通单元滤波器;同时,将权利要求17或18所述的数字滤波器,或者在与上述非对称型的滤波器系数相对应的各抽头间插入了n时钟量的延迟的数字滤波器作为基本的高通单元滤波器,任意级联连接1以上的上述低通单元滤波器以及1以上的上述高通单元滤波器而构成。
26.根据权利要求25所述的数字滤波器,其特征在于:
将权利要求14所述的上述给定的对称型的数值串作为滤波器系数的数字滤波器,或者在与由该给定的对称型的数值串构成的滤波器系数相对应的各抽头间插入n时钟量的延迟的数字滤波器,作为第2低通单元滤波器,进而任意级联连接该第2低通单元滤波器而构成。
27.根据权利要求25所述的数字滤波器,其特征在于:
将权利要求19所述的上述给定的对称型的数值串作为滤波器系数的数字滤波器,或者在与由该给定的对称型的数值串构成的滤波器系数相对应的各抽头间插入n时钟量的延迟的数字滤波器,作为第2高通单元滤波器,进而任意级联连接该第2高通单元滤波器而构成。
28.根据权利要求22所述的数字滤波器,其特征在于:
上述延迟模块,具备按照周期为基准时钟的1/n倍的第2时钟,将输入的数据保持输出的延迟电路。
29.一种数字滤波器设计用程序,其特征在于,在计算机中执行与权利要求1~7的任一项所述的滤波器设计方法相关的处理顺序。
30.一种数字滤波器设计用程序,其特征在于,令计算机具备权利要求8~10的任一项所述的各模块的功能。
CNA2003801018286A 2002-10-21 2003-10-20 数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器 Pending CN1706100A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP306315/2002 2002-10-21
JP2002306315 2002-10-21

Publications (1)

Publication Number Publication Date
CN1706100A true CN1706100A (zh) 2005-12-07

Family

ID=32105204

Family Applications (1)

Application Number Title Priority Date Filing Date
CNA2003801018286A Pending CN1706100A (zh) 2002-10-21 2003-10-20 数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器

Country Status (6)

Country Link
EP (1) EP1557946A1 (zh)
JP (1) JPWO2004036746A1 (zh)
KR (1) KR20050071604A (zh)
CN (1) CN1706100A (zh)
TW (1) TW200408813A (zh)
WO (1) WO2004036746A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101884227A (zh) * 2006-04-03 2010-11-10 Srs实验室有限公司 音频信号处理
CN101968964A (zh) * 2010-08-20 2011-02-09 北京中星微电子有限公司 一种去除语音信号中直流分量的方法及装置
US8027477B2 (en) 2005-09-13 2011-09-27 Srs Labs, Inc. Systems and methods for audio processing

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5020265B2 (ja) 2007-02-07 2012-09-05 三菱電機株式会社 乗員保護装置
TWI597939B (zh) 2009-02-18 2017-09-01 杜比國際公司 具相位偏移之複數值合成濾波器組
EP4233282A1 (en) * 2020-10-23 2023-08-30 Marvell Asia Pte, Ltd. Equalization in high-speed data channel having sparse impulse response

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE467680B (sv) * 1990-12-19 1992-08-24 Johan Hellgren Digital filterbank med minskad effektfoerbrukning
JPH05243908A (ja) * 1991-12-06 1993-09-21 Nec Corp ろ波器
JPH05235701A (ja) * 1992-02-25 1993-09-10 Nippon Steel Corp 環状畳み込みによるディジタルフィルタバンク処理方法及び装置
GB9226536D0 (en) * 1992-12-21 1993-02-17 Unilever Plc Foodstuffs and other compositions
JPH1079686A (ja) * 1996-09-02 1998-03-24 Ricoh Co Ltd ディジタル相関器

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8027477B2 (en) 2005-09-13 2011-09-27 Srs Labs, Inc. Systems and methods for audio processing
US9232319B2 (en) 2005-09-13 2016-01-05 Dts Llc Systems and methods for audio processing
CN101884227A (zh) * 2006-04-03 2010-11-10 Srs实验室有限公司 音频信号处理
CN101884227B (zh) * 2006-04-03 2014-03-26 Dts有限责任公司 音频信号处理
US8831254B2 (en) 2006-04-03 2014-09-09 Dts Llc Audio signal processing
CN101968964A (zh) * 2010-08-20 2011-02-09 北京中星微电子有限公司 一种去除语音信号中直流分量的方法及装置
CN101968964B (zh) * 2010-08-20 2015-09-02 北京中星微电子有限公司 一种去除语音信号中直流分量的方法及装置

Also Published As

Publication number Publication date
WO2004036746A1 (ja) 2004-04-29
JPWO2004036746A1 (ja) 2006-02-16
KR20050071604A (ko) 2005-07-07
TWI292044B (zh) 2008-01-01
EP1557946A1 (en) 2005-07-27
TW200408813A (en) 2004-06-01

Similar Documents

Publication Publication Date Title
US5982305A (en) Sample rate converter
CN1025398C (zh) 使用单个位及多位量化的多阶σ-δ模/数转换器
CN1706101A (zh) 数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器
CN1051415C (zh) 把1位数字数据信号转化为模拟数据的信号处理电路
CN1925323A (zh) 采样率转换方法及其电路
TW200842827A (en) Method and system for utilizing rate conversion filters to reduce mixing complexity during multipath multi-rate audio processing
CN1707954A (zh) 多线型平行处理三角积分模拟/数字转换器
CN1684405A (zh) 时钟同步器以及时钟与数据恢复装置和方法
CN1251411C (zh) 锁相环电路、时钟生成电路和时钟生成方法
CN1977455A (zh) 无抖动采样率转换
CN1835394A (zh) D类放大器
CN101060312A (zh) 音量变换装置
CN1706100A (zh) 数字滤波器的设计方法及装置、数字滤波器设计用程序、数字滤波器
CN1554147A (zh) 数字滤波器及其设计方法
US20080084343A1 (en) Methods and systems for implementing a digital-to-analog converter
CN1113469C (zh) 信号处理方法及其装置
CN1825760A (zh) 数据传输控制和采样频率转换器
CN1210871C (zh) 信号传输设备
CN1237719C (zh) 消除控制电路接收的信号中的自抖动的锁相环电路和方法
CN1792033A (zh) 数字滤波器及该设计方法、设计装置、数字滤波器设计用程序
CN1297981C (zh) 校正偏移的设备和方法
CN1442986A (zh) I/q解调器及其i/q信号抽样方法
CN101064502A (zh) 带有下采样功能的数字信号滤波装置及方法
CN1146905C (zh) 波形均衡器
US6486813B1 (en) Oversampling circuit digital/analog converter

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C02 Deemed withdrawal of patent application after publication (patent law 2001)
WD01 Invention patent application deemed withdrawn after publication