Attribute VB_Name = "modM"
'頻數統計算法
Option Explicit
'P(1 To 樣本數，1 To 2)：原始數據
'S(1 To 分組數，1 To 2)：頻段的端點
'R(1 To 分組數，1 To 2)：頻數統計結果
Public Sub FreFuzz(P, S, R)
    Dim I As Integer, J As Integer, L As Integer
    Dim M As Integer, N As Integer
    Dim P1 As Single, P2 As Single, P0 As Single
    Dim S1 As Single, S2 As Single
'根據數組求下標的上界
'M是數據的樣本數
'N是頻段的組數
    M = UBound(P, 1): N = UBound(S, 1)
'統計頻數
    For I = 1 To M
        P1 = P(I, 1): P2 = P(I, 2)
        For P0 = P1 To P2   '以1為步長，如果步長不等于1需要填“Step 步長”
            For L = 1 To N
                S1 = S(L, 1): S2 = S(L, 2)
                If P0 >= S1 And P0 < S2 Then R(L, 1) = R(L, 1) + 1
            Next L
        Next P0
    Next I
'求相對頻數
    For L = 1 To N
        R(L, 2) = R(L, 1) / M
    Next L
End Sub