Fitting an equation with x,y variables and b, d constant.

1 回表示 (過去 30 日間)
ly
ly 2016 年 11 月 17 日
コメント済み: Torsten 2016 年 11 月 21 日
Hi,
I have an equation.
x=[10,50,100,300,500,1000,1500,2000,3000];
y=[0.11,0.17,0.2,0.24,0.29,0.3,0.31,0.35,0.38];
I want to fit this equation and get b and d values.
I tried with lsqcurvefit command, but I can not convert this equation to y=function(x).
Anybody has a suggestion.

採用された回答

Torsten
Torsten 2016 年 11 月 17 日
Use lsqnonlin and define the functions f_i as
f_i = 1/a*(b-ydata(i)).^1.5-log(d./xdata(i))+0.5*log(1-ydata(i)/b)
Best wishes
Torsten.
  5 件のコメント
Walter Roberson
Walter Roberson 2016 年 11 月 21 日
lsqnonlin does not calculate Chi-Square or any other probability measure. It does not create any hypotheses about how well the model fits: it only searches for a minimum.
Torsten
Torsten 2016 年 11 月 21 日
You don't get statistics from lsqnonlin.
You will have to use tools like nlinfit in combination with nlparci and nlpredci.
Best wishes
Torsten.

サインインしてコメントする。

その他の回答 (1 件)

Walter Roberson
Walter Roberson 2016 年 11 月 18 日
a = some value
y1 = @(b, d, x) -b .* (exp(-(2/3) .* lambertw(-3 .* (b.^3 ./ a.^2).^(1/2) .* d.^3 ./ x.^3)) .* d.^2 - x.^2) ./ x.^2
y2 = @(b, d, x) -b .* (exp(-(2/3) .* lambertw(3 .* (b.^3 ./ a.^2).^(1/2) .* d.^3 ./ x.^3)) .* d.^2 - x.^2) ./ x.^2;
guessbd = rand(1,2);
fit1 = fittype(y1, 'coefficients', {'b', 'd'}, 'dependent', 'y', 'independent', x);
fit2 = fittype(y2, 'coefficients', {'b', 'd'}, 'dependent', 'y', 'independent', x);
[bd1, gof1] = fit( x, y, fit1, 'startpoint', guessbd );
[bd2, gof2] = fit( x, y, fit2, 'startpoint', guessbd );
The pair of fits is due to there being two solutions when y is expressed in terms of x, almost identical but differing in sign of the LambertW expression. You would need to check the goodness of fit results to see which was better.
  1 件のコメント
ly
ly 2016 年 11 月 21 日
Ok, I got it.
How to get Chi-Square from lsqnonlin command?

サインインしてコメントする。

カテゴリ

Help Center および File ExchangeLeast Squares についてさらに検索

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by