I=imread("lena.bmp");% 提取影象

BW1=edge(I,"sobel"); %用SOBEL運算元進行邊緣檢測

BW2=edge(I,"roberts");%用Roberts運算元進行邊緣檢測

BW3=edge(I,"prewitt"); %用prewitt運算元進行邊緣檢測

BW4=edge(I,"log"); %用log運算元進行邊緣檢測

BW5=edge(I,"canny"); %用canny運算元進行邊緣檢測

h=fspecial("gaussian’,5);

BW6=edge(I,’canny’);

subplot(2,3,1), imshow(BW1);

title(‘sobel edge check’);

subplot(2,3,2), imshow(BW2);

title(‘sobel edge check’);

subplot(2,3,3), imshow(BW3);

title(‘prewitt edge check’);

subplot(2,3,4), imshow(BW4);

title(‘log edge check’);

subplot(2,3,5), imshow(BW5);

title(‘canny edge check’);

subplot(2,3,6), imshow(BW6);

title(‘gasussian&canny edge check’);%此為用高斯濾波後Canny運算元邊緣檢測結果

（注意：程式碼中有一些標點是中文模式，若輸入程式碼後標點顯示紅色，則為中文標點，改回來就行了）