How to implement without function.?

Answers (2)

Jan
Jan on 21 Jan 2013
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.

Asked:

on 21 Jan 2013

Community Treasure Hunt

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

Start Hunting!