Discussion :

[TopCao]

Salut (c'est encore moi...)

Question de débutant : je dispose d'une série de plots(x,y,b) d'une fonction de type exp(-b).

Je voudrai faire l'air sous la courbe en chaque (x,y)
Du coup j'ai écrit la fonction suivante :

Code :

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

1
2
3
4
5
6
7
8
9
10
for x= 1:128
        for y= 1:128
            YVAL = double(diffIm(x,y,:));
            sumVAL =0;
            for i = size(b)
                sumVAL = sumVAL + YVAL(i)*b(i+1)-b(i);
            end
            fitIm(x,y,8) = sumVAL;
        end
    end

Mais évidement je ne peux pas faire b(i+1) si for dépend de size(b)

Bref, je suis coincé...

[Caro-Line]

Hum...je ne sais pas au juste ce que tu veux faire avec ton code mais je pense qu'il y a une incompréhension :

Code :

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

for i = size(b)

signifie que i prendra les valeurs retournées par size(b) (qui sera un vecteur [N M] avec N nombre de lignes de b, et M nombre de colonnes).

Je ne pense pas que ce soit ce que tu souhaites faire, si ?


Sinon ce n'est pas la fonction TRAPZ que tu essayerais de réécrire par hasard ?

[TopCao]

Je pense que size(b) est bien ce que je veux.
b prend 12 valeurs dans mon cas actuel; rangées dans un tableau 12x1
je veux que i prenne les valeurs de 1 à 12 (définies par le nombre de valeurs de b.

Je vais voir si trapz peut m'aider

[duf42]

En fait ce que tu veux c'est plutôt:

Code :

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

size(b,1)

ou encore

Code :

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

length(b)

[Caro-Line]

Oui enfin ce que je relevais c'est surtout que si b est de taille 12*1

Code :

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

for i = size(b)

est équivalent à :

Code :

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

for i = [12 1]

Donc i prendra la valeur 12 puis la valeur 1

Alors que ce que tu veux c'est :

Code :

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

for i = 1:length(b)

Yaka tester pour vérifier

[TopCao]

Merci Caroline, la fonction trapz correspondait exactement à ce que je voulais faire !