Siguiente: Análisis de Fourier Superior: Comenzando con PDL: un Anterior: Gráficos 2D

Operaciones con matrices

Para ilustrar un poco lo que se puede realizar con PDL veremos un ejemplo de cómo utilizar matrices.

Para multiplicar matrices podemos utilizar la notación siguiente

perldl> $a=pdl[ [1,2],[3,4] ]
perldl> $b=pdl[[5,6]]
perldl> $c=$a x $b
perldl> p $c
[
 [15 18]
 [35 42]
]

Para calcular la inversa necesitamos cargar previamente el módulo Slatec y utilizar el comando matinv de la siguiente forma

perldl> use PDL::Slatec
perldl> $inv=matinv($a)
perldl> p $inv
[
 [  -2    1]
 [ 1.5 -0.5]
]
perldl> p $inv x $a
[
 [1 0]
 [0 1]
]