d9r3.frm

常用的数值算法的VB程序

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   8220
   ClientLeft      =   1665
   ClientTop       =   540
   ClientWidth     =   7650
   LinkTopic       =   "Form1"
   ScaleHeight     =   8220
   ScaleWidth      =   7650
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   6120
      TabIndex        =   0
      Top             =   120
      Width           =   1215
   End
End
Attribute VB_Name = "Form1"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Private Sub Command1_Click()
    'PROGRAM D9R3
    'Driver for routine COVSRT
    MA = 10
    MFIT = 5
    Dim COVAR(10, 10), LISTA(5)
    For I = 1 To MA
        For J = 1 To MA
            COVAR(I, J) = 0#
            If I <= 5 And J <= 5 Then
                COVAR(I, J) = I + J - 1
            End If
        Next J
    Next I
    Print Tab(5); "Original matrix"
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
    'Test 1 - spread by 2
    Print Tab(5); "Test #1 - Spread by two"
    For I = 1 To MFIT
        LISTA(I) = 2 * I
    Next I
    Call COVSRT(COVAR(), MA, MA, LISTA(), MFIT)
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
    'Test 2 - reverse
    Print Tab(5); "Test #2 - Reverse"
    For I = 1 To MA
        For J = 1 To MA
            COVAR(I, J) = 0#
            If I <= 5 And J <= 5 Then
                COVAR(I, J) = I + J - 1
            End If
        Next J
    Next I
    For I = 1 To MFIT
        LISTA(I) = MFIT + 1 - I
    Next I
    Call COVSRT(COVAR(), MA, MA, LISTA(), MFIT)
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
    'Test 3 - spread and reverse
    Print Tab(5); "Test #3 - Spread and reverse"
    For I = 1 To MA
        For J = 1 To MA
            COVAR(I, J) = 0#
            If I <= 5 And J <= 5 Then
                COVAR(I, J) = I + J - 1
            End If
        Next J
    Next I
    For I = 1 To MFIT
        LISTA(I) = MA + 2 - 2 * I
    Next I
    Call COVSRT(COVAR(), MA, MA, LISTA(), MFIT)
    For I = 1 To MA
        For J = 1 To MA
            Print Tab(J * 6); Format$(COVAR(I, J), "#.0");
        Next J
    Next I
End Sub
Sub COVSRT(COVAR(), NCVM, MA, LISTA(), MFIT)
    For J = 1 To MA - 1
        For I = J + 1 To MA
            COVAR(I, J) = 0#
        Next I
    Next J
    For I = 1 To MFIT - 1
        For J = I + 1 To MFIT
            If LISTA(J) > LISTA(I) Then
                COVAR(LISTA(J), LISTA(I)) = COVAR(I, J)
            Else
                COVAR(LISTA(I), LISTA(J)) = COVAR(I, J)
            End If
        Next J
    Next I
    SWAP = COVAR(1, 1)
    For J = 1 To MA
        COVAR(1, J) = COVAR(J, J)
        COVAR(J, J) = 0#
    Next J
    COVAR(LISTA(1), LISTA(1)) = SWAP
    For J = 2 To MFIT
        COVAR(LISTA(J), LISTA(J)) = COVAR(1, J)
    Next J
    For J = 2 To MA
        For I = 1 To J - 1
            COVAR(I, J) = COVAR(J, I)
        Next I
    Next J
End Sub