从 $|00\rangle$ 出发，两步即可制备贝尔态。第一步用"分别作用"规则：$(H\otimes I)|00\rangle = (H|0\rangle)\otimes|0\rangle = (|00\rangle+|10\rangle)/\sqrt{2}$，中间态仍是乘积态（即 $|+\rangle\otimes|0\rangle$），叠加有了、纠缠还没有。第二步 CNOT 按线性分配到每个叠加分量：$\mathrm{CNOT}\,(|00\rangle+|10\rangle)/\sqrt{2} = (\mathrm{CNOT}|00\rangle+\mathrm{CNOT}|10\rangle)/\sqrt{2} = (|00\rangle+|11\rangle)/\sqrt{2} = |\Phi^+\rangle$。中间态可分解、末态不可分解（第 9 页命题兑现），纠缠恰在第二步产生。单比特门作用在多比特系统上还有两种等价算法：张量积法把 $H\otimes I$ 展成 $4\times 4$ 大矩阵作用于 $(\alpha,\beta,\gamma,\delta)^{T}$，得 $((\alpha+\gamma),(\beta+\delta),(\alpha-\gamma),(\beta-\delta))^{T}/\sqrt{2}$，适合编程模拟；线性展开法让门只作用于目标比特的基矢、逐项展开再合并同类项，手算推荐。