SymPy无法简化矩阵与其共轭矩阵的克罗内克积问题求助

你的问题出在SymPy默认的simplify函数不会主动将复指数项转换为三角形式，也没有充分利用你定义的符号实值假设来化简根式乘积。以下是针对性的解决步骤：

方法1：针对元素逐个化简，结合三角展开 #

先利用符号的实值属性简化共轭矩阵，再对克罗内克积的每个元素进行三角展开和根式化简：

import sympy as sp

# 定义符号，保留原假设
e = sp.Symbol('epsilon', real=True, domain=sp.Interval(0, 1, right_open=True))
i = sp.I
p = sp.Symbol('psi', real=True)

# 定义琼斯矩阵
J = sp.Matrix([
    [sp.sqrt(1 - e), -sp.sqrt(e) * sp.exp(i*p)],
    [sp.sqrt(e) * sp.exp(-i*p), sp.sqrt(1 - e)]
])

# 简化共轭矩阵：利用实值符号，共轭后直接化简元素
J_conj = J.applyfunc(lambda elem: sp.conjugate(elem).simplify())

# 计算克罗内克积
J_prod = sp.kronecker_product(J, J_conj)

# 对每个元素依次进行：展开复指数、化简根式乘积、转换为三角形式
J_prod_simplified = J_prod.applyfunc(
    lambda elem: sp.trigsimp(sp.powsimp(sp.expand(elem, complex=True)))
)

print(J_prod_simplified)

方法2：先将复指数转为三角形式再计算 #

如果一开始就把矩阵里的复指数展开为余弦和正弦形式，后续化简会更直接：

import sympy as sp

e = sp.Symbol('epsilon', real=True, domain=sp.Interval(0, 1, right_open=True))
i = sp.I
p = sp.Symbol('psi', real=True)

# 用三角形式定义琼斯矩阵
J = sp.Matrix([
    [sp.sqrt(1 - e), -sp.sqrt(e) * (sp.cos(p) + i*sp.sin(p))],
    [sp.sqrt(e) * (sp.cos(p) - i*sp.sin(p)), sp.sqrt(1 - e)]
])

# 计算共轭矩阵（实部不变，虚部取反）
J_conj = J.conjugate()

# 计算克罗内克积后直接化简
J_prod_simplified = sp.kronecker_product(J, J_conj).applyfunc(sp.simplify)

print(J_prod_simplified)

关键说明 #

• sp.expand(elem, complex=True)：强制展开复指数项，把exp(i*a)*exp(i*b)合并为exp(i*(a+b))，为后续转三角形式做准备。
• sp.powsimp：利用epsilon的非负性，将sqrt(x)*sqrt(x)化简为x（比如sqrt(1-e)*sqrt(1-e)变为1-e）。
• sp.trigsimp：将复指数表达式转换为余弦/正弦形式，比如exp(2i*p)转为cos(2p) + i*sin(2p)，匹配你预期的结果形式。

内容的提问来源于stack exchange，提问作者Yên Bình