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

eee

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

20
227-4.py
View File

@@ -22,12 +22,24 @@ def SqrtSolve(A,b):
for k in range(j): for k in range(j):
L[getIndexFromDownMatrix(i,j)] -= L[getIndexFromDownMatrix(i,k)]*L[getIndexFromDownMatrix(j,k)] L[getIndexFromDownMatrix(i,j)] -= L[getIndexFromDownMatrix(i,k)]*L[getIndexFromDownMatrix(j,k)]
L[getIndexFromDownMatrix(i,j)] /= L[getIndexFromDownMatrix(j,j)] L[getIndexFromDownMatrix(i,j)] /= L[getIndexFromDownMatrix(j,j)]
print(L) # 打印下三角矩阵
print("下三角矩阵 L:")
for i in range(n):
L_row = []
for j in range(n):
if j <= i:
L_row.append(L[getIndexFromDownMatrix(i,j)])
else:
L_row.append(0)
print(L_row)
# print(L)
for i in range(n): for i in range(n):
for k in range(i): for k in range(i):
b[i] -= L[getIndexFromDownMatrix(i,k)]*b[k] b[i] -= L[getIndexFromDownMatrix(i,k)]*b[k]
b[i] /= L[getIndexFromDownMatrix(i,i)] b[i] /= L[getIndexFromDownMatrix(i,i)]
# 打印 b 向量
print("y 向量:")
print(b)
for i in range(n-1,-1,-1): for i in range(n-1,-1,-1):
for k in range(i+1,n): for k in range(i+1,n):
b[i] -= L[getIndexFromDownMatrix(k,i)]*b[k] b[i] -= L[getIndexFromDownMatrix(k,i)]*b[k]
@@ -42,5 +54,5 @@ if __name__ == "__main__":
b = [10,5,4] b = [10,5,4]
print("x:") # print("x:")
print(SqrtSolve(A,b)) print("x: \n",SqrtSolve(A,b))