A= imread('F:lena.png','png');%读取图像信息
%I=imresize(A,[128,128]);
imshow(I); %显示原图
title('原图128*128');
I=im2double(I);
t=dctmtx(8);  
dct=@(x)t*x*t';  % 设置函数句柄
for i=2:12
for j=2:12
B=A(8*i-8:8*i,8*j-8:8*j);
C(i,j)=mean(mean(B));
end
end
%B=blkproc(I,[8 8],dct,'imresize'); % 图像块处理
%B=mat2cell(I,ones(8,1)*128/8,ones(8,1)*128/8);
%B =colfilt(I,[8,8],'distinct',dct);
%y1=blkproc(I,[8,8],'P1*x*P2',t,t');%进行DCT变换得到变换矩阵%放大后矩阵的大小
%y1=colfilt(I,[8,8],'P1*x*P2',t,t');
k=2;
M=round(k*8);
N=round(k*8);
%构造变换矩阵T
T1=[eye(8,8);zeros(M-8,8)];
T2=[eye(8,8);zeros(N-8,8)];
%构造目标图像频域数据
y2=blkproc(y1,[8 8],'P1*x*P2',T1,T2');
%乘以增强系数
y3=blkproc(y2,[M N],'P1.*x',k);
t1=dctmtx(8*k);
dct=@(x)t'*x*t;
y4=blkproc(y3,[M N],'P1*x*P2',t1',t1);
imshow(y4);