Выберите данные моделирования из SimData объект с использованием выражений
___ = select( возвращает запрошенные данные моделирования в указанном формате данных.simdata,query,'Format',formatValue)
Загрузить модель глюкозно-инсулиновой реакции. Для получения подробной информации о модели см. раздел «Фон» в разделе «Моделирование реакции глюкоза-инсулин».
sbioloadproject('insulindemo.sbproj','m1');
Подавление информационного предупреждения, выдаваемого во время моделирования.
warnSettings = warning('off', 'SimBiology:DimAnalysisNotDone_MatlabFcn_Dimensionless');
Смоделировать один прием пищи для обычного субъекта в течение 7 часов.
singleMeal = sbioselect(m1,'Name','Single Meal'); cs = getconfigset(m1,'active'); cs.StopTime = 7; sd = sbiosimulate(m1,singleMeal)
SimBiology Simulation Data
ModelName: Cobelli's Glucose-Insulin System
Logged Data:
Species: 15
Compartment: 0
Parameter: 24
Sensitivity: 0
Observable: 0
sbioplot(sd);
Выберите все данные видов, зарегистрированные в SimData объект sd.
[t,x,names] = select(sd,{'Type','species'});
namesnames = 15x1 cell
{'Glucose appearance.Dose' }
{'Glucose appearance.Stomach Glu Solid' }
{'Glucose appearance.Stomach Glu Tritur'}
{'Glucose appearance.Stomach Glu' }
{'Glucose appearance.Gut Glu' }
{'Glucose appearance.Plasma Glu' }
{'Glucose appearance.Plasma Glu Conc' }
{'Glucose appearance.Tissue Glu' }
{'Insulin secretion.Interstitial Ins' }
{'Insulin secretion.Portal Ins' }
{'Insulin secretion.Liver Ins' }
{'Insulin secretion.Plasma Ins' }
{'Insulin secretion.Plasma Ins Conc' }
{'Insulin secretion.Ins Delay 1' }
{'Insulin secretion.Ins Delay 2' }
Данные графика для скорости появления глюкозы и использования глюкозы, а именно Glu Appearance Rate и Glu Util.
newsd = select(sd,{'Type','parameter','name',{'Glu Appear Rate'; 'Glu Util'}})
SimBiology Simulation Data
ModelName: Cobelli's Glucose-Insulin System
Logged Data:
Species: 0
Compartment: 0
Parameter: 2
Sensitivity: 0
Observable: 0
sbioplot(newsd);
Сравните данные по концентрации глюкозы в плазме (названные виды Plasma Glu Conc) и скорость секреции инсулина (параметр, названный Ins Secr). Использовать selectbyname извлечение данных путем указания соответствующих имен.
newsd2 = selectbyname(sd,{'Plasma Glu Conc','Ins Secr'})
SimBiology Simulation Data
ModelName: Cobelli's Glucose-Insulin System
Logged Data:
Species: 1
Compartment: 0
Parameter: 1
Sensitivity: 0
Observable: 0
sbioplot(newsd2);
Выберите данные для всех видов и параметров, которые имеют Glu в своих названиях.
newsd3 = select(sd,{'Where','Name','regexp','Glu'})
SimBiology Simulation Data
ModelName: Cobelli's Glucose-Insulin System
Logged Data:
Species: 7
Compartment: 0
Parameter: 11
Sensitivity: 0
Observable: 0
newsd3.DataNames
ans = 18x1 cell
{'Stomach Glu Solid' }
{'Stomach Glu Tritur' }
{'Stomach Glu' }
{'Gut Glu' }
{'Plasma Glu' }
{'Plasma Glu Conc' }
{'Tissue Glu' }
{'Stomach Glu After Dosing'}
{'Glu Appear Rate' }
{'Glu Prod' }
{'Plasma Glu Conc Rate' }
{'Ins Dep Glu Util' }
{'Glu Util' }
{'Glu Excretion' }
{'Glu Excretion Mode' }
{'Delayed Glu Signal' }
{'Delayed Glu Signal Mode' }
{'Basal Glu Prod' }
Можно также вернуть выбранные данные в виде структуры.
sdStruct = select(sd,{'Where','Name','regexp','Glu'},'Format','struct');Восстановите параметры предупреждения.
warning(warnSettings);