Выполните скан параметра

Этот пример показывает, как выполнить скан параметра путем симуляции модели несколько раз, каждый раз изменяя значение параметра.

В модели, описанной в модели дрожжевого гетеротримерного цикла G-белка, скорость инактивации G-белка (kGd) намного ниже по мутантному штамму по сравнению со штаммом дикого типа (kGd = 0.004 от kGd = 0.11), что объясняет более высокие уровни активированного G-белка (Ga) в мутантном штамме. Для получения подробной информации о том, как изменяется уровень kGd влияет на уровень Ga, выполните скан параметров над различными значениями kGd.

Загрузите gprotein.sbproj проект, который включает переменную m1, объект модели.

sbioloadproject gprotein

Создайте вектор с пятью равномерно расположенными значениями для kGd в диапазоне от 0.001 на 0.15.

kGdValues = linspace(1e-3,0.15,5)';

Создайте SimFunction объект, где kGd - параметр входа для сканирования и Ga - наблюдаемый вид. Передайте в пустом массиве [] как последний входной параметр для обозначения отсутствия дозированных видов.

simfunc = createSimFunction(m1,{'kGd'},{'Ga'},[]);

Симулируйте модель несколько раз с различными значениями кГд. Установите время остановки равным 1000.

sd = simfunc(kGdValues,1000);

Постройте график результатов симуляции, чтобы увидеть, как изменяется уровень kGd влияет на уровень Ga.

sbioplot(sd);

Figure contains an axes. The axes with title States versus Time contains 5 objects of type line. These objects represent Run 1 - Ga, Run 2 - Ga, Run 3 - Ga, Run 4 - Ga, Run 5 - Ga.

См. также

|

Похожие темы