def func1x a b c r =a x0 2+ 2b x1 2 + c 2 printa printx0 printx1 #printx2 return rravel# 生成原始数据xx = npindices10 10z = func1xx 10 5 2 + nprandomnormalsize=100 100abcd para = cur
def func1(x, args): r = 0 for i in range(0, len(args), 3): a = args[i] b = args[i+1] c = args[i+2] r += a * x[i//3] ** 2 + 2b * x[(i//3)+1] ** 2 + c ** 2 return r.ravel()
生成原始数据
xx = np.indices([10, 10]) z = func1(xx, 10, 5, 2) + np.random.normal(size=100) / 100 abcd, para = curve_fit(func1, xx.ravel(), z) print(abcd)
原文地址: https://www.cveoy.top/t/topic/ht6E 著作权归作者所有。请勿转载和采集!