python實現(xiàn)迭代法求方程組的根過程解析

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

站在用戶的角度思考問題，與客戶深入溝通，找到合江網(wǎng)站設(shè)計與合江網(wǎng)站推廣的解決方案，憑借多年的經(jīng)驗，讓設(shè)計與互聯(lián)網(wǎng)技術(shù)結(jié)合，創(chuàng)造個性化、用戶體驗好的作品，建站類型包括：做網(wǎng)站、成都網(wǎng)站制作、企業(yè)官網(wǎng)、英文網(wǎng)站、手機端網(wǎng)站、網(wǎng)站推廣、主機域名、網(wǎng)頁空間、企業(yè)郵箱。業(yè)務(wù)覆蓋合江地區(qū)。

有方程組如下：

迭代法求解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)

結(jié)果如下：

python實現(xiàn)迭代法求方程組的根過程解析

以上就是本文的全部內(nèi)容，希望對大家的學習有所幫助，也希望大家多多支持創(chuàng)新互聯(lián)。

網(wǎng)頁標題：python實現(xiàn)迭代法求方程組的根過程解析
文章出自：http://chinadenli.net/article44/ghoohe.html

成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián)，為您提供關(guān)鍵詞優(yōu)化、虛擬主機、網(wǎng)站改版、全網(wǎng)營銷推廣、搜索引擎優(yōu)化、品牌網(wǎng)站制作

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

猜你還喜歡下面的內(nèi)容

欧美一区二区三区老妇人-欧美做爰猛烈大尺度电-99久久夜色精品国产亚洲a-亚洲福利视频一区二区

python實現(xiàn)迭代法求方程組的根過程解析