精读笔记

Problem Setting

论文标题：Optimal experimental design for passive imaging source problems（arXiv preprint / 2026-06-05）。

这篇论文实际处理的是 passive imaging source problem 的 sensor placement：观测不是一次 source 到 data 的线性响应，而是多个随机 realization 形成的 covariance/correlation data；未知量不是 source mean，而是 source covariance。作者主要落在 spatially uncorrelated source，即 Q=Diag(q) 的情形，因此反演对象是源强 q。

关键困难有两个层次。第一，数据维度平方化：m 个传感器产生 m×m covariance matrix，active inverse problem 中 m 维数据的 OED 代价在这里变成 m^2 维相关数据的代价。第二，设计变量 w 不再只是选择观测算子的若干行，而是在 F_w Diag(q) F_w^* 的左右两侧同时出现；进入 posterior precision 后近似表现为四次 w 依赖。这意味着 active OED 里常用的凸性、prior-preconditioned Hessian 低秩套路不能无缝复用。

以前方法卡住的地方不是“不会定义 A-optimality”，而是每次评估 J(w)=tr(C_post(w)) 都可能需要重新组装与设计相关的 correlation operator，并且伴随 PDE forward/adjoint solves。对于 Helmholtz 这类 PDE-governed map，优化循环中反复做这些操作基本不可接受。本文的关键矛盾是：OED 需要大量 design-dependent evaluations，但 passive imaging 的信息算子天然处于平方化的数据/算子空间；必须把设计相关部分和 PDE 相关部分拆开。

Motivation

已有 Bayesian OED for PDE inverse problems 的核心计算范式是：利用 posterior covariance update 的低秩结构，把 expensive PDE solves 放到 offline 或少量 Hessian action 中，然后 online 优化只做小矩阵运算。这套方法对 active imaging 很自然，因为数据维度随传感器线性增长，设计变量通常以 Diag(w)Γ^{-1}Diag(w) 的形式进入，A-optimal objective 在 relaxed weights 上常有较好凸性。

被动成像缺的正是一个同等级别的结构化 OED formulation。直接把 source covariance Q 作为未知会使参数维度也平方化；直接把 covariance data 向量化会使观测维度平方化；直接处理 Wishart likelihood 又会失去闭式 Gaussian posterior。作者的动机不是提出更精确统计模型，而是找到一个足够可计算、仍保留 passive imaging 本质的数据结构。

核心观察是：如果源空间非相关，则 unknown covariance 被压回 diagonal q；而 correlation forward map 的平方结构并非完全无结构，它可以通过 Schur product、Khatri-Rao product 和 Kronecker product 显式重排。换句话说，缺口在于没有把 passive covariance operator 的 product-space structure 用于 OED 压缩。

Core Idea

论文真正的核心是把被动成像 OED 的 posterior covariance update 从“巨大 correlation-data inverse problem”重组为“两层低秩 memory reuse”。第一层低秩近似压缩 PDE observation map 本身：M^{-1/2}F^TΓ^{-1/2}≈QR。这样任意设计 w 下的信息矩阵可以写成 Q L_w Q^T，L_w=R Diag(w)R^T。这里 PDE、mesh、sensor candidate geometry 的主要信息都被收进 Q,R，online 阶段只通过 L_w 注入设计。

第二层低秩近似处理 correlation data 带来的 product-space 扩张。由于 passive posterior precision 中出现 diag∘[(F^*Γ_w^†F)⊗(F^*Γ_w^†F)]∘Diag，第一层压缩后会自然出现 (Q^T * Q) 和 L_w⊗L_w。作者再对 prior-preconditioned product-space operator 做一次 QR/SVD 型压缩，得到 \hat Q,\hat R。最终 A-optimal objective 变成 tr(C0)−tr(\hat C)+tr((\hat L_w+I)^{-1}\hat C)，其中 \hat L_w=\hat R(L_w⊗L_w)\hat R^T。

本质区别在于它没有试图在 m^2 数据空间中直接优化，也没有把每个候选设计都映射回 PDE solve；它把信息流改成：PDE physics offline 编码，correlation-induced quadratic interactions offline 编码，design online 只改变低维 L_w。这个 inductive bias 很强：源必须 diagonal，信息必须能被低秩 product subspace 捕捉。但在这个假设下，它正好对准 passive imaging 的计算瓶颈。

Method

1. Bayesian passive imaging formulation：作者把 covariance estimator G 建模为 G=F Diag(q)F^*+E，并采用 Gaussian noise E∼N(0,Γ⊗Γ)。它解决的是获得闭式 posterior covariance 的问题；没有这个近似，就很难得到可优化的 A-optimal objective。核心变化是把 Wishart covariance estimation 问题转化为线性 Gaussian inverse problem。但这也是统计上最脆弱的一步。

2. Spatially uncorrelated source parameterization：把 full covariance Q 限制为 Diag(q)。它解决参数维度平方化问题，使 unknown 从 operator-valued covariance 回到 function-valued source strength。这个假设是方法可扩展的根本原因，不是普通建模便利。

3. diag/Schur product discretization：使用 piecewise constant FEM，使 φ_i φ_j=δ_ij φ_i，进而 diag(A Diag(q)B) 可表示为带 mass matrix 的 Schur product。它解决的是如何在离散层面保留 covariance diagonal 与 product operator 的代数结构。若换成高阶连续元，这个结构会明显变复杂。

4. Two-level low-rank compression：第一层压缩 observation map，第二层压缩 prior-preconditioned Khatri-Rao product space。它解决的是 online OED 中重复 PDE solve 和 dense n×n posterior update 的问题。核心变化是把 J 和 ∇J 的计算从 PDE-dependent 变成 low-rank algebra-dependent。

5. p-continuation binary design extraction：作者用 relaxed p-norm continuation 从 p=1 的非二值设计推向 binary sensor placement。它解决组合优化不可解的问题。但在本文设定下 objective 可能非凸，active OED 中的全局最优理论不能继承，因此这部分更像合理 heuristic，而不是论文最硬的理论贡献。

Key Insight / Why It Works

最重要的 insight 是：passive imaging 的“维度平方化”不是完全不可压缩的 m^2 black-box 数据扩张，而是 F 和 F 的双边作用诱导出的 structured product space。只要先找到 F 的主要 observation subspace，correlation data 的主要信息大概率落在该 subspace 的 Khatri-Rao/Kronecker 二阶组合里。因此第二层压缩不是任意降维，而是在第一层物理观测模态的 pairwise interaction space 中寻找 prior-relevant directions。

方法有效的真正原因有三点。第一，Helmholtz observation map 从源区域到有限传感器具有有限数值秩，尤其在固定频率、有限 aperture、有限候选点下，第一层压缩能捕捉绝大部分可观测源模式。第二，A-optimality 关注 posterior variance 的 trace，不要求逐项重建全部 covariance data，只需要捕捉 prior-weighted uncertainty reduction 的主方向；这使第二层低秩截断可行。第三，source diagonal 假设让 correlation operator 的二阶结构仍作用在 n 维 q 上，而不是 n^2 维 Q 上，否则低秩结构很可能不够。

最可能的核心贡献是 Theorem 3.6 这类 objective/gradient 的低秩表达：它把被动 OED 从“每步 PDE + 高维 posterior inversion”变成“预计算因子上的小矩阵微分”。这属于 memory reuse / offline-online decomposition，而不是更好的 statistical inference。p-continuation、SLSQP、随机低秩算法本身都不是新东西；它们是辅助。

需要直说的是，增益主要来自结构化 scaling，而不是更准确的被动成像统计建模。作者为了可计算性采用 Γ⊗Γ 的 Gaussian data noise，而真实 sample covariance 的 asymptotic covariance 依赖 q 和 w。也就是说，OED 优化的 utility 与真实有限样本 covariance experiment 之间存在错配。实验中设计优于随机并不意外，因为任何合理利用 PDE/geometry 的 OED 都应优于随机；更关键但文中未充分说明的是：这种简化 utility 选出的设计在真实 Wishart likelihood 或 source-dependent noise 下是否仍最优。

另一个有价值的机制性观察是 correlation data 改变了 sensor layout 的几何偏好。active data 中靠近源区的传感器常占优；correlation data 中过近的传感器可能测到高度冗余的 wave field，off-diagonal correlations 的独立信息有限，因此更分散的布设可能更好。这个 insight 值得迁移到其他 multi-sensor covariance sensing 问题：不要只最大化单点 SNR，还要设计传感器间的互信息/相关结构多样性。

Relation To Prior Work

最接近的技术谱系是 infinite-dimensional Bayesian linear inverse problem 的 A-optimal OED，尤其 Alexanderian/Petra/Stadler/Ghattas 一类 prior-preconditioned low-rank posterior update，以及 Aarset 关于 binary low-rank A-optimal designs 的 p-continuation / redundant-dominant 方法。本文在优化框架和算法外壳上明显继承这些 active imaging OED 工作。

本质差异在 passive imaging 的 correlation-data forward map。active case 的 posterior precision 大致是 F^T W F+C0^{-1}；这里变成 diag∘[(F^T W F)⊗(F^T W F)]∘Diag+C0^{-1}。设计变量进入次数更多，凸性丢失，且 product-space rank 可能平方增长。本文新增的信息是如何用 Schur/Khatri-Rao 结构把这个 passive precision update 压缩到两层低秩表示中。

看似新的部分中，Bayesian Gaussian posterior、A-optimal trace criterion、randomized low-rank approximation、Woodbury identity、p-continuation binary extraction 都是已有思想重组。实质创新在于把这些工具放到 passive covariance operator 的结构上，并给出可计算的 objective/gradient 公式。它不是 passive imaging inversion 理论的根本突破，也不是 OED 理论的全局最优突破；它是 correlation-data inverse problem 的 scalable OED formulation。

Dataset / Evaluation

评估是一个 2D Helmholtz 方程的合成 aeroacoustics-like 场景：圆形域、障碍物、源区、数百候选传感器、复值 wave field 转实值观测。任务覆盖的是固定几何、固定 PDE 类型、固定频率附近的线性 source covariance inversion。没有真实实验数据，没有多场景统计评估，也没有 3D 或多频泛化验证。

实验支持的 claim 主要是 computational scalability within this regime：offline 低秩压缩后可以高效评估 A-optimal objective/gradient，并找到比随机布设更好的 sensor designs。MAP 和 posterior samples 说明该 Bayesian formulation 在这个合成问题上能产生合理 reconstruction，但这不是 OED 核心 claim 的强证据。

评估没有充分验证的部分包括：有限样本 sample covariance 的真实 Wishart 噪声下设计是否仍好；Gaussian Γ⊗Γ 简化是否影响 sensor ranking；非凸 relaxed optimization 是否对初始化稳定；低秩截断误差是否系统影响 design；3D 场景中 offline O(n^2) product-space 操作是否仍可接受。因此实验足以证明“这个结构化算法在 2D 合成 Helmholtz 上能跑且优于随机”，但不足以证明“被动成像 OED 问题已被一般性解决”。

Limitation

最核心的限制是统计模型和计算模型都经过强约束。统计上，correlation data 来自 sample covariance，本应是 Wishart 或 asymptotic Gaussian，协方差依赖 F_w Diag(q)F_w^*+Γ，也依赖 w 和未知 q。本文用 Γ⊗Γ 作为固定 Gaussian covariance，这让 posterior 和 A-optimality 变得闭式，但也把真实实验设计中的 heteroscedastic/source-dependent uncertainty 抹掉了。文中承认这一点，但没有量化其对设计的影响。

建模上，spatially uncorrelated source 是 scalability 的关键前提。若要估计 full covariance Q，未知量变成 operator-valued，参数维度平方化回归，当前 two-level compression 很可能不够。换言之，论文解决的是 passive imaging 中最可压缩的子类，而不是一般 passive imaging source covariance problem。

优化上，J(w) 单调但不保证凸。monotonicity 只能说明加传感器不会带来负信息，不能排除局部极小或 relaxed-to-binary gap。p-continuation 在实验中工作良好，但理论上已从 guarantee 退化为 heuristic。非二值设计被当作 lower-bound-like reference 也只是经验现象。

scalability 上，方法并没有消除高维问题，而是把成本转移到 offline product-space compression。第二层 Khatri-Rao 相关操作带有 O(n^2ℓ\hatℓ) 级别瓶颈，\hatℓ 最坏可接近 ℓ^2。2D n≈几千时可行；3D、复杂几何、多频或高分辨率源区下，这个 offline 阶段可能成为新的主瓶颈。

还有一个隐含上限：低秩性取决于物理 forward map 和 prior 的 smoothing/limited aperture。若传感器数、频率、多路径复杂度或源结构使 effective rank 上升，online 低秩优势会被 \ell 和 \hat\ell 的增长吃掉。增益来源可能主要来自当前 PDE 设置的数值低秩，而非方法对任意 passive imaging 都天然 scalable。

Takeaway

• 1. 被动成像 OED 的关键不是简单把 active OED 的 data dimension 改成 m^2，而是识别 F⊗F 与 diag/Diag 之间的 product-space algebra；这个结构决定了是否可扩展。
• 2. 对 covariance/correlation data 问题，合理的 offline-online decomposition 往往需要两层压缩：先压缩物理观测模态，再压缩模态之间的二阶交互。
• 这个 insight 可迁移到 covariance sensing、array design、ambient noise imaging、多传感器统计反演。
• 3. 空间非相关源假设是本文成功的核心 inductive bias。

一句话总结

这篇论文是把 infinite-dimensional Bayesian A-optimal OED 推进到 passive correlation-data source imaging 的一次结构化 scaling 工作，真正贡献在于利用 diagonal source covariance 与 Schur/Khatri-Rao product-space 低秩压缩，把原本不可承受的被动成像传感器设计变成可离线预计算、在线小矩阵优化的问题。