calculardct.m

program that compares the DCT of two matrix

function [Mdct] = calcularDCT(M)

%Obtiene la transformada DCT de una matriz de 8x8 

Mdct=[];
M = cast(M,'int32');
for x=1:8
    for y=1:8
        if x == 1
            Cu = 1/sqrt(2);
        else 
            Cu = 1;
        end
        if y == 1
            Cv = 1/sqrt(2);
        else 
            Cv = 1;
        end        
        S=0;
        for i=0:7
            for j=0:7
                S=S+M(i+1,j+1)*cos((((2*i+1)*(x-1)*pi)/16))* ...
                    (cos(((2*j+1)*(y-1)*pi)/16));
            end
        end
        Mdct(x,y)=(1/4)*Cu*Cv*S;
    end
end
return
end