在matlab中将曲线拟合到具有多个不同系数的一组数据

时间:2016-04-17 20:26:52

标签: matlab curve-fitting

我需要在一组约100个数据点上插入一行,这组数据遵循Pacejka公式,如下所示:

Fy = Dy sin [Cy arctan {By x - Ey(By x - arctan(Byx))}] + Svy

其中Dy,Cy,By,Ey和Svy是要求解的系数。

我可以用这段代码绘制一条线,但它不是靠近数据的地方。

这是我到目前为止拟合公式的方法。如何将其更改为更接近实际的最佳拟合线?

x = SA;
y = Fy;
expr = 'D * sin(C * atan(B*x - E*(B* x - atan(B*x)))) + A';
ft = fittype(expr, 'independent', 'x');
opts = fitoptions('Method', 'NonlinearLeastSquares');
opts.StartPoint = ones(1,5);
[fitresult, gof] = fit(x, y, ft, opts)
plot(fitresult, x, y)

这是我的代码现在返回的内容 enter image description here enter image description here

1 个答案:

答案 0 :(得分:3)

就像你一样,只记下那个公式:

expr = 'D * sin(C * atan(B*x - E*(B* x - atan(B*x)))) + A';
ft = fittype(expr, 'independent', 'x');
opts = fitoptions('Method', 'NonlinearLeastSquares');
opts.StartPoint = ones(1,5);  % [A,B,C,D,E]
[fitresult, gof] = fit(x, y, ft, opts)
plot(fitresult, x, y)

编辑:

基于Wikipedia article,这是一个小例子:

% some data based on equation
B = 0.714;
C = 1.4;
D = 800;
E = -0.2;
f = @(x) D * sin(C * atan(B*(1-E)*x + E*atan(B*x)));
x = linspace(0,10,200)';  %'
y = f(x);

% add noise
yy = y + randn(size(x))*16;

% fit
%expr = 'D * sin(C * atan(B*(1-E)*x + E*atan(B*x)))';
expr = 'D * sin(C * atan(B*x - E*(B* x - atan(B*x))))';
ft = fittype(expr, 'independent', 'x', 'dependent','y');
opts = fitoptions('Method', 'NonlinearLeastSquares');
opts.StartPoint = [1 1 1000 1];  % [B C D E]
[fitresult, gof] = fit(x, y, ft, opts)

% plot
yhat = feval(fitresult, x);
h = plot(x,y,'b-', x,yhat,'r-', x,yy,'g.');
set(h, 'LineWidth',2)
legend({'y', 'yhat', 'y+noise'}, 'Location','SouthEast')
grid on

fit

结果:

fitresult = 
     General model:
     fitresult(x) = D * sin(C * atan(B*x - E*(B* x - atan(B*x))))
     Coefficients (with 95% confidence bounds):
       B =      0.5916  (0.5269, 0.6563)
       C =       1.899  (1.71, 2.089)
       D =       783.7  (770.5, 796.9)
       E =       1.172  (1.136, 1.207)
gof = 
           sse: 6.6568e+04
       rsquare: 0.9834
           dfe: 196
    adjrsquare: 0.9832
          rmse: 18.4291

请注意,使用像[1 1 1 1]这样的起点远非真正的解决方案,因此拟合将在不收敛的情况下停止(D参数的比例与其他参数B非常不同,CE)...相反,我必须从靠近[1 1 1000 1]的地方开始。