精读笔记

Problem Setting

论文标题：Generalized matrix nearness problems II（arXiv preprint / 2026）。

这篇论文实际解决的是：在给定 A,B,C 时，求满足结构约束的 X，使 A-BXC 在某种范数下尽可能小；但它关心的不是一般优化建模，而是哪些广义 nearness 问题可以被还原为可由数值线性代数直接处理的问题。

真正困难点有两个。第一，约束集合往往非凸，例如 rank、Kronecker rank、 prescribed eigenvalue / partial trace 等；通用凸优化框架天然不适用。第二，Frobenius 范数之外的正交不变范数不再给出简单的正交投影几何，经典 Eckart–Young/Mirsky 类型的“同一个截断 SVD 同时适配所有正交不变范数”的现象在 generalized setting 中会断裂。

关键矛盾是：作者希望保留前作那种“不开梯度、不做黑箱优化、靠 SVD/投影/重排拿全局解”的风格，但问题扩展后，闭式结构只在少数情形幸存。于是论文一方面最大化挖掘可闭式化约的边界，另一方面为剩余情形构造一个仍然只依赖数值线性代数的全局收敛迭代框架。

Motivation

已有路线不够的地方很明确：前作 [13] 对 Frobenius 范数下的 generalized matrix nearness 给出大量闭式解或线性收敛算法，但这些结果本质上吃的是 Hilbert 空间结构。一旦换成 nuclear norm、spectral norm 或一般 Schatten p-norm，SVD 截断不再自动给出 generalized rank-constrained 问题的统一答案。

作者的核心观察是：并不是所有扩展都会破坏结构。有些看似更复杂的问题，只是坐标组织方式变了。例如 Kronecker rank 在 rearrangement 后就是普通 rank；partial trace 条件在 rearrangement 后可解释为特征值条件；仿射项可以通过对辅助变量做最小二乘消去，转成原问题上的投影版本。

关键缺口是：需要区分“结构上仍可化约”的推广和“只能迭代处理”的推广。本文的价值在于把这条边界画得更清楚：哪些新问题仍有闭式解，哪些不存在 Mirsky 型统一解，哪些可以通过 Dykstra splitting 保证全局收敛。

Core Idea

论文最核心的思想不是提出一个新的优化器，而是持续把 generalized nearness 问题重写成已知可解的几何对象。仿射项不是新自由度，而是可被正交投影消掉的 nuisance variable；Kronecker rank 不是新 rank，而是在 rearrangement operator 下的普通矩阵 rank；partial trace 不是孤立约束，而是 rearranged matrix 的 eigenvalue condition。这个思路的本质是：通过改变矩阵条目的组织方式，把非标准结构变成标准结构。

对于无法闭式求解的 Schatten 范数版本，作者没有走通用 nonsmooth optimization 路线，而是把问题拆成两个各自可线性代数化的子问题：一个是 Frobenius generalized nearness 子问题，直接调用前作/本文闭式子程序；另一个是 Schatten norm 的 proximal 子问题，在奇异值层面显式求解。Dykstra correction 保证两个步骤不是简单交替投影，从而避免振荡并得到全局收敛。

和 prior 的本质区别在于：prior 多在固定 Frobenius 目标下扩展约束；本文把“可解性”本身作为对象，系统展示在 Kronecker 结构和 Schatten 范数下哪些部分仍可继承、哪些部分必须放弃闭式统一性。

Method

1. 仿射项消去：它解决的是目标中额外 rank-one affine term 带来的自由变量问题。必要性在于如果直接把 x 和 X 一起优化，会破坏前作现成子程序；作者对 x 做最小二乘闭式消去，把问题变成 A 和 C 右乘一个去除 e 方向的投影矩阵后的普通 generalized nearness。核心变化是把 affine nuisance 变成数据预处理。

2. Kronecker rearrangement：它解决的是 Kronecker rank 约束不可直接用普通 SVD 的问题。通过定义 R(Y⊗Z)=vec(Z)vec(Y)^T，Kronecker rank 等价于 rearranged matrix 的普通 rank。必要性在于这把非线性的 Kronecker 分解问题变回 rank-constrained approximation。核心变化是把 tensor/Kronecker 结构编码为矩阵行列重排。

3. Partial trace 到 eigenvalue：它解决的是 quantum-style partial trace 约束如何纳入 generalized nearness 框架的问题。作者证明 tr_P(X(W⊗I))=λW 等价于 R(X)vec(W)=λvec(W)。必要性在于 prescribed eigenvalue 问题已有处理路径。核心变化是把一个 matrix-valued constraint 转成标准谱约束。

4. Mirsky-type negative result：它不是算法模块，而是边界判定。作者构造例子说明 rank-constrained generalized nearness 的 minimizer 会依赖 Schatten p。必要性在于阻止读者期待“Frobenius 闭式解自动推广到所有正交不变范数”。核心变化是把问题从“找统一公式”转为“接受范数依赖，用迭代”。

5. Dykstra/proximal splitting：它解决的是一般 Schatten norm 下缺少闭式解的问题。X-step 是 Frobenius nearness，Y-step 是 Schatten norm proximal，Δ 是 Dykstra correction。必要性在于简单 alternating 可能振荡；Δ 提供了误差记忆。核心变化是把全局 nonsmooth constrained problem 分裂成两个有结构的子问题。

Key Insight / Why It Works

最重要的 insight 是：很多 generalized nearness 问题的难度不在约束本身，而在表示。Kronecker rank 看起来是双线性分解复杂度，但 rearrangement 后就是普通 rank；partial trace 看起来是量子信息里的特殊算子约束，但 rearrangement 后就是 eigenpair 条件。这类结果的本质是 latent structure alignment：找到正确坐标系后，已有的 SVD/rank/eigenvalue machinery 直接复用。

闭式解部分真正有效的原因是 Frobenius norm 对条目重排、正交变换、最小二乘消元都极其友好。作者的闭式结果基本都在利用 Frobenius 几何的等距性和正交投影性质。这里不是 scaling，也不是更好的优化技巧，而是对问题代数结构的重新编码。

Schatten norm 迭代算法有效的原因则不同。它本质是 proximal/Dykstra splitting：把目标 f(A-BXC) 看成在集合 {A-BXC:X∈S} 上最小化 f。X-step 投到结构集合附近，Y-step 降低范数，Dykstra correction 保存先前投影误差，使迭代对应一个收敛的 splitting 过程。这里的核心贡献不是 proximal operator 本身，而是把前作 Frobenius nearness solver 插入 Dykstra 框架，使大量非凸结构约束也能被统一调用。

哪部分可能只是辅助？数值实验中相对 CVX 的巨大优势，在很大程度上来自 problem-specific linear algebra 和 CVX/SCS 对这类问题不匹配，而不是算法在一般意义上优于所有优化方法。尤其对于 p=3/2，CVX 无法直接建模并不说明该算法在所有非标准 convex objective 上都有绝对优势；它说明的是专用结构求解器比通用 conic modeling 更合适。

哪部分是实质贡献？我认为有三点：Kronecker rank generalized nearness 的闭式化约；partial trace 与 eigenvalue/rank constraint 的连接；以及明确否定 Mirsky 型范数无关推广。Algorithm 3 是实用上重要的整合，但数学新意更多在“哪些结构可化约、哪些不可”的边界判断。

Relation To Prior Work

最接近的是 Li and Lim 2023 的 Generalized matrix nearness problems。本文基本沿用其哲学：不要把问题交给通用优化器，而是利用 SVD、投影、矩阵分解直接解。差异在于前作主要在 Frobenius norm + 普通矩阵乘法下组织约束，本文把边界推到 affine、Kronecker 和 Schatten norm。

与 Eckart–Young/Mirsky 路线的关系也很关键。Kronecker rank approximation 的闭式解沿用 Van Loan–Pitsianis / tilde transform / rearrangement 的传统，本质上是 Eckart–Young 在重排空间中的应用；这部分“新”更多是把已有 Kronecker approximation 思想嵌入 BXC generalized nearness。真正新增的是带 B1⊗B2、C1⊗C2 的 generalized 版本及其与前作 rank-constrained solver 的连接。

与 convex optimization / proximal methods 的关系：Dykstra iteration 和 Schatten proximal 并不是新概念，nuclear norm soft-thresholding、spectral norm proximal 也属于成熟工具。本文的新增信息在于把这些工具和 generalized matrix nearness 的闭式 Frobenius 子程序组合起来，并证明在相当广的约束类上全局收敛。

所以这篇论文属于 numerical linear algebra flavored structured optimization 谱系，而不是现代大规模一阶优化谱系。它的本质创新不是提出更通用的 optimizer，而是把一批结构优化问题重新归约到可控的线性代数 primitives。

Dataset / Evaluation

实验覆盖三类：合成 matrix recovery、系统辨识中的 nuclear norm Hankel 问题、CFAR 检测中的 spectral norm PSD 问题。任务选择是合理的，因为它们分别检验了：标准/非凸约束下的恢复精度、多个约束组合、以及实际应用中 spectral norm generalized nearness 的可用性。

实验确实支持一个核心 claim：专用线性代数算法比 CVX/CVXPY+SCS 更快、更稳，尤其在非凸约束或非标准 Schatten p-norm 下，通用 conic solver 不是合适工具。这里验证的是“结构化 solver 的工程优势”，不是某种统计泛化能力。

但 evaluation 的边界也明显。合成实验通过反向构造 A=BXC，因此更像 sanity check 和数值稳定性测试，不是噪声、病态、模型错配下的系统评估。真实应用部分虽然来自系统辨识和 CFAR 场景，但仍是随机生成矩阵/模拟设置，并非完整真实部署 benchmark。PSD 约束那组尤其微妙：理论上正是收敛保证的例外，实验只能说明在作者生成的实例上表现好，不能补足理论缺口。

总体看，实验足够说明算法比 CVX 更适合这些结构问题，但不足以证明在所有大规模/病态/真实数据场景下都鲁棒。增益主要来自结构利用和避免 conic reformulation，而不是某种黑箱优化性能突破。

Limitation

1. 闭式解的适用边界很窄：Kronecker rank 和 partial trace 的漂亮化约依赖 B,C 具有 Kronecker 可分形式。一旦换成一般 B,C，作者明确没有闭式答案。这说明方法的核心能力来自代数结构匹配，而不是普适优化能力。

2. Frobenius 几何仍是地基：即使最终目标是 Schatten norm，Algorithm 3 的 X-step 仍回到 Frobenius generalized nearness。换句话说，非 Frobenius 情形的可解性是通过 splitting 间接获得的，而不是直接掌握了 Schatten geometry 下的 constrained nearness。

3. 收敛率上限较弱：理论只给 sublinear/Fejér 型下降。文中说实验接近线性，但没有给出解释。可能存在局部误差界或强正则性，但文中未充分说明。若实际应用中子问题病态或约束集合导致非唯一，收敛行为可能明显变差。

4. 每步 SVD 成本不可忽略：对于大矩阵，Algorithm 3 的 Y-step 要 SVD；Kronecker rank 化约也可能产生维度重排后的大矩阵。论文强调数值线性代数而非通用优化，但 scalability 上限仍受 dense SVD 支配。大规模稀疏/随机化 SVD 的处理文中未充分说明。

5. PSD 约束是理论缺口：作者承认 positive semidefinite case 可能不 attain minimum，因此不在 convergence guarantee 内。CFAR 实验使用的正是这个例外场景，所以该应用的理论支撑不完整。

6. CVX 对比的解释要谨慎：CVX/SCS 表现差不意外，因为它是通用 conic modeling + first-order solver，不是为这些矩阵结构定制的强基线。若对比专门的 proximal/ADMM/partial SVD solver，增益幅度可能不同。当前实验更多证明“不要用 CVX 解这类问题”，而不是证明 Algorithm 3 是最优工程实现。

Takeaway

• 1. 最值得记住的是：Kronecker 结构问题常常不是新优化问题，而是普通矩阵问题的错误坐标表示。
• rearrangement operator 是把 latent low-rank structure 显性化的关键工具。
• 2. Mirsky 型“范数无关最优解”在 generalized nearness 中失败，这是一个重要边界结果。
• 以后处理非-Frobenius generalized matrix approximation 时，不应期待 Frobenius 解自动迁移。

一句话总结

这篇论文把 generalized matrix nearness 从 Frobenius/普通乘法体系推进到 Kronecker 结构与 Schatten 范数体系，真正贡献是划清了哪些问题可通过线性代数重排闭式化约、哪些必须用结构化 splitting 迭代求解。