我有一個使用 numpy 的 poly1d 定義的多項式:import numpy as npp = np.poly1d([-1.37970733e+07, 8.02055768e+07, -1.98442743e+08, 2.70981281e+08, -2.20611918e+08, 1.07019785e+08, -2.85873914e+07, 3.24251782e+06])然后我對這個多項式關(guān)于我稱為 的點進行級數(shù)展開radMax = 0.80868373。以下是我進行系列擴展的方法:from sympy import series, Symbol, Polyx = Symbol('x')polyExpand = p(x).series(x,radMax,3)/(p(radMax)/0.0523772)這給了我想要的展開式,我想要二階項的系數(shù)。如果運行上面的代碼,您將看到該系數(shù)為“0.2843162879664”。我想提取這個系數(shù)。
1 回答

森林海
TA貢獻2011條經(jīng)驗 獲得超2個贊
如果您想按項處理系數(shù):
from sympy import sympify
print(polyExpand.as_coefficients_dict()[sympify(f"(x - {radMax}) ** 2")])
如果您想通過索引獲取系數(shù):
coeffs = list(polyExpand.as_coefficients_dict().values())
print(coeffs[2]) # order n = 2
這些將正確輸出-0.284316287966400。你忘記了負號。
添加回答
舉報
0/150
提交
取消