python實現(xiàn)迭代法求方程組的根過程解析
更新時間:2019年11月25日 10:29:33 作者:Picassooo
這篇文章主要介紹了python實現(xiàn)迭代法求方程組的根過程解析,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下
這篇文章主要介紹了python實現(xiàn)迭代法求方程組的根過程解析,文中通過示例代碼介紹的非常詳細,對大家的學習或者工作具有一定的參考學習價值,需要的朋友可以參考下
有方程組如下:
迭代法求解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)
結果如下:
以上就是本文的全部內容,希望對大家的學習有所幫助,也希望大家多多支持腳本之家。
相關文章
JavaScript toUpperCase()方法使用詳解
這篇文章主要為大家詳細介紹了JavaScript toUpperCase()方法的使用技巧,感興趣的小伙伴們可以參考一下2016-08-08JavaScript+html5 canvas制作色彩斑斕的正方形效果
這篇文章主要介紹了JavaScript+html5 canvas制作色彩斑斕的正方形效果,實例分析了JavaScript結合html5 canvas實現(xiàn)圖形動態(tài)繪制的技巧,需要的朋友可以參考下2016-01-01讓JavaScript擁有類似Lambda表達式編程能力的方法
在前幾天的博文中我發(fā)布了一個可以自定義頁碼呈現(xiàn)方式的組件,有C#和JavaScript兩個版本。2010-09-09