lhye200/CalWay_Python
1
1
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

dfg

Browse Source
This commit is contained in:
lwj
2025-06-18 23:16:37 +08:00
parent b630b2791d
commit 97428790f3
1 changed files with 9 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

10
227-3.py
View File

@@ -19,6 +19,10 @@ def SovleRowMain(A,b):
b[i], b[row_max_index] = b[row_max_index], b[i] b[i], b[row_max_index] = b[row_max_index], b[i]
p[i], p[row_max_index] = p[row_max_index], p[i] p[i], p[row_max_index] = p[row_max_index], p[i]
print(f"第{i+1}次交换，交换行{i+1}和行{row_max_index+1}，A矩阵为：")
for row in A:
print(row)
print(f"b向量为：{b}\n")
if abs(A[i][i]) < ks: if abs(A[i][i]) < ks:
raise ValueError("A矩阵奇异，无法进行高斯消元") raise ValueError("A矩阵奇异，无法进行高斯消元")
for j in range(i + 1, n): for j in range(i + 1, n):
@@ -27,7 +31,11 @@ def SovleRowMain(A,b):
for k in range(i + 1, n): for k in range(i + 1, n):
A[j][k] -= m * A[i][k] A[j][k] -= m * A[i][k]
b[j] -= m * b[i] b[j] -= m * b[i]
print(f"系数为{-1*m}用加号")
print("消元后的A矩阵：")
for row in A:
print(row)
print(f"消元后的b向量：{b}\n")
if abs(A[n - 1][n - 1]) < ks: if abs(A[n - 1][n - 1]) < ks:
raise ValueError("A矩阵奇异，无法进行高斯消元") raise ValueError("A矩阵奇异，无法进行高斯消元")