这是一个简单的多项式方程:
b^2 + 2b + 1 = 0
我可以很容易地解决这个问题:
import numpy as np
from scipy.optimize import fsolve
eq = lambda b : np.power(b,2) + 2*b + 1
fsolve(eq, np.linspace(0,1,2))
同样,我可以求解任何项数有限的方程。但是如何求解具有无限项的方程,其给出为:
上式可以写成:
(1 - l) * (5.5 + 4.0*l + 4*l^2 + 6*l^3 + 5*l^4 + 5*l^5 + 5*l^6 + 5*l^7 + 5*l^8 + 5*l^9 + 5*l^10 ) = 5
当 n 从 1 变为 10 时。但我想解决这个问题以获得足够大的值,n
使得 LHS ~= RHS。
我知道 LHS 的值,G1 -> Ginf
但不明白我怎么能在这里计算 lambda 的值。
我尝试查看numpy 多项式函数,但找不到与此处相关的函数。