Subversion Repositories mdct

function output = dct_func(input)
 
N = 8;
 
input_s = input./2;
 
A = round(cos(pi/4)*(2^11))
B = round(cos(pi/8)*(2^11))
C = round(sin(pi/8)*(2^11))
D = round(cos(pi/16)*(2^11))
E = round(cos(3*pi/16)*(2^11))
F = round(sin(3*pi/16)*(2^11))
G = round(sin(pi/16)*(2^11))
Ce = [
   A,  A,  A,  A;
   B,  C, -C, -B;
   A, -A, -A,  A;
   C, -B,  B, -C;
];
 
Co = [
   D,  E,  F,  G;
   E, -G, -D, -F;
   F, -D,  G,  E;
   G, -F,  E, -D;
];
for i=1:N
 
   fXe = [
      input_s(i,1);
      input_s(i,2);
      input_s(i,3);
      input_s(i,4);
   ];
 
   fXo = [
      input_s(i,8);
      input_s(i,7);
      input_s(i,6);
      input_s(i,5);
   ];
   fXs = fXe+fXo;
   fXm = fXe-fXo;
 
   fYe = (Ce*fXs);
   fYo = (Co*fXm);  
   % transpose output
   output(1,i) = fYe(1,1)/2^11;
   output(2,i) = fYo(1,1)/2^11;
   output(3,i) = fYe(2,1)/2^11;
   output(4,i) = fYo(2,1)/2^11;
   output(5,i) = fYe(3,1)/2^11;
   output(6,i) = fYo(3,1)/2^11;
   output(7,i) = fYe(4,1)/2^11;
   output(8,i) = fYo(4,1)/2^11;  
end
 
return