如图公式，我现在得到了z的一维矩阵，它的顺序是z（i，j）,z（i+m，j）z(i+2m,j); 请问如何循环实现如图公式中的计算，并最终得到m和h的关系曲线。小弟初学者，以下是小弟的代码和数据文件，不知道问题出在哪里，计算结果和文献中的图不一样，求各位老师学长帮忙看看~ 感激不尽！

clear;
format long g
z=load('d:\20191106\a\a1.txt');
n=length(z);
h=0;
i=0.0195314;
m=100
for m=1:512
  for j=1:n
     for i=1:n-m
     h=h+(z(i+m,j)*-z(i,j))^2;
      end
  end
  h=1/(n*(n-m))*h;
  rr(m)=log(m*i);
  hh(m)=log(h);
end
  plot(rr,hh,'bo')

如图公式，我现在得到了z的一维矩阵，它的顺序是z（i，j）,z（i+m，j）z(i+2m,j); 请问如何循环实现如图公式中的计算，并最终得到m和h的关系曲线。小弟初学者，以下是小弟的代码和数据文件，不知道问题出在哪里，计算结果和文献中的图不一样，求各位老师学长帮忙看看~ 感激不尽！

clear;
format long g
z=load('d:\20191106\a\a1.txt');
n=length(z);
h=0;
i=0.0195314;
m=100
for m=1:512
  for j=1:n
     for i=1:n-m
     h=h+(z(i+m,j)*-z(i,j))^2;
      end
  end
  h=1/(n*(n-m))*h;
  rr(m)=log(m*i);
  hh(m)=log(h);
end
  plot(rr,hh,'bo')