diff --git a/227-3.py b/227-3.py index c56f62c..77bca4f 100644 --- a/227-3.py +++ b/227-3.py @@ -63,7 +63,7 @@ def prettyPrintMatrix(matrix): if __name__ == "__main__": # A = np.array(A, dtype=float) # b = np.array(b, dtype=float) - + #把矩阵A和b改成题干要求的##################################### A = [ [0, 3, 4], [1, -1, 1], diff --git a/227-4.py b/227-4.py index 115fa58..ad66601 100644 --- a/227-4.py +++ b/227-4.py @@ -22,7 +22,7 @@ def SqrtSolve(A,b): for k in range(j): L[getIndexFromDownMatrix(i,j)] -= L[getIndexFromDownMatrix(i,k)]*L[getIndexFromDownMatrix(j,k)] L[getIndexFromDownMatrix(i,j)] /= L[getIndexFromDownMatrix(j,j)] - + print(L) for i in range(n): for k in range(i): b[i] -= L[getIndexFromDownMatrix(i,k)]*b[k] @@ -34,7 +34,7 @@ def SqrtSolve(A,b): b[i] /= L[getIndexFromDownMatrix(i,i)] return b - +#把A,b换成题干的数值########################################### if __name__ == "__main__": # 储存下三角矩阵 a11, a21, a22, a31, a32, a33 ... A = [4,2,2,-2,-3,14] diff --git a/227-7.py b/227-7.py index 05e1cd2..51f3082 100644 --- a/227-7.py +++ b/227-7.py @@ -8,19 +8,21 @@ def ZGsolve(A,b): beta[i] = A[i][2] / A[i][1] else: beta[i] = A[i][2] / (A[i][1] - A[i][0]*beta[i-1]) - + print("beta:") + print(beta[:-1]) for i in range(n): if i == 0: b[i] = b[i] / A[i][1] else: b[i] = (b[i] - A[i][0]*b[i-1]) / (A[i][1] - A[i][0]*beta[i-1]) - + print("y:") + print(b) for i in range(n-2,-1,-1): b[i] = b[i] - beta[i]*b[i+1] return b - +#把A,b换成题干的数值########################################### if __name__ == "__main__": # 储存追赶法A矩阵 A = [