Python curve_fit 不适合正确的罪

如何解决Python curve_fit 不适合正确的罪

我使用普通的 Scipy curve_fit 代码和我自己的 csv 文档值。

from scipy.optimize import curve_fit  
from matplotlib import pyplot as plt 
import csv

x = [];
y = [];
with open('test2.csv') as csv_file:
    csv_reader = csv.reader(csv_file,delimiter=';')
    line_count = 0
    for row in csv_reader:
        v = f"{row[0]}";
        u = f"{row[1]}";
        x.append(int(u));
        y.append(float(v));
        line_count += 1;
    print(f'Processed {line_count} lines.')
  
def test(x,a,b): 
    return a * np.sin(b * x) 
param,param_cov = curve_fit(test,x,y) 
  
  
print("Sine funcion coefficients:") 
print(param) 
print("Covariance of coefficients:") 
print(param_cov)

这些值是使用 Arduino 和加速度传感器测量的。在 Excel 中，他们绘制了一个美丽的正弦。 test2.csv 文件：

1.03;3024
1.28;3027
2.16;3039
3.58;3050
4.04;3053
4.96;3064
5.04;3067
5.03;3072
4.23;3090
4.02;3093
3.88;3097
2.78;3123
2.11;3144
2.02;3147
0.99;3170
0.02;3184
-0.18;3187
-0.91;3198
-1.15;3203
-1.97;3217
-2.12;3220
-3.62;3264
-4.53;3308
-4.56;3312
-4.67;3322
-4.75;3326
-4.97;3338
-5;3341
-5.02;3345
-4.95;3348
-4.65;3362
-4.6;3366
-4.49;3371
-3.96;3385
-3.29;3406
-1.97;3439
-1.6;3450
-0.28;3476

以及 Python 脚本的输出：

Processed 38 lines.
Sine funcion coefficients:
[0.96405843 1.00058573]
Covariance of coefficients:
[[ 7.20387348e-01 -2.68335160e-05]
 [-2.68335160e-05  6.82099708e-08]]

如您所见，系数不是文件中的正弦值。 阅读有什么问题或我做错了什么？

解决方法

1. 您没有为 sin() 函数定义所有自变量。任意 sin() 函数的形式为 Asin(Bx + C)，其中 C 是相位。如果没有定义所有自变量，我认为您不会使用任何方法获得拟合。

2. 您知道 Levenberg–Marquardt 算法（由 scipy.optimize.curve_fit 使用）适用于周期函数吗？我不熟悉它，但曲线拟合器通常不能很好地处理周期函数 - 这就是在这些情况下使用傅立叶变换的原因（我建议您进行调查）。