kstwo.txt

来自「《Delphi常用数值算法集》的配书源码」· 文本 代码 · 共 45 行

procedure KSTWO(DATA1:array of real; N1:integer;
                    DATA2:array of real; N2:integer;var D, PROB:real);
label 1;
var
    J1,J2:integer;  FO1,FO2,FN1,FN2,EN1,EN2,AM,DT:real;
begin
    SORT(N1, DATA1);
    SORT(N2, DATA2);
    EN1:=N1;
    EN2:=N2;
    J1:=1;
    J2:=1;
    FO1:=0;
    FO2:=0;
    D:=0; 
1:  If (J1 <= N1) And (J2 <= N2) Then
    begin
        If DATA1[J1] < DATA2[J2] Then
        begin
            FN1:=J1 / EN1;
            If Abs(FN1 - FO2) > Abs(FO1 - FO2) Then
                AM:=Abs(FN1 - FO2)
            Else
                AM:=Abs(FO1 - FO2);
            DT:=AM;
            If DT > D Then D:=DT;
            FO1:=FN1;
            J1:=J1 + 1;
        end
        Else
        begin
            FN2:=J2 / EN2;
            If Abs(FN2 - FO1) > Abs(FO2 - FO1) Then
                AM:=Abs(FN2 - FO1)
            Else
                AM:=Abs(FO2 - FO1);
            DT:=AM;
            If DT > D Then D:=DT;
            FO2:=FN2;
            J2:=J2 + 1;
        end;
        GoTo 1;
    end;
    PROB:=PROBKS(Sqrt(EN1 * EN2 / (EN1 + EN2)) * D);
end;