有限矩证书：重构、原子界与误差迁移 ​

返回学习路线。本页使用Prony重构理路Prony 有限幂和重构Prony finite-sum reconstruction · Classical Prony method · 普罗尼方法在不同节点与非零振幅的精确有限幂和模型中，从两倍阶数的连续样本恢复湮灭递推、节点和权重，并核清模型、重根及噪声边界。、平坦Hankel矩证书理路平坦 Hankel 矩的正测度证书Flat Hankel moment certificate · Univariate flat moment certificate · 一元平坦矩证书以前一阶矩阵正定和最后Schur余量为零，构造唯一的有限正原子测度；给出循环自伴证明、全部正测度中的唯一性及半正定和噪声边界。与Christoffel原子上界理路Christoffel 函数与矩约束原子上界Christoffel function moment bound · Christoffel variational principle · 克里斯托费尔函数以固定点取值为一的最小平方多项式定义Christoffel函数，从有限矩认证原子质量上界，并区分等号测度、不达上确界、支撑限制及矩误差。完成四种交付。每项都先说明矩表是准确给定还是只在误差盒里；“重构了一个候选”“全部正代表唯一”“某点原子质量有上界”是三种不同结论。

下载标准库精确程序与完整结果JSON。程序用有理数执行消元、矩阵分解、节点证书和平方预算。根候选由任务的准确因式分解给出；它不冒充任意多项式的根求解器。含 $\sqrt{3/5}$ 与 $\sqrt{3/7}$ 的对称节点通过平方关系准确核矩，未把根舍入为小数。

任务一：同一段数据，模型识别与全部正代表唯一各需要什么？ ​

给定六个准确值

$$(y_0,\dots ,y_5)=(1,1/2,3/2,5/2,11/2,21/2).$$

第一份合同说明它们来自三个不同节点、非零振幅的幂和。恢复节点、振幅、最小递推与下一项。第二份合同不预设原子数，只说输入应当是某个正测度的矩；再给 $m_6=43/2$，认证全部正代表中的唯一性。若这个最后矩改成 $43/2-1/100$，交一份明确的拒绝见证。

最后把前四项分别换成 $0,-1,-3,-7$ 和 $1,4,12,32$。前者能否用一般Prony处理？后者能否交付两个不同节点？

完整答案 ​

三阶Hankel矩阵与右端为

$$H=\left(\begin{array}{ccc}1& 1/2& 3/2\\ 1/2& 3/2& 5/2\\ 3/2& 5/2& 11/2\end{array}\right),h=(5/2,11/2,21/2{)}^{\mathsf{T}}.$$

顺序主子式为 $1,5/4,1$，LDL对角为 $1,5/4,4/5$。它可逆且正定。解 $Hc=h$得到 $c=(0,2,1)$，所以

$$q(t)=t^3-t^2-2t=t(t+1)(t-2),y_{k+3}=2y_{k+1}+y_{k+2}.$$

三个准确根为 $-1,0,2$；用 $k=0,1,2$ 的Vandermonde系统求得权重 $1/6,1/2,1/3$，再代入全部六项。三个根不同、权重非零，故这个完整幂和的最小递推阶数确为三。零节点分量只影响 $y_0$，不能把样本起点无声地改为 $k=1$。

递推预测 $y_6=2(11/2)+21/2=43/2$。在已知三分量模型中，这份预测有模型保证；若只给一段裸数据，它只是可供下一次检查的候选值。

现在以七个数作为正矩表。最后Schur余量为

$$s=m_6-h^{\mathsf{T}}c=43/2-(11+21/2)=0.$$

所以存在且仅存在

$$\mu =\frac{1}{6}{\delta }_{-1}+\frac{1}{2}{\delta }_0+\frac{1}{3}{\delta }_2.$$

唯一性没有预设候选只含三个原子：任意正代表都必须满足 $\int q^{2}d\mu =0$，从而集中在这三个准确零点，再由前三矩固定权重。只看前六矩的Prony识别，不能单独承担这一步。

若最后矩减去 $1/100$，其余矩不动，同一个多项式给

$$q(x{)}^2=x^6-2x^5-3x^4+4x^3+4x^2,L(q^2)=-1/100.$$

一个处处非负的平方却有负矩值，直接拒绝任何正代表。无需讨论小数特征值应当怎样设阈值。

两份四项迁移分别给：

• 对 $0,-1,-3,-7$，Hankel行列式为 $-1$，多项式为 $(t-1)(t-2)$，权重为 $1,-1$。一般Prony合同接受；正测度合同拒绝，因为总质量零而一阶矩非零。
• 对 $1,4,12,32$，系统可逆却得到 $(t-2{)}^2$。普通两不同节点合同拒绝；它实际对应 $(k+1)2^k$ 的合流项。重复填两个节点2不能修复问题，因为普通振幅相加只会得到一个几何分量。

任务二：一张均匀矩表，怎样证明4/9真是可达到的上界？ ​

给定准确的 $(m_0,\dots ,m_4)=(1,0,1/3,0,1/5)$。不假设未知测度均匀，只要求它正且匹配这些矩。求 $\mu (\{0\})$的最佳二次平方证书，并交一份达到的测度。再分别追加 $m_5=0,m_6=3/25$ 与 $m_5=0,m_6=1/7$，说明唯一性结论发生什么变化。对于第二种追加，必须给真正匹配到第六阶的有限正代表，不能只交三点规则的低阶准确性。

完整答案 ​

解 $H_{2}u=(1,0,0{)}^{\mathsf{T}}$ 得 $u=(9/4,0,-15/4)$。因此

$$p(x)=1-\frac{5}{3}{x}^2,L(p^2)=1-\frac{10}{3}\cdot \frac{1}{3}+\frac{25}{9}\cdot \frac{1}{5}=\frac{4}{9}.$$

$p(0)=1$且 $p^2\ge 0$，所以所有正代表都满足 $\mu (\{0\})\le 4/9$。这份上界可达到：令 ${\rho }^2=3/5$，取

$$\nu =\frac{4}{9}{\delta }_0+\frac{5}{18}{\delta }_{-\rho }+\frac{5}{18}{\delta }_{\rho }.$$

权重和为一，奇矩为零，二、四阶矩分别为 $\frac{5}{9}{\rho }^2=1/3$ 与 $\frac{5}{9}{\rho }^4=1/5$。非零原子处 $p=0$，故等号条件完整成立。这证明的是所有满足矩约束的正测度中，零点质量的最大值；它没有证明未知测度确实把 $4/9$ 放在零。

均匀分布在零点没有原子，却也匹配同一输入表。两份代表的第六矩不同：$\nu$给 $\frac{5}{9}(3/5{)}^3=3/25$，均匀分布给 $1/7$。

在第一种追加中，三阶湮灭多项式为 $x^3-\frac{3}{5}x$，最后Schur余量为零，因此 $\nu$是全部正代表中的唯一解。在第二种追加中，余量为 $1/7-3/25=4/175>0$，该平坦证书不通过；均匀分布仍是一个准确的连续代表。

第二种追加还允许五原子代表

$$\eta =\frac{16}{45}{\delta }_0+\frac{49}{180}({\delta }_{-\sqrt{3/7}}+{\delta }_{\sqrt{3/7}})+\frac{1}{20}({\delta }_{-1}+{\delta }_1).$$

总质量为 $16/45+49/90+1/10=1$，奇矩为零，$2j$阶矩为

$$\frac{49}{90}{\left(\frac{3}{7}\right)}^j+\frac{1}{10}=\frac{1}{2j+1},j=1,2,3.$$

这次连续与有限原子代表确实分享了整张到第六阶的表。正余量既不是无解证书，也不是“没有原子”的证书。若想继续判断，需要增加信息，而不是给小正余量强行取零。

任务三：同一个方差，为什么全实线没有最大者，而紧区间有？ ​

现在只给 $(m_0,m_1,m_2)=(1,0,1/3)$。先允许支撑在整条实线上，再限制支撑在 $[-1,1]$；分别求原点质量的上确界或最大值，并明确是否达到。把后一个问题平移到支撑 $[0,2]$、均值一、二阶矩 $4/3$，求点一的最大质量，检查支撑多项式怎样变化。

完整答案 ​

全实线情况下，Christoffel最优一次多项式是常数一，上界为一。但若原点质量等于一，测度只能是 ${\delta }_0$，与二阶矩 $1/3$矛盾，所以没有达到者。

对任意 $a\ge 1$，构造

$${\mu }_a=(1-\frac{1}{3a^2}){\delta }_0+\frac{1}{6a^2}({\delta }_{-a}+{\delta }_a).$$

其总质量一、均值零，二阶矩为 $2a^2/(6a^2)=1/3$。当 $a=1,2,4,8$，原点质量分别为 $2/3,11/12,47/48,191/192$，趋于一。于是上确界准确等于一，而没有最大者。每次都保留有限的二阶矩，但沿序列趋于无穷远的节点使二阶矩不随弱极限保留。

限制支撑在 $[-1,1]$ 后，逐点不等式

$$1-x^2\ge {\mathbf{1}}_{\{0\}}(x)$$

成立，积分给质量上界 $1-m_2=2/3$。$a=1$的上述测度正好达到，因此这里的准确答案是最大值 $2/3$。当 $a>1$时，逃逸族不再符合支撑合同，不能继续用来提高这个值。

平移为 $X=Y+1$，新矩为 $E[X]=1$、$E[X^2]=4/3$。在 $[0,2]$上改用

$$1-(x-1{)}^2=2x-x^2\ge {\mathbf{1}}_{\{1\}}(x).$$

积分仍为 $2-4/3=2/3$，并由 $\frac{1}{6}{\delta }_0+\frac{2}{3}{\delta }_1+\frac{1}{6}{\delta }_2$达到。若错误地继续套 $1-x^2$，它在区间右部为负且在目标点一为零，根本不是所需的指示函数上界。支撑与目标点改变后，应重写逐点证书，而非只替换矩值。

任务四：输入只有误差盒时，还能保留哪种结论？ ​

任务一的七个中心矩各有绝对误差至多 ${10}^{-6}$，真实测度只保证非负。能否据此认证真实测度恰好三原子？构造一个落在误差盒内的四原子测度作为检验。再用固定的 $q(x)=x(x+1)(x-2)$，认证集合 $\{|q|\ge 1/10\}$的质量上界。

另一份输入以任务二的前五矩为中心，总质量准确为一，二、四阶矩各有误差至多 $1/1000$；奇数阶也可有同样误差。交付零点原子质量的确定上界，说明哪些误差实际进入预算。

完整答案 ​

令任务一的三原子测度为 $\mu$，取 $\epsilon ={10}^{-9}$，定义

$$\tilde{\mu }=(1-\epsilon )\mu +\epsilon {\delta }_3.$$

它有四个不同正原子。第 $k$ 阶矩为 ${\tilde{m}}_k=(1-\epsilon )m_k+\epsilon 3^k$，对 $0\le k\le 6$，改变量都不超过 ${10}^{-6}$；其中最大绝对改变量是第六阶的 $707.5\times {10}^{-9}$。四阶大小的Hankel矩阵具有准确行列式

$$\det H_3(\tilde{m})=144(1-\epsilon {)}^3\epsilon >0.$$

因此真实秩可以从三增长到四，中心表的平坦性不能认证整个误差盒中的准确原子数。

固定平方多项式为

$$q^2=x^6-2x^5-3x^4+4x^3+4x^2.$$

中心矩给 $L(q^2)=0$，全部非零系数绝对值和为 $1+2+3+4+4=14$。于是任何符合误差合同的正测度都满足

$$0\le \int q^{2}d{\mu }_{\mathrm{t}\mathrm{r}\mathrm{u}\mathrm{e}}\le 14\times {10}^{-6}=\frac{7}{500000}.$$

在 $|q|\ge 1/10$的集合上，$q^2\ge 1/100$，故该集合质量至多

$$\frac{7}{500000}\mathbf{\nabla }\mathbf{\cdot }\frac{1}{100}=\frac{7}{5000}.$$

这份结论控制离开三个零点附近的质量，没有声明剩余质量全部是准确原子。构造的四原子测度还给 $\int q^{2}d\tilde{\mu }=144\epsilon$，可直接核它确实服从这个预算。

第二份输入继续使用任务二已经准确归一化的 $p(x)=1-\frac{5}{3}{x}^2$。因为

$$p^2=1-\frac{10}{3}{x}^2+\frac{25}{9}{x}^4,$$

总质量没有误差，奇数系数为零，实际误差预算只包含二、四阶：

$${\mu }_{\mathrm{t}\mathrm{r}\mathrm{u}\mathrm{e}}(\{0\})\le \frac{4}{9}+\frac{10}{3}\frac{1}{1000}+\frac{25}{9}\frac{1}{1000}=\frac{811}{1800}.$$

不需要对估计矩阵的逆作未经验证的浮点解释，检查者只需核多项式、归一化和逐矩误差。若这些误差只是概率意义下的置信区间，则结论也只在联合误差事件内成立；本任务没有额外声称抽样覆盖。

交付检查 ​

提交应保留原始矩表、阶数或次数、全部适用条件、准确多项式及矩值、节点与权重的完整性，以及失败或未决出口。有限程序检查支持具体算术与迁移，不代替一般定理。准确秩、模型内唯一性、全部正代表唯一性、上确界与达到者，都应按实际证据分别写清。