simBiology exported model simulation output not in the same order as `InitialValues`

1 次查看（过去 30 天）
Paulo
Paulo 2023-2-1
评论： Paulo 2023-2-3
Hello,
I'm using a SimBiology exported model and I need to get the last value of a simulation and use it as initial condition for a following simulation (with different parameters, or dosing information), here is an example of what I need to do:
[t, y, names] = simulate(exported_model);
exported_model.InitialValues = y;
[t, y, names] = simulate(exported_model);
The problem is in the second line of the code as, surprinsingly, the order of the variables in y is not the same as the order in InitialValues. Currently, I'm performing a search for the appropriate index in names and using `exported_model.getIndex` to update each element of InitialValues. But this approach is suboptimal, it takes more than 1min to go over the entire array of initial conditions.
Is there any way to ensure that the output of model simulations come in the correct order? Or what would be the optimal approach here?
Thank you
  1 个评论
Valentina Vasic
Valentina Vasic 2023-2-2
Hallo Paulo,
I have never tried it, I usually use your method, but I never needed to do double simulation. Maybe this can give you a hand. You have to go to the Examples.
An easier and faster method I don't think there is, but here is some advice.
Best,
Valentina

请先登录，再进行评论。

请先登录，再回答此问题。

采纳的回答

Fulden Buyukozturk
Ran in:
Hello,
exported_model.InitialValues lists all of exported_model's editable objects which is by default everything (compartments, species, parameters) whereas simulate(exported_model) returns states in the order of model's statesToLog, which is by default only non-constant quantities.
If you specify editable objects and statesToLog to be the same during model export, the order of quantities should be the same. Please see below for an example.
m = sbmlimport('lotka');
% define states for editable objects and statesToLog
states = sbioselect(m,'Name',{'y1','y2'});
% set statesToLog
m.getconfigset.RuntimeOptions.StatesToLog = states;
% export model specifying editable objects
em = export(m,states);
[t,y,names] = simulate(em);
{em.ValueInfo.Name}'
ans = 2×1 cell array
{'y1'} {'y2'}
names
names = 2×1 cell array
{'y1'} {'y2'}
Fulden

更多回答（1 个）

Fangjun Jiang
Fangjun Jiang 2023-2-1
编辑：Fangjun Jiang 2023-2-1
For a Simulink simulation, there are Inputs, Outputs and States. Usually, only States need initial values.
From "exported_model.InitialValues = y", it seems that you are assigning the Outputs of last simulation to the initial States of the next simulation. That is the problem. In general, Outputs will not be the same as the States, although they could be exactly the same.
There are ways to save the values of States and utilize it for the next simulation. See "States" and "Final states" in this doc.
web(fullfile(docroot, 'simulink/gui/data-import-export-pane.html'))
web(fullfile(docroot, 'simulink/ug/state-information.html'))
web(fullfile(docroot, 'simulink/ug/saving-and-restoring-simulation-operating-point.html'))
  3 个评论
显示 1更早的评论
Fangjun Jiang
Fangjun Jiang 2023-2-1
Based on the doc link below, it seems "names" in your code is the column labels of "y". I wonder if there is a column label for model.InitialValues.
If there is and the labels match but in a different order, then you could use ismember() to shuffle them to match.
https://www.mathworks.com/help/simbio/ref/simbiology.export.model.simulate.html
Paulo
Paulo 2023-2-1
Yes, that is my current approach (something similar to it anyway), I'm getting the index from names using:
[~, i] = max(cellfun(@(x) strcmp(x, name), all_names));
I do something similar using exported_model.getIndex, and I'm able to correctly update the initial values, but it takes more than one minute to do so. And that is even after I defined an array that maps one set of indexes into the other, so I don't have to the search every time.
So I was wondering if there is a more efficient way of doing it.
Thank you,

请先登录，再进行评论。

社区

更多回答在  SimBiology Community

类别

在 Help Center 和 File Exchange 中查找有关 Import Data 的更多信息

标签

产品


版本

R2020b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by