Discussion :

[Dam2227]

Bonjour

Y'aurait-il une commande Matlab pour additionner des vecteurs colonnes et lignes de telle sorte qu'ils forment une matrice à la fin?

Par exemple, si j'ai deux vecteurs

Code :

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

1
2
3
 
a=[1 2 3];
b=[4;5;6];

je souhaiterais que en les additionnant, je finisse avec :

Code :

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

1
2
 
c=[5 6 7 ; 6 7 8 ; 7 8 9]

L'idée derrière ça est que j'ai une fonction donc chaque point est défini par une intégrale. Pour le moment, je calcule d'abord les éléments dans l'intégrale dans une boucle for et ensuite j'intègre le tout avec trapz en un seul coup. Je souhaiterais éviter la boucle for en calculant en un seul coup ma matrice, plutôt que de rentrer dans une boucle. Y'a t-il un moyen? Je peux poster un exemple si ce n'est pas clair.

Merci!

[Invité]

Bonsoir,

Avec la fonction bsxfun :

Code :

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

1
2
3
4
5
6
7
8
9
a=[1 2 3];
b=[4;5;6];
bsxfun(@plus,a,b)
 
ans =
 
     5     6     7
     6     7     8
     7     8     9

[Dam2227]

Excellent! Merci beaucoup pour la réponse rapide!

[Dam2227]

Humm, en fait, ça ne résout pas mon problème. Le truc c'est que ma fonction ne me permet pas de séparer les deux vecteurs. Je ne vais pas la poster ici parce qu'elle est trop longue, mais voici un exemple :

Code :

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

1
2
3
4
5
6
7
8
9
10
 
    t=0.5:0.1:1;
    w=0:1:10;
    lengthw=length(w);
    lengtht=length(t)
    for i=1:lengthw
        c(i,1:lengtht)=(t.^2).*sin(w(i)+t);
    end
    c
    trapz(t,c,2)

C'est une fonction qui dépend de w, et chaque point est intégré selon t.

Est-ce que je peux me passer de la boucle for dans ce cas là?

Merci!

[Invité]

Il suffit d'adapter la fonction dont le handle est passé en paramètre :
Voir Que représente le symbole @ ?

Code :

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

1
2
3
4
t=0.5:0.1:1;
w=0:1:10;
c = bsxfun(@(T,W) (T^2).*sin(W+T),t,w.')
trapz(t,c,2)

[Dam2227]

Okeche, j'ai compris.

Merci!