用Mathematica计算一类特殊矩阵的特征多项式

1、先来构造这种矩阵的定义：f[n_] := Module[{A}, A = Table[0, n, n]; (A[[#, # + 1]] = 1) & /@ Range[n - 1]; (A[[# + 1, #]] = 1) & /@ Range[n - 1]; A]

3、再定义矩阵对应的特征多项式，g[n]。g[n_] := Module[{A}, A = Ta水瑞侮瑜ble[0荑樊综鲶, n, n]; (A[[#, # + 1]] = 1) & /@ Range[n - 1]; (A[[# + 1, #]] = 1) & /@ Range[n - 1]; Expand[Det[A - x*IdentityMatrix[n]]]]

5、但是，这种无止境的枚举法，于此题没有什么益处。我们很难通过这一串式子，观察出多项式的规律。不过，这些多项式之间存在一罩脾青槊种递归规律：Reduce[ForAll[x,(g[#]-(-x*g[#-1]+m*g[#-2])//Factor)&@9==0],m]实际上，把9改为任何大于2的正整数，都有m=-1。