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