Discussion :

[tomazy4T6]

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]

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]

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]

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]

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]

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?