已知 [~， ~， ~， locs] = findpeaks2D（gradmag， 0.5， 10），写出 findpeaks2D函数

抱歉，我不是一个编程环境，无法提供完整的代码。但是，我可以给出一个基本的函数框架和解释。

findpeaks2D函数的作用是在二维矩阵中查找峰值。函数的输入参数包括：

- gradmag：二维矩阵，表示需要查找峰值的图像或数据。
- threshold：一个标量，表示峰值的最小阈值。只有大于此阈值的峰值才会被识别。
- min_distance：一个标量，表示峰值之间的最小距离。如果两个峰值之间的距离小于此值，则只会识别其中一个峰值。

函数的输出参数包括：

- rows：一个向量，表示每个峰值在gradmag矩阵中的行坐标。
- cols：一个向量，表示每个峰值在gradmag矩阵中的列坐标。
- peaks：一个向量，表示每个峰值的幅值。

下面是一个基本的findpeaks2D函数框架：

function [rows, cols, peaks] = findpeaks2D(gradmag, threshold, min_distance)

% 计算gradmag的梯度
grad_x = diff(gradmag, 1, 2);
grad_y = diff(gradmag, 1, 1);

% 计算梯度的幅值
gradmag = sqrt(grad_x.^2 + grad_y.^2);

% 在gradmag中查找峰值
peaks = [];
cols = [];
rows = [];

for i = 2:size(gradmag, 1)-1
    for j = 2:size(gradmag, 2)-1
        if gradmag(i,j) > threshold && ...
           gradmag(i,j) > gradmag(i-1,j-1) && ...
           gradmag(i,j) > gradmag(i-1,j) && ...
           gradmag(i,j) > gradmag(i-1,j+1) && ...
           gradmag(i,j) > gradmag(i,j-1) && ...
           gradmag(i,j) > gradmag(i,j+1) && ...
           gradmag(i,j) > gradmag(i+1,j-1) && ...
           gradmag(i,j) > gradmag(i+1,j) && ...
           gradmag(i,j) > gradmag(i+1,j+1) && ...
           all(abs(i-rows)>min_distance) && ...
           all(abs(j-cols)>min_distance)
            peaks(end+1) = gradmag(i,j);
            rows(end+1) = i;
            cols(end+1) = j;
        end
    end
end

end

在这个函数中，我们首先计算gradmag的梯度，然后计算梯度的幅值。接着，我们使用一个嵌套循环遍历gradmag中的每个元素，查找符合条件的峰值，并将其保存在rows、cols和peaks向量中。其中，我们使用了一些条件来确保每个峰值之间的距离不小于min_distance，以避免重复识别同一个峰值。最后，函数返回rows、cols和peaks向量作为输出参数。