Can the Symbolic Math Toolbox Be Used to Prove Matrix Identities?
2 vues (au cours des 30 derniers jours)
Afficher commentaires plus anciens
Paul
le 30 Sep 2023
Commenté : Walter Roberson
le 1 Oct 2023
With the following defintions
syms s
A = symmatrix('A',3);
B = symmatrix('B',[3 1]);
C = symmatrix('C',[1 3]);
L = symmatrix('L',[3 1]);
I = sym(eye(3));
LHS = inv(s*I - A + L*C)*B + inv(s*I - A + L*C)*L*C*inv(s*I - A)*B;
RHS = inv(s*I - A)*B;
I'd like to prove the following identity
ident = LHS == RHS
Is there any other way to have the software verify that ident is always true besides subbing in symbolic values for the matrix elements?
isAlways(symmatrix2sym(ident))
0 commentaires
Réponse acceptée
Walter Roberson
le 30 Sep 2023
No.
5 commentaires
Walter Roberson
le 1 Oct 2023
Although children does not appear to be defined for symmatrix you can work around that using internal calls that are deliberately no longer documented.
syms s
A = symmatrix('A',3);
B = symmatrix('B',[3 1]);
C = symmatrix('C',[1 3]);
L = symmatrix('L',[3 1]);
I = sym(eye(3));
LHS = inv(s*I - A + L*C)*B + inv(s*I - A + L*C)*L*C*inv(s*I - A)*B;
RHS = inv(s*I - A)*B;
ident = LHS == RHS
feval(symengine, 'op', LHS, 0)
feval(symengine, 'op', LHS).'
feval(symengine, 'op', RHS, 0)
feval(symengine, 'op', RHS).'
whos
Walter Roberson
le 1 Oct 2023
syms s
A = symmatrix('A',3);
B = symmatrix('B',[3 1]);
C = symmatrix('C',[1 3]);
L = symmatrix('L',[3 1]);
I = sym(eye(3));
LHS = inv(s*I - A + L*C)*B + inv(s*I - A + L*C)*L*C*inv(s*I - A)*B;
RHS = inv(s*I - A)*B;
ident = LHS == RHS
feval(symengine, 'coerce', ident, 'DOM_STRING')
Plus de réponses (0)
Voir également
Catégories
En savoir plus sur Conversion Between Symbolic and Numeric dans Help Center et File Exchange
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!