| 24小時(shí)熱門(mén)版塊排行榜 |
| 1 | 1/1 | 返回列表 |
| 查看: 1169 | 回復(fù): 0 | |||
| 【懸賞金幣】回答本帖問(wèn)題,作者xjzeng將贈(zèng)送您 10 個(gè)金幣 | |||
xjzeng木蟲(chóng) (小有名氣)
|
[求助]
python解非線性方程的問(wèn)題
|
||
|
廢話不多說(shuō),上代碼 方法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) 程序運(yùn)行,返回的是一個(gè)conditionset對(duì)象,說(shuō)明不能解這個(gè)方程。 方法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) 運(yùn)行結(jié)果報(bào)錯(cuò) 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種方法能解出一個(gè)根 [303.3],但是方程有3個(gè)根。如果把“result=fsolve(f,[300])”這句改成“result=fsolve(f,[320])“,可以得到第二個(gè)根 [323.7];如果“result=fsolve(f,[300])”這句改成“result=fsolve(f,[360])“,可以得到第三個(gè)根 [349.4]。 如果把“result=fsolve(f,[300])”這句改成“result=fsolve(f,[300,320.360])“,則會(huì)報(bào)錯(cuò) typeerror: only size-1 arrays can be converted to python scalars 各位大佬,如何才能把三個(gè)根一次性解出來(lái)呢? |
| 1 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料工程085601調(diào)劑求老師收留 +9 | 強(qiáng)木木木 2026-03-07 | 11/550 |
|
|---|---|---|---|---|
|
[考研] 308求調(diào)劑 +8 | 倘若起風(fēng)了呢 2026-03-05 | 10/500 |
|
|
[考研] 283求調(diào)劑,材料、化工皆可 +5 | 蘇打水7777 2026-03-11 | 5/250 |
|
|
[考研] 求調(diào)劑 +5 | 18880831720 2026-03-11 | 5/250 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +6 | jiaanl 2026-03-11 | 6/300 |
|
|
[考研] 0856化工原理 +5 | z2839474511 2026-03-10 | 5/250 |
|
|
[考研] 工科0856專碩化學(xué)工程269能調(diào)劑嗎 +9 | 我想讀研11 2026-03-10 | 9/450 |
|
|
[考研] 086000生物與醫(yī)藥319分求調(diào)劑 +4 | Tolkien 2026-03-07 | 8/400 |
|
|
[考博] 26申博求助 +3 | 跳躍餅干 2026-03-10 | 4/200 |
|
|
[考研] 材料工程專碩調(diào)劑自薦信,初試323 +6 | 李白26 2026-03-07 | 6/300 |
|
|
[考研] 求調(diào)劑! +3 | 朔朔話 2026-03-09 | 3/150 |
|
|
[考研] 數(shù)二英二309分請(qǐng)求調(diào)劑 +3 | dtdxzxx 2026-03-09 | 4/200 |
|
|
[考研] 0703求調(diào)劑 +10 | c297914 2026-03-06 | 11/550 |
|
|
[考研] 考研求調(diào)劑,工科,總分304 +4 | 矮子不想長(zhǎng)高了 2026-03-08 | 7/350 |
|
|
[考研] 材料化工求調(diào)劑 +7 | 博斯特525 2026-03-06 | 8/400 |
|
|
[考研] 一志愿武漢理工085601,初試301分,請(qǐng)問(wèn)能調(diào)劑到湖北嗎 +3 | 肖yang 2026-03-06 | 3/150 |
|
|
[考研]
|
程晴之 2026-03-06 | 6/300 |
|
|
[考研] 哈爾濱理工大學(xué)2026年研究生調(diào)劑,材料科學(xué)與化學(xué)工程學(xué)院研究生調(diào)劑 +3 | xinliu866 2026-03-06 | 3/150 |
|
|
[考研] 復(fù)試調(diào)劑 +5 | 呼呼?~+123456 2026-03-05 | 5/250 |
|
|
[考研] 紡織、生物、化學(xué)、材料等專業(yè) +3 | Eember. 2026-03-05 | 7/350 |
|