用于产生高质量声音信号的解码器和编码器

一般地说，本发明是关于诸如音乐信号等声音信号的高质量、低 比特速率(bit-rate)的数字变换编码。本发明通过对每个取样 声音段自适应地选择最佳变换(transform)、窗函数及变换块长 (block length)来实现变换代码中时间分辨力 (resolution)和频率分辨力之间的最佳权衡(trade off)。

本发明适用于所有离散正交变换。变换的正交性保证了利用正/ 逆变换能够实现严格的信号重建(Signal reconstruction)。 因此，本发明允许对正交变换自适应地选择块长而不丢失信息，就是 说，在没有量化误差(quantization errors)的情况下，能 够利用本发明的解码部分严格地恢复原始信号。

然而，在本发明的最佳实施例中却将自适应技术用于非正交变 换。自适应块长选择技术保持了变换的下述特点：(1)在没有系数 量化误差的情况下完全消除混淆(aliasing)；(2)临界采样 (Critical Sampling)。选择了一个变换相位项 (transform phase term)从而消除了时间域混淆(time doma in aliasing)。

在信号处理领域的人们当中存在相当大的兴趣去发现用最少量信 息来适当地表达给定信号的方法。通过减少所需信息，可以使信号在 具有较低带宽的信道(Communi cation Channels)上传送。 对于数字技术，最少信息需求与最少二进位需求是同义语。

有两个因素限制了位需求(bit requirement)的降低。第 一，带宽为W的信号可以由采样率不小于2·W的采样序列(Series of Samples)来精确表示。第二，采样信号可以是连续取值范围 中的任何一个值，对它的量化会使信号的代表值中引入不精确量，这 个不精确量与量化步长或分辨力成正比。这些不精确量叫作量化误 差，它们与所能得到的代表信号采样量化值的二进制位数成反比。

如果将编码技术应用于满带宽，则全部量化误差(它们表现为噪 声)将均匀地分扩散在这一带宽。可以应用于频谱中选定部分的技术 能够限制量化噪声的频率扩散(spread)。两种这类技术是“子 频段编码(Subband Coding)”和“变换编码(trans form Coding)”。通过使用这些技术，能够在量化噪声特别有害的那 些特定频段内，以较小步长量化那个频段，从而减小了量化误差。

子频段编码可以由多个数字带通滤波器来实现。变换编码可以由 若干个模仿一组数字带通滤波器的由时间域到频率域的离散变换中的 任何一种来实现。虽然离散变换与数字滤波器相比更易于实现，而且 需要较少的计算能力和硬件，但它们缺少设计灵活性，因为由一个变 换系数所代表的一个带通滤波器“频率仓”(frequency bin) 有一个均一的带宽。与此相反，一组带通滤波器能够设计成具有不同 的子频带宽。然而，变换系数能够组合在一起来确定其带宽为单个变 换系数带宽的倍数的“子频段”。这里用术语“子频段”来指整个信 号频带中的各选定部分，而不论是由子频段编码或是由变换编码来实 现的。由变换编码实现的子频段是由一组一个或多个相邻变换系数或 者说频率仓来确定的。一个变换编码频率仓的带宽决定于编码的采样 率和每个信号采样块中的采样数目(离散变换长度)。

子频段滤波器有两个特点，它们对于要达到高质量音乐信号编码 的给定主观(Subjective)性能的数字处理系统所需要的位速率 特别具有关键性。第一个特性是滤波器通带和阻带之间的区间(即过 渡带)的宽度。第二个特性是阻带的衰减水平。这里所用的滤波器的 “选择性”这一名词是指滤波器响应曲线在其过渡带的陡度 (Steepness)(过渡带滚降(rolloff)陡度)以及在阻带 的衰减水平(阻带的抑制深度)。

这两个滤波器特性对于编码的主观性能是极端重要的，这是因为 人耳表现出来的频率分析性能犹如具有可变中心频率的高度非对称调 谐(tuned)滤波器的频率分析性能。人耳调谐滤波器的频率分辨 力在整个音频范围内随频率而变化。在频率500赫兹以下，人耳能 分辨出频率靠近的几个信号，但是随着频率增大到可听 (audibiloty)极限，信号频率相差增大才能分辨开。这种听觉 滤波器的有效带宽称作临界带宽(Critical band)。

临界带宽的一个重要性质是在临界带宽内音质掩蔽效应最强烈地 表现出来一一在临界带宽内的一个主信号能压低该临界带宽内任何位 置的其他信号的可听能力。频率在临界带宽之外的信号不会这样强烈 地被掩蔽。一个主信号不仅可以掩蔽与掩蔽信号同时发生的其他信号 分量，而且可以掩蔽在掩蔽信号之前和之后发生的其他信号分量。在 一个临界带宽内的前掩蔽效应和后掩蔽效应的持续时间决定于掩蔽信 号的强度，但通常前掩蔽效应的持续时间比后掩蔽效应的持续时间短 得多。对此，一般可参阅声音工程手册(Audio Engineering Handbook)，K.B.Benson编，San Francisco(旧金 山)：McGraw-Hill出版，1988 ，第1.40-1.42和 4.8-4.10页。

音质掩蔽通过转换编码更容易实现，只要整个可听频谱范围内的 子频段带宽大约是在频谱中同一部分的人耳临界频带的一半。这是因 为人耳的临界频带具有可变的中心频率来适应声音刺激，而子频带和 转换编码一般具有固定的子频带中心频率。为了最佳地利用音质掩蔽 效应(PsyChoacoustic-masking effects)，由于一个 主信号的存在而造成的失真信号(distortion artifact)都 应该限制在含有该主信号的子频段内。如果子频段的带宽为临界频带 的一半左右或者小于一半，而且滤波器的选择性足够高，那么即使信 号频率靠近子频段通带带宽的边缘，也有可能对不希望要的失真信号 进行有效的掩蔽。如果子频段带宽大于一个临界频段的一半，则有可 能由于主信号使人耳的临界频带偏离编码子频带，从而使一些处在人 耳临界带宽之外的不希望的失真信号不能被掩蔽。这一效应在人耳临 界频带较窄的低频段是最有害的。

变换编码的性能和比特速率要求还依赖于其他因素，包括信号采 样块长、变换编码误差、以及消除假频(aliasing cancellation)。

                    块长(block length)

信号取样块长影响变换编码的时间分辨力和频率分辨力。对于较 长的块长，时间分辨力一般较差，因为量化误差会使变换编码/解码 系统“模糊”(Smear)了在整个信号采样块长度上的一个采样信号 的各频率分量。在由逆变换恢复的信号中的时间畸变对于持续时间短 的事件(event)或瞬时信号是最容易听得出来的。瞬变信号的时间 畸变将表现为瞬变前和瞬变后的振鸣(ringing)。

对于较短的块长其频率分辨力较差，因为此时离散变换频率仓 (bin)较宽而滤波器选择性较低(降低了过渡带滚降率并减低了 阻带抑制水平)。这种滤波器性能的降低会造成在邻近的频率仓内出 现不希望的变换系数来响应所需信号。这些不需要的分量是所谓“边 瓣泄漏(Sidelobe Leakage)”效应的结果。

这样，短块长可以造成标称(nominal)滤波器带宽在较低频 段或全部频段上起过人耳临界带宽。即使标称子频段带宽比人耳临界 带宽要窄，其滤波器性能的降低却表现为一个宽过渡带和/或差的阻 带抑制水平，可以在人耳临界带宽之外造成显著的信号分量。另一方 面，长块长可能造成差的时间分辨力，使变换编码造成出现在人耳时 间音质掩蔽间隔之外的信号畸变分量。

                    变换代码误差

离散变换不能产生绝对精确的一组频率系数，因为这些离散变换 只作用于有限长的一段信号。严格地讲，离散变换产生的是时间域信 号的时间-频率表示，並不是它的真正的频率域表示，后者需要无 限变换长度。然而，这里为便于讨论，将把离散变换的输出称作频率 域表示。从效果上看，离散变换假定采样信号只含有其周期为有限采 样段长度的约数(Submultiple)的那些频率分量。这等效于假 定有限长信号是周期性的。当然，这一假定通常并不成立。假定的周 期性造成有限时间段边缘的不连续性，它使变换产生虚假的高频分 量。

减小这种效应的一种技术是在变换之前通过对信号采样加权使其 时间段边缘附近的采样为零或接近于零，借以减小其边缘不连续性。 在时间段中央的采样一般为不变化地通过，即加权因子为1。这一加 权函数称作“分析窗”。窗的形状直接影响滤波器的选择性，降低过 渡带滚降率和阻带抑制水平。某些窗的形状所引起的这种滤波器性能 降低要大大小于其他窗形状。

这里使用的术语“分析窗”只是指实行正变换之前使用的窗函 数。如下面将要讨论的，在本发明中使用的分析窗可以由合成窗设计 考虑(Synthesis window design consid erations) 来约束。所以，作为本领域技术中一般使用的术语“分析窗”的设计 和性能特性可能与本发明中应用的这种分析窗不同。

分析窗是一个时间域函数。如果不提供窗效应的补偿，则恢复的 或“合成的(Synthesized)”信号将会依其分析窗的形状而产生 相应的畸变。称作“重迭加(Overlap-add)”的一种补偿方法 是本领域所熟知的。这种方法要求编码器去变换输入信号采样的重迭 块。通过小心地设计分析窗，使两个相邻窗在其重迭部分相加为单位 1，于是窗效应被完全补偿。

当采用某些类型的变换，例如离散付立叶变换(DFT)时，这 项补偿技术增加了表示该信号所需二进位的数目，因为在重迭时间段 的信号部分必须变换和传送两次，即在两个重迭信号采样块各有一 次。这种变换不是临界采样的。术语“临界采样的(Critically Sampled)”是指该变换在其时间区间上给出的频率系数数目与它 接收的输入信号采样数目相同。因此，对于非临界采样变换，希望所 设计的窗的重迭区间尽可能小。

某些变换(包括本发明最佳实施例中使用的那一种)要求逆变换 的合成输出加窗。合成窗(Synthesis window)用于对每个合 成信号块整形。所以，合成信号将是受过分析窗和合成窗调制过的。 这种两步调制在数学上类似于用这样一个窗来调制原始信号一次，该 窗的形状等于分析窗和合成窗的采样点逐点乘积。所以，为了用重迭 相加来补偿窗畸变，两个窗的设计必须使得两个窗的乘积在重迭相加 区间之和为单位1。

尽管没有唯一的标准可用来评价一个窗的好坏，但是，如果与这 个窗合用的滤波器的选择性被认为是“好”的，通常就认为这个窗是 “好”的，所以，一对设计得好的分析/合成窗能减少边瓣泄漏。

减少边瓣泄漏是重要的，因为泄漏会使变换造成的频谱系数错误 地表示成滤波器通带外的信号分量频率。这种错误表示 (misrepresentation)是一种畸变，称作“混淆 (aliasing)。”

                       消除假频

奈奎斯特(Nyquist)定理指出，当采样点之间的间隔不大于 信号最高频分量的周期的一半时，该信号能由离散采样精确地恢复。 当采样率低于这个奈奎斯特采样率时，较高频分量被错误地表示成较 低频分量。这种较低频分量是真正分量的“假频(alias)”分 量。

子频段滤波器和有限长数值变换都不是完美的带通滤波器。通带 和阻带之间的过渡不是无限陡的，而且在阻带中信号的衰减不是无限 大。结果，即使将一个带通滤波后的信号按照由通带截止频率给出 的奈奎斯特采样率来采样，那么能够通过该滤波器的超过滤波器标称 频率的信号(当然被衰减了)将不能被可靠地表示出来。

有可能设计出分析滤波器和合成滤波器从而使假频畸变被逆变换 自动消掉。二次镜象滤波器在时间域具有这种特性。某些变换编码技 术，包括本发明的最佳实施例中使用的那种，也能消掉假频畸变。

随着采样块长变短，压低变换编码中假频畸变的可听后果 (audible Consequences)就变得更加困难。如前面解释过 的，较短的采样块会降低滤波器性能：通带宽度增大，通带到阻带的 过渡陡度降低，以及阻带抑制能力减弱。结果，假定信号以滤波器标 称截止频率作为奈奎斯特采样率进行采样，那么在截止频率以上的频 率分量以较小衰减通过滤波器从而使假频混淆变得更加显著。

即使在原理上正变换和逆变换可以消除假频畸变，但量化误差将 使逆变换不能完全消除假频畸变。残余假频畸变将是可听见的，除非 这种畸变在音质上(Psychoacoüstically)被掩蔽掉。然而， 以短采样块来变换时，一些变换频率仓(bin)的通带可能大于听 觉临界频带，特别是在低频段，其人耳临界频带有最大的分辨力。结 果，假频畸变可能未被掩蔽。减小这种畸变的一种途径是在有问题的 子频段中增大量化精度，但这会增加所需的比特速率。

本发明的一个目的是采用一种编码/解码技术和方法对宽带声音 信息(特别是音乐)进行数字处理，它能以96千位/秒(kbs) 的低编码速率提供出高主观质量的声音。

本发明的目的之二是所提供的这种编码/解码技术和方法适于宽 频带声音信号的高质量传送、存贮和再生，这里的再生质量适合于例 如广播声音链路。

本发明的目的之三是使再生质量的主观感觉如同从密纹唱片得到 的一样好。

本发明的目的之四是所提供的这种编码/解码装置和方法所采用 的数字处理系统只需要小量空间来存贮编码信号。

本发明的目的之五是提供来选择最佳变换块长度。

本发明的目的之六是用于对每个采样声音信号段根据信号的特性 (例如是否是瞬变信号)从一组窗函数中选择出最佳的分析窗函数。

本发明的目的之七是用于对每个采样声音信号段根据其信号特性 (例如是否是瞬变信号)来从一组离散变换函数中选出最佳离散变 换。

本发明的目的之八是对于用低比特速率离散变换来处理宽频带声 音信号特别是音乐信号时产生的畸变提供时间上的音质掩蔽。

关于本发明的上述目的及其他目的的进一步细节将在本发明的文 件中给出，特别是在下文中的“本发明详细描述”一节中给出。

根据本发明的原理，一个编码器提供宽带声音信息的数字编码功 能。宽带声音信号被采样、量化和以N个采样点为一组组成时间域信 号采样块。一个信号分析器(analyzer)对当前信号采样块进行 数值分析以确定适当的变换类型、块长和窗类型，以使编码保真度 (fidelity)最佳。在一个具体实施例中，有一个瞬变信号检测 器来确定在当前信号采样块中是否存在任何瞬变信号需要更短的块长 度以避免可听见的瞬时畸变。具有这类瞬变信号的采样块被再分割成 最佳长度采样子块，以保持足够的变换频率选择性，并保证在音质上 掩蔽掉由瞬变信号引起的编码畸变。每个采样块，不论是正常长度的 还是减小了长度的，都被一个根据瞬变信号检测器输出而选定的分析 窗函数来调制。然后，根据瞬变信号检测器输出所选定的离散正变 换，响应经分析窗加权的时间域采样块，产生出频率域频谱分量。确 定编码器所用信号采样块长度、分析窗函数和正变换类型所需的信息 作为编码变换系数的伴随信息被传送给解码器。

也是根据本发明的原理，一个解码器提供了由本发明的编码器编 码的数字编码宽带声音信号的高质量再生信号。将确定所选信号采样 块长度、分析窗函数和正变换类型的伴随信息从编码信号中提取出 来。将这一信息用于建立逆变换长度，并报告合成窗函数和离散逆变 换类型的选择。通过一个具有与产生频率域频谱分量的编码器中的特 性相反特性的离散变换，响应其频率域频谱分量而产生时间域信号采 样块。在使用离散变换的本发明实施例中需要各采样块由合成窗调 制。如果使用的话，该解码器合成窗和编码器分析窗的形状应使两个 相邻重迭采样块的组合窗函数之和为单位1。将相邻采样块重迭相加 以消除窗调制的加权效应和恢复时间域信号的数字再生信号，然后再 把它转换成高质量模拟输出。

再根据本发明的技术，一个编码/解码系统提供了宽频带声音信 息的数字编码和高质量再生。在该系统的编码器部分，宽带声音信号 被采样、量化和按N个采样点长度将信号分组构成时间域信号采样 块。一个信号分析器对当前信号采样块进行数值分析以确定适当的变 换类型、块长和窗类型，以使编码保真度最佳。在一个具体实施例 中，有一个瞬变信号检测器来确定在当前信号采样块中是否存在任何 瞬变信号需要更短的块长度以避免可听见的瞬时畸变。将具有这类瞬 变信号的采样块再分割成最佳长度采样子块，以保持足够的变换频率 选择性，并保证从音质上掩蔽掉由瞬变信号引起的编码畸变。每个采 样块，不论是正常长度还是减小了长度的，都被一个根据瞬变信号检 测器输出而选定的分析窗函数来调制。然后，根据瞬变信号检测器输 出而选定的离散正变换，响应经分析窗加权的时间域采样块，从而产 生出频率域频谱分量。将确定编码器所用信号采样块长度、分析窗函 数和正变换类型所需的信息作为编码变换系数的伴随信息传送给解码 器。系统的解码器部分经过存贮装置或传输路径接收编码器的数字输 出。解码器部分从编码信号中提取确定所选定的信号采样块长度、分 析窗函数及正变换类型的伴随信息。将这一信息用于建立逆变换长 度，并报告合成窗函数和离散逆变换类型的选择。通过一个与产生频 率域频谱分量的编码器中的特性相反的离散变换，响应其频率域频谱 分量，产生时间域信号采样块。在使用离散变换的本发明实施例中需 要各采样块由合成窗调制。如果使用的话，该解码器合成窗和编码器 分析窗的形状应使两个相邻重迭采样块的组合窗函数之和为单位1。 将相邻采样块被重迭相加以消除窗调制的加权效应和恢复时间域信号 的数字再生信号，然后再把它转换成高质量模拟输出。

根据本发明的技术，在一个具体实施例中有一个瞬变信号检测 器，当需要以较短的信号长度来保证在音质上掩蔽由大瞬变信号造成 的畸变信号时，该检测器便相应地选择比本发明通常使用的最大块长 要短的信号采样块长度。让数字化信号采样通过一个高通滤波器并被 分成子块。将每个子块中的峰值幅度与前一个子块中的峰值幅度相比 较。如果相邻峰值振幅的变化没有超过预先给定的阈值，则编码器使 用最大信号采样块长度。如果相邻子块的振幅变化超过了预定的阈 值，则选择较短的块长度。

对于每个信号采样块，选定一个分析窗函数，它适合于信号采样 块的信号成分和由瞬变信号检测器选定的块长度。如前面讨论过的， 这个分析窗函数用于调制信号采样块内的采样，以改善离散变换的总 体性能。在本发明的编码器的一个实施例中，分析窗是由同一类型而 不同长度的多个窗函数中选出的。

在本发明中可以采用其形状随着块长度或瞬时特性的函数而改变 的窗函数，而不偏离本发明的范围和目的。例如，本发明的一个实施 例可以选择这样一个窗，它对于含有瞬变信号的信号采样块提供较高 的阻带抑制能力，其代价是过通带滚降较浅(Shallower)。然 而，窗函数的选择不能违反约束条件，例如重送相加特性的约束。

在本发明的编码器的一个实施例中，离散正变换响应经分析窗加 权的时间域来样块产生频率域频谱分量。离散变换的长度置成等于由 瞬变信号检测器选定的信号采样块长度。

任何正交的时域到频域变换都能使用。但本发明的一个最佳实施 例使用了一个非正交离散变换，它等效于一个修正的离散余弦变换 (DCT)和一个修正的离散正弦变换(DST)的交替使用。在另 一个实施例中，由单一的修正的离散余弦变换(DCT)来实现其非 正交离散变换。

精通本领域技术的人应该理解，可以根据选定的块长或瞬变信号 特性选定相应的不同正交离散变换。可以从多个离散变换中适当地选 择一个离散正交变换，而不偏离本发明的范围或精神。

例如，本发明的一个实施例能够正常使用一个通常的DCT(离 散余弦变换)并选择一个通常的DST(离散正弦变换)用于宽有瞬 时脉冲的信号采样块。在当今技术中已知，低频分量的编码精度对于 大多数音乐信号实现给定的主观质量而言是很足够的。通常的DCT 对低频信号具有选越的编码性能。另一方面，较高频率分量的编码精 度对于瞬变信号实现给定的主观质量水平是很重要的。通常的DST 对于较高频率信号具有优越的编码性能。所以，本发明的一个实施例 中可以适应地选择最适于每个信号段编码的变换。

在编码器和解码器的一个最佳实施例中，所用采样率为44.1 千赫兹。尽管该采样率不是临界采样率，44.1千赫兹是一个适当 的采样率，而且它是方便的，因为它也是密纹唱片采用的采样率。另 一个实施例采用48千赫兹采样率。在采用44．1千赫兹采样率的 最佳实施例中，标称频率响应延伸到15千赫兹，时间域采样块的最 大长度为1024个采样点。可以适当地选择较短的长度：512、 256、或128个采样点。在本发明的一个最佳实施例中，在适于 专业广播使用的主观质量水平上的音乐编码可以用低达96Kbs (包括错误校正代码等附加信息)的几种位速率来实现。也可以使用 产生不同水平信号质量的其他速率而不偏离本发明的基本精神。

在下面的“本发明详细描述”一节将结合附图更详细地说明本发 明的各个特点及其最佳实施例。

图1a和1b是说明本发明基本功能结构的功能块图。

图2a到图2e中的方块图给出本发明一个实施例的硬件构成。

图3a和3b中的方块图更详细地给出本发明的一个双通道实施例 中处理机的串行通信接口。

图4a到图4e为假设的图形表示，用以说明一个时间域信号分组 成一系列重迭的、加窗的时间域信号采样块。

图5a到5d为假设的图形表示，用以说明由E-TDAC变换产 生的时间域假频畸变。

图6a到图6g为假设的图形表示，用以说明借助E-TDAC变 换信号合成过程中的重迭相加来抵消时间域假频信号。

图7q是假设的图形表示，用以说明由一对信号采样块构成单道系 统中的框架。

图7b是假设的图形表示，用以说明由一对信号采样块构成双通 道系统中的框架。

图8a是带有一个瞬时脉冲的模拟音乐信号段的图形表示。

图8b是合成音乐信号段的图形表示，它是用一个固定的信号采 样块长512采样点，通过编码/解码系统从一个信号采样块恢复而 来的。

图8c是合成音乐信号段的图形表示，它是由本发明的一个实施 例采用固定帧排列(Eixed-Frame Alignment)技术以最大 信号采样块长512采样点恢复得来的。

图9的示意图表示一个递归滤波器，它实现了本发明编码器部分 的瞬时信号检测器中所使用的高通滤波器。

图10是一个假设的图形表示，说明瞬变信号检测器如何将一个 信号采样块的一半分成一组更短长度的子块。

图11是一个假设的图形表示，说明瞬变信号检测器如何在一个 层次结构(hierarchical)子帧内识别每一个子块中的举振辐 值。

图12是一个二进制树的假设的图形表示，说明由瞬变信号检 测器的比较器部分构成的节点。

图13是一个流程图，给出本发明一个实施例中使用的瞬变信号 检测器的比较器部分的逻辑。

图14是一个假设的图形表示，说明由瞬变信号检测器的比较器 部分对二进制树的修正过程。

图15是一个假设的图形表示，说明由瞬变信号检测器根据所示 二进制树选定的子块长度。

图16是一个音乐波形的图形表示，它是按本发明的瞬变信号检 测器部分所选定的子块长度区间迭加而成的。

图17是时间域信号采样块的假设图形表示。

图18是时间域信号采样块的又一个假设图形表示，说明当假定 块内信号为周期性时由离散变换引起的在采样块边缘的不连续性。

图19a是一张功能块图，给出函数X(t)被函数W(t)调 制，结果生成函数Y(t)。

图19b到19d是说明时间域信号采样块被一个分析窗调制情 况的进一步假设图形表示。

图20是适用于本发明的一个最佳实施例的一族分析-合成窗的 图形表示。

图21a至21c是一系列信号采样块的假设图形表示，用以说 明固定帧排列(Fixed-Framc Alignment)技术。

图22a用图形表示出使用完全对称的只作分析的窗的一组滤波 器的过渡带滚降和阻带抑制，并与只用这种窗的前半部的一组滤波器 的过渡带滚降和阻带抑制相比较。

图22b用图形表示出使用完全对称的只作分析的窗的一组滤波 器的过渡带滚降和阻带抑制与使用非对称的只作分析的窗的一组滤波 器的过渡带滚降和阻带抑制的比较。

图23a是一个假设图形表示，说明由E-TDAC变换采用通 常的相应项造成的时间域假频分量信号的时间逆转(reversal) 区间。

图23b是一个假设图形表示，说明由E-TDAC变换采用固 定框架排列技术所要求的相位项造成的时间域假频分量信号的时间逆 转区间。

图23c是一个假设图形表示，说明采用为消除一个信号采样子 块中的时间域假频信号所需要的相位项时，由E-TDAC造成的时 间域假频分量信号的时间逆转区间。

图24是一种桥变换(bridge transform)的假设图形表 示，说明时间域假频分量信号的时间逆转区间。

图25是一系列信号采样块的假设图形表示，说明扩展的固定帧 排列技术。

图26a至26f是流程图，说明本发明的扩展固定帧排列技术 实施例所需要的帧控制逻辑。    

图27是一系列可变长度信号采样块的假设图形表示，用以说明 动态帧排列(Dynamic-Framc Alignmcnt)技术。

图28的假设图形表示显示出相邻加窗块的重迭-相加性质。

图29a至29e为假设图形表示，说明一个时间域信号被分成 一系列重迭的和加权的时间域信号采样块，特别是为O-TDAC变 换使用。

图30a至30d中的假设图形表示说明由O-TDAC变换造 成的时间域假频畸变。

图31a至31g中的假设图形表示说明在O-TDAC变换的 信号合成过程中通过重迭-相加得到的时间域假频信号抵消。

图32的假设图形表示说明为消除时间域假频信号由固定帧排列 所要求的O-TDAC变换相位项。

图33a和33b是E-TDAC和O-TDAC变换的变换系 数带宽的假设图形表示。

表Ⅰ给出用于瞬变信号检测器第一部分的aKHj高通滤波器的 系数。

表Ⅱ给出用于瞬变信号检测器第一部分的4KHz高通滤波器的 系数。

表Ⅲ给出瞬变信号检测器第4部分中的增高/衰减(attack／ decay)阈值，用于确定构成瞬时信号状态所需的信号振幅变化 量。

1．本发明的硬件实现

图1a和1b给出本发明的基本功能结构。图1a中所示本发明 的编码器部分包括时间域信号输入102、低通滤波器104、信号 采样器和量化器106、信号采样缓存108、瞬变信号检测器110 (它选择信号采样块长度)、分析窗调制器112(它用分析窗函数 调制数字化时间域信号块)、数字滤波器组116(它把采样和量化 的信号变换成频率系数)、帧(frame)控制器114(它根据选 定的信号采样块长度来控制分析窗调制器和滤波器组)、量化器118 (它根据所要求的主观信号质量和音质效应来对变换系数编码)、以 及格式器(formatter)120(它把编码的频率系数和选定的信 号采样块长度组合成二进制位流，供传输或存贮之用)。图1a中显 示出一个传输路径122，然而应该理解，其编码信号可以被存贮起 来供其后使用。

图1b所示本发明的解码器部分包括：编码二进位流信号输入 132、格式解析器(deformatter)134(它从组合二进位 流中提取每个编码频率系数和选定的信号采样块长度)、线性化器 (Linearixer)136(它把每个编码系数转换成线性取值的变 换系数)、帧控制器140(它根据所提取的信号采样块长度来控制 反滤波器组和合成窗调制器)、逆字滤波器组138(它将变换系数 变换成时间域信号块)、合成窗调制器142(它用选定的合成窗函 数来调制每个合成的时间域信号块)、信号块重迭相加器144(它 恢复时间域信号的数字表示)、模拟信号发生器145、低通滤波器 148、以及模拟信号输出150。

A．处理硬件

图2a-2c和图3a-3b绘出本发明的基本硬件结构。经验 研究表明：常规的整数变换计算必须以至少20个有效二进位的精度 来完成才能达到所宣称的性能指标。

本发明的一个单通道方案(采用44.1KHz采样率或48 KHq采样率)的最佳实施例的实际实施中利用了一个16位模-数 转换器(ADC)，其周期时间不超过20微秒，用于量化输入的时 间域信号。每个16位数字化采样用于构成24位字的16个高有效 位，该24位字将用于其后的计算。以20.5MHz无等待形态运 行的Motorola DSP56001 24位数字信号处理器 (DSP)用于完成所需要的计算和控制编码与解码处理过程。静态 随机存取存贮器(RAM)为DSP提供程序和数据存贮器。周期时 间不大于20微秒的16位数-模特换器(DAC)用于从解码的数 字信号产生模拟信号。

图2a所示编码器硬件结构由下列部件构成：模拟信号输入 200、低通滤波器(LPF)200A、ADC201、DSP 202、静态RAM203、可抹挣可编程只读存贮器(EPROM) 204、及编码串行信号输出206。LPF200A保证输入信号 为有限带宽。ADC201将输入信号数字化(采样和量化)，变成 串行16位字流。DSP202按收并暂存数字化采样串行流。检测 是否存在任何瞬变信号、选择信号来样块长度、将采样分组构成采样 块、用分析窗函数调制采样块、完成将采样块变换到频率域所需的计 算、对变换系数编码、将代码字和信号采样块长度格式组合到数据流 中，并通过串行数据通道206将编码信号传送出去。DSP的程序 和数据工作区存贮于3个24KB(千字节)静态RAM组203 中，它们组织成了3组，各含有8，192个24位字。DSP要求 快速访问时间的程序存贮器，这可以在RAM中实现，比在可编程 ROM中实现要更廉价。因此，EPRM204以压缩格式存贮程 序和静态数据，当编码器第一次加电时DSP将压缩格式的数据分解 成可使用形式送入RAM203。

图2b和2c给出了2个DSP接口的更详细结构。图2b给出 DSP202、ADC201及串行数据通道206的数据通信接 口。时间发生器202A为编码器产生接收计时、字同步、和传送计 时信号。线SCO计时控制数字化输入信号采样的串行二进位流从 ADC201经母线SDR进入DSP202。线SCI向ADC和 DSP提供字同步信号，用于标志每个16位字的起始。线SCK计 时控制编码信号串行位流沿着线STD从DSP进入串行数据通道 206 。

图2c给出存贮器选址接口。Motorola Dsp56001的存贮 器被分成四段：程序ROM、程序RAM、X数据、及Y数据。当线 PS为低时只有程序RAM203或EPROM204被选通，但它 们被映射(mapped)到各自的地址空间。转换器205C允许 DSP202根据地址线A15的状态来选择BAM或EPROM。 当DSP置A15为高时，转换器205C置RAM203和 EPROM204的片选通(Chip-Select)(CS)线为 低。当DSP202置A15为低时，转换器205C置RAM 203和EPROM204的CS线为高。当线CS为高而线PS为 低时，选通程序RAM203。

当DSP使线DS为低和线XY为高时，选通RAM203的X 数据库。当DSP使线Xy为低和线CS为高时，RAM203的y 数据库被选通。

图2d给出解码器硬件结构，它包括：编码的串行信号输入路经 207、DSP208、静态RAM209、EPROM210、 DAC212、LPF213A、以及模拟信号输出213。DSP 208接受和暂存编码信号、将信号格式解析成编码变换系数和信号 采样块长、完成将系数变换到时间域所需的计算。将系数分组构成时 间域块、用合成窗函数调制这些块(如果在编码器/解码器中使用的 离散变换不需要使用合成窗，则不需进行合成窗调制)、将这些块重 迭相加构成数字采样的时间域序列、以及将串行位流形式的数字采样 传送给DAC212。DSP的程序和数据工作区存贮于3个24个 KB静态RAM209组中，它们被组织成3组，各有8，192个 24位字。EPROM210以压缩格式存贮程序和静态数据中，当 解码器第一次加电时DSP将这些压缩格式数据分解成可用形式送入 RAM209。DAC212产生一个模拟信号，对应于来自DSP 的串行数据流。LPF213A保证信号输出213不含有由于编 码/解码过程造成的任何虚假的高频分量。

图2e给出DSP208、串行信号输入通道207及DAC 212之间的串行通信接口。时间信号发出器208A使用锁相环电 路从编码串行位输入信号中提取出定时基准信号，为解码器产生接收 定时信号、字同步、及传送定时信号。线SCO定时控制串行位信号 沿线SRD进入CSP208。线SCK计算控制编码数字信号采样 的串行位流沿线STD从DSP208进入DAC212。线SC2 向DAC和DSP提供字同步信号，它标志每个16位字的起始。 DSP208和内存贮器地址总线之间的接口的实现方式与前文中对 编码器的同类接口的实现方式相同。参见图2C。    

基本硬件结构可以被修正。例如，一个以27MHz无等待状态 运行的，Motorola DSP56001如果采用简单量化技术，便 能够实现双通道编码器或解码器。

编码器的双通道实施例需要有LPF200A和和200B、 ADC201A和201B，按图3a所示方式联接。在DSP和 ADC部件之间是接口的运行方式与上文描述的单通道编码器的接口运 行方式类似。定时信号发生器202A以字同步信号速率一串的速率 向DSP的线SC2提供一个附加信号，用以控制多路转换器202 B和向DSP指出两个ADC中当前是哪一个在发送数字化数据。

一个解码器的双通道实施例需要DAC212A和212B、 LPF213A和213B，按图3b所示方式联接。DSP和DAC 部件之间的接口运行方式与上文描述的单通道解码器接口的运行方式 类似。定时信号发生器208A以字同步信号速率的一串的速率向 DSP的线SC1提供一个附加信号，用于控制信号分离器 (dcmultiplexer)208B并向DSP指出两个DAC中哪一 个当前正在接收数字数据。

再有，可以使用专门化的硬件来实现某些功能，例如窗调制或快 速向立叶变换(FFT)。整个编码器/解码器可以在一个专门设计 的集成电路中实现。对于精通本门技术的人而言，显然有许多其他可 能的实现方法。

B．输入信号采样和缓存

在本发明的一个最佳实施例中，图1a所示的信号采样器和量化 器106是一个模一数转换器，它把输入信号量化成16位，然后在 右边补上8个零，构成一个24位二进制整数形式。所有其后的变换 计算都是以24位二进制整数的算术运算来实现的。模拟输入信号被 图1a中框104所示低通滤波器(LPF)限制其宽带最高为15 KHz(对于20KHz带宽的编码器，带宽限制在20KHz)。

具有至少是密纹唱片(CD)质量的音乐信号，除了具有其他特 性外，其带宽超过15KHz。根据奈奎斯特定理，已知15KHz 带宽的信号采样率必须不低于30KHz 。对本发明的一个实施例选 用了44。1KHz的采样率，因为这是在CD的应用中使用的速 率，而且这一选择简化了在这类应用中使用本发明所必须的装置， (这一采样率也支持本发明的另一种20KHz带宽的实施例)。

也可以利用其他采样率，例如48KHz是通常许多专业性声音 系统应用中使用的速率。如果选用另一个采样率，那么相邻变换系数 之间的频率间隔将要改变，代表希望的信号带宽所需要的系数个数也 将改变。对于精通本门技术的人而言，采样率变化对实现本发明所造 成的全部影响是显而易见的。

Ⅱ．本发明的最佳实施

A．分析滤波器族(bank)-正变换

实现本发明的细节受所选择的离散变换的影响，该离散变换用于 实现图1a中框116所代表的数字滤波器族的功能。几种变换技术 中的任何一个都可以用来实现这滤波器族。在本发明的一个最佳实施 例中使用的变换技术首先见Princen和Bradley的“基于在时 域抵消假频信号的分析/合成滤波器族设计”IEEE Trans. on Acoust，Speech，Signal Proc(声学、话言、 号处理)，Vol(卷)ASSP-34，1986，1153- 1161页。这一技术是一个偶迭加临界采样单边分析合成系统 (evenly-stacked cnticallq sampled single- sideband analqsis-synthesis system)的时间域 等效技术。在此地将这一变换称作“偶送加时间域抵消假频” (Evenlq-Stacked Time-Domain Aliasing Cancellation)(E-TDAC)技术。TDAC变换的另一 种形式可以用于本发明的另一实施例。该技术在Princen、 Johnson和Bradleq的下列文中描述过：“基于时间域消除 假频信号的子频带/变换编码用滤波器族的设计”，ICASSP 1987 Conf.Proc.，1987年5月，2161-64页 这另一种变换是“迭加临界采样单边带分析合成系统”的时间域等 效技术。在此地将这一变换称作“迭加时间域抵消假频”技术(O -TDAC)。

在下面充分描述了E-TDAC实施例之后，将讨论使用O- TDC变换和用快速付立叶变换(FFT)来实现离散付立叶变换 (DFT)的本发明实施例。

E-TDAC使用的变换函数等效于交替使用修正的离散余弦变 换(DCT)和修正的离散正弦变换(DST)。方程1所示DCT 和方程2所示DST是： 这里K=频率系数序号，

n=输入信号采样序号，

N=信号采样块长度，

m=E-TDAC的相位项，

x(n)=输入信号x(t)在采样点n的量化值，

C(k)=DCT第K个系数，

S(k)=DST第K个系数。

E-TDAC变换对于每个信号采样块交替产生两组频谱系数或 变换块之一。这些变换块的形式是 这里i=信号采样块序号，

C(K)=DCT系数(见方程1)，

S(K)=DST系数(见方程2)。

Princen指出，利用适当的相位分量m(见方程6)和精心设 计的一对分析合成窗，E-TDAC技术能从固定长余弦和正弦变 换块的交替序列中精确地恢复输入信号。这一序列的形式是

{C(k)}0，{S(k)}1，{C(k)}2，{S(k)}3，…(5) 这里每个变换块代表一个时间域信号采样块。这一处理过程见图 4a-4e、5a-5d及6a-6g。

参考图4a，可以看到量化输入信号X(t)被分组构成块。一 组块经过窗函数Wc(示于图4b)调制，产生信号Xc(t)(示于 图4d)。信号Xc(t)输入到DCT。另一组采样输入信号X(t) 的块以二分之一块长与第一组重迭，被窗函数Ws(示于图4c)加 窗(这个窗函数与Wc完全相同，但在时间上移动了二分之一块长)， 产生出信号Xs(t)(图4e所示)，然后被送到DST。

如果只是交替使用DCT和DST变换块，会使变换块的被去挣 的那一半中所包含的信息丢失。这一损失产生了时间域假频分量，但 这一畸变可以通过在方程1和2中选择适当的相位项m、应用正变换 去重迭时间域信号块以及通过重迭和相加由逆变换恢复的相相邻的时 间域信号采样块的办法来消除。

方程1和2中的相位项m控制时间域假频畸变的相移。图5a- 5b和6a-6g以图说明了这一畸变。从逆DCT恢复的信号Yc (t)示于图5a。图5b说明了被恢复的信号由两个分量组成：加 窗的原始信号(实线)及时间域假频畸变(虚线)。图5c和5d给 出由逆DST恢复的信号Ys(t)的类似信息。为了消除这一假频 畸变並精确地恢复原始时间域信号，E-TOAC要求其假频信号如 下所述：对于DCT，时间域假频分量的构成是：采样信号的前一半 在时间上以来样块四分之一点为参考点逆转，采样信号的后一半在时 间上以采样块四分之三点为参考点逆转。对于DST，假频分量与 DCT的相似，只是其振恢的符号相反。见图5b和5d。为消除假 频信号所需相位项是 $m=\frac{\frac{N}{2}+1}{2}--------\left(6\right)$ 这里N=信号采样块长度。

E-TDAC也要求将一对精心设计的分析-合成窗应用于重迭 的信号采样块。信号采样块必须有100％重迭，即给定块的50％ 被前一块重迭，而同一块的另50％被它的后续块重迭。图6a- 6g说明了信号采样块的重迭和造成的假频畸变的抵消。图6a和 6d所示信号Yc(t)和Ys(t)是从逆DCT和DST恢复而 来的，它们分别被窗函数Wc(t)和Ws(t)调制(如图6b和 6e所示)，产生出示于图6c和图6f的信号Yc(t)和Ys (t)。当这些加窗信号的重迭块相加时，假频分量被抵消，所产生 的信号Y(t)(示于图6g)是原始输入信号X(t)的精确重建 信号。

下面更详细地讨论窗设计和合成过程中使用的重迭-相加。在这 一点上注意到下述情况就是够了：略挣变换块的二分之一使所需位速 率减半，但100％窗重迭(这是在信号合成过程中E-TOAC所 要求的)使所需位速率加倍。结果，E-TDAC是临界采样的，并 对所需位速率是有中性影响。

E-TDAC所用的计算算法是快速付立叶变换(FFT)。参 见Cooleq和Tukeq的“复付立叶级数的机器计算算法”，Math。 Comput.卷19，1965，297-301页。单次FFT能用 来同时完成两个等长度信号采样块的修正DCT和修正DST，其作 法是规定每一块分别为一个单一复变换的实分量和虚分量。这一技术 利用了这一事实：FFT是一个复变换，而两个输入信号采样块只含 有实数值采样。通过将这些变换因式分解成一个FFT和一个复常数 阵到(arraq)的汞积，DCT系数从FFT中产生，作为一组实部 值，而DST系数由一组虚部值表示。所以，一个信号采样块的DCT 能与另一个长度信号采样块的DST同时计算，只需一个FFT后接 受数阵列相汞的相加。

使用一个FFT同时计算两个变换的基本技术是在本门技术中众 所周知的，并由Brigham在下列书中描述：速快付立叶变换， Englewood Cliffs，NJ，Prentice-Hall，Inc.， 1974。关于同时计算E-TDAC变换所需的修正的DCT和修 正的DST的附加信息可在Lookabaugh的博士论文中找到：“话 言的可变速率和自适应频率域失量化”Stanford，CA． Stanford Universitq，PhD Thesis，1988年6月。

在本发明的单通道方案的一个实施例中，两个相邻的重迭的最大 长度信号采样块被存放在缓存区并一同传送到一对DCT/DST块 中。如图7a所示，这个块对构成了一帧(frarne)。然后将这一 帧中的两块量化和格式化供传输或存贮。

在双通道系统中，通过处理来自两个通道中每一个的最大长度信 号采样块，可以实现同时计算两个块的变换：从一个通道产生一个 CCT块，从第二个通道产生一个DST块。这对块构成一帧。见图 7b。对于给定通道的编码块在DCT和DST之间交替变换(见表 正式5)，并且其类型总是与另一通道块的相反。在一帧中的两块被 一起变换和格式化。

然而，对于本发明的一个最佳实施例，正变换的同时处理更复 杂，因为本发明可以自适应地选择信号采样块长度。如下文中将要讨 论的，块长的适应性变化也可以要求改变正变换相位项(见方程6) 及变换长度。可以用来切换块长度的基本方法将在下文中描述。

B．瞬变信号检测

图1a中框110代表的瞬变信号检测器监测输入信号振辐的快 速变化，并在检测到足够大的瞬变信号时选定短的信号采样块长度。 如上文所讨论的，量化误差将使离散变换编码模糊了整个信号采样块 区间内瞬变信号的频谱分量。例如，比较图8a、8b和8c。在这 些图中的信号图比例尺使得清楚地显示出瞬变信号前面的低水平信 号。瞬变信号的最大振辐超过了图的范围，因而未显示出来。图8a 给出变换编码前的原始信号。图8b给出使用512点采样块长度从 非适应性块长度变换编码恢复出来的信号。瞬变信号的高频分量调制 了整个块内的音乐信号，但在瞬变信号之前最明显。

在使用短信号采样块的编码中这一畸变是无害的，因为该畸变包 含在较短的块中，并且在音质上被瞬变信号本身掩蔽。当瞬变信号检 测器确定需要用较短的块来保证瞬变信号畸变产物有暂时音质掩蔽 时，瞬变信号检测器可以选择较高的时间分辨力，其代价是由于自适 应选择较短块长而使变换滤波器族的选择性变差。如图8C所示，对 于本发明的一个移到256点采样块的实施例，其恢复得到的信号中 瞬变信号之前的高频瞬变分量的延续时间要短于图8b所示使用 512点采样块来编码的情况。

通常，在包含瞬变信号的信号采样子块内，滤波器族的选择性降 低不会造成损害，因为这种较差的滤波器族性能将被瞬变信号本身所 掩蔽。然而，在瞬变信号之前或之后的较短子块中编码性能的降低是 可以被听出来的，所以，最大限度地减小其编码使用较短块长的时间 区间是重要的。

瞬变信号检测器由四部分组成：第一部分是一个高通滤波器 (HPF)，它从瞬变信号检测处理过程中排除了较低频率分量。HPF 是用四阶送归数字滤波器来实现的。图9给出这个滤波器的示意图， 它可以表示成方程7和8所给出的一对差分方程。其等效的传递函数 于方程9。 Wn=b01Xn+b11Xn-1+b21Xn-8-a11Wh-a-

   a21 Wn-2                                       (7) Yn=b02Wn+b12Wn-1+b22Wn-2-a12Yn-1-

   a22 Yn-2                                         (8) 于使用1024点采用最大块长的一个本发明最佳实施例，HPF 的载止频率是2KHq。这个载止频率是通过听音试验建立的。表Ⅰ 给出差分方程系数值和传递函数值。一个使用512个采样点最大块 长的本发明实施例将使用截止频率为4KHq的高通滤波器(HPF)。 表Ⅱ给出这一滤波器截止频率时的系数。

瞬变信号检测器的第二部分将高通滤波后的信号采样分段，构成 子块的层次结构子帧(hierachicol subframe)。在本发明 的一个实施例中，使用了三级层次结构。参见图10。1024采样 的最大块长的二分之一被置于第一级。在第二级，这个二分之一信号 采样块又分段成二个子块，每个为256点采样。在第三级，这个二 分之一信号采样块分段成四个子块，每个为128个采样点。一般而 言，子块长度可表示成 这里Nj=第j级子块的长度，

N=最大长度信号采样块的长度，

M=在该分段层次结构中的层次级别数目。 尽管O级没有被瞬变信号检测器所使用，但从概念上讲，它代表最大 长度信号采样块。固此N0=N。

瞬变信号检测器的第三部分是一个峰值振幅检测器。对于在当前 子帧内的所有层次结构级别中的每个子块，识别出具有最大振幅的信 号采样。子块峰值振幅表示为 这里Xn=N/2采样长度的子块中第n个信号采样

j=子帧层次结构级别序号

k=在级别j中的子块序号

N=最大长度采样块的长度。 符号Pjo用来代表当前子帧中第j级第一个子块之前紧邻的前一个子 帧中第j级最后一个子块的峰值振幅。参见图11。例如，在前一子 帧中的P34表示成当前子帧中的P30。

瞬变信号检测器的第四部分是一个阀值比较器。这一部分的第一 级防止瞬变信号检测器在所通过的信号中只含有低振幅信号时改变块 长度。这一功能的实现办法是：对于其峰值P11不超过一个“默阈值 (Silence threshold)”Ts的所有信号采样块，将块长置为 其最大值。这一阈值表示为

Ts=0.003Xs    (12) 这里Xs=信号采样Xn的饱和值。这是一个信号采样在超出其编码表 示范围之前所能得到的最大值。例如，对于24位二进制整数的二进 制补码表示，其饱和值是8，388，607。

如果峰值P11不超过默阈值，则瞬变信号检测器的剩余部分被旁 路。否则，如果存在瞬变信号，则比较器的第二级一个较小的信号采 样块长度来适应于采样信号中存在的瞬变信号。如果首先给出一个概 念性描述，会对比较器第二级和第三级所完成的处理有更好的理解。 从概念上讲，第二级计算出一套该层次结构子帧的每一级相邻子块峰 值振幅之比。这些比值表示成 这里j=子帧层次结构级别序号

k=在第j级的子块序号

Pjk=子块峰值振幅(见方程11)。 将每个比值Rjk与两个阈值TAj和TDj相比较。TAj表示第j级的 增长阈值(attack fhreshold)，它总量大于1。TDj代表第 j级的衰减(decay)阈值，它总是小于1。表Ⅲ给出这两个阈值 之值。如果峰值振幅比向上或向下超过(Cross)了其中的一个阈 值，则瞬变信号检测器选择一个较短的信号采样块长度。这个自适应 过程由第三级完成，将在下面几段中描述。

比较器的第三级构成一个二进制值节点树，如图12所示。在树 结构中的每个节点代表一个在当前子帧中可能被选中的一个可能的子 块长度。这个树中每个节点之值由下列表达式确定： 这里j=子帧层次结构级序号，

k=第j级内的子块序号，

Rjk=子块峰值振幅比(见方程13)。 然后，修改树的节点。从树的“叶”节点开始(在本发明的最佳实施 例中为第3级)，对于每个等于1的节点，它的父节点和所有侧节点 也都是1。

比较器第二级和第三级的实际实施与上文中的概念性描述的差别 仅仅在于数学表达式的形式不同。应该承认，当Pj(k-1)为0时在 方程13中表示的比值没有定义。这一问题的解决办法是避免除法， 重新定义表达式14为下面的形式： 这里j=子帧层次结构级序号，

k=第j级内的子块序号，

Pjk=子块峰值振幅(见方程11)。 图13中的流程图说明了构成和修改结构树的逻辑。图14给出一个 结构树在修改前和修改后的实例。

比较器的第四级从“树”上削去所有的零节点。如果该结构树只 是由零节点组成，瞬变信号检测器将指出：如果可能的话帧控制器 (图1a所示框114代表的)将被允许使用最大长度采样块。否 则，余下的叶节点确定子块长度。对于第j级叶节点的子块长度是方 程10中确定的Nj。例如参考图14，经过删简过程之后留下来的 叶节点是t31、t33和t22。在本发明的最佳实施例中，这些子块中 每一个的长度将是图15所示的。图16提供了为含有一个瞬变信号 波形的信号段所选定的子块长度的实例。

请注意，对于使用E-TDAC变换的本发明最佳实施例，由瞬 变信号检测器构成的树的叶节点所代表的子块长度之和总是等于N／2 采样点数。这个总数等于每次通过瞬变信号检测器过程中向前推进的 信号采样数目。因此，瞬变信号检测器只需对每个信号采样处理一 次。

本领域的专家应该理解，可以使用其他瞬变信号检测器和子块长 度选择方法(Scheme)而不偏离本发明的范围和精神。

再有，在下文中描述的本发明的一个最佳实施例中，是根据由本 发明的瞬变信号检测器部分所指定的信号采样块长度来选择窗函数和 变换类型的。窗函数或变换类型可以直接由瞬变信号特性来选择而不 偏离本发明的范围或目的。

C．加窗

因为变换假定块中信号是周期性的，所以除非信号采样块被修 正，否则一个离散变换将错误地造成不存在的频谱分量。参见图 17。这些变换误差是由于图18所示块边缘的不连续性造成的。可 以对这些不连续处进行光滑最大限度地减小这种效应。图19a至 19d中以图来说明一个块如何被修正或加权，从而使块边缘附近的 采样值接近于零。图19a所示乘法器电路以图19c所示加权函数 来调制图19b所示的采样输入信号x(t)。结果信号示于图19 d。这一过程由图1a中的框112代表。这个加权函数称作分析 窗，是对信号采样块的采样点逐个相乘，它已是相当多研究工作的主 题，因为它的形状对数字滤波器的性能有深广影响。例如，可参考 Harris的文章“论使用离散付立叶变换进行谐合分析时窗的应 用”，Proc.IEEE，卷66，1978，51-58页。简单 地说，对于给定的阻带抑制水平，一个好的窗允许较陡的过渡带滚 降，并允许对该窗的调制效应进行校正。下文中将更详细地讨论窗的 设计。

窗的选择也受到由瞬变信号检测器选定的信号采样块长度的影 响。因为编码要适应于根据瞬变信号状态而定的信号采样块长度，编 码也必须选择具有适当形状和长度的分析窗。根据由瞬变信号检测器 选定的信号采样块长度来选择分析窗的过程是由图1a中框114所 代表的帧控制器完成的。下面将更详细地讨论帧控制器。

应该注意的是，在瞬变信号检测器要求减小块长时，为子块选择 的窗可能不是最佳的，因为它可能是不对称的，也可能在它的两个边 缘不趋于零。如上文所讨论的那样，使用一个欠佳(Suboptimal) 的窗将会降低正变换和逆变换中的滤波器选择性。尽管如此，由于下 述两个理由，还必须选择这种欠佳的窗：(1)在本发明的最佳实施 例中使用的正变换和逆变换对窗的形状有一定要求，以便能消除时间 域假频畸变，(2)组合分析/合成窗的形状必须保证相邻重迭窗在 重迭区间之和为单位1。

对于一个给定的窗函数W(下文中将更详细讨论一个最佳窗函 数)，记号Wab将用于代表一个“满(full)”窗，它开始于 零或接近零，在Na/2采样点之后上升到值1或接近1，在Nb／2 采样点之后其值降到零或接近零。窗Wab的总长度可以表示为 这里记号Nx=在第X级的子块长度(见方程10)。如果a等于b， 则窗Wab是对称的。与这一记号相一致，用于调制最大长度信号采 样块的满长度窗将记为Woo。

部分窗W以零或接近零值开始，在Na/2采样点之后上升到一 个等于1或接近1的终止值，该部分窗用记号Wa+表示。以1或接 近于1的值开始并在Na/2采样点之后终止于零或接近零值的部分 窗由记号Wa-表示。

图20中给出一组窗。窗W00是由窗W0+和W0-联合构成 的曲线代表。窗W32是由窗W3+和W2-联合代表。

D．帧控制

图1a中的框114代表的过程是根据瞬变信号检测器的输出来 控制选择分析窗和选择变换。这些过程控制处理一个信号帧所需的活 动，将在下文中描述。还将简要讨论用单个FFT同时完成两个变换 时的若干考虑。

1．固定帧排列(Fixed-Frame Alignement)

在一个使用E-TDAC变换的本发明最佳实施例中，固定帧排 列(FFA)不仅允许对一对满长度信号采样块同时处理变换，而且 允许对子块对同时处理变换。FFA产生不变长度的帧，所以帧边界 的排列是固定的或者说是同步的(Synchronous )。然而，将FFA 限定到为最大信号采样块长度二分之一的选择的子块长度。所以，供 FFA使用的瞬变信号检测器要比上文描述的检测器简单，其简单之 处在于它的层次结构子帧只由一级构成。

图21a至21c图形表示在一个单通道系统中的一系列加窗信 号采样块。块A和块B代表已由单个FFA同时变换过的两块。如果 不是由于发生了瞬变信号，那么块c和块D本来也已经同时变换过 了。

在图21a中给出一个瞬变信号，它出现在块G的后一半(子块 Cb)和块D的前一半(子块Da)当中。编码器只需对块C和块D 中这四个子块转移到N/2(512)点变换。注意，对于这些子块 使用的窗与最大长度来样块使用的相同，即为W00：用于子长Ca 和Da的窗W00的前一半是W0+，用于子块Cb和Db的窗W00 的后一半是W0+和W0-是欠佳的(suboptimal)，因为它们 在两个边缘不趋于零。如图22a所示，使用欠佳窗会使滤波器族的 选择性明显降低。虽然在含有瞬变信号的子块内较差的编码器性能一 般将被瞬变信号本身所掩蔽，但在瞬变信号之前和之后的较短长度子 块中这种性能降低仍可能被听到。所以，希望在瞬变状态平息之后尽 可能快地恢复使用最大长度信号采样块。

用于E-TDAC抵消假频的FFA需要两个不同的相位项。 图23a给出的两个相邻重迭满长度信号采样块，是在利用相邻块重 迭相加以抵消时间域假频之前，从逆修正DST恢复得来的。被恢复 的信号由两个分量组成：加窗的原始信号(实线)和时间域假频畸变 (虚线)。假频分量是加窗原始信号的时间逆转复制品，但时间逆转 发生在两个分开的区域。E-TDAC变换的相位项m(见等式1和 2)控制这两个区域之间的边界位置。对于正常(normal)E- TDAC，其边界位于信号采样块的中点。所需相位项示于等式6。

图23b给出的两个相邻重迭变换块是从本发明的一个实施例在 其移到较小子块长度时利用FFA恢复的。满长度块已从逆修正DCT 恢复了。半长度块已从逆修正DST恢复。在DCT块中的假频分量 与前述相同。但是，如果要利用相邻块重迭相加来消去假频信号，则 第一个DST子块中的假频分量必须是从一端到另一端完全时间逆转 的。对于所有子块都要求其假频分量是一端到另一端逆转的。适当的 相位项是 $m=\frac{1}{2}------\left(16\right)$ 可以说明，相位项可写成更一般的形式： 这里φ=时间逆转区域之间边界的位置。这个公式代表的情况比 FFA所要求的更具有一般性，但它需要下文中讨论的其他帧控制技 术。

图23c说明一个信号采样块被分成两个子块并被两个窗调制。 右侧的子块和窗的长度是N/4采样点。这一子块内的两区域之间的 边界位于该子块从右侧(或者说尾部边缘)算起N/8采样点处。在 每个N/4采相子块区域内引起假频分量时间逆转所需要的相位项是 $m=\frac{\frac{N}{8}+1}{2}++\left(18\right)$ 这里N=最大信号采样块长度。

对于图21a所示实例，一个单一的FFT能够同时处理子块 Ca的DCT和子块Da的DST。类似地，单一FFT能同时处理子 块Cb的DCT和子块Db的DST。能对块E和F恢复最大长度块的 同时处理。

图21b所示实例只与图21a所示实例稍有不同。图21b中 所示的瞬变信号在块D的后一半(子块Db)中发生。在这种情况 下，由单个FFT来同时处理变换要求在一个比图21a所示情况要 长的区间内使用较短的子块长度。因此，由于滤波器族选择性降低所 引起的较差的编码性能要影响到较长的区间。此时可对子块对 Ca/Da、Cb/Db、Ea/Fa及Eb/Fb进行同时处理DCT/DST 块对。可在块G和H恢复最大长度块的同时处理。为消除假频信号E -TDAC所要求的各子块相位项与上面公式16给出的相同。

图21c所示实例与图21b所示相同，但处理顺序稍有不同。 对图21c所示实例，满长度块C是单独变换的。子块对Da／Ea和 Db/Eb是同时变换的。能对块F和G恢复最大长度块的同时处理， 然而，现在这个块对的前一个块是进行DST变换。这一顺序允许本 发明的一个实施例更快地恢复最大块长编码，然而需要一个更快的处 理器，因为对块C到块下进行变换所需要的计算量比使用固定长度块 同时变换的编码序列所需要的计算量大20％左右。

2．扩展的固定帧排列

本发明的一个最佳实施例使用的一种技术叫做扩展的固定帧排列 (Enhanced Fixed-Frame Alignment)(EFA)，它 产生一个固定的或同步的帧排列，但允许使用的子块长度为最大信号 采样块长度的2的幂次约数。在本发明的一个最佳实施例中，EFA 技术可以使用的信号采样块长度为1024，512，256及128 采样点。(如下文中所示，在某些限定情况下将使用64个采样点的 块长。)然而，在实际情况中这个限制并不是严格的，只是因为用 FFT实现离散变换时对于块长为2的幂次的情况更为有效。

在能够充分描述EFA之前，有必要引入“桥变换”(bridge trangform)的概念。桥变换是为从一个信号采样块长度向另一个 采样块长度过渡来搭桥的一种变换。例如，如图24所示，假定要求 本发明去处理一个N/2采样点的子块，接着再处理另一个N/4采 样点的子块。这时对每个子块分别完成变换将是可能的。然而，注意 到N/2采样子块要加窗W0+，而N/4子块要加窗W1-。两个窗 显然都是欠佳的，会引起滤波器选择性的显著降低。例如可参见图 22a。

一个桥变换通过避免使用Wj+和Wj-之类窗来使滤波器选择性 的降低程度最小。对于图24所示实例，桥变换允许一个单一变换使 用一个W01窗，而不需要两个变换各自分别使用W0+和W1-窗。 从效果上讲，桥变换完成了3N/4采样长度的一个单一子块的变 换。对W01窗的滤波器响应曲线示于图22b。

图24也给出了为消除时间域假频由桥变换所需要的时间逆转区 域。两个时间逆转区域的边界位于含有要进行桥变换的信号采样的那 两个子块的彼此边缘处。这个边界位于两个搭桥子块的右手边缘起算 N/4采样点处。因此，根据公式17，对于N/2到N/4的桥变 换所需相位项为 $m=\frac{\frac{N}{4}+1}{2}$ 这里N=最大信号采样块长度。

图24所示单个桥变换的计算可以利用一个FFT对3个N/4 长子块计算其变换，然后再进行重新组合操作。这一技术是本领域中 公知的。可参见Oppenheim和Schafer的“书”数字信号处理” Engleloood Cliffs，N.J.：Prentice-Hall， Inc.，1975，307-314页。FFT加上重组合操作 也能用于同时处理两个E-TDAC桥变换。为实现同时处理所需要 的预处理和后处理与Brigham和Lookabaugh针对正常满长度 E-TDAC变换所描述的桥变换需要的预处理和后处理相同。然而 要注意的重要一点是：在E-TDAC中的同时处理只是对于具有相 同长度和相同相位项的修正DCT和修正DST才是可能的。

利用同时处理允许本发明的一个实施例采用较慢的处理器。然 而，对于本发明使用ETA的单道通方案，同时处理几乎没有得到什 么好处。对于使用同时处理的实施例，所选择的处理器应该能够对付 最坏情况的信号瞬变事件。全部满长度块的分别变换只给处理器稍微 增加了一点负担。所以，与利用EFA和同时处理来实现本发明所需 要的处理器相比。利用EFA而不用同时处理来实现本发明的单通道 方案只需要稍微快一点的处理器。

对于双通道系统，如果两个通道总是用同时一块长来编码，那么 用同时处理是有好处的。这可以通过从同一个二进制树来为道通选择 子块长度的办法实现。这个树可以通过对两个单独树的相应节点进行 布尔“或”组合来构成，这里的两个单独的二进树是通过它们各自的 瞬变信号检测器来处理各通信号构成的。

下文中对EFA的解释适合于本发明的单通道实施例，并假定不 进行同时处理。

图26a至26f所示流程图说明了能用于控制EFA的“帧控 制过程”的逻辑流程。图26a给出这一过程的全貌。其中框2600 代表的初始化步骤示于图26b。这一过程的主体由四个主要部分构 成。

图26c所示帧控制的第一步是设置子块长度。检查一个队列， 以确定是否在瞬变信号检测器的二进制树上有任何节点尚未被处理。 如果没有，瞬变信号检测器则去处理下一个N/2信号来样点，如上 述方式构成二进制树，并把每个叶节点的级别放入队列。从队列中取 出第一个条目，并指定为J1。这个值为“当前”子块规定了所希望 的长度。符号J0代表“前一个”子块的树节点级别。如果不曾检测 到瞬变信号，而且瞬变信号检测器已经指出可以使用N采样点块长， 即J1=0，那么就必须测试标志“PARTIAL”(部分)以确 定是否允许编码器恢复满长度处理。如果不可以，则迫使J1的取值 为1 ，置子块长度为N/2样点。

图26d所示帧控制的第二部分确定所需分析窗长度、变换长 度。及E-TDAC相位项。所选窗是Wxy，这里x是由值J0建立 的，而y是由值J1建立的。变换长度和相位项的确定方式与图26b 所示情况类似。

第三部分示于图26e。这一部分通过检验标志“TYPE(类 型)”来选择所需变换的类型。这一个二进位标志指出是否应完成修 正的DCT或修正的DST。

帧控制的第四部分示于图26f。如果J0和J1双双为零，那 么便知道刚刚完成了满长度变换。于是清除“PARTIAL(部 分)”和“BRIDGE(桥)”两个标志。否则，标志“BRIDGE (桥)”被触发。正在清除/触发步骤之前，这个标志指出是否曾使 用过正常变换或桥变换作为当前变换。标志“TYPE(类型)”被 触发，且当前子块长度J1被够入J0。

如图26g所示，如果标志“BRIDGE(桥)”是清除的， 则处理过程继续其第一部分，得到下一个子块长度。否则，PARTIAL (部分)标志被触发，处理过程继续其第二部分。

图25给出单通道系统中加窗信号采样块序列。所示瞬变信号发 生在块C的后一半和块D的前一半中。瞬变载号检测器为这一子帧建 立的二进制树示于图14。为每个子块指定的相应长度适于图15。 编码器对块C的处理是首先以W03窗进行9N/16采样长度的桥 变换。其次，使用W33窗完成N/8(12)采样长度的一个变 换。再其次，以W32窗使用长度3N/16采样的桥变换。在块C中 的剩余N/8采样留给下一个子帧，与块E的开始部分一起进行桥变 换。    

编码器对块D的前一半进行二个长度为N/8(128)采样的 变换和一个N/4(256)采样长度的变换。一个W33窗用于两个 N/8采样子块，一个W22窗用于N/4采样子块。块D的后一半用 一个N/2长度变换加W11窗。在子块D使用的窗都是最佳的。

使用桥变换和非对称窗允许E-TDAC利用由公式17计算的 相位项来消掉时间域假频信号。这一个采样块序列所需要的时间逆转 区域及相应的相位项示于图25。

3．动态帧排列(Dynamic-Frame Alignment)

进行E-TDAC变换的动态帧排列能够选择任意子块长度L， 只要L是偶数而且大于或等于4。然而由FFT来实现的本发明实施 例只使用最大信号采样块长度的2的幂次约数作为子块长度。这是因 为对于变换长度为2的幂次的情况完成FFT最为有效。对于大多数 情况，DFA不允许对子块对进行同时的变换处理，但DFA的确改 善了编码器性能，以其允许在瞬变状态过去之后编码器立即返回最大 长度信号取样块而使编码器性能超过了FFA或EFA所可能达到的 编码器性能。

DFA与EFA很相似。这两种技术之间的主要区别在于：与 EFA不同，DFA不限制于保持固定的或同步的帧排列。DFA可 以自由地截断当前块帧并立即转回满长度信号采样块。与此相反， EFA必须保证使用适当的子块长度来结束当前帧，然后才能恢复满 长度块的处理。

DFA帧控制过程只在少数几个方面与EFA帧控制不同。将参 考图26a、26c及26f所示流程图来解释这些差别。第一，标 志PARTIAL(部分)不象图26a中框2660所示那样被触 发，而是永远是零。第二，不需要象图26c中框2611所示那样 进行测试来确定树节点队列是否为“空”，因为这个队列永远是空 的。第三，不象图26f中框2644所示那样触发标志BRIDGE (桥)，而是永远把该标志置零。

DFA瞬变信号检测器只有一点与EFA检测器不同：对于每次 重复，DFA检测器并不总是前进N/2个新采样点。假定瞬变信号 检测器对当前子块选定了要转向(Shift)子块长度N/8采样， 那么在处理下一个子帧之前瞬变信号检测器只前进N/8信号采样 点。这样，与EFA的情况不同，DFA瞬变信号检测器可能检验一 个给定信号采样不只一次。

图27给出一个单通道系统中的加窗信号采样块序列。图中给出 一个瞬变信号，发生在块C的后半段。下文中的讨论假定瞬变信号检 测器选择了移向长度为N/8或者说128个采样的两个子块。

编码器首先开始用W03窗长度为9N/16采样的桥变换来处 理块C。块C是一个被截断的块，它以长度为N/8(128)采样 的变换使用W33窗结束。长度为9N/16采样的桥变换加W30窗 用于块E之中，再返回满长度块处理。

编码器用两个长度N/8(128)采样的变换来处理被截断的 块D。在被截断的块D之后立即开始满长度块F。请注意块帧边界已 移动了满块长的3/4。

F．系数量化和格式化

由正变换产生的频率系数一般不适于低位速率传输或有效存贮。 多种量化技术可以使用，利用信号的冗余性和不相干性 (irrelevancy)来减少所需位数。例如，可参考Jayant和 Noll的“波形的数字编码”，Englewood Cliffs，N.J： prentice-Hall，Inc.，1984，563-576页。这 一过程由图1a中框118代表。

如果一个信号成份能被接收器只利用信号中的其余部分便可以预 测出来或提供出来，那么这个信号成份便是冗余的(redundant)。 在当今技术中公知的利用信号冗余的技术有线性预测编码、自适应预 测编码和熵(entropy)编码。

一个信号成份如果被略掉了确不降低所接收信号的可察觉质量， 则这个信号成份是可有可无的(irrelevant)。已有技术中有几 种可以利用信号的可有可无性，包括对数量化、浮点表示、以及基于 频率成份量或音质掩蔽效应的自适应位分配(bit allocation)。

格式化过程把量化的变换系数与信号采样块长度组合在一起以供 传输或存贮。这一过程由图1a中的框120代表。在本发明的一个 最佳实施例中，一个具体的子块长度序列由所有可能的子块长度序列 表中的索引号来代表。例如，表中一个条目将代表子块长度N/2、 N/8、N/8及N/4采样序列。另一个条目代表8个N/8采样 子块构成的序列。

帧同步位和错误检验/纠错代码可以使用，供格式化数据传输时 需要。数据库指针或关键字可以加入，供格式化数据存贮时需要。现 在格式化帧已经准备好，可以沿图1a所示路经122传送或存贮。

DFA技术给信号格式化增加了更大的复杂性。DFA格式化的 一个优选方案是利用循环缓存区。当代表变换系数的适当数量的数据 位与附加信息已经放入缓存区时，那些数据可以是格式化的，加上帧 同步位，然后传输或存贮。这允许格式化帧有一个不变长度。固此， 使用DFA时的格式化帧不再等效于变换块帧。DFA格式化帧不保 证含有相同数量的变换块。固此，由解码器接收的变换块长度信息的 准确性是至关重要的，因而应该采取某种办法(例如几余代码或错误 校正代码)来加以保护。变换块长度信息出错会使解码器不能对那个 变换块及随后的所有块进行适当的解码，直到重新建立变换块同步为 止。

              G．格式解析和系数线性化

信号通道132或者接收被传送来的信号或从存贮器取出信号， 当这样接收到数字化的编码信号时，便开始了格式解析 (deformatting)过程。格式解析过程由图1b中的框134代 表。格式解析过程提取解码器使用的量化变换系数和信号采样块长 度。利用在编码器中对系数进行量化过程的逆过程来改变换系数转换 成线性表示形式。这一过程由图1b中框136代表。线性采样块长 度用于选择适当的合成窗函数和设置逆变换长度。

              H．合成滤波器族-逆变换

图1b中框138代表一族合成滤波器，它们将格式解析和线性 化处理所恢复的频率域系数变换时间域信号采样块。图1a中的分析 滤波器族116所用变换的逆变换实现了合成滤波器族138。本发 明的这一实施例中使用的E-TDAC技术采用的逆变换是交替应用 修正的逆DCT和修正的逆DST。每个变换的长度设置成等于从格 式化信号中提取出来的每个信号采样块长度。所需相位项的计算采用 的处理过程与上文描述的本发明编码器部分使用的相同。

因为变换系数中有一半从传输或存贮中略掉了(见表达式3和 4)，所以必须为逆变换重新建立那些系数。失掉的DCT系数可以 从所能得到的DCT系数按公式19所示重新建立。失掉的DST系 数可以按公式20所示重新建立。逆修正的DST由公式21表示， 逆修正DST由公式22表示。 这里k=变换系数序号，

n=信号采样序号， K=变换系数总数， N=采样块长度， m=E-TDAC所用相位项(见公式6)， (k)=量化的第k个DCT系数， (k)=量化的第k个DST系数， (n)=恢复的量化信号x(n)。

图4a-4e及6a-6g说明了分析-合成滤波器族的变换过 程。分析滤波器族将时间域信号变换成交替出现的DCT和DST块 序列。逆变换对每隔一个块应用逆DCT，而对这些块中的另一半应 用逆DST。如图5a-5d所示。所恢复的信号中含有假频畸变。 这一畸变是在其后进行的时间域块重迭-相加过程(图1b中框144 代表这一过程)中被消去的。下文中讨论这-重迭-相加过程。

逆变换的计算是由FFT算法完成的。在逆变换中采用了正变换 中使用的同样技术，以允许使用单个FFT同时计算一对逆修正 DCT和逆修正DST。

图1b中框140所代表的帧控制过程控制子块对的排列及逆 O-TDAC变换为消除时间域假频信号所需相位项的选择。这一过 程与本发明编码器部分使用的帧控制过程相同。

Ⅰ．合成窗

图6a-6g说明了通过对相邻时间域信号采样块进行重迭-相 加来消除时间域假频信号的过程。正如Princen导出的那样，为消 除时间域假频畸变，E-TDAC变换要求应用一个与分析窗完全相 同的合成窗并进行相邻块的重迭-相加。每个块重迭100％；50％ 被前一块重迭，50％被后续块重迭。图16中框142代表了合成 窗调制过程。分析窗和合成窗的设计必须满足下列条件：当相邻的乘 积窗被重迭时，该分析一合成乘积窗总是相加为单位1。这个重迭- 相加过程由图1b中的框144代表，并示于图6a-6g。分别从 逆DCT和DST恢复的信号Ya(t)和Ys(t)示于图6a和 6d。每个信号被分成一系列块。每个信号块被图6b和图6e所示 合成窗函数调制。结果生成的信号块Yc(t)和Ys(t)示于图 6c和6f。以二分之一块长度重迭的两个信号相加，产生出信号 Y(t)，示于图6g。信号Y(t)是原始输入信号的精确重建。

如图28所示，在块k和块k+1之间的重迭区内某些时刻not 上的信号采样是由这两个块的每一个中的采样来代表的。在两个加窗 块重迭一相加之后，可以看到在时刻not处的恢复信号采样为取自 加窗块k和k+1的采样之和，它可由下式表示：

x(not)=WPk(not)·x(not)+WPk+1(not)·x(not )

                                                                  (23) 这里WPk(not)=WAk(not)·WSk(not )=(WAk(not))2，

WAk(not)=块k中在时刻not的分析窗，

WSk(not)=块k中在时刻not的合成窗，

WAk(not)=WSk(not)，这是E-TDAC变换所要求的。

如果两个相邻乘法窗在其窗重迭区间之和等于单位1，则乘法窗 的调制效应被抵消。所以，如果对块k和块k+1之间的重迭区间内 所有时间采样，满足

WPk(nt)+WPk+1(nt)=1，    (24) 则信号x(nt)可以被精确恢复。

虽然具有适当的重迭-相加性质的任何窗都可用作导出分析-合 成窗对的基础，但在本发明的一个最佳实施例中使用的窗是从Kaiser -Bessel窗导出的，其α值在4至7范围内。参见公式25。 这里α=Kaiser-Bessel   α因子， n=窗采样序号， N=以采样点数表示的窗长度， $I_0=\underset{k=0}{\overset{\infty }{\Sigma }}\frac{{(x/2)}^k}{k!}$ 窗的导出是用一个长度等于块长减重迭区间的矩形窗与一个 Kaiser-Bessel窗进行褶积。参见公式26。这一公式可以简 化成公式27。 n=乘法窗采样序号 V=窗重迭区间内采样数， N=乘法窗的预期长度， W(n)=长度为V+1的开始窗函数， WP(n)=长度为N的导出乘法窗， 对于使用E-TDAC变换的本发明最佳实施例，分析窗和合成 上面导出的乘法窗WP(n)取平方根得到的。分析和合成窗有 长度1024、512、256和128采样点，有100％窗 或N/2采样)。分析窗函数示于公式28。 N∈(128，256，512，1024)，

W(n)长度为N+1的Kaiser-Bessel函数；α因子

取质范围是4至7。

J．信号输出

图1b中框146代表一个常规的数-横转换器，它响应数字输 入，产生一个可变电压横拟信号。数字信号取自重迭-相加过程产生 的24位整数字的16个最高有效二进制位。横拟输出应由一个通带 宽度为15KHz(对20KHz编码器为20KHz)的低通滤波 器来滤波。该滤波器示于图1b中的框148。所恢复的适于放大的 声音信号可以图1b中所示输出终端150得到。

Ⅱ．本发明的另一个O-TDAC实施例

本发明的另一个实施例利用的变换称作“奇迭时间域抵消假频 (Oddly-Stacked Time-Domain Aliasing Cancellation)(O-TDAC)。下面的描述讨论本发明的 E-TDAC和O-TDAC方案在实现中的差别。

A．正变换

O-TDAC利用的变换函数是修正的离散余弦变换(DCT)， 示于等式29。 这里k=频率系数序号， n=输入信号采样序号， N=采样块长度， m=O-TDAC所用相位项(见等式6)， x(n)=输入信号x(t)在采样n的量化值， C(k)=第k个DCT系数。 O-TDAC变换产生一组谱系数或变换块，其形式是： 这里i=信号采样块序号，

C(k)=DCT系数(见等式29)。

所用的计算算法是快速付立叶变换(FFT)。与E-TDAC 方案不同，O-TDAC的实现不能用单一FFT来同时变换两个信 号采样块。然而，利用与E-TDAC方案中采用的“先乘-变换- 后乘处理”相类似的技术可以把变换的复杂性从N2减小为 N·log3N计算。前乘步骤将信号采样x(n)的实数值序列转变成 一个复数，其值序列是用复数函数 $e^{-\mathrm{j\pi }\frac{n}{N}}-------\left(31\right)$ 来乘信号采样点。这里， n=输入信号采样序号， N=采样块长度。

由FFT实现的离散付立叶变换将修改过的信号采样转换成一组 变换系数。因为FFT是复数变换，修正过的信号采样组的实部和虚 部能被同时变换。最后，后乘(Postmultiply)步骤得到DCT 系数。这一过程由下面的等式32和33表示： 这里

n=输入信号采样序号，

N=采样块长度，

k=频率系数序号，

m=O-TDAC的相位项(参见等式6)，

R(k)=系数X*(k)的实部，

Q(k)=系数X*(k)的虚部，

C(k)=第k小DCT系数。

在本发明的单通道方案最佳实施例中，每个信号采样块用FFT 分别变换。一个块构成一帧。在双通道系统中，从两个通道中的每一 个得到的信号采样块由FFT处理变换的DCT1/DCT2块对。 这个块对构成一帧。

Princeh指出适当的相位分量m(见式b)和精心设计的分析 -合成窗对，O-TDAC技术能从形如

{C(K)}0，{C(K)}1，{C(K)}2

{C(K)}3，…    (34) 的一组余弦变换块精确地恢复输入信号。

O-TDAC变换和假频信号消除过程与E-TDAC变换的情 况很相似，这一过程示于图29a-29e、30a-30d及 31a-31g。其主要区别是假频信号分量的形式。对于采样信号 块的前一半，假频分量是输入信号相对于采样编四分之一点处的时间 逆转图象，但它的振幅的符号与输入信号的符号相反。对于采样信号 块的后一半，其假频信号是相对于采样块四分之三点的时间逆转，而 其符号不变。参见图30b和30d。

分析窗和合成窗的设计和使用都与E-TDAC的情况完全相 同。参见图31a-31g。上文中对E-TDAC讨论的帧排列技 术在这里也同样有效。

在O-TDAC中使用适应性信号采样块长度时，与-TDAC 的情况所需要的考虑相比，并不需要再作其他额外的考虑，只是FFA 技术所用的变换相位项的选择这一点除外。图32中给出消除假频信 号所需要的相位项，图中显示出信号采样块为C至E，这是由逆 O-TDAC变换恢复而来的一组三个重迭的满长度变换块。根据 FFA技术，块D和E已被分段成为子块。

为了利用重迭-相加来实现假频信号的消除，在子块Da和Do 中的假频分量必须在时间上逆转，并且其振幅要从加窗原始信号那里 翻转得来。请注意，在图31a中这一逆转/翻转特性是对于O- TDAC在左手区由假频分量产生出来的。所以，逆转区之间的边界 必须放在该子块右手边缘处。由式17，我们知道子块Da和Ea所 需相位项是 $m=\frac{1}{2}-------\left(35\right)$

在子块Db和Eb中的假频分量只是在时间上逆转。这一特性对 O-TDAC是在右手区表现出来的。所以，逆转区之间的边界中必 须放在该子块的左手边缘。根据式17，所需相位项应是 $m=\frac{\frac{N}{2}+1}{1}------\left(36\right)$ 这里N=最大信号采样块长度。

对于O-TDAC使用FFA的确定相位项的一般规律如下。对 于每个满长度块中的第一子块，其变换使用的相位项如式35所示。 对于所有的满长度块以及在每个满长度块中的第二子块，其相位项如 式36所示。

B．逆变换

DCT系数中有一半从传输或存贮中略掉，然后利用式37所示 关系从所能得到的DCT块中重新建立起来。逆修正DCT示于式 38。 这里k=变换系数序号，

n=信号采样序号，

K=变换系数数， N=采样块长度， m=E-TDAC的相位项(见式6)， (k)=第k个量化DCT系数， (n)=恢复的量化系数X(n)。

在应用O-TDAC实现逆变换时，利用与正变换中使用的相类 似的“前乘-变换-后乘”处理，使计算复杂性从N2减少到N· log2N。这一过程将实数值DCT系数变成一组修正的复数值系 数，利用逆FFT同时变换修正系数的实部和虚部(IFFT)，再 从“后乘”得到时间域信号，如下列等式所示。 这里 m=O-TDAC的相位项(见式6)， N=采样块长度， k=频率系数序号， n=输入信号采样序号， r(n)=采样X*(n)的实部， q(n)=采样X*(n)的虚部， (n)=恢复的量化信号X(n)。

其后进行的加窗、重迭-相加、及信号输出处理与上文中描述的 本发明E-TDAC实施例的情况完全相同。

C．自适应使用E-TDAC和O-TDAC变换

由上面的讨论，一个熟悉本门技术的人应该理解，自适应地在 E-TDAC和O-TDC之间进行选择也是可能的，选择分析窗 和合成窗的过程与上面描述的相同。为消除时间域假频信号所需的相 位项的计算也将遵循上面讨论的相同原则。

变换选择可以用来改善瞬变信号的编码器性能或为其他感兴趣的 信号特征来改善编码器性能。例如，通过在这两个TDAC变换之间 来回改变，编码器/解码器系统能借助自适应地移动频率仓的中心频 率来改善单调(single-tone)信号的编码。这是可能的，因为 对这两种变换各自的频率仓带宽有差异。如图33a所示，E- TOAC变换仓零(直流分量)和仓N/2有带宽等于所有其他仓带 宽的二分之一。然而，对于O-TDAC变换，每个频率仓是等间隔 的。见图33b。在本发明的一个最佳实施例中使用1024采样点 的块长，采样率为44.1KHq，有可能例如将仓1的频率从 21.5Hq移到43Hq，这只需从E-TDAC变换转移到O- TDAC变换就可以实现。

Ⅲ．本发明的另一种DFT实施例

对于大多数应用，E-TDAC是最佳选择，然而，E-TDAC 编码器所需信号处理资源(resources)比以离散付立叶变换 (DFT)为基础的编码器所需资源要多。使用DFT实现编码器需 要较少存贮器和较低处理速度。本发明的DFT实施例的基本结构与 图1a和图1b所示相同，只是没有使用合成窗调制器142。

下面的描述讨论了本发明的DFT变换方案和E-TDAC变换 方案的区别。

A．瞬变信号检测器和帧控制

上文中描述的关于E-TDAC实施例的瞬变信号检测器逻辑在 概念上与DFT实施例的相同，只是其检测器分析整个块的采样而不 是只分析半个块的采样。固此，帧控制过程尽管在概念上与E- TDAC的相似，但一次选择各子块长度是针对一个满长度块，而不 是只对半块子帧长度来选择各子块长度。

B．加窗

因为E-TDAC变换与DFT变换之间的差异，其分析窗与 E-TDAC编码器使用的分析窗不同。选择一个能减少输入信号采 样块重迭量的窗设计是很重要的，因为DFT不是临界采样的，也就 是说，对于重迭区间的信号采样，其传输速率或数据存贮量都是要加 倍的。使用100％重迭的分析窗会使DFT编码器需要将近两倍于 E-TDAC编码器所需的位速率。

这样，在本发明的一个DFT实施例中使用的满长度块优选窗显 示出在大区间上有一倍增益，从而把块重迭长度由N/2(512) 采样减小到16个采样。这一减小降低了数字滤波器的阻带抑制水 平，但与E-TDAC编码器的情况相比它只引起数据速率的适度增 加。

从DFT窗族中导出一个窗的方式与E-TDA实施例中导出 汞积窗WP(n)所用方式类似。对于DFT实施例，核(Kernel) Kaiser-Bessel函数长度为17采样点，其α固子之值在 1.5至3范围内。参见式25。DFT分析窗是由Kaiser- Bessel窗与一个矩形窗被积得到的。对于窗W00，该矩形窗的长 度是1008采样点(块长1024点减去重迭区间的16个点)。 类似地，对窗W33其矩形窗的长度是112采样(或者说128- 16)。将这些值代入式27，得到的这两个窗函数是： 这里W(n)=Kaiser-Bessel函数，其长度为17采样点， 其α因子取值区间是1.5至3。

C．分析滤波器族-正变换

DFT实现其滤波器族，表达式是： X(k) x(n)e 对0≤n＜N    (43) 这里 k=频率系数序号， n=输入信号采样序号， N=信号采样块最大长度， x(n)=输入信号x(t)在采样点n处的量化值， X(k)=第k个变换频率系数。

D．合成窗族-逆变换 DFT实现其逆滤波器族，表达式是： x(n)= X(k)e 对0≤K＜N(44) 这里 k=频率系数序号， n=输入信号采样序号， N=信号采样块最大长度， x(n)=在采样点n处的恢复的最化信号x(t)。 X(k)=第k个变换频率系数。

E．合成窗 如上所述，DFT不需要使用合成窗。

F．其他变换

可以使用任何正交变换而不偏离本发明的精神和目的。利用上面 对于DFT所讨论的各种考虑，根据瞬变信号的存在或任何其他信号 特性来自适应地选择任一正变变换也是可能的。例如，本发明的一个 实施例可以在正常情况下使用式45所示通常DCT，但选择式46 所示通常DST来编码含有瞬变信号的子块。其逆DCT和逆DST 分别示于式47和48。 X(k) (k) x(n)cos (45) X(k) x(n)sin (46) x(n) α(k)X(k)cos (47) x(n) α(k)X(k)sin (48) 这里K=频率系数序号，

n=输入信号采样序号，

N=最大信号采样块长度，

x(n)=输入信号x(t)在采样点n处的量化值，

X(k)=第k个变换频率系数， a(k)= 对＜＞0

1    否则。

   表：

       a1j           a2j  ai1    -1.641731     0.678759  ai2    -1.809110    0.853238

       b0j           b1j         b2j  bi1    0.830277     -1.659936  0.830277  bi2    0.941365     -1.878653  0.941365

表Ⅰ  2KHz高通滤收器系数

       a1j          a2j  ai1    -1.266438  0.416893  ai2    -1.517035  0.709569

       b0j          b1j        b2j  bi1   0.671333    -1.340665  0.671333  bi2   0.852569    -1.690380  0.852569

表Ⅱ  4KHz高通滤收器系数

                      比值阀 子帧级别号    子块长度   TAj      TDj

1            512     5.0   0.2000

2            256     10.0  0.1000

3            128     13.3  0.0752

表Ⅲ 瞬变信号检测器增大一衰减阀值