利用matlab实现图像的角点检测

function [y1y2rc]=harris（image）
image=rgb2gray（image）;
image=double（image）/255;
[heightwidth]=size（image）;
X=[-1 0 1; -1 0 1; -1 0 1]; %X方向梯度算子
Y=X‘; %Y方向梯度算子
fx=conv2（imageX‘same‘）; %X方向滤波
fy=conv2（imageY‘same‘）;
fx2=fx.^2;
fy2=fy.^2;
fxy=fx.*fy;
clear fx;
clear fy;
g=fspecial（‘gaussian‘[55]2）;
A=conv2（fx2g‘same‘）;
B=conv2（fy2g‘same‘）;
C=conv2（fxyg‘same‘）;
R=zeros（heightwidth）;
Result=zeros（heightwidth）;
Rmax=0.0;
for i=1:height
    for j=1:width
        M=[A（ij） C（ij）
           C（ij） B（ij）];
       R（i