mathematica 计算分块矩阵的特征值
问题描述：
mathematica 计算分块矩阵的特征值
比如A,B,C,D,均为n阶方阵,如何计算矩阵{{A,B},{C,D}的特征值?我直接输入总是报错啊.
比如，我输入
a = ({{1,2},
{3,4}});
b = ({{1,1},
{2,4}});
c = ({{1,4},
{3,2}});
d = ({{1,3},
{4,2}});
AD = ({
{a,b},
{c,d}
});
Det[AD]
输出提示为
Det::matsq:Argument {{{{1,2},{3,4}},{{1,1},{2,4}}},{{{1,4},{3,2}},{{1,3},{4,2}}}} at position 1 is not a non-empty square matrix.>>
哪里有错呢？
答
你这样输入不是分块矩阵了,而是一个由四个二姐矩阵组成的二阶矩阵.你应该这样：a = ({{1,2},{3,4}});b = ({{1,1},{2,4}});c = ({{1,4},{3,2}});d = ({{1,3},{4,2}});AD = ArrayFlatten[({{a,b},{c,d}})];Det[AD]...