二分法求函数极小值画图代码python
时间: 2024-05-15 18:11:35 浏览: 98
二分法是一种基本的数值计算方法,用于求解单调函数的极值。对于一个单峰函数 $f(x)$,可以通过二分法求解其极小值。
下面是使用 Python 画单峰函数的代码,以 $f(x) = x^2\sin(5x)$ 为例:
```python
import numpy as np
import matplotlib.pyplot as plt
def func(x):
return x**2 * np.sin(5*x)
x = np.linspace(-2, 2, 1000)
y = func(x)
plt.plot(x, y)
plt.show()
```
下面是使用 Python 实现二分法求解单峰函数的极小值的代码:
```python
import numpy as np
def func(x):
return x**2 * np.sin(5*x)
def bisection(a, b, tol=1e-6):
"""
二分法求解单峰函数的极小值
:param a: 左端点
:param b: 右端点
:param tol: 精度要求
:return: 极小值的近似解
"""
while (b - a) > tol:
m = (a + b) / 2
if func(m-1e-6) > func(m) < func(m+1e-6):
return m
elif func(m-1e-6) > func(m+1e-6):
a = m
else:
b = m
return (a + b) / 2
result = bisection(-2, 2)
print(result)
```
以上代码实现了一个函数 `bisection`,它接收左右端点 `a` 和 `b`,以及精度要求 `tol`,返回单峰函数在 `[a, b]` 区间内的极小值。
阅读全文