编写一个利用巴特沃斯低通滤波器进行图像去噪得的Matlab程序。记得先进行加躁

在Matlab中，我们可以使用巴特沃斯低通滤波器对图像进行去噪处理。首先，你需要加载图像并添加噪声，然后应用滤波器去除噪声。以下是基本步骤：

% 1. 加载图像
img = imread('your_image.jpg'); % 替换为你要处理的图像文件名

% 2. 添加噪声 (这里以高斯噪声为例)
noiseLevel = 0.05; % 噪声水平
noisy_img = img + noiseLevel * randn(size(img));

% 3. 定义滤波器参数 (例如，巴特沃斯滤波器通常需要截止频率和阶数)
cutoff_freq = 0.1; % 截止频率，可以根据需求调整
filter_order = 4; % 滤波器阶数

% 4. 创建滤波器
[num, den] = buttord(cutoff_freq, 0.95, 'low', filter_order); % 巴特沃斯设计函数
b = fir1(num, den, 'type', 'low'); % 构建巴特沃斯滤波器系数

% 5. 应用滤波器
filtered_img = imfilter(noisy_img, b, 'replicate'); % 使用imfilter函数，'replicate'表示边缘复制边界处理

% 6. 显示原始图像、噪声图像和去噪后的图像
subplot(1, 3, 1), imshow(img);
title('Original Image');
subplot(1, 3, 2), imshow(noisy_img);
title('Noisy Image');
subplot(1, 3, 3), imshow(filtered_img);
title('Filtered Image');

%