my_dct.m

其中

function C=my_DCT(A,k)
s=size(A);
for m=1:s(1)
    if rem(m,8)==0
        for n=1:s(2)
            if rem(n,8)==0
                B((m-8+1):m,(n-8+1):n)=dct2(A((m-8+1):m,(n-8+1):n));
                B((m-k+1):m,(n-8+1):n)=0;
                B((m-8+1):m,(n-k+1):n)=0;
            end
        end
    end
end
%
for m=1:128
    if rem(m,8)==0
        for n=1:128
            if rem(n,8)==0
                C((m-8+1):m,(n-8+1):n)=idct2(B((m-8+1):m,(n-8+1):n));
            end
        end
    end
end
% 
C=uint8(C);