delphi - 计算赫斯特指数

Question

你好 =) 提前抱歉我的英语不好我的任务是通过线性回归方法计算赫斯特指数。我有解决方案的文字描述。它看起来很简单，但我总是得到超出范围 0..1 的值。通常，值为 1.9 或类似值。有时它会得到接近于零的负值。我已经检查了代码大约一千次，但没有发现错误。

var
max_z,min_z,x_m:real; //max and min of cumulative sum and mean value of X for every Tau
st,ss,sst,st2 :real;
Al, Herst: real;
x_vr:array of double;   //a piece of array with length=tau
i, j, nach: integer;
begin
    //file opening and getting values of X array are in another function
    nach:=3;    //initial value of tau
    Setlength(ln_rs,l-nach); //length of  ln(R/S) array
    Setlength(ln_t,l-nach);  //length of  ln(tau) array
    Setlength(r,l-nach);   //length of  R array
    Setlength(s,l-nach);   //length of S array


      //Let's start
    for tau:=nach to l do  //we will change tau
    begin
        Setlength(x_vr,tau+1); //set new local array (length=tau)
        for i:=0 to length(x_vr)-1 do
            x_vr[i]:=x[i];

        x_m:=Mean(x_vr);    //mean value
        Setlength(y,tau+1);   //length of array of difference from mean value
        Setlength(z,tau+1);   //length of array of cumulative sum

        for i:=0 to tau do
            y[i]:=x_vr[i]-x_m;      //difference from mean value

        z[0]:=y[0];
        for i:=1 to tau do      //cumulative sum
             for j :=i downto 0 do
                z[i]:=z[i]+y[j];

        max_z:=z[0];
        for i:=1 to tau do        //max of cumulative sum
            max_z:=max(max_z,z[i]);

        min_z:=z[0];
        for i:=1 to tau do        //min of cumulative sum
            min_z:=min(min_z,z[i]);

        r[tau-nach]:=max_z-min_z;    //R value
        s[tau-nach]:=0;
        for i:=0 to tau do
            s[tau-nach]:=power(y[i],2)+s[tau-nach];         //S value

        s[tau-nach]:=sqrt(s[tau-nach]/(tau+1));

        //new array values
        ln_rs[tau-nach]:=Ln(R[tau-nach]/S[tau-nach]);   // ln(R/S)
        ln_t[tau-nach]:=ln(tau);                        // ln (tau)

    end;    //End of calculating

    //Method of Least squares
    for i:=0 to length(ln_rs)-1 do  
        st:=st+ln_t[i];

    st:=(1/length(ln_rs))*st;

    for i:=0 to length(ln_rs)-1 do
        ss:=ss+ln_rs[i];

    ss:=(1/length(ln_rs))*ss;

    for i:=0 to length(ln_rs)-1 do
        sst:=sst+ln_t[i]*ln_rs[i];

    sst:=(1/length(ln_rs))*sst;

    for i:=0 to length(ln_rs)-1 do
        st2:=st2+ln_t[i]*ln_t[i];

    st2:=(1/length(ln_rs))*st2;


    Herst:=(sst-st*ss)/(st2-st*st);      //coefficient of approximal function
    al:=ss-st*Herst;

谢谢大家=)

附注

 for tau:=nach to l do

有 L，而不是 1。L 是 X 数组的长度。当 l=nach 时，L>nach 总是除了最后一步之外。

P.P.S。它有效，伙计们。但值(value)观不对。他们走出了射程。也许，算法有错误。或者也许我跳过了某些步骤。

最后更新

这很神秘，但我只改变了计算数组 Z 的方法，它开始正常工作......谢谢大家 =)

Answer 1

我看到的第一件事:

nach := 3;    
for tau := nach to l do  //w

这算数了。并且因为 nach>1，所以不会执行该循环体。

如果您希望倒计时。使用 downto 变体。倒计时:

for tau := nach downto l do  //w