Discussion :

[zizo08]

Bonjour,

J'ai besoin de votre aide svp.
Mon problème est le suivant: je dois coder une fonction qui me permet de calculer la complexité définit par l'expression suivante:
D(P(i)|P(i+1))=sum(P(i)*log(p(i)/p(i+1));
ses probabilités sont calculées à partir d'un signal audio x.
Pour résoudre ce problème, j'écris ce code:

Code :

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

1
2
3
4
5
6
7
w=166;
for i=1:100
  y=x((i-1)*w+1:i*w);  
hist_norm(i)=hist(y)./sum(hist(y)); 
comp(i)=sum(hist_norm(i).*log(hist_norm(i)./hist_norm(i+1)));
complexity=[complexity comp(i)];
end

mais ça n'a pas marché.

Y a t-il quelqu'un qui peut m'aider?

Merci beaucoup en avance.

[Caro-Line]

mais ça n'a pas marché.

Qu'est-ce qui n'a pas "marché" ?

Message d'erreur ? Résultat incorrect ? ...

[tug83]

Je ne sais pas si c'est que tu cherchez mais tu peux déterminer la complexité d'un programme dans MATLAB avec le mlint:

Code :

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

mlint('filename','-cyc')

[zizo08]

en fait, le message d'erreur est indiqué au niveau de cette ligne:

Code :

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

hist_norm(i)=hist(y)./sum(hist(y));

et indique que 'A[I]=B should be the same dimension'.
comment pourrais je résoudre ce problème.

le x que je traite est un vecteur construit à partir d'un signal audio.

merci pour votre aide.

[zizo08]

j'ai voulu ajouter une autre clarification.
le programme que je code dois implémenter cette expresion.
D(P(i)|P(i+1))=sum(P(i)*log(p(i)/p(i+1));
mais je ne suis pas certaine de la façon dont je clacule les probabilités.

merci une autre fois

[Caro-Line]

Ce ne serait pas plutôt ça le message d'erreur :
In an assignment A(I) = B, the number of elements in B and I must be the same.

Merci de faire un copier-coller exact des messages d'erreur, c'est plus facile pour trouver le problème.

A priori hist(y)./sum(hist(y)) est un vecteur et tu essayes de le mettre dans un scalaire, donc forcément MATLAB n'aime pas.

Pour le reste je ne peux pas vraiment t'aider, mes connaissances en probas remontant à loin loin loin.

[zizo08]

est ce qu'il y a pas une façon pour enregistrer cette quantité
hist(y)./sum(hist(y))
parce que pour moi je dois calculer le différence de deux quantités successives=> je dois changeer la valeur de y pour chaque i.


je suis désolée pour les erreurs que je commise.

merci pour votre aide

[Caro-Line]

Tu peux essayer avec un tableau de cellules.

[Jerome Briot]

De toute façon, peut importe ce que tu calcules, ton code ne fonctionnera pas :

Code :

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

1
2
3
4
5
6
7
for i=1:100
 
hist_norm(i)=hist(y)./sum(hist(y)); 
comp(i)=sum(hist_norm(i).*log(hist_norm(i)./hist_norm(i+1)));

end

hist_norm(i+1) ne sera définit qu'à l'itération suivante.

Exemple plus simple :

Code :

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

1
2
3
4
for n=1:5
   a(n) = n;
   b(n) = a(n)+a(n+1);    
end

Code :

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

1
2
3
4
??? Attempted to access a(2); index out of bounds because numel(a)=1.
 
Error in ==> Untitled at 3
   b(n) = a(n)+a(n+1);

[zizo08]

merci

la boucle commence à partir de i=1.
J'ai définit hist_norm(0)=1 en dehors de la boucle.

[Jerome Briot]

J'ai définit hist_norm(0)=1 en dehors de la boucle.

Comme tu veux, mais ça ne marchera pas plus parce que :

1. l'indexage commence à 1 avec MATLAB.
2. ce n'est pas i-1 que tu utilises dans la boucle, mais i+1

[zizo08]

oui vous avez tout à fait raison.
comment pourrais je faire donc?

merci pour l'aide

[Jerome Briot]

Il faut calculer hist_norm et comp séparément.

Je te suggère par exemple d'utiliser 2 boucles FOR-END distinctes

Code :

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

1
2
3
4
5
6
7
8
 
for n = 1:5
   a(n) = n;
end
 
for n=1:4
   b(n) = a(n)+a(n+1);
end

[zizo08]

merci beaucoup pour votre aide.
c'est bon maintenant.

merci