gm(1,1).m

来自「many matlab script ,i colect many things」· M 代码 · 共 29 行

q=0.6;
n=4;
x0=[1120,1150,1450,1710];
g1=exp(-2/(n+1));g2=exp(2/(n+1));g3=0;
for i=2:n
    g4=x0(i-1)/x0(i);
    if (g4<g1)|(g4>g2)
        g3=1+g3;
    end
end
g3
x1(1)=x0(1);c=0;d=0;e=0;f=0;
for i=2:n
    x1(i)=x0(i)+x1(i-1);
    z1(i)=q*x1(i)+(1-q)*x1(i-1);
    z2(i)=z1(i)*z1(i);
    c=z1(i)+c;
    d=x0(i)+d;
    e=z1(i)*x0(i)+e;
    f=z2(i)+f;
end
a=(c*d-(n-1)*e)/((n-1)*f-c*c);
b=(d*f-c*e)/((n-1)*f-c*c);
g=b/a;
m=-1*n*a;
y=(x0(1)-g)*exp(m)+g;
x=y-x1(n);
x