load kernel;

load goldhill-sample; III = imread('goldhill.bmp');

CC = imresize(I,2,'bicubic');

size=256;

II=zeros(size,size);

g=ones(M*M,1);

num1=128/M; num2=128/M;

for t1=1:1:num1
    for t2=1:1:num2
        for i=1:1:M
            for j=1:1:M
                g((i-1)*M+j) = I((t1-1)*M+i,(t2-1)*M+j);
            end
        end
        temp=H*g;
        for i=1:1:N
            for j=1:1:N
                II((t1-1)*N+i,(t2-1)*N+j)=temp((i-1)*N+j);
            end
        end
    end
end

figure;imshow(CC,[]);
figure;imshow(II,[]);
max=II(1,1); min=max;
max1=CC(1,1); min1=max1;
for i=1:1:size
    for j=1:1:size
        if II(i,j)>max
            max=II(i,j);
        end
        if II(i,j)<min
            min=II(i,j);
        end       
        if CC(i,j)>max1
            max1=CC(i,j);
        end
        if CC(i,j)<min1
            min1=CC(i,j);
        end       
    end
end
In=(II-min)./(max-min);
Ic=(CC-min1)./(max1-min1);
%imwrite(In,'LenaM.bmp','bmp');
%imwrite(Ic,'LenaC.bmp','bmp');

% Calculate the errors
temp=ones(size,size);
for i=1:1:size
    for j=1:1:size
        temp(i,j)=III(i,j);
    end
end
sum=0;sum1=0;
for i=1:1:size
    for j=1:1:size
        sum=sum+(II(i,j)-temp(i,j))^2;
        sum1=sum1+temp(i,j)^2;
    end
end
MSE=sum/sum1;
fprintf('Markov___MSE = %4.3e \n',MSE);

sum=0;
for i=1:1:size
    for j=1:1:size
        sum=sum+(CC(i,j)-temp(i,j))^2;
    end
end
MSE=sum/sum1;
fprintf('Cubic____MSE = %4.3e',MSE);