- `(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`.
Encoding conditional assignments in SimBiology
23 views (last 30 days)
Show older comments
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.
0 Comments
Accepted Answer
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:
Thus, `param1` will be assigned the value of `a + b` if `a > 0`, and `a + c` if `a <= 0`.
Hope it helps!
0 Comments
More Answers (1)
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.
0 Comments
See Also
Categories
Find more on Perform Sensitivity Analysis in Help Center and File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!