Why does solve ignore if a variable is assumed real or positive?

5 Ansichten (letzte 30 Tage)
Prabhakar
Prabhakar am 18 Jan. 2011
When I use 'solve' on a complicated equation with variables assumed to be real or positive, these assumptions are sometimes ignored.
syms a positive;
solve('5*ln(a) - 32*exp(a/16) + 8*a^2 + 4*a^4=0','a')
ans = 0.22609034746836648974751899999113*i - 1.3900380484673922332261785846175

Antworten (1)

Prabhakar
Prabhakar am 18 Jan. 2011
This behaviour occurs because the solver cannot find an analytical solution, and switches to the numeric solver. The numeric solver ignores any assumptions made on variables. There are a couple of workarounds to this issue:
1. You can force the result to follow the desired assumption, for example:
abs(solve('5*ln(abs(a)) - 32*exp(abs(a)/16) + 8*abs(a)^2 + 4*abs(a)^4=0','a') )
2. You can access the numeric solver directly, giving the input it requires to find the desired result:
evalin(symengine, 'numeric::fsolve(5*ln(a) - 32*exp(a/16) + 8*a^2 + 4*a^4=0,a=0..infinity)')
  2 Kommentare
Walter Roberson
Walter Roberson am 20 Jan. 2011
Interestingly, Maple's solve() is not able to find a symbolic solution for the version with abs().
Christopher Creutzig
Christopher Creutzig am 24 Jan. 2011
Note that with the assumption that a be positive, your abs(a) is automatically simplified to a and your first command returns an incorrect solution (try subs('5*ln(a) ... 4*a^4', 'a', ans), I get -1.63.)

Melden Sie sich an, um zu kommentieren.

Tags

Community Treasure Hunt

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

Start Hunting!

Translated by