离散余弦变换DCT

离散余弦变换 DCT

先说结论，DCT 变换就是输入信号为实偶函数的 DFT 变换。

对于一个长度为 N 的实数离散信号： x [ n ] = { x [ 0 ] , x [ 1 ] , . . . , x [ N − 1 ] } x[n]=\{x[0],x[1],...,x[N-1]\} x[n]={x[0],x[1],...,x[N−1]}。其 DFT 可以表示为：
$$X[k]=\sum _{n=0}^{N-1}x[n]e^{-in\frac{2\pi }{N}k}\text{\hspace{1em}(1)}$$

先将其扩大为原来的 2 倍，则新的信号$x^{\prime }[n]$ 为：
$$\begin{gathered} x^{\prime }[n]=x[n],where0\le n\le N-1 \\ {x}^{\prime }[n]=x[-n-1],where-N\le n\le -1 \end{gathered}$$
由于其并不关于 $n=0$ 对称，因此，需要将其右移 $\frac{1}{2}$ 个单位。其对应的 DFT 改变成：
$$X^{\prime }[k]=\sum _{n=-N+\frac{1}{2}}^{N-\frac{1}{2}}{x}^{\prime }[n-\frac{1}{2}]e^{-in\frac{2\pi }{2N}k}$$
进一步：
$$X^{\prime }[k]=\sum _{n=-N+\frac{1}{2}}^{N-\frac{1}{2}}{x}^{\prime }[n-\frac{1}{2}]\cos (n\frac{2\pi }{2N}k)-\sum _{n=-N+\frac{1}{2}}^{N-\frac{1}{2}}{x}^{\prime }[n-\frac{1}{2}]\sin (n\frac{2\pi }{2N}k)$$
因为 $x^{\prime }[n-\frac{1}{2}]$ 是偶函数，$\sin (n\frac{2\pi }{2N}k)$ 是奇函数，${\sum }_{n=-N+\frac{1}{2}}^{N-\frac{1}{2}}{x}^{\prime }[n-\frac{1}{2}]\sin (n\frac{2\pi }{2N}k)=0$，则：
$$\begin{gathered} X^{\prime }[k]=\sum _{n=-N+\frac{1}{2}}^{N-\frac{1}{2}}{x}^{\prime }[n-\frac{1}{2}]\cos (n\frac{2\pi }{2N}k) \\ =2\sum _{n=\frac{1}{2}}^{N-\frac{1}{2}}{x}^{\prime }[n-\frac{1}{2}]\cos (n\frac{2\pi }{2N}k) \end{gathered}$$
令 $m=n-\frac{1}{2}$，最终：
$$X^{\prime }[k]=2\sum _{m=0}^{N-1}{x}^{\prime }[m]\cos (\frac{(m+\frac{1}{2})\pi k}{N})$$
在实际的工程中，为了在 DCT 变换变成二维矩阵运算的形式时，将该矩阵正交化以便于进一步的计算，那么，可以乘以一个系数 $\sqrt{\frac{1}{2N}}$ （$k=0$ 时除外）：
$$X^{\prime }[k]=\sqrt{\frac{2}{N}}\sum _{m=0}^{N-1}{x}^{\prime }[m]\cos (\frac{(m+\frac{1}{2})\pi k}{N})$$