CN1407720A - 可变抽样率的递归数字滤波器 - Google Patents

可变抽样率的递归数字滤波器 Download PDF

Info

Publication number
CN1407720A
CN1407720A CN02129721.5A CN02129721A CN1407720A CN 1407720 A CN1407720 A CN 1407720A CN 02129721 A CN02129721 A CN 02129721A CN 1407720 A CN1407720 A CN 1407720A
Authority
CN
China
Prior art keywords
coefficient
digital filter
factor
new
ratio
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.)
Granted
Application number
CN02129721.5A
Other languages
English (en)
Other versions
CN100592627C (zh
Inventor
K·M·费尔古森
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.)
Tektronix Inc
Original Assignee
Tektronix Inc
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 Tektronix Inc filed Critical Tektronix Inc
Publication of CN1407720A publication Critical patent/CN1407720A/zh
Application granted granted Critical
Publication of CN100592627C publication Critical patent/CN100592627C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related 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/0294Variable filters; Programmable filters
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/04Recursive filters
    • H03H17/0416Recursive filters with input-sampling frequency and output-delivery frequency which differ, e.g. extrapolation; Anti-aliasing

Landscapes

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

Abstract

本发明提供可变抽样率的递归数字滤波器,它根据各个相当的抽样率改变自适应数字滤波器的滤波系数。最初针对特定的抽样率来校准自适应数字滤波器,所述特定的抽样率确定要采用的特定的滤波系数。当改变抽样率时,确定经校准的或者初始的抽样率与新抽样率之比,并根据该比值计算系数因子。系数因子用于修改对应于经校准的抽样率的初始滤波系数,从而产生对应于新抽样率的新系数。把新系数应用于自适应数字滤波器以提供所需的信号处理而不需要再抽样或再校准,而且不改变频率响应。

Description

可变抽样率的递归数字滤波器
技术领域
本发明涉及数字滤波,更具体地说,涉及在一定的抽样率范围上维持恒定的频率响应的可变抽样率的递归数字滤波器。
技术背景
诸如抽样率转换、抖动/漂移补偿和/或测量、图像质量测量等若干信号处理应用所共有的问题是设法在一定的抽样率范围上维持恒定的频率响应。在基于人类视觉模型的图像质量测量中,其中用不同的格式、显示器(帧频、空间分辨率)和观看距离(在人眼的视场内不同的空间抽样率)来观看测试和参考视频序列,保持效率和精确度是重要的。如2001年5月16日提交的序列号为09/858775的共同未决的美国专利申请中所描述的,维持自适应滤波器的空间和时间频率响应,使人们可以模拟接近显示器(例如2倍的屏幕高度)和远离(例如10倍的屏幕高度)时观看视频序列之间的视觉差异,而不需要例如再抽样。
现有方法需要人类视觉模型的再抽样或再校准。校准可能是非常耗时的,所以不希望对每个有差异的观看距离进行校准。再抽样涉及到对数据再抽样两次,一次是针对人类视觉模型的输入,然后是针对中间人类视觉模型图像,这刚好在同时涉及参考和测试视频的计算(求差)之前。这需要额外的处理和引入内插不准确度。
所需的是这样一种可变抽样率的递归数字滤波器,它在一定的抽样率范围内维持恒定的频率响应,而不需要对每个新的抽样率再抽样或再校准。
发明内容
因此,本发明提供可变抽样率的递归数字滤波器,它根据各个相当的抽样率改变自适应数字滤波器的滤波系数。最初针对特定的抽样率来校准自适应数字滤波器,所述特定的抽样率确定要采用的特定的滤波系数。当改变抽样率时,确定经校准的或者初始的抽样率与新抽样率之比,并根据该比值计算系数因子。系数因子用于修改对应于经校准的抽样率的初始滤波系数,从而产生对应于新抽样率的新系数。把新系数应用于自适应数字滤波器,以便提供所需的信号处理而不需要再抽样或再校准,而且不改变频率响应。
结合所附权利要求书和附图来阅读以下详细描述,本发明的目的、优点及其它新颖的特征是显而易见的。
附图说明
图1是结合了根据本发明的可变抽样率递归数字滤波器的基于人类视觉模型的图像质量测量装置的框图。
图2是用于根据本发明的可变抽样率递归数字滤波器的可变相当抽样率系数转换器的框图。
图3是说明根据本发明的可变抽样率递归数字滤波器在不同抽样率下的恒定频率响应的曲线图。
具体实施方式
参照图1,作为利用可变抽样率递归数字滤波器的信号处理装置的实例,所示的基于人类视觉模型的图像质量测量装置10具有作为输入的参考视频序列和测试(退化的)视频序列。把视频序列存储并显示在各个显示装置12、14上,而且由人类视觉模型16来处理。最初由如上文引用的未决美国专利申请所描述的自适应数字滤波器18、20处理各个通道、即参考和测试通道,然后才由其余的模型元件22处理,从而提供图像质量测量输出。每个数字滤波器18、20具有相当观看距离或空间抽样率和帧频或时间抽样率控制器24、26,它们向数字滤波器提供可变系数,从而独立地将其调到特定的观看距离和/或帧抽样率。
为了给出对输入和输出进行再抽样的等效响应,用于数字滤波器18、20的控制信号或新系数是这样导出的:取在各个抽样率的双线性变换之比。按照常规,双线性变换B把z-域(空间)中的滤波器极点转换到s-域(频率),其中抽样率变化是线性的。将S-域中的极点乘以抽样率之比,然后又转换回z-域。
B{Pz}→Ps; B-1{Ps *R}→Pz′=K*Pz其中B是双线性变换,Pz是z-域中的滤波器极点,Ps是变换到s-域中的滤波器极点,R是初始抽样率和新抽样率之比,K是常数或z因子,而Pz′是变换回z-域中的新滤波器极点。
给定初始或经校准的抽样率以及相应的滤波器系数,就可通过把初始滤波器系数乘以如下确定的K或z因子而得出对应于新抽样率的新滤波器系数:
z因子(z,R):=(1/z){(z(1+R)+(1-R))/(z(1-R)+(1+R))}其中z=初始z域系数,而R=抽样率之比(旧/新)。此等式如图2所示,可以用硬件或软件来实现。
考虑以下形式的滤波器的递归一阶差分方程:
yn+1:=yn+c(xn-yn); yn=yn-1+c(xn-1-yn-1)=(1-c)yn-1+cxn-1
方程的z-变换是:
H(z,c):=c/(z-1+c)
如果直流增益保持为1:c/(1-1+c)=1;则z=1且H(z,c)=1,而滤波器在(1-c)具有极点,也就是说,对于c=0.3,极点=0.7。
频率响应由下式给定:
F(c,ω):=H(e-jω,c)和ωn=π(n/N),其中n/N在0到1之间变化。
以下是在基于双线性变换的时间再抽样之后匹配滤波器合成频率响应的实例。所示滤波器是由两个低通滤波器之间的差异产生的带通滤波器。如果初始抽样率T0是1/60,而各个低通滤波器系数是Cts0=0.6和Ctc0=0.48,则当抽样率变为T1=1/24,则比值R变为0.4,即T0/T1=(1/60)/(1/24)=24/60=0.4。于是,新系数变为Cts1:=z因子(Cts0,R)*Cts0,Ctc1:=z因子(Ctc0,R)*Ctc0,即Cts1=0.818而Ctc1=0.754。新抽样率ω1=ω0 *R。则带通滤波器的频率响应由对应于初始抽样率的Ft0(ω):=|F(Ctc0,ω)2-F(Cts0,ω)2|和对应于新抽样率的Ft1(ω):=|F(Ctc1,ω)2-F(Cts1,ω)2|给出。当适当地缩放并画在一起时,如图3所示,其中实线是组合滤波器在初始抽样率的频率响应,而虚线是组合滤波器在新抽样率下的频率响应,显然,利用双线性变换为滤波器提供了所需的恒定频率响应结果。
因此,本发明利用双线性变换来确定作为初始系数和初始抽样率与新抽样率之比的函数的系数因子,从而提供了可变抽样率递归数字滤波器,所述系数因子被用来针对新抽样率调整初始系数,以便当应用于可变抽样率递归数字滤波器时,在不同抽样率下的频率响应是恒定的。

Claims (8)

1.一种可变抽样率递归数字滤波器,它包括:
自适应数字滤波器,它具有作为输入的初始系数,该系数定义自适应数字滤波器如何在初始抽样率下处理数字信号;以及
为所述自适应数字滤波器确定作为所述初始系数和所述初始抽样率与新抽样率之比的函数的新系数的装置,当在所述新抽样率下处理所述数字信号时,使得所述自适应数字滤波器的所述频率响应在两个抽样率下都是恒定的。
2.如权利要求1所述的可变抽样率递归数字滤波器,其特征在于所述确定装置包括:
计算作为所述初始系数与所述比值的函数的系数因子的装置,以及
用所述系数因子来修改所述初始系数以产生所述新系数的装置。
3.如权利要求2所述的可变抽样率递归数字滤波器,其特征在于所述计算装置执行以下公式:
z因子(z,R):=(1/z){(z(1+R)+(1-R))/(z(1-R)+(1+R))}其中z因子(z,R)是所述系数因子,z是所述初始系数,而R是所述比值。
4.如权利要求2或3所述的可变抽样率递归数字滤波器,其特征在于:所述修改装置包括用于把所述初始系数乘以所述系数因子而产生所述新系数的装置。
5.一种对数字信号进行可变抽样率递归数字滤波的方法,它包括以下步骤:
为处理所述数字信号的自适应数字滤波器确定初始抽样率和相应的初始系数;
确定用于处理所述数字信号的、作为所述初始系数和所述初始抽样率与新抽样率之比的函数的新系数;以及
把所述新系数应用于所述自适应数字滤波器,使得所述自适应数字滤波器在所述两个抽样率下具有恒定的频率响应。
6.如权利要求5所述的方法,其特征在于所述新系数确定步骤包括以下步骤:
计算作为所述初始系数和所述比值的函数的系数因子;以及
用所述系数因子修改所述初始系数以产生所述新系数。
7.如权利要求6所述的方法,其特征在于所述计算步骤包括以下步骤:
求出以下公式:
z因子(z,R):=(1/z){(z(1+R)+(1-R))/(z(1-R)+(1+R))}其中z因子(z,R)是所述系数因子,z是所述初始系数,而R是所述比值。
8.如权利要求6或7所述的方法,其特征在于:所述修改步骤包括把所述初始系数乘以所述系数因子而产生所述新系数的步骤。
CN02129721A 2001-08-08 2002-08-08 可变抽样率的递归数字滤波器及对应的方法 Expired - Fee Related CN100592627C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US09/925546 2001-08-08
US09/925,546 US6976045B2 (en) 2001-08-08 2001-08-08 Variable sample rate recursive digital filter

Publications (2)

Publication Number Publication Date
CN1407720A true CN1407720A (zh) 2003-04-02
CN100592627C CN100592627C (zh) 2010-02-24

Family

ID=25451883

Family Applications (1)

Application Number Title Priority Date Filing Date
CN02129721A Expired - Fee Related CN100592627C (zh) 2001-08-08 2002-08-08 可变抽样率的递归数字滤波器及对应的方法

Country Status (5)

Country Link
US (1) US6976045B2 (zh)
EP (1) EP1289143B1 (zh)
JP (1) JP4159318B2 (zh)
CN (1) CN100592627C (zh)
DE (1) DE60230372D1 (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101258682B (zh) * 2005-07-07 2011-05-18 德克萨斯仪器股份有限公司 自动输入误差恢复电路和用于递归数字滤波器的方法
CN102647168A (zh) * 2010-10-28 2012-08-22 矽统科技股份有限公司 动态滤波装置及其方法
CN102694528A (zh) * 2011-03-23 2012-09-26 美国亚德诺半导体公司 异步采样率转换器中抽取率的自适应控制方法和装置
CN101800845B (zh) * 2010-01-06 2013-05-22 北京中科大洋科技发展股份有限公司 一种可调节的基于光滑曲线的频域滤波器设计方法

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7158591B2 (en) * 2001-05-09 2007-01-02 Signum Concept, Inc. Recursive resampling digital filter structure for demodulating 3G wireless signals
US7405747B2 (en) * 2004-04-14 2008-07-29 Tektronix, Inc. Extended K2T measurement of video
JP4563310B2 (ja) * 2005-12-01 2010-10-13 パナソニック株式会社 無線受信機
CN100463367C (zh) * 2007-03-12 2009-02-18 北京航空航天大学 一种章动频率自动跟踪的高通数字滤波方法
US8355567B2 (en) 2008-12-12 2013-01-15 Tektronix, Inc. Method and apparatus for implementing moving image color appearance model for video quality ratings prediction
US20130103283A1 (en) * 2011-10-19 2013-04-25 Delphi Technologies, Inc. Digital filter for an event driven engine controller having a variable sample rate
US8917932B2 (en) 2011-12-29 2014-12-23 Tektronix, Inc. Method of predicting perceptual differences between memory and directly viewed colors in video
JP2016511588A (ja) * 2013-02-19 2016-04-14 テクトロニクス・インコーポレイテッドTektronix,Inc. 画像アピアランス・フレームワーク及びデジタル画像制作及び表示用アプリケーション
WO2014179489A1 (en) * 2013-05-01 2014-11-06 Starkey Laboratories, Inc. Adaptive feedback cancellation coefficients based on voltage
CN103873016B (zh) * 2014-02-24 2016-09-21 北京遥测技术研究所 一种任意采样率变换数字重构滤波器系数设计方法
US10045050B2 (en) * 2014-04-25 2018-08-07 Vid Scale, Inc. Perceptual preprocessing filter for viewing-conditions-aware video coding
DE102017212715B3 (de) * 2017-07-25 2019-01-31 Robert Bosch Gmbh Verfahren zum Verarbeiten von kontinuierlichen Sensorsignalen und Sensorsystem
CN112865753B (zh) * 2020-12-28 2023-04-07 珠海格力电器股份有限公司 一种滤波系数调整方法、装置、存储介质及滤波器

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
NL8503478A (nl) 1985-12-18 1987-07-16 Philips Nv Coefficienten generator met fase vergrendelde lus voor een filterinrichting met niet-rationele verhouding tussen ingangs- en uitgangsbemonsterfrequentie.
JPH07307629A (ja) * 1994-05-11 1995-11-21 Asahi Kasei Micro Syst Kk トーンコントロール回路
US5928313A (en) * 1997-05-05 1999-07-27 Apple Computer, Inc. Method and apparatus for sample rate conversion
US6134570A (en) * 1997-05-12 2000-10-17 Level One Communications, Inc. Method and apparatus for efficient implementation of a multirate LMS filter
KR100244767B1 (ko) * 1997-06-25 2000-02-15 전주범 디지탈 자기 기록/재생 시스템의 선택적 동기/비동기 부분 응답 채널 데이터 검출 장치
US6310919B1 (en) * 1998-05-07 2001-10-30 Sarnoff Corporation Method and apparatus for adaptively scaling motion vector information in an information stream decoder
DE19829289C2 (de) 1998-06-30 2001-12-06 Siemens Ag Verfahren zur Berechnung der Koeffizienten eines nichtrekursiven digitalen Filters
US6531970B2 (en) * 2001-06-07 2003-03-11 Analog Devices, Inc. Digital sample rate converters having matched group delay

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101258682B (zh) * 2005-07-07 2011-05-18 德克萨斯仪器股份有限公司 自动输入误差恢复电路和用于递归数字滤波器的方法
CN101800845B (zh) * 2010-01-06 2013-05-22 北京中科大洋科技发展股份有限公司 一种可调节的基于光滑曲线的频域滤波器设计方法
CN102647168A (zh) * 2010-10-28 2012-08-22 矽统科技股份有限公司 动态滤波装置及其方法
CN102647168B (zh) * 2010-10-28 2015-01-14 矽统科技股份有限公司 动态滤波装置及其方法
CN102694528A (zh) * 2011-03-23 2012-09-26 美国亚德诺半导体公司 异步采样率转换器中抽取率的自适应控制方法和装置
CN102694528B (zh) * 2011-03-23 2015-05-06 美国亚德诺半导体公司 异步采样率转换器中抽取率的自适应控制方法和装置

Also Published As

Publication number Publication date
EP1289143B1 (en) 2008-12-17
JP2003133909A (ja) 2003-05-09
EP1289143A2 (en) 2003-03-05
EP1289143A3 (en) 2005-05-04
DE60230372D1 (de) 2009-01-29
JP4159318B2 (ja) 2008-10-01
US6976045B2 (en) 2005-12-13
US20030031281A1 (en) 2003-02-13
CN100592627C (zh) 2010-02-24

Similar Documents

Publication Publication Date Title
CN100592627C (zh) 可变抽样率的递归数字滤波器及对应的方法
Moller et al. Evaluation and design of filters using a Taylor series expansion
US20070094317A1 (en) Method and system for B-spline interpolation of a one-dimensional signal using a fractional interpolation ratio
Valimaki et al. Fractional delay filter design based on truncated Lagrange interpolation
US5473555A (en) Method and apparatus for enhancing frequency domain analysis
RU2628219C2 (ru) Вещественное-в-комплексное преобразование с малой задержкой в банках фильтров с перекрытием для частично комплексной обработки
JP2005519269A (ja) デジタル周波数応答補正器及び随意応答発生器の方式
US20070053420A1 (en) Method, apparatus, and program for designing digital filters
US7542053B2 (en) Methods and apparatus for re-scaling image by variable re-scaling factor
US6748409B1 (en) Data interpolating system
CN1112802C (zh) 图像格式转换器及其方法
US6766338B1 (en) High order lagrange sample rate conversion using tables for improved efficiency
JP2002199353A (ja) 情報信号処理装置、情報信号処理方法、画像信号処理装置およびそれを使用した画像表示装置、並びに情報提供媒体
JP2003116103A (ja) ディジタル映像処理装置及び方法
AU2007219336B8 (en) Method and apparatus for abitrary ratio image reduction
US7446820B2 (en) Methods and apparatus for re-scaling image by variable re-scaling factor using B-spline interpolator
JPH06294830A (ja) 周波数領域の解析を強化するための方法と装置
Hsue et al. Blind deconvolution of symmetric noncausal impulse responses using two-sided linear prediction
US20050156625A1 (en) Asynchronus sampling rate conversion
CN1901617A (zh) 使输出图像质量最优化的自适应图像缩放器及其方法
WO2006106713A1 (ja) 最悪値誤差尺度を最小化する信号近似システム
JP3224381B2 (ja) データサンプル間の補間方法及び装置
JPH09214747A (ja) デジタル鮮鋭化/ぼかしフィルタの生成方法
US20040254968A1 (en) Interpolation for waveform and vector displays
CN118136042B (zh) 基于iir频谱拟合的频谱优化方法、系统、终端及介质

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20100224

Termination date: 20170808

CF01 Termination of patent right due to non-payment of annual fee