How to numerically differentiate provided data?
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Lingbai Ren
le 15 Sep 2021
Commenté : Lingbai Ren
le 24 Sep 2021
I have the measurements of x with corresponding y displacement lengths:
x = [0,0.375,0.75,1.125,1.5,1.875,2.25,2.625,3];
y = [0,-0.2571,-0.9484,-1.9689,-3.2262,-4.6414,-6.1503,-7.7051,-9.275];
and dy/dx = theta(x) >> theta is the slope
My question is how to numerically differentiate the provided data for displacement y(x)
and the hint is I can choose formulas of any error order.
It's a question combined both math and numerical method, can any one give some help?
0 commentaires
Réponse acceptée
Walter Roberson
le 15 Sep 2021
x = [0,0.375,0.75,1.125,1.5,1.875,2.25,2.625,3];
y = [0,-0.2571,-0.9484,-1.9689,-3.2262,-4.6414,-6.1503,-7.7051,-9.275];
theta = gradient(x,y)
plot(x, y, x, theta)
legend({'x', 'theta'})
3 commentaires
Walter Roberson
le 15 Sep 2021
x = [0,0.375,0.75,1.125,1.5,1.875,2.25,2.625,3];
y = [0,-0.2571,-0.9484,-1.9689,-3.2262,-4.6414,-6.1503,-7.7051,-9.275];
orders = 1:7;
for order = orders
p = polyfit(x, y, order);
predict = polyval(p, x);
plot(x, predict, 'displayname', "order = " + order);
hold on
err(order-min(orders)+1) = sum((predict - y).^2);
end
xlabel('x'); ylabel('y predicted')
legend show
figure(2)
plot(orders, err);
xlabel('polynomial order'); ylabel('sse')
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Loops and Conditional Statements dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!