computing a interpolating polynomial for some function
显示 更早的评论
I am trying to understand how to use the commands polyfit and polyval so i am trying to interpolate some function, I wrote something so simple but even it is simple it gives me error
x=linspace(-1,1,20);
y= @(x) 1/(1+16*(x^2));
p = polyfit (x,y,2);
plot(x,p)
and this is the error message
Error using polyfit (line 48)
The first two inputs must have the same number of elements.
Error in polyfit_polyval (line 7)
p = polyfit (x,y,2);
采纳的回答
x = linspace(-1,1,20);
yfunc = @(x) 1./(1+16*x.^2);
y = yfunc(x);
...
or
x=linspace(-1,1,20);
y=1./(1+16*x.^2);
...
9 个评论
Anas Gharsa
2022-1-26
Torsten
2022-1-26
Yes, the dot after 1./ and the dot in x.^2.
Look up elementwise multiplication and division.
Anas Gharsa
2022-1-26
Anas Gharsa
2022-1-26
Torsten
2022-1-26
You didn't specify the degree of the fitting polynomial.
Anas Gharsa
2022-1-26
okay my bad !!!!!! but why the graph looks like this
Torsten
2022-1-26
It's the best you can get for a quadratic polynomial as interpolant.
Anas Gharsa
2022-1-26
In your original code, the main problem was that y was a function handle but polyfit requires its first two inputs to be numeric arrays of the same size. Note that Torsten defined the function not as y but as yfunc then evaluated yfunc with x as the input to generate y. So the data that got passed into polyfit as its second input was the result of evaluating that function at the points in x rather than the function.
The use of the element-wise operators allowed the function to be evaluated on an array and return an array with the same size. The matrix form of the operator would require the inputs to be square matrices.
A = [1 2; 3 4]
y = A.*A % each element of y is the square of the corresponding element of A
z = A*A % this perform matrix multiplication instead of element-wise
更多回答（0 个）
类别
在 帮助中心 和 File Exchange 中查找有关 Creating and Concatenating Matrices 的更多信息
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!
Translated by 
