d1r10.frm

矩阵特征值的求解过程之二

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   4800
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   4680
   LinkTopic       =   "Form1"
   ScaleHeight     =   4800
   ScaleWidth      =   4680
   StartUpPosition =   3  'Windows Default
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   2760
      TabIndex        =   0
      Top             =   4200
      Width           =   1695
   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 D1R10
    'Driver program for routine CHOBSB,CHODCM
    N = 5
    Dim A(5, 5), B(5), C(5, 5), D(5), T(5), X(5)
    '输入已知的方程组的系数矩阵
    A(1, 1) = 1.4: A(1, 2) = 2.1: A(1, 3) = 2.1: A(1, 4) = 7.4: A(1, 5) = 9.6
    A(2, 1) = 1.6: A(2, 2) = 1.5: A(2, 3) = 1.1: A(2, 4) = 0.7: A(2, 5) = 5#
    A(3, 1) = 3.8: A(3, 2) = 8#:  A(3, 3) = 9.6: A(3, 4) = 5.4: A(3, 5) = 8.8
    A(4, 1) = 4.6: A(4, 2) = 8.2: A(4, 3) = 8.4: A(4, 4) = 0.4: A(4, 5) = 8#
    A(5, 1) = 2.6: A(5, 2) = 2.9: A(5, 3) = 0.1: A(5, 4) = 9.9: A(5, 5) = 7.7
    '输入已知的方程组的右端向量B
    B(1) = 4#
    B(2) = 9.3
    B(3) = 8.4
    B(4) = 0.4
    B(5) = 4.1
    Print
    Print Tab(5); "已知的方程组的右端向量"
    Print Tab(14); Format$(B(1), "##.##")
    Print Tab(14); Format$(B(2), "##.##")
    Print Tab(14); Format$(B(3), "##.##")
    Print Tab(14); Format$(B(4), "##.##")
    Print Tab(14); Format$(B(5), "##.##")
    For I = 1 To N
        For J = 1 To N
           C(I, J) = A(I, J)
        Next J
    Next I
    Call CHODCM(C(), N, D(), T())
    For I = 1 To N
        X(I) = B(I)
    Next I
    Call CHOBSB(C(), N, D(), X())
    Print
    Print Tab(5); "计算出的方程组的解"
    Print Tab(10); Format$(X(1), "#.####E+00")
    Print Tab(10); Format$(X(2), "#.####E+00")
    Print Tab(10); Format$(X(3), "#.####E+00")
    Print Tab(10); Format$(X(4), "#.####E+00")
    Print Tab(10); Format$(X(5), "#.####E+00")
    '将计算出的解X乘以系数矩阵，以验证计算结果正确
    For L = 1 To N
        B(L) = 0#
        For J = 1 To L
            B(L) = B(L) + A(J, L) * X(J)
        Next J
        For J = L + 1 To N
            B(L) = B(L) + A(L, J) * X(J)
        Next J
    Next L
    Print
    Print Tab(5); "计算出的解乘以系数矩阵的结果"
    Print Tab(14); Format$(B(1), "##.##")
    Print Tab(14); Format$(B(2), "##.##")
    Print Tab(14); Format$(B(3), "##.##")
    Print Tab(14); Format$(B(4), "##.##")
    Print Tab(14); Format$(B(5), "##.##")
End Sub
Sub CHODCM(A(), N, D(), T())
    For I = 1 To N
        Sum = A(I, I)
        For J = 1 To I - 1
            T(J) = A(J, I)
            For K = 1 To J - 1
                T(J) = T(J) - T(K) * A(J, K)
            Next K
            If D(J) = 0# Then
                If T(J) <> 0# Then
                    Print "No Cholesky decomposition"
                Else
                    A(I, J) = 1#
                End If
            Else
                A(I, J) = T(J) / D(J)
            End If
            Sum = Sum - T(J) * A(I, J)
        Next J
        D(I) = Sum
    Next I
End Sub
Sub CHOBSB(A(), N, D(), B())
    For I = 1 To N
        Sum = B(I)
        For J = 1 To I - 1
            Sum = Sum - A(I, J) * B(J)
        Next J
        B(I) = Sum
    Next I
    For I = N To 1 Step -1
        If D(I) = 0# Then
            Print "  Singular Matrix"
            Exit Sub
        Else
            Sum = B(I) / D(I)
        End If
        For J = I + 1 To N
            Sum = Sum - A(J, I) * B(J)
        Next J
        B(I) = Sum
    Next I
End Sub