d13r6b.txt

常用的数值算法的VB程序

Private Sub Command1_Click()
    'PROGRAM D13R6B
    'Driver for routine TPTEST
    'Compare two correlated distributions vs. two
    'uncorrelated distributions
    NPTS = 500
    EPS = 0.01
    NSHFT = 10
    ANOISE = 0.3
    Dim DATA1(500), DATA2(500), DATA3(500)
    IDUM& = -5
    Print
    Print Tab(5); "         Correlated:             Uncorrelated:"
    Print Tab(5); "Shift     T      Probability      T     Probability"
    OFFSET = (NSHFT / 2) * EPS
    For J = 1 To NPTS
        GAUSS = GASDEV(IDUM&)
        DATA1(J) = GAUSS
        DATA2(J) = GAUSS + ANOISE * GASDEV(IDUM&)
        DATA3(J) = GASDEV(IDUM) + ANOISE * GASDEV(IDUM&)
    Next J
    Call AVEVAR(DATA1(), NPTS, AVE1, VAR1)
    Call AVEVAR(DATA2(), NPTS, AVE2, VAR2)
    Call AVEVAR(DATA3(), NPTS, AVE3, VAR3)
    For J = 1 To NPTS
        DATA1(J) = DATA1(J) - AVE1 + OFFSET
        DATA2(J) = DATA2(J) - AVE2
        DATA3(J) = DATA3(J) - AVE3
    Next J
    For I = 1 To NSHFT + 1
        SHIFT = I * EPS
        For J = 1 To NPTS
            DATA2(J) = DATA2(J) + EPS
            DATA3(J) = DATA3(J) + EPS
        Next J
        Call TPTEST(DATA1(), DATA2(), NPTS, T1, PROB1)
        Call TPTEST(DATA1(), DATA3(), NPTS, T2, PROB2)
        Print Tab(5); Format$(SHIFT, "#.#0");
        Print Tab(12); Format$(T1, "#.###0");
        Print Tab(24); Format$(PROB1, "#.###0");
        Print Tab(38); Format$(T2, "#.###0");
        Print Tab(52); Format$(PROB2, "#.###0")
    Next I
End Sub