%Buzási Bence
%TXEVVQ

function [kernel_inverse_matrix] = TXEVVQ_kernel_inverse_creator()

n1=8;
n2=8;
u = 1:8;
v = 1:8;

[u1,v1] = meshgrid(u,v);
n = @(x) ((x~=1)*1) + ((x==1)*(1/sqrt(2)));
kernel = @(x,y) 1/4 .* n(u1) .* n(v1) .* cos( ((2.*(x-1)+1).*(u1-1).*pi) ./ ...
(2* n1)) .* cos( ((2.*(y-1)+1).*(v1-1).*pi) ./ (2* n2));

figure(3)
kernel_inverse_matrix = zeros(8,8,64);
for x = 1 : n1
    for y = 1 : n2
       kernel_inverse_matrix(:,:,(x-1)*8+y) = kernel(x,y); 
       subplot(8,8,(x-1)*8+y)
       imshow(uint8(kernel_inverse_matrix(:,:,(x-1)*8+y)*128+128))
    end
end

end