Bonjour, voici mon probleme:
J'ai un vecteur Al et un vecteur Zt. Je dérive Al par rapport à Zt: j'obtiens dAldZ. J'aimerai à partir du vecteur dAldZ retouver le vecteur Al (avec la fonction cutrapz ou autre?).

Voici ma méthode:
Code : Sélectionner tout - Visualiser dans une fenêtre à part
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
Al = [1.9e-2 7.31e-3 1.6e-3 3.23e-4 6.68e-5 1.65e-5 5.22e-6 1.48e-6 3.48e-7];
Zt = 0:10:80;
 
% interpolation sur une grille Z_in plus fine:
Z_in=0:80;
Al_in = interp1(Zt,Al,Z_in,'pchip')
 
% calcul de la derivee:
dAldZ=diff(Al_in)./diff(Z_in);
 
% cumtrapz:
Al2=cumtrapz(Z_in(1:end-1),dAldZ);
 
% je trace pour verifier:
plot(Zt,Al,'o',Z_in,Al_in,'-',Z_in(1:end-1),Al2,'r')
Et ça ne marche pas! je ne omprends pas mon erreur! Merci de m'aider