1112112我玩玩

2112

import numpy as np  # 本次仅需要用到这一个库

# a = input('需要计算的方程组为：')
a = np.array([[0.01, 2, -0.5, -5], [-1, -0.5, 2, 5], [5, -4, 0.5, 9]])  # 双括号   a即要计算的
i = a.shape[0]  # 得到a的行数
j = a.shape[1]  # 得到a的列数
print(i)
print(j)
X = []  # 建立一个X空列表，便于放入所得解
print(X)

for z in range(0, j - 2):  # 当z从第1列到j-1列循环，每一次循环中，都是#针对所在列，先排序，后消元（以下对应这两部分） 2.z还代表循环次数
    for x in range(z, i):  # x从z行到第最后一行循环，每次循环扔掉上次循#环的第一行，因为每次的第一行都用来存储每次的最大
        n = np.argmax(abs(a[x:i, z]))  #n为剩下的几行中，绝对值最大的值#所在行的索引
        print(n + x)  # 索引再加x，因为n为剩下的几行中的索引，而我们
#接下来的交换针对所有行
        a[[x, n + x], :] = a[[n + x, x], :]  # 数组特殊的交换数值方式
        print(a)

    for y in range(1 + z, i):  # 之所以从1+z，开始，最开始的z循环不仅仅#是列，也是循环次数，每循环一次，向下移动一列
        k = a[z, z] / a[y, z]
        print(k)
        a[y,] = a[z,] - k * a[y,]
print(a)

w = j - 1
g = 0
for q in range(i - 1, -1, -1):  # i为矩阵行数  这里是从最后一行开始向上#循环行！
    w = w - 1
    h = 2 - q  #3元，所以2

    if h == 0:
        X.append(a[q, j - 1] / a[q, w])
        print(X)
    else:
        for n in range(0, h):  # e+1决定让n循环几次
            g = g + a[q, 2 - n] * X[n]
            print(g)
        X.append((a[q, j - 1] - g) / a[q, w])  # 公式
        print(X)
        g = 0