一维离散余弦变换（DCT)及其反变换（IDCT）公式的一种推导

前言

  取一维离散Fourier变换与反变换的公式为
$$\begin{array}{c}F\left(u\right)=c\sum _{n=0}^{N-1}f\left(n\right)e^{-j\frac{2\pi un}{N}},u=0,1,\cdots ,N-1\\ f\left(n\right)=d\sum _{u=0}^{N-1}F\left(u\right)e^{j\frac{2\pi un}{N}},n=0,1,\cdots ,N-1\\ cd=\frac{1}{N}\end{array}$$
相关文章可参考一维离散傅里叶变换(DFT)和逆变换(IDFT)公式的一种推导

DCT公式初步推导

  基于公式
$$F\left(u\right)=c\sum _{n=0}^{N-1}f\left(n\right)e^{-j\frac{2\pi un}{N}},u=0,1,\cdots ,N-1$$

一般形式

  令$$f\left(n\right)=g\left(n\right)e^{-j\frac{p\pi u}{qN}},p,q\ne 0$$其中$g\left(n\right)$通过个数为$tN\left(t\in {\mathbb{Z}}_{\ge 1}\right)$的一个空间域（时域）实数采样序列：
$$g\left(0\right),g\left(1\right),\cdots ,g\left(N-1\right),g\left(N\right),\cdots ,g\left(tN-1\right),\forall n,g\left(n\right)\in \mathbb{R}$$
且当$t>1$时， ∀ n ∈ { N , N + 1 , ⋯   , t N − 1 } ,   g ( n ) = 0 \forall n\in \left\{ N,N+1,\cdots ,tN-1 \right\},\text{ }g\left( n \right)=0 ∀n∈{ N,N+1,⋯,tN−1}, g(n)=0，即在拥有前$N$个采样点的基础上新增$\left(t-1\right)N$个零采样点（无效采样点），则有：
$$\begin{array}{rl}& F\left(u\right)=c\sum _{n=0}^{tN-1}\left[\left(g\left(n\right)e^{-j\frac{p\pi u}{qN}}\right)\cdot e^{-j\frac{2\pi un}{tN}}\right]\\ & =c\sum _{n=0}^{tN-1}\left[g\left(n\right)\cdot \left(e^{-j\frac{p\pi u}{qN}}\cdot e^{-j\frac{2\pi un}{tN}}\right)\right]\\ & =c\sum _{n=0}^{tN-1}g\left(n\right)e^{-j\frac{\left(pt+2nq\right)u\pi }{qtN}}\\ & =c\sum _{n=0}^{tN-1}g\left(n\right)\left(\cos \frac{\left(pt+2nq\right)u\pi }{qtN}-i\sin \frac{\left(pt+2nq\right)u\pi }{qtN}\right)\\ & =c\sum _{n=0}^{N-1}g\left(n\right)\left(\cos \frac{\left(pt+2nq\right)u\pi }{qtN}-i\sin \frac{\left(pt+2nq\right)u\pi }{qtN}\right)\end{array}$$
  对$F\left(u\right)$取实部，得到计算$C\left(u\right)$的DCT公式
$$C\left(u\right)=\mathrm{Re}\left(F\left(u\right)\right)=c\sum _{n=0}^{N-1}g\left(n\right)\cos \frac{\left(pt+2nq\right)u\pi }{qtN},u=0,1,\cdots ,N-1$$

标准形式

  在DCT公式一般形式中，令$$p=1,q=t=2$$则有
$$\begin{array}{rl}& C\left(u\right)=c\sum _{n=0}^{N-1}g\left(n\right)\cos \frac{\left(2+2\times n\times 2\right)u\pi }{4N}\\ & =c\sum _{n=0}^{N-1}g\left(n\right)\cos \frac{\left(2n+1\right)u\pi }{2N},u=0,1,\cdots ,N-1\end{array}$$

IDCT公式初步推导

  基于公式
$$f\left(n\right)=d\sum _{u=0}^{N-1}F\left(u\right)e^{j\frac{2\pi nu}{N}},n=0,1,\cdots ,N-1$$

一般形式

  令
$$F\left(u\right)=C\left(u\right)e^{j\frac{p\pi u}{qN}}$$
其中$C\left(u\right)$通过个数为$tN\left(t\in {\mathbb{Z}}_{\ge 1}\right)$的一个频率域实数离散序列：
$$C\left(0\right),C\left(1\right),\cdots ,C\left(N-1\right),C\left(N\right),\cdots ,C\left(tN-1\right),\forall u,C\left(u\right)\in \mathbb{R}$$
且当$t>1$时， ∀ u ∈ { N , N + 1 ,   ⋯   ,   t N − 1 } , C ( u ) = 0 \forall u\in \left\{ N,N+1,\text{ }\cdots ,\text{ }tN-1 \right\},C\left( u \right)=0 ∀u∈{ N,N+1, ⋯, tN−1},C(u)=0，相当于在拥有前$N$个离散点的基础上新增$(t-$