[论文解读] On Solving a Generalized Chinese Remainder Theorem in the Presence of Remainder Errors

Guangwu Xu|arXiv (Cornell University)|Aug 30, 2014

Digital Filter Design and Implementation参考文献 7被引用 10

一句话总结

本文提出了两种高效算法，用于求解带有余数错误的广义中国剩余定理（CRT），在先前工作的基础上，通过简化商的计算并引入对受损余数的直接方法，实现了改进。关键贡献在于证明了误差界 τ = maxᵢ minⱼ gcd(Mᵢ,Mⱼ)/4 是紧确的，即无法进一步改进，并通过反例证明了在此界限之外会出现可识别性失败。

ABSTRACT

In estimating frequencies given that the signal waveforms are undersampled multiple times, Xia et. al. proposed to use a generalized version of Chinese remainder Theorem (CRT), where the moduli are $M_1, M_2, \cdots, M_k$ which are not necessarily pairwise coprime. If the errors of the corrupted remainders are within $τ=\sds \max_{1\le i\le k} \min_{\stackrel{1\le j\le k}{j eq i}} \frac{\gcd(M_i,M_j)}4$, their schemes can be used to construct an approximation of the solution to the generalized CRT with an error smaller than $τ$. Accurately finding the quotients is a critical ingredient in their approach. In this paper, we shall start with a faithful historical account of the generalized CRT. We then present two treatments of the problem of solving generalized CRT with erroneous remainders. The first treatment follows the route of Wang and Xia to find the quotients, but with a simplified process. The second treatment considers a simplified model of generalized CRT and takes a different approach by working on the corrupted remainders directly. This approach also reveals some useful information about the remainders by inspecting extreme values of the erroneous remainders modulo $4τ$. Both of our treatments produce efficient algorithms with essentially optimal performance. Finally, this paper constructs a counterexample to prove the sharpness of the error bound $τ$.

研究动机与目标

• 提供关于秦九韶广义CRT的历史忠实叙述，特别是其“大衍求一术”与“大衍求一法”的方法。
• 在存在余数误差的情况下，通过简化商的计算过程，改进现有的鲁棒CRT方法。
• 提出一种新颖的直接方法，直接作用于受损余数，无需预先估计商，并通过模 4τ 的极值分析揭示其结构特性。
• 证明误差界 τ = maxᵢ minⱼ gcd(Mᵢ,Mⱼ)/4 的紧确性，表明其在一般情况下无法被改进。
• 证明所提出的算法仅通过一次CRT计算即可实现渐近最优性能。

提出的方法

• 通过简化算法步骤，在保持正确性的前提下，简化了王和夏的商查找过程。
• 引入一种直接方法，作用于模 4τ 的受损余数，利用极值分析推断真实余数的信息。
• 利用模数 Mᵢ = d·mᵢ 的广义CRT模型结构，基于模 d 下受损余数的平均值推导出重构公式。
• 应用扩展欧几里得算法高效计算CRT系数，确保计算效率。
• 基于 a = floor(d/2) 的取值进行分类讨论，以处理不同的误差配置，并确保重构误差有界。
• 构造一个反例，使用两个不同的解 N₁ 和 N₂，在误差大小为 d/4 时产生完全相同的受损余数，从而证明该界限是紧确的。

实验结果

研究问题

• RQ1 在不牺牲正确性或性能的前提下，是否可以简化鲁棒CRT中商的计算步骤？
• RQ2 是否可能直接从受损余数中重构解，而无需估计商？
• RQ3 在广义CRT设定下，唯一重构仍可能成立的最大误差界限是多少？
• RQ4 误差界 τ = maxᵢ minⱼ gcd(Mᵢ,Mⱼ)/4 是否可以改进，或其本身已是最优？
• RQ5 通过分析模 4τ 下受损余数的极值，可以揭示余数的哪些结构性特征？

主要发现

• 所提出的第一个算法简化了商的计算，仅通过一次CRT计算即实现了渐近最优性能。
• 第二个算法直接作用于受损余数，并通过模 4τ 的极值分析揭示了有用信息。
• 误差界 τ = maxᵢ minⱼ gcd(Mᵢ,Mⱼ)/4 被证明是紧确的，反例表明当误差超过 d/4 时重构会失败。
• 构造了一个反例，其中两个不同的解 N₁ 和 N₂ 在误差大小为 d/4 时产生完全相同的受损余数，证明了在此界限之外存在不可识别性。
• 在所有情况下，重构误差严格小于 τ，确认该界限在给定假设下既是必要也是充分的。
• 结果表明，简化模型中的误差界 d/4 是最优的，即使在附加假设下也无法进一步改进。