deter.de

微分进化法的pascal程序

Function Determination ( n_pop  : integer;tempui: popsizematrix ) : Double;
   VAR    suma1,suma2,suma3  : Double;
Begin
   suma1:=0.0;
   suma2:=0.0;
   suma3:=0.0;


   For i:=1 To XY Do
       ydata_obl[i]:=tempui[n_pop,1]+tempui[n_pop,2]*xdata[i]+tempui[n_pop,3]*xdata[i]*xdata[i];


   For i:=1 To XY Do
       suma2:= suma2+ ydata[i];
   suma2:=suma2/XY;

   For i:=1 To XY Do
   Begin
      suma1:= suma1+ ((ydata[i]-ydata_obl[i]) * (ydata[i]-ydata_obl[i]) );
      suma3:= suma3+ ((ydata[i]-suma2) * (ydata[i]-suma2) );
   End;

   Determination:= suma1/suma3;
End;