Discussion :

[stage75]

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 :

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

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

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

[Nathaniel_etudiant]

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

Code :

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

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

[stage75]

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

Code :

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

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

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

[Jean Dumoncel]

Salut,

il ne manquait pas grand chose pourtant :

Code :

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

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.

[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.

[stage75]

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 :

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

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 :

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

inner matrix dimensions must agree????

[Nathaniel_etudiant]

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

[Nathaniel_etudiant]

Dans la boucle il faut écrire

Code :

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

1
2
 
for i=1:1:71;

et non

Code :

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

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..

[Jean Dumoncel]

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

Code :

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

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 :

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

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

[stage75]

Dans la boucle il faut écrire

Code :

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

1
2
 
for i=1:1:71;

et non

Code :

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

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?

[Nathaniel_etudiant]

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

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

Code :

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

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;

[stage75]

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

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

Code :

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

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