DCT一维变换原理及实现,二维dct变换,matlab

clear; clc; X=round(rand(32)*100) %产生随机矩阵 A=zeros(32); for i=0:31 for j=0:31 if i==0 a=sqrt(1/32); else a=sqrt(2/32); end A(i+1,j+1)=a*cos(pi*(j+0.5)*i/32); end end Y=A*X*A' ; %DCT变换 YY=dct2(X); Y8=Y(1:8,1:8); M=mean(Y8); N=mean(M);