用matlab写出边缘检测的代码

边缘检测是图像处理中常用的一种技术，Matlab提供了多种边缘检测方法，比如Sobel、Prewitt、Canny等，下面是一个Sobel边缘检测的Matlab代码示例：

% 读取图像
img = imread('lena.jpg');
% 转为灰度图像
gray_img = rgb2gray(img);
% Sobel算子
Sobel_x = [-1 0 1; -2 0 2; -1 0 1];
Sobel_y = Sobel_x';
% 垂直方向和水平方向卷积
Gx = conv2(double(gray_img), Sobel_x, 'same');
Gy = conv2(double(gray_img), Sobel_y, 'same');
% 边缘梯度幅值
G = sqrt(Gx.^2 + Gy.^2);
% 边缘梯度方向
theta = atan2(Gy, Gx);
% 非极大值抑制
[nms_G, nms_theta] = non_max_suppression(G, theta);
% 双阈值检测
low_th = 0.1 * max(nms_G(:));
high_th = 0.3 * max(nms_G(:));
[edges, ~] = hysteresis_thresholding(nms_G, nms_theta, low_th, high_th);
% 显示结果
imshow(edges);

其中，`non_max_suppression`和`hysteresis_thresholding`是辅助函数，用于进行非极大值抑制和双阈值检测。完整代码可以参考以下链接：

https://github.com/zhixuan-lin/matlab-edge-detection/blob/master/edge_detection.m