Clear Filters
Clear Filters

Encoding conditional assignments in SimBiology

23 views (last 30 days)
Is it possible to define conditional assignments in SimBiology, i.e.
if a > 0
param1 = a + b
else
param1 = a + c
end
I am aware of events but it is insuffient as it lacks the 'else' expression.

Accepted Answer

Ayush Singh
Ayush Singh on 15 Jul 2024
Hey emjey,
In SimBiology, conditional assignments like the one you described are not directly supported in the same way as typical programming languages. However, you can achieve similar functionality using repeated assignment rules. Since events do not support an `else` clause directly, you will need to set up your model in a way that effectively mimics the conditional logic.
You can achieve this using repeated assignment rules with logical conditions. SimBiology allows you to define rules that continuously update the value of a parameter based on other parameters or species in the model.
You can follow below steps to achieve similar functionality
1. Define Parameters:
- Define the parameters 'a', 'b', 'c', and 'param1' in your model.
2. Create a Rule:
- Add a repeated assignment rule to your model that uses logical conditions.
Example:
Suppose you have parameters 'a', 'b', 'c', and 'param1'. You can define a repeated assignment rule for 'param1' using logical conditions as follows:
% Create a SimBiology model
model = sbiomodel('conditionalAssignmentModel');
% Add parameters
a = addparameter(model, 'a', 1);
b = addparameter(model, 'b', 2);
c = addparameter(model, 'c', 3);
param1 = addparameter(model, 'param1', 0);
% Add a repeated assignment rule
% The rule will use logical conditions to mimic the if-else structure
addrule(model, 'param1 = (a > 0) * (a + b) + (a <= 0) * (a + c)', 'repeatedAssignment');
% Simulate the model (optional)
cs = getconfigset(model);
cs.StopTime = 10;
simData = sbiosimulate(model);
% Plot the results (optional)
sbioplot(simData);
Below is the explanation for above
The rule `param1 = (a > 0) * (a + b) + (a <= 0) * (a + c)` works as follows:
  • `(a > 0)` evaluates to `1` if `a` is greater than `0`, and `0` otherwise.
  • `(a <= 0)` evaluates to `1` if `a` is less than or equal to `0`, and `0` otherwise.
  • When `a > 0`, `(a > 0) * (a + b)` will be `a + b` and `(a <= 0) * (a + c)` will be `0`.
  • When `a <= 0`, `(a > 0) * (a + b)` will be `0` and `(a <= 0) * (a + c)` will be `a + c`.
Thus, `param1` will be assigned the value of `a + b` if `a > 0`, and `a + c` if `a <= 0`.
Hope it helps!

More Answers (1)

Arthur Goldsipe
Arthur Goldsipe on 15 Jul 2024
Hi emjey,
I see this problem solved in one of two ways:
1) Take advantage of true/false getting converted to 1/0 in mathematical operations. This is the approach suggested in the other answer. Another altenative for the rule expression is "param1 = a + b*(a>0) + c*(a<=0)". One disadvantage of this approach is that the rule expressions can easily get complicated and hard to debug, especially for more complicated "if, else if, else" expressions.
2) You can put any complicated logic, include if/else in a helper function that you call from the repeated assignment. That way, you can use if/else to write the code the most natural way. You ca find an example of that approach here. The main downside to this is that you must share this function file with anyone else who wants to use the model, since it isn't automatically bundled with the model or project file.

Categories

Find more on Perform Sensitivity Analysis in Help Center and File Exchange

Products


Release

R2022a

Community Treasure Hunt

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

Start Hunting!