Numerical derivative central difference
4 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Arctgx
le 24 Mar 2020
Réponse apportée : Ameer Hamza
le 25 Mar 2020
Hello, i try to write a code about central difference. But something is wrong. My code is below. What is wrong?
function derivativecentral(numericalvalue, step, n)
centraltable=[0 -1 0 1 0;0 1 -2 1 0;-1 2 0 -2 1;1 -4 6 -4 1];
x=numericalvalue-2*step: step:numericalvalue+2*step;
m=2;
if (n/2)==floor(n/2); m=1; coeff = centraltable(n,:);
der=sum(coeff.*fonk(x)) / (m*step^n);
display(der);
end
end
function [f] = fonk(x) f=x*exp(x); end
At the end of code, i am writing derivativecentral(0,0.1,1)
But there is nothing.
0 commentaires
Réponse acceptée
Ameer Hamza
le 25 Mar 2020
The code displays nothing because the condition did not become true. Try the following code with a different input
derivativecentral(0,0.1,4)
function derivativecentral(numericalvalue, step, n)
centraltable=[0 -1 0 1 0;0 1 -2 1 0;-1 2 0 -2 1;1 -4 6 -4 1];
x=numericalvalue-2*step: step:numericalvalue+2*step;
m=2;
if (n/2)==floor(n/2); m=1; coeff = centraltable(n,:);
der=sum(coeff.*fonk(x)) / (m*step^n);
display(der);
end
end
function [f] = fonk(x)
f=x.*exp(x); % <---- this is syntax for element-wise multiplication
end
0 commentaires
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Simscape Multibody 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!