如何用Mathematica求解相合转换矩阵剖析.docx

Shanglong Ning University of Chinese Academy of Sciences  DATE \@ dddd, MMMM d, y 星期一, 六月 5, 17 ALL RIGHTS RESERVED 如何用Mathematica求解相合转换矩阵 Mathematica介绍 Mathematica是一款功能强大的数学软件，在各个学科中均有广泛应用。如果你对此感兴趣，可以参考维基百科或百度对此的介绍。废话不多说，下面进入正题：如何用Mathematica求解相合转换矩阵？ 相合矩阵的定义： 在数学上，所谓矩阵相合是指，在某一域上的矩阵A和B，如果存在一可逆矩阵（非奇矩阵）P使得PTAP=B,则称A与B矩阵相合。相合关系是一种等价关系，因为它满足自反性，对称性与传递性。 下面举一个例子说明在已知矩阵A，B的情况下如何求解转换矩阵P。 A=012210379 ， B=41438-42在mathematica中操作如下： 1、首先定义已知矩阵A,B与未知矩阵P(注意加了;可以使结果不显示出来，一切都请使用英文标点)，其中a11,a12,a13,a21,a22,a23,a31,a32,a32为未知数 2、然后开始解方程，请注意这里矩阵相乘用dot(.)然后Flatten是必须的。本质上讲，这里是把它换成线性方程组求解，这里是核心