用matlab写出用逆谐波均值滤波消除图片椒盐噪声的代码并且每行写上注释

以下是用MATLAB实现逆谐波均值滤波器去除图像椒盐噪声的代码，每行都有注释：

function [outputImg] = nhmean_filter(inputImg, order, Q)
% inputImg: 输入图像
% order: 滤波器阶数
% Q: 逆谐波均值滤波器的 Q 值

% 获取输入图像的大小
[row, col] = size(inputImg);

% 初始化输出图像
outputImg = zeros(row, col);

% 对于每个像素，用逆谐波均值滤波器处理
for i = 1:row
    for j = 1:col
        % 计算滤波器对于当前像素的局部窗口
        window = inputImg(max(i-order,1):min(i+order,row), max(j-order,1):min(j+order,col));
        % 计算窗口中非零像素的 Q 次方和
        numerator = sum(window.^Q);
        % 计算窗口中非零像素的 Q 次方和的倒数
        denominator = sum(window.^(Q+1));
        % 计算当前像素的输出值
        outputImg(i,j) = numerator / denominator;
    end
end

% 将输出图像转换为无符号 8 位整数类型
outputImg = uint8(outputImg);
end

使用方法：

% 读取图像
inputImg = imread('lena_salt_noise.png');

% 调用逆谐波均值滤波器函数
outputImg = nhmean_filter(inputImg, 3, 1.5);

% 显示输入图像和输出图像
figure;
subplot(1,2,1);
imshow(inputImg);
title('输入图像');
subplot(1,2,2);
imshow(outputImg);
title('输出图像');