load samplematrix32

K=zeros(N*N,N*N);

for i1=1:1:N  
    for i2=1:1:N
        for j1=1:1:N
            for j2=1:1:N
                %K((i1-1)*N+j1,(i2-1)*N+j2)=(exp(-((i1-i2)^2+(j1-j2)^2)/2.646^2)+0.2)/1.2;                 
                K((i1-1)*N+j1,(i2-1)*N+j2)=0.98^(abs(i1-i2))*0.98^(abs(j1-j2));                
            end
        end
    end
end

H=inv(SM*K*SM');

H=SM'*H;

H=K*H;

save Kernel H M N;