新数字信号处理 第2版 教学课件 张小虹 4数字信号处理5.pptVIP

* 在线教务辅导网： 教材其余课件及动画素材请查阅在线教务辅导网 QQ:349134187 或者直接输入下面地址： §4.8 离散余弦变换 在一般有限长正交变换中，其正反变换式可写为 式中?k(n)称基序列或变换核，它们相互正交，即 (4.8-1) A(k)= ? n=0 ? N?1 x(n)?k (n) (4.8-2) N 1 k=0 ? N?1 x(n) = A(k)?k(n) (4.8-3) ?k (n) ?m(n)= ? N 1 n=0 ? N?1 C 0 m=k m?k 如果C =1 是归一化正交变换 DFT就是这类变换之一。在DFT中，基序列是复周期序 序列。 也是复序列。实际上有若干种实数基序列（变换核）的 列WN 。在这种变换中，即使x(n) 是实序列，其变换X(k) 变换，使得x(n) 是实序列时，其变换A(k)=X(k) 也是实 压缩方面特别有用十分重要。 的特点，在数字信号处理的应用，尤其是语音、图像 DCT与DFT关系密切。由于DCT具有能量集中（压缩） 离散余弦变换(DCT)就是一种实序列的正交变换，且 4.8.1、离散余弦变换的定义 DCT的变换形式如(4.8-1) (4.8-2)式，其中变换核(基序 列)是余弦函数。 即与DFT隐含着周期性相似，DCT同时隐含周期对称 因为余弦函数既是周期的又是对称的， 使(4.8-2)式中 的x(n)在区间0≤ n ≤N?1外的延伸也是周期对称的。 性。由x(n) 的展开不同，有数种DCT的定义。 例 讨论x(n)的展开。 3 2 1 1/4 1/8 1/2 1 0 x(n) n n 1/4 1/2 1/4 1/8 1/2 1 0 1 2 3 4 5 1 1/2 1/4 … 6 7 8 9 ? x1(n) ? x1(n)的周期为2N?2=6 ，对应DCT-1 n 1/8 1 1/4 1/2 1/4 1/8 1/2 1 0 1 2 3 4 5 6 7 8 9 10 11 1 1/2 1/4 … ? x2(n)的周期为2N=8 ，对应DCT-2 ? x2(n) 3 2 1 1/4 1/8 1/2 1 0 x(n) n n ? x3(n) ? x3(n)的周期为4N=16 ，对应DCT-4 0 1 2 3 4 5 6 7 8 1 8 ? 1 4 ? 1 2 ? ? 1 1 8 1 4 1 2 13 14 15 9 10 11 12 1 8 ? 1 4 ? 1 2 ? 1 8 1 4 1 2 1 3 2 1 1/4 1/8 1/2 1 0 x(n) n 0 1 2 3 4 5 6 7 8 9 10 11 1 8 ? 1 4 ? 1 2 ? ?1 1 8 1 4 1 2 12 13 14 15 1 8 ? 1 4 ? 1 2 ? 1 8 1 4 1 2 1 ?1 1 n ? x4(n) ? x4(n)的周期为4N=16 ，对应DCT-4 3 2 1 1/4 1/8 1/2 1 0 x(n) n 对于给定的实序列x(n) ，0≤ n ≤N?1 ，DCT定义为 由x(n) 的展开不同，有数种DCT的定义，只讨论最常用 的一种，即DCT-2。 Xc2 (k) =DCT[x(n)] (4.8-5) 式中 c(k)= k=0 1 ? 2 1?k?N?1 1 (4.8-4) c(k) x(n)cos 2N (2n+1)k? n=0 ? N?1 2 ? N = 反变换 系数是归一化正交变换所需要的 x(n)=IDCT[Xc(k)] 在本节中 Xc2 (k)或X c(k)为简便记为X (k)。 1?n?N?1 (4.8-6) c(k)Xc(k)cos k=0 ? N?1 2 ? N = 2N (2n+1)k? ┇ 变换可以用矩阵表示为 Xc=CNx 其中 CN是N×N的变换矩阵，例N=8 由DCT的定义计算Xc(k) Xc(0)= [x(0)+ x(1)+ ???+x(N?1) ] 1 ? N Xc(1)= x(0)cos +x(1)cos + ???+x(N?1)cos 2 N 2N ? 2N 3? 2N (2N?1)k? ? Xc(N?1)= x(0)cos +x(1)cos