【超详解】量子克拉美罗下界QCRB：时延估计精度极限推导+MATLAB多场景绘图实战

文章目录

- [第一章 QCRB基础认知与应用场景](#第一章 QCRB基础认知与应用场景)
- - [1.1 从经典估计理论到量子精度极限](#1.1 从经典估计理论到量子精度极限)
  - [1.2 时延估计QCRB公式的物理含义](#1.2 时延估计QCRB公式的物理含义)
- [第二章时延估计QCRB的完整推导逻辑](#第二章时延估计QCRB的完整推导逻辑)
- - [2.1 量子Fisher信息与下界的对应关系](#2.1 量子Fisher信息与下界的对应关系)
  - [2.2 高斯脉冲下的QCRB公式推导过程](#2.2 高斯脉冲下的QCRB公式推导过程)
- [第三章 MATLAB基础仿真：单参数影响曲线绘制](#第三章 MATLAB基础仿真：单参数影响曲线绘制)
- - [3.1 脉宽τ_p对时延估计下界的影响](#3.1 脉宽τ_p对时延估计下界的影响)
  - [3.2 脉冲数量n对精度极限的缩放规律](#3.2 脉冲数量n对精度极限的缩放规律)
- [第四章多场景QCRB对比仿真与结果分析](#第四章多场景QCRB对比仿真与结果分析)
- - [4.1 不同量子态的精度极限对比](#4.1 不同量子态的精度极限对比)
  - [4.2 信道损耗对QCRB的劣化效应分析](#4.2 信道损耗对QCRB的劣化效应分析)
- [第五章工程应用中的参考与常见误区](#第五章工程应用中的参考与常见误区)
- - [5.1 脉宽定义与公式适配注意事项](#5.1 脉宽定义与公式适配注意事项)
  - [5.2 QCRB与实际实验精度的差距](#5.2 QCRB与实际实验精度的差距)
  - [5.3 测量次数与光子数的区分](#5.3 测量次数与光子数的区分)

第一章 QCRB基础认知与应用场景

1.1 从经典估计理论到量子精度极限

经典参数估计理论中，克拉美罗下界（CRB）给出了无偏估计量所能达到的最小方差，是统计推断中的精度基准。当估计场景进入量子尺度，待测参数编码在量子态上，经典CRB不再适用，此时需要用量子克拉美罗下界（QCRB）描述物理定律允许的终极估计精度。

QCRB的核心价值在于，它不依赖具体测量方案，仅由量子态本身的性质决定，是所有量子测量方案都无法突破的精度下限。在光学位移测量、时延估计、相位传感、磁力计等精密测量领域，QCRB是评估方案性能、判断是否实现量子增强的核心标尺。

与时域测量直接相关的场景中，最典型的就是光脉冲时延估计------无论是光梳测距、光纤传感还是量子雷达，最终都可以归结为对脉冲到达时间的估计，而QCRB直接决定了这套系统理论上能达到的最高时间分辨率。

1.2 时延估计QCRB公式的物理含义

对于高斯型光脉冲的时延估计问题，QCRB的最简形式可以写成：

σ t = γ τ p n (1) \sigma_{t}=\gamma \frac{\tau_{p}}{\sqrt{n}} \tag{1} σt=γn τp(1)

公式中每个参数的物理意义明确：

σ t \sigma_t σt：时延估计量的标准差，代表估计精度极限，数值越小精度越高；
τ p \tau_p τp：光脉冲的时域半高全宽，也就是常说的脉宽；
n n n：独立重复测量的脉冲数量，也就是累计测量次数；
γ \gamma γ：无量纲系数，由输入量子态的类型、测量方案的优劣决定，是区分经典测量与量子增强测量的核心参数。

从公式可以直接读出两个基础规律：

第一，估计精度与脉宽成正比，脉冲越窄，理论上能达到的时延精度越高，这也是超短脉冲激光在精密测量中广泛应用的原因；

第二，估计精度与测量次数的平方根成反比，这是统计估计的通用缩放规律，量子增强也无法突破这一统计极限，只能通过降低 γ \gamma γ系数来整体抬升精度基线。

经典相干态探测的极限下， γ \gamma γ存在一个固定的经典极限值；而利用压缩态、纠缠态等非经典量子态，可以让 γ \gamma γ小于经典极限，实现所谓的"超分辨"量子增强测量。

第二章时延估计QCRB的完整推导逻辑

2.1 量子Fisher信息与下界的对应关系

QCRB的推导起点是量子Fisher信息 （Quantum Fisher Information, QFI）。对于一个携带待估计参数 θ \theta θ的量子态 ρ ( θ ) \rho(\theta) ρ(θ)，量子Fisher信息 F Q F_Q FQ定量描述了该量子态中包含的关于参数 θ \theta θ的信息量。

量子克拉美罗下界的通用表达形式为：

Var ( θ ^ ) ≥ 1 n F Q (2) \text{Var}(\hat{\theta}) \geq \frac{1}{n F_Q} \tag{2} Var(θ^)≥nFQ1(2)

其中 Var ( θ ^ ) \text{Var}(\hat{\theta}) Var(θ^)是无偏估计量 θ ^ \hat{\theta} θ^的方差， n n n为独立测量次数。该式表明，估计方差的下限由量子Fisher信息和测量次数共同决定，信息越丰富、测量次数越多，理论精度越高。

对于纯态 ∣ ψ ( θ ) ⟩ |\psi(\theta)\rangle ∣ψ(θ)⟩的情况，量子Fisher信息可以简化为：

F Q = 4 ( ⟨ ∂ θ ψ ∣ ∂ θ ψ ⟩ − ∣ ⟨ ∂ θ ψ ∣ ψ ⟩ ∣ 2 ) (3) F_Q = 4\left( \left\langle \partial_\theta \psi | \partial_\theta \psi \right\rangle - \left| \left\langle \partial_\theta \psi | \psi \right\rangle \right|^2 \right) \tag{3} FQ=4(⟨∂θψ∣∂θψ⟩−∣⟨∂θψ∣ψ⟩∣2)(3)

这个式子的物理意义是：参数变化引起的量子态在希尔伯特空间中的变化速率，速率越快，态对参数越敏感，估计精度就越高。

2.2 高斯脉冲下的QCRB公式推导过程

时延估计中，待估计参数是时间偏移 t t t，脉冲的时域波函数会整体平移 t t t。对于变换极限高斯脉冲，时域电场包络为高斯分布，对应的频谱同样是高斯分布。

我们将时延 t t t作为待估参数，脉冲的频域波函数可以表示为 E ( ω ) e − i ω t E(\omega) e^{-i\omega t} E(ω)e−iωt，代入纯态量子Fisher信息公式计算后，可以得到频域形式的QFI：

F Q = 4 ∫ − ∞ + ∞ ( ω − ⟨ ω ⟩ ) 2 ∣ E ( ω ) ∣ 2 d ω = 4 ⋅ Var ( ω ) (4) F_Q = 4 \int_{-\infty}^{+\infty} (\omega - \langle \omega \rangle)^2 |E(\omega)|^2 d\omega = 4 \cdot \text{Var}(\omega) \tag{4} FQ=4∫−∞+∞(ω−⟨ω⟩)2∣E(ω)∣2dω=4⋅Var(ω)(4)

也就是说，时延估计的量子Fisher信息，等于脉冲频谱角频率方差的4倍。频谱越宽，频率方差越大，对时延的敏感度就越高。

对于半高全宽为 τ p \tau_p τp的高斯时域脉冲，其角频率的半高全宽 Ω p \Omega_p Ωp满足傅里叶变换极限关系：

τ p ⋅ Ω p = 4 ln ⁡ 2 (5) \tau_p \cdot \Omega_p = 4\ln2 \tag{5} τp⋅Ωp=4ln2(5)

高斯分布的方差与半高全宽的关系为 σ ω = Ω p / ( 2 2 ln ⁡ 2 ) \sigma_\omega = \Omega_p / (2\sqrt{2\ln2}) σω=Ωp/(22ln2 )，代入方差计算后可得：

Var ( ω ) = σ ω 2 = 2 ln ⁡ 2 τ p 2 (6) \text{Var}(\omega) = \sigma_\omega^2 = \frac{2\ln2}{\tau_p^2} \tag{6} Var(ω)=σω2=τp22ln2(6)

将其代入量子Fisher信息表达式，再代入QCRB通用公式，最终整理后就可以得到式(1)的形式。其中系数 γ \gamma γ包含了量子态性质、测量效率等修正项：理想相干态平衡零拍测量下对应固定的经典极限值；引入压缩态等非经典资源后，等效频谱方差会增大， γ \gamma γ随之减小，精度极限进一步降低。

第三章 MATLAB基础仿真：单参数影响曲线绘制

3.1 脉宽τ_p对时延估计下界的影响

首先分析固定测量次数下，脉冲宽度与时延精度极限的定量关系。我们分别仿真经典相干态和理想压缩态两种场景，绘制 σ t \sigma_t σt随 τ p \tau_p τp的变化曲线。

以下是完整MATLAB代码，可直接运行生成高清曲线：

matlab 复制代码

%% 图1：脉宽对时延QCRB的影响曲线
clear; clc; close all;

% 参数设置
n_pulse = 100;          % 测量脉冲数
tau_p_list = linspace(10, 1000, 200);  % 脉宽范围 10fs ~ 1000fs

% 不同量子态的gamma系数
gamma_coherent = 0.28;  % 相干态经典极限
gamma_squeezed = 0.12;  % 10dB压缩态理想值

% 计算QCRB
sigma_coherent = gamma_coherent * tau_p_list / sqrt(n_pulse);
sigma_squeezed = gamma_squeezed * tau_p_list / sqrt(n_pulse);

% 绘图
figure('Color','w','Position',[100 100 600 400]);
plot(tau_p_list, sigma_coherent, 'b-', 'LineWidth', 2);
hold on;
plot(tau_p_list, sigma_squeezed, 'r--', 'LineWidth', 2);

% 格式设置
xlabel('脉冲半高全宽 \tau_p (fs)','FontSize',12);
ylabel('时延估计精度极限 \sigma_t (fs)','FontSize',12);
legend('相干态经典极限','10dB压缩态理想QCRB','Location','northwest');
grid on;
set(gca,'FontSize',11);
title('固定测量次数下脉宽与时延精度极限的关系','FontSize',13);

运行代码后可以得到线性增长的曲线：脉宽越窄，精度极限越高；相同脉宽下，压缩态的QCRB显著低于相干态，体现了量子增强的优势。需要注意的是，这里的 γ \gamma γ取值对应理想无损耗场景，实际实验中会因损耗、探测效率等因素有所上升。

3.2 脉冲数量n对精度极限的缩放规律

第二个基础仿真是验证 1 / n 1/\sqrt{n} 1/n 的统计缩放规律。固定脉冲宽度，观察累计测量脉冲数从1到10000变化时，精度极限的变化趋势，采用双对数坐标可以更直观地看出幂律关系。

对应MATLAB代码如下：

matlab 复制代码

%% 图2：测量脉冲数对QCRB的影响（双对数坐标）
clear; clc; close all;

% 参数设置
tau_p = 100;    % 固定脉宽100fs
n_list = logspace(0, 4, 200);  % 脉冲数 1 ~ 10000

gamma_coherent = 0.28;
gamma_squeezed = 0.12;

sigma_coherent = gamma_coherent * tau_p ./ sqrt(n_list);
sigma_squeezed = gamma_squeezed * tau_p ./ sqrt(n_list);

% 绘图
figure('Color','w','Position',[100 100 600 400]);
loglog(n_list, sigma_coherent, 'b-', 'LineWidth', 2);
hold on;
loglog(n_list, sigma_squeezed, 'r--', 'LineWidth', 2);

xlabel('独立测量脉冲数 n','FontSize',12);
ylabel('时延估计精度极限 \sigma_t (fs)','FontSize',12);
legend('相干态经典极限','10dB压缩态理想QCRB','Location','northeast');
grid on;
set(gca,'FontSize',11);
title('固定脉宽下测量次数与精度极限的缩放关系','FontSize',13);

双对数图中曲线呈现斜率为-0.5的直线，完美对应 1 / n 1/\sqrt{n} 1/n 的幂律关系。这也说明，无论是经典方案还是量子增强方案，都遵循统计测量的基本规律；量子增强的作用是将整条曲线向下平移，在相同测量次数下获得更高精度。

第四章多场景QCRB对比仿真与结果分析

4.1 不同量子态的精度极限对比

不同的输入量子态对应不同的 γ \gamma γ系数，也就对应不同的精度极限。我们选取四种典型场景进行对比：热光态、相干态、10dB压缩态、理想NOON态，以柱状图的形式直观展示精度差异。

MATLAB绘图代码如下：

matlab 复制代码

%% 图3：不同量子态下的时延QCRB对比
clear; clc; close all;

% 参数设置
tau_p = 100;    % 脉宽100fs
n_pulse = 100;  % 测量100个脉冲

% 各量子态的gamma系数
state_names = {'热光态','相干态','10dB压缩态','理想NOON态'};
gamma_list = [0.39, 0.28, 0.12, 0.07];

sigma_list = gamma_list * tau_p / sqrt(n_pulse);

% 绘图
figure('Color','w','Position',[100 100 600 420]);
bar_colors = [0.6 0.6 0.6; 0.2 0.4 0.8; 0.8 0.2 0.2; 0.5 0 0.5];
b = bar(sigma_list, 0.6);
b.FaceColor = 'flat';
for i = 1:length(gamma_list)
    b.CData(i,:) = bar_colors(i,:);
end

% 数值标注
for i = 1:length(sigma_list)
    text(i, sigma_list(i)+0.2, sprintf('%.2f fs',sigma_list(i)),...
        'HorizontalAlignment','center','FontSize',10);
end

set(gca,'XTickLabel',state_names,'FontSize',11);
ylabel('时延估计精度极限 \sigma_t (fs)','FontSize',12);
title('100fs脉宽、100次测量下不同量子态的QCRB对比','FontSize',13);
ylim([0 max(sigma_list)*1.2]);
grid on;

从结果可以看出，从经典热光到相干态已经有明显的精度提升；而压缩态和NOON态等非经典资源，可以进一步将精度极限推低数倍。需要注意的是，NOON态等多光子纠缠态虽然理论 γ \gamma γ极低，但制备难度大、抗损耗能力差，实际工程中很少直接使用。

4.2 信道损耗对QCRB的劣化效应分析

实际应用场景中，光脉冲在传输路径上不可避免会遇到损耗，比如光纤传输、自由空间衰减、探测器效率不足等。损耗会让量子态退化为混合态，导致量子Fisher信息下降， γ \gamma γ系数上升，最终劣化估计精度。

对于衰减信道，透过率为 η \eta η（0<η≤1）时，压缩态的等效 γ \gamma γ系数满足如下近似关系：

γ ( η ) = γ 0 η \gamma(\eta) = \frac{\gamma_0}{\sqrt{\eta}} γ(η)=η γ0

其中 γ 0 \gamma_0 γ0为无损耗时的理想系数。我们绘制损耗从0到20dB变化时，相干态与压缩态的QCRB变化曲线，对比两者的抗损耗能力。

对应MATLAB代码：

matlab 复制代码

%% 图4：信道损耗对QCRB的影响
clear; clc; close all;

% 参数设置
tau_p = 100;
n_pulse = 100;
loss_dB = linspace(0, 20, 200);  % 损耗 0~20dB
eta = 10.^(-loss_dB/10);         % 透过率

gamma0_coh = 0.28;   % 相干态无损耗gamma
gamma0_sq = 0.12;    % 压缩态无损耗gamma

% 损耗下的gamma与sigma
gamma_coh = gamma0_coh ./ sqrt(eta);
gamma_sq = gamma0_sq ./ sqrt(eta);

sigma_coh = gamma_coh * tau_p / sqrt(n_pulse);
sigma_sq = gamma_sq * tau_p / sqrt(n_pulse);

% 绘图
figure('Color','w','Position',[100 100 600 400]);
plot(loss_dB, sigma_coh, 'b-', 'LineWidth', 2);
hold on;
plot(loss_dB, sigma_sq, 'r--', 'LineWidth', 2);

xlabel('信道损耗 (dB)','FontSize',12);
ylabel('时延估计精度极限 \sigma_t (fs)','FontSize',12);
legend('相干态','10dB压缩态','Location','northwest');
grid on;
set(gca,'FontSize',11);
title('信道损耗对时延估计QCRB的劣化效应','FontSize',13);

仿真结果显示，损耗会让两种方案的精度都同步劣化；损耗越大，压缩态相对于相干态的绝对优势会被削弱。在实际长距离传感、测距场景中，必须将损耗纳入考量，否则会严重高估量子增强的实际收益。

第五章工程应用中的参考与常见误区

5.1 脉宽定义与公式适配注意事项

QCRB公式中的 τ p \tau_p τp默认为时域半高全宽，但不同文献、不同仿真中可能采用均方根脉宽、1/e脉宽等不同定义，这是套用公式时最容易出错的地方。

三种常见脉宽的换算关系明确：

半高全宽 τ FWHM \tau_{\text{FWHM}} τFWHM与均方根脉宽 τ rms \tau_{\text{rms}} τrms： τ FWHM = 2 2 ln ⁡ 2 ⋅ τ rms ≈ 2.355 τ rms \tau_{\text{FWHM}} = 2\sqrt{2\ln2} \cdot \tau_{\text{rms}} \approx 2.355 \tau_{\text{rms}} τFWHM=22ln2 ⋅τrms≈2.355τrms
半高全宽与1/e脉宽 τ 1 / e \tau_{1/e} τ1/e： τ FWHM = 2 ln ⁡ 2 ⋅ τ 1 / e ≈ 1.177 τ 1 / e \tau_{\text{FWHM}} = \sqrt{2\ln2} \cdot \tau_{1/e} \approx 1.177 \tau_{1/e} τFWHM=2ln2 ⋅τ1/e≈1.177τ1/e

如果拿到的脉宽数据是均方根定义，直接代入式(1)计算，得到的精度极限会比真实值偏小2倍以上，造成对系统性能的误判。工程计算前必须先统一脉宽定义，再代入公式。

5.2 QCRB与实际实验精度的差距

QCRB是理想条件下的理论极限，不代表实际实验中一定能达到。实际测量精度劣于QCRB是常态，主要来源包括：测量方案不是最优方案、存在额外经典噪声、探测器存在暗噪声、模式匹配不完美、脉冲存在啁啾偏离变换极限等。

评估一套量子测量方案的性能时，正确的做法是将实际测得的精度与QCRB对比，计算两者的比值。比值越接近1，说明测量方案越接近最优；如果实验精度优于QCRB，要么是参数定义出错，要么是实验结果存在问题。

5.3 测量次数与光子数的区分

很多资料中会将式(1)里的 n n n解释为光子数，这是典型的概念混淆。公式中的 n n n是独立重复测量的次数，也就是发送了多少个独立脉冲；而单个脉冲内的平均光子数，是通过 γ \gamma γ系数来影响精度的。

光子数增加会提升信噪比，等效于降低 γ \gamma γ值，但不会改变 1 / n 1/\sqrt{n} 1/n 的缩放规律。在光子数受限的场景（比如弱光测距、量子雷达），需要同时考虑单脉冲光子数对 γ \gamma γ的影响，以及累计脉冲数对统计精度的提升，不能直接用总光子数代替 n n n代入公式。

你在量子精密测量仿真或者时延估计系统设计中，遇到过哪些参数定义混淆、理论与实测不符的问题？欢迎留言交流。

【超详解】量子克拉美罗下界QCRB：时延估计精度极限推导+MATLAB多场景绘图实战

文章目录

第一章 QCRB基础认知与应用场景

1.1 从经典估计理论到量子精度极限

1.2 时延估计QCRB公式的物理含义

第二章 时延估计QCRB的完整推导逻辑

2.1 量子Fisher信息与下界的对应关系

2.2 高斯脉冲下的QCRB公式推导过程

第三章 MATLAB基础仿真：单参数影响曲线绘制

3.1 脉宽τ_p对时延估计下界的影响

3.2 脉冲数量n对精度极限的缩放规律

第四章 多场景QCRB对比仿真与结果分析

4.1 不同量子态的精度极限对比

4.2 信道损耗对QCRB的劣化效应分析

第五章 工程应用中的参考与常见误区

5.1 脉宽定义与公式适配注意事项

5.2 QCRB与实际实验精度的差距

5.3 测量次数与光子数的区分

第二章时延估计QCRB的完整推导逻辑

第四章多场景QCRB对比仿真与结果分析

第五章工程应用中的参考与常见误区