Parametric solutions to system of equations inequations via solve()

3 Ansichten (letzte 30 Tage)
Arthur
Arthur am 8 Nov. 2014
Kommentiert: Arthur am 13 Nov. 2014
Given a square matrix A in n x n symbolic variables a_{i,j}, I want a parametric solutions for the a_{i,j} constrained to a system of equations and inequalities.
For instance, given
syms a, b, c, d, x
A = [[a, b]; [c, d]]
e(x) = charpoly(A, x)
I would like to find parametric solutions for all a, b, c, d that satisfy, for example
e(1)==0 and a >= 0, b >= 0, c >= 0, d >= 0.
Maple has functionality like this ("Solve semi-algebraic", or something) and I would like to be able to do this with Matlab.
I attempted to use solve() (see attached screenshot), but, even in this simpler case I consider there, I get the warning and the empty solution, so I'm guessing that I'm not understanding the solve() function correctly.
Thank you.

Melden Sie sich an, um diese Frage zu beantworten.

Akzeptierte Antwort

Stefan Wehmeier
Stefan Wehmeier am 11 Nov. 2014
I see that you are using R2014a. Then your only option is to solve for one variable
solve(e(1) == 0, a)
and take this as a parameterization with which you can play by plugging in values for two variables, and solving result>=0 for the third.
If you upgrade to R2014b, just
solve([e(1) == 0, a>=0, b>=0, c>=0, d>=0], [a, b, c, d])
works and gives you some solutions.
  1 Kommentar
Arthur
Arthur am 13 Nov. 2014
Alas, there is no student edition of R2014b. Strange that this feature would be unavailable until that release.
Anyway, thanks for your answer.

Melden Sie sich an, um zu kommentieren.

Weitere Antworten (0)

Kategorien

Mehr zu Equation Solving finden Sie in Help Center und File Exchange

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by