CN1407720A - 可变抽样率的递归数字滤波器 - Google Patents
可变抽样率的递归数字滤波器 Download PDFInfo
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/0294—Variable filters; Programmable filters
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
- H03H17/00—Networks using digital techniques
- H03H17/02—Frequency selective networks
- H03H17/04—Recursive filters
- H03H17/0416—Recursive 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所述的方法,其特征在于:所述修改步骤包括把所述初始系数乘以所述系数因子而产生所述新系数的步骤。
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)
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)
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)
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 |
-
2001
- 2001-08-08 US US09/925,546 patent/US6976045B2/en not_active Expired - Fee Related
-
2002
- 2002-07-16 DE DE60230372T patent/DE60230372D1/de not_active Expired - Lifetime
- 2002-07-16 EP EP02254993A patent/EP1289143B1/en not_active Expired - Lifetime
- 2002-07-30 JP JP2002221582A patent/JP4159318B2/ja not_active Expired - Fee Related
- 2002-08-08 CN CN02129721A patent/CN100592627C/zh not_active Expired - Fee Related
Cited By (6)
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 |