答案就在qz算法的原始论文中，你可以参考https://epubs.siam.org/doi/abs/10.1137/0710024?journalCode=sjnaam第251页的第5步，Fortran代码也可以在这篇论文的旧版本中找到。

我遇到了同样的问题，最终转向了scipy提供的低级LAPACK函数。
在matlab中，这些量可以被找到为

[AA,BB,Q,Z,V,W] = qz(A,B)

在python中，同样的结果可以通过

from scipy.linalg import qz
from scipy.linalg.lapack import zggev

# matlab uses complex output by default

# python uses real

AA, BB, Q, Z = qz(A, B, 'complex')
_, _, V, W, _, _ = zggev(A, B)

不同版本的LAPACK ggev函数取决于数据类型（真实的与复数，单精度与双精度）。文档here以及其他许多地方。