Discussion :

[adnino]

Bonjour,
comment pourrais-je calculer l'air d'un signal qui se présente ainsi sachant que temps d’incrémentation est de (tinc 1e-9 s)
Merci à avous

[shaiHulud]

Bonjour,
Via une intégrale. Voir les fonctions trapz, integral.

[adnino]

j'ai essayé de faire ça

Code :

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

q = integral(S_inv,0,(tmax/tinc))

mais ça donne l'erreur

Error using integral (line 82)
First input argument must be a function handle.

merci

[shaiHulud]

integral() est faite dans le cas ou tu as une vraie fonction. Toi, tu as juste une grille de valeurs, utilise trapz.

[adnino]

Ca fera aussi une intégrale rectangulaire aussi ?

Trapz fait ce calcul
Elle fait l'integral de chaque valeur du vecteur du signal X le pas d'incrementation ?

[adnino]

Parfait,
Un grand merci
Et si je veux calculer l’énergie de ce signal ?

[shaiHulud]

Il faut que tu me rappelles la définition de l'énérgie

[adnino]

E=∫x(t)²dt % integrale de 0 à tmax%
Merci

[shaiHulud]

bah il suffit de passer les valeurs au carré à trapz -> trapz (S_inv.^2)

[adnino]

J'ai pris un exemple de test mais je ne trouve pas la même valeur,

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
t=(0:10);
 
y=t.*2;
plot(t,y);
 
Energie=trapz (t,y.^2) 
----% resultat matlab---%
Energie =
 
        1340
-----%résultat par calcul m'a donné
111.11

je ne trouve pas la même valeur

[Invité]

Bonjour,

Si tu calcules : y=t.*2; (déjà fois 2 et non exposant ?)
Puis prends le carré de ceci trapz (t,y.^2) = trapz(t,(t.*2).^2), tu ne risques pas de trouver la bonne valeur

Ensuite, théoriquement on a
Mais tu veux je présume calculer l'aire de ta courbe non ?

[adnino]

bonjour,
j'ai fait ce petit programme
Je pense que je suis arrivé mais il y a une difference importante entre les deux résultat,
Il y aurait pas une autre fonction réduit cet erreur (difference)

Code :

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

1
2
3
4
5
6
7
8
t=(0:10);
 
y=t.^2;
plot(t,y);
 
Energie=trapz (t,y.^2) 
 
surface=trapz(t,y)

------------resultat matlab
Energie =

20333


surface =

335
------------calcul théorique------
Surface=333.3
Energie=20000

[Invité]

Tu ne réponds pas à ma question :

Mais tu veux je présume calculer l'aire de ta courbe non ?

qui plus est ton code n'a presque pas changé
Là tu calcules au mieux l'intégrale

[adnino]

Oui ,je voudrais calculer l'aire et l'énergie du signal

[Invité]

Donc prend TON signal, et non t.^2.

[adnino]

Vous voulez dire ca ?
Energie=trapz (t,y.^2)

[Invité]

Je veux dire :

Code :

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

Energie=trapz (t,la_variable_qui_contient_les_donnees_de_ton_signal.^2)

et non

Code :

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

1
2
y = n_importe_quoi_d_autre;
Energie=trapz (t,y.^2)

[adnino]

Exactement, Bravo et merci infiniment