读出错误给每个比特配一个 $2\times 2$ 混淆矩阵（行是真值、列是读数），readout_error = {q: [p(0→1), p(1→0)]}，输出分布等于各比特边缘概率乘上混淆矩阵。设 GHZ 理想分布 $P(000) = P(111) = 0.5$，只给 q0（位串最右位）配 {0: [0.10, 0.02]}，四个非零概率可全部手算：真值 0 分支 $P(000) = 0.5\times(1-0.10) = 0.4500$、$P(001) = 0.5\times 0.10 = 0.0500$（0 被读成 1，翻转的是最右位）；真值 1 分支 $P(111) = 0.5\times(1-0.02) = 0.4900$、$P(110) = 0.5\times 0.02 = 0.0100$。验证 $0.45+0.05+0.01+0.49 = 1.00$——混淆矩阵行和为 1（$0.90+0.10$ 与 $0.02+0.98$），总概率守恒。NoisySimulator(backend_type="density_matrix", readout_error={0: [0.10, 0.02]}) 的运行输出与手算逐位一致；限制是密度矩阵路径不支持 simulate_shots（抛 NotImplementedError: Density matrix is not supported for shot simulation.），端到端含噪采样改用 dummy 后端。