修改一下代码使其能运行：function [outputImg] = saltpepper(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

function [outputImg] = saltpepper(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));
% 计算当前像素的输出值
if denominator == 0
outputImg(i,j) = inputImg(i,j);
else
outputImg(i,j) = numerator / denominator;
end
end
end

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