跪求基于matlab对sem图像进行边缘检测，各个算子程序怎么写？

im=imread('h1.jpg')

I=rgb2gray(im)

bw1=edge(I,'sobel')

bw2=edge(I,'roberts')

bw3=edge(I,'canny')

bw4=edge(I,'prewitt')

figure(2),subplot(2,2,1)imshow(bw4)title('prewitt算子效果图')

subplot(2,2,2)imshow(bw1)title('sobel算子效果图')

subplot(2,2,3)imshow(bw2)title('roberts算子效果图')

subplot(2,2,4)imshow(bw3)title('canny算子效果图')

具体如下：

1、第一步，打开matlab软件，出现如下界面，见下图，转到下面的步骤。

2、第二步，完成上述步骤后，敲入命令“clear；clc ”来清理工作空间，见下图，转到下面的步骤。

3、第三步，完成上述步骤后，敲入命令“syms   x”来定义一个符号变量，见下图，转到下面的步骤。

4、第四步，完成上述步骤后，敲入命令“f1 =  atan（x）”来定义一个函数，见下图，转到下面的步骤。

5、第五步，完成上述步骤后，敲入命令“df1  = diff（f1，x）”来求解导数函数的符号解，见下图，转到下面的步骤。

6、第六步，完成上述步骤后，敲入以下指令以查看处理结果，见下图，转到下面的步骤。

7、第七步，完成上述步骤后，查看派生结果，见下图。这样，就解决了这个问题了。

comp=[ ]

定义了空矩阵，

comp=[comp sem]

这句是在矩阵comp的后面，每次循环就加上一个sem，你可以理解吗？

我的腾讯联系方式，散坝而要龄要伞流午

---