Solving dynamic system of equations with parameters

7 views (last 30 days)
Hello,
I would like to solve a system of equations with parameters. An example (with trivial math) is attached. The system and the variables I would like to solve it for are dependent on a user input. In my example that is the varialbe "length" which can take an arbitrary positive natural number.
I can solve it by explicitly listing the variables "a(1)", etc. in the solve function call. Since the varialbe vector "a" is userdependent I cant use that in general though.
The second function call to solve results in warnings and wrong answers, because solve assumes that the second input "a" is an equation array and not the variables to be solved for.
Any suggestions to solve the system for, say, a random (positive and natural number) "length"?
Thanks a lot in advance!
length = 3;
a = sym('a', [length, 1]);
equations = sym('something', [length, 1]);
syms p
for runIdx = 1:length
equations(runIdx) = a(runIdx) - p == 0;
end
answers_correct_but_static = solve(equations, a(1), a(2), a(3))
answers_incorrect_dynamic = solve(equations, a)

Accepted Answer

A Jenkins
A Jenkins on 15 Jan 2014
Due to the way solve choses the variables to solve for ( symvar ), you may not need to specify the variables at all. For example, if you use 'x' as your variable name instead of 'a', this works in 2013b:
length = 3;
x = sym('x', [length, 1]);
equations = sym('something', [length, 1]);
syms p
for runIdx = 1:length
equations(runIdx) = x(runIdx) - p == 0;
end
answers_correct_dynamic = solve(equations)
  2 Comments
Sebastian
Sebastian on 15 Jan 2014
Edited: Sebastian on 15 Jan 2014
Hugh. Your suggestion works indeed.
I'm still confused though. I was aware of the fact that solve uses symvar to select the variables to solve for. Hence I used "a" as the variable name, because I figured it is selected first from symvar. Simply calling
answers = solve(equations);
In my above example doesn't work though.. How comes symvar picks the x's over the p? Or in other words: How robust is this solution?
A Jenkins
A Jenkins on 15 Jan 2014
Try this:
symvar(equations,length)
The documentation says it prefers variables "alphabetically closest to x."
Someone from TMW would have to comment on robustness, but I guess I agree that it seems natural to me that if you aren't sure, x,y,z are a good guess for variables and a,b,c are usually some kind of constant.

Sign in to comment.

More Answers (0)

Community Treasure Hunt

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

Start Hunting!