test_main.m

code efficace dct of with matlab

clc;clear;close all;
disp('EFFICIENT 8X8 DISCRETE COSINE TRANSFORM');
disp('In order to obtain the source code please visit');
disp(' ');
disp('http://www.advancedsourcecode.com/efficientdctccode.asp');
disp(' ');
disp('DCT and IDCT for arbitrary size images are available on request.');
disp('For any question please email me [e-mail]');
disp(' ');
disp('Luigi Rosa');
disp('Via Centrale 35');
disp('67042 Civita Di Bagno');
disp('L''Aquila - ITALY');
disp('mobile [phone]');
disp('email [e-mail]');
disp('website http://www.advancedsourcecode.com');
disp(' ');

a = 255*rand(8,8);

% NOTE: if you work with uint8 images you have to use dct_uint8 and
% idct_uint8

y1 = dct2(a);
y2 = dct_double(a);

figure,mesh(abs(y1-y2)),title('Difference between DCT2 and C-code');


y3 = idct2(y1);
y4 = idct_double(y1);

figure,mesh(abs(y3-y4)),title('Difference between IDCT2 and C-code');