[chinafjzjh]

clear
close all
I = imread('F:\a.jpg'); %注意输入的图像
subplot(2,2,1);
imshow(I);
title('原图像');
[m,n]=size(I);
J(m,n)=0;
J=double(J);
X=imnoise(I,'gaussian');Y=double(X); %加入噪声
subplot(2,2,2);
imshow(X);
title('加噪图像1');
J=J+Y/10;
X=imnoise(I,'gaussian');Y=double(X);
subplot(2,2,3);
imshow(X);
title('加噪图像2');
J=J+Y/10;
for i=1:14 %再循环14次
X=imnoise(I,'gaussian');
Y=double(X);
J=J+Y/10;
end
subplot(2,2,4);
imshow(mat2gray(J));
title('16幅噪声图像的平均结果');