fandct.m

自己手写的图像的FFT变换和DCT变换的MATLAB代码

%**********************************************
%************* 二维的DCT变换  %****************
%************************************
function uu=fandct(y) 
[wideth,height]=size(y);
temp1=zeros(1,height);
%**********************************************
%*************  进行行的的DCT变换  %****************
%************************************
    for i=1:wideth
        for j=1:height
            temp1(j)=y(i,j);            
        end;
        temp1=iyDCT(temp1,height);
        for j=1:height
            uu(i,j)=temp1(j);
        end;
    end;
%**********************************************
%*************  进行列的DCT变换  %****************
%************************************
     temp2=zeros(1,height);
    for j=1:height
        for i=1:wideth
            temp2(i)=uu(i,j);            
        end;
        temp2=iyDCT(temp2,wideth);
        for i=1:wideth
            uu(i,j)=temp2(i);
        end;
    end;
end