sbiosteadystate

Найдите устойчивое состояние модели SimBiology

Синтаксис

[success, variant_out]
= sbiosteadystate(model)

[success, variant_out]
= sbiosteadystate(model, variant_in)

[success, variant_out]
= sbiosteadystate(model, variant_in, scheduleDose)

 [success, variant_out, model_out]
= sbiosteadystate(model,___)

 [success, variant_out, model_out, exitInfo]
= sbiosteadystate(model,___)

[___] = sbiosteadystate(___, Name,Value)

Описание

пример

[success, variant_out] = sbiosteadystate(model) попытки найти устойчивое состояние модели SimBiology^®, model. Функция возвращает success, который является true если устойчивое состояние было найдено, и SimBiology Variant object, variant_out, со всеми непостоянными разновидностями, отсеками и параметрами модели, имеющей установившиеся значения. Если устойчивое состояние не было найдено, то success false и variant_out содержит последние значения, найденные алгоритмом.

пример

[success, variant_out] = sbiosteadystate(model, variant_in) применяет альтернативные значения количества, сохраненные в объекте варианта SimBiology, variant_in, к модели прежде, чем попытаться найти установившиеся значения.

пример

[success, variant_out] = sbiosteadystate(model, variant_in, scheduleDose) применяет объект дозы расписания SimBiology, scheduleDose, или вектор из доз расписания к соответствующим количествам модели прежде, чем попытаться найти значения устойчивого состояния. Только дозы во время = 0 позволены, то есть, время дозы каждого объекта дозы должно быть 0. Задавать дозу, не задавая различное, набор variant_in к пустому массиву, [].

пример

[success, variant_out, model_out] = sbiosteadystate(model,___) также возвращает модель SimBiology, model_out это - копия входа model с набором состояний к установившемуся решению, которое было найдено. Кроме того, model_out имеет все первоначальные отключенные правила присвоения.

пример

[success, variant_out, model_out, exitInfo] = sbiosteadystate(model,___) также возвращает информацию о выходе о расчете устойчивого состояния.

пример

[___] = sbiosteadystate(___, Name,Value) дополнительные опции использования заданы одним или несколькими Name,Value парные аргументы.

Примеры

свернуть все

Найдите устойчивое состояние простой модели регуляции генов

Скрипт Open Live Script

В этом примере показано, как найти устойчивое состояние простой модели регуляции генов, где продукт белка от перевода управляет записью.

Загрузите демонстрационный проект SimBiology, содержащий модель, m1. Модель имеет пять реакций и четыре разновидности.

sbioloadproject('gene_reg.sbproj','m1')

Отобразите реакции модели.

m1.Reactions

ans = 
   SimBiology Reaction Array

   Index:    Reaction:                    
   1         DNA -> DNA + mRNA            
   2         mRNA -> mRNA + protein       
   3         DNA + protein <-> DNA_protein
   4         mRNA -> null                 
   5         protein -> null

Вычисление устойчивого состояния пытается найти значения устойчивого состояния непостоянных количеств. Чтобы узнать, какие количества модели являются непостоянными в этой модели, используйте sbioselect.

sbioselect(m1,'Where','Constant*','==',false)

ans = 
   SimBiology Species Array

   Index:    Compartment:    Name:          Value:    Units:  
   1         unnamed         DNA            50        molecule
   2         unnamed         DNA_protein    0         molecule
   3         unnamed         mRNA           0         molecule
   4         unnamed         protein        0         molecule

Существует четыре разновидности, которые не являются постоянными, и начальные суммы трех из них обнуляются.

Используйте sbiosteadystate найти значения устойчивого состояния для тех непостоянных разновидностей.

[success,variantOut] = sbiosteadystate(m1)

success = logical
   1

variantOut = 
   SimBiology Variant - SteadyState (inactive)

   ContentIndex:     Type:        Name:             Property:           Value:
   1                 compartment  unnamed           Capacity            1
   2                 species      DNA               InitialAmount       8.79024
   3                 species      DNA_protein       InitialAmount       41.2098
   4                 species      mRNA              InitialAmount       1.17203
   5                 species      protein           InitialAmount       23.4406
   6                 parameter    Transcription.k1  Value               0.2
   7                 parameter    Translation.k2    Value               20
   8                 parameter    [Binding/Unbin... Value               0.2
   9                 parameter    [Binding/Unbin... Value               1
   10                parameter    [mRNA Degradat... Value               1.5
   11                parameter    [Protein Degra... Value               1

Начальные суммы всех разновидностей модели были установлены в установившиеся значения. DNA сохраненная разновидность начиная с общего количества DNA и DNA_protein равно 50.

Можно также использовать вариант, чтобы сохранить альтернативные начальные суммы и использовать их во время вычисления устойчивого состояния. Например, вы могли установить начальное количество ДНК к 100 молекулам вместо 50.

variantIn = sbiovariant('v1');
addcontent(variantIn,{'species','DNA','InitialAmount',100});
[success2,variantOut2,m2] = sbiosteadystate(m1,variantIn)

success2 = logical
   1

variantOut2 = 
   SimBiology Variant - SteadyState (inactive)

   ContentIndex:     Type:        Name:             Property:           Value:
   1                 compartment  unnamed           Capacity            1
   2                 species      DNA               InitialAmount       12.7876
   3                 species      DNA_protein       InitialAmount       87.2124
   4                 species      mRNA              InitialAmount       1.70502
   5                 species      protein           InitialAmount       34.1003
   6                 parameter    Transcription.k1  Value               0.2
   7                 parameter    Translation.k2    Value               20
   8                 parameter    [Binding/Unbin... Value               0.2
   9                 parameter    [Binding/Unbin... Value               1
   10                parameter    [mRNA Degradat... Value               1.5
   11                parameter    [Protein Degra... Value               1

m2 = 
   SimBiology Model - cell 

   Model Components:
     Compartments:      1
     Events:            0
     Parameters:        6
     Reactions:         5
     Rules:             0
     Species:           4
     Observables:       0

Поскольку алгоритм нашел устойчивое состояние, третий выход m2 модель устойчивого состояния, где значения непостоянных количеств были установлены к значениям устойчивого состояния. В этом примере начальные суммы всех четырех разновидностей были обновлены к значениям устойчивого состояния.

m2.Species

ans = 
   SimBiology Species Array

   Index:    Compartment:    Name:          Value:     Units:  
   1         unnamed         DNA            12.7876    molecule
   2         unnamed         DNA_protein    87.2124    molecule
   3         unnamed         mRNA           1.70502    molecule
   4         unnamed         protein        34.1003    molecule

Входные параметры

свернуть все

`model` — Модель SimBiology
Объект модели SimBiology

Модель SimBiology в виде SimBiology Model object.

`variant_in` — Вариант SimBiology
различный объект

Вариант SimBiology в виде a Variant object. Альтернативные значения количества, сохраненные в варианте, применяются к модели прежде, чем найти устойчивое состояние.

`scheduleDose` — Информация о дозах
Объект дозы расписания SimBiology

Информация о дозах в виде SimBiology планирует объект дозы. Доза должна быть шариком, то есть, не должно быть никакой задержки или время администрирования для дозы. Другими словами, его LagParameterName и DurationParameterName свойства должны быть пустыми, и время дозы (Time свойство), должен быть 0. Для получения дополнительной информации о том, как создать дозу шарика, смотрите Дозы Создания Программно.

Аргументы в виде пар имя-значение

Задайте дополнительные разделенные запятой пары Name,Value аргументы. Name имя аргумента и Value соответствующее значение. Name должен появиться в кавычках. Вы можете задать несколько аргументов в виде пар имен и значений в любом порядке, например: Name1, Value1, ..., NameN, ValueN.

Пример: 'AbsTol',1e-6 задает, чтобы использовать значение абсолютной погрешности ^10–6.

`'Method'` — Метод, чтобы вычислить устойчивое состояние
`'auto'` (значение по умолчанию) | `'simulation'` | `'algebraic'`

Метод, чтобы вычислить устойчивое состояние modelВ виде разделенной запятой пары, состоящей из 'Method' и вектор символов 'auto', 'simulation', или 'algebraic'. Значение по умолчанию ('auto') поведение должно использовать 'algebraic' метод сначала. Если тот метод неудачен, функция использует 'simulation' метод.

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

Для алгебраического метода функция вычисляет устойчивое состояние путем нахождения, что корень потока функционирует алгебраически. Для нелинейных моделей этот метод требует Optimization Toolbox™. Для получения дополнительной информации см. Алгебраический Метод.

Примечание

Устойчивое состояние, возвращенное алгебраическим методом, как гарантируют, не совпадет с тем, найденным методом симуляции. Алгебраический метод быстрее, поскольку он не включает симуляции, но метод симуляции может смочь найти устойчивое состояние, когда алгебраический метод не мог.

Пример: 'Method','algebraic'

`'AbsTol'` — Абсолютная погрешность, чтобы обнаружить сходимость
`1e-8` (значение по умолчанию) | положительный, действительный скаляр

Абсолютная погрешность, чтобы обнаружить сходимость в виде разделенной запятой пары, состоящей из 'AbsTol' и положительный, действительный скаляр.

Когда вы используете алгебраический метод, абсолютная погрешность используется, чтобы задать настройки оптимизации и обнаружить сходимость. Для получения дополнительной информации см. Алгебраический Метод.

Когда вы используете метод симуляции, абсолютная погрешность используется, чтобы определить сходимость при нахождении решения для устойчивого состояния прямым интегрированием можно следующим образом: $(‖ \frac{d \vec{S}}{d t} ‖ < A b s T o l) o r (‖ \frac{d \vec{S}}{d t} ‖ < R e l T o l * ‖ \vec{S} ‖)$ , где $\vec{S}$ вектор из непостоянных разновидностей, параметров и отсеков.

`'RelTol'` — Относительная погрешность, чтобы обнаружить сходимость
`1e-6` (значение по умолчанию) | положительный, действительный скаляр

Относительная погрешность, чтобы обнаружить сходимость в виде разделенной запятой пары, состоящей из 'RelTol' и положительный, действительный скаляр. Этот аргумент пары "имя-значение" используется для simulation метод только. Алгоритм сходится и сообщает об устойчивом состоянии, если алгоритм находит состояния модели прямым интегрированием, таким что $(‖ \frac{d \vec{S}}{d t} ‖ < A b s T o l) o r (‖ \frac{d \vec{S}}{d t} ‖ < R e l T o l * ‖ \vec{S} ‖)$ , где $\vec{S}$ вектор из непостоянных разновидностей, параметров и отсеков.

`'MaxStopTime'` — Максимальная сумма времени симуляции, чтобы взять прежде, чем завершить работу без устойчивого состояния
100000 (значение по умолчанию) | положительное целое число

Максимальная сумма времени симуляции, чтобы взять прежде, чем завершить работу без устойчивого состояния в виде разделенной запятой пары, состоящей из 'MaxStopTime' и положительное целое число. Этот аргумент пары "имя-значение" используется для simulation метод только.

`'MinStopTime'` — Минимальное количество времени симуляции, чтобы взять перед поиском устойчивого состояния
1 (значение по умолчанию) | положительное целое число

Минимальное количество времени симуляции, чтобы взять перед поиском устойчивого состояния в виде разделенной запятой пары, состоящей из 'MinStopTime' и положительное целое число. Этот аргумент пары "имя-значение" используется для simulation метод только.

Выходные аргументы

свернуть все

`success` — Отметьте, чтобы указать, найдено ли устойчивое состояние модели
`true` | `false`

Отметьте, чтобы указать, найдено ли устойчивое состояние модели, возвращено как true или false.

`variant_out` — Вариант SimBiology
различный объект

Вариант SimBiology, возвращенный как различный объект. Вариант включает все разновидности, параметры и отсеки модели с непостоянными количествами, имеющими установившиеся значения.

`model_out` — Модель SimBiology в устойчивом состоянии
объект модели

Модель SimBiology в устойчивом состоянии, возвращенном как объект модели. model_out копия входа model, с непостоянными разновидностями, параметрами и набором отсеков к установившимся значениям. Кроме того, model_out имеет все первоначальные отключенные правила присвоения. Симуляция модели в устойчивом состоянии требует, чтобы первоначальные правила присвоения были неактивны, поскольку эти правила могут изменить значения в variant_out.

Примечание

Если вы решаете фиксировать variant_out к входу model это имеет первоначальные правила присвоения, затем model как ожидают, не будет в устойчивом состоянии, потому что правила тревожат систему, когда вы симулируете model.
model_out в устойчивом состоянии только если симулирован без любых доз.

`exitInfo` — Выйдите из информации о расчете устойчивого состояния
символьный вектор

Выйдите из информации о расчете устойчивого состояния, возвращенном как вектор символов. Информация содержит различные сообщения для соответствующих выходных условий.

Steady state found (simulation) – Устойчивое состояние найдено с помощью метода симуляции.
Steady state found (algebraic) – Устойчивое состояние найдено с помощью алгебраического метода.
Steady state found (unstable) – Нестабильное устойчивое состояние найдено с помощью алгебраического метода.
Steady state found (possibly underdetermined) – Устойчивое состояние то есть, возможно, не асимптотически устойчивый найдено с помощью алгебраического метода.
No Steady state found – Никакое устойчивое состояние не найдено.
Optimization Toolbox (TM) is missing – Метод установлен в 'algebraic' поскольку нелинейные модели и Optimization Toolbox отсутствует.

Больше о

свернуть все

Метод симуляции

sbiosteadystate симулирует модель до MaxStopTime. Во время симуляции функция аппроксимирует градиент с помощью конечного дифференцирования (правая разность) в зависимости от времени, чтобы обнаружить устойчивое состояние.

Алгебраический метод

sbiosteadystate попытки найти устойчивое состояние модели алгебраически путем нахождения корня потока функционируют v. Функция потока включает уравнения реакции, правила скоростей, и алгебраические уравнения, то есть, v(X,P) = 0, где X и P являются непостоянными количествами и параметрами модели. Таким образом, массовое сохранение, наложенное уравнениями реакции, уважают.

Для нелинейных моделей, sbiosteadystate использование fmincon (Optimization Toolbox), чтобы получить исходное предположение для корня. Решение, найденное fmincon затем улучшен fsolve (Optimization Toolbox). Обнаружить сходимость, sbiosteadystate использует абсолютную погрешность ('AbsTol'). Другими словами, OptimalityTolerance, FunctionTolerance, и StepTolerance опции соответствующей оптимизационной функции установлены в 'AbsTol' значение.

Для линейных моделей, sbiosteadystate находит корни функции потока v путем решения линейной системы, заданной уравнениями реакции и сохранения. Для линейных моделей нет никакого уровня или алгебраических уравнений.

Введенный в R2016a

Документация

sbiosteadystate

Синтаксис

Описание

Примеры

Найдите устойчивое состояние простой модели регуляции генов

Входные параметры

`model` — Модель SimBiology
Объект модели SimBiology

`variant_in` — Вариант SimBiology
различный объект

`scheduleDose` — Информация о дозах
Объект дозы расписания SimBiology

Аргументы в виде пар имя-значение

`'Method'` — Метод, чтобы вычислить устойчивое состояние
`'auto'` (значение по умолчанию) | `'simulation'` | `'algebraic'`

`'AbsTol'` — Абсолютная погрешность, чтобы обнаружить сходимость
`1e-8` (значение по умолчанию) | положительный, действительный скаляр

`'RelTol'` — Относительная погрешность, чтобы обнаружить сходимость
`1e-6` (значение по умолчанию) | положительный, действительный скаляр

`'MinStopTime'` — Минимальное количество времени симуляции, чтобы взять перед поиском устойчивого состояния
1 (значение по умолчанию) | положительное целое число

Выходные аргументы

`success` — Отметьте, чтобы указать, найдено ли устойчивое состояние модели
`true` | `false`

`variant_out` — Вариант SimBiology
различный объект

`model_out` — Модель SimBiology в устойчивом состоянии
объект модели

`exitInfo` — Выйдите из информации о расчете устойчивого состояния
символьный вектор

Больше о

Метод симуляции

Алгебраический метод

Смотрите также

Документация SimBiology

Поддержка

Документация

sbiosteadystate

Синтаксис

Описание

Примеры

Найдите устойчивое состояние простой модели регуляции генов

Входные параметры

model — Модель SimBiology Объект модели SimBiology

variant_in — Вариант SimBiology различный объект

scheduleDose — Информация о дозах Объект дозы расписания SimBiology

Аргументы в виде пар имя-значение

'Method' — Метод, чтобы вычислить устойчивое состояние 'auto' (значение по умолчанию) | 'simulation' | 'algebraic'

'AbsTol' — Абсолютная погрешность, чтобы обнаружить сходимость 1e-8 (значение по умолчанию) | положительный, действительный скаляр

'RelTol' — Относительная погрешность, чтобы обнаружить сходимость 1e-6 (значение по умолчанию) | положительный, действительный скаляр

'MinStopTime' — Минимальное количество времени симуляции, чтобы взять перед поиском устойчивого состояния1 (значение по умолчанию) | положительное целое число

Выходные аргументы

success — Отметьте, чтобы указать, найдено ли устойчивое состояние модели true | false

variant_out — Вариант SimBiology различный объект

model_out — Модель SimBiology в устойчивом состоянии объект модели

exitInfo — Выйдите из информации о расчете устойчивого состояния символьный вектор

Больше о

Метод симуляции

Алгебраический метод

Смотрите также

Документация SimBiology

Поддержка

`model` — Модель SimBiology
Объект модели SimBiology

`variant_in` — Вариант SimBiology
различный объект

`scheduleDose` — Информация о дозах
Объект дозы расписания SimBiology

`'Method'` — Метод, чтобы вычислить устойчивое состояние
`'auto'` (значение по умолчанию) | `'simulation'` | `'algebraic'`

`'AbsTol'` — Абсолютная погрешность, чтобы обнаружить сходимость
`1e-8` (значение по умолчанию) | положительный, действительный скаляр

`'RelTol'` — Относительная погрешность, чтобы обнаружить сходимость
`1e-6` (значение по умолчанию) | положительный, действительный скаляр

`'MinStopTime'` — Минимальное количество времени симуляции, чтобы взять перед поиском устойчивого состояния
1 (значение по умолчанию) | положительное целое число

`success` — Отметьте, чтобы указать, найдено ли устойчивое состояние модели
`true` | `false`

`variant_out` — Вариант SimBiology
различный объект

`model_out` — Модель SimBiology в устойчивом состоянии
объект модели

`exitInfo` — Выйдите из информации о расчете устойчивого состояния
символьный вектор