码迷,mamicode.com
首页 > 其他好文 > 详细

给定函数关系,求解参数

时间:2020-07-12 18:48:59      阅读:49      评论:0      收藏:0      [点我收藏+]

标签:相对   imp   UNC   array   运用   code   相关   print   简单   

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

这应该算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)

 

给定函数关系,求解参数

标签:相对   imp   UNC   array   运用   code   相关   print   简单   

原文地址:https://www.cnblogs.com/zxixiu/p/13289057.html

(0)
(0)
   
举报
评论 一句话评论(0
登录后才能评论!
© 2014 mamicode.com 版权所有  联系我们:gaon5@hotmail.com
迷上了代码!