polyfit3

Fit polynomial to data, with new features as forced to zero coefficients and data weighting.
3.1K Downloads
Updated 16 Dec 2010

View License

POLYFIT3(X,Y,N,NUL,W) finds the coefficients of a polynomial P(X) of degree N that fits the data, P(X(I))~=Y(I), in a least-squares sense. Any of the coefficients can be forced to be zero, and data can be weighted.

NUL is a vector with coefficients forced to be zero in usual order.

W is a vector containing the weights of data.

POLYFIT3(X,Y,N,NUL) assumes that all weights of data are equal to 1.
POLYFIT3(X,Y,N,[],W) allows to weight data without forcing any coefficient to be zero.

The original usage is still available.

Tested under MATLAB 6.0.0.

DEMO (for not weighted data).
x=linspace(0,11,12);
noise=(rand(size(x))-0.5);
y=3-x.^3+2*x.^5+noise;
[p1,S] = polyfit3(x,y,5)
w=ones(size(x));
[p2,S] = polyfit3(x,y,5,[2 4 5],w)
close all;
figure;
plot(x,y,'.r'); hold on;
x=linspace(0,11);
plot(x,polyval(p1,x),'b');
plot(x,polyval(p2,x),'k');

DEMO (For weighted data).
x=linspace(-5,5,11);
noise=(rand(size(x))-0.5);
y=x.^4+noise; y(11)=y(11)+100;
[p1,S] = polyfit3(x,y,4)
w=ones(size(x)); w(11)=0.1;
[p2,S] = polyfit3(x,y,4,[2 3 4 5],w)
close all;
figure;
plot(x,y,'.r'); hold on;
x=linspace(-5,5,11);
plot(x,polyval(p1,x),'b');
plot(x,polyval(p2,x),'k');

Cite As

Antoni J. Canós (2024). polyfit3 (https://www.mathworks.com/matlabcentral/fileexchange/4262-polyfit3), MATLAB Central File Exchange. Retrieved .

MATLAB Release Compatibility
Created with R12
Compatible with any release
Platform Compatibility
Windows macOS Linux
Categories
Find more on Elementary Math in Help Center and MATLAB Answers

Community Treasure Hunt

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

Start Hunting!
Version Published Release Notes
1.1.0.0

Error fixed in the calculation of S.normr and S.R. Now p and structure S created by POLYFIT3 can be used by polyval and prediction error estimates DELTA is correctly obtained.

1.0.0.0