修改型离散余弦转换

修改型离散余弦转换（modified discrete cosine transform，下称缩写：MDCT）是一种衍生自傅立叶转换的转换。MDCT以第四型离散余弦转换（DCT-IV）为基础，但具有重叠转换的性质：它在较大资料集合中，用于处理连续资料区块，其中当前资料区块的后半段与下一个资料区块的前半段相重叠。这种重叠情形有助于避免资料区块边界产生多余资料，除了具有离散余弦转换的能量压缩特性外，也使MDCT尤其有利于数据压缩方面应用。借由以上优势，MDCT在音讯压缩领域中成为应用最为广泛的有损数据压缩手段，应用于当前大多数音讯编码格式，如MP3、杜比数位、Ogg Vorbis、AAC等等。
MDCT是由Princen、Johnson和Bradley于1987年提出，其基本原理即时域混叠消除法（time-domain aliasing cancellation，缩写：TDAC）是由Princen和Bradley于1986年提出。其他类似的转换还有如以离散正弦转换为基础的修改型离散正弦转换（modified discrete sine transform，缩写：MDST）以及其他较少使用的转换，例如以其他不同类型的DCT或DCT/DST的组合为基础的MDCT。
定义.
作为一种重叠转换，MDCT与其他傅立叶相关的转换比起有些不一样，因为它的输出数为输入数的一半（而非相同）。特别的是，它是一个线性函数"F" : R2"N" -> R"N"（其中R代表实数集合）。根据公式：
将2"N"个实数"x"0, ..., "x"2"N"-1转换成"N"个实数"X"0, ..., "X""N"-1。（在这个转换前面的正规化系数，单位不为一定且根据不同的解释方法而有所改变。下面只考虑MDCT跟IMDCT的正规化乘积。）
反转换.
MDCT之反转换为离散余弦反转换（inverse modified discrete cosine transform，缩写：IMDCT）。
MDCT由于输入数与输出数并不相同，乍看之下应不可逆。然而，透过加入造成错误并且必须要回复原始资料的子带重叠区块中重叠的IMDCT部份，MDCT可具备完全可逆性。这种技术称为时域混叠消除法。
IMDCT根据公式：
将"N"个实数"X"0, ..., "X""N"-1转换成2"N"个实数"y"0, ..., "y"2"N"-1。（如同第四型离散余弦转换为一正交转换，其逆转换具有与正转换一样的形式）。具有一般窗式正规化（usual window normalization,如下所示）的窗式MDCT（windowed MDCT），其逆转换的正规化系数必须乘2（即变成2/"N"）。
计算.
虽然直接应用MDCT公式的计算复杂度为O（"N"2），但是透过递回的方式将计算分解化，可以将计算复杂度降到O（"N" log "N"），并且得出一样的结果，就如同快速傅立叶转换一样。另外也可以透过其他转换计算MDCT。通常一个离散傅立叶转换（或快速傅立叶转换）或是一个离散余弦转换包含的前置处理跟后置处理的复杂度为O("N")。此外，如下所述，任何第四型离散余弦转换的演算法都能衍生出一种方法来计算偶数项的MDCT和IMDCT。
窗函数.
在一般的信号压缩应用中，转换特性可透过利用将窗函数"w""n"（"n" = 0, ..., 2"N"-1）乘以MDCT中的"x""n"及IMDCT中的"y""n"作进一步的改良，为了避免在"n" = 0 and 2"N"的边界点产生不连续并且使函数于这些点可以变得平滑而趋近零（即window the data "before" the MDCT and "after" the IMDCT）。原则上，"x"和"y"可以有不同的窗函数，而窗函数再资料区块转换时可以改变（尤其是当两个不同尺寸的资料区块结合在一起时）。但为简单起见，我们考虑给相同尺寸的资料区块应用的窗函数。 
转换仍然是可逆的（即TDAC依然可以正常作用），考虑一个对称的窗"w""n" = "w"2"N"-1-"n"，只要"w"满足Princen-Bradley condition：
各种不同应用的窗函数都是一样的，例如： 
为应用于MP3和MPEG - 2的AAC的窗函数，而 
为Vorbis所应用。此外AC-3与MPEG - 4 AAC均采用了Kaiser-Bessel derived (KBD) window。值得注意的是，应用在MDCT的窗函数与应用于其他类型的信号分析的窗函数是不同的。因为它们必须满足Princen – Bradley condition。原因之一，MDCT应用了两次窗函数，包括MDCT（分析）和IMDCT（合成）。
第四型离散余弦转换与原始TDAC的关系.
由定义可以看出，一个偶数点"N"的MDCT基本上就相当于一个输入位移了"N"/2且两个N点资料区块同时转换的第四型离散余弦转换。而再更进一步的检视这个等效关系，则重要的特性例如TDAC可以很容易地推导出来。为了更明确的定义与第四型离散余弦转换的关系，我们必须了解第四型离散余弦转换分别对应到偶数点跟奇数点的边界条件：偶数点在左边边界，（约在"n"=–1/2），奇数点在右边边界（约在"n"="N"–1/2）诸如此类（而不是DFT那样的周期性边界）。这些遵从以下所述：
和
因此，如果输入为一组长度为"N"的"x"阵列，我们可以想见这个阵列会扩增至（"x", –"x""R", –"x", "x""R", ...），其中"x""R"代表"x"的反向序列。考虑一个有2"N"个输入与"N"个输出的MDCT。我们将输入分成4组区块（"a", "b", "c", "d"），每组大小为"N"/2。如果我们把这些资料位移"N"/2（即MDCT的定义中的+"N"/2项），然后（"b", "c", "d"）扩增到"N"点第四型离散余弦转换的输入结尾。因此，我们必须将这些资料根据上述边界条件“叠”（“fold”）回来。是故，一个具有2"N"点输入（"a", "b", "c", "d"）的MDCT，就确实的相当于一个具有"N"点输入（–"c""R"–"d", "a"–"b""R"）的第四型离散余弦转换，其中"R"的意义与前述相同。（以这种方式，任何用来计算第四型离散余弦转换的演算法都可以很直观的套用在MDCT上）。同样地，上述IMDCT公式恰好正是第四型离散余弦转换的1/2，其中输出位移了"N"/2且扩增长度为2"N"（透过边界条件）。而第四型离散余弦逆转换可以简单的给定如上述之输入（–"c""R"–"d", "a"–"b""R"）。当透过边界条件完成位移跟扩增的动作后，我们可以得到： 
（因此有一半的IMDCT输出是多余的。）
我们现在明白了TDAC的运作原理。考虑计算一个MDCT的子带，有50％的重叠，2"N"个资料区块（"c", "d", "e", "f"）。则得到的IMDCT会跟上述类似：("c"–"d""R", "d"–"c""R", "e"+"f""R", "e""R"+"f") / 2。当加入这些到先前由重叠到一半所产生的IMDCT时，反向的项被消除了，并获得一简单的已回复的原始资料（"c", "d"）。
TDAC的起源.
”时间域混叠消除法”的起源现在很清楚了。逻辑上第四型离散余弦转换中根据边界条件所扩增的输入资料造成资料发生混叠的现象，正如同频谱中奈奎斯特频率会与较低的频率发生混叠一样。因此，当"c"–"d""R"及其他被加进来会有准确的符号（right sign）提供作消除之用。对于奇数点"N"（事实上很少使用），"N"/2不是整数，所以MDCT并不是由第四型离散余弦转换透过一个简单的位移就可以互相对应。在这种情况下，额外的位移一半的样本，代表著MDCT/IMDCT等同于第三/二型离散余弦转换，其分析类似上述。
窗式MDCT的TDAC.
如上所述，TDAC性质已被证明适用于普通MDCT，显示加入子带资料区块的IMDCT于重叠的一半部份可回复原始的资料。推导这项窗式MDCT逆性质较上述稍微复杂一些，叙述如下。根据前述当formula_8和formula_9经过MDCT与IMDCT，并且加上其重叠一半时，我们得到formula_10，即原始的资料。 
现在我们假设我们将MDCT的输入和IMDCT的输出两者均乘上一个长度为2"N"的窗函数。如前所述，我们假设一个对称的窗函数其形式为formula_11其中"w"跟"z"是长度"N"/2的向量，而"R"代表反向。如此则Princen – Bradley condition可以写成：formula_12，其中乘法和加法为点对点的运算，或等效成formula_13（将"w"跟"z"反向）。 
因此，我们对formula_14作MDCT取代对formula_8作MDCT（所有的乘法跟加法均为点对点）。当完成IMDCT再乘以窗函数（点对点相乘），最后"N"点的一半，成为：
（请注意，我们不再乘1/2，因为IMDCT的正规化与窗式的例子中的因数2有所不同）。同样的，对于产生formula_9的窗式MDCT和窗式IMDCT，是在其前"N"点的一半： 
当我们这两个各半的结果加总起来，我们得到： 
即回复的原始资料。
音讯编码应用.
就MP3格式而言，MDCT并不直接用于处理音讯讯号，而是用于处理32波段多相正交滤波器（polyphase quadrature filter，缩写：PQF）阵列输出端信号。MDCT输出结果会再经由一条混叠削减公式作后处理，以减弱PQF当中典型的混叠。这种转换与滤波器阵列组合，被称作混合滤波器阵列或子带MDCT。
AAC则相反，通常纯粹采用MDCT；其中唯一例外是索尼的罕见格式MPEG – 4 AAC - SSR——既采用MDCT，又在其后用到四波段PQF阵列。自适应听觉转换编码（adaptive transform acoustic coding，缩写：ATRAC）利用运用MDCT的堆叠型正交镜像滤波器（quadrature mirror filter，缩写：QMF）。