Sympy: Exp (-iHt): ошибка Mutabledensematrix

Question

Этот вопрос следует за этим ( Sympy: Exp (-iHt) ). Я перепробовал много глупостей с этой первой проблемой: без комментариев.

Моя настоящая проблема заключается в следующем:

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)

Я получаю эту ошибку: Ошибка типа: нечитаемый тип: 'MutableDenseMatrix'

Как я могу решить это?

Спасибо за ответ.

Answer 1

Ваш код должен работать, но любой способ обойти это - вычислить экспоненту самостоятельно из формы Иордана:

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                          ℯ      ⎦