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方向?yàn)V波
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