| 1 | 1/1 | 返回列表 |
| 查看: 1166 | 回復: 0 | ||
| 【懸賞金幣】回答本帖問題,作者xjzeng將贈送您 10 個金幣 | ||
xjzeng木蟲 (小有名氣)
|
[求助]
python解非線性方程的問題
|
|
|
廢話不多說,上代碼 方法1: from sympy import * t=symbols('t') eq=eq(10**20*exp(-12000/t)/(0.01+10**14*exp(-12000/t))-18700*(t-300)) #eq=[10**20*exp(-12000/t)/(0.01+10**14*exp(-12000/t))-18700*(t-300)]#改成這樣結(jié)果也是一樣的 result=nonlinsolve(eq,t) print(result) 程序運行,返回的是一個conditionset對象,說明不能解這個方程。 方法2: import numpy as np from scipy.optimize import fsolve def f(t): return 10**20*np.exp(-12000/t)/(0.01+10**14*np.exp(-12000/t))-18700*(t-300) result=fsolve(f,[300]) print(result) 運行結(jié)果報錯 minpack.error: result from function call is not a proper array of floats. 方法3 import math from scipy.optimize import fsolve def f(t): return 10**20*math.exp(-12000/t)/(0.01+10**14*math.exp(-12000/t))-18700*(t-300) result=fsolve(f,[300]) print(result) 第3種方法能解出一個根 [303.3],但是方程有3個根。如果把“result=fsolve(f,[300])”這句改成“result=fsolve(f,[320])“,可以得到第二個根 [323.7];如果“result=fsolve(f,[300])”這句改成“result=fsolve(f,[360])“,可以得到第三個根 [349.4]。 如果把“result=fsolve(f,[300])”這句改成“result=fsolve(f,[300,320.360])“,則會報錯 typeerror: only size-1 arrays can be converted to python scalars 各位大佬,如何才能把三個根一次性解出來呢? |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研]
|
簡木ChuFront 2026-03-09 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 一志愿南大化學339分求調(diào)劑,四六級已過,有比賽,有文章 +7 | Gallantzhou 2026-03-07 | 7/350 |
|
|
[考研] 考研調(diào)劑,一志愿山東大學材料與化工,328分,政治51 +5 | 關(guān)你西紅柿929 2026-03-08 | 6/300 |
|
|
[考研] 086000生物與醫(yī)藥319分求調(diào)劑 +4 | Tolkien 2026-03-07 | 4/200 |
|
|
[考研] 0703化學求調(diào)劑 +4 | 很老實人 2026-03-09 | 4/200 |
|
|
[基金申請] 第二個面上項目的獲批幾率如何? +3 | wawm9296 2026-03-04 | 5/250 |
|
|
[考研] 求調(diào)劑,數(shù)一英一274分 +4 | 小菲會努力 2026-03-08 | 4/200 |
|
|
[考研] 一志愿山東大學105500藥學專碩,總分302求調(diào)劑 +5 | 五維天空 2026-03-04 | 11/550 |
|
|
[考研] 0703化學求調(diào)劑,總分309分,一志愿華南師范 +3 | 花與葉@ 2026-03-08 | 3/150 |
|
|
[考研] 理學287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
[考研] 誠求調(diào)劑,323分 有專利/科研/實習經(jīng)歷 +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
[考研] 278求調(diào)劑 +5 | Gale1314 2026-03-06 | 5/250 |
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
[考研] 材料277分求調(diào)劑 +13 | 飯飯星球 2026-03-04 | 14/700 |
|
|
[考研] 求調(diào)劑 +5 | danyyyy 2026-03-04 | 5/250 |
|
|
[考研] 求材料調(diào)劑 +4 | berdmond 2026-03-05 | 4/200 |
|
|
[考研] 316求調(diào)劑 +3 | 林小星發(fā)大財 2026-03-04 | 3/150 |
|
|
[考研] 293求調(diào)劑 +3 | 是樂渝哇 2026-03-04 | 3/150 |
|
|
[考研] 085601 材料305分求助 +4 | 泡泡郵件 2026-03-03 | 6/300 |
|