How to implement without function.?

1 次查看(过去 30 天)
Lalit Patil
Lalit Patil 2013-1-21
polyfit(x,y,1)
How to implement the same thing without function..?

回答(2 个)

Walter Roberson
Walter Roberson 2013-1-21

Jan
Jan 2013-1-21
I assume that this is a homework question. There are good reasons to look for a faster implementation than Matlab's polyfit, but you do not ask for speed.
At first you need to construct the Vandermonde matrix. There are a lot of different algorithms to get it, and Google or WikiPedia know more details.
At second the least squares problem must be solved. V \ y with the Vandermonde matrix V and the y values is sufficient. A QR decomposition is more stable, so I suggest to try it and ask again in case or problems.
Be sure to use properly scaled inputs for x: A mean of 0 and a standard deviation of 1 is a fair choice to support the numerical stability of the algorithm. It is easy to consider the scaling afterwards.

类别

Help CenterFile Exchange 中查找有关 Mathematics 的更多信息

Community Treasure Hunt

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

Start Hunting!

Translated by