python用fsolve、leastsq對(duì)非線性方程組求解
背景:
實(shí)現(xiàn)用python的optimize庫(kù)的fsolve對(duì)非線性方程組進(jìn)行求解??梢钥吹竭@一個(gè)問(wèn)題實(shí)際上還是一個(gè)優(yōu)化問(wèn)題,也可以用之前擬合函數(shù)的leastsq求解。下面用這兩個(gè)方法進(jìn)行對(duì)比:
代碼:
from scipy.optimize import fsolve,leastsq
from math import sin,cos
def f(x):
x0 = float(x[0])
x1 = float(x[1])
x2 = float(x[2])
return [
5*x1+3,
4*x0*x0 - 2*sin(x1*x2),
x1*x2-1.5
]
x0 = [1,1,1]
result = fsolve(f,x0)
print("===================")
print()
print("求解函數(shù)名稱(chēng):",fsolve.__name__)
print("解:",result)
print("各向量值:",f(result))
#擬合函數(shù)來(lái)求解
h = leastsq(f,x0)
print("===================")
print()
print("求解函數(shù)名稱(chēng):",leastsq.__name__)
print("解:",h[0])
print("各向量的值:",f(h[0]))
結(jié)果:
===================
求解函數(shù)名稱(chēng): fsolve
解: [-0.70622057 -0.6 -2.5 ]
各向量值: [0.0, -9.126033262418787e-14, 5.329070518200751e-15]
===================
求解函數(shù)名稱(chēng): leastsq
解: [-0.70622057 -0.6 -2.5 ]
各向量的值: [0.0, -2.220446049250313e-16, 0.0]
以上就是本文的全部?jī)?nèi)容,希望對(duì)大家的學(xué)習(xí)有所幫助,也希望大家多多支持腳本之家。
相關(guān)文章
python讀取.mat文件及將變量存為.mat文件的詳細(xì)介紹
這篇文章主要給大家介紹了關(guān)于python讀取.mat文件及將變量存為.mat文件的詳細(xì)介紹,?mat文件是matlab的數(shù)據(jù)存儲(chǔ)的標(biāo)準(zhǔn)格式,文中通過(guò)代碼介紹的非常詳細(xì),需要的朋友可以參考下2024-06-06
PyCharm 2021.2 (Professional)調(diào)試遠(yuǎn)程服務(wù)器程序的操作技巧
本文給大家分享用 PyCharm 2021 調(diào)試遠(yuǎn)程服務(wù)器程序的過(guò)程,通過(guò)圖文并茂的形式給大家介紹的非常詳細(xì),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友參考下吧2021-08-08
推薦技術(shù)人員一款Python開(kāi)源庫(kù)(造數(shù)據(jù)神器)
今天小編給大家推薦一款Python開(kāi)源庫(kù),技術(shù)人必備的造數(shù)據(jù)神器!非常不錯(cuò),對(duì)大家的學(xué)習(xí)或工作具有一定的參考借鑒價(jià)值,需要的朋友參考下吧2020-07-07
Django 中自定義 Admin 樣式與功能的實(shí)現(xiàn)方法
這篇文章主要介紹了Django 中自定義 Admin 樣式與功能的實(shí)現(xiàn)方法,本文圖文并茂給大家介紹的非常詳細(xì),具有一定的參考借鑒價(jià)值,需要的朋友可以參考下2019-07-07

