d11r1.frm

方程组的矩阵解法

VERSION 5.00
Begin VB.Form Form1 
   Caption         =   "Form1"
   ClientHeight    =   3195
   ClientLeft      =   60
   ClientTop       =   345
   ClientWidth     =   5865
   LinkTopic       =   "Form1"
   ScaleHeight     =   3195
   ScaleWidth      =   5865
   StartUpPosition =   3  'Windows Default
   Begin VB.CommandButton Command1 
      Caption         =   "Command1"
      Height          =   375
      Left            =   3960
      TabIndex        =   0
      Top             =   2640
      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 D11R1
    'Driver for routine MOMENT
    PI = 3.14159265
    NPTS = 10000
    NBIN = 100
    NDAT = NPTS + NBIN
    Dim DATA(100100)
    I = 1
    For J = 1 To NBIN
        X = PI * J / NBIN
        NLIM = CInt(Sin(X) * PI / 2# * NPTS / NBIN)
        For K = 1 To NLIM
            DATA(I) = X
            I = I + 1
        Next K
    Next J
    Print
    Print Tab(5); "Moments of a sinusoidal distribution"
    Call MOMENT(DATA(), NPTS, AVE, ADEV, SDEV, VAR, SKEW, CURT)
    Print
    Print Tab(30); "Calculated    Expected"
    Print Tab(5); "Mean : "; Tab(32); Format$(AVE, "#.###0");
    Print Tab(44); Format$(PI / 2#, ".####")
    Print Tab(5); "Average Deviation : ";
    Print Tab(32); Format$(ADEV, ".###0"), Tab(44); 0.5708
    Print Tab(5); "Standard Deviation : ";
    Print Tab(32); Format$(SDEV, ".###0"), Tab(44); 0.6837
    Print Tab(5); "Variance : ";
    Print Tab(32); Format$(VAR, ".###0"), Tab(44); 0.4674
    Print Tab(5); "Skewness : ";
    Print Tab(32); Format$(SKEW, "##.###0"), Tab(44); 0#
    Print Tab(5); "Kurtosis : ";
    Print Tab(32); Format$(CURT, "##.###0"), Tab(44); -0.8062
End Sub
Sub MOMENT(DATA(), N, AVE, ADEV, SDEV, VAR, SKEW, CURT)
    If N <= 1 Then
        Print " 'N must be at least 2'"
      Exit Sub
    End If
    S = 0#
    For J = 1 To N
        S = S + DATA(J)
    Next J
    AVE = S / N
    ADEV = 0#
    VAR = 0#
    SKEW = 0#
    CURT = 0#
    For J = 1 To N
        S = DATA(J) - AVE
        ADEV = ADEV + Abs(S)
        P = S * S
        VAR = VAR + P
        P = P * S
        SKEW = SKEW + P
        P = P * S
        CURT = CURT + P
    Next J
    ADEV = ADEV / N
    VAR = VAR / (N - 1)
    SDEV = Sqr(VAR)
    If VAR <> 0# Then
        SKEW = SKEW / (N * SDEV ^ 3)
        CURT = CURT / (N * VAR ^ 2) - 3#
    Else
        Print " no skew or kurtosis when zero variance"
    End If
End Sub