python實現迭代法求方程組的根過程解析-創新互聯

這篇文章主要介紹了python實現迭代法求方程組的根過程解析,文中通過示例代碼介紹的非常詳細，對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下

成都創新互聯公司專注于網站建設|成都網站維護|優化|托管以及網絡推廣，積累了大量的網站設計與制作經驗，為許多企業提供了網站定制設計服務，案例作品覆蓋玻璃貼膜等行業。能根據企業所處的行業與銷售的產品，結合品牌形象的塑造，量身制作品質網站。

有方程組如下：

迭代法求解x，python代碼如下：

import numpy as np
import matplotlib.pyplot as plt
 
A = np.array([[8, -3, 2], [4, 11, -1], [6, 3, 12]])
b = np.array([[20, 33, 36]])
 
# 方法一：消元法求解方程組的解
result = np.linalg.solve(A, b.T)
print('Result:\n', result)
 
# 方法二：迭代法求解方程組的解
B = np.array([[0, 3/8, -2/8], [-4/11, 0, 1/11], [-6/12, -3/12, 0]])
f = np.array([[20/8, 33/11, 36/12]])
error = 1.0e-6
steps = 100
xk = np.zeros((3, 1)) # initialize parameter setting
errorlist = []
for k in range(steps):
  xk_1 = xk
  xk = np.matmul(B, xk) + f.T
  print('xk:\n', xk)
  errorlist.append(np.linalg.norm(xk-xk_1))
  if errorlist[-1] < error:
    print('iteration: ', k+1)
    break
 
# 把誤差畫出來
x_axis = [i for i in range(len(errorlist))]
plt.figure()
plt.plot(x_axis, errorlist)

本文名稱：python實現迭代法求方程組的根過程解析-創新互聯
地址分享：http://m.kartarina.com/article0/ejooo.html

成都網站建設公司_創新互聯，為您提供品牌網站制作、定制開發、網站排名、網頁設計公司、網站維護、網站建設

廣告

聲明：本網站發布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網站立場，如需處理請聯系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經允許不得轉載，或轉載時需注明來源： 創新互聯