% Script A1DCT
% Checks out DCT8 and IDCT8
clc
% The DCT8 matrix...
C8 = cos((pi/16)*((0:7)'*(1:2:15)));
% Correctness...
x = randn(8,5);
errDCT8  = norm(DCT8(x)-C8*x);
errIDCT8 = norm(IDCT8(DCT8(x))-x);
fprintf('|| DCT8(x)-C8*x ||       = %10.3e\n',errDCT8)
fprintf('|| IDCT8(DCT8(x)) - x || = %10.3e\n',errIDCT8)








