矩阵求逆标准算法.docxVIP

矩阵求逆标准算法 矩阵求逆标准算法(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(), n a()) Dim a() n = UBou nd(qa, 1) ReDim na(n, n) ReDim a(n, 2 * n) For i = 1 To n Forj = 1 To n a(i, j) = qa(i, j) Next j Next i For i = 1 To n For j = n + 1 To 2 * n If j - i = n The n 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 The n For q = i To n If a(q, i) 0 The n 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 The n MsgBox 此矩阵不可逆： If q n The n 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 The n 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 j Next i For i = n To 1 Step -1 If a(i, i) = 0 The n For q = i - 1 To 1 Step -1 If a(q, i) 0 The n 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 The n 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 j Next i For i = 1 To n Forj = 1 To n n a(i, j) = a(i, j + n) Next j Next i End Sub 调用示例：下面代码随机产生一个 10*10的矩 阵，并求逆，打印于窗体 Private Sub Comma nd1_Click() Dim a(10, 10), b() Cis Ran domize For i = 1 To 10 For j = 1 To 10 a(i, j) = In t(R nd * 100) Print a(i, j); Next j Print Next i Print Call jzqn(a(), b()) For i = 1 To 10 For j = 1 To 10 Print Format(b(i, j), 0.000), Next j Print Next i End Sub 矩阵运算是数值运算中经常碰到的， 砖头”抛岀多天，尚未 引岀玉来”，我自己再来个补充吧！ 矩阵求逆上面给出的程序，虽然可以使用，但远不完善，更不精炼。下面将其修改一下，例如：使用 IIF ()函数 简化判断分支语句，将 约化”过程合并，添加一个矩阵无逆的判断， ……。但还是属于小打小闹的修修补 补，希望 诸位能挑岀程序中