Discussion :

[tubaas]

Bonjour.
c'est peut-être plus un problème de calcul numérique mais je ne trouve pas la solution.
j'ai des données issues d'une fft donc avec des abcisses linéaires mais que je trace avec semilogx.
je m'intéresse alors à l'aire sous la courbe.
pour donner une idée du problème, voilà un exemple simple que je n'arrive pas à résoudre.

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

1
2
3
4
5
6
 
x=(20:2000);
y=zeros(2000-19,1);
y(1:200-19)=1;
y(200-20:2000-19)=-1;
semilogx(x,y)

je cherche donc à définir une fonction somme sur y qui donne 0.
j'ai pensé à pondérer (mettre un poids plus fort pour les x petits) mais je n'ai pas réussi à régler le problème
merci si vous avez des solutions

[Jerome Briot]

Code :

Sélectionner tout - Visualiser dans une fenêtre à part

trapz(log(x),y)

Non ?

[tubaas]

merci
j'étais parti plus sur une méthode des rectangles sans trop m'en rendre compte.
je recherchais le principe pourtant fort simple que la base du rectangle (ou du trapèze) est [log10(x(i+1))-log10(x(i))], facteur que je recherchais et non log(x(i))
en tout cas trapz est plus commode.