矩阵求逆标准算法.doc

矩阵求逆标准算法（VB）源码 2006-11-29 13:49 类别：默认 ?? 本程序依据矩阵初等变换的基本原理编写，算法较为繁琐，但易于理解适合VB初学者。 ?? 本程序适合任何（n*n）的矩阵求逆，对于不可逆矩阵有提示信息，并结束程序 ?? 本程序在XP，VB6.0下调试通过 ?? 本程序由本人原创，请慎用。如有疑问，或调试有误，请联系本人QQ?? 本程序可在VB6.0内任何地方用call jzqn(qa(),na()))语句调用 ?? 其中qa()是输入的矩阵数组，调用此函数后 na()为返回的逆矩阵数组 ?? 注意：调用本程序前不要声明na()的维数，仅用dim na()即可。 ?? 请不要试图对一个病态矩阵求逆、否则计算结果未必是你想要的 ?? 病态矩阵是指行列式计算结果极其接近于零的矩阵 Public Sub jzqn(qa(), na())Dim a()n = UBound(qa, 1)ReDim na(n, n)ReDim a(n, 2 * n)? For i = 1 To n??? For j = 1 To n????? a(i, j) = qa(i, j)??? Next j? Next iFor i = 1 To n?For j = n + 1 To 2 * n? If j - i = n Then? a(i, j) = 1? Else? a(i, j) = 0? End If? Next j? Next i?For i = 1 To n? If a(i, i) = 0 Then??? For q = i To n??? If a(q, i) 0 Then?????? For w = i To 2 * n???????? zj = a(i, w)???????? a(i, w) = a(q, w)???????? a(q, w) = zj?????? Next w?????? Exit For??? End If??? Next q??? If q n Then MsgBox 此矩阵不可逆: Exit Sub? End If? For k = 2 * n To i Step -1? a(i, k) = a(i, k) / a(i, i)? Next k? For j = i + 1 To n? If a(j, i) 0 Then???? For k = 2 * n To i Step -1??? a(j, k) = a(j, k) / a(j, i) - a(i, k)??? Next k??? End If?Next jNext iFor i = n To 1 Step -1? If a(i, i) = 0 Then??? For q = i - 1 To 1 Step -1??? If a(q, i) 0 Then?????? For w = i To 2 * n???????? zj = a(i, w)???????? a(i, w) = a(q, w)???????? a(q, w) = zj?????? Next w?????? Exit For??? End If??? Next q? End If? For k = 2 * n To i Step -1? a(i, k) = a(i, k) / a(i, i)? Next k? For j = i - 1 To 1 Step -1? If a(j, i) 0 Then? xxx = a(j, i)? For k = 2 * n To 1 Step -1??? a(j, k) = a(j, k) / xxx - a(i, k)??? Next k??? End If?Next jNext iFor i = 1 To n?For j = 1 To n?? na(i, j) = a(i, j + n)Next jNext iEnd Sub ?? 调用示例：下面代码随机产生一个10*10的矩阵，并求逆，打印于窗体 Private Sub Command1_Click()Dim a(10, 10), b()ClsRandomizeFor i = 1 To 10?? For j = 1 To 10???? a(i, j) = Int(Rnd * 100)??? Print a(i, j);?? Next j?? PrintNext iPrintCall jzqn(a(), b()) For i = 1 To 10?? For j = 1 To 10