如图公式,我现在得到了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')
|
|