如何快速消除下式中的KroneckerDelta，目前我只会半自动求解，利用Block对每一项分别令KroneckerDelta下的角标相等进行化简，但是这样觉得很是愚蠢，但是查看了Documentation之后并没有找到快速消除Kroneckerdelta的方法。
希望得到指点。
待求式子如下：
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -q2] KroneckerDelta[
p2, -q1] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] +
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -q1] KroneckerDelta[
p2, -q2] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -p2] KroneckerDelta[
q1, -q2] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q2] KroneckerDelta[p1, -q1] KroneckerDelta[-p2,
k - p1 - p2 - q1 - q2] Sinh[Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q1] KroneckerDelta[p1, -q2] KroneckerDelta[-p2,
k - p1 - p2 - q1 - q2] Sinh[Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] +
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q2] KroneckerDelta[-p1,
k - p1 - p2 - q1 - q2] KroneckerDelta[p2, -q1] Sinh[
Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2]
希望得到指点。
待求式子如下：
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -q2] KroneckerDelta[
p2, -q1] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] +
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -q1] KroneckerDelta[
p2, -q2] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] +
Subscript[\[Lambda], q1] - Subscript[\[Lambda], q2])
Cosh[Subscript[\[Lambda], k]] KroneckerDelta[k,
k - p1 - p2 - q1 - q2] KroneckerDelta[p1, -p2] KroneckerDelta[
q1, -q2] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q2] KroneckerDelta[p1, -q1] KroneckerDelta[-p2,
k - p1 - p2 - q1 - q2] Sinh[Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] -
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q1] KroneckerDelta[p1, -q2] KroneckerDelta[-p2,
k - p1 - p2 - q1 - q2] Sinh[Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2] +
1/2 E^(-Subscript[\[Lambda], p1] - Subscript[\[Lambda], p2] -
Subscript[\[Lambda], q1] +
2 (Subscript[\[Lambda], p1] + Subscript[\[Lambda], q1]) -
Subscript[\[Lambda], q2])
KroneckerDelta[-k, -q2] KroneckerDelta[-p1,
k - p1 - p2 - q1 - q2] KroneckerDelta[p2, -q1] Sinh[
Subscript[\[Lambda], k]] Sinh[
Subscript[\[Lambda], -k + p1 + p2 + q1 + q2]] Subscript[h, p1,
p1 + p2, p2] Subscript[h, q1, q1 + q2, q2]
