%裝載原始圖像
load sinsin;
subplot(2,2,1);image(X);colormap(map);title('含噪圖像');
%對圖像加噪聲并顯示出含噪圖像
init=2788605800;
randn('seed',init);
X=X+10*randn(size(X));
subplot(2,2,2);image(X);colormap(map);title('含噪圖像');
%應用中值濾波進行圖像平滑處理
[p,q]=size(X);
for i=2:p-1
    for j=2:q-1
        Xtemp=0;
        for m=1:3
            for n=1:3
                Xtemp=Xtemp+X(i+m-2,j+n-2);
            end
        end
        Xtemp=Xtemp/9;
        X1(i,j)=Xtemp;
    end
end
%顯示結果
subplot(2,2,3);image(X1);colormap(map);
title('平滑后圖像');