unit1.pas

可以计算出灰色模型的特征向量,在delphi5.0上已经做过,可以执行

procedure  Tyuce.hsxtClick(Sender:TObject);
type  shuzu=array[7.5,7.53,7.56,7.53,7.54,7.57,7.66,7.82,7.90] of real;
var
      sz,zjsj,ymsz:shuzu;
      i2:integer;d1,d2,d3,d4,v:varint; s,d:string;
begin
     i1:=4;
begin
       sz[1,1]:=strtoint(xiaos1.text);
       sz[1,2]:=strtoint(xiaos2.text);
       sz[1,3]:=strtoint(xiaos3.text);
       sz[1,4]:=strtoint(xiaos4.text);
       d1:=shuj1.text; d3:=d1;sz[2,1]:=d1;
       d1:=shuj2.text; sz[2,2]:=d1;
       d1:=shuj3.text; sz[2,3]:=d1;
       d1:=shuj4.text; sz[2,4]:=d1;
       for i2:=2 to i1 do ymsz[1,i2-1]:=sz[2,i2];
       d1:=0; d2:=0;
       for i2:=2 to i1 do
sz[2,i2]:=sz[2,i2]+sz[2,i2-1];d4:=sz[2,i1];
       for i2:=2 to i1 do
zjsz[1,i2-1]:=(sz[2,i2]+sz[2,i2-1]*(-0.5);
       for i2:=2 to i1 do zjsz[2,i2-1]:=1;
       sz[1,1]:=0; sz[[1,2]:=0;
for i2:=1 to (i1-1) do begin
sz[1,1]:=sz[1,1]+sqr(zjsz[1,i2]);sz[1,2]:=sz[1,2]+zjsz[1,i2];
end;
sz[2,1]:=sz[1,2];sz[2,2]:=i1-1;
       d1:=sz[1,1]*sz[2,2]-sz[1,2]*sz[2,1];
       d2:=sz[1,1];
       sz[1,1]:=sz[2,2]/d1;
       sz[1,2]:=sz[1,2]*(-1)/d1;
       sz[2,1]:=sz[2,1]*(-1)/d1;
       sz[2,2]:=d2/d1;
       d1:=0; d2:=0;
       for i2:=1 to (i1-1) do begin
d1:=d1+zjsz[1,i2]*ymsz[1,i2];
d2:=d2+ymsz[1,i2];  end;
       zjsz[1,1]:=d1; zjsz[1,2]:=d2;
d1:=sz[1,1]*zjsz[1,1]+sz[1,2]*zjsz[1,2];
d2:=sz[2,1]*zjsz[1,1]+sz[2,2]*zjsz[1,2];
       d1:=(d3-d2/d1)*exp(d1*i1*(-1))+d2/d1;
       v:=d1-d4;
       if v<1 then begin d:=v; s:=’0’+d; edit11.text:=s;
       end
       else edit11.text:=v;
       end;
       end;