Curve fitting with fails misarebly
9 vues (au cours des 30 derniers jours)
I'm trying to fit a curve in matlab using fit in command line. The input data is:
And the target function is
Y = 1/(kappa*X.^a)
I am using fittype, fitoptions, and fit as follow:
model1 = fittype('1/(kappa*x.^pow)');
opt1 = fitoptions(model1);
opt1.StartPoint = [1e-5 -2];
[fit1,gof1] = fit(X,Y.^-1,model1,opt1)
I get results with rsquare of roughly -450 which are vaguely in the same direction as the measurement.
A plot of Y as a function of X in a log log scale, yield an almost linear plot. The resulting fit however, has a smaller slope (closer to 0) and is higher than any of the data points (almost twice as high as highest data point). It is needless to say that the line doesn't cross any of the data point.
How can I generate better quality of data fitting?
When I try to fit the results to
I get a good result with pow being roughly -1.5.
But if I try to fit to:
I get a poor result. Why is that so?
Richard Willey le 9 Jan 2012
Curve Fitting Toolbox includes a lot of nice code that will automatically choose "good" starting points for nonlinear regressions for supported model types. You're much better off using the library of standard models rather than hard coding your own model.
One of the nice things about cftool is that you can use the tool to generate a fit and then (automatically) create a function that with replicate this analysis without any need to open the GUI.
I attached some code that I auto generated from cftool.
The function had no trouble generating a good fit. I used a "power" model which looks to be the same as your custom equation.
function [fitresult, gof] = createFit1(X, Y)
% Create a fit.
% Data for 'untitled fit 1' fit:
% X Input : X
% Y Output: Y
% fitresult : a fit object representing the fit.
% gof : structure with goodness-of fit info.
% See also FIT, CFIT, SFIT.
% Auto-generated by MATLAB on 09-Jan-2012 09:08:03
%%Fit: 'untitled fit 1'.
[xData, yData] = prepareCurveData( X, Y );
% Set up fittype and options.
ft = fittype( 'power1' );
opts = fitoptions( ft );
opts.Display = 'Off';
opts.Lower = [-Inf -Inf];
opts.StartPoint = [1.92785231569241e-005 -1.5782710696971];
opts.Upper = [Inf Inf];
% Fit model to data.
[fitresult, gof] = fit( xData, yData, ft, opts );
% Plot fit with data.
figure( 'Name', 'untitled fit 1' );
h = plot( fitresult, xData, yData );
legend( h, 'Y vs. X', 'untitled fit 1', 'Location', 'NorthEast' );
% Label axes
xlabel( 'X' );
ylabel( 'Y' );
Andrew Newell le 8 Jan 2012
I got a nice looking fit using cftool with the type of fit chosen to be "Power".
EDIT: I think your mistake is using Y.^-1 in the last line. Try this code:
model1 = fittype(@(kappa,pow,x) kappa*x.^pow);
opt1 = fitoptions(model1);
opt1.StartPoint = [2e-5 -1.5];
[fit1,gof1] = fit(X,Y,model1,opt1)