Bonjour à tous !
J'ai deux courbes :
y1<=>T(z) = -2.5000e-07 *z^2 +0.0334 z
y2<=>T(z) -5.0000e-07*z^2 + 0.0584 z
Je dois faire une approximation du gradient de la température dT/dz en fonction de z avec la méthode des différences finies pour trouver une dérivée mais je dois avoir un soucis dans mon code :
Pour la première dérivée:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8 h=1000 X=0:h:50000 f = y1(X) ydérivée=diff(f)/h f2=y2(X) ydérivée2=diff(f2)/h plot(X,ydérivée,'b',X,ydérivée2,'g')
Et pour la dérivée seconde
Mais ce que je rencontre depuis tout à l'heure ce n'est que des messages d'erreurs^^
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4 yseconde1=diff(ydérivée)/h yseconde2=diff(ydérivée2)/h plot(X,yseconde1,'k',X,yseconde2,'r')
si qqn pouvait m'aider svp
Partager