%Script per la sperimentazione sulla sfocatura di immagini
A=imread('gatto.png'); % Carico l’immagine
A=double(A);
A=A/255;
[m,n,~]=size(A);
[B,psf]=sfoca(A); % Applico la sfocatura
B = B + randn(m,n, 3) * 0.02; % Aggiungo rumore
C=recover(B, psf); % Applico il metodo del gradiente coniugato
D=recover_reg(B,psf); % Applico il metodo del gradiente coniugato con regolarizzione
imwrite(B,'gatto_sfocato.png');
imwrite(C,'gatto_rec.png');
imwrite(D,'gatto_rec_reg.png');