computing a interpolating polynomial for some function

4 vues (au cours des 30 derniers jours)
Anas Gharsa
Anas Gharsa le 26 Jan 2022
Commenté : Steven Lord le 26 Jan 2022
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);

Réponse acceptée

Torsten
Torsten le 26 Jan 2022
Modifié(e) : Torsten le 26 Jan 2022
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 commentaires
Anas Gharsa
Anas Gharsa le 26 Jan 2022
i understand now!! Thank you once again
Steven Lord
Steven Lord le 26 Jan 2022
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]
A = 2×2
1 2 3 4
y = A.*A % each element of y is the square of the corresponding element of A
y = 2×2
1 4 9 16
z = A*A % this perform matrix multiplication instead of element-wise
z = 2×2
7 10 15 22

Connectez-vous pour commenter.

Plus de réponses (0)

Catégories

En savoir plus sur Interpolation dans Help Center et File Exchange

Produits


Version

R2021b

Community Treasure Hunt

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

Start Hunting!

Translated by