En analyse numérique, la méthode des différences finies est une technique courante de recherche de solutions

approchées d'équations aux dérivées partielles qui consiste à résoudre un système de relations (schéma
numérique) liant les valeurs des fonctions inconnues en certains points suffisamment proches les uns des autres.
On va utiliser cette méthode pour calculer les dérivées d’ordre N numériquement.

En utilisant la méthode des différences finis on va montrer que la dérivée en un point xi s’exprime en fonction
des points voisins xi 1 et xi 1

Une discrétisation des opérateurs différentiels (dérivées premières, secondes, ………..) peut être obtenue par
les formules de Taylor. La formulation de Taylor-Young est préférable dans son utilisation simple

h² h3 hn
f  x  h   f ( x)  f '( x)h  f ''( x) f '''( x)  .......  f ( n ) ( x)  O(h N 1 ) Mais il est préférable d’utiliser
2! 3! n!
''  dx  ² '''  dx  ( n )  dx 
3 n

la notation : y  x  dx   yi  yi dx  yi  yi  .......  yi
'

2! 3! n!
1) Dérivée Première :

Dérivée première :

yi 1  yi
A droite : y ( xi  dx)  yi 1  yi 1  yi dx 
 yi 
' '

dx
yi  yi 1
A gauche : y ( xi  dx)  yi 1  yi 1  yi dx 
 yi  ' '

dx
dx ² 
I ) y ( xi  dx)  yi 1  yi  yi' dx  yi''
2!  I  II y y
Centrée :    yi'  i 1 i 1
dx ²  2dx
II ) y ( xi  dx)  yi 1  yi  yi' dx  yi''
2! 
2) dérivée seconde :

A droite :

dx 2 dx 2
y( xi  dx)  yi 1  yi  yii dx  yi''  yi  yii dx  yi''
2! 2

4dx 2
y( xi  2dx)  yi  2  yi  2 y dx  y
i
i  yi  2 yii dx  2 yi''dx 2
''
i
2!
yi  2  2 yi 1  yi
yi  2  2 yi 1   yi  dx 2 yi'' 
 yi'' 
dx 2
A gauche :

dx 2 '' dx
2
y ( xi  dx)  yi 1  yi  y dx  y
i
i
''
i  yi  yi dx  yi
i

2! 2
2
4dx
y ( xi  2dx)  yi  2  yi  2 yii dx  yi''  yi  2 yii dx  2 yi''dx 2
2!
y  2 yi 1  yi
yi  2  2 yi 1   yi  dx 2 yi'' 
 yi''  i  2
dx 2
 dx 2 dx3 
I ) yi 1  yi  yi' dx  yi''  yi'''
2 9  I  II y  y  2 yi
   yi''  i 1 i 1
Centrée : dx 2 3
dx  dx ²
II ) yi 1  yi  yi' dx  yi''  yi'''
2 9 

3 ) Dérivée tiers :

A droite on a :

dx 2 ''' dx
3
y (x  dx)  yi 1  yi  y dx  y
'
i
''
i  yi
2 9
2
4dx 4dx3
y (x  2dx)  yi  2  yi  2 yi' dx  yi''  yi'''
2 3
2
9dx 9dx3
y (x  3dx)  yi 3  yi  3 yi' dx  yi''  yi'''
2 2
 yi  3 yi 1  3 yi  2  3 yi 3
A partir de ces 3 relations on obtient : yi 
'''

dx3
A gauche on a :

dx 2 dx3
y ( x  dx)  yi 1  yi  yi' dx  yi''  yi'''
2 9
4
y (x  2dx)  yi  2  yi  2 yi'  2 yi''dx 2  dx 3 yi'''
3
9 9
y (x  2dx)  yi 3  yi  3 yi' dx  yi'' dx 2  yi'''dx 3
2 2
yi  3 yi 1  3 yi  2  yi 3
A partir de ces 3 relations on obtient : yi 
'''

dx3
yi 3 yi  2 yi 1 yi yi 1 yi  2 yi 3
'
ygauche 1 -1

1 -1
'
ydroite
''
ygauche -1 -2 1

1 -2 -1
''
ydroite
'''
ygauche 1 -3 3 -1

1 -3 3 -1
'''
ydroite

x=1:0.5:6;
y=[9;15.375;26;41.625;63;90.875;126;169.125;221;282.375;354];
d=size(y);
N=d(1);
dx=0.5;
% deriveé premier yp
for i=1:N;
if (i==1)
yp(i)= (y(i+1)-y(i))/dx ;
else
yp(i)= (y(i)-y(i-1))/dx ;
end
end
%deriveé seconde ys
for i=1:N;
if (i<3)
ys(i)= (y(i)-2*y(i+1)+y(i+2))/(dx^2) ;
else
ys(i)= (y(i-2)-2*y(i-1)+y(i))/(dx^2) ;
end
end
%deriveé 3 yt
for i=1:N;
if (i<4)
yt(i)= (-y(i)+3*y(i+1)-3*y(i+2)+y(i+3))/(dx^3) ;
else
yt(i)= (-y(i-3)+3*y(i-2)-3*y(i-1)+y(i))/(dx^3) ;
end
end
y