编写matlab程序，实现下图所示的图像处理过程:(a) Pa=Pb=0.1的脉冲噪声；(b)3×3的中值滤波器；(c) 第二次中值滤波器处理；(d)第三次中值滤波器处理

在MATLAB中实现上述图像处理过程，你需要分步骤进行：

1. **添加脉冲噪声**:
使用`imnoise`函数可以添加椒盐噪声。假设你的原始图像为`img`，添加0.1的脉冲噪声的代码如下：

   Pa = 0.1; % 高斯噪声比例
   Pb = 0.1; % 脉冲噪声比例
   noisy_img = imnoise(img, 'salt & pepper', Pa, Pb);

2. **3x3中值滤波**:
中值滤波器可以使用`medfilt2`函数，将图像传递给它，并指定过滤器大小为3x3：

   med_filt_size = 3;
   filtered_img_1 = medfilt2(noisy_img, [med_filt_size, med_filt_size]);

3. **多次中值滤波**:
你可以重复这个步骤，对第一次过滤后的图像再次进行中值滤波：

   filtered_img_2 = medfilt2(filtered_img_1, [med_filt_size, med_filt_size]);
   filtered_img_3 = medfilt2(filtered_img_2, [med_filt_size, med_filt_size]);

4. **显示各步骤结果**:
对每个处理后的图像保存并查看，可以使用`imshow`函数：

   figure;
   subplot(2, 2, 1), imshow(img, 'InitialMagnification', 'fit');
   title('Original Image');
   subplot(2, 2, 2), imshow(noisy_img, 'InitialMagnification', 'fit');
   title(['Noisy Image (Pa=' num2str(Pa) ', Pb=' num2str(Pb) ')']);
   subplot(2, 2, 3), imshow(filtered_img_1, 'InitialMagnification', 'fit');
   title('First Median Filter');
   subplot(2, 2, 4), imshow(filtered_img_3, 'InitialMagnification', 'fit');
   title('Final Result after Third Filter');