要估計函數的局部極值,可以使用SciPy庫中的optimize模塊中的minimize函數。下面是一個簡單的示例:
import numpy as np
from scipy.optimize import minimize
# 定義一個函數
def func(x):
return x**2 + 2*x + 1
# 使用minimize函數找到函數的局部極小值
result = minimize(func, x0=0)
print(result.x) # 打印局部極小值的位置
print(result.fun) # 打印局部極小值的函數值
在這個示例中,我們定義了一個簡單的二次函數func(x) = x**2 + 2*x + 1
,然后使用minimize函數找到該函數的局部極小值。x0
參數指定了搜索的起始點,默認為0。最后,我們打印出找到的局部極小值的位置和函數值。