Using polyfit to fit power function including the initial point x=0

Question

Laura 2014-10-9

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/158077-using-polyfit-to-fit-power-function-including-the-initial-point-x-0

评论： Matt J 2014-10-10

I want to find the equation for data points that are given below

 x= [ 0 0.0005 0.001 0.005 0.01 0.05 0.1];
 y=[0.43 0.47 0.51 0.77 1 1.9 2.44];

Because the first x component is zero hence I could not include the first data set if I fit it as power function.

Is there a way to fit it as y=y(0) +a*x^m?

Thank you.

1 个评论
显示 -1更早的评论隐藏 -1更早的评论

Matt J 2014-10-10

Is m unknown? Does it have to be integer-valued?

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Star Strider 2014-10-9

2
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/158077-using-polyfit-to-fit-power-function-including-the-initial-point-x-0#answer_154618

在 MATLAB Online 中打开

You can do a power fit with fminsearch:

   x= [ 0 0.0005 0.001 0.005 0.01 0.05 0.1];
   y=[0.43 0.47 0.51 0.77 1 1.9 2.44];
   xc = linspace(min(x),max(x));
f = @(b,x) b(1).*x.^b(2);
SSE = @(b,f,x,y) sum((y-f(b,x)).^2);
B = fminsearch(@(b) SSE(b,f,x,y),  [1;1]);
figure(1)
plot(x, y, 'pr')
hold on
plot(xc, f(B,xc), '-g')
hold off
grid
legend('Data', 'Power Function Fit', 'Location', 'NW')
text(0.05, 1.2, sprintf('\\itf\\rm(\\itx\\rm) = %.1f \\itx\\rm^{%.2f}',B))

producing:

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 2

the cyclist 2014-10-9

1
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/158077-using-polyfit-to-fit-power-function-including-the-initial-point-x-0#answer_154615

在 MATLAB Online 中打开

I assume that you mean you want to fit

y == y0 + a * x^m

and you want to estimate values for y0, a, and m.

If you have the Statistics Toolbox, you can use the nlinfit function to do that sort of fit.

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

Answer 3

Chad Greene 2014-10-9

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/158077-using-polyfit-to-fit-power-function-including-the-initial-point-x-0#answer_154614

在 MATLAB Online 中打开

x=[0 0.0005 0.001 0.005 0.01 0.05 0.1];
y=[0.43 0.47 0.51 0.77 1 1.9 2.44];
plot(x,y,'ko')
hold on
linearFit = polyfit(x,y,1); 
xfit = linspace(0,.1,100);
yfitLinear = linearFit(1)*xfit + linearFit(2); 
plot(xfit,yfitLinear,'b')
secondOrderFit = polyfit(x,y,2); 
yfitSecondOrder = secondOrderFit(1)*xfit.^2 + secondOrderFit(2)*xfit + secondOrderFit(3); 
plot(xfit,yfitSecondOrder,'r')
legend('original data','linear fit','second order fit','location','southeast')
legend boxoff
box off