CN105849802B - 用于正交镜像滤波的方法和装置 - Google Patents

用于正交镜像滤波的方法和装置 Download PDF

Info

Publication number
CN105849802B
CN105849802B CN201480071503.6A CN201480071503A CN105849802B CN 105849802 B CN105849802 B CN 105849802B CN 201480071503 A CN201480071503 A CN 201480071503A CN 105849802 B CN105849802 B CN 105849802B
Authority
CN
China
Prior art keywords
samples
array
sample
filter
extracted
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
CN201480071503.6A
Other languages
English (en)
Other versions
CN105849802A (zh
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.)
Samsung Electronics Co Ltd
Original Assignee
Samsung Electronics Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Samsung Electronics Co Ltd filed Critical Samsung Electronics Co Ltd
Publication of CN105849802A publication Critical patent/CN105849802A/zh
Application granted granted Critical
Publication of CN105849802B publication Critical patent/CN105849802B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS TECHNIQUES OR SPEECH SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING TECHNIQUES; SPEECH OR AUDIO CODING OR DECODING
    • G10L19/00Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis
    • G10L19/02Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders
    • G10L19/0204Speech or audio signals analysis-synthesis techniques for redundancy reduction, e.g. in vocoders; Coding or decoding of speech or audio signals, using source filter models or psychoacoustic analysis using spectral analysis, e.g. transform vocoders or subband vocoders using subband decomposition
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F9/00Arrangements for program control, e.g. control units
    • G06F9/06Arrangements for program control, e.g. control units using stored programs, i.e. using an internal store of processing equipment to receive or retain programs
    • G06F9/30Arrangements for executing machine instructions, e.g. instruction decode
    • G06F9/34Addressing or accessing the instruction operand or the result ; Formation of operand address; Addressing modes
    • G06F9/355Indexed addressing
    • G06F9/3552Indexed addressing using wraparound, e.g. modulo or circular addressing
    • HELECTRICITY
    • H03ELECTRONIC CIRCUITRY
    • H03HIMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00Networks using digital techniques
    • H03H17/02Frequency selective networks
    • H03H17/0248Filters characterised by a particular frequency response or filtering method
    • H03H17/0264Filter sets with mutual related characteristics
    • H03H17/0272Quadrature mirror filters

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • Audiology, Speech & Language Pathology (AREA)
  • Human Computer Interaction (AREA)
  • Computational Linguistics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Compression, Expansion, Code Conversion, And Decoders (AREA)
  • Signal Processing For Digital Recording And Reproducing (AREA)

Abstract

一种执行正交镜像滤波器(QMF)合成滤波的方法包括:在包括调制的QMF子带采样的第一阵列中,在要丢弃的采样的位置处记录对应于当前时隙的新的采样。该方法还包括:从第一阵列中提取采样以移除在相邻子带之间的混叠,通过执行模运算确定对应于提取的采样的滤波器系数,以及通过使用提取的采样和滤波器系数合成移除混叠的时域采样。

Description

用于正交镜像滤波的方法和装置
技术领域
这里公开的一个或多个实施例涉及用于使用正交镜像滤波器(quadraturemirror filter,QMF)组执行分析滤波或合成滤波的方法和装置。
背景技术
数字音频编码方法不仅在移动设备的领域而且在家庭音频系统的领域都是非常重要的因素。在根据ISO MPEG音频标准的音频编码算法中,MPEG层3(MP3)可用于广播和多媒体内容。例如AAC,HE-AAC v1和HE-AAC v2是MP3的改进,它们以相对较少的比特率提供高品质的音乐。根据提供多声道的MPEG环绕编码方法,消费者可能不仅体验高品质的声音而且体验现场的和生动的多声道声音。最近,包括在杜比睿波(Dolby Pulse)和MS10中的HE-AAC被广泛用作家用音频编码方法。上述声音编码方法都使用基于QMF的子带编码技术。因为人类的听觉系统基于音频信号的频率识别声音,所述基于QMF的子带编码技术在音频和声音信号的压缩中非常有效。但是,基于QMF的子带编码需要大量的计算的负担以实现其效果。
发明内容
技术问题
基于QMF的子带编码需要大量的计算的负担以实现其效果。
解决方案
这里公开的一个或多个实施例包括用于减少正交镜像滤波器(QMF)滤波的复杂性并迅速执行QMF滤波的方法和装置。
有益效果
可使用一种编码方法,在其中,资源被有效地分配,使得在对人耳相对不敏感的高频率范围内的频带信号由相对小数量的比特表示,并且在对人耳相对非常敏感的低频率范围内的频带信号由相对大数量的比特表示,从而减少编码比特的数量。
附图说明
通过结合附图的以下各实施例的描述,这些和/或其他方面将变得明了和更容易理解,其中:
图1示出根据本公开的实施例的QMF滤波系统;
图2示意示出在根据本公开的实施例的QMF合成滤波方法和传统方法之间的比较;
图3是用于解释QMF合成滤波算法的流程图;
图4描述记录采样阵列的过程;
图5描述根据图4的算法提取的采样和滤波器系数;
图6描述根据本公开的实施例的QMF合成滤波方法;
图7描述根据本公开的实施例的QMF合成滤波方法;
图8描述根据本公开的实施例的记录采样阵列的过程;
图9描述根据本公开的实施例的从采样阵列提取采样的过程;
图10示出根据本公开的实施例的滤波器系数阵列;
图11描述根据本公开的实施例的选择滤波器系数的过程;
图12描述根据本公开的实施例提取的采样和滤波器系数;
图13是根据本公开的实施例的用于解释QMF合成滤波方法的流程图;
图14是根据本公开的实施例的用于解释记录采样阵列的方法的流程图;
图15是根据本公开的实施例的用于解释记录采样阵列的方法的流程图;
图16是根据本公开的实施例的用于解释从采样阵列提取采样的方法的流程图;
图17是根据本公开的实施例的用于解释从采样阵列提取采样的方法的流程图;
图18是根据本公开的实施例的用于解释选择滤波器系数的方法的流程图;
图19是根据本公开的实施例的用于解释合成时域采样的方法的流程图;以及
图20是根据本公开的实施例的合成滤波装置的框图。
具体实施方式
这里公开的一个或多个实施例包括用于减少正交镜像滤波器(QMF)滤波的复杂性并迅速执行QMF滤波的方法和装置。
额外的方面部分将在随后的描述中阐述,部分从描述中将明了,或者可以通过公开的实施例的实践而领悟。
根据本公开的一个或多个实施例,执行正交镜像滤波器(QMF)合成滤波的方法可以包括:在包括调制的QMF子带采样的第一阵列中在要丢弃的采样的位置处记录对应于当前时隙的新的采样;从第一阵列中提取采样以移除在相邻子带之间的混叠;通过使用模运算确定对应于提取的采样的滤波器系数;以及通过使用提取的采样和滤波器系数合成已移除混叠的时域采样。
对应于当前时隙的新的采样的记录位置可以通过模运算在先前时隙中的采样记录位置来确定。新的采样可以被记录在从先前时隙中的采样记录位置偏移N数量的新的采样的位置处。当在先前时隙中的采样记录位置是第一阵列的开始时,可以从第一阵列的最后第N位置起记录新的采样。在第一阵列中记录新的采样可以包括:通过逆变换变换到基带的K数量的复子带采样,产生对应于K数量的子带的N数量的实数采样;以及在要丢弃采样的位置处记录N数量的实数采样。
可以根据当前时隙是偶数还是奇数确定要提取的采样的位置。在当前时隙是偶数时从第一阵列提取的采样的位置与在当前时隙是奇数时从第一阵列提取的采样的位置可以相互排斥。当第一阵列划分为多个部分时,对于自然数n,在当前时隙是偶数时可以在第(4n-3)和第4n部分中提取采样,以及在当前时隙是奇数时可以在第(4n-2)和第(4n-1)部分中提取采样。
确定滤波器系数可以包括:通过使用当前时隙和合成滤波器的级(degree)执行模运算;以及基于模运算的结果从第二阵列选择滤波器系数。相同的滤波器系数在第二阵列中可以被记录两次,并且该相同的滤波器系数被记录的位置可以按照合成滤波器的级彼此分离。从第二阵列选择的滤波器系数的位置可以从自先前的时隙选择的滤波器系数的位置起移位1。
可以通过累计提取的采样与确定的滤波器系数的相乘的结果来合成该时域采样。可以根据时隙的变化递归地执行记录、提取、确定和合成。
根据本公开的一个或多个实施例,一种执行正交镜像滤波器(QMF)合成滤波的方法可以包括:通过根据时隙的变化在第一阵列中循环采样记录位置来在第一阵列中记录新的采样,考虑采样记录位置的循环来确定要从第一阵列中提取的采样的位置,以及基于在确定的位置提取的采样合成对于当前时隙的时域采样。
在第一阵列中记录新的采样可以包括:确定在先前时隙在第一阵列中记录的N数量的采样的领先(leading)采样的位置是否是第一阵列的开始,并且如果在先前时隙的领先采样的位置是第一阵列的开始,则从第一阵列的最后位置的第N位置起记录新的采样。如果在先前时隙的领先采样的位置不是第一阵列的开始,则可以从自领先采样的位置减少N的位置起记录新的采样。
在确定要从第一阵列中提取的采样的位置中,位置这样确定使得从与在先前时隙中提取的采样的位置不同的位置中提取采样。在确定从第一阵列中提取的采样的位置中,要提取的采样的位置可以根据当前时隙是偶数还是奇数来确定。
该方法可以还包括选择对应于从记录滤波器系数的第二阵列提取的采样的系数。在通过按照预定间隔划分提取的采样获得的M数量的组的每个中,分别对应于第n采样和第n+1采样的第一滤波器系数和第二滤波器系数在第二阵列中可以按2M彼此分离。在第二阵列中选择系数中,可以在一位置选择对于M数量的组的第m组的第n采样的系数,该位置通过将2*M*(n-1)+m-1和当前时隙被M除的余数求和获得。
根据本公开的一个或多个实施例,一种执行正交镜像滤波器(QMF)合成滤波的方法可以包括:将经调制的子带采样记录在具有大小M*N的第一阵列中根据时隙的改变循环的具有大小N的记录部分中,在具有大小M*N的第二阵列中选择预定(例如M*N/2)数量的滤波器系数,其中作为一组具有连续的M数量的滤波器系数的预定(如N)数量的系数组的每个被连续记录两次,以及通过使用从第一阵列所提取的预定(例如M*N/2)数量的采样和从第二阵列中选择的(如M*N/2)数量的滤波器系数合成对于当前时隙的时域采样。例如,M和N可以是正数或整数,大于零。例如,M和N可以是计数号码。
M可以对应于合成滤波器的级而N可以是QMF子带数的量两倍。
在经调制的子带采样的记录中,如果在先前时隙中在第一阵列记录的采样的领先采样的位置i是0,则经调制的子带采样可以从位置[(M-1*N)]起记录,以及如果在先前时隙中在第一阵列记录的采样的领先采样的位置i不是0,则经调制的子带采样可以从位置[i-N]起记录。
在滤波器系数的选择中,可以在第二阵列的位置[m+2*M*n+MOD(l,M)]选择滤波器系数,其中m是大于或等于0并且小于M的整数,n是大于或等于0并且小于N/2的整数,1是当前时隙,以及MOD是模运算。
该方法还可以包括从第一阵列提取(M*N/2)数量的采样,其中,如果当前时隙是奇数,则可以在第一阵列的位置[2*N*n+k+N/2]提取(M*N/2)数量的采样,其中,k是大于或等于0且小于N的整数,并且n是大于或等于0且小于N/2的整数。
该方法可以进一步包括从所述第一阵列中提取采样和在第三阵列中存储提取的采样,其中第二阵列是第三阵列和第一阵列的大小的两倍。
根据本公开的一个或多个实施例,一种非暂时性计算机可读记录介质可以在其上记录用于执行任何在此公开的方法的一个或多个程序。
根据本公开的一个或多个实施例,用于执行正交镜像滤波器(QMF)合成滤波的装置可以包括存储器,存储第一阵列和第二阵列,其中在第一阵列中记录调制的子带信号的采样,以及在第二阵列中记录滤波器系数;调制单元(调制器),将子带信号的采样调制到基带的K数量的子带并且在第一阵列中存储调制的采样;和合成滤波器,从第一阵列提取采样以移除在相邻子带之间的混叠,通过使用模运算从第二阵列选择对应于提取的采样的滤波器系数,以及通过使用提取的采样和选择的滤波器系数合成时域采样,其中调制单元通过根据时隙的变化在第一阵列中循环采样记录位置来记录对应于当前时隙的新的采样。
该合成滤波器可以根据当前时隙是偶数还是奇数确定要从第一阵列提取的采样的位置。相同的滤波器系数在第二阵列中可以被记录两次,并且该相同的滤波器系数被记录的位置可以按照合成滤波器的级彼此分离。该装置可以还包括解码器,其解码编码的多声道比特流并且输出子带的信号。
根据本公开的一个或多个实施例,一种执行正交镜像滤波器(QMF)合成滤波的方法可以包括:在调制的QMF子带采样的第一阵列的第一部分的位置记录对应于当前时隙的新的采样,所述第一部分的位置包括要丢弃的采样,从第一阵列中提取采样,确定对应于提取的采样的滤波器系数,并通过使用所提取的采样和滤波器系数合成时域采样。新的采样被记录在第一部分的位置,并且例如同时地从该第一部分丢弃采样,无需转移位于限定第一阵列的第一阵列的剩余部分的采样。
随着时隙的值增加,新的采样可以记录在第一阵列的第二部分中,第二部分与第三部分相邻,在第三部分中在先于新的采样在第二部分中记录的时隙的时隙中记录采样。
可以通过选择在第一时隙的部分来从第一阵列中提取采样,第一时隙中的部分不同于在第二时隙中选择的部分(从中提取来自第一阵列的采样)。
滤波器系数可以这样确定,对于第一时隙通过从第一组滤波器系数中选择滤波器系数,而对于跟随第一时隙的第二时隙,通过从第二组滤波器系数中选择滤波器系数,并从第一组滤波器系数选择部分的滤波器系数。第一组滤波器系数和第二组滤波器系数可以相同。
发明模式
现在将详细地参考各实施例,其示例在附图中示出,其中相似的参考标号始终指代相似的元件。在这方面,所公开的实施例可具有不同的形式并且不应被解释为限于在此阐述的描述。因此,本实施例仅仅如下通过参照附图来描述,以解释本公开的各方面。如这里所使用的,术语“和/或”包括一个或多个相关所列的项目的任意组合和所有组合。诸如“至少一个”的表达在元素列表的前面时,修该整个列表的元素,而不修改该列表的单个元素。
图1示出根据本公开的实施例的QMF滤波系统1。参照图1,QMF滤波系统1可以包括分析滤波器组101、编码器102、解码器111以及合成滤波器组112。分析滤波器组101和编码器102构成分析滤波设备10,而解码器111和合成滤波器组112构成合成滤波设备11。分析滤波设备10和合成滤波设备11可以是独立的设备。
分析滤波器组101可对输入的时域的语音信号X执行分析滤波,并输出K数量的子带信号S0~K-1 103。编码器102可编码子带信号S0~K-1的一个或多个,并输出比特流。在这样做时,可使用这样的编码方法,在其中资源被有效地分配,使得在对人耳相对不敏感的高频率范围内的频带信号由相对小数量的比特表示,并且在对人耳相对非常敏感的低频率范围内的频带信号由相对大数量的比特表示,从而减少编码比特的数量。
解码器111可接收经编码的比特流,并进行编码过程的逆过程,并由此产生子带信号S'0~K-1 113,其与编码之前的子带信号S0~K-1非常类似。子带信号S0~K-1与信号S'0~K-1不相同的原因是,在通过编码器102的编码过程中子带信号S0~K-1通过除去对人的听力系统不重要的信号分量来进行压缩。合成滤波器组112可以对每个子带信号S'0~K-1执行合成滤波以将每个子带信号S'0~K-1变换为可由人听到的时域的信号Y。编码器102和解码器111可以根据实施的编解码器的类型采用各种方法。在下面的描述中,将主要探讨QMF分析或合成滤波过程。
图2示意示出在根据本公开的实施例的QMF合成滤波方法和传统方法之间的比较。由于图2示意性示出与传统方法比较的概念,本公开的具体实施例可以通过参考以下描述来理解。
在左边的合成滤波器组20示出了传统的合成滤波方法,而在右边的合成滤波器组21示出根据本公开的实施例的合成滤波方法。在根据合成滤波器组20的合成滤波方法中,调制201,滤波202和存储器更新203可以被递归执行。
在根据依据本公开的实施例的右边的合成滤波器组21的合成滤波方法中,虽然调制211和滤波212被递归执行,但是存储器更新可以不执行。为了省略存储器更新,调制201和滤波202可以被新定义为调制211和滤波212,这将在下面进行详细说明。为了检查根据在传统合成滤波方法中的存储器更新203的计算量,传统合成滤波器组20的操作首先被讨论如下。
在调制201的过程中,被变换为基带的每个子带信号S'0~K-1被频带偏移(band-shifted)到原始频带。每个子带信号S'0~K-1的调制采样被记录在采样阵列V,其由通过等式1表示。
【等式1】
在等式1,K表示子带的总数,l表示时隙索引,n表示采样数目,N表示相对于时隙l的新调制采样的总数,而S'k,l表示在k子带信号中时隙l的子带采样。例如,N可以等于K x2。为了防止在术语上混淆,定义有术语“子带采样”、“子采样”或“子带信号的采样”都由S'k,l指示,并且术语“采样”表示或对应于根据等式1调制的子带采样V[n]。除非另有定义,在以下描述中上述变量以相同的方式被使用。
fk,n是逆离散余弦变换(IDCT)函数并且可以根据编解码器的类型使用各种方法。例如,fk,n可以由等式2表示。
【等式2】
滤波202可以表示或对应于带通滤波,其移除在采样阵列V中相邻子带之间的混叠,同时只留下k子带信号。带通滤波可以通过等式3和4来定义。
【等式3】
等式3表示或对应于从采样阵列V中提取采样的过程。当采样阵列V中N/2数量的采样为一组时,根据等式3从第一到第四采样组提取N数量的采样。按类似于等式3的方式从采样阵列V的其他采样组中提取采样。从采样阵列V中提取的采样被记录在提取阵列g中。
【等式4】
其中0≤n<N/2
等式4表示或对应于累积所提取的采样和滤波系数的相加结果,并相对于时隙1合成时域采样的过程。在等式4中,c表示存储滤波器系数的滤波器系数阵列。滤波器系数可以被称为窗口系数。M表示滤波202的级。例如,滤波202可通过使用有限脉冲响应(FIR)滤波器来执行,而M可以表示FIR滤波器的级。当N是128时,时域采样的总共64个可以被合成。合成的结果被存储在阵列y中。
当相对于该时隙l的滤波202被完成时,存储器更新203被执行以执行相对于下一个时隙l+1的调制201和滤波202。存储器更新203可以由等式5定义。
【等式5】
V[n]=V[n-N],其中N≤n<N·M
根据等式5,在采样阵列V中N*(M-1)数量的采样被移位N。这里,符号“*”表示乘法运算。读和写每个对于一个采样的一次移位执行一次。根据存储器更新203,在更新之前记录在位置[(M-1)N]到[MN-1]的采样被丢弃,而记录在位置[(M-2)N]到[(M-1)N-1]的采样被移位到该丢弃采样的位置上。在以下的说明中,为方便起见,在位置[0]到[N]被表示为[0~N]。
当存储器更新203完成时,相对于下一个时隙l+1调制201和滤波202被递归执行。
图3是用于解释QMF合成滤波算法的流程图。图3示出当滤波器的级是10(M=10)和子带的数量是64以及新的采样数为128(K=64,其中N=2*K=128)时的QMF合成滤波算法。参照图3,操作A301、A302-A304和A305分别表示调制201,滤波202和存储器更新203。
在操作A301,新的采样被记录在V[0~127]。X[k][l]对应于等式1的S'k,l。换言之,子带采样X[k][l]表示到第k子带信号的第l时隙的采样。在操作A301,通过调制X[k][l]产生新的采样,并记录在V[0~127]中。
在其中参考图4的操作A301中记录新的采样的部分中,当l=0时,新的采样被记录在部分400中,当l=1时,新的采样被记录在部分401中,而当l=M时,新的采样被记录在部分402中。因此,新的采样总是被记录在V[0~127]中。
如图4所示,根据传统技术,在采样阵列V中,其中记录新的采样的部分40和其中丢弃现有的采样的部分41(也见部分410)不管时隙的改变总被维持。这是因为存在于部分42(也见部分420)的采样被移位N。
在操作A302,六百四十(640)个采样被提取并记录在提取阵列g。操作A302是相应于等式3的过程。
参照图4,从采样阵列V提取的采样用阴影示出。换句话说,当l=0,M数量(例如,十)的采样组V2M-1,V2M-4,V2M-5,V2M-8,V2M-9,...,V4,V3和V0被提取。每个采样组包括64个采样。另外,当l=1时,M数量的采样组V2M+1,V2M-2,V2M-3,V2M-6,V2M-7,...,V6,V5及V2被提取。因此,当时隙改变时,提取的采样的值被改变,但提取的采样的位置[0~N/2-1][3N/2~5N/2-1]...[MN-N/2~MN-1]没有改变。因此,采样总在相同部分提取而不管时隙如何。
在操作中A303,提取的采样和滤波系数彼此相乘并且其结果被记录在滤波器阵列w中。在操作A304,六十四(64)个时域采样通过累加滤波器阵列W而输出。操作A303和A304是对应于等式4的过程。
参考图4和5讨论操作A303和A304。图5示例性示出等式4的结果。图5示出等式4的部分结果,按顺序从(l=0,n=0,m=0)开始到(l=0,n=0,m=1)...(l=0,n=0,m=9),(l=0,n=1,m=0),(l=0,n=1,m=1)…(l=0,n=63,m=9)和(l=1,n=0,m=0)。例如,在等式4中,m可以是最下面例程的变量而l可以是最上面例程的变量。等式4可以通过参考图5而容易被理解,在图5中当n相同时,滤波器系数的位置根据m增加+1。此外,当m相同时,滤波器系数根据n中的增加而增加+M。如可以看到的,在时隙l的变化不影响滤波器系数的位置。
在图5所示的结果参照图4描述。阴影的系数组是从采样阵列V中提取的采样,并且表示在等式4中的提取阵列g。图4的箭头表示滤波器系数。在同一高度的水平方向中箭头的运动表示值m的变化。另外,在垂直方向中箭头的移动表示值n的变化。可以看出,尽管时隙由l=0改变至l=M,乘以提取的采样的滤波器系数的位置被无变化地保持。
在操作中A305,记录在位置[1151~128]的采样被顺序地移位。参照图4,当时隙从l=0改变至l=1时,存在于采样阵列V[(M-1)N~MN-1]的两个采样组V0和V1被丢弃。根据移位,存在于采样阵列V[(M-2)N~(M-1)N-1]的两个采样组V2和V3被记录在V[(M-1)N~MN-1]。
根据操作A305,丢弃的采样的位置总是固定为相对于全部时隙变化的常数。换言之,位于[(M-1)N~MN-1]的采样总是被丢弃。
根据操作A305,(M-1)N次的存储器读取和写入被执行以移位(M-1)N数量的采样。例如,当在处理器中分别按4个周期和6个周期执行读和写时,上面的过程按10(M-1)N个周期执行。对于实际的标准编解码器,如HE-AAC v2,杜比睿波或MS10DDT,因为M是10,N是128,用于对于每个声道的存储器更新过程的周期是368640(10*9*128*32)。因此,需要最小17Mcycle用于立体声的情况,并且需要51.84Mcycle用于5.1环绕声道的情况。来看上述计算量对总计算量的比率,对于HE-AAC v2的解码器,例如,用于滤波器的存储器更新的周期占总计算量的约32%以及占只有QMF滤波处理的约55%。
如下所述,模寻址方法用于在存储器中存储输入数据,以使得在公开的实施例中存储器更新过程可以省略,另外,滤波处理被重新配置或修改,以适应上述方法使得存储器更新的复杂性被除去。
图6描述了根据本公开的实施例的QMF合成滤波方法。除非另外提及,下面描述的合成滤波方法的实施例可以理解为根据本公开的实施例的合成滤波设备11的操作。
首先,当时隙l为0(l=0)时,对应于为0的时隙l(l=0)的N数量的新的采样被记录在采样阵列V的部分600中。当时隙l改变为(l=1)时,位于部分610的采样被丢弃,并且同时N数量的新的采样被记录在部分601中,其对应于部分610并且占据所丢弃的采样的相同位置。换句话说,位于位置[M(N-1)~MN-1]的采样组V0,1被丢弃,而同时采样组V2M,2M+1记录在位置[M(N-1)~MN-1]。因此,可以理解根据时隙改变,丢弃的采样被新的采样取代,该新的采样对应于等于1的当前时隙l(l=1)。换句话说,按照覆写方法对应于当前时隙的新的采样被记录在丢弃的采样的位置处。不同于传统的方法,在本实施例中,位于部分620的采样没有移位。当在l=0和l=1的采样阵列V相互比较时,可以看到采样组V2到V2M-1被维持不变。
在该示例实施例中,由于新的采样被记录而没有移位采样,故可以省去根据等式5和操作305的存储器更新203。
当时隙l为2(l=2)时,存在于部分611的采样组V2,3被丢弃,同时采样组V2M+2,2M+3被记录。当时隙l为M(l=M)时,存在于部分603的采样组被丢弃,同时采样组V4M-1,4M-2被记录在部分602中,其对应于部分611。
根据依据示例实施例的采样记录方法,在采样阵列V的采样记录部分根据时隙的改变来循环。由于记录部分被循环,本公开的采样阵列V操作上类似于图8的圆形阵列80。
参照图8的圆形阵列80,圆形阵列80的大小对应于MN并且可以划分为总共M数量的部分(如10)。每个部分的大小为N(如128)。圆形阵列80的开始位置[0]和最后位置[MN](如1280)被认为是相同的位置。当时隙l为0(l=0)时,N数量的采样被记录在部分811。当时隙l为1(l=1)时,N数量的采样被记录在部分812。当时隙l为2(l=2)时,N数量的采样被记录在部分813。因此,随着时隙改变,记录部分按照虚线指示的方向循环。在示例实施例中,领先采样的记录位置可以用于识别记录部分的方法。领先采样表示在一个记录部分中记录的采样中首先记录的采样。例如,当l=0时,领先采样821的记录位置是[0]并且在领先采样821之后的(N-1)数量的采样(如127)被顺序记录在记录位置[1~N-1](如0,然后1,2,3到127)。其中记录采样的方向通过由实线指示的箭头来指示。因此,其中领先采样的记录位置被循环的方向和其中记录N数量的采样的方向可以彼此相对。
根据示例实施例,对应于当前时隙的新采样的记录位置或记录部分可以根据模寻址(modulo addressing)来确定。同样,可根据模寻址来确定将根据时隙的改变被丢弃的采样的位置。这是因为,要丢弃的采样的位置和新的采样的位置根据时隙的改变被循环。
在当前时隙中要记录的新的采样中领先采样的记录位置il可以由等式6表示。
【等式6】
il-1是在先前时隙记录的采样中的领先采样的记录位置。虽然等式6表示领先采样的记录位置il为分支语句,这基本上是指模寻址。换句话说,这是因为,当相对于前一个时隙的记录位置il是采样阵列V的开始位置[0]时,相对于当前时隙的记录位置il回到最后的第N位置。当相对于前一个时隙的记录位置不是开始位置[0]时,N数量的新采样是从自位置il-1减少N的位置开始记录。当等式6由模函数表达时,它可通过等式7来表示,但本公开不限于此。
【等式7】
il=MOD(M-MOD(l,M),M)×N
=MOD(MN-MOD(il-1+N,MN),MN)
=(M-1-MOD(l-1,M))×N
新采样的记录位置可通过不同于领先采样的记录位置的方法来表达。例如,新的采样的记录位置可按照例如使用最后的采样、第n个采样、记录部分的开始和结束位置、记录部分的标识符等(它们不对应于领先采样)的各种方法来表示。
根据示例实施例的调制过程可以由等式8来表示。
【等式8】
其中0≤n<N
参考等式8,N数量的新采样可以从通过等式6或7确定的记录位置[il]起顺序记录。新的采样表示或对应于通过等式2的函数调制的子带采样S′k,l。S′k,l例如可以是复子带采样。同样,仅等式8计算的值的实数可以存储在采样阵列V。换句话说,通过等式2和8变换到基带的K数量的复子带采样被逆变换并且产生对应于K数量的子带的N数量的实数采样。通过等式6,7和8可以将产生的采样记录在采样阵列V中。
可以提取采样阵列V的采样以除去混叠。然而,由于在采样阵列V中采样的记录位置和顺序被改变,传统的采样提取方法可能不适用于示例实施例。根据示例实施例,从采样阵列V提取采样中,考虑到根据时隙的改变的采样记录位置的循环来提取采样。
根据图4的传统方法和图6的公开的示例实施例之间的阴影部分的比较,可以看到,虽然在图4和6中提取的采样组的数量相同,但是顺序和提取位置是不同的。当时隙1为0(1=0)时,在图4和6中采样的提取顺序与位置相同。然而,当时隙l为1至M-1(l=1至M-1)时,在图4和6中采样的提取顺序与位置相当不同。
根据示例实施例,要从采样阵列提取的采样的位置可以根据当前时隙是否为奇数还是偶数来确定。在图6,当时隙是偶数时提取的采样的位置和当时隙是奇数时提取的采样的位置是排斥的。每当时隙改变时要提取的采样的位置被反转。
当采样阵列V通过采样的N/2的间隔划分时,如果当前时隙是偶数,则采样可以在第1,第4,第5,第8,第9,第12,第13,…,和第2M部分中提取。如果当前时隙是奇数,则采样可以在第2,第3,第6,第7,第10,第11,…,第(2M-2)和第(2M-1)部分中提取。换句话说,在当前时隙是偶数时,可以在第(4n-3)和第4n部分中提取采样。在当前时隙是奇数时,可以在第(4n-2)和第(4n-1)部分中提取采样。n表示小于或等于M/2的自然数。
等式9表示根据示例实施例提取采样的过程。
【等式9】
如果l=偶数;
否则,l=奇数;
g[Nn+k]=V[2Nn+k+N/2],其中0≤k<N
根据等式9,从采样阵列V提取(M*N/2)数量的采样。图12示出根据等式9的采样提取的部分结果。图9反映在图6的图示、等式9和图12的结果之间的关系。当对等式9应用情况(l=0,n=1)时,可以从采样阵列提取(1)V[0~N/2-1],(2)V[3N/2~4N/2-1],(3)V[4N/2~5N/2-1],和(4)V[7N/2~8N/2-1]。如图9所示,提取的采样(1)~(4)可以顺序记录在提取阵列g中。
当根据示例实施例的提取阵列是g1而根据常规方法的提取阵列是g2时,除了其中时隙是M的整数倍的情况,g1[n]≠g2[n]。相应地,乘以g1[n]的滤波器系数不同于乘以g2[n]的滤波器系数。
这是因为,不同于g2,考虑到采样阵列的采样记录位置的循环而产生提取阵列g1。因此,根据示例实施例,记录对应于提取阵列g1的滤波器系数的滤波器系数阵列c′通过新的方法来定义,其不同于传统的滤波器系数阵列c。
等式10可定义根据示例实施例的滤波器系数阵列c′。
【等式10】
c′[m+2Mn]=c′[m+2Mn+M]=c[m+Mn]其中0≤n<N/2
其中0≤m<M
参考等式10和图10来描述根据示例实施例的滤波器系数阵列。图10示出当M=10和N=128时的传统的滤波器系数阵列c(t002)和滤波器系数阵列c′(t001)两者。
滤波器系数阵列c(t002)的大小可以与提取阵列g的大小MN/2相同。然而,滤波器系数阵列c′(t001)的大小是MN,这是提取阵列g和滤波器系数阵列c(t002)每个的大小的两倍。
滤波器系数阵列c′(t001)的大小是滤波器系数阵列c(t002)的大小的两倍的原因是滤波器系数阵列c(t002)的滤波器系数在滤波器系数阵列c′(t001)中记录两次。相同的滤波器系数在滤波器系数阵列c′(t001)中记录两次并且用于记录相同的滤波器系数的位置分离合成滤波器的级M。
当滤波器系数阵列c(t002)划分为N/2数量的滤波器系数组时,第一系数组(t01)在滤波器系数阵列c′(t001)中被连续地记录两次(t02)。同样,来自第二系数组到最后系数组的各组在滤波器系数阵列c′(t001)中被连续记录两次,并依此类推。
根据示例实施例,可以根据模寻址确定对应于提取的采样的滤波器系数。等式11表示根据模寻址选择滤波器系数的方法。
【等式11】
c′[m+2Mn+MOD(l,M)]其中0≤n<N/2其中0≤m<M-l
根据等式11,滤波器系数在位置[m+2MN+MOD(1,M)]处被选择,[m+2MN+MOD(1,M)]通过将2*M*(n-1)+m-和当前时隙除以M的余数求和获得。MOD(1,M)表示或对应于使用合成滤波器的级M和当前时隙1的模运算。
图11示出当在等式11中M=10和N=128时在滤波器系数阵列c′中选择滤波器系数的示例过程。参照图11,当n=0时,滤波器系数相对于全部时隙1在部分t100中被选择。当n=1时,滤波器系数相对于全部时隙1在部分t101中被选择。但是,当时隙1增大1时,选择的滤波器系数的位置增加+1。
首先,当1=0且n=0时,对于0≤m<M选择c′[0~9]。这里,c′[0~9]=c[0~9]。当n增加1时,对于0≤m<M选择c′[20~29]。换句话说,当n增加1时,滤波器系数的位置增加[2M]。这是因为相同的滤波器系数在滤波器系数阵列c′中按照间隔M记录两次。这里,c′[20~29]=c[10~19]。
当1=0时,所选择的滤波器系数可以是c[0~9]、c[10~19],这和当1=0时选择的滤波器系数相同。这是因为,当1=0时,在采样阵列V中采样的记录位置和提取位置与根据传统方法的那些相同。然而,随着时隙1被改变,在公开的实施例中,与根据传统方法的那些不同,改变了在采样阵列V中采样的记录位置和提取位置,因此与传统方法不同的滤波器系数被选择。
当l=1且n=0时,对于0≤m<M选择c′[1~10]。这里,c′[1~9]=c[1~9]以及c′[10]=c[0]。当n增加1时,对于0≤m<M选择c′[21~30]。这里,c′[21~29]=c[11~19]以及c′[20]=c[10]。如图4所示,根据传统方法,按照c[0~9],c[10~19]…的顺序选择滤波器系数。根据本公开的实施例,滤波器系数的选择顺序变为c[1~9],c[0],c[11~19],c[10]…。当1=2时,滤波器系数的选择顺序是c[2~9],c[0~1],c[12~19],c[10~11]。换句话说,从在先前时隙选择的滤波器系数的位置移位1的位置处选择滤波器系数。
这是因为,根据公开的实施例,随着在采样阵列V中采样的记录位置和提取位置根据模寻址循环改变,滤波器系数的选择位置根据模寻址循环改变。
等式12表示通过累计由等式9提取的采样与由等式11选择的滤波器系数的相乘的结果的合成时域采样的过程。
【等式12】
其中0≤n<N/2
再参考图6,根据等式11,提取采样的位置和选择的滤波器系数被示出。
当时隙1=0改变为l=1时,采样对于n=0和m=0在V[N/2]提取,以及记录在g[1],然后g[1]乘以c'[1]。采样对于n=0和m=1在V[2N/2]提取,以及记录在g[N/2],然后g[N/2]乘以c'[2]。可以重复相同的过程直到n=0和m=9,然后n改变为1(n=1)。采样对于n=1和m=0可以在V[N/2+1]提取,以及记录在g[2],然后g[2]乘以c'[2M+1]。当相同的过程被重复直到n=64和m=9时,时隙l改变为2(l=2)。
图6所示的在合成滤波过程中对于变量n产生滤波器系数可以如图7中所示。图7直观地示出了循环根据时隙的改变选择的滤波器系数的位置。
上面的说明描述了根据公开的示例实施例的在合成滤波设备11中执行的QMF合成滤波方法。根据本公开的实施例的上述等式的至少一部分可以通过如图12的先前计算的查找表来体现。例如,关于根据时隙改变的循环采样的记录位置、提取采样的位置和所选择的滤波器系数的位置的信息的至少一部分可以作为查找表存储在合成滤波设备11中并随后使用。
图13是根据本公开的实施例的用于解释QMF合成滤波方法的流程图。上述实施例的任何多余的描述将在此省略。因此,可以参考上述说明而没有任何特殊的评论。在以下的说明中,第一、第二和第三阵列分别与采样阵列V、滤波器系数阵列c'和提取阵列g(对应或)具有相同的含义。
参照图13,合成滤波设备11可以在第一阵列中记录新的采样(A1305)。合成滤波设备11可以将对应于当前时隙l的新的采样记录在第一阵列中要被丢弃的采样的位置处,该第一阵列包括调制的QMF子带采样。
合成滤波设备11可以记录新的采样,同时根据时隙l的改变循环在第一阵列中的采样记录位置il。例如,合成滤波设备11可以如等式6到8所示模计算在先前的时隙l-1的采样记录位置il-1,并确定对应于当前时隙l的新采样的记录位置il
根据等式6,合成滤波设备11可以在从先前时隙l-1中的采样记录位置il-1偏移N数量的新的采样的位置处记录新的采样。然而,当先前时隙l-1中的采样记录位置il-1是第一阵列的开始位置[0]时,从第一阵列的最后第N位置[(M-1)N]起记录新的采样。
下面具体参考图14描述操作A1305。
合成滤波设备11可确定采样记录位置il-1是否是第一阵列的开始位置[0](A1405)。根据实施例,il-1可以是在先前时隙i-1中记录在第一阵列中的N数量的采样的领先采样的位置。合成滤波设备11可确定在先前时隙的领先采样的位置il-1是否是第一阵列的开始位置[0]。
如果在先前时隙的领先采样的位置il-1是第一阵列的开始位置[0],则合成滤波设备11可将当前时隙的采样记录位置il更新为(M-1)*N(即,il=(M-1)*N)(A1415)。
如果在先前时隙的领先采样的位置il-1不是第一阵列的开始位置[0],则合成滤波设备11可将当前时隙的采样记录位置il更新为il-1-N(即,il=il-1-N)(A1410)。
根据等式8,合成滤波设备11可以逆变换被变换到基带的K数量的复子带采样,以产生对应于K数量的子带的N数量的实数采样(A1420)。合成滤波设备11可以在要丢弃的采样的位置上记录N数量的实数采样(A1425)。操作A1420可以早于操作A1405进行。可替代地,操作A1420可以与操作A1405同时进行,或在操作A1405之后进行。
如上所述,当M对应于合成滤波器的级而N对应于加倍的QMF子带的数量时,合成滤波设备11可以将经调制的子带采样记录在具有大小M*N的第一阵列中根据时隙的改变循环的具有大小N的记录部分中。
返回参考图13,合成滤波设备11可以从第一阵列提取采样(A1310)。例如,采样可以被提取以除去相邻子带之间的混叠。合成滤波设备11可通过在第一阵列中循环采样记录位置来确定要从第一阵列提取的采样的位置。要被提取的采样的数目可以是M*N/2,它是第一阵列的大小的一半。
合成滤波设备11可确定对应于提取的采样的滤波器系数(A1315)。合成滤波设备11例如可以通过使用模操作确定滤波器系数。合成滤波设备11可以在具有大小M*N的第二阵列中选择(M*N/2)数量的滤波器系数,其中作为一组具有连续的M数量的滤波器系数的N数量的系数组的每个被连续记录两次。
合成滤波设备11可以通过使用所提取的采样和滤波器系数合成时域采样(A1320)。其中混叠被除去的合成的时域采样可以是对应于当前时隙的(N/2)数量的采样。合成滤波设备11可以根据时隙的改变递归地执行上述操作。
图15是根据本公开的实施例的用于解释记录采样阵列的方法的流程图。上述实施例的任何多余的描述将在此省略。因此,可以参考上述说明而没有任何特殊的评论。同样,本公开不局限于图15的实施例并且记录采样阵列的其他方法可以落入主张的主题的范围内。
首先,合成滤波设备11可以初始化等式8的各种变量(A1505)。例如,时隙索引l、采样索引n、子带索引k以及记录位置索引都被设置为0(l=0,n=0,k=0,i=0)(A1505)。术语“temp”可以指一个变量,用于临时存储在调制过程的中间输出的值,其也可以设置为0(temp=0)。
合成滤波设备11可以根据等式8执行操作使得temp=temp+S'k,l*fk,n(A1510)。合成滤波设备11确定是否k<K-1(A1515),并将k增加1直到k=K-1(A1520),从而执行操作使得temp=temp+S'k,l*fk,n。
当k=K-1,(即,k不小于K-1)时,合成滤波设备11可确定记录位置索引i是否是0(i=0)(A1525)。如果i≠0,则合成滤波设备11可设置记录位置索引i为i-N(i=i-N)(A1530)。如果i=0,则记录位置索引i为(M-1)*N(i=M-1)*N)(A1535)。
合成滤波设备11可以在V[i+n]中记录temp值(临时值)并且初始化temp(如temp=0)(A1540)。
合成滤波设备11可确定采样索引n是否小于N(n<N)(A1545)。合成滤波设备11可以将n增加1直到n=N并且初始化k(例如n=n+1,k=0),并且可以对应于下一采样索引重复调制操作(A1550)。
当n=N时,(即,n不小于N),合成滤波设备11可确定当前时隙l是否小于最大时隙L(A1560)。如果l<L,则合成滤波设备11改变当前时隙l至下一个时隙,初始化n和k(例如,n=0,k=0),并且递归地执行相对于下一时隙的调制操作A1510直到l=L(A1555)。当l=L时,(即,l大于L-1),合成滤波设备11可确定记录采样阵列的方法完成。
图16是根据本公开的实施例的用于解释从采样阵列提取采样的方法的流程图。上述实施例的任何多余的描述将在此省略。因此,可以参考上述说明而没有任何特殊的评论。同样,本公开不局限于图16的实施例并且从采样阵列提取采样的其他方法可以落入主张的主题的范围内。
合成滤波设备11可确定当前时隙是偶数还是奇数(A1605)。在当前时隙是奇数时,合成滤波设备11可以提取在第一阵列的第(4n-2)和第(4n-1)部分的采样(A1610)。在当前时隙是偶数时,合成滤波设备11可以提取在第一阵列的第(4n-3)和第4n部分的采样(A1615)。例如,合成滤波设备11可以根据等式9提取采样。根据等式9,提取(M*N/2)数量的采样。在当前时隙是奇数时,在位置[2*N*n+k+N/2]处提取采样。如上所述,k是大于或等于0且小于N的整数,并且n是大于或等于0且小于N/2的整数。
在当前时隙是偶数时从第一阵列中提取采样的位置,以及在当前时隙是奇数时从第一阵列中提取采样的位置可以是相互排斥的。例如,合成滤波设备11可确定从第一阵列提取的采样的位置,使得采样可以在不同于在先前时隙提取采样的位置的位置处提取。
合成滤波设备11可以在第三阵列记录所提取的采样(A1620)。
图17是根据本公开的实施例的用于解释从采样阵列提取采样的方法的流程图。上述实施例的任何多余的描述将在此省略。因此,可以参考上述说明而没有任何特殊的评论。同样,本公开不局限于图17的实施例并且从采样阵列提取采样的其他方法可以落入主张的主题的范围内。
合成滤波设备11可以根据等式9初始化各变量(A1705)。例如,合成滤波设备11可以设置时隙索引l、采样索引n和子带索引k的每个为0。
合成滤波设备11可确定l是否是偶数(A1710)。如果l不是偶数,则合成滤波设备11可以将k加1(A1725),直到k=N(即,k不小于N)(A1720),并且可以根据等式提取采样从而g[Nn+k]=V[2Nn+k+N/2](A1715)。
当k=N(即,k不小于N)时,合成滤波设备11可确定是否n<N/2(A1745)。合成滤波设备11可以将n增加1直至n=N/2,初始化k,以及递归地执行操作A1715中的采样提取(A1755)。
当n等于N/2(即,n不小于N/2),合成滤波设备11可确定时隙是否小于最大时隙L(A1790)。合成滤波设备11增加时隙直到l=L,并且初始化n和k(例如,n=0,k=0),以返回到操作A1710(A1760)。当l=L时,(即,l大于L-1),合成滤波设备11可确定从采样阵列中提取采样的方法完成。
当l被确定为偶数时,合成滤波设备11可确定是否k<N/2(A1730)。合成滤波设备11可以将k增加1直到k=N/2(即,k不大于N/2),(A1740),并且可以根据等式提取采样,使得g[Nn+K]=V[2Nn+k+N](A1735)。当k=N/2(即,k不小于N/2)时,合成滤波设备11可确定是否k<N(A1770)。合成滤波设备11可以将k增加1直到k=N(即,k不小于N),(A1775),并且可以根据等式提取采样,使得g[Nn+k]=V[2Nn+k+N/2](A1765)。当k=N时,(即,k不小于N)时,合成滤波设备11可确定是否n<N/2(A1780)。合成滤波设备11将n增加1直至n=N/2,(即,n不小于N/2),并可初始化k(例如k=0),以返回到操作A1730(A1785)。当n为N/2(即,n不小于N/2)时,合成滤波设备11返回到操作A1790。
图18是根据本公开的实施例的用于解释选择滤波器系数的方法的流程图。上述实施例的任何多余的描述将在此省略。因此,可以参考上述说明而没有任何特殊的评论。同样,本公开不局限于图18的实施例并且选择滤波器系数的其他方法可以落入主张的主题的范围内。
首先,合成滤波设备11可以根据等式11初始化各变量(A1805)。例如,合成滤波设备11可以设置时隙索引l、采样索引n和滤波器级索引m的每个为0。
合成滤波设备11可以在第二阵列的位置[m+2Mn+MOD(l,M)]选择系数(A1810)。合成滤波设备11可以根据等式11选择系数。合成滤波设备11可以通过使用当前时隙和滤波器的级M进行模运算并且基于模运算的结果从第二阵列中选择滤波器系数。
如上所述,相同的滤波器系数在第二阵列中记录两次。以滤波器级M,记录相同的滤波系数的位置被彼此分离。因此,合成滤波设备11可以在具有大小M*N的第二阵列中选择(M*N/2)数量的滤波器系数,其中作为一组具有连续的M数量的滤波器系数的N数量的系数组的每个被连续记录两次。
合成滤波设备11可以通过将m增加1选择M数量的滤波器系数(A1840)。当m=M-1,(即,m不小于M-1)时,合成滤波设备11可确定是否n<N/2(A1825)。合成滤波设备11可以将n增加1直到n=N/2,(即,n不大于N/2),可以初始化m(例如,m=0),(A1830),并重复A1810的滤波系数的选择。如上所述,在通过预定的间隔划分从第一阵列提取的采样的M数量的组的每个中,分别对应于第n采样和第(n+1)采样的第一滤波器系数和第二滤波器系数在第二阵列中可以按2M彼此分离。
当n=N/2,(即,n不小于N/2)时,合成滤波设备11可确定时隙是否小于最大时隙L(A1835)。合成滤波设备11将时隙增加1直到时隙等于L(即,l大于L-1),初始化m和n(如n=0,m=0),并返回到操作A1810(A1845)。当时隙被增加1时,从第二阵列中选择的滤波器系数的位置可以从在先前的时隙中选择的滤波器系数的位置起移位1。当l=L,(即,l大于L-1)时,合成滤波设备11可确定选择滤波器系数的方法完成。
图19是根据本公开的实施例的用于解释合成时域采样的方法的流程图。上述实施例的任何多余的描述将在此省略。因此,可以参考上述说明而没有任何特殊的评论。同样,本公开不局限于图19的实施例并且合成时域采样的其他方法可以落入主张的主题的范围内。
虽然图18和19中示出的某些过程是相同或相似的,但是将从图18中选择的滤波器系数乘以采样并且累计乘法的结果的过程被进一步添加到图19的过程中。
合成滤波设备11可以根据等式12初始化各变量(A1905)(例如,l=0,n=0,m=0)。术语“temp”可以指一个变量,用于临时存储在时域采样的合成过程的中间输出的值,其也可以被初始化(例如temp=0)。
合成滤波设备11可以根据等式合成中间输出值,从而temp=temp+c'[m+2Mn+MOD(l,M)]*g[Nm/2+n](A1910)。合成滤波设备11可以通过将m增加1(A1920)来重复操作A1910直到m=M-1(即,m不小于M-1)(A1915)。当m=M-1,(即,m不大于M-1)时,合成滤波设备11可在在时域采样阵列y1[n]中存储temp的值并且可初始化temp(例如,temp=0)(A1950)。
合成滤波设备11可以通过将n增加1递归执行上述过程(A1930)。当n=N/2(即,n不小于N/2)时,合成滤波设备11可确定当前时隙是否是等于最大时隙L(A1935)。合成滤波设备11可以将时隙增加1(A1945),直到l=L(A1935),并且可以递归地执行上述过程,其中n和m再次初始化(例如,n=0,m=0)(A1945)。当l=L,(即,l大于L-1)时,合成滤波设备11可确定合成时域采样的方法完成。
图20是根据本公开的实施例的合成滤波装置11的框图。上述实施例的任何多余的描述将在此省略。因此,可以参考上述说明而没有任何特殊的评论。同样,图20所示的全部元件不都是本质的元件。可以在图20中进一步提供其他常见元件。常见元件的图示被省去以避免使得本公开的主题不清楚。
参考图20,合成滤波设备11可以包括解码器2001、调制单元2002、合成滤波器2003和存储器2004。解码器2001可解码经编码的比特流并且可以输出多个子带的信号。例如,比特流可以是多声道音频比特流。存储器2004可以存储第一阵列和第二阵列,其中子带信号的调制采样被记录在第一阵列中而滤波器系数被记录在第二阵列中。上述第三阵列也可被存储在存储器2004中。
调制单元(调制器)2002可以将子带信号的采样调制到基带的k数量的子带,并且可以在第一阵列中存储调制的采样。调制单元2002可以通过根据时隙的改变在第一阵列中循环采样记录位置来记录对应于当前时隙被调制的新采样。
合成滤波器2003可以包括多个滤波器,例如,其中的一些或全部可平行安置。为了消除在相邻子带的混叠,合成滤波器2003可提取来自第一阵列的采样,从第二阵列中选择对应于通过使用模运算提取的采样的滤波器系数,并通过使用提取的采样和选择的滤波器系数合成时域采样。合成滤波器2003可根据当前时隙是偶数或奇数确定从第一阵列提取的采样的位置。在这样做时,相同的滤波器系数可以在第二阵列中被记录两次,并且该相同的滤波器系数被记录的位置按照合成滤波器2003的级彼此分离。
合成滤波设备11可以通过各种设备实现,其可包括多媒体内容再现装置、音频系统、家庭影院系统、扩音器、移动终端、电视机、计算机、个人数字助理(PDA)、平板计算机、便携式媒体播放器(PMP)、智能电话机、或其他类型的电子设备,并且可以被包括作为上述设备的一部分。
在分析滤波过程中,时域采样被输入,从而混叠被移除,时域采样变换成K数量的子带,并且执行修正的离散余弦变换(MDCT变换)。根据实施例,存储器更新过程可通过在记录时域采样的阵列上执行模寻址而被省略。由于合成滤波过程和分析滤波过程具有相互对应的结构,本领域的普通技术人员可以理解根据上述合成滤波的实施例的分析滤波过程。
如上所述,根据一个或多个本公开的上述实施例的QMF滤波方法,QMF滤波可通过减小QMF滤波的复杂性而快速地执行,并且可以减少用于QMF滤波的功耗。
此外,本公开的上述实施例的各方面还可以通过在例如计算机可读介质的介质中/上的计算机可读代码/指令来实现,以控制至少一个处理元件来实现任何上述的实施例。该介质可以对应于允许计算机可读代码的存储和/或发送的任何介质/媒体。
计算机可读代码可以按各种方式在介质上记录/传送,其中介质的例子包括记录介质,诸如磁存储介质(例如,ROM,软盘,硬盘等)和光学记录介质(例如,CD-ROM或DVD)和传输介质,如因特网传输介质。因此,该介质可以这样包括或携带信号或信息的被定义并且可测量结构,诸如根据本公开的一个或多个实施例的携带比特流的设备。
例如,用于实现这里公开的实施例的计算机可读代码和/或程序指令可被记录在非暂时性计算机可读介质中。介质还可以包括(单独的或与程序指令组合的)数据文件,数据结构等等。非暂时性计算机可读媒体的例子包括磁介质,诸如硬盘、软盘和磁带;光学介质,诸如CD ROM盘和DVD;磁光介质,诸如光盘;以及专门配置来存储和执行程序指令的硬件设备,诸如只读存储器(ROM)、随机存取存储器(RAM)、快闪存储器、USB存储器等。程序指令的示例包括机器代码(诸如由编译器产生的)和文件两者,该文件包含可由计算机使用解释器执行的高级代码。所述程序指令可以由一个或多个处理器来执行。所描述的硬件设备可被配置为充当一个或多个软件模块,以便执行上述实施例的操作,或反之亦然。此外,非暂时性计算机可读存储介质可以在通过网络连接的计算机系统中分发并且计算机可读代码或程序指令可以分散的方式被存储和执行。此外,计算机可读存储介质也可以在至少一个专用集成电路(ASIC)或现场可编程门阵列(FPGA)中实施。
根据上述示例实施例的QMF滤波系统、合成滤波装置、以及使用QMF滤波系统和/或合成滤波装置或其部分实现的方法,可以使用一个或多个处理器。例如,处理设备可以使用一个或多个通用或专用计算机来实现,并且可以包括例如一个或多个处理器、控制器和算术逻辑单元、中央处理单元(CPU)、图形处理单元(GPU)、数字信号处理器(DSP)、微型计算机、现场可编程阵列、可编程逻辑单元、专用集成电路(ASIC)、微处理器、或能够以定义的方式响应并执行指令的任何其他设备。
如本文公开的,存储器可以存储一个或多个阵列和/或数据库。数据库可以包括数据的收集和可以被存储在例如存储设备(如存储器)的支持数据结构。例如,存储器可以实施为存储介质,诸如非易失性存储器设备,诸如只读存储器(ROM)、可编程只读存储器(PROM)、可擦除可编程只读存储器(EPROM)和闪速存储器,USB驱动器,易失性存储器设备,诸如随机存取存储器(RAM),硬盘,软盘,蓝光盘,或光介质,例如CD ROM盘和DVD,或它们的组合。然而,存储器的例子不限于上面的描述,并且存储器可以由其它各种设备和如将由本领域技术人员可以理解的结构来实现。
流程图的每个块可以表示单元、模块、段或代码部分,其包括用于实现指定的逻辑功能(多个功能)的一个或多个可执行指令。还应当指出的是,在一些替代实现方式中,块中示出的功能可能会不按顺序。例如,连续示出的两个块实际上可以基本上并发(同时)执行,或者这些块有时可以以相反的顺序执行,这取决于所涉及的功能。
应当理解,这里描述的示例实施例应该仅按描述性的意思考虑,而不是为了限制的目的来考虑。在每个实施例的特征或方面的描述通常应该被认为是可用于在其他实施例中的其他类似的特征或方面。
虽然已经参考附图描述本公开的一个或多个实施例,但是本领域的那些普通技术人员将理解,这里可以在形式和细节上进行各种改变而不脱离如随后的权利要求中定义的本公开的精神和范围。

Claims (12)

1.一种执行正交镜像滤波器(QMF)合成滤波的方法,该方法包括:
在包括调制的QMF子带采样的第一阵列中,在要从第一阵列丢弃的采样的位置处记录对应于当前时隙的新的采样;
从第一阵列中提取采样;
通过执行模运算确定对应于提取的采样的滤波器系数;以及
使用至少一个处理器,通过使用提取的采样和滤波器系数合成时域采样,
其中确定滤波器系数包括:
通过使用当前时隙和合成滤波器的级执行模运算,以及基于模运算的结果从第二阵列选择滤波器系数,
其中相同的滤波器系数在第二阵列中被记录两次,并且该相同的滤波器系数被记录的位置按照合成滤波器的级彼此分离。
2.根据权利要求1所述的方法,其中,在第一阵列中记录新的采样中,对应于当前时隙的新的采样的记录位置通过模运算在先前时隙中的采样记录位置来确定。
3.根据权利要求1所述的方法,其中,在第一阵列中记录新的采样中,新的采样被记录在从先前时隙中的采样记录位置偏移N数量的新的采样的位置处。
4.根据权利要求3所述的方法,其中,在第一阵列中记录新的采样中,当在先前时隙中的采样记录位置是第一阵列的开始时,从第一阵列的最后第N位置起记录新的采样。
5.根据权利要求1所述的方法,其中,在第一阵列中记录新的采样包括:
通过逆变换变换到基带的K数量的复子带采样产生对应于K数量的子带的N数量的实数采样;以及
在要丢弃采样的位置处记录N数量的实数采样。
6.根据权利要求1所述的方法,其中,在从第一阵列中提取采样中,根据当前时隙是偶数还是奇数确定要提取的采样的位置。
7.根据权利要求1所述的方法,其中,在从第一阵列中提取采样中,在当前时隙是偶数时从第一阵列提取的采样的位置与在当前时隙是奇数时从第一阵列提取的采样的位置相互排斥。
8.根据权利要求1所述的方法,其中,在从第一阵列中提取采样中,当第一阵列划分为多个部分时,对于自然数n,在当前时隙是偶数时在第(4n-3)和第4n部分中提取采样,以及在当前时隙是奇数时在第(4n-2)和第(4n-1)部分中提取采样。
9.根据权利要求1所述的方法,其中从第二阵列选择的滤波器系数的位置从自先前的时隙选择的滤波器系数的位置起移位1。
10.根据权利要求1所述的方法,其中,在合成时域采样中,通过累计提取的采样与确定的滤波器系数的相乘的结果来合成该时域采样。
11.根据权利要求1所述的方法,其中根据时隙的变化递归地执行所述记录、提取、确定和合成。
12.一种用于执行正交镜像滤波器(QMF)合成滤波的装置,该装置包括:
存储器,用于存储记录调制的子带信号的采样的第一阵列和记录滤波器系数的第二阵列;
调制器,用于将子带信号的采样调制到基带的K数量的子带,并且在第一阵列中存储调制的采样;和
合成滤波器,用于从第一阵列提取采样,使用模运算从第二阵列选择对应于提取的采样的滤波器系数,以及通过使用提取的采样和选择的滤波器系数合成时域采样,
其中调制器通过根据时隙的变化在第一阵列中循环采样记录位置来记录对应于当前时隙的新的采样,
其中所述合成滤波器进一步配置为:
通过使用当前时隙和合成滤波器的级执行模运算,以及基于模运算的结果从第二阵列选择滤波器系数,
其中相同的滤波器系数在第二阵列中被记录两次,并且该相同的滤波器系数被记录的位置按照合成滤波器的级彼此分离。
CN201480071503.6A 2013-10-28 2014-10-28 用于正交镜像滤波的方法和装置 Expired - Fee Related CN105849802B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
KR10-2013-0128645 2013-10-28
KR1020130128645A KR102244613B1 (ko) 2013-10-28 2013-10-28 Qmf 필터링 방법 및 이를 수행하는 장치
PCT/KR2014/010168 WO2015065002A1 (en) 2013-10-28 2014-10-28 Method and apparatus for quadrature mirror filtering cross-reference to related applications

Publications (2)

Publication Number Publication Date
CN105849802A CN105849802A (zh) 2016-08-10
CN105849802B true CN105849802B (zh) 2019-12-10

Family

ID=52996386

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201480071503.6A Expired - Fee Related CN105849802B (zh) 2013-10-28 2014-10-28 用于正交镜像滤波的方法和装置

Country Status (5)

Country Link
US (1) US9812140B2 (zh)
EP (1) EP3050053B1 (zh)
KR (1) KR102244613B1 (zh)
CN (1) CN105849802B (zh)
WO (1) WO2015065002A1 (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10180829B2 (en) * 2015-12-15 2019-01-15 Nxp Usa, Inc. System and method for modulo addressing vectorization with invariant code motion
CN109584277A (zh) * 2018-12-07 2019-04-05 上海应用技术大学 一种基于二次搜索的核相关滤波跟踪方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE465393B (sv) * 1990-01-16 1991-09-02 Ericsson Telefon Ab L M Adressprocessor foer en signalprocessor
US6199039B1 (en) * 1998-08-03 2001-03-06 National Science Council Synthesis subband filter in MPEG-II audio decoding
US8015368B2 (en) * 2007-04-20 2011-09-06 Siport, Inc. Processor extensions for accelerating spectral band replication

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5436940A (en) 1992-06-11 1995-07-25 Massachusetts Institute Of Technology Quadrature mirror filter banks and method
US5526464A (en) 1993-04-29 1996-06-11 Northern Telecom Limited Reducing search complexity for code-excited linear prediction (CELP) coding
JP3349421B2 (ja) * 1998-01-30 2002-11-25 三洋電機株式会社 音声データの圧縮/伸長装置及びデジタルフィルタ
US6446037B1 (en) 1999-08-09 2002-09-03 Dolby Laboratories Licensing Corporation Scalable coding method for high quality audio
US7333929B1 (en) 2001-09-13 2008-02-19 Chmounk Dmitri V Modular scalable compressed audio data stream
US8117248B2 (en) 2005-02-28 2012-02-14 Hitachi Global Storage Technologies Netherlands B.V. Digital filter instruction and filter implementing the filter instruction
US7676374B2 (en) * 2006-03-28 2010-03-09 Nokia Corporation Low complexity subband-domain filtering in the case of cascaded filter banks
EP4325723A3 (en) * 2006-10-25 2024-04-17 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Apparatus and method for generating time-domain audio samples
KR20080076691A (ko) * 2007-02-14 2008-08-20 엘지전자 주식회사 멀티채널 오디오신호 복호화방법 및 그 장치, 부호화방법및 그 장치
EP2144171B1 (en) * 2008-07-11 2018-05-16 Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. Audio encoder and decoder for encoding and decoding frames of a sampled audio signal
KR101599884B1 (ko) * 2009-08-18 2016-03-04 삼성전자주식회사 멀티 채널 오디오 디코딩 방법 및 장치
US9443534B2 (en) * 2010-04-14 2016-09-13 Huawei Technologies Co., Ltd. Bandwidth extension system and approach
US8886523B2 (en) * 2010-04-14 2014-11-11 Huawei Technologies Co., Ltd. Audio decoding based on audio class with control code for post-processing modes
PL2676268T3 (pl) 2011-02-14 2015-05-29 Fraunhofer Ges Forschung Urządzenie i sposób przetwarzania zdekodowanego sygnału audio w domenie widmowej

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SE465393B (sv) * 1990-01-16 1991-09-02 Ericsson Telefon Ab L M Adressprocessor foer en signalprocessor
US6199039B1 (en) * 1998-08-03 2001-03-06 National Science Council Synthesis subband filter in MPEG-II audio decoding
US8015368B2 (en) * 2007-04-20 2011-09-06 Siport, Inc. Processor extensions for accelerating spectral band replication

Also Published As

Publication number Publication date
CN105849802A (zh) 2016-08-10
EP3050053A1 (en) 2016-08-03
KR102244613B1 (ko) 2021-04-26
EP3050053B1 (en) 2019-12-11
WO2015065002A1 (en) 2015-05-07
US9812140B2 (en) 2017-11-07
US20150120306A1 (en) 2015-04-30
EP3050053A4 (en) 2017-03-15
KR20150048502A (ko) 2015-05-07

Similar Documents

Publication Publication Date Title
JP5265682B2 (ja) デジタルコンテンツの符号化および/または復号
CN101401305B (zh) 利用复调制滤波器组的滤波器、抽头产生器及滤波方法
DK2337224T3 (en) Filter unit and method for generating subband filter pulse response
KR100778349B1 (ko) 이산값의 시퀀스를 갖는 신호 처리 장치 및 방법
KR100776235B1 (ko) 변환 표시로의 변환 또는 변환 표시의 역변환을 위한 장치및 방법
JP2010537245A5 (zh)
TWI575962B (zh) 部份複數處理之重疊濾波器組中的低延遲實數至複數轉換
JP3814611B2 (ja) 時間離散オーディオサンプル値を処理する方法と装置
RU2323469C2 (ru) Устройство и способ для обработки, по меньшей мере, двух входных значений
JP2014524045A (ja) サンプルレートスケーラブル可逆オーディオコーディング
KR20210003507A (ko) 오디오 코딩을 위한 잔차 신호 처리 방법 및 오디오 처리 장치
US8433584B2 (en) Multi-channel audio decoding method and apparatus therefor
CN105849802B (zh) 用于正交镜像滤波的方法和装置
EP3329486A1 (en) Method and apparatus for generating from an hoa signal representation a mezzanine hoa signal representation
JP6089878B2 (ja) 直交変換装置、直交変換方法及び直交変換用コンピュータプログラムならびにオーディオ復号装置
JP2007178684A (ja) マルチチャンネルオーディオ復号装置
JP6907341B2 (ja) ステレオ信号処理方法及び装置
EP2784776B1 (en) Orthogonal transform apparatus, orthogonal transform method, orthogonal transform computer program, and audio decoding apparatus
US20210166701A1 (en) Device and method for encoding / decoding audio signal using filter bank
Kim Lossless wideband audio compression: Prediction and transform
TW202334938A (zh) 正交鏡像濾波器域中之沉浸式音訊及視訊服務空間重建濾波器庫
RU2643662C2 (ru) Устройство и способ для обработки звукового сигнала с использованием комбинирования в диапазоне перекрытия
Kim Lossless Wideband Audio Compression
JP2010537552A (ja) デジタル信号ソーティングに対するシステム、方法、および/または装置
CN113470616A (zh) 语音处理方法和装置以及声码器和声码器的训练方法

Legal Events

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

Granted publication date: 20191210

Termination date: 20211028

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