Pelembutan Citra (MATLAB)

function pelembutan_citra
%menggunakan filter lolos rendah
%syarat: semua koefisien filter harus positif
%syarat: jumlah semua koefisien harus sama dengan 1

A=imread('cameraman.tif');
m=[1/9 1/9 1/9;1/9 1/9 1/9;1/9 1/9 1/9];
%m=[1/10 1/10 1/10; 1/10 1/5 1/10;1/10 1/10 1/10];
%m=[1/16 1/8 1/16; 1/8 1/4 1/8; 1/16 1/8 1/16];
d=A
[baris kolom]=size(A);
for i=2:baris-1;
    for j=2:kolom-1;
        d(i,j)=A(i-1,j-1)*m(1,1)+A(i-1,j)*m(1,2)+A(i-1,j+1)*m(1,3)+A(i,j-1)*m(2,1)+A(i,j)*m(2,2)+A(i,j+1)*m(2,3)+A(i+1,j-1)*m(3,1)+A(i+1,j)*m(3,2)+A(i+1,j+1)*m(3,3);
        if d(i,j)>255
            d(i,j)=255;
        else if d(i,j)<0
                d(i,j)=0;
            end
        end
    end
end
subplot(2,2,1), imshow(A);
subplot(2,2,2), imhist(A);
subplot(2,2,3), imshow(d);
subplot(2,2,4), imhist(d);