close?all
clear?all
clc
%-----------原始數據以及參數設置---------------%
Y=imread（‘D:\matlab\bin\tv.BMP‘）;
N=length（Y）;
figure
subplot（131）
imshow（Y）;
title（‘original?picture‘）;
Y=im2double（Y）;
imwrite（Y‘D:\matlab\bin\LENA.jpg‘）
xlabel（‘方差=??0‘）
T=?dctmtx（N）;
T1=dctmtx（8）;
p=0.01;
pause
%-------------directly?DCT---------------------%
Y1=T*Y*T‘;
%-------------directly?IDCT---------------------%
Y1=T‘*Y1*T;
subplot（132）
imshow（Y1）;
title（‘after?directly?DCT‘）;
t=mean（var（Y1-Y））;
xlabel（t）
pause
%----------8*8?block?