乐筑天下

搜索
欢迎各位开发者和用户入驻本平台 尊重版权,从我做起,拒绝盗版,拒绝倒卖 签到、发布资源、邀请好友注册,可以获得银币 请注意保管好自己的密码,避免账户资金被盗
查看: 14|回复: 0

[VBA]各位老大我改了一段矩阵求逆的程序~请帮忙改正一下谢谢

[复制链接]

1

主题

1

帖子

1

银币

初来乍到

Rank: 1

铜币
5
发表于 2006-5-30 11:26:00 | 显示全部楼层 |阅读模式
Dim N As Integer
Dim ch As Double, Z As Double, zz As Double
Dim P As Integer
Dim qn(1 To 6, 1 To 6) As Double
Dim nb(1 To 6, 1 To 6) As Double '返回的逆阵
'对mc(1to 6,1 to 6)求逆

Dim E(1 To 6, 1 To 12) As Double
'应用的原理是: [ A I]经过矩阵行变换后成为[ I ivt_A], E阵(N,2N)阶 既 [ A I]
For i = 1 To 6
  For j = 1 To 6
E(i, j) = mc(i, j)
  Next j
Next i
For i = 1 To 6
  For j = 7 To 12
If i = j - 6 Then
E(i, j) = 1
Else
E(i, j) = 0
End If
  Next j
Next i
'以上步骤实现用A,I初始化E阵

For j = 1 To 6
If E(j, j) = 0 Then
For i = 1 To 6
If E(i, j)  0 Then Exit For
Next i
For P = 1 To 12
ch = E(i, P)
E(i, P) = E(j, P)
E(j, P) = ch
Debug.Print E(i, P)
Next P '将对角线元素为0的行与对应列上元素非零的行,进行行变换
End If

Z = E(j, j)
For P = 1 To 12
E(j, P) = E(j, P) / Z '将矩阵E (A,I) A矩阵每行都除以对角线元素,使对角线元素变为1
Next P
For i = 1 To 6
If i  j And E(i, j)  0 Then
zz = E(i, j)
For P = 1 To 12
E(i, P) = E(i, P) - E(j, P) * zz
Next P
End If
Next i
Next j

For i = 1 To 6
For j = 1 To 6
nb(i, j) = E(i, j + 6)

Next j
Next i

  ' qn(i, j) = nb(i, j) '返回求出的A阵的逆阵B阵
回复

使用道具 举报

发表回复

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

  • 微信公众平台

  • 扫描访问手机版

  • 点击图片下载手机App

QQ|关于我们|小黑屋|乐筑天下 繁体中文

GMT+8, 2025-7-5 19:46 , Processed in 1.137113 second(s), 65 queries .

© 2020-2025 乐筑天下

联系客服 关注微信 帮助中心 下载APP 返回顶部 返回列表