Discussion :

[clem3900]

Bonjour,

je voudrais calculer une intégrale en utilisant la fonction quad, mais je ne vois pas comment exprimer ma fonction T, en fonction de ma variable lambda...
Voici mon programme:

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
function T = integral(????)
 
fichier='XXX.asc';
S=load(fichier);
 
 
i=2:152;
 
for k=1:1024
 
    Q= S(k,1);
    R = sum(S(k,i));
    T(k,:) = [Q R];
end
 
plot(T);

En fait j'ai un tableau de 1024x152, dont la première colonne contient l'ensemble des lambda (ils correspondent aux S(k,1)) servant pour l'axe des abscisses. Et pour chaque ligne je calcule la somme, et je trace le résultat de cette somme en fonction des lambda.
La courbe que j'obtiens est bonne (car j'ai la même sous excel), mais sur l'axe des abscisses j'ai les valeurs de k au lieu de lambda ( je ne pense pas que se soit normal).
Et du coup je ne vois pas comment exprimer ma fonction T en fonction de ma variable lambda, car celle-ci n'apparait pas dans le code vu qu'il s'agit de valeurs de mon tableau...

Si quelqu'un a une idée.

Merci

[TRSRT]

En premier, il me semble que ta boucle for est inutile.

Tu peux écrire:

Code :

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

1
2
    Q= S(:,1);
    R = sum(S(1:1024,2:152),2);

Pour le reste, je pense que le problème est dans ta manière de définir le plot.
Ecris:

Code :

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

Plot(Q,R)

et cela doit marcher

Pour le calcul de l'intégral, je n'ai jamais utilisé quad mais trapz(X,Y) qui doit fonctionner ici il me semble

[clem3900]

j'ai changé mon programme comme tu l'as dit, ce qui fonctionne très bien (merci). Par contre pour la fonction trapz matlab ne m'affiche aucune valeur...
Mon programme est le suivant:

Code :

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

1
2
3
4
5
6
7
8
fichier='XXX.asc';
S=load(fichier);
 
   Q= S(:,1);
   R = sum(S(1:1024,2:152),2);
   plot(Q,R);
 
Z = trapz(R,Q);

[TRSRT]

Corrige moi si je me trompe mais c'est bien en fonction de lambda que tu veux intégrer non?
(En gros calculer l'aire sous la courbe dessinée juste avant?)

Donc ce serait plutôt trapz(Q,R) non?

[clem3900]

Comme j'étais pas sûr j'ai essayé les deux ^^
mais dans les deux cas ça ne m'affiche rien...

[TRSRT]

ton Z n'a aucune valeur?

Moi j'ai fait cela pour tester:

Code :

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

1
2
3
4
5
S=[1 4 7 10;2 5 8 11;3 6 9 12]
Q= S(:,1)
 R=sum(S(1:3,2:4),2)
plot(Q,R)
trapz(Q,R)

Et il me dit bien que l'intégrale vaut 48, ce qui correspond à la valeur calculée à la main

[clem3900]

...
Non je n'ai rien moi, dans ma commane window j'ai:

Code :

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

>>

et c'est tout...
Peut être que le problème vient de la taille de S.
S fait 1024 par 152, peut être que c'est trop grand pour utiliser trapz, non?

[clem3900]

je suis stupide...
J'avais mis

Code :

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

;

après trapz...
Merci pour votre aide!!

[TRSRT]

Bon bah parfait, problème résolu, content d'avoir pu t'aider