function [x h2] = BlurImageAddNoise (s, sigma_blur, sigma_noise)

h = zeros(512);
for x=1:512 
    for y=1:512 
        h(x,y) = exp(-((x-256)^2+(y-256)^2)/sigma_blur^2); 
    end;
end;
h = h./sum(h(:));
h2=circshift(h,[256 256]);
x = ifft2(fft2(s).*fft2(h2)) + sigma_noise*randn(512);