📄 form05.frm
字号:
VERSION 5.00
Object = "{8AE029D0-08E3-11D1-BAA2-444553540000}#3.0#0"; "VSFLEX3.OCX"
Begin VB.Form Form05
AutoRedraw = -1 'True
BackColor = &H80000013&
Caption = "贝叶斯判别分析"
ClientHeight = 6255
ClientLeft = 60
ClientTop = 450
ClientWidth = 10485
FillColor = &H8000000F&
ForeColor = &H8000000F&
Icon = "Form05.frx":0000
LinkTopic = "Form1"
ScaleHeight = 6255
ScaleMode = 0 'User
ScaleWidth = 5664.506
StartUpPosition = 3 '窗口缺省
Begin vsFlexLib.vsFlexArray vsFlexArray1
Height = 1695
Left = 3600
TabIndex = 14
Top = 240
Width = 6495
_Version = 196608
_ExtentX = 11456
_ExtentY = 2990
_StockProps = 228
Appearance = 1
ConvInfo = 1413783674
ExtendLastCol = -1 'True
End
Begin VB.TextBox Text4
Height = 375
Left = 1800
TabIndex = 2
Top = 1320
Width = 855
End
Begin VB.CommandButton Command2
Caption = "分 析"
Height = 495
Left = 1560
TabIndex = 9
Top = 3240
Width = 1095
End
Begin VB.CommandButton Command1
Caption = "重 置"
Height = 495
Left = 360
TabIndex = 8
Top = 3240
Width = 1095
End
Begin VB.Frame Frame2
BackColor = &H80000013&
Caption = "结果显示"
Height = 2295
Left = 360
TabIndex = 7
Top = 3840
Width = 9855
Begin VB.TextBox Text5
BackColor = &H80000013&
Height = 1935
Index = 0
Left = 120
MultiLine = -1 'True
ScrollBars = 3 'Both
TabIndex = 10
Top = 240
Width = 9615
End
End
Begin VB.TextBox Text3
Height = 1455
Index = 1
Left = 3600
MultiLine = -1 'True
ScrollBars = 3 'Both
TabIndex = 6
Top = 2160
Width = 6495
End
Begin VB.Frame Frame1
BackColor = &H80000013&
Caption = "数据输入方式"
Height = 1095
Left = 360
TabIndex = 3
Top = 1920
Width = 2295
Begin VB.OptionButton Option2
BackColor = &H80000013&
Caption = "使用文本框"
Height = 255
Left = 480
TabIndex = 5
Top = 720
Width = 1455
End
Begin VB.OptionButton Option1
BackColor = &H80000013&
Caption = "使用表格"
Height = 255
Left = 480
TabIndex = 4
Top = 240
Value = -1 'True
Width = 1335
End
End
Begin VB.TextBox Text2
Height = 375
Left = 1800
TabIndex = 1
Top = 720
Width = 855
End
Begin VB.TextBox Text1
Height = 375
Left = 1800
TabIndex = 0
Top = 120
Width = 855
End
Begin VB.Label Label1
BackColor = &H80000013&
Caption = "样品总数"
Height = 495
Left = 480
TabIndex = 13
Top = 120
Width = 975
End
Begin VB.Label Label2
BackColor = &H80000013&
Caption = "指标总数"
Height = 495
Left = 480
TabIndex = 12
Top = 720
Width = 975
End
Begin VB.Label Label3
BackColor = &H80000013&
Caption = "类别总数"
Height = 495
Left = 480
TabIndex = 11
Top = 1320
Width = 975
End
End
Attribute VB_Name = "Form05"
Attribute VB_GlobalNameSpace = False
Attribute VB_Creatable = False
Attribute VB_PredeclaredId = True
Attribute VB_Exposed = False
'*************************************************************************************************
'
'
' 贝叶斯(Bayes)判别分析
'
'
'
' 函数功能部分
'
'**************************************************************************************************
Option Base 1
Dim Y() As String, x() As Double, m%, n%, k%, k1() As Integer, k2() As Integer
Sub putout(A() As Double, m1%, n1%)
'****************************************************************************************
' 功能 ： 输出矩阵
' 参数 ： A为待输出的矩阵；m1为a矩阵的行；n1为a矩阵的列
' 返回值： 无返回值
'****************************************************************************************
For i = 1 To m1
For j = 1 To n1
Print Tab(j * 22); A(i, j);
Next j
Next i
Print
End Sub
Sub chengfa(a1() As Double, a2() As Double, F() As Double)
'****************************************************************************************
' 功能 ： 求两矩阵的乘积
' 参数 ： a1，a2为待求乘积的两矩阵（a1的列数必须与a2的行数相同，否则结果无效）
' 返回值： F返回乘积结果的矩阵，行数与a1的相同，列数与a2的相同
'**************************************************************************************
If UBound(a1, 2) <> UBound(a2, 1) Then
MsgBox ("要相乘的第一矩阵的列数和第二矩阵的行数必须相同")
Exit Sub
Else
For i = 1 To UBound(a1, 1)
For j = 1 To UBound(a2, 2)
F(i, j) = 0 '原来f矩阵可能被存入非零数据，所以要逐个赋0
For kk = 1 To UBound(a1, 2)
F(i, j) = F(i, j) + a1(i, kk) * a2(kk, j)
Next kk, j, i
End If
End Sub
Sub pingjun(A() As Double, m1%, qd%, zd%, pj() As Double)
'****************************************************************************************
' 功能 ： 求m行矩阵从第qd列到第zd列的平均
' 参数 ： A为待求矩阵；m为矩阵的行数；qd为起始列下标；zd终止列下标
' 返回值： pj返回最终结果，是一m为列向量
'****************************************************************************************
ReDim pj(m1, 1) As Double '原来的pj列向量可能已存入非零数据，故重新声明使其数据均为零
For i = 1 To m1
For j = qd To zd
pj(i, 1) = pj(i, 1) + A(i, j)
Next j
pj(i, 1) = pj(i, 1) / (zd - qd + 1)
Next i
End Sub
Sub licha(A() As Double, m1%, qd%, zd%, lc() As Double)
'****************************************************************************************
' 功能 ： 求m行矩阵从第qd列到第zd列的离差
' 参数 ： A为待求矩阵；m为矩阵的行数；qd为起始列下标；zd终止列下标
' 返回值： lc返回最终结果，是一m1为列向量
'****************************************************************************************
ReDim lc(m1, m1) As Double '原来的lc矩阵可能已存入非零数据，所以要重新声明使其数据均为零
Dim pj1() As Double: ReDim pj(m, 1)
Call pingjun(A(), m1, qd, zd, pj1())
For r = 1 To m1
For t = 1 To m1
For j = qd To zd
sum = sum + (A(r, j) - pj1(r, 1)) * (A(t, j) - pj1(t, 1))
Next j
lc(r, t) = sum: sum = 0
Next t, r
End Sub
Sub zhuanzhi(A() As Double, B() As Double)
'****************************************************************************************
' 功能 ： 求矩阵的转置
' 参数 ： A为待求矩阵
' 返回值： B返回最终结果即A的转置
'****************************************************************************************
Dim m1%, n1%
m1 = UBound(A, 1): n1 = UBound(A, 2)
ReDim B(n1, m1)
For i = 1 To n1
For j = 1 To m1
B(i, j) = A(j, i)
Next j, i
End Sub
Sub MRinv(n1 As Integer, A() As Double, mtxA() As Double, flag As Boolean)
'****************************************************************************************
' 功能 ： 用全选主元高斯－约当法实现矩阵求逆
' 参数 ： A为待求矩阵；n为其阶数
' 返回值： flag返回False表示求逆失败（无逆），返回True表示求逆成功结果存入mtxA矩阵
'****************************************************************************************
ReDim nIs(n1) As Integer, nJs(n1) As Integer
Dim i As Integer, j As Integer, kk As Integer
Dim d As Double, p As Double
For i = 1 To n1
For j = 1 To n1
mtxA(i, j) = A(i, j)
Next j, i
' 全选主元，消元
For kk = 1 To n1
d = 0#
For i = kk To n1
For j = kk To n1
p = Abs(mtxA(i, j))
If (p > d) Then
d = p
nIs(kk) = i
nJs(kk) = j
End If
Next j
Next i
' 求逆失败
If (d + 1# = 1#) Then
flag = False
Exit Sub
End If
If (nIs(kk) <> kk) Then
For j = 1 To n1
p = mtxA(kk, j)
mtxA(kk, j) = mtxA(nIs(kk), j)
mtxA(nIs(kk), j) = p
Next j
End If
If (nJs(kk) <> kk) Then
For i = 1 To n1
p = mtxA(i, kk)
mtxA(i, kk) = mtxA(i, nJs(kk))
mtxA(i, nJs(kk)) = p
Next i
End If
mtxA(kk, kk) = 1# / mtxA(kk, kk)
For j = 1 To n1
If (j <> kk) Then mtxA(kk, j) = mtxA(kk, j) * mtxA(kk, kk)
Next j
For i = 1 To n1
If (i <> kk) Then
For j = 1 To n1
If (j <> kk) Then mtxA(i, j) = mtxA(i, j) - mtxA(i, kk) * mtxA(kk, j)
Next j
End If
Next i
For i = 1 To n1
If (i <> kk) Then mtxA(i, kk) = -mtxA(i, kk) * mtxA(kk, kk)
Next i
Next kk
For kk = n1 To 1 Step -1 ' 调整恢复行列次序
If (nJs(kk) <> kk) Then
For j = 1 To n1
p = mtxA(kk, j)
mtxA(kk, j) = mtxA(nJs(kk), j)
mtxA(nJs(kk), j) = p
Next j
End If
If (nIs(kk) <> kk) Then
For i = 1 To n1
p = mtxA(i, kk)
mtxA(i, kk) = mtxA(i, nIs(kk))
mtxA(i, nIs(kk)) = p
Next i
End If
Next kk
flag = True ' 求逆成功
End Sub
'*********************************************************************'
'以下为求矩阵行列式的函数 '
' '
'其中m()为要计算的矩阵 '
' '
'*********************************************************************'
Public Function Determinant(A() As Double, mm() As Double) As Double
Dim i, j, kk, row, order, r, c As Long
Dim Pivot, Pivot2 As Single
Dim temp() As Double
Determinant = 1
row = UBound(A, 1)
ReDim mm(row, row)
If UBound(A, 2) <> row Then MsgBox "这不是方阵": Exit Function
ReDim temp(1 To row)
For i = 1 To row
For j = 1 To row
mm(i, j) = A(i, j)
Next j, i
For i = 1 To row
Pivot = 0
For j = i To row
For kk = i To row
If Abs(mm(kk, j)) > Pivot Then
Pivot = Abs(mm(kk, j))
r = kk: c = j
End If
Next kk
Next j
If Pivot = 0 Then Determinant = 0: Exit Function
If r <> i Then
order = order + 1
For j = 1 To row
temp(j) = mm(i, j)
mm(i, j) = mm(r, j)
mm(r, j) = temp(j)
Next j
End If
If c <> i Then
order = order + 1
For j = 1 To row
temp(j) = mm(j, i)
mm(j, i) = mm(j, c)
mm(j, c) = temp(j)
Next j
End If
Pivot = mm(i, i)
Determinant = Determinant * Pivot
For j = i + 1 To row
Pivot2 = mm(j, i)
If Pivot2 <> 0 Then
For kk = 1 To row
mm(j, kk) = mm(j, kk) - mm(i, kk) * Pivot2 / Pivot
Next
End If
Next
Next
Determinant = Determinant * (-1) ^ order
End Function
Private Sub Command1_Click()
For i = 1 To n
For j = 1 To m + 1
vsFlexArray1.row = i
vsFlexArray1.Col = j
vsFlexArray1.Text = ""
Next j, i
Text3(1).Text = ""
End Sub
Private Sub Command2_Click()
Dim Y() As String, x() As Double, m%, n%, k%, k1() As Integer, k2() As Integer, q%
n = Val(Text1.Text)
m = Val(Text2.Text)
k = Val(Text4.Text)
ReDim k1(k): ReDim k2(k): ReDim x(m, n)
Dim xx() As Double: ReDim xx(m + 1, n)
If Text3(1).Enabled = True Then
Y = Split(Text3(1).Text, ",")
For i = 0 To UBound(Y)
xx(i Mod (m + 1) + 1, i \ (m + 1) + 1) = Val(Y(i))
Next i
Call putout(xx(), m + 1, n)
For i = 1 To n
For j = 1 To m
x(j, i) = xx(j, i)
Next j, i
⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -