Discussion :

[george216]

Chers amis, je vous souhaite tout d'abord une année 2010 pleine de bonheur et de succès.
Ma discussion concerne qques manipulations d'un nombre très élevé de colonnes d'une matrice. En fait, j'ai une matrice A de format(100,20000). Chaque colonne de A est composée de -1, 0 et 1. les différents passages de -1 à 1 définissent des expansions, pour chaque colonne. La durée de chaque expansion est la soustraction de l'indice de 1 celui de (-1) (Bien évidemment, les éléments 1 et -1 qui interviennent pour cette expansion). De même, les différents passages de 1 à -1 constituent des récessions, pour chaque colonne.
La durée de chaque récession est la soustraction de l'indice de (-1) celui de 1. Je voudrais dire par indice le numéro d'observations dans la colonne.
Pour mieux éclaircir, je donne l'exemple suivant:
soit le vecteur u comportant 10 observations:
[1 0 0 0 0 -1 0 -1 0 1 ]'
Dans ce vecteur u il y a une récession:le passage du premier élément du vecteur (qui est 1) au 6ème élément (qui est -1) et la durée de cette récession est 6-1=5.
Dans ce vecteur u il y a une expansion: le passage du 8ème élément qui est -1 au 10ème élément qui est 1. La durée de cette expansion est 10-8=2.
Y a-t-il quelqu'un qui peut m'aider par un programme qui peut déterminer pour chaque colonne de ma matrice A les différentes expansions, récessions ainsi que leurs durées. SVP, j'ai besoin de ça dans ma thèse. Merci d'avance.
George

[Jerome Briot]

As-tu au moins essayé de coder une simple boucle FOR-END ?

[phryte]

Bonjour.
Une idée à creuser :

Code :

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

1
2
col=[1 0 0 0 0 -1 0 -1 0 1]';
r=find(col==-1)-find(col==1)

[george216]

J'ai bien réflichi à la boucle for end et à la fonction "find". Le problème est que les colonnes n'ont pas le même nombre d'expansions et de récessions et par là, elles n'ont pas le même nombre de durées. D'ailleurs la fonction find donne les différents nombres de durée de toutes les colonnes de la matrice A dans un vecteur colonne.

[phryte]

Tu peux tester une seule colonne :

Code :

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

1
2
A=[1 0 0 0 0 -1 0 -1 0 1;0 0 1 0 0 0 -1 0 0 1]';
r=find(A(:,1)==-1)-find(A(:,1)==1)

[george216]

Est-il possible de faire cette opération sur toutes les colonnes de A? je voudrais avoir une vue d'ensemble sur les durées pour toutes ces colonnes.

[phryte]

Bonjour.
Je n'ai pas réussi à me passer de la boucle for !

Code :

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

1
2
3
4
5
A=[1 0 0 0 0 -1 0 -1 0 1;0 0 1 0 0 0 -1 0 0 1]';
for n=1:2
r(:,n)=find(A(:,n)==-1)-find(A(:,n)==1);
end
r