问题描述
试图找到矩阵的逆。下面的代码有效,逆结果正确,但输出中的一些零以负零的形式出现。不知道问题是什么,任何帮助将不胜感激。谢谢
# Importing NumPy Library
import numpy as np
import sys
n=3
a = np.zeros((n,2*n))
#a = np.eye(3,dtype = float)
print(a)
b = [[1,-1,0],[0,-4,2],-2]]
print (b)
c = np.concatenate((b,a),axis =1)
print(c)
for i in range(n):
for j in range(n):
a[i][j] = b[i][j]
# Augmenting Identity Matrix of Order n
for i in range(n):
for j in range(n):
if i == j:
a[i][j+n] = 1
# Applying Guass Jordan Elimination
for i in range(n):
if a[i][i] == 0.0:
sys.exit('Divide by zero detected!')
for j in range(n):
if i != j:
ratio = a[j][i]/a[i][i]
for k in range(2*n):
a[j][k] = a[j][k] - ratio * a[i][k]
**OUTPUT**
[[ 1. -0.25 -0.25]
[-0. -0.25 -0.25]
[-0. -0. -0.5 ]]
解决方法
暂无找到可以解决该程序问题的有效方法,小编努力寻找整理中!
如果你已经找到好的解决方法,欢迎将解决方案带上本链接一起发送给小编。
小编邮箱:dio#foxmail.com (将#修改为@)