matlab边缘检测代码

A=imread（‘yanghuawumoli.png‘）; 
B=rgb2gray（A）; 
C=im2bw（B）; 
subplot（121）; 
imshow（C）; 
[mn]=size（C）; 
W=zeros（334）;   
W（::1）=[-1 -1 -1; 2 2 2; -1 -1 -1]; 
W（::2）=[-1 -1 2; -1 2 -1; 2 -1 -1]; 
W（::3）=[-1 2 -1; -1 2 -1; -1 2 -1]; 
W（::4）=[2 -1 -1; -1 2 -1; -1 -1 2]; 
D=zeros（mn）; 
R=zeros（41）; 
M=zeros（33）; 
for i=2:m-1  
    for j=2:n-1 
        if（C（ij））  
            M（::）=C（i-1:i+1 j-1:j+1）;