python解非线性方程的问题

2020-10-16 01:00:53 646 阅读 AI 写同款 · GEO 优化AI 选题数据监测一键分发 ›

用CentOS 7安装cadence搭建适合IC Design的科研环境（四）——IC617、MMSIM151、calibre2015安装过程step by step 以下操作都在root权限下执行。建立安装目录建议把这种大型软件装在/opt目录下。mkdir -p /opt/cadence/installscape mkdir -p /opt/mentor/Calibre2015其中-p参数是直接建立父目录。因为IC617和MMISM151会自己建立安装目录，所以在这里没有建相关目录。而calibre则是默认安装在安装包所在的目录下，所以需要提前建立。解压安阅读详情

废话不多说，上代码
方法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)]#改成这样结果也是一样的
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)

运行结果报错
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

各位大佬，如何才能把三个根一次性解出来呢？

...全文