x=imread('mj.jpg');
figure(5);imagesc(xgray);colormap(gray)
title('original')
x=double(x);
r=x(:,:,1);g=x(:,:,2);b=x(:,:,3);
xgray=0.2126*r+0.7152*g+0.0722*b;
xza=zeros(600,600);
for i=1:75
for j=1:75
xz=xgray(i*8-7:i*8,j*8-7:j*8);
xd=double(xz);
xc=xd-128;
y=dct(dct(xc')');
p=1;
q=p*8./hilb(8);
yq=round(y./q);
ydq=yq.*q;
xdq=idct(idct(ydq')');
xe=xdq+128;
xf=uint8(xe);
xza(i*8-7:i*8,j*8-7:j*8)=xf;
end
end
figure(1);imagesc(xza);colormap(gray)
title('p=1')
for i=1:75
for j=1:75
xz=xgray(i*8-7:i*8,j*8-7:j*8);
xd=double(xz);
xc=xd-128;
y=dct(dct(xc')');
p=2;
q=p*8./hilb(8);
yq=round(y./q);
ydq=yq.*q;
xdq=idct(idct(ydq')');
xe=xdq+128;
xf=uint8(xe);
xza(i*8-7:i*8,j*8-7:j*8)=xf;
end
end
figure(2);imagesc(xza);colormap(gray)
title('p=2')
for i=1:75
for j=1:75
xz=xgray(i*8-7:i*8,j*8-7:j*8);
xd=double(xz);
xc=xd-128;
y=dct(dct(xc')');
p=4;
q=p*8./hilb(8);
yq=round(y./q);
ydq=yq.*q;
xdq=idct(idct(ydq')');
xe=xdq+128;
xf=uint8(xe);
xza(i*8-7:i*8,j*8-7:j*8)=xf;
end
end
figure(3);imagesc(xza);colormap(gray)
title('p=4')
for i=1:75
for j=1:75
xz=xgray(i*8-7:i*8,j*8-7:j*8);
xd=double(xz);
xc=xd-128;
y=dct(dct(xc')');
p=10;
q=p*8./hilb(8);
yq=round(y./q);
ydq=yq.*q;
xdq=idct(idct(ydq')');
xe=xdq+128;
xf=uint8(xe);
xza(i*8-7:i*8,j*8-7:j*8)=xf;
end
end
figure(4);imagesc(xza);colormap(gray)
title('p=10')