Guassian Kernel

 

m = 5; n = 5;
sigma = 1;
[h1, h2] = meshgrid(-(m-1)/2:(m-1)/2, -(n-1)/2:(n-1)/2);
hg = exp(- (h1.^2+h2.^2) / (2*sigma^2));
hg = hg / (2 * pi * (sigma ^2));
h = hg ./ sum(hg(:));