MATLAB叠加高斯白噪声

clear;
I = imread('.\img\Lena.tif');  % 读取一个灰度图或彩色图
I = im2double(I);   % 将uint8转化成double并归一化

% 避免叠加噪声后像素值超过[0 255]
I = imadjust(I,[0 1],[0.35 0.65]);
I = I*255;

% 叠加高斯白噪声
[M,N] = size(I);
NOISE = 0+10*randn(M,N); % 均值为0，标准差为10（方差为100）
I1 = I + NOISE;

figure(1);
subplot(1,2,1);imshow(I,[0 255]);title('原图');
subplot(1,2,2);imshow(I1,[0 255]);title('叠加噪声后');

此外，可以用以下语句验证所生成噪声的均值和方差是否满足要求 

var_NOISE = var(NOISE(:))
mean_NOISE = mean(NOISE(:))