我正在试图找出,如果我能找到一个矩阵指数,以某种形式的封闭形式的符号矩阵:a = symbols(['a'+str(1+k) for k in range(n)], real=True);我猜他们试图通过乔兰·布洛克获得矩阵的力量,而我的结构对此并不好。我怎么能继续下去呢?A**2是可以计算的,但是A**1和A**0的总和不能计算吗?!
我想创建一个块带对角线矩阵使用sympy。基本上,矩阵应该看起来像R D N 0 0 00 0 R D N 00 0 0 R D
该矩阵随后将通过numpy通过sympy.lambdify我确实找到了一种使用sympy.kronecker_product构造矩阵的方法,但是它非常慢。我看了一下为什么我的程序这么慢,发现kronecker.py是罪魁祸首。这是一个