function fdct = fast_dct2D(img)
%快速DCT變換，其根據文獻《一中簡便快速的DCT算法及其硬件實現》
%文獻中提到這種DCT變換可以將復雜的乘法計算分解為簡單的移位和加法計算
%從而提升計算速度

%%  程序
T=dctmtx(8);                           % 這里用了MATBLE里面的函數生成矩陣
B=blkproc(I,[8 8], 'P1*x*P2',T,T');    % 根據文獻， 做DCT快速變換
I2=blkproc(B2,[8 8], 'P1*x*P2',T',T);
figure,imshow(I2);                     % 重建后的圖像如圖3-4所示