Discussion :

[madchemiker]

Bonsoir,

J'ai établi un modèle simulink dépendant de beaucoup de paramètres. Afin de fitter certains paramètres inconnus, je compte me caler sur des manips. J'ai donc crée une fonction dans un fichier.m ou je calcule un écart à une valeur souhaitée.

Allure de ma fonction

Code :

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

1
2
3
4
5
6
7
8
9
10
11
function ecart=Optimisation(X)
 
assignin ('base','deltaGmoinsa',X(1));
assignin ('base','deltaGmoinsd',X(2));
assignin ('base','rDEG',X(3));
 
sim('fitting',3600000);
 
R0=resistance.signals.values(100);
Rf=resistance.signals.values(length(resistance.signals.values));
ecart=(Rf-R0-0.8e-3)^2;

Ensuite j'appelle cette fonction dans la toolbox Optimisation et je recherche le minimum de ma fonction. Le problème c'est que quand Simulink prend un jeu de valeurs faisant 'planter' la simulation, je perds tout.

On m'a vaguement dit "Fais un try ...catch".

Merci mais je ne sais pas quoi en faire.

Qqun aurait'il une piste pour passer outre le plantage et redémarrer avec de nouvelles données?

Merci beaucoup

[Jerome Briot]

Je ne suis pas sûr de répondre exactement à ta question mais voici un exemple très simple.

Le code suivant :

Code :

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

1
2
i=-2;
X(i) = rand;

ne fonctionne pas est renvoie le message :

Code :

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

??? Attempted to access (-2); index must be a positive integer or logical.

=> voir la : Subscript indices must either be real positive integers or logicals.

On peut donc trouver une alternative comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
i=-2;
 
while 1
    try
        X(i) = rand;
        break
    catch
        i=i+1;
    end
end

Ce code (très théorique) fonctionne.

Pour voir ce qui se passe exactement, il suffit de modifier un peu le CATCH comme ceci :

Code :

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

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
i=-2;
 
while 1
    try
        X(i) = rand;
        break
    catch ME
 
        disp(ME.message)
 
        i=i+1;
    end
end
 
X

Ce qui retourne :

Code :

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

1
2
3
4
5
6
7
8
 
Attempted to access (-2); index must be a positive integer or logical.
Attempted to access (-1); index must be a positive integer or logical.
Attempted to access (0); index must be a positive integer or logical.
 
X =
 
    0.6557

Les erreurs sont donc bien retournées, mais le programme ne s'arrête pas pour autant.