python计算矩阵方程_矩阵方程求解器Python

大多数象征性的工作都可以手工完成：X^2 - CX - I = 0

-> X^2 + 2EX - I = 0 // sub C = -2E

-> X^2 + 2EX + E^2 - I = E^2 //add E^2 to both sides (i.e., complete the square)

-> (X + E)^2 = E^2 + I //simplify and add I to both sides

-> X+E = +/-(E^2 + I)^(1/2) //take square root (now we may have more than one answer)

-> X = -E +/- (E^2 + I)^(1/2) //subtract E from both sides

矩阵平方根可能是，也可能不是你想象征性地解决的问题。SymPy当然会让你象征性地表示它，但事实证明它无法用数字来计算它，在我的尝试中(在MinGW64上的Python3中)。在

你的矩阵C是对称的，所以我们可以检查平方根下的项(即1/2)是否有明确的计算公式。一些初步事实：两个对称矩阵的和与差又是对称的。在

给定对称矩阵A和B，则AB是对称的当且仅当A和B交换。在

每个实对称矩阵都可以通过实正交相似性对角化。在对于任何可对角化矩阵A=VDV^(-1)，则A^(1/2) = VD^(1/2)V^(-1)。在

从C开始，我们问，E^2+I是可对角化的(因此有一个简单的矩阵平方根的显式公式)？

C是对称的，E = -(1/2)C；标量乘法不会改变{}的对称性，因为它会影响每个细胞；因此{}是对称的。

E^2 = (E * E)通勤，所以{}是对称的。

最后，I是对称的，所以(E^2 + I)是对称的。在

因此，可以使用上述(4)对角化的平方根。对角矩阵的平方根是通过取对角线上元素的平方根来计算的。在这里你可能会遇到另一个问题，如果这些元素是否定的，你的答案将是复杂的。对于每个平方根，也可能有多个答案，可能会给你一些答案来考虑。这很可能就是为什么SymPy没有给出一个数字答案。在