关于MATLAB图像相关性分析的问题

 对于下边的程序，它的功能是对读入的图像随机的选取1000对像素点，然后对像素点进行相关性分析：
clear all;
a=imread('lenac.bmp');
for k=1:1000
    i=fix(rand*250);
    j=fix(rand*250);
    x(k)=a(i+1,j+1);
    y(k)=a(i+2,j+1);
end
sum=0.0;
for i=1:1000
    sum=sum+x(i);
end
ex=sum/1000;
sum=0.0;
for i=1:1000
    sum=sum+y(i);
end
ey=sum/1000;
sum=0.0;
for i=1:1000
    sum=sum+x(i)*x(i);
end
ex2=sum/1000;
sum=0.0;
for i=1:1000
    sum=sum+y(i)*y(i);
end
ey2=sum/1000;
dx=ex2-ex*ex;
dy=ey2-ey*ey;
sum=0.0;
for i=1:1000
    sum=sum+x(i)*y(i);
end
exy=sum/1000;
cov=exy-ex*ey;
xg=cov/(sqrt(dx)*sqrt(dy))

这里有两个错误：首先系统会提示：sqrt使用错误，原因是dx和dy是unit8类型的。还有一个就是计算错误，你会发现不管你怎么加sum后边的值都一直会是255。
上边两处错误的原因其实都是出在unit8类型上，默认情况下，系统读入图像后储存的类型就是unit8。此时我们直接对储存图像的数据进行操作，则后边程序中出现的所有的变量就都是unit8类型的，由于unit8类型的最大值就是255，所以sum会一直停留在255不变，不管你加多少次。
要改正这种情况有以下两