Результаты моделирования печати на одном рисунке
sbioplot( выводит на график каждый прогон моделирования из sd)sd, a SimData объект или массив объектов на том же рисунке. График - это временной график каждого состояния в sd. На рисунке также показано иерархическое отображение всех прогонов как различных узлов в дереве, и можно выбрать, какой прогон будет отображаться.
sbioplot( печать результатов моделирования путем вызова дескриптора функции sd,fcnHandle,xArgs,yArgs,Name,Value)fcnHandle с входами sd, xArgs, и yArgsи использует дополнительные параметры, заданные одним или несколькими аргументами пары имя-значение. Например, можно указать x-метку и y-метку графика. xArgs и yArgs должны быть массивами ячеек или строковыми векторами имен состояний для построения графика.
Постройте график данных о добыче и хищнике из стохастически смоделированной модели лотки с использованием пользовательской функции (plotXY).
Загрузите модель. Установите тип решателя SSA для выполнения стохастического моделирования и установите время остановки равным 3.
sbioloadproject lotka; cs = getconfigset(m1); cs.SolverType = 'SSA'; cs.StopTime = 3; rng('default') % For reproducibility
Задание количества прогонов и использование sbioensemblerun для моделирования.
numRuns = 2; sd = sbioensemblerun(m1,numRuns);
Постройте график данных моделирования. По умолчанию sbioplot показывает график времени каждого вида для каждого прогона.
sbioplot(sd);
Постройте график выбранных состояний друг против друга; при этом постройте участок популяции добычи по сравнению с популяцией хищника. Используйте функцию plotXY (показано в конце этого примера) для построения графика смоделированных данных y1 (добыча) по сравнению с данными y2 (хищник). Укажите функцию как дескриптор функции.
Если для этого примера используется файл сценария в реальном времени, plotXY функция уже включена в конец файла. В противном случае необходимо определить plotXY в конце файла .m или .mlx или добавьте его в качестве файла по пути MATLAB.
sbioplot(sd,@plotXY,{'y1'},{'y2'},'xlabel','y1','ylabel','y2','title','Prey versus Predator');
Определение функции plotXY
sobioplot принимает дескриптор функции для функции с сигнатурой:
function [handles,names] = functionName(sd,xArgs,yArgs).
plotXY функция строит график двух выбранных состояний друг против друга. Первый вход sd - данные моделирования (SimBiology SimData объект или вектор объектов). В этом конкретном примере xArgs представляет собой клеточный массив, содержащий имя вида, подлежащего построению на оси x, и yArgs представляет собой клеточный массив, содержащий имя второго вида, подлежащего построению на оси y. Однако входные значения xArgs и yArgs можно использовать любым способом в пользовательской функции печати. Функция возвращает handlesмассив дескрипторов функций для линейных графиков, и namesмассив ячеек символьных векторов, показанных на узлах, которые являются дочерними узлами узла Run в иерархическом отображении.
function [handles,names] = plotXY(sd,xArgs,yArgs) % Select simulation data for each state from each run. xData1 = selectbyname(sd(1),xArgs); xData2 = selectbyname(sd(2),xArgs); yData1 = selectbyname(sd(1),yArgs); yData2 = selectbyname(sd(2),yArgs); % Plot the species against each other. fH1 = plot(xData1.Data,yData1.Data); fH2 = plot(xData2.Data,yData2.Data); % The first output, handles, is a two-dimensional array of handles of the line plots. It must be of size M x N, % where M is the number of line plots for each run and N is the number of runs. handles = [fH1,fH2]; % The second output, names, must be a one-dimensional cell array of character vectors. % Its length must be equal to the number of rows in handles, and the texts are displayed on the % nodes that are children of a Run node. names = {'y1 vs y2'}; end