Subversion Repositories mdct

x = [

   139,144,149,153,155,155,155,155;

   144,151,153,156,159,156,156,156;

   150,155,160,163,158,156,156,156;

   159,161,162,160,160,159,159,159;

   159,160,161,162,162,155,155,155;

   161,161,161,161,160,157,157,157;

   162,162,161,163,162,157,157,157;

   162,162,161,161,163,158,158,158;

   %255,144,149,153,155,155,155,155;

   %255,151,153,156,159,156,156,156;

   %255,155,160,163,158,156,156,156;

   %255,161,162,160,160,159,159,159;

   %254,160,161,162,162,155,155,155;

   %254,161,161,161,160,157,157,157;

   %254,162,161,163,162,157,157,157;

   %254,162,161,161,163,158,158,158;

];

% level shift input colors. This will range DC coeff to boundary 9 bits

xn = x-128;

%Cm = [

%   C4,  C4,  C4,  C4,  C4,  C4,  C4,  C4;

%   C1,  C3,  C5,  C7, -C7, -C5, -C3, -C1;

%   C2,  C6, -C6, -C2, -C2, -C6,  C6,  C2;

%   C3, -C7, -C1, -C5, -C5,  C1,  C7, -C3;

%   C4, -C4, -C4,  C4,  C4, -C4, -C4,  C4;

%   C5, -C1,  C7,  C3, -C3, -C7,  C1, -C5;

%   C6, -C2,  C2, -C6, -C6,  C2, -C2,  C6;

%   C7, -C5,  C3, -C1,  C1, -C3,  C5, -C7;

%   ];

%Cz = [

%   A,  A,  A,  A,  0,  0,  0,  0;

%   B,  C, -C, -B,  0,  0,  0,  0;

%   A, -A, -A,  A,  0,  0,  0,  0;

%   C, -B,  B, -C,  0,  0,  0,  0;

%   0,  0,  0,  0,  D,  E,  F,  G;

%   0,  0,  0,  0,  E, -G, -D, -F;

%   0,  0,  0,  0,  F, -D,  G,  E;

%   0,  0,  0,  0,  G, -F,  E, -D;

%   ];

Ye = [

   0;

   0;

   0;

   0;

];

Yo = Ye;

Y = [

   zeros(1,8);

   zeros(1,8);

   zeros(1,8);

   zeros(1,8);

   zeros(1,8);

   zeros(1,8);

   zeros(1,8);

   zeros(1,8);

];

Y = dct_func(xn);

Z = dct_func(Y);

Z

%1/2*Cm*x;

Y;

Q = [ 16,11,10,16,24,40,51,61;

   12,12,14,19,26,58,60,55;

   14,13,16,24,40,57,69,56;

   14,17,22,29,51,87,80,62;

   18,22,37,56,68,109,103,77;

   24,35,55,64,81,104,113,92;

   49,64,78,87,103,121,120,101;

   72,92,95,98,112,100,103,99

];

N=8;

% quantization

for i = 0:N-1

   for j = 0:N-1

      Zq(i+1,j+1) = round( Z(i+1, j+1)/Q(i+1,j+1));

   end

end

Zq;