阻尼振动曲线拟合python
时间: 2023-12-06 10:38:16 浏览: 255
python-曲线拟合-原理-代码.doc
为了拟合阻尼振动曲线,可以使用Python中的SciPy库。具体步骤如下:
1. 导入必要的库
```python
import numpy as np
from scipy.optimize import curve_fit
import matplotlib.pyplot as plt
```
2. 定义阻尼振动函数
```python
def func(x, a, b, c, d):
return a * np.exp(-b * x) * np.cos(c * x + d)
```
3. 准备数据
```python
x = np.linspace(0, 10, 100)
y = 2.5 * np.exp(-0.3 * x) * np.cos(1.2 * x + 0.5) + np.random.normal(size=100)
```
4. 拟合数据
```python
popt, pcov = curve_fit(func, x, y)
```
5. 绘制拟合曲线
```python
plt.plot(x, y, 'b-', label='data')
plt.plot(x, func(x, *popt), 'r-', label='fit')
plt.legend()
plt.show()
```
其中,`popt`是拟合函数的参数,`pcov`是协方差矩阵。
阅读全文