value = [];
for w = 0:5:50
x = 1:50;
y = roundn(0.1*x+2+w*rand(1,size(x,2)),-2);
plot(x,y,'bo');
value = [value sum((x(1,:)-mean(x(1,:))).*(y(1,:)-mean(y(1,:))))/(size(x,2))];
value = roundn(value,-2);
end
value = [];
coef = [];
for w = 0:5:50
x = 1:50;
y = x+w*rand(1,size(x,2));
plot(x,y,'bo');
coef = [coef sum((x(1,:)-mean(x(1,:))).*(y(1,:)-mean(y(1,:))))/(std(x)*std(y)*(size(x,2)))];
coef = roundn(coef,-2);
end
参考地址:https://blog.youkuaiyun.com/xiwan0902/article/details/72280352/