Ajouter des résultats de simulation dans une variable

**tomazy4T6** · 07/01/2016, 15h04

Bonjour à toutes et à tous,

Je travaille sous Matlab et SimHYDRAULICS.
Le script suivant permet de lancer des simulations de manière bouclée.
Lorsque une simulation est terminée, je demande à Matlab de stocker la variable F.signals.values (à n lignes et 1 colonne) dans ma variable effort, puis d'enregistrer cette variable.
Chaque colonne de la variable effort doit contenir les valeurs de F.signals.values de chaque itération.
De plus, à chaque itération, la variable F.signals.values change de taille.

La première itération fonctionne puis à la deuxième, Matlab, affiche : "Subscripted assignment dimension mismatch".
Je ne crée pas la variable effort car je ne connais pas sa dimension avant de lancer la simulation.

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
17
 
clear jj
 
rel_tol=5e-5;
max_step=5e-5;
 
for jj=1:1:taille-1
 
opt=simset('Solver','ode15s','MaxStep',max_step,'RelTol',rel_tol);
NomModele='mdlSA19SImhydroV1try4_carto';
sim(NomModele,[Tdeb Tfin(jj,1)],opt)
effort(:,jj)=F.signals.values(:,1);
vaR=effort(:,jj);
filename = 'fichier_effort';
save(filename, 'vaR');
 
end

Ce que j'aimerai obtenir, c'est une variable effort sous forme de tableau avec dans chaque colonne les valeurs de la variable F.signals.values à chaque itération.
Pourriez vous s'il vous plaît m'aider?

**Gooby** · 07/01/2016, 15h09

Bonjour

Le problème est que tu veux concaténer (par exemple)

ceci:
1
2
3

avec ceci
5
5
5
5

Problème de dimension donc.

La question est: comment souhaites tu remplir les "trous" de ta matrice?

Deux solutions, réadapter la taille de ta matrice dynamiquement ou utiliser un tableau de cellules.

**tomazy4T6** · 07/01/2016, 15h48

Bonjour Gooby,

Voilà, je ne t'ai pas précisé que je voulais remplir les trous manquants avec des zéros.

Exemple:
la variable F.signals.values contient pour chaque itération les valeurs suivantes:
ité 1:[1 2 3 4 5 6 7]
ité 2:[1 3 5 7 9]
ité 3:[2 4 6]
... etc

Je voudrais pouvoir sauvegarder dans un fichier .mat une seule variable se nommant "var" et contenant le tableau suivant:

ité n° 1 2 3 ...etc

1 1 2
2 3 4
3 5 6
4 7 0
5 9 0
6 0 0
7 0 0

**Gooby** · 07/01/2016, 15h59

Dans ce cas, il va falloir un peu adapté ton code pour faire 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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
 
clear jj
 
rel_tol=5e-5;
max_step=5e-5;
 
for jj=1:1:taille-1
 
opt=simset('Solver','ode15s','MaxStep',max_step,'RelTol',rel_tol);
NomModele='mdlSA19SImhydroV1try4_carto';
sim(NomModele,[Tdeb Tfin(jj,1)],opt)
 
% Il faut ici adapté les tailles de matrices pour que la concaténation puisse avoir lieu
if ( size(F.signals.values(:,1),1) > size(effort,1) )
    myDiff = size(F.signals.values(:,1)) - size(effort,1);
    effort(end+1:end+myDiff,:) = zeros(myDiff,size(effort,2)); % on complète effort avec des lignes de zéros
    effort(:,jj) = F.signals.values(:,1);
elseif ( size(F.signals.values(:,1),1) < size(effort,1) )
     myDiff = size(effort,1) - size(F.signals.values(:,1));
     temp = cat(1,F.signals.values(:,1),zeros(myDiff,1)); % on crée un vecteur temporaire de même longueur que effort (donc F.signals.values complétés par des zéros)
     effort(:,jj) = temp;
else
    effort(:,jj)=F.signals.values(:,1);
end
 
vaR=effort(:,jj);
filename = 'fichier_effort';
save(filename, 'vaR');
 
end

**tomazy4T6** · 07/01/2016, 16h26

Je viens de lancer le code en utilisant ton script.
Bon là où ça coince, c'est qu'il ne connaît pas la variable "effort" (Undefined function or
variable 'effort').

Si je déclare la variable "effort" avant la boucle conditionnelle, le code tourne sans erreur mais les résultats ne sont pas bons.
Je vous transmet un bout du résultat:

12776,5600000000 0 0 0 12776,5600000000
12776,5600000000 0 0 0 12776,5600000000
12776,5600000000 0 0 0 12776,5600000000
12776,5600000002 0 0 0 12776,5600000002
12776,5600000002 0 0 0 12776,5600000002
12776,8695936674 0 0 0 12776,8695936674

Pensez-vous qu'il soit plus facile d'utiliser un tableau de cellule dans mon cas?

**Gooby** · 07/01/2016, 17h00

Moi, de base, j'aurai utiliser un tableau de cellules. Ensuite, dans ton cas, on peut s'en sortir en facilement en faisant 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
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
 
clear jj
 
rel_tol=5e-5;
max_step=5e-5;
 
for jj=1:1:taille-1
 
opt=simset('Solver','ode15s','MaxStep',max_step,'RelTol',rel_tol);
NomModele='mdlSA19SImhydroV1try4_carto';
sim(NomModele,[Tdeb Tfin(jj,1)],opt)
 
% Il faut ici adapté les tailles de matrices pour que la concaténation puisse avoir lieu
if ( jj == 1 )
    effort(:,jj)=F.signals.values(:,1); % Comme ça ton tableau s'initialise à jj = 1
elseif ( size(F.signals.values(:,1),1) > size(effort,1) )
    myDiff = size(F.signals.values(:,1)) - size(effort,1);
    effort(end+1:end+myDiff,:) = zeros(myDiff,size(effort,2)); % on complète effort avec des lignes de zéros
    effort(:,jj) = F.signals.values(:,1);
elseif ( size(F.signals.values(:,1),1) < size(effort,1) )
     myDiff = size(effort,1) - size(F.signals.values(:,1));
     temp = cat(1,F.signals.values(:,1),zeros(myDiff,1)); % on crée un vecteur temporaire de même longueur que effort (donc F.signals.values complétés par des zéros)
     effort(:,jj) = temp;
else
    effort(:,jj)=F.signals.values(:,1);
end
 
vaR=effort(:,jj);
filename = 'fichier_effort';
save(filename, 'vaR');
 
end

Ensuite, pour le fait que le "résultat n'est pas bon", c'est difficile à comprendre avec le morceau que tu nous donnes. De façon simple, qu'obtiens-tu?

**tomazy4T6** · 07/01/2016, 17h18

Vous avez raison, le tableau de cellule est plus pratique. J'essaierai de mon côté de faire fonctionner le code avec cette méthode.

Je veux cependant y arriver avec l'autre méthode; celle pour laquelle vous m'avez aidé.
Quand je lance le code tel quel, il m'affiche: Undefined function or variable 'effort.

Il faut donc déclarer la variable effort en amont. Celle-ci est égale au vecteur présent dans la variable issue de la simulation: F.signals.values.
J'ai rajouté la ligne 3 après la ligne 2 de lancement de la simu.

Code :

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

1
2
3
 
sim(NomModele,[Tdeb Tfin(jj,1)],opt)
effort=F.signals.values;

J'ai ajouté cela hasardeusement pour que le code tourne mais cela ne fonctionne pas.
A savoir que la variable "effort" donne le tableau suivant (j'en transmet que le début):
12776,5600000000 0 0 0 12776,5600000000
12776,5600000000 0 0 0 12776,5600000000
12776,5600000000 0 0 0 12776,5600000000
12776,5600000002 0 0 0 12776,5600000002
12776,5600000002 0 0 0 12776,5600000002
12776,8695936674 0 0 0 12776,8695936674

Les colonnes 1 et 5 sont identiques (elles devraient être différentes car les simulations 1 et 5 ont des vitesses d'entrée différentes).
De plus les colonnes 2,3 et 4 sont entièrement nulles.
Il faut donc trouver comment déclarer la variable "effort" pour que le code que vous avez écrit puisse la reconnaître.

Je réfléchis de mon côté également. Je suis sûr que c'est tout bête.

**Gooby** · 07/01/2016, 17h43

Reprends mon code du post précédent, j'ai déclaré la variable effort à la première itération de jj

**tomazy4T6** · 07/01/2016, 18h15

Autant pour moi Gooby.
J'étais sur l'ancien code.
Je vous remercie bien.
Je vais essayer de le refaire de mon côté sans vous plagier.

**Gooby** · 07/01/2016, 18h26

Si tu veux comprendre le comportement du code (pourquoi le résultat n'est pas celui que tu attends), je te conseille de mettre un point d'arrêt, de l'exécuter ligne par ligne et de vérifier que le contenu de tes variables soient bons. Là, je vois pas pourquoi cela ne donne pas ce que tu veux, j'ai peut être fait une erreur mais je ne vois pas laquelle.