1. 课程导入与仿真环境搭建

量子噪声信道与量子态保真度仿真，核心目标只有一句话：把“量子态在噪声中如何退化”这件事，从抽象公式变成可计算、可观察、可比较的数字与曲线。在理想情况下，一个量子态经过信道后保持不变；但在真实物理系统中，退相干、能量弛豫、控制误差都会让输出态偏离输入态。保真度正是衡量这种偏离程度的标尺——它取值在 0 到 1 之间，越接近 1 表示输出态与目标态越接近。本课程将围绕“构造量子态 → 施加噪声信道 → 计算输出态 → 用保真度量化偏差 → 分析参数影响”这条主线展开，最终使你能够独立搭建噪声模型、批量扫描参数并解读仿真结果。

整体学习路径可以概括为四个层次。第一层是数学与物理基础：密度矩阵、迹、部分迹、Kraus 算符与量子信道映射，这些是后续所有代码背后的语言。第二层是典型噪声信道建模：比特翻转、相位翻转、去极化、振幅阻尼、相位阻尼，每一种信道都有明确的物理来源和 Kraus 表示。第三层是保真度计算与对比：区分纯态保真度与一般态保真度公式，理解为什么不同噪声对同一初始态的影响差异巨大。第四层是参数扫描与结果分析：以噪声强度 p、时间 t、阻尼率 γ 等为自变量，绘制保真度衰减曲线，并从中读出信道的“破坏性”特征。这四层不是割裂的，而是层层递进：没有正确的信道算符，保真度计算就失去意义；没有保真度对比，参数扫描就只是数字堆砌。

工欲善其事，必先利其器。本课程使用 Python 生态中的 QuTiP（Quantum Toolbox in Python）作为主仿真框架，它内置了量子态、算符、信道、保真度等对象，能显著减少重复造轮子的成本。推荐使用 Python 3.10 及以上版本，通过虚拟环境隔离依赖，避免与系统其他包冲突。安装命令如下：

python -m venv qnoise-env
# Windows 激活
qnoise-env\Scripts\activate
# macOS / Linux 激活
source qnoise-env/bin/activate

pip install --upgrade pip
pip install qutip numpy matplotlib scipy

安装完成后，必须做一次最小可用性验证，而不是装完就假定成功。下面这段代码同时检查版本、构造一个单量子比特态、施加一个去极化信道并计算保真度，是后续所有实验的“冒烟测试”：

import qutip
import numpy as np

print("QuTiP 版本:", qutip.__version__)

# 构造 |0> 态
psi0 = qutip.basis(2, 0)
rho0 = psi0 * psi0.dag()

# 去极化信道，噪声强度 p = 0.2
p = 0.2
rho_out = (1 - p) * rho0 + p * qutip.qeye(2) / 2

# 保真度：纯态与一般态
F = qutip.fidelity(rho0, rho_out)
print("输出态:\n", rho_out)
print("保真度:", F)

若环境正常，你会看到 QuTiP 版本号、一个 2×2 密度矩阵，以及约 0.9 的保真度。这里有一个值得注意的细节：去极化信道把 |0⟩ 态以概率 p 替换为完全混合态 I/2，因此保真度随 p 线性下降，p=0.2 时理论值为 0.9。这个简单例子已经包含了本课程的全部要素——初始态、信道映射、输出态、保真度指标。后续章节只是把信道换得更“真实”、把初始态换得更“复杂”、把单点计算扩展为参数扫描。

为了让仿真结果可复现，建议在脚本开头固定随机种子，并统一约定：密度矩阵用 qutip.Qobj 表示，保真度统一调用 qutip.fidelity，噪声强度参数统一记为 p（0 到 1）或 γ（非负）。此外，绘图时横轴统一为噪声参数、纵轴统一为保真度，便于跨信道横向比较。下面这张图展示了本节的核心工作流与各环节的输入输出关系，可作为你搭建自己实验脚本时的结构参照。

这张图里有两个容易被忽视但很关键的点。其一是第③步的 Kraus 表示：任何完全正的保迹映射都可以写成 ρ_out = Σᵢ Kᵢ ρ₀ Kᵢ†，且满足 Σᵢ Kᵢ†Kᵢ = I。理解这一点，你就能自己推导去极化、振幅阻尼等信道的算符，而不必死记公式。其二是第⑥步的判读逻辑：保真度曲线下降得越快，说明该信道在相同噪声强度下对量子态的破坏越强；若曲线趋于某个非零平台，则说明信道存在稳态，这在振幅阻尼信道中尤为典型——长时间演化后系统趋于 |0⟩，对 |0⟩ 初始态的保真度反而回升。

最后给出一个可直接执行的配置检查清单，建议在开始每章实验前逐项确认：Python 版本不低于 3.10；qutip.__version__ 能正常打印；qutip.fidelity 对相同态返回 1.0；密度矩阵的迹约等于 1（可用 rho.tr() 检查）；绘图后端可正常弹出窗口或保存 PNG。这五项全部通过，说明你的仿真环境已经就绪，可以进入后续的信道建模与保真度计算环节。