Fit a polynomial function
조회 수: 10 (최근 30일)
이전 댓글 표시
Does someone know how it is possible to fit a polynomial function whent the x value is a vector? In other words, if we want to fit a polynomial function with output data y and input parameters x where x=[x1,x2,x3,....,xn]. Because until now the only thing that I have found is only if x is a single parameter.
댓글 수: 0
채택된 답변
the cyclist
2015년 6월 21일
편집: the cyclist
2015년 6월 21일
The easiest way to do this, if you have the Statistics Toolbox, is to use the fitlm function. Here is some example code:
load carsmall
tbl = table(Weight,Acceleration,MPG,'VariableNames',{'Weight','Acceleration','MPG'});
% % The formula from the example
% formula = 'MPG ~ Weight + Acceleration';
% A formula with 2nd- and 3rd-order polynomials
formula = 'MPG ~ Weight^3 + Acceleration^2';
lm = fitlm(tbl,formula)
A few things to notice:
This code is based on the example in that documentation.
I have included -- but commented out -- the formula in which Weight and Acceleration are included to first order.
The model I fit includes Acceleration to 2rd order and Weight to 3rd order. Notice that the way I specified the model, MATLAB automatically included the lower-order terms (including the intercept).
You could also include cross terms like Acceleration*Weight, but I did not.
Remember that a "linear" model is one that is linear in the coefficients. Even though this model includes terms like Acceleration^2, it is still linear because the coefficient of that term will be linear.
댓글 수: 6
dpb
2015년 6월 21일
But, again, ML has to have some way to "know" what you intend and with multiple variables the possibilities are endless so there's no attempt to try to generalize beyond some lower-order special cases that are fairly frequent but you kept dancing around the question of whether you did or did not, think interaction terms would be important for starters...
추가 답변 (2개)
Image Analyst
2015년 6월 21일
See my polyfit demo, attached below this image it creates
댓글 수: 4
Image Analyst
2015년 6월 21일
Attached is an example. Consider x1 to be the horizontal direction, and x2 to be the orthogonal (vertical) direction. It fits the data (models it) to a 4th order polynomial in both directions. For each (x1, x2) pair, I have a value f(x1,x2) which is the intensity of the image. Then I fit a 2D 4th order polynomial surface to those values.
dpb
2015년 6월 21일
There are several regression and curve fitting routines if you have the Statistics and/or Curve Fitting toolboxes; if you don't you can use the "backslash" operator that will do a least squares solution to an overdetermined system. In this case you write the explicit model as the design matrix (not forgetting to include the column of Ones for the intercept term, of course) and a vector of the observation values as
X=[ones(length(y),1) x1 x2 ... xN];
and solve as
c=X\y;
where each xi is the (column) vector of the values of the associated independent variable and y is your vector of observations.
For example, if you were to try to fit a model of the form z=f(x,y) with the cross term, you could write
X=[ones(size(z)) x x.*y y];
c=X\z;
Again, the above assumes column vectors x,y,z are the two independent variables and z is the observation. c will be the coefficients of the the model
z=c(1) + c(2)*x + c(3)*x.*y + c(4)*y;
댓글 수: 0
참고 항목
카테고리
Help Center 및 File Exchange에서 Polynomials에 대해 자세히 알아보기
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!