Discussion :

[eddyleh]

Bonjour,

Je tente actuellement de modéliser le niveau d'eau dans un bassin soumis à un débit d'entrée, aux précipitations et à l'évaporation.

Mon bassin possède trois composantes: un chenal profond de 3m, une prairie avec une hauteur max de submersion de 0.4 m et des digue sur le pourtour de 1.5 m de haut.

J'ai exprimé l'ensemble de la topographie à l'aide de forme géométriques simples et le tout en fonction de la hauteur. Dès lors lorsque je run mon modèle, le bassin se remplis: d'abord le chenal puis, l'eau se déverse dans la prairie et si tout est plein, l'eau se monte le long des digues.

Pour cela j'utilise une boucle "for" qui parcours mes données et des boucles if:

Si V<V_chenal : calcul de la hauteur, du volume et de la surface dans le chenal

Si V>V_chenal et V<V_chenal+V_prairie: calcul de la hauteur, du volume et de la surface dans la prairie

Si V>V_chenal+V_prairie: calcul de la hauteur, du volume et de la surface dans la prairie.

Le problème est que dans certains cas, entre deux itérations, le volume diminue. Par exemple le niveau d'eau de la prairie s'abaisse et retourne dans le chenal.

A ce niveau, le code matlab n'arrive pas à retourner dans le bloc chenal précédent et me donne un message d'erreur.

Des solutions pour renvoyer mon code au début???

Merci d'avance

[shaiHulud]

Bonjour,

Impossible en matlab de renvoyer un code au début (pas de goto).

Il faut transformer tes GOTO en boucles. Par exemple

ce que tu voudrais (LABEL et GOTO n'existent pas en matlab)

Code :

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

1
2
3
4
5
 
<LABEL 1>
x= traitement1()
if x; <GOTO 1>; end
traitement2()

est équivalent à

Code :

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

1
2
3
4
5
6
 
x= true;
while x
   x= traitement1();
end
traitement2();

[tachmou]

Salut,

Mon commentaire peut soit t'être utile, soit t'embrouiller mais bon (si tu ne vois pas de quoi je parle, ignore ce message) : si tu cherches à faire un code chiadé, je te conseille d'utiliser l'orienté objet... mais si tu préfères développer un simple code de calcul, utilise la méthode de VV33D.

Cordialement,

[le fab]

je vais répondre à coté, mais vu ton problème je ferais plutôt appel à simulink pour le modéliser