Question

我想要拟合一个数据集，由n - 维向量data给出（值介于-1和1之间，并且具有相应的x值linspace(0,9,n)），到

给出的参数曲线

x = cos(t)/sqrt(1-(a/b)^2)
y = b*sin(t)-a

代表a < b和0 < t < pi。该曲线是椭圆的上半部分，它在x轴下方a下降，并且在-1和1处与x轴相交。我适合参数a,b。< / p>

在y(x)表单上编写此函数似乎不可能（编辑：不是这样，请参阅注释 - 但是，我仍然想知道以下问题的答案：）。 ..可以从参数形式完成Matlab的拟合吗？怎么样？

感谢。

Answer 1

您可以用y = f(x)格式编写函数，然后使用fitnlm估算非线性模型的参数。

fh = @(x,a,b)b*sin(acos(x*sqrt(1 - (a/b)^2))) - a % Define y = f(x)
x = [-1:0.01:1]; % Generate some x values
yval = fh(x,1,2); % Calculate y values for the x 
yval_ = yval + 0.1*rand(1,201); % Add artifical noise

nonLinMdl = fitnlm(x,yval_,@(p,x)fh(x,p(1),p(2)),[1.5, 2.5])

nonLinMdl = 


Nonlinear regression model:
    y ~ y(x,p1,p2)

Estimated Coefficients:
          Estimate       SE       tStat       pValue  
          ________    ________    ______    __________

    p1    0.71909     0.053201    13.516    5.8209e-30
    p2     1.7668     0.055472    31.849    4.6625e-80


Number of observations: 201, Error degrees of freedom: 199
Root Mean Squared Error: 0.032
R-Squared: 0.988,  Adjusted R-Squared 0.988
F-statistic vs. zero model: 6.36e+04, p-value = 5.87e-280

请注意，计算的参数与我们使用的(1,2)不同。对于不同的参数集，函数形状是相同的这一事实可能需要做些什么。你也可以尝试不同的起点。

ycalc = fh(x,nonLinMdl.Coefficients.Estimate('p1'),nonLinMdl.Coefficients.Estimate('p2'))
plot(x,yval_,x,ycalc)

我无法上传文件，因为我在受限制的环境中，但我绘制了两个系列，模型计算值非常接近输入数据。

是否可以在Matlab中拟合参数曲线？

1 个答案: