Discussion :

[mostyx]

Bonjour à tous,

J'aimerais réussir à faire une boucle for avec ce type d'équation

x(a un moment T+dt) = x (au moment T)*8

Pour un x allant de 1 a 100
avec un dt petit
Voici mon début de code mais cela ne me donne rien

Code :

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

1
2
3
4
5
6
7
8
9
 
dt=0.01
%Temperature sans pertes
 
u(1)=0;
for t=1:dt:100
t(i+1)=t(i)*8
u(i+1)=u(i)+t1(i)*dt 
end

Merci de votre aide !

[FR119492]

Salut!
Deux remarques:
Tout d'abord, tu n'as pas initialisé T(1).
Ensuite, pour intégrer une équation différentielle, tu utilises la méthode d'Euler, dont on a dit: "Admirez, mais n'appliquez pas".
Jean-Marc Blanc

[mostyx]

Mais en fait je ne veux pas du tout intégrer

je veux simplement tracer le graphique d'une courbe qui avance en fonction de sa valeur précédente

[Jean Dumoncel]

Bonjour,

ta demande et ton code ne sont pas clairs du tout... pourrais-tu réexpliquer ce que tu souhaites faire?

[mostyx]

Quand je tape cela :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
 
clc
clear all;
dt=1
 
u=0;
x=0;
for w=1:100
u=u+1
x=x+1
end
plot(x,u)

Mon plot ne marche pas que dois-je faire?

[Jean Dumoncel]

Mon plot ne marche pas

ben si... tu affiches un point de coordonnées (100,100), en lui-même le code n'est pas faux, il ne génère pas de message d'erreur (même si j'ai bien compris qu'il ne réalise pas ce que tu souhaites faire).

que dois-je faire?

Nous ré-expliquer ce que tu souhaites faire précisément?

[tubaas]

salut

Mon plot ne marche pas

A mon avis il "marche" seulement tu essaye de tracer uns calaire donc c'est un unique point aux coordonnées (100,100)
tu veux peut-être faire quelque chose comme cela

Code :

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

1
2
3
4
5
6
7
8
9
10
maxlim=100;
u=zeros(1,maxlim);
x=u;
u(1)=0;
x(1)=0;
for w=2:maxlim
u(w)=u(w-1)+1;
x(w)=x(w-1)+1;
end
plot(x,u)

ensuyite la bonne question est celle de Magelan : que cherches-tu à faire exactement ? Si c'est la résolution d'un équation différentielle il existe des schémas suivant le type d'équation.
Tu as des cours sur internet comme celui-ci ou encore celui-là pour les équations aux dérivées partielles.

[mostyx]

En fait , ce que je cherche réellement a faire , c'est mettre cette équation
dans matlab.

et d'en obtenir le graphique

Ou plutot : Tb(t+Δt)= Tb(t) (PcombΔt + 4α + √( Pcomb2Δt 2 + 8α PcombΔt + 4α2))/6α

Le code en est a ca :

dt=0.1;
t=1:dt:2;
for i=1:50
t(i+1)=t(i)+dt ;
T(i)=T(1)*(Pcarbu*dt+4*alpha+(Pcarbu^2*dt^2+8*alpha*Pcarbu*dt+ 4*alpha^2)^0.5)/(6*alpha);
end

alpha , Pcarbu et dt sont des valeurs définie au début.
Mais malheureusement , je recois en rouge des erreurs me disant que les vecteurs ne sont pas de même longueur

Je vérifie cela je tape donc length(T) et length(t) et je vois qu'il sont effectivement différent :s
Comment remedier a cela?


Ps mes question précédente avait pour but de savoir si je comprenais bien comment fonctionnait les boucles )

[Jean Dumoncel]

Ps mes question précédente avait pour but de savoir si je comprenais bien comment fonctionnait les boucles )

Tes boucles for sont mal construites mais comme on ne savait pas ce que tu voulais faire, c'est difficile de te répondre...

Bon j'essaie de résumer :
tu as donc un vecteur T pour différentes valeurs de t (par exemple pour t=1, T(1)=... pour t=2, T(2)=..., etc). Et tu voudrais calculer T(t+dt) avec par exemple dt=0.1 en appliquant ta formule, tu obtiendrais alors T(1.1), T(2.1), etc...
C'est bien ça?

[mostyx]

Merci de votre aide mais j'ai réussi !
Ce fut fastidieux mais maintenant c'est fait!

Sinon tu ne saurais pas par hasard comment lancer un fichier matlab (.m) depuis matlab une commande ?
Genre launch(superfichierdelamort.m) ?

[Jean Dumoncel]

Merci de votre aide mais j'ai réussi !
Ce fut fastidieux mais maintenant c'est fait!

Par curiosité (et pour donner un peu de sens à ce fil) tu aurais pu nous montrer ton code...

Sinon tu ne saurais pas par hasard comment lancer un fichier matlab (.m) depuis matlab une commande ?
Genre launch(superfichierdelamort.m) ?

Dans Matlab? il suffit de taper le nom de ton fichier.

Code :

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

>> superfichierdelamort

[mostyx]

Voici la solution

Code :

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

1
2
3
4
5
6
7
8
9
10
11
%Temperature sans pertes
T(1)=293;
t(1)=0 
dt=0.1;
 
Tb(1)=293
for i=1:30000  %Pour 2000 secondes
t(i+1)=t(i)+dt ;
T(i+1)=T(i)*(Pcarbu*dt+4*alpha+(Pcarbu^2*dt^2+8*alpha*Pcarbu*dt+ 4*alpha^2)^0.5)/(6*alpha);
end
plot(t,T);

Et merci pour l'execution d'un sous fichier