本程序利用奇异值分解对3通道彩色图像进行压缩分解，具体步骤如下： 压缩过程： 1. 选取子图像大小K值，把图像分解成M×M个子图像，IMG(s),s=1,2,…, M2，其中M=N/K，原始图像IMG大小为N×N。 2. 计算这M2个子图像的平均值average，对每幅子图像减去均值图像得到新图像。 3. 计算相关矩阵R，其元素定义为 。 4. 计算R的特征值与特征向量，计算每幅子图像与最大特征向量的内积，便得到编码，即压缩后的图像。

继 Discrete Cosine Transform: Algorithms, Advantages, Applications (K. R. Rao and P. Yip) 之后的一本关于离散余弦变换(DCT)的快速实现算法的最新专著. 本书第五章“Integer Discrete Cosine/Sine Transforms”关于DCT整数近似实现写的非常详细和精彩, 很方便大家编程实现. 值得一提的是本书的第一作者Vladimir Britanak 一直在从事DCT, 尤其是改进的离散余