myidct2.m

function P = myidct2(C)
% 2-D离散反余弦变换
% 输入矩阵P 输出变换结果矩阵C
N = size(C,1);
% 根据公式构造变换矩阵D
D = (1:2:2*N-1).*(0:1:N-1)';
D = cos(D*pi/(2*N));
D(1,:)=sqrt(0.5);
D = D*sqrt(2/N);
% 计算变换结果
P = D'*C*D;
end