Python解方程組 scipy.optimize.fsolve()函數(shù)如何求解帶有循環(huán)求和的方程式
最近做的課題需要求解方程組，方程組還挺復(fù)雜，未知參數(shù)比較多
簡(jiǎn)單粗暴介紹fsolve()函數(shù)用法：
假如要解方程組
x*y=6
x+y=5
先變換成右邊為0的形式
x*y-6
x+y-5
func函數(shù)中直接替換要解的方程就可
import scipy
from scipy.optimize import fsolve
def func(i):
x, y = i[0], i[1]
return [ # 這里寫要求解的方程組式子，變成等于0的形式
x*y-6,
x+y-5 ]
if __name__ == '__main__':
r = scipy.optimize.fsolve(func, [0, 0])
print(r)結(jié)果
[2. 3.]
再看看有參數(shù)的解法：
如果有參數(shù)a,b需要傳入
axy-6
x+by-5
import scipy
from scipy.optimize import fsolve
def func(i, params):
x, y = i[0], i[1]
a, b = params
return [ # 這里寫要求解的方程組式子，變成等于0的形式
a*x*y-6,
x+b*y-5 ]
if __name__ == '__main__':
a = 2
b = 0.5
params = [a, b]
r = scipy.optimize.fsolve(func, [0, 0], args=params)
print(r)結(jié)果
[0.32055053 9.35889894]
近似解
更復(fù)雜的方程式帶有循環(huán)求和：
附三角函數(shù)寫法
import math math.sin(theta)
帶求和符號(hào)的方程求解（這里舉一個(gè)簡(jiǎn)單的例子好理解一些）
這里用sum來(lái)計(jì)算和
import scipy
from scipy.optimize import fsolve
def func(i, params):
x, y = i[0], i[1]
a, b = params
list_e1 = [x * a[j][0] + y * b[j][0] - 5 for j in range(2)]
list_e2 = [x * a[j][1] - y * b[j][1] - 2 for j in range(2)]
return [ # 這里寫要求解的方程組式子，變成等于0的形式
sum(list_e1),
sum(list_e2)
]
if __name__ == '__main__':
a = [[3, 4], [3, 4]]
b = [[1, 3], [4, 1]]
params = [a, b]
r = scipy.optimize.fsolve(func, [0, 0], args=params)
print(r)結(jié)果
[0.9375 0.875 ]
到此這篇關(guān)于Python解方程組 scipy.optimize.fsolve()函數(shù)如何求解帶有循環(huán)求和的方程式的文章就介紹到這了,更多相關(guān)Python scipy.optimize.fsolve()函數(shù)內(nèi)容請(qǐng)搜索腳本之家以前的文章或繼續(xù)瀏覽下面的相關(guān)文章希望大家以后多多支持腳本之家！
相關(guān)文章
Python numpy多維數(shù)組實(shí)現(xiàn)原理詳解
這篇文章主要介紹了python numpy多維數(shù)組實(shí)現(xiàn)原理詳解,文中通過(guò)示例代碼介紹的非常詳細(xì)，對(duì)大家的學(xué)習(xí)或者工作具有一定的參考學(xué)習(xí)價(jià)值,需要的朋友可以參考下2020-03-03
詳解利用Python scipy.signal.filtfilt() 實(shí)現(xiàn)信號(hào)濾波
這篇文章主要介紹了詳解利用Python scipy.signal.filtfilt() 實(shí)現(xiàn)信號(hào)濾波，小編覺(jué)得挺不錯(cuò)的，現(xiàn)在分享給大家，也給大家做個(gè)參考。一起跟隨小編過(guò)來(lái)看看吧2019-06-06
python?web?開發(fā)之Flask中間件與請(qǐng)求處理鉤子的最佳實(shí)踐
Flask作為輕量級(jí)Web框架,提供了靈活的請(qǐng)求處理機(jī)制,中間件和請(qǐng)求鉤子允許開發(fā)者在請(qǐng)求處理的不同階段插入自定義邏輯,實(shí)現(xiàn)諸如身份驗(yàn)證、日志記錄、數(shù)據(jù)預(yù)處理等功能,本文將詳細(xì)介紹Flask的請(qǐng)求處理鉤子、g對(duì)象以及中間件模式的使用方法和最佳實(shí)踐,感興趣的朋友一起看看吧2025-05-05
詳解如何為eclipse安裝合適版本的python插件pydev
這篇文章主要介紹了詳解如何為eclipse安裝合適版本的python插件pydev，pydev是一款優(yōu)秀的Eclipse插件，具有一定的參考價(jià)值，感興趣的小伙伴們可以參考一下2018-11-11
Python 實(shí)現(xiàn)數(shù)據(jù)庫(kù)更新腳本的生成方法
下面小編就為大家?guī)?lái)一篇Python 實(shí)現(xiàn)數(shù)據(jù)庫(kù)更新腳本的生成方法。小編覺(jué)得挺不錯(cuò)的，現(xiàn)在就分享給大家，也給大家做個(gè)參考。一起跟隨小編過(guò)來(lái)看看吧2017-07-07
這篇文章主要介紹了Python元字符的用法實(shí)例解析，具有一定借鑒價(jià)值,需要的朋友可以參考下2018-01-01
Python讀寫Redis數(shù)據(jù)庫(kù)操作示例
Redis是一個(gè)開源的非關(guān)系型數(shù)據(jù)庫(kù)，它采用C語(yǔ)言編寫，是一個(gè)key-value存儲(chǔ)系統(tǒng)，它存儲(chǔ)的value類型很多，包括string（字符串），list（鏈表），set（集合），zset（有序集合），hash（哈希）2014-03-03
在本篇文章里小編給大家分享的是關(guān)于python進(jìn)行矩陣運(yùn)算的方法及實(shí)例代碼，需要的朋友們可以學(xué)習(xí)下。2020-06-06
Python構(gòu)建企業(yè)級(jí)離線包倉(cāng)庫(kù)的實(shí)戰(zhàn)指南與腳本
在內(nèi)網(wǎng)環(huán)境與安全要求日益嚴(yán)格的今天,構(gòu)建離線Python包倉(cāng)庫(kù)已成為Python項(xiàng)目管理的必備技能,本文將為大家介紹完整的腳本構(gòu)建過(guò)程,希望對(duì)大家有所幫助2025-09-09
解決遇到:PytorchStreamReader failed reading zip&n
本文針對(duì)"PytorchStreamReaderfailedreadingziparchive:failedfindingcentral"錯(cuò)誤提出解決方案,包括檢查文件完整性、文件路徑,嘗試更新PyTorch版本,檢查壓縮文件格式,代碼問(wèn)題,或?qū)で蠹夹g(shù)支持等,希望這些經(jīng)驗(yàn)?zāi)芙o遇到同樣問(wèn)題的人一個(gè)參考2024-09-09