如何计算具有某些已知系数的POLY_FIT?

时间:2015-11-08 20:19:57

标签: data-fitting idl-programming-language

让我们举个例子:

我们使用对应于已知多项式f(x)= 0.25-x + x2的X和Y数据。使用POLY_FIT计算二次多项式拟合返回精确系数(在机器精度范围内)。

; Define an 11-element vector of independent variable data:
X = [0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0]
; Define an 11-element vector of dependent variable data:
Y = [0.25, 0.16, 0.09, 0.04, 0.01, 0.00, 0.01, 0.04, 0.09, $
   0.16, 0.25]
; Define a vector of measurement errors: 
measure_errors = REPLICATE(0.01, 11)
; Compute the second degree polynomial fit to the data:
result = POLY_FIT(X, Y, 2, MEASURE_ERRORS=measure_errors, $
   SIGMA=sigma)
; Print the coefficients:
PRINT, 'Coefficients: ', result
PRINT, 'Standard errors: ', sigma

此示例打印,

Coefficients:    0.250000    -1.00000    1.00000

Standard errors:    0.00761853    0.0354459    0.0341395

按预期工作,但假设我已经知道系数0.25,我怎样才能将该系数传递给POLY_FIT?或者我可能必须使用其他一些FIT功能?

来源:http://www.exelisvis.com/docs/POLY_FIT.html

1 个答案:

答案 0 :(得分:1)

你看过MPFIT了吗?它有关键字来传递有关参数的信息。 Docs