[plot] La fonction ne relie pas mes points

j'ai une question à savoir que je veux faire une moyenne sur les colonnes de la matrice Pt0 (une 600x9) d'un certain type à savoir que je veux moyenner des valeurs espacées de 14 intervalles (la valeur à considérer est située chaque 15 lignes de la précédente) mais je ne vois pas comment je peux faire avec la fonction sum car je peux spécifier si je suis sur une ligne ou une colonne mais je ne vois pas comment sauter lors de la sommation. J'obtiendrai donc une matrice 40*9 .

En combinant ce que je t'ai montré dans mon message précédent, et mon tout premier:

Citation:

---

Envoyé par Winjerome

Tu peux faire les deux cumulés: sum( sum(Pt0(k:k+14,:) ).
Une solution sans boucle serait d'utiliser RESHAPE sur Pt0 avant d'appliquer SUM dessus.

---

Citation:

---

Envoyé par Winjerome

Code:

---

1 2 3

for k=1:15:size(Pt0,1) ... end

---

Ou version sans boucle:

Code:

---

reshape(Pt0,15,[])

---

---

Sachant que maintenant, tu as une dimension en plus.

avec la fonction reshape y a t'il moyen de spécifier qu'elle me mette que les valeurs par pas de 15 sur la première colonne et ainsi de suite

ex: première colonne je veux avoir les valeurs Pt0(1,1) Pt0(15,1) Pt0(30,1) ...
deuxième colonne les valeurs Pt0(2,1) Pt0(16,1) Pt0(31,1)


j'ai lu la notice dans le help que cette fonction organise en lisant colonne par colonne


merci beaucoup pour l'aide winjerome

Que fais-tu de la deuxième dimension?

Petite rectification:

Code:

---

1:15:600 = 1 16 31 46 61 ...

---

Donc

Citation:

---

Envoyé par hittie

ex: première colonne je veux avoir les valeurs Pt0(1,1) Pt0(16,1) Pt0(31,1) ...
deuxième colonne les valeurs Pt0(2,1) Pt0(17,1) Pt0(32,1)

---

Je t'ai déjà donné tous les éléments.
Avec Pt0 de taille 600x1, on avait reshape(Pt0,15,[]) ou reshape(Pt0,15,40) vu que 600 = 15*40.
Reste à ajouter une dimension, et appliquer mon petit cours sur SUM ;)

j'ai rajouté cette ligne de code

Code:

---

G=reshape(Pt0,[15 40 9]);

---

et je voulais savoir si un tableau sur 3 dimensions se traite de la même façon dans le cas d'une sommation qu'un tableau 2 D

Exactement la même façon:
sum(maVariable,3) somme sur la troisième dimension
sum(maVariable,4) somme sur la quatrième dimension
etc...

winjerome j'ai avancé et je voulais savoir si

Code:

---

1 2 3

for i=1:9 H(:,:,i)=sum(G(:,:,i),2)/40; end

---

correspond bien à la moyenne sur les lignes de la matrice G (une 15*40*9) (j'ai l'impression que oui car la dimension correspond )

Oui, tu fais bien la moyenne sur chaque ligne (tu peux aussi utiliser MEAN qui fonctionne de la même façon en ce qui concerne le maniement des dimensions).
Par contre:

Code:

---

1 2 3

for i=1:9 H(:,i)=sum(G(:,:,i),2)/40; %ou mean(G(:,:,i),2); end

---

Serait plus juste, vu que tu obtiens un tableau 2D.
Ou sans boucle: H = mean(G,2);, voire H(:,:) = mean(G,2); pour obtenir le même tableau 2D

ok merci winjerome (en fait je voulais garder les 9 tableaux séparés).

Sinon j'ai une question : j'essaie de faire une multiplication entre deux matrices :

Code:

---

1 2 3

for i=1:9 H(:,:,i)=AreaMatrix(1,:,i)*G(:,:,i); end

---

AreaMatrix est une 1x15x9 et G une 15x40x9 donc normalement la multiplication devrait s'effectuer sauf que j'obtient comme erreur mtimes (probleme avec la multiplication)...

AreaMatrix étant comme tu dis un tableau 1x15x9 et G un 15x40x9, on aurait: AreaMatrix(1,:,i) de taille 1x15 et G(:,:,i) de taille 15x40, donc normalement pas de problème...
As-tu vérifié les dimensions des deux?

en fait j'avais oublié que j'avais une autre boucle qui n'avait pas été mise en commentaire.

Sinon j'essaie de faire la moyenne en divisant la somme des valeurs sur une ligne en divisant par l'aire totale de la ligne (l'aire infinitésimale contenue dans chaque colonne x le nombre d'élément par ligne)

Code:

---

1 2 3 4 5

for i=1:9 for j=1:15 H2(:,:,i)=sum(AreaMatrix(1,:,i)*G(:,:,i),2)/(AreaMatrix(1,j,i)*40); end end

---

mais je n'arrive pas car j'obtient une 1x1x9 au lieu d'une 15x1x9
mais

Même genre de remarque que dans mon message #28: chaque ligne

Code:

---

sum(AreaMatrix(1,:,i)*G(:,:,i),2)/(AreaMatrix(1,j,i)*40);

---

est un simple scalaire, il n'y a donc pas lieu de faire H2(:,:,i).
Réfléchi: tu boucles sur i et j ;)

oui je voulais reprendre la même expression que pour la matrice H1 (moyenne arithmétique) sauf que je n'arrive pas à voir comment je peut diviser la somme par la valeur de AreaMatrix correspondante à la ligne mutlipliée par le nombre d'éléments de la ligne (voilà pourquoi j'avais mis l'indice j)