两种快速DCT算法的矩阵分解与分析.pdf

( ) 第 31卷 第 2 期 湘潭师范学院学报 自然科学版 Vo l. 3 1 N o. 2 　　　　　　 ( ) 　　　　　 2009年 6月 J o u rn a l o f X ia n g ta n N o rm a lU n ive rs ity N a tu ra l S c ien ce Ed itio n J un e. 2009 两种快速 DCT算法的矩阵分解与分析 1 2 2 刘 　琼 ,欧阳万里 , 肖创柏 ( 1. 湖南涉外经济学院 ,湖南 长沙 4 10205; 2. 北京工业大学 计算机学院 ,北京 100 124 ) 摘 　要 :利用矩阵乘法理论来分析 Loeffler DCT算法和 Feig DCT算法 。通过使用矩阵分解的表示形式 ,指出了两种算 法的区别与联系 ,这种矩阵分解的表示形式和分解过程有利于对算法的理解和进一步提出更好的快速算法 。 关键词 : DCT;矩阵分解 ;快速算法 中图分类号 : TP39 1. 4 1　　　文献标识码 : A 　　　文章编号 : 167 1 - 023 1 (2009) 02 - 0024 - 03 ( ) 在图像信号处理 ,特别是如 JPEG,M PEG1/2 /4和 H. 263等图像和视频压缩标准中 , 8 ×8 离散余弦变换 DCT 是不可 或缺但又是计算复杂度非常高的。许多快速算法的提出就是为了降低算法的计算复杂度 [ 1 - 5 ] 。文献 [ 1 ]中提出的 Loeffler DCT算法由于其规整和简单的特性而被广泛使用 ,文献 [ 1 ]使用数据流图描述了该算法 ;另一个广泛应用的算法是 Feig 的 2D - Scaled DCT算法 。一般来说 , Feig算法运算复杂度最低 。然而 , Feig没有说明怎样获得该算法 , 只给出了结果 。文献 [ 6 ]中对 Feig的 Scaled - DCT算法进行了 3种不同形式的拓展 。本文通过使用矩阵分解的表示形式 ,指出在不同的分解阶 段 ,两个算法都乘上了一个对角阵 。这种矩阵分解的表示形式和分解过程有利于对 Feig 的 2D - Scaled DCT算法的理解并 通过矩阵乘法理论对其做进一步改进 。 1　8 点 DCT的初步分解模型 N 点 D C T - II变换可以定义为 : N - 1 π ( ) 2 n 2k + 1 ( ) y = c co s x 。 1 n n ∑ 4N k k =0 其中, n, k = 0, 1, 2,