我正在嘗試在 sympy 中計算哈密頓時間演化。H 是一個矩陣（哈密頓矩陣）。例如 ：Matrix([[1, 2], [2, 2]]t 是一個符號：t=Symbol('t')如何輕松計算 A = exp (-iHt) ？我嘗試了很多事情但沒有成功：計算特征值，特征向量......
2 回答
叮當貓咪
TA貢獻1776條經驗 獲得超12個贊
應該能夠做這樣的事情：
import sympy
t = sympy.symbols('t')
H = sympy.Matrix([[1, 2], [2, 2]])
A = sympy.exp(-sympy.I * H * t) # sympy.I is imaginary constant
A按預期返回一個復雜的 sympy 矩陣。
森欄
TA貢獻1810條經驗 獲得超5個贊
這是你想要的嗎？
In [1]: M = Matrix([[1, 2], [2, 2]])
In [2]: (I*M*t).exp()
Out[2]:
? 3???t √17???t 3???t √17???t √17???t 3???t √17???t 3???t 3???t √17???t 3???t √17???t √17???t 3???t √17???t 3???t?
? ───── + ─────── ───── + ─────── - ─────── + ───── - ─────── + ───── ───── + ─────── ───── + ─────── - ─────── + ───── - ─────── + ─────?
? 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 ?
? 4?√17?? + 68?? + 13?√17?? + 85?? 2?√17?? + 18?? - 18?? - 2?√17?? ?
? ──────────────────────────────────────────────────────────────────────────────────────────────── ───────────────────────────────────────────────────────────────────────────────────────────────?
? 17?√17 + 153 17 + 9?√17 ?
? ?
? ? 3???t √17???t 3???t √17???t √17???t 3???t √17???t 3???t? 3???t √17???t 3???t √17???t √17???t 3???t ?
? ? ───── + ─────── ───── + ─────── - ─────── + ───── - ─────── + ─────? ───── + ─────── ───── + ─────── - ─────── + ───── ?
? ? 2 2 2 2 2 2 2 2 ? 2 2 2 2 2 2 ?
?-?- 18?√17?? - 34?? + 34?? + 18?√17?? ? √17?? + 9?? + 8?? ?
?─────────────────────────────────────────────────────────────────────────────────────────────────────── ──────────────────────────────────────────────────────────────── ?
? 17?√17 + 153 √17 + 17 ?
添加回答
舉報
0/150
提交
取消