FindFit用于隐式函数 - mathematica

时间:2017-04-05 13:19:17

标签: wolfram-mathematica implicit test-data

我想将一些测试数据放入一些隐式函数中。

我想将一些参数拟合到椭圆方程f(x,y)= a其中a是已知变量。我的测试数据和函数更复杂,但是我得到的数据点多于变量。有点不可能将我想要拟合的方程式转换为像f(x)= y这样的显式形式。因此我附加了一些代码来获得基本思想。

Test = {{0, 1}, {0.1, 0.9}, {1.1, 0}};

Ftest = a*x^2 + b*y^2

FindFit[Test, Ftest == 2, {a, b}, {x, y}];

但是这会导致错误:坐标数(1)不等于变量数\ (2)。 >>

1 个答案:

答案 0 :(得分:2)

您可以将其设为最小二乘最小化:

data = {{0, 1}, {0.1, 0.9}, {1.1, 0}}
Ftest[x_, y_] := a*x^2 + b*y^2
fit = FindMinimum[ Total[(Ftest @@@ data - 2)^2] , {a, b}] 
ContourPlot[ (Ftest[x, y] /. fit[[2]]) == 2 , {x, 0, 1.5}, {y, 0, 
  1.5}, Epilog -> {Red, Point /@ data}]

enter image description here

使用你需要解决的拟合函数,最终得到:

fit = NonlinearModelFit[data, Sqrt[2 - a*x^2]/Sqrt[b], {a, b}, x]

Plot[fit[x], {x, 0, 1.2}, Epilog -> {Red, Point /@ data}, 
 AspectRatio -> 1] 

enter image description here