给定函数关系，求解参数

Posted 2020-12-17 zxixiu

最近忙着搞论文仿真相关的东西，从参考师兄师姐的仿真到运用到自己的仿真里面，有些东西因为仿真材料的不同，相关的参数的也是不同的，只能看理论根据自己学的东西来弄清楚，这些参数的数值了。

这应该算python的一个简单应用吧，根据实验得出的数据，拟合出相对应的参数，虽然比较简单，但是这个感觉还是很棒的，感受到了自己学的东西力量和成就感，棒(?•??•?)??

import numpy as np
from scipy.optimize import curve_fit
def func(x,a,b):
    return a/pow((1+b*np.square(x)),1/3)
x=[0,
40,
80,
120,
160,
200,
240,
280,
320,
360,
400,
440,
480,
520,
560,
600,
640,
680,
720,
760,
800,
840,
880,
920,
960,
1000,
1040,
1080,
]
x=np.array(x)
num=[1632.92791,
1631.30402,
1630.24111,
1624.56044,
1622.82435,
1615.90951,
1611.07326,
1604.99695,
1596.1984,
1586.0594,
1578.11119,
1568.84614,
1556.78799,
1551.49705,
1539.58063,
1526.34147,
1518.16296,
1503.63059,
1493.09595,
1476.38461,
1465.37757,
1454.28195,
1439.80273,
1430.43139,
1418.14295,
1403.0496,
1389.00735,
1376.90197,
]
y=np.array(num)
popt,pcov=curve_fit(func,x,y)
print(popt)
a=popt[0]
b=popt[1]
yvals=func(x,a,b)
print("popt:",popt)
print("a",a)
print("b",b)