# computing a interpolating polynomial for some function

3 visualizzazioni (ultimi 30 giorni)
Anas Gharsa il 26 Gen 2022
Commentato: Steven Lord il 26 Gen 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);
##### 0 CommentiMostra -2 commenti meno recentiNascondi -2 commenti meno recenti

Accedi per commentare.

### Risposta accettata

Torsten il 26 Gen 2022
Modificato: Torsten il 26 Gen 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 CommentiMostra 7 commenti meno recentiNascondi 7 commenti meno recenti
Anas Gharsa il 26 Gen 2022
i understand now!! Thank you once again
Steven Lord il 26 Gen 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

Accedi per commentare.

### Categorie

Scopri di più su Polynomials in Help Center e File Exchange

R2021b

### Community Treasure Hunt

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

Start Hunting!

Translated by