如何自己写medfilt2函数

medfilt2函数是MATLAB中用于进行二维中值滤波的函数，可以通过以下步骤自己编写：

1. 定义输入和输出参数。输入参数为待滤波的二维矩阵A和滤波器的大小K，输出参数为滤波后的二维矩阵B。

2. 根据滤波器大小K计算出滤波器的中心点位置。

3. 对于输入矩阵A中的每个像素，以其为中心位置，根据滤波器大小K确定其邻域内的像素值，并将这些像素值存储到一个一维数组中。

4. 将存储像素值的一维数组进行排序，然后取其中位数作为滤波器输出值。

5. 将所有像素的滤波器输出值按照原始像素的位置放入输出矩阵B中。

以下是一个示例代码，实现了一个3x3的中值滤波器：

function B = medfilt2(A,K)
% A: 输入矩阵
% K: 滤波器大小
% B: 输出矩阵

% 计算滤波器中心位置
center = floor(K/2);

% 通过循环对每个像素进行中值滤波
for i = 1:size(A,1)
    for j = 1:size(A,2)
        % 确定滤波器邻域内的像素值
        row_min = max(i-center,1);
        row_max = min(i+center,size(A,1));
        col_min = max(j-center,1);
        col_max = min(j+center,size(A,2));
        neighborhood = A(row_min:row_max,col_min:col_max);
        
        % 对邻域内的像素值进行排序，并取中位数作为输出值
        sorted = sort(neighborhood(:));
        B(i,j) = sorted(floor(numel(sorted)/2)+1);
    end
end

需要注意的是，这个示例代码只实现了一个简单的中值滤波器，如果需要处理更复杂的情况，可能需要进行更多的优化和调整。