2. 量子态描述与密度矩阵基础
在量子态层析中，我们最终要回答的问题是：给定一组测量数据，如何反推出量子系统真实的状态。要回答这个问题，首先必须有一套能同时描述“纯态”和“混态”的数学语言，这就是密度矩阵。本节从纯态与混态的区别出发，建立密度矩阵的定义与性质，再借助布洛赫球获得几何直觉，最后落到量子测量公设与期望值计算，为后续层析算法打下基础。
纯态与混态：两种不同的“不确定性”
纯态是可以用一个态矢量 $|\psi\rangle$ 完全描述的量子态。例如单量子比特的 $|0\rangle$、$|1\rangle$、$(|0\rangle+|1\rangle)/\sqrt{2}$ 都是纯态。对纯态做同样的测量，结果的不确定性来自量子力学的相干叠加，是“内禀”的。
混态则不同。它表示系统以概率 $p_i$ 处于某个纯态 $|\psi_i\rangle$，而我们并不知道具体是哪一支。这种不确定性来自经典意义上的“信息缺失”。例如一个量子比特以 50% 概率处于 $|0\rangle$、50% 概率处于 $|1\rangle$，它不是叠加态 $(|0\rangle+|1\rangle)/\sqrt{2}$，而是一个混态。两者的测量统计在特定基下可能相似，但在其它基下会明显不同，这正是层析要区分的关键。
密度矩阵的定义与性质
为了统一描述纯态和混态，引入密度算符：
ρ = Σ_i p_i |ψ_i⟩⟨ψ_i|, Σ_i p_i = 1, p_i ≥ 0纯态是它的特例，此时只有一项 $p=1$。密度矩阵具有三条核心性质：
- 厄米性：$\rho = \rho^\dagger$，保证测量期望值为实数。
- 单位迹：$\mathrm{Tr}(\rho)=1$，对应概率归一。
- 半正定：$\rho \ge 0$，所有本征值非负，保证任意测量结果概率非负。
由此可定义纯度 $\mathrm{Tr}(\rho^2)$。纯态满足 $\mathrm{Tr}(\rho^2)=1$，混态则 $\mathrm{Tr}(\rho^2)<1$。对单量子比特，密度矩阵可写成
ρ = (I + r_x σ_x + r_y σ_y + r_z σ_z) / 2其中 $\vec r=(r_x,r_y,r_z)$ 是布洛赫矢量，$\sigma_{x,y,z}$ 为泡利矩阵。纯态对应 $|\vec r|=1$，混态对应 $|\vec r|<1$，完全混态对应 $\vec r=0$。
量子测量公设与期望值计算
量子测量由一组测量算符 $\{M_m\}$ 描述，指标 $m$ 对应可能的测量结果，且满足完备性关系 $\sum_m M_m^\dagger M_m = I$。若系统处于密度矩阵 $\rho$，得到结果 $m$ 的概率为
P(m) = Tr(M_m† M_m ρ)测量后系统状态坍缩为
ρ_m = (M_m ρ M_m†) / P(m)对于投影测量，$M_m$ 是投影算符 $P_m$，此时 $P(m)=\mathrm{Tr}(P_m\rho)$。可观测量 $A$ 的期望值为
⟨A⟩ = Tr(ρ A)以单量子比特为例，若 $\rho=(I+\vec r\cdot\vec\sigma)/2$，测量 $\sigma_z$ 时，$P(+1)=\mathrm{Tr}(P_0\rho)=(1+r_z)/2$，$P(-1)=(1-r_z)/2$。因此只要统计 $|0\rangle$ 和 $|1\rangle$ 的出现频率，就能估计 $r_z$。类似地，测量 $\sigma_x$、$\sigma_y$ 可分别得到 $r_x$、$r_y$。这正是量子态层析的基本思路：通过在不同基下测量，把密度矩阵的未知参数逐个“投影”出来。
从测量频率到密度矩阵：一个可操作的示例
假设我们在 Z 基下对 1000 个单量子比特副本测量，得到 700 次 $|0\rangle$、300 次 $|1\rangle$。则
r_z ≈ (700 - 300) / 1000 = 0.4再在 X 基下测量，得到 600 次 $|+\rangle$、400 次 $|-\rangle$，则
r_x ≈ (600 - 400) / 1000 = 0.2在 Y 基下测量，得到 550 次 $|+i\rangle$、450 次 $|-i\rangle$，则
r_y ≈ (550 - 450) / 1000 = 0.1于是 $\vec r=(0.2,0.1,0.4)$，$|\vec r|\approx 0.458<1$，说明这是一个混态。代入密度矩阵表达式可得
ρ ≈ 1/2 [ [1+0.4, 0.2-0.1i],
[0.2+0.1i, 1-0.4] ]
= [ [0.7, 0.1-0.05i],
[0.1+0.05i, 0.3] ]可以验证 $\mathrm{Tr}(\rho)=1$，$\rho=\rho^\dagger$，且本征值约为 0.771 和 0.229，均非负，满足密度矩阵要求。纯度 $\mathrm{Tr}(\rho^2)\approx 0.771^2+0.229^2\approx 0.647<1$，与混态判断一致。