get_matrix() 把整条线路读成酉矩阵 $U$。读法：第 $j$ 列就是 $U|j\rangle$，即输入基态 $|j\rangle$ 演化出的态矢；列编号采用 little-endian 约定 $j = q_0 + 2q_1 + \ldots$，q[0] 是态矢索引的最低位。Bell 线路的 $4\times 4$ 矩阵非零元均为 $\pm 0.7071$（$1/\sqrt{2}$ 的四位小数舍入）。第 0 列 $(0.7071, 0, 0, 0.7071)^{T} = (|00\rangle+|11\rangle)/\sqrt{2} = |\Phi^+\rangle$；第 1 列对应 $j = 1$（q0=1、q1=0），H(q0) 给出 $(|00\rangle-|10\rangle)/\sqrt{2}$，其中 q0=1 的分量触发 CNOT 翻转 q1，末态 $(|00\rangle-|11\rangle)/\sqrt{2} = |\Phi^-\rangle$，与矩阵第 1 列 $(0.7071, 0, 0, -0.7071)^{T}$ 逐元一致——该线路把整个计算基映成 Bell 基。两个限制：矩阵规模 $2^n\times 2^n$ 使元素数随 $n$ 每加 1 乘 4，$n = 12$ 时约 $1.7\times 10^7$ 个元素，故约 12 比特以内适用；含测量的线路没有酉矩阵表示（测量不可逆），强行提取抛 NotMatrixableError: Measured circuits have no unitary matrix。