question matlab calcul

Bonjour,

je suis débutant en Matlab, j'ai un parametre L qui varie de 1 à 5 avec un pas de temps de 0.5 et je veux calculer "source" à chaque valeur de L tel que:

source = L*exp( -L/25), le programme que j'ai fait et il marche pas et le suivant:

Code:

---

1 2	L=1:0.5:5; source = L*exp( -L/25)

---

mais me signal erreur que le matrix et pas valide!!t

tu peux faire une boucle sur les éléments de ta colonne L

Code:

---

1 2 3 4 5 6 7 8

T=[0:1:50].'; %définition d'un vecteur colonne "temps" L=[1:0.5:5].'; % tu obtiens une matrice colonne de 9 éléments source=zeros(1,9).'; % tu définis ta matrice source et l'initialise for i=1:1:9 source=L(i)*exp(-L(i)/25); % calcul de la valeur de "source" pour les différents L end; plot(T,source); % tracé de "source" en fonction du temps

---

Voilà, il y a surement plus simple, si c'est le cas des membres avisés de ce forum de te le diront :ccool:

Citation:

---

Envoyé par Nathaniel_etudiant

tu peux faire une boucle sur les éléments de ta colonne L

Code:

---

1 2 3 4 5 6 7 8

T=[0:1:50].'; %définition d'un vecteur colonne "temps" L=[1:0.5:5].'; % tu obtiens une matrice colonne de 9 éléments source=zeros(1,9).'; % tu définis ta matrice source et l'initialise for i=1:1:9 source=L(i)*exp(-L(i)/25); % calcul de la valeur de "source" pour les différents L end; plot(T,source); % tracé de "source" en fonction du temps

---

Voilà, il y a surement plus simple, si c'est le cas des membres avisés de ce forum de te le diront :ccool:

---

merciii bcp c'est ce que je cherche même si ya plus simple c'est trop tard ùnt ta déja donné la réponse :D

Salut,

il ne manquait pas grand chose pourtant :

Code:

---

1 2	L=1:0.5:5; source = L.*exp( -L/25)

---

L'opérateur .* permet d'effectuer une multiplication termes à termes et non une multiplication matricielle.

Hum, je pense m'être trompé pour la définition de "source". Il faudrait écrire "source(i)" à la place de "source". De plus il faut insérer la ligne "plot(T,source(i))" dans la boucle for.

Citation:

---

Envoyé par Nathaniel_etudiant

Hum, je pense m'être trompé pour la définition de "source". Il faudrait écrire "source(i)" à la place de "source". De plus il faut insérer la ligne "plot(T,source(i))" dans la boucle for.

---

voilà ce que j'ai fait dans monpetit programme:

Code:

---

1 2 3 4 5 6

T = [300:50:3800].'; source=zeros(1,71).'; for i=[1:1:71].'; R(i) =6.38e14*(1/T(i))*exp(-0.384e5/T(i)); end; plot(T,R(i))

---

il me donne erreur :

Code:

---

inner matrix dimensions must agree????

---

tu peux aller voir la FAQ ici
http://matlab.developpez.com/faq/?pa...ons_must_agree

Dans la boucle il faut écrire

Code:

---

1 2	for i=1:1:71;

---

et non

Code:

---

1 2	for i=[1:1:71].';

---

de plus met la commande "plot" dans la boucle, et pas en dehors. tu peux rajouter la commande "hold on;" dans la boucle, à la suite de la commande "plot".

Je ne sais pas si c'est ce que tu veux faire..

Pourquoi mettez-vous des .' à chaque fois que vous créez un vecteur?

Code:

---

1 2 3 4 5 6

T = [300:50:3800]; R=zeros(1,71); for i=1:71; R(i) =6.38e14*(1/T(i))*exp(-0.384e5/T(i)); end; plot(T,R)

---

Ou plus simplement (sans boucle) :

Code:

---

1 2 3

T = [300:50:3800]; R =6.38e14.*(1./T).*exp(-0.384e5./T); plot(T,R)

---

Citation:

---

Envoyé par Nathaniel_etudiant

Dans la boucle il faut écrire

Code:

---

1 2	for i=1:1:71;

---

et non

Code:

---

1 2	for i=[1:1:71].';

---

de plus met la commande "plot" dans la boucle, et pas en dehors. tu peux rajouter la commande "hold on;" dans la boucle, à la suite de la commande "plot".

Je ne sais pas si c'est ce que tu veux faire..

---

oui c'est ce que je veux faire mais me signal tjrs erreur: inner matrix dimensiosn must agree !! et si tu as remarqué que a matrix elle a les même dimensiosn si je sais bien compter!! non?

Magelan: oui c'est inutile, ":" est l'opérateur colonne, pas besoin de faire la transposée avec ".' ". Autant pour moi :roll:

Avec ce code là, je n'ai pas d'erreur.

Code:

---

1 2 3 4 5 6 7 8

T = 300:50:3800; R=zeros(1,71); for i=1:1:71; R(i) = 6.38e14*(1/T(i))*exp(-0.384e5/T(i)); plot(T,R(i)); hold on; end;

---

Citation:

---

Envoyé par Nathaniel_etudiant

Magelan: oui c'est inutile, ":" est l'opérateur colonne, pas besoin de faire la transposée avec ".' ". Autant pour moi :roll:

Avec ce code là, je n'ai pas d'erreur.

Code:

---

1 2 3 4 5 6 7 8

T = 300:50:3800; R=zeros(1,71); for i=1:1:71; R(i) = 6.38e14*(1/T(i))*exp(-0.384e5/T(i)); plot(T,R(i)); hold on; end;

---

---

OK merci tout le monde