\( \DeclareMathOperator{\abs}{abs} \newcommand{\ensuremath}[1]{\mbox{$#1$}} \)
(%i3) |
A
:
matrix([
1,
2],[
0,
-
1])
$
B : matrix([ 0, - 1],[ 1, 3]) $ print( A, "y", B) $ |
Como base que constituyen de un su sistema generador, ortogonalizamosla:
(%i6) |
u1
:
A
$
u2 : B -(( mat_trace( transpose( u1). B)) / mat_trace(( transpose( u1). u1))) * u1 $ print( "Base ortogonal:", "{", u1, ",", u2, "}") $ |
El ejercicio nos pide la traza de la matriz producto:
(%i7) | mat_trace((( u1). u2)) ; |
Observemos que si elegimos u1 de otra forma el resultado será distinto:
(%i10) |
u1
:
B
$
u2 : A -(( mat_trace( transpose( u1). A)) / mat_trace(( transpose( u1). u1))) * u1 $ mat_trace((( u1). u2)) ; |
Created with wxMaxima.
Consideremos las matrices que debemos ortogonalizar: