DCT算法的MATLAB代码实现

function [F]=DCT（A）
   A=im2double（A）;
   [mn]=size（A）;
   F=zeros（mn）;
   S=zeros（1m*n）;
   k=2;
   for u=1:m
       for v=1:n
           for a=1:m
               for b=1:n
                   S（11）=A（11）*cos（（（2*1+1）*u*pi）/（2*n））*cos（（（2*1+1）*v*pi）/（2*m））;
                   if k==2
                       S（1k）=A（ab+1）*cos（（（2*a+1