2024-04-19 21:14:35 发布
网友
这个问题跟在这个问题后面。我尝试了许多愚蠢的事情来解决第一个问题:无可奉告
我真正的问题是:
t=sympy.Symbol('t') H=sympy.Matrix([[3, 0, 0, 0], [0, 1 , 2, 0], [0, 2, 2, 0], [0, 0, 0, 4]]) A = sympy.exp(-sympy.I * H * t)
我得到了这个错误:TypeError:unhabable类型:“MutableDenseMatrix”
我怎样才能解决它
谢谢你的回答
您的代码应该可以工作,但解决方法是自己从Jordan表单计算指数:
In [11]: P, J = H.jordan_form() In [12]: P*(-I*t*J).exp()*P.inv() Out[12]: ⎡ -3⋅ⅈ⋅t ⎤ ⎢ℯ 0 0 0 ⎥ ⎢ ⎥ ⎢ ⎛3 √17⎞ ⎛3 √17⎞ ⎛3 √17⎞ ⎛3 √17⎞ ⎥ ⎢ -ⅈ⋅t⋅⎜─ + ───⎟ -ⅈ⋅t⋅⎜─ - ───⎟ -ⅈ⋅t⋅⎜─ + ───⎟ -ⅈ⋅t⋅⎜─ - ───⎟ ⎥ ⎢ ⎝2 2 ⎠ ⎝2 2 ⎠ ⎝2 2 ⎠ ⎝2 2 ⎠ ⎥ ⎢ 4⋅√17⋅ℯ 2⋅ℯ 16⋅ℯ 16⋅ℯ ⎥ ⎢ 0 - ───────────────────── + ───────────────────────────────────── ─────────────────────── - ────────────────────── 0 ⎥ ⎢ ⎛ √17 1⎞ ⎛ 1 √17⎞ ⎛ 2 2 ⎞ ⎛ √17 1⎞ ⎛ 1 √17⎞ ⎥ ⎢ 17⋅⎜- ─── - ─⎟ ⎜- ─ + ───⎟⋅⎜- ───────── + ─────────⎟ (-17 + √17)⋅⎜- ─── - ─⎟ ⎜- ─ + ───⎟⋅(√17 + 17) ⎥ ⎢ ⎝ 2 2⎠ ⎝ 2 2 ⎠ ⎜ √17 1 1 √17⎟ ⎝ 2 2⎠ ⎝ 2 2 ⎠ ⎥ ⎢ ⎜ - ─── - ─ - ─ + ───⎟ ⎥ ⎢ ⎝ 2 2 2 2 ⎠ ⎥ ⎢ ⎥ ⎢ ⎛3 √17⎞ ⎛3 √17⎞ ⎛3 √17⎞ ⎛3 √17⎞ ⎥ ⎢ -ⅈ⋅t⋅⎜─ + ───⎟ -ⅈ⋅t⋅⎜─ - ───⎟ -ⅈ⋅t⋅⎜─ + ───⎟ -ⅈ⋅t⋅⎜─ - ───⎟ ⎥ ⎢ ⎝2 2 ⎠ ⎝2 2 ⎠ ⎝2 2 ⎠ ⎝2 2 ⎠ ⎥ ⎢ 2⋅√17⋅ℯ ℯ 8⋅ℯ 8⋅ℯ ⎥ ⎢ 0 ───────────────────── - ─────────────────────── - ───────────────── + ───────────────── 0 ⎥ ⎢ 17 2 2 -17 + √17 √17 + 17 ⎥ ⎢ - ───────── + ───────── ⎥ ⎢ √17 1 1 √17 ⎥ ⎢ - ─── - ─ - ─ + ─── ⎥ ⎢ 2 2 2 2 ⎥ ⎢ ⎥ ⎢ -4⋅ⅈ⋅t⎥ ⎣ 0 0 0 ℯ ⎦
您的代码应该可以工作,但解决方法是自己从Jordan表单计算指数:
相关问题 更多 >
编程相关推荐