How do I get subs to behave as desired for symbolic expressions ... does collect change an expression?

1 visualización (últimos 30 días)
I am trying to simplify a large symbolic expression and am running into a problem.
The problem is represented by the following code:
syms a b x
f = a*x + b*x + a*b*x*x
subs(f,a+b,'c')
subs(collect(f,x),a+b,'c')
Here, I see two different behaviors from subs. In the first, the substitution of a+b with c is ignored and the original expression is returned. In the second, the desired behavior occurs and the resulting expression replaces a*x+b*x with c*x. Is there something I am missing -- how can I get the first subs to work?
In a related issue, I am finding that the follwing code yields baffling results:
isequal(f,f)
isequal(f,collect(f,x))
The first equality check returns true -- as expected since f should equal itself. However, the second equality check returns false. Is this expected behavior, and if so, how does someone check for equality between expressions?

Respuesta aceptada

Walter Roberson
Walter Roberson el 17 de Ag. de 2019
collect(f,x) is a different expression with a different internal symbol than f is.
isAlways(f==collect(f,x))
It is tempting to use
simplify(f-collect(f,x))
but in sufficiently complicated expressions, simplify() will not be able to prove that the value is 0. For example, if you rewrite complex values into exp and sincos, then the transform is too advanced for simplify() to detect equality.
  3 comentarios
Brandon Laflen
Brandon Laflen el 22 de Ag. de 2019
Thanks Walter, these are useful tips. Especially, using subs(f,a,c-b) makes it work as desired when I follow with a simplify.
Do you happen to know if isAlways(f==collect(f,x)) is guaranteed to always return true, regardless of f?
Walter Roberson
Walter Roberson el 22 de Ag. de 2019
No, there is a known case where it will fail: f = sym(nan) then isAlways will return false because NaN == NaN is considered to be false.
I do not know if there are expressions that are too complicated for isAlways to prove.
You should watch out for the cases that isAlways cannot prove; see https://www.mathworks.com/help/symbolic/isalways.html

Iniciar sesión para comentar.

Más respuestas (0)

Productos


Versión

R2017a

Community Treasure Hunt

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

Start Hunting!

Translated by