CN100379154C - 非均匀性取样的限频信号的重建方法 - Google Patents

非均匀性取样的限频信号的重建方法 Download PDF

Info

Publication number
CN100379154C
CN100379154C CNB018167640A CN01816764A CN100379154C CN 100379154 C CN100379154 C CN 100379154C CN B018167640 A CNB018167640 A CN B018167640A CN 01816764 A CN01816764 A CN 01816764A CN 100379154 C CN100379154 C CN 100379154C
Authority
CN
China
Prior art keywords
centerdot
omega
equation
sampling
sigma
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.)
Expired - Fee Related
Application number
CNB018167640A
Other languages
English (en)
Other versions
CN1468433A (zh
Inventor
H·强森
P·罗温伯格
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.)
Infineon Technologies AG
Original Assignee
Infineon Technologies AG
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 Infineon Technologies AG filed Critical Infineon Technologies AG
Publication of CN1468433A publication Critical patent/CN1468433A/zh
Application granted granted Critical
Publication of CN100379154C publication Critical patent/CN100379154C/zh
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/06Continuously compensating for, or preventing, undesired influence of physical parameters
    • H03M1/0617Continuously compensating for, or preventing, undesired influence of physical parameters characterised by the use of methods or means not specific to a particular type of detrimental influence
    • H03M1/0626Continuously compensating for, or preventing, undesired influence of physical parameters characterised by the use of methods or means not specific to a particular type of detrimental influence by filtering
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/06Continuously compensating for, or preventing, undesired influence of physical parameters
    • H03M1/08Continuously compensating for, or preventing, undesired influence of physical parameters of noise
    • H03M1/0836Continuously compensating for, or preventing, undesired influence of physical parameters of noise of phase error, e.g. jitter
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03MCODING; DECODING; CODE CONVERSION IN GENERAL
    • H03M1/00Analogue/digital conversion; Digital/analogue conversion
    • H03M1/12Analogue/digital converters
    • H03M1/1205Multiplexed conversion systems
    • H03M1/121Interleaved, i.e. using multiple converters or converter parts for one channel
    • H03M1/1215Interleaved, i.e. using multiple converters or converter parts for one channel using time-division multiplexing

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Analogue/Digital Conversion (AREA)
  • Complex Calculations (AREA)
  • Superconductors And Manufacturing Methods Therefor (AREA)
  • Glass Compositions (AREA)
  • Transition And Organic Metals Composition Catalysts For Addition Polymerization (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
  • Wire Bonding (AREA)
  • Measurement Of Radiation (AREA)

Abstract

本发明涉及一种非均匀性取样的限频模拟信号xa(t)的重建方法与装置,该非均匀性取样信号包括N个子序列xk(m),k=0,1,...,N-1,N≥2,其是以1/(MT)的取样速率依据xk(m)=xa(nMT+tk)取样而获得,其中M为整数,而tk=kMT/N+Δtk,Δtk不为零。本发明包括自该N个子序列xk(m)形成新序列y(n),以使y(n)至少含有与x(n)=xa(nt)相同的信息,即,xa(t)在一低于ω0的频率区域中以1/T的取样速率取样,ω0是预定限制频率,这通过下列步骤完成:(i)以因子M上取样该N个子序列xk(m)的每一项,k=0,1,...,N-1,而M是正整数;(ii)以对应的数字滤波器过滤该上取样的N个子序列xk(m)的每一序列,k=0,1,...,N-1;及(iii)添加该N个数字滤波子序列以形成y(n)。对应数字滤波器较佳为分段延迟滤波器,且具有一在频带|ωT|≥ω0T中的频率响应Gk=ake(-jωsT),k=0,1,...,N-1,ak是常数,且s等于d+tk,而d为整数。

Description

非均匀性取样的限频信号的重建方法
技术领域
本发明大体上涉及取样领域,更特别地,涉及非均匀性取样的限频信号的重建方法与装置,涉及时间交错的模数转换器(ADC)中的时间偏移的补偿方法与装置,及关于一种用于执行所述重建方法的计算机程序产品。
相关技术与发明背景说明
均匀性取样时,通过在t=nT,-∞<n<∞,等距离取样一模拟信号xa(t)而获得一序列x(n),即,x(n)=xa(nt),T是取样周期,如图1a所示。在此情况下,二连续取样间隔之间的时间恒为T。另一方面,在非均匀性取样时,二连续取样间隔之间的时间依取样间隔而定。本发明处理的情况是其中样本可以分为N个子序列xk(m),k=0,1,2...,N-1,其中xk(m)通过在t=nMT+tk以1/(MT)的取样速率做取样xa(t)而获得,即,xk(m)=xa(nMT+tk),M是正整数。这一取样方式以N=2且M=2而示于图1b。这样的非均匀性取样信号例如发生在时间偏移误差所导致的时间交错的模数转换器(ADC)。
产生的问题是如何自xk(m)形成新序列y(n),而使得y(n)恰等于或大约(某种程度地)等于x(n)。对于传统的时间交错的模数转换器而言,N=M,且理想上,tk=kT。在此情况下,简单地通过使xk(m)交错就可获得y(n)=x(n)。然而实际上,时间偏移误差导致tk不恰等于kT,时间偏移误差将频叠(aliasing)成分引入Y(ejωT),Y(ejωT)是y(n)的傅立叶转换。这意味着y(n)≠x(n),因此y(n)中的信息不再与x(n)中的信息相同。
应注意,众所周知,如果各tk不同,从而使全部样本在时间上是分离的,则xa(t)由xk(m)中的样本所唯一地确定。同样众所周知的是如何使用模拟内插功能而自xk(m)得到xa(t)。然而在实际应用时,这些功能不容易(如果有丝毫可能性的话)达成,于是需要其它的解决方案。
发明概述
因此,本发明的一个目的是分别提供非均匀性取样的限频模拟信号xa(t)的重建方法与非均匀性取样的限频模拟信号xa(t)的重建装置的运作方法,该非均匀性取样信号包括N个子序列xk(m),k=0,1,...,N-1,N≥2,其是以1/(MT)的取样速率依据xk(m)=xa(nMT+tk)经由取样而获得,其中M为正整数,T为采样周期,而tk=kMT/N+Δtk,Δtk不为零,其可以自该N个子序列xk(m)形成新序列y(n),而使得y(n)至少含有与x(n)=xa(nt)相同的信息,即,xa(t)在一低于ω0(可能包含ω0)的频率区域中以1/T的取样速率取样,ω0是预定限制频率。
本发明的又一目的是分别提供这样的方法与装置的运作方法,其有效、快速、简单且成本低。
本发明的再一目的是分别提供此方法与装置的运作方法,其可以减少噪声,例如,数字化噪声。
除其它目的外,上述目的分别通过一方法与一装置的运作方法而获得,所述方法与装置的运作方法执行下列步骤:
(i)以因子M上取样(upsampling)N个子序列xk(m)的每一项,k=0,1,...,N-1;
(ii)以一对应的数字滤波器过滤上取样的N个子序列xk(m)的每一项,k=0,1,...,N-1;及
(iii)添加N个数字滤波子序列以形成y(n)。
较佳地,该对应的数字滤波器是分段延迟(fractional delay)滤波器,且具有一在频带|ωT|≤ω0T中的频率响应Gk(ejωT)=akce-jω(t k +dT),k=0,1,...,N-1,ak为常数,d为整数。
如果ω0T是小于π的固定值,从而原始模拟信号包括一比ω0更高的频率的频率成分,则取得区域性最佳重建,即,y(n)含有与x(n)=xa(nT)相同的信息,即,只有在频率区域|ω|≤ω0中以1/T的取样速率取样的xa(t)。区域性最佳重建特别重视过取样的(oversampled)系统,其中较低频成分载有基本信息,而较高频成分含有要由数字和/或模拟滤波器除去的所不希望的成分(例如噪声)。
此处,分段延迟滤波器具有在频带ω0T<|ωT|≤π中的频率响应Gk=akAk(ejωT),k=0,1,...,N-1,其中Ak(ejωT)是任意复数函数。
另一方面,如果ω0不包含原始模拟信号的频率成分(即,ω0T包含直到π为止的全部频率),则取得最佳重建,即,y(n)等于x(n)。
在任一情况下,都产生二种不同的情形:(1)2K0+1=N与(2)2K0+1<N,其中K0得自于
Figure C0181676400101
以用于区域性最佳重建,其中分别地,
Figure C0181676400102
应读作是大于或等于 M ( ω 0 T + ω 1 T ) 2 π 的最小整数,而[-ω1,ω1]是所述限频模拟信号xa(t)所处的频带,且通过
K0=M-1
而用于最佳重建。
在状况(1),ak计算如下
a=B-1c,
a是ak的向量形式,且由下式计算得出
a=[a0a1...aN-1]T
B-1是B的逆矩阵,且B得自于
B = u 0 - K 0 u 1 - K 0 . . . u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) . . . u N - 1 - ( K 0 - 1 ) . . . u 0 K 0 u 1 K 0 . . . u N - 1 K 0
其中 u k = e - j 2 π MT t k
而c是
c=[c0c1...c2K0]T
其中
c k = M , k = K 0 0 , k = 0,1 , . . . , 2 K 0 , k ≠ K 0
在状况(2),各ak计算如下
a = B ^ - 1 c ^
a定义为
a=[au afix]T
其中au与afix含有(2K0+1)个未知数ak与L=N-2K0-1个固定常数ak
Figure C0181676400111
Figure C0181676400112
的逆矩阵,
Figure C0181676400113
得自于
B ^ = B S
其中B得自于
B = u 0 - K 0 u 1 - K 0 . . . u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) . . . u N - 1 - ( K 0 - 1 ) . . . u 0 K 0 u 1 K 0 . . . u N - 1 K 0
其中 u k = e - j 2 π MT t k
S得自于
S=[SzSd],其中
S Z = 0 0 . . . 0 0 0 . . . 0 . . . 0 0 . . . 0
并且
Sd=diag[11...1]
Figure C0181676400118
c ^ = [ c a fix ] T
其中c得自于
c=[c0c1...c2K0]T
其中
c k = M , k = K 0 0 , k = 0,1 , . . . , 2 K 0 , k ≠ K 0
故L个ak可任意地选择。优选将其选择为零(在此情况下对应的频道被除去)或者是M/N(在此情况下任何数字化噪声可减至最小)。
本发明的其它目的是提供一种在时间交错的模数转换器(ADC)系统的运作方法——其包括多个模数转换器(ADC)——中的时间偏移的补偿方法,以及提供模数转换器系统本身。
因此,所提供的这样的方法与模数转换器系统运作方法分别包括了上述的方法与装置的运作方法,其中N个子序列xk(m),k=0,1,2...,N-1,N≥2中的每一序列由一个对应的模数转换器取样。
本发明的再一目的是提供一种计算机程序产品,其用于非均匀性取样的限频模拟信号的重建。
此目的是通过一计算机程序产品而实现,该计算机程序产品可加载到一数字信号处理装置的内部存储器,包括软件码部分,其在该产品于该装置上运行的时候,用于执行上述任何方法。
本发明的一优点为,可以产生完全或部分重建的数字信号,不需要应用很复杂且难以实行的模拟内插功能。
通过本发明的实施例的下列详细说明,将可明白本发明的其它特征及其优点。
附图简单说明
通过以下所提供的本发明较佳实施例的详细说明与附图1~6,将可更完整地了解本发明,所给附图仅用于说明,因此并不会限制本发明。
图1a示意说明均匀性取样,其中通过在t=nT,-∞<n<∞,等距离取样一模拟信号xa(t)而获得一序列x(n),即,x(n)=xa(nt);图1b示意说明非均匀性取样,其中样本分为二个子序列xk(m),k=0,1,而xk(m)通过在t=n2T+tk以1/(2T)的取样速率取样xa(t)而获得,即,xk(m)=xa(n2T+tk)。
图2示意说明一均匀性取样器与数字转换器。
图3示意说明一上取样器。
图4示意说明一混合模拟/数字滤波器组模数转换器系统。
图5示意说明一解析滤波器组系统,其用于产生xk(m),k=0,1,2...,N-1,其中xk(m)是在时间间隔t=nMT+tk取样xa(t)而获得的N个子序列。
图6示意说明图4的系统中的上取样与合成组的多相代表。
实施例详细说明
在下列说明中,为了解释而非限制,提出了特定细节,以供完整了解本发明。然而,对本领域技术人员显而易见的是,本发明可通过不同于这些特定细节的其它变例而实行。在其它情况下,略去公知的方法与装置的详细说明,以免让不需要的细节混淆本发明的说明。
本发明考虑非均匀性取样的限频信号的重建问题。此问题例如发生在时间偏移误差所导致的时间交错的模数转换器(ADC)。为了精确起见,我们处理下列状况:已知有N个子序列xk(m),k=0,1,2...,N-1,其是依据xk(m)=xa(nMT+tk),以1/(MT)的取样速率来取样限频模拟信号xa(t)而获得,如何自xk(m)形成新序列y(n),以使y(n)恰等于或大约(某种程度地)等于x(n)=xa(nT),即,以1/T的取样速率取样的xa(t)。为此目的,我们在本发明中建议使用N频道数字合成滤波器组。整个系统可以看成传统时间交错的模数转换器的概括化,而前者是将其简化而成的特殊状况。我们证明,使用适当的理想合成滤波器,所建议的系统可以达成y(n)=x(n)。然而,这些合成滤波器不适于以实际的数字滤波器加以近似。所以,我们也考虑y(n)≠x(n)但其中y(n)与x(n)在低频区域含有相同信息的状况。我们证明,整个系统可于|ωT|≤ω0T时达成Y(ejωT)=X(ejωT),Y(ejωT)与X(ejωT)分别为y(n)与x(n)的傅立叶转换,而ω0是预定限制频率,再次地,关于适当的理想合成滤波器,其在此情况下可以通过实际的数字滤波器加以近似。此方案对于(略微)过取样的模数转换器系统很有用,其可以忍受频带ω0T<|ωT|≤π中的频叠。理想合成滤波器是全通滤波器,其一般具有不同的增益常数。我们分析使用实际的滤波器来近似理想滤波器的效果。
此说明的其余部分的概要如下。首先,扼要重述均匀性取样、上取样与混合模拟/数字滤波器组,其中后者可以在分析非均匀性取样系统时方便地使用。下文中处理非均匀性取样与重建。其后,考虑时间交错的模数转换器及其概括化。接续部分分别涉及误差分析与使噪声作数字转换。最后,提供一方程式表,该方程式是在上文中所提到的。
均匀性取样、上取样与滤波器组
在图2中,均匀性取样与数字转换由均匀性取样器与数字转换器代表。忽略数字转换,通过对于全部的n,均匀地在时间间隔nT取样模拟输入信号xa(t)而得到输出序列x(n),请参看本说明结尾的方程式表中的方程式(1)。这里,T是取样周期,而fsample=1/T是取样频率,x(n)与xa(t)的傅立叶转换根据泊松求和公式而相关,请参看方程式(2)。
图3的上取样器用于以因子M增加取样频率。与较低速率有关的取样周期及取样频率,此处分别标示为T1与fsample,1,显然与方程式(3)中的T及fsample有关。输出序列y(n)得自于方程式(4),而如方程式(5),y(n)与x(m)的傅立叶转换彼此相关。
考虑图4的系统,我们将它称为混合模拟/数字滤波器组或滤波器组ADC。此系统使用一模拟分析滤波器组、均匀性取样器与数字转换器、及一数字合成滤波器组。取样与数字转换发生于解析滤波器的输出端,因为T1=MT,故取样频率为1/T1=fsample/M。在滤波器组ADC中,取样与数字转换二者于是以低取样速率fsample/M执行。
忽略图4所示系统中数字转换。输出序列y(n)的傅立叶转换可借助于以上关系容易地获得,参看方程式(6),其中Xk(ejMωT)得自于方程式(7)。方程式(6)可以重写为方程式(8),其中Vp(jω)得自于方程式(9)。
考虑图2与4所示的系统,X(ejωT)与Y(ejωT)分别得自于方程式(2)与方程式(8)。请记住,一取样信号的光谱恒为周期性,其周期为2π(2π周期)。于是,X(ejωT)显然为2π周期。只要全部Gk(ejωT)为2π周期,则这对于Y(ejωT)而言亦为真。因此,在-π≤ωT≤π间隔中考虑X(ejωT)与Y(ejωT)就足够了。现在,我们将处理二种不同类型的重建。
最佳重建:如果方程式(10)对某非零常数c与整数常数d成立,则图4的系统具有最佳重建(PR)。在时域中,我们有最佳重建的状况y(n)=cx(n-d)。亦即c=1时,y(n)只是x(n)的移位形式。自方程式(2)、(8)与(10),我们看到,如果对于-∞≤r≤∞,方程式(11)成立,则可以获得最佳重建。
区域性最佳重建:令x(n)与y(n)如方程式(12)所示而分解,对应的傅立叶转换得自于方程式(13)与(14),其中ω0T<π。如果方程式(15)——或等效地,方程式(16)——对某非零常数c与整数常数d成立,则图4的系统具有区域性最佳重建(RPR)。在时域中,我们有区域性最佳重建的状况,ylow(n)=cxlow(n-d)。亦即c=1时,ylow(n)只是xlow(n)的移位形式。然而,y(n)不是x(n)的移位形式,即,y(n)≠cx(n-d)。自方程式(2)、(8)与(16),我们看到,如果对于-∞≤r≤∞,满足方程式(17),则可以获得区域性最佳重建。区域性最佳重建系统重视过取样系统,其中xlow(n)承载基本信息,而xhigh(n)含有要由数字和/或模拟滤波器除去的所不希望的成分(例如噪声)。
限频状况:当Xa(jω)限频时,于-π≤ωT≤π的间隔,在方程式(2)与(8)的求和中,只有数目有限的项需要处理。我们考虑二种不同的状况。
状况A(最佳重建):令xa(t)依据方程式(18)而限频。在此情况下,达到以1/T的有效取样频率取样而无频叠的奈奎斯特(Nyquist)判据。因此若避免频叠进入频带-π≤ωT≤π,则可以保留xa(t)。
首先考虑图2中的x(n)。根据方程式(2),当Xa(jω)依据方程式(18)而限频时,显然,我们在区域-π≤ωT≤π中无频叠。其次,考虑图4中的y(n)。在区域-π≤ωT≤π中,Xa(jω)依据方程式(18)而限频,易证明,我们只需要考虑方程式(8)中的2K0+1项,p=-K0,-(K0-1),...,K0,其中K0得自于方程式(19)。
如果方程式(20)成立,其中K0得自于方程式(19),则现在获得最佳重建。故在此情况下,只要图4中的系统具有最佳重建,则xa(t)可自x(n)与y(n)而保留。
状况B(区域性最佳重建状况):令xa(t)依据方程式(21)而限频,且如方程式(22)所示而分解,其中对应的傅立叶转换得自于方程式(23)、(24)与(25)。
在此情况下,只要避免频叠进入频带-ω0T≤ωT≤ω0T,则不能保留xa(t),但可以保留xa,low(t)。
首先考虑图2中的x(n)。在区域-π≤ωT≤π中,Xa(jω)依据方程式(21)与(25)而限频,显然,我们只需要考虑方程式(2)中的3项,r=-1,0,1。此外,在区域-ω0T≤ωT≤ω0T中,ω0得自于方程式(25)。易于证明,我们仅需考虑一项:r=0。也就是说,进入此频带的频叠是自动避免的。其次,考虑图4中的y(n)。在区域-π≤ωT≤π中,Xa(jω)依据方程式(21)与(25)而限频,易证明,我们只需要考虑方程式(8)中的2K0+1项:p=-K0,-(K0-1),...,K0,而K0得自于方程式(26),其中
Figure C0181676400161
代表大于或等于x的最小整数。此外,在区域-ω0T≤ωT≤ω0T中,ω0得自于方程式(25)。易证明,我们只需要考虑方程式(8)中的2K0+1项:p=-K0,-(K0-1),...,K0,而K0得自于方程式(27)。
如果满足方程式(28),其中K0得自于方程式(27),且A(jω)为某任意函数,则现在可以获得区域性最佳重建。故在此情况下,只要图4中的系统具有区域性最佳重建,则可以自x(n)与y(n)保留xa.low(t)。
非均匀性取样与重建
令xk(m),k=0,1,...,N-1,是在时间间隔t=nMT+tk经由取样而获得的N个子序列,即如从方程式(29)所得出的。就M=N=2而言,xa(t)依据图1b进行取样。
如果依据方程式(30)来选择图4中的这些解析滤波器,则子序列xk(m)可以通过自该解析滤波器将输出信号取样而获得。在此情况下,解析滤波器组如图5所示。
结合方程式(9)与(30),我们得到方程式(31)。
以下显示如何在限频状况A与B(请参看前文)中选择合成滤波器,从而分别获得最佳重建与区域性最佳重建。
状况A(最佳重建状况):在此情况下,xa(t)依据方程式(18)而限频。令Gk(ejωT)为得自于方程式(32)的2π周期滤波器。由方程式(31)与(32),可以获得方程式(33)。就最佳重建而言,需要令得自于方程式(33)的Vp(jω)满足方程式(20)。即,如果满足方程式(34),即获得最佳重建。
状况B(区域性最佳重建状况):在此情况下,xa(t)依据方程式(21)而限频。令Gk(ejωT)得自于方程式(35)的2π周期滤波器,其中Ak(ejωT)为某任意复数函数。由方程式(31)与(35),我们获得方程式(36),其中A(jω)得自于方程式(37)。
就区域性最佳重建而言,需要使得由方程式(36)所给出的Vp(jω)满足方程式(28)。即,再次地,如果满足方程式(34),即获得区域性最佳重建。
其次,考虑如何计算ak。就最佳重建与区域性最佳重建(状况A与B)二者而言,必须满足方程式(34)。此方程式可以写为如同方程式(38)的矩阵形式,其中B是依据方程式(39)的(2K0+1)×N阶矩阵,uk得自于方程式(40)。此外,分别依据方程式(41)与(42),a为有N个元素的列向量而c是有2K0+1元素的列向量,其中T代表转置(无复数共轭部分)。ak为未知数,而ck得自于方程式(43)。
方程式(38)是有2K0+1个方程式的线性系统,其中有N个未知参数ak。因此,如果2K0+1≤N,则方程式(38)有解。我们将二种不同的状况加以区别。
状况1:2K0+1=N。在此情况下,未知数的数目与方程式的数目相等。在以下命题所述的条件下,可以唯一确定此状况中的ak
命题1:如果B与c分别得自于方程式(39)与(42),2K0+1=N,且tk≠tm+MTr,k≠m,r∈Z,则存在唯一的a,其满足方程式(38),且唯一的ak也满足方程式(34)。此外,a中的全部ak均是实数值常数。
证明:我们首先证明存在唯一的解。既然2K0+1=N,则B是N×N阶方阵。如果B非奇异,则a由方程式(44)唯一地确定,其中B-1是B的逆矩阵。因此在所述的条件下,表明B非奇异就足够了。为此目的,我们首先要看到,得自于方程式(39)的B可以写成方程式(45),其中A得自于方程式(46),而C则是依据方程式(47)的对角线矩阵。
矩阵A是范德蒙德(Vandermonde)矩阵。所以,A非奇异的必要与充分条件为,uk是不同的,即,uk≠um,k≠m,而由于方程式(40),这是与tk≠tm+MTr,k≠m,r∈Z相同的条件。此外,因为B的行列式为det B=det A det C,且|det C|=1,我们获得方程式(48)所给出的关系。即,若当且唯当A为非奇异,B为非奇异。这就证明在所述条件下,B为非奇异,且恒存在a的唯一解。
为证明a中的ak是实数值常数,我们进行如下。假设我们具有满足方程式(34)的唯一值ak。利用方程式(40),则方程式(34)同样可以写为方程式(49),其中x*代表x的复数共轭部分。由方程式(49),我们得到方程式(50)。这表明,值ak *也满足方程式(34)。然而,因为ak是唯一的,即得出它们必为实数值。
状况2:2K0+1<N。在此情况下,未知数的数目超过了方程式的数目。所以,我们可以在ak中加入L=N-2K0-1个额外的线性约束条件,而仍然满足方程式(34)。此处,我们将自己限制于这样一种状况:其中对于k=N-L+1,N-L+2,...,N,L个ak固定于某些常数。这种状况涵盖具有偶数频道的传统时间交错的模数转换器。因为L个ak可以自由选择,所以,在将要除去对应频道的状况下,我们当然可以将它们设定为零。如此一来,不需要考虑具有偶数频道的状况。然而,以下我们将可看到,可能值得考虑这些状况,为的是减少在整个系统的输出端的数字化噪声。
待解的线性方程式系统在此处可写为如同方程式(51)的矩阵形式,而分别依据方程式(52)、(53)与(54),
Figure C0181676400181
是N×N阶矩阵,a与
Figure C0181676400182
是具有N元素的列向量,其中B为得自于方程式(39)的(2K0+1)×(2K0+1)阶矩阵,au与afix分别含有a的(2K0+1)个未知数与L个固定常数,c为得自于方程式(43)的具有(2K0+1)个元素的列向量,S为得自于方程式(55)的L×N阶矩阵,其中Sz是得自于方程式(56)的L×(2K0+1)阶空矩阵,Sd是L×L对角线矩阵,其中对角线元素等于一,请参看方程式(57)。
如同状况1,状况2中的ak可以在以下命题所述条件下唯一地确定。命题2:如果
Figure C0181676400183
Figure C0181676400184
分别得自于方程式(52)与(54),方程式(53)中的afix含有L个实数值固定常数,2K0+1<N,且tk≠tm+MTr,k≠m,r∈Z,则存在唯一的a,其满足方程式(51),因而亦存在唯一的一组ak,其满足方程式(34)。此外,a中的全部ak均是实数值常数。
证明:本证明接续命题1的证明。为了证明存在性与唯一性,证实在所述条件下是非奇异就足够了,其原因为a是由方程式(58)唯一确定的。
为了证明
Figure C0181676400186
的非奇异性,我们观察到,它的行列式得自于方程式(59),其中
Figure C0181676400187
是自B将k=N-L+1,N-L+2,...,N共L行消去而得的(2K0+1)×(2K0+1)阶子矩阵,即,如方程式(60)所示。我们从命题1的证明知道, det B ~ ≠ 0 , 于是,在所述条件下, det B ^ ≠ 0 . 这样就证明了是非奇异的,且唯一的解恒存在。a中的ak是实数值的证明方式与命题1相同。
时间交错的模数转换器及其概括化
本节考虑传统时间交错的模数转换器及其概括化。首先考虑N=M而tk得自于方程式(61)与(62)的状况。
此外,令合成滤波器Gk(ejωT)是通过使方程式(32)中的ak=1,k=0,1,...,M-1,c=1且d=0而得出,即,如方程式(63)所示。由方程式(31)与(63),我们得到方程式(64)。
由此,获得最佳重建。在此情况下,我们具有一个时间交错的模数转换器。此处,输出序列y(n)是通过使xk(m)交错而得出。
然而,实际上Δtk不再恰等于零。若Δtk已知,如果N为奇数且2K0+1=N,则ak可以依据方程式(44)计算,而如果2K0+1<N,则依据方程式(58)计算。在此情况下,不能达到最佳重建,其原因为N=M且最佳重建需要使K0=M-1。于是,既不能满足2K0+1=N且不能满足2K0+1<N。另一方面,可以获得区域性最佳重建。就此状况而言,产生下列问题:已知N=M与K0,则我们可以允许且仍然可以获得区域性最佳重建的ω0T的最大值为何?易于证实,为了达成区域性最佳重建,我们必须满足方程式(65)。如果2K0+1=N,则我们得到方程式(66)。
其次,考虑N≠M而tk得自于方程式(67)与(68)的状况。此外,令合成滤波器Gk(ejωT)是通过使方程式(32)中的ak=M/N,k=0,1,...,N-1,c=1且d=0而得出,即,如方程式(69)所示。由方程式(31)与(69),我们得到方程式(70)。
由此,获得最佳重建。在此情况下,我们具有一系统,其可以视为时间交错的模数转换器的概括化。然而,在此情况下,我们不再能够通过使xk(m)交错而获得输出序列。
再次地,实际上,Δtk不再恰等于零。若Δtk已知,如果N为奇数且2K0+1=N,则ak可以依据方程式(44)计算,而如果2K0+1<N,则依据方程式(58)计算。与M频道的状况相反,此处,我们在N频道的状况下通过分别依据方程式(19)与(27)来选择K0,且当然通过选择N,从而使得2K0+1≤N,则可以达成最佳重建与区域性最佳重建二者。为了达成区域性最佳重建,就已知的M与K0而言,ω0T必须再次满足方程式(65)。如果2K0+1=N,则我们得到方程式(71)。因此,通过增加频道数目,我们在一更宽的频率区域获得区域性最佳重建。
误差与噪声解析
其次,提供误差分析。更精确地说,当B与a分别为B+ΔB与a+Δa所取代的时候,我们导出a与c中误差的界限。在有关数字化噪声方面,a中的误差是重点所在,在下文中将可明白。c中的误差则告诉我们任何实际的滤波器必须与理想合成滤波器有何等近似,以满足c的某些指定允许误差。
我们将要利用方程式(72)所定义的L范数以用于一具有元素xi的N×1(1×N)阶向量x,及方程式(73)所定义的L范数以用于一具有元素xik的N×N阶矩阵X。
a的误差:首先,考虑2K0+1=N的状况1。首先假设对于tk=dkT与ak,我们有Ba=c。其次,假设tk=dkT与ak分别为tk=dkT+Δtk与ak+Δak所取代,而c保持固定。由此得到方程式(74)。矩阵ΔB为依据方程式(75)的N×N阶矩阵,其中Δbpk与Δtpk分别得自于方程式(76)与(77)。
现在,如果满足方程式(78),则可以证实方程式(79)成立。由方程式(75)~(77),我们可得方程式(80)。
我们有B=AC,并从而得到B-1=C-1A-1。此外,因为此处的A是DFT矩阵,其逆矩阵A-1是IDFT矩阵;因此,‖A‖=1。我们也得到‖C-1=1,其原因为显然C-1是具有对角线元素uk k 0的对角线矩阵,其中uk得自于方程式(40)。于是,我们得到方程式(81),其与方程式(80)一起,导出方程式(82)。利用方程式(79)~(82),且假设‖ΔB‖‖B-1<<1,我们最后得到方程式(83)。
其次,考虑2k0+1<N的状况2。此状况比状况1略为困难,其原因为我们通常不能以DFT矩阵与对角线矩阵的积表示
Figure C0181676400201
。然而,如果我们将自己限制在时间交错的模数转换器及其概括化,则显然我们可以将方程式(51)重写为方程式(84),其中B′是依据方程式(85)的N×N阶矩阵,而uk得自于方程式(40),且c′为依据方程式(86)的具有N个元素ck的列向量。
显然,我们可以将B′表示为一DFT矩阵与一对角线矩阵的积。所以我们将以和状况1相同的结果作为结束,即,以方程式(83)为界。
c的误差:假设tk=dkT与ak时,我们得到Ba=c。现在假设tk=dkT与ak分别为tk=dkT+Δtpk与ak+Δak所取代。这样得出方程式(87),我们由该方程式可得到方程式(88)。依次地,由方程式(88),我们得到方程式(89)。利用方程式(39)与方程式(75)~(77),我们最后得到方程式(90),其在设计合成滤波器Gk(z)时很有用。
由上述,请记住理想滤波器应具有在重点频率范围上的频率响应ake-jωt k(如果在方程式(32)与(35)中c=1且d=0的话)。实际上,Gk(z)当然可以只近似于理想响应。我们可以将Gk(z)的频率响应表示为方程式(91),其中Δak(ωt)与Δtpk(ωt)分别是理想大小与相位响应的偏移量。已知c的允许误差与方程式(90)及(91),于是可以容易地设计Gk(z),以使需求得到满足。
为了分析在图4的系统输出端的噪声变化,方便的是依据图6,用其所谓多相实现方式(polyphase realization)来表示合成滤波器组。输出序列y(n)是通过将yi(m),i=0,1,...,M-1交错而获得。输出y(n)的转换函数得自于方程式(92),其中Y(z)得自于方程式(93),而X(z)、Y(z)、与G(p)(z)分别在方程式(94)、(95)与(96)中定义。Gik(z)是依据方程式(97)的Gk(z)的多相成分。
通常,在噪声分析中,数字转换误差经模型化而成为静止的白噪声(white noise)。令xk(m),k=0,1...,N-1为具有零平均值与变化量σxk 2的不相关的白噪声源。既然G(p)(z)描述一线性与非时变系统,则输出yi(m),i=0,1,...,M-1亦为具有零平均值的静止的白噪声。然而,yi(m)的变化量,此处标为σyi 2(n),大体上是不同的,即使当σxk 2相等时亦然。输出yi(m)也可以是相关的。所以,输出噪声y(n)大体上将不静止。因此,它的变化量,此处标为σy 2(n),是随时间变化的。因为显然方程式(98)成立,故它又是周期性的,周期为N。
我们在方程式(99)中定义输出端的平均数字化噪声。已知合成滤波器Gk(z)及它的多相成分Gik(z),则可在方程式(100)中计算(σy 2)ax
现在,令合成滤波器得自于方程式(101),且依据方程式(102),全部输入变化量σxk 2相等。结合方程式(100)~(102),我们就得到方程式(103)。
现在发生的一问题是如何选择ak,以使得自于方程式(103)的(σy 2)av在同时达成最佳重建或区域性最佳重建的约束下减至最小。让我们将问题视为由方程式(104)所定义。方程式(104)中的约束是为了获得最佳重建或区域性最佳重建所必须满足的约束之一。因为ak之和为M,故方程式(104)中待减至最小的目标函数可以重写为方程式(105)。因此,获得方程式(104)中对于ak=M/N,k=0,1,...,N-1的解,其中(σy 2)ax的最小值如方程式(106)所示。
这表明,选择ak=M/N以用于时间交错的模数转换器及它们的概括化使得输出端的平均数字化噪声减至最小。
实际上Δtk不再恰为零,其意味着ak为ak+Δak所取代。如果Δak小(且ak>0),则平均数字化噪声在此情况下得自于方程式(107)。以ak=M/N,我们得到方程式(108)。数量得自于方程式(83)。
本发明已经考虑使用数字滤波器组的非均匀性取样的限频信号的重建问题。整个系统可视为传统的时间交错的模数转换器的概括化,而前者简化为其特殊状况。通过将时间交错的模数转换器概括化,可以消除实际上由于时间偏移误差所引入的误差。我们考虑最佳重建(PR)与区域性最佳重建(RPR)系统二者,且表明了如何通过适当选择(理想)数字滤波器而获得上述系统。
非均匀性取样的限频信号的重建方法可以任何适当的数字信号处理装置上实行,例如专用硬件或计算机。在后者的状况,所述方法的执行要通过一包括软件代码部分的计算机程序产品,其加载到一适当装置的内部存储器内。
显然,本发明可通过很多方式加以变动。这种变动不应视为偏离本发明的领域。对本领域技术人员显而易见的所有这种变动均应认为包含于所附权利要求的范围内。
方程式表展示在后续各页。
方程式表
x(n)=x3(t)|t=nT  -∞≤n≤∞(1)
X ( e jωT ) = 1 T Σ r = - ∞ ∞ X a ( jω - j 2 πr T ) - - - ( 2 )
T1=MT, f sample , 1 = f sample M - - - ( 3 )
Figure C0181676400233
Y ( e jωT ) = X ( e jω T 1 ) = X ( e jMωT ) - - - ( 5 )
Y ( e jωT ) = Σ k = 0 N - 1 G k ( e jωT ) X k ( e jMωT ) - - - ( 6 )
X k ( e jMωT ) = X k ( e jω T 1 )
= 1 T 1 Σ p = - ∞ ∞ H k ( jω - j 2 πp T 1 ) X a ( jω - j 2 πp T 1 )
= 1 MT Σ p = - ∞ ∞ H k ( jω - j 2 πp MT ) X a ( jω - j 2 πp MT ) - - - ( 7 )
Y ( e jωT ) = 1 T Σ p = - ∞ ∞ V p ( jω ) X a ( jω - j 2 πp MT ) - - - ( 8 )
V p ( jω ) = 1 M Σ k = 0 N - 1 G k ( e jωT ) H k ( jω - j 2 πp MT ) - - - ( 9 )
Y(ejωT)=ce-jdωTX(ejωT),|ωT|≤π      (10)
V p ( jω ) = c e - jdωT , p = rM , | ω | ≤ π / T 0 , p ≠ rM , | ω | ≤ π / T - - - ( 11 )
x(n)=xlow(n)+xhigh(n)
y(n)=ylow(n)+yhigh(n)
X(ejωT)=Xlow(ejωT)+Xhigh(ejωT)(13)
Y(ejωT)=Ylow(ejωT)+Yhigh(ejωT)
Xlow(ejωT)=0,ω0T<|ωT|≤π
Xhigh(ejωT)=0,|ωT|≤ω0T
Ylow(ejωT)=0,ω0T<|ωT|≤π(14)
Yhigh(ejωT)=0,|ωT|≤ω0T
Y(ejωT)=ce-jdωTX(ejωT),|ωT|≤ω0T    (15)
Ylow(ejωT)=ce-jdωTXlow(ejωT),|ωT|≤π(16)
V p ( jω ) = c e - jdωT , p = rM , | ω | ≤ ω 0 0 , p ≠ rM , | ω | ≤ ω 0 - - - ( 17 )
X3(jω)=0,|ω|≥π/T                     (18)
K0=M-1                            (19)
V p ( jω ) = c e - jdωT , p = 0 , | ω | ≤ π 0 , | p | = 1 , 2 , . . . , K 0 , 0 ≤ ω ≤ π - - - ( 20 )
X2(jω)=0,|ω|≥ω1              (21)
xa(t)=xa,low(t)+xa,high(t)      (22)
Xa(jω)=X3,low(jω)+Xa,high(jω)(23)
Xlow(jω)=0,|ω|>ω0            (24)
Xhigh(jω)=0,|ω|≤ω0,|ω|≥ω1
0<ω0<ω1,ω01≤2π/TT       (25)
Figure C0181676400252
Figure C0181676400253
V p ( j&omega; ) = c e - jd&omega;T , p = 0 , | &omega; | &le; &omega; 0 A ( j&omega; ) , p = 0 , &omega; 0 < | &omega; | &le; &pi; / T 0 , | p | = 1,2 , . . . , K 0 , | &omega; | &le; &omega; 0 - - - ( 28 )
xk(m)=x(nMT+tk),k=0,1,...,N-1(29)
H k ( s ) = e s t k , k = 0,1 , . . . , N - 1 - - - ( 30 )
V p ( j&omega; ) = 1 M &Sigma; k = 0 N - 1 G k ( e j&omega;T ) e j ( &omega; - 2 &pi;p MT ) t k - - - ( 31 )
G k ( e j&omega;T ) = a k c e - j&omega; ( t k + dT ) , | &omega;T | < &pi; - - - ( 32 )
V p ( j&omega; ) = 1 M c e - jd&omega;T &Sigma; k = 0 N - 1 a k e - j 2 &pi;p MT t k - - - ( 33 )
&Sigma; k = 0 N - 1 a k e - j 2 &pi;p MT t k = M , p = 0 0 , | p | = 1,2 , . . . , K 0 - - - ( 34 )
G k ( e j&omega;T ) = a k c e - j&omega; ( t k + dT ) , | &omega;T | &le; &omega; 0 T a k A k ( e j&omega;T ) , &omega; 0 T < | &omega;T | &le; &pi; - - - ( 35 )
V p ( j&omega; ) = 1 M c e - jd&omega;T &Sigma; k = 0 N - 1 a k e - j 2 &pi;p MT t k , | &omega;T | &le; &omega; 0 T A ( j&omega; ) , &omega; 0 T < | &omega;T | &le; &pi; - - - ( 36 )
A ( j&omega; ) = 1 M &Sigma; k = 0 N - 1 a k A k ( e j&omega;T ) e j ( &omega; - 2 &pi;p MT ) t k - - - ( 37 )
Ba=c    (38)
B = u 0 - K 0 u 1 - K 0 . . . u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) . . . u N - 1 - ( K 0 - 1 ) . . . u 0 K 0 u 1 K 0 . . . u N - 1 K 0 - - - ( 39 )
u k = e - j 2 &pi; MT t k - - - ( 40 )
a=[a0a1...aN-1]T    (41)
c = [ c 0 c 1 . . . c 2 K 0 ] T - - - ( 42 )
c k = M , k = K 0 0 , k = 0,1 , . . . , 2 K 0 , k &NotEqual; K 0 - - - ( 43 )
a=B-1c              (44)
B=AC                (45)
A = 1 1 . . . 1 u 0 u 1 . . . u N - 1 . . . u 0 2 K 0 u 1 2 K 0 . . . u N - 1 2 K 0 - - - ( 46 )
c = diag [ u 0 - K 0 u 1 - K 0 . . . u N - 1 - K 0 ] - - - ( 47 )
det A &NotEqual; 0 &DoubleLeftRightArrow; det B &NotEqual; 0 (48)
det A = 0 &DoubleLeftRightArrow; det B = 0
&Sigma; k = 0 N - 1 a k = M , p = 0 (49)
&Sigma; k = 0 N - 1 a k [ u k p ] * = &Sigma; k = 0 N - 1 a k u k p = 0 , p = 1,2 , . . . , K 0
&Sigma; k = 0 N - 1 a k * = M , p = 0 (50)
&Sigma; k = 0 N - 1 a k * u k p = &Sigma; k = 0 N - 1 a k * [ u k p ] * = 0 , p = 1,2 , . . . , K 0
B ^ a = c ^ - - - ( 51 )
B ^ = B S - - - ( 52 )
a=[auafik]T    (53)
c ^ = c a fix T - - - ( 54 )
S=[SzSd]       (55)
S z = 0 0 . . . 0 0 0 . . . 0 . . . 0 0 . . . 0 - - - ( 56 )
Sd=diag[1 1...1]                  (57)
a = B ^ - 1 c ^ - - - ( 58 )
det B ~ = det B ~ &Pi; l = 0 L - 1 S d , H = det B ~ - - - ( 59 )
B ~ = u 0 - K 0 u 1 - K 0 . . . u 2 K 0 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) . . . u 2 K 0 - ( K 0 - 1 ) . . . u 0 K 0 u 1 K 0 . . . u 2 K 0 K 0 - - - ( 60 )
tk=dkT+Δtk,k=0,1,...,M-1    (61)
dk=k,k=0.1,...,M-1
                                   (62)
Δtk=0,k=0,1,...,M-1
Gk(ejωT)=e-jkωT,|ωT|<π      (63)
V p ( j&omega; ) = 1 M &Sigma; k = 0 M - 1 e - &Integral; 2 &pi;pk M = 1 , p = 0 0 , p &NotEqual; 0 - - - ( 64 )
&omega; 0 T &le; 2 &pi; ( K 0 + 1 ) M - &omega; 1 T - - - ( 65 )
&omega; 0 T &le; &pi; ( M + 1 ) M - &omega; 1 T = &pi; M + &pi; - &omega; 1 T - - - ( 66 )
tk=dkT+Δtk,k=0,1,...,N-1    (67)
d k = kM N , k=0,1,...,N-1
                              (68)
Δtk=0,k=0,1,...N-1
G k ( e j&omega;T ) = M N e - jMk&omega;T N , | &omega;T | < &pi; - - - ( 69 )
V p ( j&omega; ) = 1 N &Sigma; k = 0 N - 1 e - j 2 &pi;pk N = 1 , p = 0 0 , p &NotEqual; 0 - - - ( 70 )
&omega; 0 T &le; &pi; ( N + 1 ) M - &omega; 1 T - - - ( 71 )
‖x‖=max|xi|,0≤i≤N-1   (72)
| | X | | &infin; = max &Sigma; k = 0 N - 1 | x ik | , 0 &le; i &le; N - 1 - - - ( 73 )
(B+ΔB)(a+Δa)=ε.             (74)
&Delta;B = &Delta; b - K 0 , 0 &Delta; b - K 0 , 1 . . . &Delta; b - K 0 , N - 1 &Delta; b - ( K 0 - 1 ) , 0 &Delta; b - ( K 0 - 1 ) , 1 . . . &Delta; b - ( K 0 - 1 ) , N - 1 . . . &Delta; b K 0 , 0 &Delta; b K 0 , 1 . . . &Delta; b K 0 , N - 1 - - - ( 75 )
&Delta; b pk = e j 2 &pi;pk M ( e j&Delta; t pk - 1 ) - - - ( 76 )
&Delta; t pk = 2 &pi;p MT &Delta; t k - - - ( 77 )
‖ΔB‖·‖B-1<1            (78)
| | &Delta;a | | &infin; | | a | | &infin; &le; | | &Delta;B | | &infin; &CenterDot; | | B - 1 | | &infin; 1 - | | &Delta;B | | &infin; &CenterDot; | | B - 1 | | &infin; - - - ( 79 )
| | &Delta;B | | &infin; = max &Sigma; k = 0 N - 1 | &Delta; b pk | &ap; max &Sigma; k = 0 N - 1 | &Delta; t pk | - - - ( 80 )
&le; N ( N - 1 ) &pi; M max { | &Delta; t k | &tau; } .
‖B-1≤‖C-1·‖A-1=1.(81)
| | &Delta;B | | &infin; &CenterDot; | | B - 1 | | &infin; &le; N ( N - 1 ) &pi; M &CenterDot; max { | &Delta; t k | T } . - - - ( 82 )
| | &Delta;a | | &infin; &le; | | a | | &infin; N ( N - 1 ) &pi; M &CenterDot; max { | &Delta; t k | T } . - - - ( 83 )
Bta=c1                            (84)
B &prime; = u 0 - K 0 u 1 - K 0 . . . u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) . . . u N - 1 - ( K 0 - 1 ) . . . u 0 N - K 0 - 1 u 1 N - K 0 - 1 . . . u N - 1 N - K 0 - 1 - - - ( 85 )
c k = M , k - K 0 0 , k = 0,1 , . . . N - 1 , k &NotEqual; K 0 - - - ( 86 )
(B+ΔB)(a+Δa)=c+Δc                                          (87)
Δc=BΔa+ΔBa+ΔBΔa                                          (88)
‖Δc‖≤‖B‖‖Δa‖+‖ΔB‖‖a‖+‖ΔB‖‖Δa‖(89)
‖Δc‖≤Nmax{|Δak|}+Nmax{|Δtpk|}max{|ak|}
+Nmax{|Δtpk|}max{|Δak|}                                      (90)
=N(max{|Δak|}+max{|Δtpk|}max{|ak})
G k ( e j&omega;T ) = [ a k + &Delta; a k ( &omega;T ) ] e - j [ &omega; t k + &Delta; t pk ( &omega;T ) ] - - - ( 91 )
Y ( z ) = &Sigma; i = 0 M - 1 z - i Y i ( z M ) - - - ( 92 )
Y(=)=G(P)(z)X(z)                                             (93)
X(z)=[X0(z)X1(z)...XN-1(z)]T                                  (94)
Y()=[Y0(z)Y1(z)...YM-1(z)]T                                   (95)
G ( p ) ( z ) = G 00 ( z ) G 01 ( z ) . . . G 0 , N - 1 ( ) G 10 ( z ) G 11 ( z ) . . . G 1 , N - 1 ( z ) . . . G M - 1 , 0 ( z ) G M - 1 , 1 ( z ) . . . G M - 1 , N - 1 ( z ) - - - ( 96 )
G k ( z ) = &Sigma; i = 0 M - 1 z - i G ik ( z M ) - - - ( 97 )
&sigma; y 2 ( nM + i ) = &sigma; y i 2 - - - ( 98 )
( &sigma; y 2 ) av = 1 M &Sigma; i = 0 M - 1 &sigma; y i 2 - - - ( 99 )
( &sigma; y 2 ) av = 1 M &Sigma; i = 0 M - 1 &sigma; y i 2 = 1 M &Sigma; i = 0 M - 1 &Sigma; k = 0 N - 1 &sigma; x k 2 &Sigma; n = - &infin; &infin; | g ik ( n ) | 2
= 1 M &Sigma; k = 0 N - 1 &sigma; x k 2 &Sigma; i = 0 M - 1 &Sigma; n = - &infin; &infin; | g ik ( n ) | 2
= 1 M &Sigma; k = 0 N - 1 &sigma; x k 2 &Sigma; n = - &infin; &infin; | g k ( n ) | 2
= 1 M &Sigma; k = 0 N - 1 &sigma; x k 2 1 2 &pi; &Integral; - &pi; &pi; | G k ( e j&omega;T ) | 2 d&omega;T - - - ( 100 )
G k ( e j&omega;T ) = a k e - jMk&omega;T / N , | &omega;T | < &omega; c T 0 , &omega; c T &le; | &omega;T | &le; &pi; - - - ( 101 )
&sigma; x k 2 = &sigma; x 2 . - - - ( 102 )
( &sigma; y 2 ) av = &omega; c T M&pi; &sigma; x 2 &Sigma; k = 0 N - 1 a k 2 - - - ( 103 )
min imize &Sigma; k = 0 N - 1 a k 2 subject to &Sigma; k = 0 N - 1 a k = M - - - ( 104 )
&Sigma; k = 0 N - 1 a k 2 = ( &Sigma; k = 0 N - 1 a k ) 2 + &Sigma; k = 0 N - 2 &Sigma; q = k - 1 N - 1 ( a k - a q ) 2
= M 2 + &Sigma; k = 0 N - 2 &Sigma; q = k + 1 N - 1 ( a k - a q ) 2 - - - ( 105 )
( &sigma; y 2 ) av , min = M &omega; c T N&pi; &sigma; x 2 - - - ( 106 )
( &sigma; y 2 ) av = &omega; c T M&pi; &sigma; x 2 &Sigma; k = 0 N - 1 ( a k + &Delta; a k ) 2
= &omega; c T M&pi; &sigma; x 2 &Sigma; k = 0 N - 1 ( a k 2 + 2 a k &Delta; a k )
&le; &omega; c T M&pi; &sigma; x 2 &Sigma; k = 0 N - 1 ( a k 2 + 2 a k | &Delta; a k | max ) - - - ( 107 )
( &sigma; y 2 ) av &le; &omega; c T M&pi; &sigma; x 2 ( M 2 N + 2 M | &Delta; a k | max )
= M &omega; c T N&pi; &sigma; x 2 ( 1 + 2 N | &Delta; a k | max M )
= ( &sigma; y 2 ) av , min ( 1 + 2 N | &Delta; a k | max M ) - - - ( 108 )

Claims (19)

1.一种非均匀性取样的限频模拟信号xa(t)的重建方法,该非均匀性取样信号包括N个子序列xk(m),k=0,1,...,N-1,N≥2,其是以1/(MT)的取样速率依据xk(m)=xa(nMT+tk)取样而获得,其中M为正整数,T为采样周期,而tk=kMT/N+Δtk,Δtk不为零,该方法特征在于以下步骤:
-自该N个子序列xk(m)形成新序列y(n),使得y(n)至少含有与x(n)=xa(nT)相同的信息,即,在一低于ω0的频率区域中,ω0是预定限制频率,以1/T的取样速率使xa(t)通过下列步骤完成取样:
(i)以因子M上取样该N个子序列xk(m)的每一项,k=0,1,...,N-1;
(ii)以一对应的数字滤波器过滤该上取样的N个子序列xk(m)的每一项,k=0,1,...,N-1;及
(iii)添加该N个数字滤波的子序列,以形成y(n)。
2.根据权利要求1所述的方法,其中该对应的数字滤波器是分段延迟滤波器,且具有在频带|ωT|≤ω0T中的频率响应 G k ( e j&omega;T ) = a k ce - j&omega; ( t k + dT ) ,k=0,1,...,N-1,ak为常数,d为整数,c为 c = [ c 0 , c 1 &CenterDot; &CenterDot; &CenterDot; c 2 K 0 ] T .
3.根据权利要求2所述的方法,其中该对应的分段延迟滤波器具有在频带ω0T<|ωT|≤π中的频率响应Gk=akAk(ejωT),k=0,1,...,N-1,其中Ak(ejωT)是任意复数函数。
4.根据权利要求2至3中任一项所述的方法,其中选择ak使其满足
&Sigma; k = 0 N - 1 a k e - j 2 &pi;p MT t k = M , p = 0 0 , | p | = 1,2 , . . . , K 0
其中K0得自于
K 0 = [ M ( &omega; 0 T + &omega; 1 T ) 2 &pi; ] - 1
其中
Figure C018167640002C5
应读作是大于或等于
Figure C018167640002C6
的最小的
整数,且[-ω1,ω1]是该限频模拟信号xa(t)所处的频带。
5.根据权利要求4所述的方法,其中ak计算如下,
a=B-1c,
a是ak的向量形式,且得自于
a=[a0 a1...an-1]T
B-1是B的逆矩阵,且B得自于
B = u 0 - K 0 u 1 - K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u N - 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u 0 K 0 u 1 K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 K 0
其中
u k = e - j 2 &pi; MT t k
而c为
c = [ c 0 , c 1 &CenterDot; &CenterDot; &CenterDot; c 2 K 0 ] T ,其中当2K0+1=N时,ck按照如下计算:
c k = M , k = K 0 0 , k = 0,1 , . . , 2 K 0 , k &NotEqual; K 0 .
6.根据权利要求4所述的方法,其中ak计算如下,
a = B ^ - 1 c ^ ,
a定义为
a=[au afix]T
其中au与afix含有(2K0+1)个未知数ak以及L=N-2K0-1个固定常数ak,而
Figure C018167640003C6
Figure C018167640003C7
的逆矩阵,得自于下式,
B ^ = B S
其中B由下式得出
B = u 0 - K 0 u 1 - K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u N - 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u 0 K 0 u 1 K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 K 0
其中
u k = e - j 2 &pi; MT t k ,且
S得自于
S=[SZ Sd],其中
S z = 0 0 &CenterDot; &CenterDot; &CenterDot; 0 0 0 &CenterDot; &CenterDot; &CenterDot; 0 &CenterDot; &CenterDot; &CenterDot; 0 0 &CenterDot; &CenterDot; &CenterDot; 0 ,
Sd=diag[11...1],且
Figure C018167640004C3
c ^ = c a fix T ,其中
c得自于
c = [ c 0 , c 1 &CenterDot; &CenterDot; &CenterDot; c 2 K 0 ] T ,且其中当2K0+1<N时,ck按照如下计算:
c k = M , k = K 0 0 , k = 0,1 , . . , 2 K 0 , k &NotEqual; K 0 .
7.根据权利要求6所述的方法,其中L=N-2K0-1个ak由ak=M/N,k=N-L+1,N-L+2,...,N来计算。
8.根据权利要求6所述的方法,其中L=N-2K0-1个ak由ak=0,k=N-L+1,N-L+2,...,N来计算。
9.根据权利要求1至2中任一项所述的方法,其中N=M。
10.根据权利要求1至2中任一项所述的方法,其中N≠M。
11.根据权利要求9所述的方法,其中ω0是依据下式选择:
&omega; 0 T &le; 2 &pi; ( K 0 + 1 ) M - &omega; 1 T
K0得自于
K 0 = [ M ( &omega; 0 T + &omega; 1 T ) 2 &pi; ] - 1 ,
其中应读作是大于或等于
Figure C018167640005C2
的最小的
整数,且[-ω1,ω1]是该限频模拟信号xa(t)所处的频带。
12.根据权利要求10所述的方法,其中ω0是依据下式选择:
&omega; 0 T &le; 2 &pi; ( K 0 + 1 ) M - &omega; 1 T
K0得自于
K 0 = [ M ( &omega; 0 T + &omega; 1 T ) 2 &pi; ] - 1 ,
其中
Figure C018167640005C5
应读作是大于或等于
Figure C018167640005C6
的最小的
整数,且[-ω1,ω1]是该限频模拟信号xa(t)所处的频带。
13.根据权利要求1所述的方法,其中该对应的数字滤波器是分段延迟滤波器,且具有在频带|ωT|≤π中的频率响应 G k ( e j&omega;T ) = a k ce - j&omega; ( t k - dT ) , k=0,1,...,N-1,ak为常数,d为整数,c为 c = [ c 0 , c 1 &CenterDot; &CenterDot; &CenterDot; c 2 K 0 ] T ,且因此所形成
的该新序列y(n)恰等于x(n)。
14.根据权利要求13所述的方法,其中选择ak,使其满足
&Sigma; k = 0 N - 1 a k e - j 2 &pi;p MT t k = M , p = 0 0 , | p | = 1,2 , . . . , K 0
其中K0得自于
K0=M-1。
15.根据权利要求14所述的方法,其中ak计算如下,
a=B-1c,
a是ak的向量形式,且得自于
a=[a0 a1...aN-1]T
B-1是B的逆矩阵,且B得自于
B = u 0 - K 0 u 1 - K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u N - 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u 0 K 0 u 1 K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 K 0
其中
u k = e - j 2 &pi; MT t k
而c为
c = [ c 0 , c 1 &CenterDot; &CenterDot; &CenterDot; c 2 K 0 ] T ,其中当2K0+1=N时,ck按照如下计算:
c k = M , k = K 0 0 , k = 0,1 , . . , 2 K 0 , k &NotEqual; K 0 .
16.根据权利要求14所述的方法,其中ak计算如下,
a = B ^ - 1 c ^ ,
a定义为
a=[au afix]T
其中au与afix含有(2K0+1)个未知数ak以及L=N-2K0-1个固定常数ak,而
Figure C018167640006C6
Figure C018167640006C7
的逆矩阵,
Figure C018167640006C8
得自于下式,
B ^ = B S
其中B由下式得出
B = u 0 - K 0 u 1 - K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 - K 0 u 0 - ( K 0 - 1 ) u 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u N - 1 - ( K 0 - 1 ) &CenterDot; &CenterDot; &CenterDot; u 0 K 0 u 1 K 0 &CenterDot; &CenterDot; &CenterDot; u N - 1 K 0
其中
u k = e - j 2 &pi; MT t k ,且
S得自于
S=[SZ Sd],其中
S Z = 0 0 &CenterDot; &CenterDot; &CenterDot; 0 0 0 &CenterDot; &CenterDot; &CenterDot; 0 &CenterDot; &CenterDot; &CenterDot; 0 0 &CenterDot; &CenterDot; &CenterDot; 0 ,且
Sd=diag[11...1],且
Figure C018167640007C2
c ^ = c a fix T ,其中
c得自于
c = [ c 0 , c 1 &CenterDot; &CenterDot; &CenterDot; c 2 K 0 ] T ,且其中当2K0+1<N时,ck按照如下计算:
c k = M , k = K 0 0 , k = 0,1 , . . . , 2 K 0 k &NotEqual; K 0 .
17.根据权利要求16所述的方法,其中L=N-2K0-1个ak由ak=M/N,k=N-L+1,N-L+2,...,N来计算。
18.根据权利要求13至17中任一项所述的方法,其中N>M。
19.根据权利要求1、2或13中任一项所述的方法,其中该N个子序列xk(m)在上取样以前经过数字转换。
CNB018167640A 2000-10-02 2001-09-28 非均匀性取样的限频信号的重建方法 Expired - Fee Related CN100379154C (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
SE00035493 2000-10-02
SE0003549A SE518092C2 (sv) 2000-10-02 2000-10-02 Förfarande resp. digital signalbehandlingsanordning för rekonstruktion av olikformigt samplade signaler.

Publications (2)

Publication Number Publication Date
CN1468433A CN1468433A (zh) 2004-01-14
CN100379154C true CN100379154C (zh) 2008-04-02

Family

ID=20281267

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB018167640A Expired - Fee Related CN100379154C (zh) 2000-10-02 2001-09-28 非均匀性取样的限频信号的重建方法

Country Status (9)

Country Link
US (1) US6476754B2 (zh)
EP (1) EP1330824B1 (zh)
CN (1) CN100379154C (zh)
AT (1) ATE386323T1 (zh)
AU (1) AU2001290483A1 (zh)
DE (1) DE60132797T2 (zh)
SE (1) SE518092C2 (zh)
TW (1) TW498623B (zh)
WO (1) WO2002029806A1 (zh)

Families Citing this family (20)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004039021A1 (en) * 2002-10-25 2004-05-06 The Trustees Of Columbia University In The City Of New York Time encoding and decoding of a signal
US6980147B2 (en) * 2003-04-07 2005-12-27 Photon Products, Inc. Channelized analog-to-digital converter
WO2004112298A2 (en) * 2003-05-27 2004-12-23 The Trustees Of Columbia University In The City Of New York Multichannel time encoding and decoding of a signal
US7233270B2 (en) * 2005-01-28 2007-06-19 Realtek Semiconductor Corp. Receiver capable of correcting mismatch of time-interleaved parallel ADC and method thereof
US7741982B2 (en) * 2005-02-04 2010-06-22 Signal Processing Devices Sweden Ab Estimation of timing errors in a time-interleaved analog to digital converter system
US7138933B2 (en) * 2005-04-26 2006-11-21 Analog Devices, Inc. Time-interleaved signal converter systems with reduced timing skews
DE102005060874A1 (de) * 2005-12-20 2007-06-21 Robert Bosch Gmbh Verfahren und Vorrichtung zur echtzeitfähigen Rekonstruktion eines nicht-äquidistant abgetasteten, analogen Signals
DE102006023697B4 (de) * 2006-05-19 2008-02-07 Infineon Technologies Ag Verfahren zum Decodieren, Decodierer, Codierer-Decodierer-System und Wandler
WO2008151137A2 (en) * 2007-06-01 2008-12-11 The Trustees Of Columbia University In The City Of New York Real-time time encoding and decoding machines
WO2009006405A1 (en) * 2007-06-28 2009-01-08 The Trustees Of Columbia University In The City Of New York Multi-input multi-output time encoding and decoding machines
US20100135446A1 (en) * 2008-12-03 2010-06-03 Electronics And Telecommunications Research Institute Digital-intensive rf receiver
US8825415B2 (en) 2008-12-16 2014-09-02 Signal Processing Devices Sweden Ab Methods and apparatuses for estimation and compensation on nonlinearity errors
CN102414988B (zh) * 2009-04-29 2014-03-19 香港大学 用于从不均匀样本重构均匀样本的方法或结构
KR100964378B1 (ko) * 2009-06-23 2010-06-25 한국전자통신연구원 디지털 수신기
US8532238B2 (en) * 2010-02-08 2013-09-10 Electronics And Telecommunications Research Institute Subsampling based receiver using frequency selective noise canceller
JP2011166773A (ja) * 2010-02-11 2011-08-25 Korea Electronics Telecommun サブサンプリング技法を利用する受信器のデジタル処理構造
WO2012109407A1 (en) 2011-02-09 2012-08-16 The Trustees Of Columbia University In The City Of New York Encoding and decoding machine with recurrent neural networks
US9209844B2 (en) 2012-01-11 2015-12-08 Electronics And Telecommunications Research Institute Subsampling receiver using interstage off-chip RF band pass filter
CN105659501A (zh) * 2013-10-18 2016-06-08 华为技术有限公司 一种模拟信息转换设备和方法
WO2017053635A1 (en) * 2015-09-25 2017-03-30 Google Inc. Increased sampling in non-uniform sampling analog-to-digital converters

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5610825A (en) * 1994-11-08 1997-03-11 Johnson; William J. Method and apparatus for the display of digitized analog signal loss
CN1247414A (zh) * 1998-08-25 2000-03-15 索尼公司 数字/模拟转换装置和重放装置
JP2000173186A (ja) * 1998-12-01 2000-06-23 Stmicroelectronics Srl 雑音指数を減少させた読取チャンネルにおける時間インタ―リ―ブ型デジタル信号処理

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4792787A (en) * 1987-02-04 1988-12-20 The United States Of America As Represented By The Secretary Of The Navy Wide dynamic range analog-to-digital converter using linear prediction
CA2282776C (en) * 1999-09-17 2004-10-05 Timothy John Borer Method and apparatus for high speed analog to digital conversion

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5610825A (en) * 1994-11-08 1997-03-11 Johnson; William J. Method and apparatus for the display of digitized analog signal loss
CN1247414A (zh) * 1998-08-25 2000-03-15 索尼公司 数字/模拟转换装置和重放装置
JP2000173186A (ja) * 1998-12-01 2000-06-23 Stmicroelectronics Srl 雑音指数を減少させた読取チャンネルにおける時間インタ―リ―ブ型デジタル信号処理

Also Published As

Publication number Publication date
ATE386323T1 (de) 2008-03-15
US6476754B2 (en) 2002-11-05
DE60132797D1 (de) 2008-03-27
WO2002029806A1 (en) 2002-04-11
US20020075177A1 (en) 2002-06-20
SE0003549D0 (sv) 2000-10-02
AU2001290483A1 (en) 2002-04-15
CN1468433A (zh) 2004-01-14
DE60132797T2 (de) 2009-02-12
TW498623B (en) 2002-08-11
EP1330824B1 (en) 2008-02-13
EP1330824A1 (en) 2003-07-30
SE518092C2 (sv) 2002-08-27
SE0003549L (sv) 2002-04-03

Similar Documents

Publication Publication Date Title
CN100379154C (zh) 非均匀性取样的限频信号的重建方法
Gopinath et al. Optimal wavelet representation of signals and the wavelet sampling theorem
Blu et al. Quantitative Fourier analysis of approximation techniques. I. Interpolators and projectors
Unser et al. A generalized sampling theory without band-limiting constraints
Wei et al. Lattices sampling and sampling rate conversion of multidimensional bandlimited signals in the linear canonical transform domain
Blu et al. Approximation error for quasi-interpolators and (multi-) wavelet expansions
Aldroubi et al. Wavelets on irregular grids with arbitrary dilation matrices and frame atoms for L2 (Rd)
Prendergast et al. Reconstruction of band-limited periodic nonuniformly sampled signals through multirate filter banks
Bui et al. Translation-invariant denoising using multiwavelets
Van De Ville et al. Isotropic polyharmonic B-splines: Scaling functions and wavelets
Daubechies et al. Factoring wavelet transforms into lifting steps
US5838377A (en) Video compressed circuit using recursive wavelet filtering
Blu et al. Wavelets, fractals, and radial basis functions
Evangelista et al. Discrete frequency warped wavelets: Theory and applications
Shao et al. Order tracking by discrete Gabor expansion
US5495432A (en) Apparatus and method for sampling rate conversion
Yadav et al. Orthogonal ramanujan sums, its properties, and applications in multiresolution analysis
JP2975691B2 (ja) 離散コサイン変換を用いる画像のサンプルレート変換のシステムおよび方法
Mitra et al. New methods of digital ladder realization
Mohan et al. A novel approach for satellite image resolution enhancement
Entezari et al. Quasi‐interpolation on the Body Centered Cubic Lattice
Abdelnour Symmetric wavelets dyadic sibling and dual frames
Marks Restoration of continuously sampled band-limited signals from aliased data
CN110243473B (zh) 干涉式高光谱仪器的光谱局部插值方法
Nguyen et al. Design of two-band critically sampled rational rate filter banks with multiple regularity orders and associated discrete wavelet transforms

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
ASS Succession or assignment of patent right

Owner name: INFINEON TECHNOLOGIES AG

Free format text: FORMER OWNER: TELEFONAKTIEBOLAGET LM ERICSSON (PUBL)

Effective date: 20041231

C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20041231

Address after: Munich, Germany

Applicant after: Infennian Technologies AG

Address before: Stockholm, Sweden

Applicant before: Telefonaktiebolaget LM Ericsson

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: 20080402

Termination date: 20160928

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