VB.Net矩阵求秩函数

 Public Function Math_Matrix_Rank(ByVal K(,) As Integer, ByVal error_ As Integer, GetList As Integer) As Integer  '返回矩阵的秩.
        Dim n As Integer = GetList
        Dim m As Integer = K.Length \ n
        Dim i As Integer = 0
        Dim i1 As Integer
        Dim j As Integer = 0
        Dim j1 As Integer
        Dim temp1 As Single
        If m > n Then '保证m≤n
            i = m
            m = n
            n = i
            i = 1
        End If
        m -= 1
        n -= 1
        Dim temp(m, n) As Single
        If i = 0 Then
            For i = 0 To m
                For j = 0 To n
                    temp(i, j) = K(i, j)
                Next
            Next
        Else
            For i = 0 To m
                For j = 0 To n
                    temp(i, j) = K(j, i)
                Next
            Next
        End If
        If m = 0 Then
            i = 0
            While i <= n
                If K(0, i) <> 0 Then
                    Return 1
                End If
                i += 1
            End While
            Return 0
        End If
        Dim error0 As Double
        If error_ = -1 Then
            error0 = System.Math.Pow(0.1, 10)
        Else
            error0 = System.Math.Pow(0.1, error_)
        End If
        i = 0
        While i <= m '保证误差可控制
            j = 0
            While j <= n
                If temp(i, j) <> 0 Then
                    error0 *= temp(i, j)
                    i = m
                    Exit While
                End If
                j += 1
            End While
            i += 1
        End While
        Dim error1 As Double
        For i = 0 To m '消0过程
            j = 0
            While j <= n
                If temp(i, j) <> 0 Then
                    Exit While
                End If
                j += 1
            End While
            If j <= n Then
                i1 = 0
                While i1 <= m
                    If temp(i1, j) <> 0 And i1 <> i Then
                        temp1 = temp(i, j) / temp(i1, j)
                        error1 = System.Math.Abs((temp(i, j) - temp(i1, j) * temp1)) * 100 '误差控制。因为有时候temp(i, j) - temp(i1, j) * (temp(i, j) / temp(i1, j))≠0
                        error1 += error0
                        For j1 = 0 To n
                            temp(i1, j1) = temp(i, j1) - temp(i1, j1) * temp1
                            If System.Math.Abs(temp(i1, j1)) < error1 Then
                                temp(i1, j1) = 0
                            End If
                        Next
                    End If
                    i1 += 1
                End While
            End If
        Next
        i1 = 0 '作为返回值的临时变量
        For i = 0 To m
            For j = 0 To n
                If temp(i, j) <> 0 Then
                    i1 += 1
                    Exit For
                End If
            Next
        Next
        Return i1
    End Function

推荐文章：那些年，做的几个应用