MATLAB 图像处理：添加高斯噪声并绘制时域频域图

以下是通过 MATLAB 绘制加入高斯噪声图像的时域及频域图的示例代码：

% 读取原始图像
I = imread('lena.png');

% 添加高斯噪声
sigma = 10; % 噪声标准差
noise = sigma * randn(size(I)); % 高斯噪声
J = im2double(I) + noise; % 添加噪声

% 绘制时域图像
figure;
subplot(1,2,1); imshow(I); title('原始图像');
subplot(1,2,2); imshow(J); title('加噪图像');

% 绘制频域图像
F1 = fft2(I); % 原始图像的傅里叶变换
F2 = fft2(J); % 加噪图像的傅里叶变换
F1 = fftshift(F1); % 将零频移到中心
F2 = fftshift(F2);
S1 = log(1 + abs(F1)); % 原始图像的频谱
S2 = log(1 + abs(F2)); % 加噪图像的频谱
figure;
subplot(1,2,1); imshow(S1, []); title('原始图像频谱');
subplot(1,2,2); imshow(S2, []); title('加噪图像频谱');

运行上述代码，将会显示原始图像、加噪图像的时域图像和傅里叶变换后的频域图像，如下所示：