求逆.frm

计算半方差

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   3195
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   4680
   LinkTopic       =   "Form1"
   ScaleHeight     =   3195
   ScaleWidth      =   4680
   StartUpPosition =   3  '窗口缺省
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   735
      Left            =   1800
      TabIndex        =   0
      Top             =   2160
      Width           =   1695
   End
End
Attribute VB_Name = "Form1"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
Option Base 1
Private Sub Command1_Click()
Dim b(6, 6) As Double: Dim q(6, 6) As Double: Dim e As Double
Dim i, j, k As Integer
b(1, 1) = 2.5: b(1, 2) = 7.739: b(1, 3) = 9.999: b(1, 4) = 7.656: b(1, 5) = 5.939: b(1, 6) = 1#
b(2, 1) = 7.739: b(2, 2) = 2.5: b(2, 3) = 8.677: b(2, 4) = 6.381: b(2, 5) = 7.196: b(2, 6) = 1#
b(3, 1) = 9.999: b(3, 2) = 8.677: b(3, 3) = 2.5: b(3, 4) = 7.656: b(3, 5) = 9.206: b(3, 6) = 1#
b(4, 1) = 7.656: b(4, 2) = 6.381: b(4, 3) = 7.656: b(4, 4) = 2.5: b(4, 5) = 4.936: b(4, 6) = 1#
b(5, 1) = 5.939: b(5, 2) = 7.196: b(5, 3) = 9.206: b(5, 4) = 4.936: b(5, 5) = 2.5: b(5, 6) = 1#
b(6, 1) = 1#: b(6, 2) = 1#: b(6, 3) = 1#: b(6, 4) = 1#: b(6, 5) = 1#: b(6, 6) = 0#
For i = 1 To 6
   For j = 1 To 6
     q(i, j) = 0
   Next j
    q(i, i) = 1
Next i
 For k = 1 To 6
    e = b(k, k)
   For j = 1 To 6
    b(k, j) = b(k, j) / e: q(k, j) = q(k, j) / e
   Next j
   For i = 1 To 6
     If k = i Then GoTo 30
    For j = 1 To 6
      If j > k Then Let b(i, j) = b(i, j) - b(i, k) * b(k, j)
    q(i, j) = q(i, j) - b(i, k) * q(k, j)
    Next j
30: Next i
 Next k
 For i = 1 To 6
   For j = 1 To 6
    Print Format(q(i, j), "##.000") & "   ";
   Next j
   Print
  Next i
    
 End Sub