From: Subject: Math::Matrix - Multiply and invert Matrices Date: Fri, 5 Aug 2005 16:42:21 +0900 MIME-Version: 1.0 Content-Type: multipart/related; type="text/html"; boundary="----=_NextPart_000_0009_01C599DC.A6ACBB40" X-MimeOLE: Produced By Microsoft MimeOLE V6.00.2900.2670 This is a multi-part message in MIME format. ------=_NextPart_000_0009_01C599DC.A6ACBB40 Content-Type: text/html; charset="shift_jis" Content-Transfer-Encoding: quoted-printable Content-Location: http://www4.kcn.ne.jp/~felix/Softwares/MathPPM/Math-Matrix.html Math::Matrix - Multiply and invert = Matrices
Location : Home > Softwares = > Perl = Math=20 Packages
Title : Math::Matrix

3D"Toolbox

Math::Matrix - Multiply and invert=20 Matrices

=96=BC=8F=CC

Math::Matrix - =8Ds=97=F1=82=CC=90=CF=82=C6=8Bt=8Ds=97=F1

=8AT=97v

use Math::Matrix;

=90=E0=96=BE

=81@=88=C8=89=BA=82=CC=83=81=83\=83b=83h=82=AA=97=98=97p=82=C5=82=AB=82= =E9=81B

new

=81@=8EQ=8F=C6=82=B3=82=EA=82=BD=94z=97=F1=82=C6=93=AF=82=B6=92=B7=82=B3= =82=CC=8Ds=97=F1=82=F0=90=B6=90=AC=82=B7=82=E9=81B=94z=97=F1=82=CD=83R=83= s=81[=82=B3=82=EA=82=E9=81B=83G=83=89=81[=82=CC=8F=EA=8D=87 = <undef> =82=F0=95=D4=82=B7=81B

        $a =3D new =
Math::Matrix ([rand,rand,rand],
                               [rand,rand,rand],
                               [rand,rand,rand]);

=81@<new> = =82=F0=83=81=83\=83b=83h=82=C6=82=B5=82=C4=8C=C4=82=D1=8Fo=82=B5=82=BD=8F= =EA=8D=87=81A=8Ae=97v=91f=82=F0 0 = =82=C5=96=9E=82=BD=82=B5=82=BD=8Ds=97=F1=82=F0=95=D4=82=B7=81B

clone

=81@=82=DC=82=C1=82=BD=82=AD=93=AF=82=B6=8Ds=97=F1=82=F0=90=B6=90=AC=82= =B7=82=E9=81B

	$b =3D $a->clone;

size

=81@=8Ds=97=F1=82=CC=83T=83C=83Y=82=F0=8C=88=92=E8=82=B7=82=E9=81B

= 
	($m, $n) =3D $a->size;

concat

=81@=93=AF=82=B6=8Ds=82=CC=90=94=82=CC=8Ds=97=F1=82=F0=98A=8C=8B=82=B7= =82=E9=81B=8C=8B=89=CA=82=CD=90V=82=B5=82=A2=8Ds=97=F1=82=C5=95=D4=82=B7=82= =AA=81A=83G=83=89=81[=82=CC=8F=EA=8D=87=82=CD <undef> = =82=F0=95=D4=82=B7=81B 

	$b =3D new Math::Matrix =
([rand],[rand],[rand]);
	$c =3D $a->concat($b);

transpose

=81@=93]=92u=8Ds=97=F1=82=F0=95=D4=82=B7=81B=82=B1=82=EA=82=CD=8Ds=97=F1= =82=CC=8Ds=82=C6=97=F1=82=CC=92l=82=AA=8C=F0=8A=B7=82=B3=82=EA=82=BD=82=E0= =82=CC=82=C5=82=A0=82=E9=81B

multiply

=81@=82P=82=C2=82=DF=82=CC=8Ds=97=F1=82=CC=8Ds=82=CC=92=B7=82=B3=82=C6= =82Q=82=C2=82=DF=82=CC=8Ds=97=F1=82=CC=97=F1=82=CC=92=B7=82=B3=82=AA=88=EA= =92v=82=B7=82=E9=82=C6=82=AB=82=C9=82Q=82=C2=82=CC=8Ds=97=F1=82=CC=90=CF=82= =F0=8B=81=82=DF=82=E9=81B=90=CF=82=CC=8C=8B=89=CA=82=F0=95=D4=82=B7=82=AA= =81A=83G=83=89=81[=82=CC=8F=EA=8D=87=82=CD <undef>=20 =82=F0=95=D4=82=B7=81B

solve

=81@=8Ds=97=F1=82=C5=97^=82=A6=82=E7=82=EA=82=BD=97=F5=97=9D=92m=95=FB= =92=F6=8E=AE=82=F0=89=F0=82=AD=81B=8Ds=82=CC=90=94=82=CD=97=F1=82=CC=90=94= =82=E6=82=E8=82=E0=91=E5=82=AB=82=AD=82=C8=82=AF=82=EA=82=CE=82=C8=82=E7=82= =C8=82=A2=81B=95=CF=90=94=82=AA=8C=DD=82=A2=82=C9=93=C6=97=A7=82=C5=82=C8= =82=AF=82=EA=82=CE=82Q=94=D4=96=DA=88=C8=8D~=82=CC=95=CF=90=94=82=C9=91=CE= =82=B7=82=E9=8CW=90=94=82=CD 0=20 =82=C6=82=C8=82=E9=81B=89=F0=82=F0=8A=DC=82=DE=8Ds=97=F1=82=F0=95=D4=82=B7= =82=AA=81A=83G=83=89=81[=82=CC=8F=EA=8D=87=82=CD <undef> = =82=F0=95=D4=82=B7=81B 





multiply_scalar


=81@=8Ds=97=F1=82=CC=83X=83J=83=89=81[=94{=82=F0=8Cv=8EZ=82=B7=82=E9=81=
B 
  $a->multiply_scalar(2);  scale matrix by factor 2





add


=81@=93=AF=8E=9F=82=CC=82Q=82=C2=82=CC=8Ds=97=F1=82=CC=98a=82=F0=8B=81=
=82=DF=82=E9=81B


substract


=81@add($other->negative)=82=CC=97=AA=8BL=81B


equal


=81@=82Q=82=C2=82=CC=8Ds=97=F1=82=AA=93=99=82=B5=82=A2=82=A9=94=DB=82=A9=
=82=F0=94=BB=92f=81B=8A=DB=82=DF=8C=EB=8D=B7=82=C9=92=8D=88=D3=82=B9=82=E6=
=81B


slice


=81@=97=F1=81icolumns=81j=82=F0=92=8A=8Fo=82=B7=82=E9=81B 
	=
a->slice(1,3,5);





determinant


=81@=8Ds=97=F1=82=CC=8Ds=97=F1=8E=AE=82=F0=8B=81=82=DF=82=E9=81B


dot_product


=81@=82Q=82=C2=82=CC=83x=83N=83g=83=8B=82=CC=93=E0=90=CF=82=F0=8B=81=82=
=DF=82=E9=81B


absolute


=81@=83x=83N=83g=83=8B=82=CC=90=E2=91=CE=92l=82=F0=8B=81=82=DF=82=E9=81=
B


normalizing


=81@=83x=83N=83g=83=8B=82=F0=90=B3=8BK=89=BB=82=B7=82=E9=81B


cross_product


=81@=83x=83N=83g=83=8B=82=CC=8AO=90=CF=82=F0=8B=81=82=DF=82=E9=81B
=


print


=81@=8Ds=97=F1=82=F0=95W=8F=80=8Fo=97=CD=82=C9=8Fo=97=CD=82=B7=82=E9=81=
B=92=C7=89=C1=93I=82=C8=83p=83=89=83=81=81[=83^=82=AA=82=A0=82=EA=82=CE=8D=
s=97=F1=82=CC=91O=82=C9=8Fo=97=CD=82=B3=82=EA=82=E9=81B





=8Eg=97p=97=E1


        use Math::Matrix;

        srand(time);
        $a =3D new Math::Matrix ([rand,rand,rand],=20
                         [rand,rand,rand],=20
                         [rand,rand,rand]);
        $x =3D new Math::Matrix ([rand,rand,rand]);
        $a->print("A\n");
        $E =3D $a->concat($x->transpose);
        $E->print("Equation system\n");
        $s =3D $E->solve;
        $s->print("Solutions s\n");
        $a->multiply($s)->print("A*s\n");






=92=98=8E=D2


Ulrich Pfeifer =
<pfeifer@ls6.informatik.uni-dortmund.de>
Brian J.=20
Watson <bjbrew@power.net>
Matthew Brett=20
<matthew.brett@mrc-cbu.cam.ac.uk>



  
  

    Math::Matrix - Multiply and invert=20
      Matrices 





Update : 2003/04/30



  
  

    
      
3D"Toolbox

    Update :=20
      2003/04/3


------=_NextPart_000_0009_01C599DC.A6ACBB40
Content-Type: image/png
Content-Transfer-Encoding: base64
Content-Location: http://www4.kcn.ne.jp/~felix/img/Toolbox.png

iVBORw0KGgoAAAANSUhEUgAAAIIAAAA3CAIAAABPZgVsAAAAAXNSR0IArs4c6QAAAARnQU1BAACx
jwv8YQUAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAAt5JREFU
eF7tml2WgjAMhZkt6c54cTE8uA+Oy9BN6CIUWwht6V9KIZGJx6c5MFzul6Rp6l9zbuTDwIEBg3zJ
HRAGLBxgIYI8GMkFCAYWDrAQQR6M5AIEAwsHWIggD0ZyAYKBhQMsRJAHI7kAwcDCARYiyIORXIBg
YOEACxHkwUguQDCwcICFCPJgJBcgGFg44BXR3p7v0s/9SnOA0T2U4kdnvdGpu6s/U6nKZey57tT2
r1II6pW70+4kJrvdp181hQUe8irkCKieDYIhNwNMErh7Ln1GtZJswBcDHIY59yNVSzAwwECzGB5y
bQimiJENz1uLy6RNV8V/i4Em6kMsBQOLnEhiePXt/m10fgHAmWgUpWQ24PaA07ar7K4miSG1EYLX
KRUwbrZcW2DjogQEK/lWGLJaW9ObZ39pmrK7hqBbjWHc7hQLgBvNbbxnIxxoI9lgUPMGtAvTlKIC
hrI4mASAcrP6efp7KgzxomzEi4UBdRcqGwJDp6FgWOmIFfDbGM7NOJLT0zfvy/gWOuuuVRhggIbB
4MgeBAiGb12d01+5CeTsFdI3gq2NQa3Eo4ZlUQol2eZrQ6pTrJ0N77f5xN0wnBunxfo2RS4GO0TK
R3v5DWtZeSm7y8oG1X0BiR0xfAM60mLEw/GA2WAeMOyMQQU45Pfcjzt9wXK1OywGHX37Y/APoVNT
58Ni0IeAO2OI7nticwfBMJ34rmtYF0u0PhV3RyOh6iQY6mBY5gE4bpYppg1r/V00rIuZRWn86Ua1
XTQ83xouwVOohhmp0eZsW52ZUgpDUE+lmdKPYyhzAfarix/16PRHzApVnCKu135PAiJFKTJrgc4V
tzbMa06qA1uuTvG0WHneMLyP3RRBW+Lp4n1K4CSg8LzBfV/XHy2j0nlD/nGSXIlyAJsNcv0mDmzy
T1GBIBcPDggGFg6wECEJIRhYOMBChGSDYGDhAAsRkg2CgYUDLERINggGDg58ACbUQ7InrI1DAAAA
AElFTkSuQmCC

------=_NextPart_000_0009_01C599DC.A6ACBB40
Content-Type: image/png
Content-Transfer-Encoding: base64
Content-Location: http://www4.kcn.ne.jp/~felix/img/greenTip.png

iVBORw0KGgoAAAANSUhEUgAAAAkAAAAHCAIAAABV+fA3AAAAAXNSR0IArs4c6QAAAARnQU1BAACx
jwv8YQUAAAAgY0hSTQAAeiYAAICEAAD6AAAAgOgAAHUwAADqYAAAOpgAABdwnLpRPAAAABhJREFU
GFdjZDBmwA2AcrgQTgmQefTVBwCwAwr3/S7ywgAAAABJRU5ErkJggg==

------=_NextPart_000_0009_01C599DC.A6ACBB40
Content-Type: text/css;
	charset="iso-2022-jp"
Content-Transfer-Encoding: quoted-printable
Content-Location: http://www4.kcn.ne.jp/~felix/base.css

.engLine {
	FONT-FAMILY: Arial
}
.engLine_w {
	COLOR: white; FONT-FAMILY: Arial
}
.engLine_dg {
	COLOR: darkgreen; FONT-FAMILY: Arial
}
.engTitle {
	FONT-WEIGHT: bold; FONT-SIZE: 24pt; COLOR: darkgreen; FONT-FAMILY: =
Times New Roman
}
.engTitle2 {
	FONT-WEIGHT: bold; FONT-SIZE: 20pt; COLOR: darkgreen; FONT-FAMILY: =
Times New Roman
}
.engTitle3 {
	FONT-WEIGHT: bold; FONT-SIZE: 16pt; COLOR: darkgreen; FONT-FAMILY: =
Times New Roman
}
.jpTitle {
	FONT-WEIGHT: bold; FONT-SIZE: 24pt; COLOR: darkgreen; FONT-FAMILY: MS =
Gothic
}
.jpTitle2 {
	FONT-WEIGHT: bold; FONT-SIZE: 20pt; COLOR: darkgreen; FONT-FAMILY: MS =
Gothic
}
.jpTitle3 {
	FONT-WEIGHT: bold; FONT-SIZE: 16pt; COLOR: darkgreen; FONT-FAMILY: MS =
Gothic
}
.sigWord {
	FONT-WEIGHT: bold; COLOR: darkgreen
}
S.prgCode {
	FONT-SIZE: 10.5pt; FONT-FAMILY: Courier,Arial
}

------=_NextPart_000_0009_01C599DC.A6ACBB40--