高数复习--什么是DCT

  

  离散余弦变换英语:discrete cosine transform, DCT)是与傅里叶变换相关的一种变换,类似于离散傅里叶变换,但是只使用实数。离散余弦变换相当于一个长度大概是它两倍的离散傅里叶变换,这个离散傅里叶变换是对一个实偶函数进行的(因为一个实偶函数的傅里叶变换仍然是一个实偶函数),在有些变形里面需要将输入或者输出的位置移动半个单位(DCT有8种标准类型,其中4种是常见的)。

最常用的一种离散余弦变换的类型是下面给出的第二种类型,通常我们所说的离散余弦变换指的就是这种。它的逆,也就是下面给出的第三种类型,通常相应的被称为"反离散余弦变换","逆离散余弦变换"或者"IDCT"。

有两个相关的变换,一个是离散正弦变换,它相当于一个长度大概是它两倍的实奇函数的离散傅里叶变换;另一个是改进的离散余弦变换,它相当于对交叠的数据进行离散余弦变换。

  离散余弦变换,尤其是它的第二种类型,经常被信号处理和图像处理使用,用于对信号和图像(包括静止图像和运动图像)进行有损数据压缩。这是由于离散余弦变换具有很强的"能量集中"特性:大多数的自然信号(包括声音和图像)的能量都集中在离散余弦变换后的低频部分,而且当信号具有接近马尔可夫过程的统计特性时,离散余弦变换的去相关性接近于K-L变换(Karhunen-Loève变换——它具有最优的去相关性)的性能。

例如,在静止图像编码标准JPEG中,在运动图像编码标准MJPEG和MPEG的各个标准中都使用了离散余弦变换。在这些标准制中都使用了二维的第二种类型离散余弦变换,并将结果进行量化之后进行熵编码。这时对应第二种类型离散余弦变换中的n通常是8,并用该公式对每个8x8块的每行进行变换,然后每列进行变换。得到的是一个8x8的变换系数矩阵。其中(0,0)位置的元素就是直流分量,矩阵中的其他元素根据其位置表示不同频率的交流分量。

一个类似的变换, 改进的离散余弦变换被用在高级音频编码,Vorbis和MP3音频压缩当中。

离散余弦变换也经常被用来使用谱方法来解偏微分方程,这时候离散余弦变换的不同的变量对应着数组两端不同的奇/偶边界条件。

  公式定义:

形式上来看,离散余弦变换是一个线性可逆函数F:R^n
ightarrow R^n其中R实数集,或者等价的说一个n	imes n方阵。离散余弦变换有几种变形的形式, 它们都是根据下面的某一个公式把n个实数x_{n-1}}x_0,ldots ,x_{n-1}变换到另外n个实数f_0,ldots ,f_{n-1}的操作。

DCT-I

f_m = frac{1}{2} (x_0 + (-1)^m x_{n-1}) 
     + sum_{k=1}^{n-2} x_k cos left[frac{pi}{n-1} m k 
ight]

有些人认为应该将x_{0}x_{n-1}乘以{sqrt {2}},相应的将f_0f_{n-1}乘以frac{1}{sqrt{2}}。这样做的结果是这种DCT-I矩阵变为了正交矩阵(再乘一个系数的话),但是这样就不能直接和一个实偶离散傅里叶变换对应了。

一个n=5的对实数abcde的DCT-I型变换等价于一个8点的对实数abcdedcb(偶对称)的DFT变换,结果再除以2(对应的,DCT-II~DCT-IV相对等价的DFT有一个半个抽样的位移)。需要指出的是,DCT-I不适用于n<2的情况(其它的DCT类型都适用于所有的整数n)。

所以,DCT-I暗示的边界条件是:x_k相对于k=0点偶对称,并且相对于k=n-1点偶对称; 对f_{m}的情况也类似。

DCT-II

f_m =
   sum_{k=0}^{n-1} x_k cos left[frac{pi}{n} m left(k+frac{1}{2}
ight) 
ight]

DCT-II大概是最常用的一种形式,通常直接被称为DCT。

有些人更进一步的将f_0再乘以frac{1}{sqrt{2}}(参见下面的DCT-III型的对应修改)。这将使得DCT-II成为正交矩阵(再乘一个系数的话),但是这样就不能直接和一个有半个抽样位移的实偶离散傅里叶变换对应了。

所以,DCT-II暗示的边界条件是:x_k相对于k=-frac{1}{2}点偶对称,并且相对于k=n-frac{1}{2}点偶对称; 对f_{m}相对于m=0点偶对称,并且相对于m=n点奇对称。

DCT-III

f_m = frac{1}{2} x_0 +
   sum_{k=1}^{n-1} x_k cos left[frac{pi}{n} left(m+frac{1}{2}
ight) k 
ight]

因为这是DCT-II的逆变换(再乘一个系数的话),这种变形通常被简单的称为逆离散余弦变换。

有些人更进一步的将x_{0}再乘以{sqrt {2}}(参见上面的DCT-II型的对应修改),这将使得DCT-III成为正交矩阵(再乘一个系数的话),但是这样就不能直接和一个结果有半个抽样位移的实偶离散傅里叶变换对应了。

所以,DCT-III暗示的边界条件是:x_k相对于k=0点偶对称,并且相对于k=n点奇对称; 对f_{m}相对于m=-frac{1}{2}点偶对称,并且相对于m=n-frac{1}{2}点偶对称。

DCT-IV

f_m =
   sum_{k=0}^{n-1} x_k cos left[frac{pi}{n} left(m+frac{1}{2}
ight) left(k+frac{1}{2}
ight) 
ight]

DCT-IV对应的矩阵是正交矩阵(再乘一个系数的话)。

一种DCT-IV的变形,将不同的变换的数据重叠起来,被称为改进的离散余弦变换

DCT-IV暗示的边界条件是:x_k相对于k=-frac{1}{2}点偶对称,并且相对于k=n-frac{1}{2}点奇对称;对f_{m}类似。

DCT V~VIII

上面提到的DCT I~IV是和偶数阶的实偶DFT对应的。原则上,还有四种DCT变换(Martucci, 1994)是和奇数阶的实偶DFT对应的,它们在分母中都有一个{displaystyle n+{frac {1}{2}}}的系数。但是在实际应用中,这几种变型很少被用到。

最平凡的和奇数阶的实偶DFT对应的DCT是1阶的DCT(1也是奇数),可以说变换只是乘上一个系数{displaystyle a}a而已,对应于DCT-V的长度为1的状况。

反变换

DCT-I的反变换是把DCT-I乘以系数{displaystyle {frac {2}{n-1}}}。 DCT-IV的反变换是把DCT-IV乘以系数frac{2}{n}。 DCT-II的反变换是把DCT-III乘以系数frac{2}{n},反之亦然。

离散傅里叶变换类似,变化前面的归一化系数仅仅是常规而已,改变这个系数并不改变变换的性质。例如,有些人喜欢在DCT-II变换的前面乘以sqrt{frac{2}{n}},这样反变换从形式上就和变换更相似,而不需要另外的归一化系数。

每日名言: 自知者不怨人,知命者不怨天,怨人者穷,怨天者无志-----------荀子

参考文档:

1 https://zh.wikipedia.org/wiki/%E7%A6%BB%E6%95%A3%E4%BD%99%E5%BC%A6%E5%8F%98%E6%8D%A2 

原文地址:https://www.cnblogs.com/dylancao/p/9144700.html