Large Signal Model
Foundations of realization, refinement and predictive closure
Komorebi P.
This mathematical body contains 87 numbered propositions with written proofs. The complete publication also includes its proof register, formal companion, glossary and references. The formal companion records 87 complete mappings, no partial mappings, and 0 propositions without a formalization. The exact scope and accepted execution evidence are recorded separately in the formal coverage map.
1. Three objects that must remain distinct
An observation is an element \(x\) of a declared space \(X\). Examples include a finite real signal, a complex field on a finite grid, or an array of categorical labels. Its type includes shape, scalar field or palette, and boundary convention. A zero-extended finite sequence and a periodic sequence are different objects even when their displayed samples agree.
An archive is a lossless representation of the accepted observation. In a mathematical model this means an injective analysis map with an explicit left inverse. In an implementation, conversion and floating-point error must be accounted for separately. An archive of a measured signal does not recover physical variables that were never observed.
A predictive state is a summary \(z=\gamma(x)\), or more generally a summary of an observation history. It is sufficient only relative to a specified family of future quantities, actions and dynamics. There is no implication from exact storage of \(x\) to sufficiency of \(\gamma(x)\).
The historical model's macro, meso and micro vectors are nested projections of a feature vector. They are not, merely by being nested, state variables on which the signal dynamics close. The distinction is retained in the present theory.
The term large denotes the structured hypothesis family: representations, resolutions, typed operators, compositions, memory mechanisms and their parameters. It is not a theorem about capacity, sample complexity or scale. Such a theorem needs a declared grammar, a complexity measure and a learning protocol.
2. Exact realization and atlas consistency
2.1 A left inverse is the reconstruction contract
Let \(X\) and \(C_a\) be finite-dimensional Hilbert spaces. A linear atlas consists of analysis and synthesis maps
The coefficients representing actual signals form \(\mathcal C_a=\operatorname{im}A_a\). This distinction matters for redundant representations, where \(\mathcal C_a\) is a proper subspace of the ambient coefficient space.
Proposition 1. Exact realization and consistency. The analysis map is injective, \(\Pi_a=A_aS_a\) is an idempotent with image \(\mathcal C_a\), and \(\Pi_ac=c\) exactly when \(c\) is a realizable coefficient vector. Moreover,
The first inequality is understood for nonzero \(X\). No assertion of orthogonality follows from a left inverse alone.
Proof. If \(A_ax=0\), then \(x=S_aA_ax=0\). Also \(\Pi_a^2=A_a(S_aA_a)S_a=\Pi_a\). Its image is contained in \(\mathcal C_a\) and it fixes every \(A_ax\), so its image is exactly \(\mathcal C_a\). Both norm inequalities follow from the operator-norm inequality applied to \(S_a\). \(\square\)
For a unitary atlas \(U_a\), the synthesis is \(U_a^*\) and the condition number is one. For a learned invertible atlas, exact invertibility is not enough: the synthesis norm controls amplification of coefficient errors. An arbitrarily ill-conditioned exact transform is not a stable representation.
2.2 Several exact atlases do not create new observations
For unitary \(U_a\) and positive weights \(w_a\), put
Then \(\mathcal S\mathcal A=I\) and \(\|\mathcal Ax\|^2=c\|x\|^2\), as proved in [K1]. Here \(\mathcal A\mathcal S=c^{-1}\mathcal A\mathcal A^*\) is the orthogonal projection onto the consistent coefficient subspace.
Each complete atlas contains the same observed information because it has a left inverse. Multiple atlases offer different sparsity, locality, conditioning and operator descriptions. They do not provide independent observations of the world. In particular, multiplying likelihoods from deterministic recodings as though they were independent measurements generally counts the same evidence more than once.
Categorical one-hot encoding is a separately typed injective map of sets. It is not a unitary map on the scalar-valued signal space. Its exact decoder is valid on one-hot arrays, not on arbitrary coefficient vectors. The linear tight-frame theorem does not apply to it without a separately specified linear embedding and inner product. This separation is already explicit in the capstone's categorical section.
3. Coordinate transport and one underlying operator
For a nonempty atlas family, define transport on realizable coefficients from atlas \(a\) to atlas \(b\) by
For an endomorphism \(F:X\to X\), define its coordinate realization by \(F_a=A_aFS_a\) on \(\mathcal C_a\). In this paragraph \(F\) may be nonlinear. For a partial map, every formula is restricted to the transported admissible domain. Maps between different signal types use separate input and output atlases.
Proposition 2. Coherent transport and operator realization. On the consistent coefficient spaces,
Conversely, suppose maps \(G_a:\mathcal C_a\to\mathcal C_a\) satisfy \(G_bT_{ba}=T_{ba}G_a\) for every pair. Then there is a unique map \(F:X\to X\) with \(G_a=A_aFS_a\) for every \(a\).
Proof. The transport identities follow by inserting \(S_bA_b=I\). The operator identity follows because either path is \(A_bFS_a\). For the converse, choose a reference atlas \(r\) and put \(F=S_rG_rA_r\). Compatibility gives \(G_aA_a=A_aF\). Every consistent coefficient is \(A_ax\), which proves the asserted realization. Applying \(S_a\) proves uniqueness. \(\square\)
An independently fitted family of atlas predictors need not satisfy this proposition. Its compatibility residual is a well-defined quantity,
After synthesis it becomes an object-space discrepancy:
Thus \(\|S_bd_{ba}(x)\|\leq\|S_b\|\,\|d_{ba}(x)\|\). Zero discrepancy establishes that the chart descriptions agree, not that their common prediction is correct. All charts can agree on the same wrong operator.
Redundant coefficient outputs must also be tested for consistency. Synthesizing an arbitrary coefficient vector always produces some signal, but reanalysis returns \(\Pi_ac\), not necessarily the original \(c\). Reconstruction and coefficient consistency are different checks.
4. Packets, refinement and transfer blocks
4.1 Packets belong to a declared resolution
Within one orthogonal atlas, let \(E_i\) be pairwise orthogonal coordinate selectors with \(\sum_iE_i=I\). Pull them back by a unitary map:
A packet is the component \(P_ix\) together with its projector, domain, atlas and resolution. A bare coefficient block without this information is not a fully specified packet.
Projectors from different atlases need not commute and must not be treated as one larger mutually orthogonal family. Weighted frame packets can sum to the identity while failing to be idempotent. These are different structures.
Proposition 3. Complete refinement preserves reconstruction. Suppose a finite rooted tree has root projector \(I\), and every internal projector is the sum of pairwise orthogonal child projectors. Every complete antichain \(\Lambda\), meeting each root-to-leaf path exactly once, satisfies
Proof. Replace an internal node by the sum of its children. Repeated replacement reaches any complete antichain and preserves the total projector. Different selected nodes have disjoint descendant branches, hence orthogonal ranges. The norm identity is Pythagoras. \(\square\)
This is a conservation law for information under complete refinement. Deleting a child, storing only its energy, or imposing a coefficient threshold changes the contract. For a retained subset \(J\) of an orthogonal decomposition, the lost signal is exactly \(\sum_{i\notin J}P_ix\), and its squared norm is the sum of the omitted packet energies.
The finite best-basis recursion
minimizes an additive cost over the declared finite tree [K1, R1]. The recurrence itself does not make a nonadditive predictive objective additive. If a prediction cost couples two leaves, a scalar cost at each node is generally insufficient for global optimization. Packet construction and cost evaluation must also be included in an implementation's complexity account.
4.2 Operations, not just signals, decompose
Let \(L:X\to Y\) be linear, with orthogonal input projectors \(P_i\) and output projectors \(Q_j\). Its transfer blocks are \(L_{ji}=Q_jLP_i\).
Proposition 4. Transfer reconstruction and refinement. One has
If \(P_i=\sum_rP_{i,r}\) and \(Q_j=\sum_sQ_{j,s}\) are complete refinements, then
For another linear map \(M:Y\to Z\) and output projectors \(R_k\),
Proof. Insert the appropriate identity decomposition on the input, output or intermediate space. Use \(Q_j^2=Q_j\) in the composition formula. \(\square\)
This is the algebraic reason to make transfer blocks first-class objects. The same operation can be expressed at several resolutions, while the coarse and fine expressions remain linked by exact identities. Approximations acquire explicit omitted-block terms rather than an unnamed representation error.
4.3 Nonlinear operations need interactions
For a nonlinear map, \(F(\sum_iP_ix)\) generally differs from \(\sum_iF(P_ix)\). Applying a linear transfer formula to such a map silently loses interactions.
Proposition 5. Polynomial packet interactions. Let
where each \(B_m:X^m\to Y\) is multilinear. Then
Proof. Substitute \(x=\sum_iP_ix\) in every argument and distribute by multilinearity. \(\square\)
For example, \((u+v)^2=u^2+2uv+v^2\). A grammar containing only independent packet squares misses \(2uv\). For nonpolynomial maps, one needs direct nonlinear semantics or a local expansion with a stated remainder bound. The word algebraic does not remove that obligation.
5. Predictive closure is the central contract
Let \(\gamma:X\to Z=\gamma(X)\) be a summary and let \(F_u:X\to X\) describe a deterministic transition for action \(u\). The reduced description closes when there is \(\bar F_u:Z\to Z\) such that
This is an equality of functions on the declared domain, not a fit statistic.
Proposition 6. Exact quotient criterion. Such a reduced transition exists if and only if
for every admissible \(x,x'\) and every action. When it exists it is unique on \(\gamma(X)\).
Proof. Necessity follows by applying \(\bar F_u\) to equal summaries. For sufficiency, define \(\bar F_u(z)=\gamma(F_u(x))\) for any representative with \(\gamma(x)=z\). The stated condition makes this well-defined and also proves uniqueness. \(\square\)
This criterion also applies when \(X\) is a latent physical state and \(\gamma\) is an observation map. Two physically different states can produce identical images and different futures. Exact image reconstruction does not resolve that hidden-state ambiguity. A history, a belief state or additional observations may be needed.
For a stochastic kernel, the corresponding condition is equality of the distributions of the next summary for states with the same current summary. Section 7 makes this finite criterion explicit.
5.1 The obstruction in the linear case
Suppose \(X\) is a Hilbert space, \(P\) is an orthogonal projection, \(Q=I-P\) and \(x_{t+1}=Lx_t\). The retained state is \(z_t=Px_t\). Then
Proposition 7. Closure defect and compositional defect. The retained state closes for every initial signal if and only if \(PLQ=0\). If it does not close, the error in using only \(PLP\) is
For two linear operators \(L\) and \(M\),
Consequently, compressing two operations separately can discard a path that leaves the retained subspace and later returns.
Proof. If \(PLQ=0\), the displayed recurrence depends only on \(z_t\). Conversely, apply Proposition 6 to \(x\) and \(Px\), which have the same retained state, to obtain \(PLQx=0\) for every \(x\). The norm bound is submultiplicativity. The compositional identity follows by inserting \(I=P+Q\) between \(L\) and \(M\). \(\square\)
The relevant quantity is not omitted energy alone. A very small omitted packet can have a large effect if its transfer into retained packets is strongly amplified. Conversely, an energetic packet can be irrelevant to a specified prediction if every corresponding return block is zero.
A residual is therefore typed. Reconstruction residual, coefficient-consistency residual, predictive residual, closure residual and statistical uncertainty should not be collapsed into one scalar without a specified objective.
6. Constructing the smallest closed linear refinement
There is a constructive finite solution when the dynamics are a known finite family of linear maps \(L_u\) on a \(d\)-dimensional space. It identifies which additional observables are required, rather than merely requesting a larger state vector.
Let the initial summary be \(C_0x\). Regard its rows as linear observables. Starting from \(W_0=\operatorname{row}(C_0)\), define
Proposition 8. Minimal linear predictive completion. The sequence stabilizes after at most \(d-\dim W_0\) strict enlargements. Its limit \(W_*\) is the smallest row space containing \(W_0\) and invariant under right multiplication by every \(L_u\). If \(R\) is any full-row-rank matrix whose rows form a basis of \(W_*\), then
for suitable matrices \(H\) and \(B_u\). Therefore \(z=Rx\) retains the original summary and has exact reduced dynamics \(z_{t+1}=B_uz_t\). Every linear summary with those two properties has row space containing \(W_*\).
Proof. The spaces are nested and have dimension at most \(d\). Equality at one step means invariance under every \(L_u\), so it is a fixed point. Containment in any invariant row space containing \(W_0\) follows by induction. The matrix identities express inclusion of the relevant rows in the basis \(R\). Conversely, if \(C_0=HR\) and \(RL_u=B_uR\), the row space of \(R\) is invariant and contains \(W_0\), so it contains every \(W_n\). \(\square\)
For one transition matrix this reduces to the span of the rows of
Its annihilator is the largest \(L\)-invariant subspace inside \(\ker C_0\), the standard unobservable-subspace construction [R2]. With several actions, observables pulled back through action words are required, but the linear span still stabilizes in finite dimension.
This gives a precise role for obstruction-driven refinement: add observable directions created by the dynamics until the chosen prediction variables close. It does not guarantee that the completion is small. In some systems \(W_*=X^*\) and exact memoryless closure needs the full state.
The theorem assumes the operators and relevant state are available. Applying it to estimated operators constructs closure for those estimates. It is not a proof that the real system has the same closure. If only \(C_0x\) is measured and the full \(x\) is unavailable, the augmented state \(Rx\) is not automatically accessible.
7. Joint closure of prediction and Bayesian inference
The finite categorical case admits an equally explicit construction that joins dynamics with Paper VI's inference refinement. This is useful because a state can be closed for prediction but too coarse for a new observation, or sufficient for an update but not for the next transition.
Let \(\Omega\) be a finite state set. Let \(K_u(i,j)\) be known row-stochastic transition matrices. Let \(\ell_y(i)>0\) be a finite family of observation likelihoods. A partition \(\pi\) has aggregation matrix
Distributions are column vectors. Prediction is \(p\mapsto K_u^\top p\), and a Bayesian update is
The coarse state must commute with both operations:
Fix a mandatory initial partition \(\pi_0\) with at least two blocks. This preserves distinctions that the intended application requires and excludes the trivial one-state quotient.
Proposition 9. Coarsest joint finite closure. First split \(\pi_0\) by the likelihood vectors \((\ell_y(i))_y\). Then repeatedly split each current block by the signatures
where \(B\) ranges over the current partition. The process terminates after at most \(|\Omega|-|\pi_{\mathrm{init}}|\) strict refinement rounds. The final partition is the unique coarsest refinement of \(\pi_0\) on which all declared predictions and Bayesian updates close for every prior.
Proof. At a stable partition, states in a common block have the same transition probability into every coarse block. Those probabilities define \(\bar K_u\) and give the first commutation identity, by checking it on point-mass priors. Likelihood constancy on each block defines \(\bar\ell_y\) and gives the second identity by direct summation.
For minimality, every admissible partition refining \(\pi_0\) has at least two blocks. By Proposition 12 below, exact updating for all priors requires each fine likelihood to be constant within each block. Thus it refines the first likelihood split. Suppose it refines a current partition. Each current block is a union of its blocks. Exact prediction makes the probability of reaching each such union constant within any admissible block, so it also refines the next split. Induction proves it refines the final partition. Each strict refinement increases the number of blocks, which proves termination. \(\square\)
The transition condition is the finite lumpability criterion [R4], expressed directly as a commutation identity. The joint construction adds the explicitly declared likelihood distinctions. Its output is a quotient of a known finite model, not an automatically learned model of an unknown world.
The accompanying rational example starts with six states. Likelihoods create three blocks, transition distinctions force two further rounds of splitting, and the final model has five states. Two original states remain genuinely equivalent for every action and observation in the declared family. The example verifies both commuting identities exactly.
8. Memory is an exact alternative to adding state
Refinement and memory are two ways to represent the influence of unresolved components. This can be derived rather than merely proposed as an architectural option.
For fixed linear dynamics and an orthogonal split, write
and define the restricted blocks
Thus \(z_{t+1}=Az_t+Bq_t\) and \(q_{t+1}=Cz_t+Dq_t\).
Proposition 10. Exact projected memory and a conditional tail bound. For every \(t\geq0\),
If \(\|D\|\leq\rho<1\) and \(\|z_s\|\leq Z\) over the relevant history, omitting terms with \(k\geq H\) changes the one-step right-hand side by at most
Dropping the initial unresolved-state term adds at most \(\|B\|\rho^t\|q_0\|\).
Proof. Iteration of the unresolved recurrence gives
Substitute this into the retained recurrence and change indices. Bound the omitted terms by the geometric series. \(\square\)
The memory kernel is \(BD^kC\). It records a path from retained variables into unresolved variables, evolution there, and return. The first kernel \(BC\) is precisely the return path appearing in the two-step composition defect. The term \(BD^tq_0\) is deterministic when \(q_0\) is specified. It becomes uncertain only relative to incomplete information or a probability model for \(q_0\).
This is the elementary linear instance of projection-memory reduction, associated with Mori-Zwanzig theory [R3]. A finite hand-chosen power-law or exponential memory feature is an approximation class for such effects, not a general derivation of the exact kernel.
The tail bound compares right-hand sides using the same resolved history. Replacing that history by recursively generated predictions introduces additional rollout error. Stability of that recursion needs its own bound. Also, a spectral-radius condition alone does not imply \(\|D\|\leq\rho\) in the norm already being used, especially for nonnormal operators.
9. Statistical refinement needs its own theorem
Nested feature vectors are algebraically nested when the smaller vector can be recovered from the larger by a fixed map. Predictive risk is a separate matter.
Proposition 11. Ideal statistical benefit of refinement. Let a target \(Y\) take values in a finite-dimensional Hilbert space with \(\mathbb E\|Y\|^2<\infty\). Let \(\mathcal F\subseteq\mathcal G\) be information sets and define \(m_{\mathcal F}=\mathbb E[Y\mid\mathcal F]\) and \(m_{\mathcal G}=\mathbb E[Y\mid\mathcal G]\). Then
Proof. Expand \(Y-m_{\mathcal F}=(Y-m_{\mathcal G})+(m_{\mathcal G}-m_{\mathcal F})\). The cross term has zero conditional expectation given \(\mathcal G\). \(\square\)
Thus finer information cannot worsen the optimal squared-error predictor under a fixed distribution. It can worsen an estimated predictor with finite data, restricted function classes, changed regularization or shifted distributions. The theorem supplies the ideal target for refinement, not a guarantee for every fitting procedure.
An exact archive makes later reanalysis possible without reconstructing discarded observations from a lossy state. It does not eliminate estimation error or distribution shift.
10. Bayesian coarsening and the information-loss packet
This section imports a specific existing Komorebi foundation rather than inventing another inference heuristic. Let \(p_i\) be a prior on a finite fine state set and let \(B_a\) be coarse blocks. Define
On zero-mass blocks the value of \(r_a\) is immaterial. Updating finely and then aggregating is exactly updating \(P\) with the effective likelihood \(r(p)\). The effective likelihood generally depends on the fine prior inside each coarse block, information that \(P\) alone does not retain.
Proposition 12. Bayesian closure and its exact defect. Suppose there are at least two coarse blocks and the fine and declared coarse likelihoods are strictly positive. Updating and coarsening commute for every prior if and only if
with one common \(c>0\). For a fixed prior, put \(\delta=r(p)-\bar\ell\), \(\bar b=\mathcal B_{\bar\ell}(P)\) and \(Z_r=r(p)^\top P\). Then
The fixed-prior identity also holds for nonnegative likelihoods when both normalizing evidences are positive.
Proof. Constancy up to one global factor gives commutation by summing the fine posterior. Conversely, compare priors supported on one state in each of two distinct coarse blocks. Equality of coarse posterior odds forces the ratio of those states' likelihoods to equal the ratio of the coarse likelihoods. Varying the representatives proves constancy within each block and the common factor across blocks. For the defect, use \(Z_r=\bar\ell^\top P+P^\top\delta\) and subtract the two normalized vectors. \(\square\)
With only one coarse block, both posteriors aggregate to one regardless of the likelihood, so the necessity criterion is false there. This boundary is explicitly retained from the corrected Paper VI statement [K2].
The matrix in the defect formula annihilates \(\bar\ell\), recording the irrelevant common-scale direction. It does not make every likelihood discrepancy irrelevant. The construction in Proposition 9 refines until both the likelihood and transition defects vanish for a specified finite family.
10.1 Parameter inference is not program truth
Within one declared regression model, a normal-inverse-gamma prior gives a standard conditional Bayesian parameter posterior. The software's program weights instead use penalized held-out predictive scores and then may refit coefficients on all available training records. These are finite-candidate predictive weights. Without a corresponding generative derivation, they should not be conflated with the marginal-likelihood posterior over programs under one fixed prior.
Every candidate must be scored on the same target event and base measure. Densities on differently sized coordinate spaces are not directly comparable. An invertible nonlinear coordinate change requires the appropriate Jacobian. A redundant deterministic embedding induces a measure on its consistent subspace, not an independent full-dimensional density on all coefficient entries.
The law of total variance separates within-candidate uncertainty from disagreement of candidate means. It does not supply a probability that the truth lies outside the grammar. An explicit outside-grammar component needs its own prior and predictive distribution. Assigning an arbitrary leftover mass does not solve open-world inference.
11. What a typed operator language must guarantee
A program has a source type, target type, admissible domain, evaluation rule and declared approximation status. A nonlinear partial operation is not a total matrix merely because its inputs are stored in a vector. The domain of a composite is
Composition must preserve those types and domains. Coordinate changes must preserve the represented operation as in Proposition 2. Packet refinement must preserve exact transfer semantics as in Proposition 4, or carry an explicit approximation residual. Nonlinear interactions must retain their cross-packet terms as in Proposition 5.
A grammar can be closed under a specified set of constructors without containing every possible physical law. It can also be expressive while computationally intractable. Universal discovery, practical search and semantic closure are distinct claims.
For a known operator, Proposition 8 can build its minimal linear summary. For an unknown operator, learning must address identifiability: which different candidates make identical predictions under the available observations and actions? More coordinate systems alone do not distinguish candidates that agree on the entire observed signal. Additional interventions or data may be needed.
The foundation should therefore treat a learned hypothesis as an operator with an evidence record, not just a matrix with a high normalized weight. Its record includes fitting assumptions, admissible inputs, compared alternatives, unresolved equivalences and the representation in which its residual is measured.
12. A bounded infinite-resolution bridge
Finite reconstruction identities do not alone prove convergence of an unbounded or infinite-dimensional model. A useful bounded extension can nevertheless be established directly.
Proposition 13. Strong fixed-horizon refinement limit. Let \(X\) be a Hilbert space, \(L:X\to X\) bounded, and \(P_n\) orthogonal projections converging strongly to \(I\). Put \(L_n=P_nLP_n\). For every fixed \(x\in X\) and every fixed nonnegative integer \(k\),
Proof. Since \(\|P_n\|\leq1\),
Also \(\|L_n\|\leq\|L\|\). Induct on \(k\), using
The first term tends to zero by the uniform norm bound and induction, and the second by strong convergence on the fixed vector \(L^kx\). \(\square\)
If estimated operators \(\widehat L_n\) are uniformly bounded and \((\widehat L_n-L_n)x\to0\) for every fixed \(x\), the same argument gives fixed-horizon convergence for \(\widehat L_n\). These are explicit assumptions, not conclusions of the finite learning procedure.
This result is strong convergence, not operator-norm convergence and not a bound uniform in time. On \(\ell^2\), finite coordinate projections converge strongly to \(I\) while \(\|I-P_n\|=1\) for every \(n\). Unbounded operators require domain control. Adaptive data-dependent partitions and fitted coefficients require additional consistency arguments. These remain separate mathematical tasks.
13. A complete small signal example
Consider the rational linear signal dynamics
The spatial archive stores \(x_t\). A second atlas uses the orthogonal Haar matrix
Its coefficients \(c_t=Ux_t\) reconstruct as \(x_t=U^*c_t\). The same dynamics in that atlas are \(ULU^*\), not a separate physical law. Thus coordinate transport preserves the full dynamics.
Now retain only the first spatial coordinate. With \(P=\operatorname{diag}(1,0)\), the retained scalar \(z_t\) and discarded scalar \(q_t\) obey
The exact first-coordinate sequence begins
The compressed one-step matrix predicts \(z_2=1/4\) instead of \(1/2\). The lost return path contributes exactly
There are two exact repairs. One is state refinement: start with \(C_0=(1,0)\) and add \(C_0L=(1/2,1)\). The resulting row space is the full two-dimensional space, so it closes. The other is memory:
since \(q_0=0\). This is the exact same system expressed without retaining \(q_t\) as a separate current-state variable. At \(t=2\), it yields
To see a nontrivial reduced closure, adjoin a third independent coordinate with multiplier \(5/6\). The smallest closed refinement of the first coordinate then contains the first two coordinates, not the third. Exact prediction of the chosen summary does not require reconstructing every component of the enlarged signal.
Finally, low energy is not the same as low predictive importance. For
omitting the second coordinate loses norm \(1/1000\), but changes the next first coordinate by exactly one. The correct certificate is the product of omitted magnitude and transfer amplification.
14. Connection to observable learning
Exact realization, complete decomposition and operator transport provide the representation laws. Propositions 6–10 separate predictive closure, constructive refinement and memory. The remaining requirement for learning is access to a sufficient predictive state. Part II addresses a finite linear observable sector; it does not infer its rank or coverage assumptions from a successful fit.
\newpage
Part II. Learning, stability and coherent expansion
The first stage constructed predictive closures for declared operators. The second stage asks how those constructions relate to what a learner can actually observe. A closure that depends on an unavailable hidden coordinate is a mathematical description, not yet an accessible state. A fitted recurrence can be accessible but wrong outside the span represented by its data. A one-step approximation can be accurate and still fail under repeated composition.
The results below connect these distinctions. They give an exact learning result for a finite linear observable sector, deterministic perturbation bounds, a propagation law for typed defects, a sufficient stability condition for truncated memory, controlled nonlinear interaction bounds, coherent atlas reconciliation and posterior accounting for an expanding grammar. These are finite or explicitly bounded results. They do not assume that every unknown process has a short linear state.
The observable-state perspective has classical antecedents in predictive state representations [R5] and observable-operator learning [R7]. Finite invariant observable spaces also underlie Koopman reductions [R6]. The present development uses elementary statements proved below and identifies their specific place in Komorebi's realization and obstruction calculus. It does not import those papers' sample-complexity guarantees into this model.
15. From hidden state to an accessible predictive state
Let a finite-dimensional autonomous linear process have hidden state \(x_t\in\mathbb R^d\) and observed output \(y_t\in\mathbb R^p\):
The observation map need not be injective. For a positive integer \(h\), define the observation map
The coordinates of \(O_hx_s\) are future observations relative to time \(s\). They are not available at time \(s\) merely because this matrix is written down. At time \(t\geq h-1\), however, the same vector is the observed history
Thus a forward observation map becomes a causal state after a declared warm-up interval. The time index is part of the contract.
Proposition 14. Exact causal delay closure. The following statements are equivalent:
- There is a linear map \(D_h:V_h\to V_h\) satisfying \(O_hL=D_hO_h\).
- \(\ker O_h\subseteq\ker(O_hL)\).
- \(\operatorname{rank}O_{h+1}=\operatorname{rank}O_h\).
When these statements hold, \(D_h\) is unique on \(V_h\). The history state obeys \(z_{t+1}=D_hz_t\), and its last output block is \(y_t\). Its iterates therefore determine every subsequent observed output. If a positive integer \(d\) bounds the hidden dimension, \(h=d\) is sufficient for closure, though a smaller positive \(h\) may suffice.
Proof. The equivalence of the first two statements is the linear factorization criterion: define \(D_h(O_hx)=O_hLx\) and require independence from the representative \(x\). The rows of \(O_hL\) are \(CL,\ldots,CL^h\). All except the final block already occur in \(O_h\). The kernel condition is therefore equivalent to \(CL^h\) being in the row space of \(O_h\), which is exactly the rank equality. It also implies \(O_hLx\in V_h\) because this vector is \(O_h\) applied to \(Lx\). Uniqueness follows because every element of \(V_h\) has a representative. The history recurrence follows directly. Cayley-Hamilton expresses \(L^d\) as a linear combination of lower powers, so \(CL^d\) is in the row space of \(O_d\). \(\square\)
This does not recover hidden variables in \(\ker O_h\). It identifies the quotient needed for the declared outputs. An independent hidden component can remain unobservable forever without preventing exact output prediction.
The rank equality here is a property of the true observation maps. A plateau in a finite empirical data matrix is not the same statement. It may reflect unexcited directions rather than a closed observation space. Section 17 makes the difference exact.
For several actions \(L_u\), a fixed delay vector does not automatically have one autonomous transition matrix. Products of the actual action matrices determine its observation map, which can change with the action history. One needs an action-aware state, a common closed observable family, or a separately justified belief-state construction. Proposition 14 is not silently generalized to arbitrary controlled or stochastic dynamics.
16. What finite data identify
Suppose candidate states and their successors have been recorded as columns of matrices \(Z_-,Z_+\in\mathbb R^{m\times N}\). Their construction may use output histories. No hidden-state coordinates need be supplied to the fitting procedure.
Write \(W=\operatorname{im}Z_-\). Exact interpolation of these columns is a question about a linear map on \(W\), not yet about every physically realizable state.
Proposition 15. Data-span identification and its boundary. A linear map \(B\) on the ambient space with \(BZ_-=Z_+\) exists if and only if
Its restriction to \(W\) is unique. It is an endomorphism of \(W\) if and only if \(\operatorname{im}Z_+\subseteq W\). A particular ambient extension is
where \(\dagger\) denotes the Moore-Penrose inverse. Every other fitted extension agrees with \(B_0\) on \(W\). If \(W\) is a proper subspace of the true realizable state space \(V\), the data alone do not identify the operator on \(V\setminus W\).
Proof. Necessity follows by applying \(B\) to every linear dependence among the input columns. Under the kernel condition, define \(B(Z_-a)=Z_+a\). Two choices of \(a\) differ by a kernel vector and give the same answer. This defines and uniquely determines the map on \(W\). It maps \(W\) into \(W\) exactly under the image condition. The pseudoinverse formula gives that map because \(a-Z_-^\dagger Z_-a\in\ker Z_-\). If \(W\subsetneq V\), choose a nonzero linear functional \(\varphi\) on \(V\) vanishing on \(W\) and a nonzero \(v\in V\). The operators \(B\) and \(B+v\varphi\) agree on every observed column and differ elsewhere. \(\square\)
The three checks are different. The kernel test checks consistency of a deterministic linear fit. The image test checks whether the observed span is closed under the fitted transitions. Equality \(W=V\) is a coverage condition. Neither of the first two checks establishes it.
For a known fixed feature space of dimension \(r\), coverage is witnessed by \(r\) independent input states. For an unknown larger feature space, the same rank only gives a lower bound on dimension. A model card must say which interpretation is being used.
If several action-specific operators act on a common declared \(r\)-dimensional state, each action's paired data must determine that operator on the required span. Pooling all actions into one high-rank matrix does not identify an action that has never been tried on some direction.
16.1 Perturbations have an explicit amplification factor
A noisy rank test is not an exact certificate. Even after the state space is fixed, fitting can be ill-conditioned. The following deterministic bound isolates that issue without introducing a probability model for noise.
Work in a fixed \(r\)-dimensional coordinate system, where \(Z_-\in\mathbb R^{r\times N}\) has full row rank. Let
The matrix \(R\) records deviations from the declared linear recurrence on these samples. The matrices \(E\) and \(F\) record input and successor perturbations. All norms in the next proposition are spectral norms.
Proposition 16. Fixed-space identification bound. Put \(s=\sigma_{\min}(Z_-)>0\). If \(\|E\|<s\), then \(\widetilde Z_-\) has full row rank, and
satisfies the exact error identity and bound
Proof. For every unit \(v\in\mathbb R^r\), \(\|\widetilde Z_-^\top v\|\geq\|Z_-^\top v\|-\|E^\top v\|\geq s-\|E\|\). Thus the perturbed matrix has full row rank and its pseudoinverse has norm at most \((s-\|E\|)^{-1}\). Since \(\widetilde Z_-\widetilde Z_-^\dagger=I_r\), subtract \(B\) from the fitted expression and collect terms. The bound follows by submultiplicativity and the triangle inequality. \(\square\)
This is a bound for a fixed state space and a fixed coordinate basis. If the basis or its dimension is also estimated from noisy data, subspace-estimation error must be added. No such error is hidden inside \(R\) without a derivation. Correlated and adversarial perturbations are allowed by the displayed deterministic inequality, but there is no sample-complexity claim until their magnitudes are bounded under a sampling model.
The denominator explains why numerical rank alone is insufficient. Full rank can coexist with severe noise amplification. The reference calculations include examples attaining the \(1/s\) dependence and a rank-collapse example at the excluded boundary \(\|E\|=s\).
17. Why finite agreement is not unrestricted discovery
Proposition 17. Finite-prefix indistinguishability. For every integer \(T\geq0\), there are two finite-dimensional autonomous rational linear systems with identical scalar observations at times \(0,\ldots,T\) and different observations at time \(T+1\). Both transition matrices can be chosen nilpotent.
Proof. One system has scalar output identically zero. For the other, take dimension \(T+2\), initial state \(x_0=e_{T+2}\), observation row \(C=e_1^\top\), and shift \(Le_j=e_{j-1}\) for \(j\geq2\), \(Le_1=0\). Then \(CL^tx_0=0\) for \(0\leq t\leq T\) and \(CL^{T+1}x_0=1\). The shift is rational and nilpotent. \(\square\)
An invisible delayed component of this kind can also be adjoined to any fixed base model. The combined output agrees with the base output on the entire inspected prefix and then differs. This is not an argument against learning. It identifies why a claimed guarantee must state a dimension bound, a model class, a coverage condition, an intervention protocol or a probabilistic assumption.
Even nilpotent operators have only zero eigenvalues, so stability language alone does not remove this ambiguity. A finite rank plateau cannot certify the absence of a delayed mode when arbitrarily large hidden dimensions are admitted.
18. A complete output-only identification example
Return to the rational system from Section 13, now with an additional independent hidden coordinate:
The third coordinate never influences the observed output. The first coordinate alone is not closed, but
has rank two and \(O_3\) has the same rank. Proposition 14 gives an exact two-dimensional delay state. Choose \(x_0=(1,0,q)^\top\) for any \(q\). The first four observations are
A learner receives those four outputs and the declared delay length \(h=2\). It does not receive \(L\), \(C\), \(O_2\) or \(q\). From consecutive observed windows it forms
Since \(\det Z_-=1/4\), these windows span the declared two-dimensional state. Exact arithmetic gives
The learned recurrence is
It predicts \(y_4=61/144\), which was not used to identify the coefficients. Proposition 14 and the verified rank condition on the true observation maps establish why this recurrence continues for every initial state of this particular system. The four samples alone establish the fitted recurrence only within the declared closed two-state model class. Proposition 17 forbids confusing those two conclusions.
The reference script separates the roles. One routine fits only the observed samples. A different check, with access to the generating system, verifies the theoretical assumptions and compares subsequent outputs. This prevents hidden access to the true operator from being presented as learning.
The same example connects back to the exact memory kernel. Section 13 expressed the output through a decaying convolution of its history. Here the history is summarized by a two-component recursive state. Neither representation stores the independent third coordinate. Memory and state realization describe the same observable behavior in two different ways.
19. A single propagation law for typed defects
A reconstruction error, an operator-fit error and a refinement error have different origins. Once each has been translated into an error in the state used by the next operation, their propagation can be analyzed by the same identity.
Allow the representation to change with time. Let \(\gamma_t:X\to Z_t\) be the declared representation at step \(t\), with a norm on each \(Z_t\). Let \(x_{t+1}=F_t(x_t)\) be the true transition and let \(\widehat f_t:Z_t\to Z_{t+1}\) be the implemented step. Its defect is
This definition already includes a representation switch. It is not restricted to one fixed feature vector.
Proposition 18. Compositional defect propagation. Suppose \(\widehat f_t\) is \(\lambda_t\)-Lipschitz on a domain containing both \(\gamma_t(x_t)\) and the predicted state \(\widehat z_t\). Suppose \(\|r_t(x_t)\|\leq\delta_t\) and \(\widehat z_{t+1}=\widehat f_t(\widehat z_t)\). Then, with \(e_t=\|\gamma_t(x_t)-\widehat z_t\|\),
Empty products are one. If the output \(h_t(x)\) has a decoder \(H_t\) with Lipschitz constant \(\kappa_t\) and \(\|h_t(x_t)-H_t(\gamma_t(x_t))\|\leq\eta_t\), then
Proof. Add and subtract \(\widehat f_t(\gamma_t(x_t))\) to obtain \(e_{t+1}\leq\delta_t+\lambda_te_t\). Repeated substitution yields the product-sum formula. Add and subtract \(H_t(\gamma_t(x_t))\) for the output inequality. \(\square\)
For a fixed linear representation \(R\) and an estimated transition \(\widehat B_u\), the defect is \((RL_u-\widehat B_uR)x\). For a reference reduced operator \(B_u\), it splits exactly as
The first term is representational nonclosure. The second is operator-estimation error. A fit statistic for the second cannot certify that the first vanishes. Proposition 16 supplies a conditional bound on the second in its declared fixed-space setting.
If \(\lambda_t\leq\rho<1\) and \(\delta_t\leq\delta\), then
This is a genuinely recursive prediction bound. It is not merely a one-step comparison made using the true history at every step. If the Lipschitz constants exceed one, the same theorem remains valid but can permit error growth.
19.1 Stability under changing operations
For linear candidate operators, one sufficient way to obtain a common contraction norm is to exhibit a positive-definite matrix \(M\) and \(0\leq\rho<1\) with
Then \(\|z\|_M=(z^\top Mz)^{1/2}\) gives \(\|\widehat B_uz\|_M\leq\rho\|z\|_M\) for all actions, directly from the displayed inequality. Proposition 18 therefore applies under arbitrary switching among them. This is the common quadratic Lyapunov construction from control theory [R8], used here as an explicit certificate rather than inferred from a list of eigenvalues.
Individual spectral stability does not suffice. The matrices
are both nilpotent, but \(BA=\operatorname{diag}(0,4)\). Alternating them amplifies an appropriate state by four every two steps. Even a single stable nonnormal matrix can produce large transient amplification in a specified norm.
If a learned policy chooses actions from predicted states while the reference process uses actions chosen from true states, the two action sequences may differ. Proposition 18 then needs a bound for that action mismatch, or a Lipschitz bound for the full closed-loop map. The result is not an automatic control-policy guarantee.
20. Truncating memory without losing the prediction bound
Section 8 bounded the omitted tail of the exact linear memory expression using the same resolved history. We now close the separate recursion argument.
Write the exact relation as
For negative indices set \(z_s=0\), solely as a convention for missing prehistory in this initialized formula. A length-\(H\) approximate model, with \(H\geq0\), is
An empty sum is zero. Let \(d_t\) be the exact right-hand side minus this truncated right-hand side evaluated on the true history.
Proposition 19. Uniform prediction bound for a stable finite-memory truncation. Suppose
Let \(E_0\) bound the initialization errors at indices \(-H,\ldots,0\). Then the recursively generated predictions satisfy
If \(\|z_t\|\leq Z\) for \(t\geq0\) and \(\|D\|\leq\rho<1\), a valid choice of defect bound is
The last term can be removed if the exact initialized forcing \(f_t\) is retained in the approximate recurrence.
Proof. Subtract the two recurrences and bound every retained history error by the largest of the errors at indices \(t-H,\ldots,t\). This gives \(e_{t+1}\leq a\max_{t-H\leq s\leq t}e_s+\Delta\). The ball of radius \(E=\max(E_0,\Delta/(1-a))\) is invariant because \(aE+\Delta\leq E\). Induction proves the bound. The estimate for \(\Delta\) separates coefficient errors, the geometric kernel tail and the initialized forcing, using Proposition 10. \(\square\)
This norm-sum condition is sufficient, not necessary. A stable recurrence with cancellation can fail this test. Failure to find this certificate is not proof of instability.
For the scalar example with \(q_0=0\), take exact retained coefficients and \(z_0=1\). All outputs lie in \([0,1]\): the full nonnegative memory coefficients sum to at most \(7/8\), so induction bounds every new output by the previous maximum. Thus \(Z=1\) is valid. With
one obtains
The right-hand side is a theorem for every time in this initialized example. The figure compares it with a separately computed finite-horizon error. The computed maximum is not substituted for the uniform bound.
21. Nonlinear packet truncation with an explicit remainder
Polynomial interactions were expanded exactly in Proposition 5. A computational model may retain only some ordered interaction tuples. It then needs a bound on what was omitted, rather than an assumption that packetwise evaluation is exact.
Let \(P_i\) be an orthogonal resolution of a finite-dimensional real Hilbert space. Suppose
with bounded multilinear maps. For each ordered tuple \(\boldsymbol i=(i_1,\ldots,i_m)\), let \(\beta_{m,\boldsymbol i}\) bound the norm of \(B_m\) restricted to \(\operatorname{im}P_{i_1}\times\cdots\times\operatorname{im}P_{i_m}\). Let \(\mathcal I_m\) be the retained tuple set, and let \(F_{\mathcal I}\) retain the constant term and exactly those interactions.
Proposition 20. Packet-interaction truncation bound. Put \(a_i=\|P_ix\|\). Then
On the ball \(\|x\|\leq R\), a valid Lipschitz bound for \(F_{\mathcal I}\) is
Proof. Subtract the retained terms from the exact multilinear expansion. Bound each omitted term by its restricted multilinear norm and sum. For the Lipschitz estimate, telescope the difference between two \(m\)-fold evaluations by replacing one argument at a time. Each changed argument has norm at most \(\|x-y\|\) after projection, and each unchanged argument has norm at most \(R\). There are \(m\) terms for each retained tuple. \(\square\)
If the declared polynomial itself approximates a nonpolynomial map with remainder at most \(\varepsilon_R\) on the ball, add \(\varepsilon_R\) to the defect bound. Without such a remainder estimate, the result is only a theorem for the polynomial, not for the unexpanded physical law.
For a retained state \(z=Px\) and omitted state \(q=(I-P)x\), dropping every term that uses \(q\) gives the coarser majorant
The tuple-level result is finer because it can distinguish a harmless omitted packet from a strongly coupled one. Ordered tuples must be counted with their multiplicities unless symmetry has explicitly combined them.
A prediction guarantee follows by inserting these defect and Lipschitz bounds into Proposition 18, provided the true and predicted trajectories remain in the declared domain. That domain condition cannot be inferred from a small initial error. One possible construction is to assume that the true map preserves a closed convex set \(K\) and to project every approximate output back onto \(K\). Euclidean projection is nonexpansive and fixes the true output, so it does not increase the one-step approximation bound. This is a declared modification of the predictor, not an identity for an unmodified recurrence.
For the scalar quadratic \((u+v)^2\), retaining only \(u^2+v^2\) omits \(2uv\). The packet bound is \(2|u||v|\) and is exact in magnitude. Small energy in one packet matters in proportion to both its amplitude and the cross-interaction coefficient.
22. Learning several charts without inventing several worlds
The exact atlas identities constrain the meaning of a fitted model, but independently trained predictors need not satisfy them. Here is a simple coherent reconciliation with a precise objective.
Fix a finite-dimensional represented space \(X\). It may be the observed signal space or a predictive quotient, but it must be the same space for all atlases in this section. Let \(A_a:X\to C_a\) and \(S_a:C_a\to X\) satisfy \(S_aA_a=I\). Let \(G_a:C_a\to C_a\) be linear fitted maps. Pull them back to object space:
For positive weights summing to one, define
Proposition 21. Coherent reconciliation and coordinate-stable transport. The family \(\widetilde G_a\) is consistent and satisfies every atlas intertwining identity. The operator \(\overline L\) is the unique minimizer of
over linear maps on \(X\), where \(\|\cdot\|_F\) is the Frobenius norm. For any target operator \(L_*\) and any chosen operator norm,
Further, the norm \(\|c\|_a=\|S_ac\|\) on \(\mathcal C_a\) makes every exact transport \(T_{ba}=A_bS_a\) an isometry. A represented transition \(A_bLS_a\) has exactly the same induced Lipschitz constant in these norms as \(L\) has in object space.
Proof. Coherence follows from Proposition 2 because every chart realizes the same \(\overline L\). Expanding the squares gives
The weighted cross term vanishes, proving the unique minimizer. The target bound follows from the triangle inequality. The pullback expression is a norm because \(S_a\) is inverse to \(A_a\) on \(\mathcal C_a\). Finally, \(\|T_{ba}c\|_b=\|S_bA_bS_ac\|=\|S_ac\|\). The same identity for a transition, and the bijection between \(X\) and \(\mathcal C_a\), establish equality of Lipschitz constants. \(\square\)
This reconciliation is ordinary weighted least squares in a common semantic space. It does not prove that consensus is correct, reduce statistical error under arbitrary dependence, or turn deterministic recodings into independent evidence. The output must also be reanalyzed into a realizable coefficient vector. An independently fitted redundant output can have an additional inconsistency residual \((I-A_aS_a)G_aA_ax\) before reconciliation.
For Euclidean coefficient errors, conditioning remains important. If \(A_a\) has smallest singular value at least \(m_a>0\) and largest singular value at most \(M_a\), its canonical synthesis \(A_a^\dagger\) has norm at most \(1/m_a\), and \(\|T_{ba}\|\leq M_b/m_a\). An arbitrary left inverse can behave much worse away from the consistent coefficient space, so the canonical-synthesis qualification matters.
A learned atlas can therefore be parameterized within prescribed singular-value bounds, with synthesis tied to its inverse or pseudoinverse. This gives a construction contract. It is not a claim that every unconstrained training procedure respects those bounds.
When atlases change over time but represent the same state, using the pullback norms avoids multiplying arbitrary chart condition numbers at every step. Exact transport telescopes to the direct transport. Numerical synthesis errors and genuinely inconsistent predictions remain defects and must still be bounded.
23. Refinement cannot manufacture missing information
An atlas change is an invertible recoding of the same represented information. Predictive refinement usually adds information. Treating these as identical operations would undo the central distinction of the foundation.
Proposition 22. The upward-lift obstruction. Let \(\gamma:X\to Z\) be a coarse representation and \(\eta:X\to W\) a proposed refined representation. There is a deterministic map \(J:\gamma(X)\to\eta(X)\) satisfying \(\eta=J\gamma\) if and only if
For linear representations, this is equivalent to \(\ker\gamma\subseteq\ker\eta\). If two states share a coarse representation but their refined representations differ, any predicted refined value \(w\) has worst-case error at least
Proof. A lift must be constant on each coarse fiber. Conversely, the stated condition defines it unambiguously on every fiber. The linear version is the kernel factorization criterion. The lower bound follows from the triangle inequality between the two refined values through \(w\). \(\square\)
For a strict refinement, the coarse state is recoverable from the refined one but not conversely. In the linear case this typically means \(\ker\eta\subsetneq\ker\gamma\), exactly the opposite of the condition required for an upward lift.
For example, \(\gamma(x_1,x_2)=x_1\) and \(\eta(x_1,x_2)=(x_1,x_2)\). The states \((0,1)\) and \((0,-1)\) have the same coarse state, but any common point initialization has worst-case error at least one in Euclidean norm. Padding a vector with zeros is a choice, not information recovery.
A refinement step therefore needs one of three explicitly distinguished inputs: an archive that already contains the relevant observations, new observations, or a set-valued or probabilistic state retaining uncertainty about the missing details. An exact archive of measured outputs can initialize a longer history when those samples were stored. It cannot recover a latent coordinate that never affected any measurement.
For adaptive prediction, an attempted lift \(J_t\) contributes its own initialization defect
It belongs in Proposition 18's initial error or in the defect of the composed transition. It cannot be set to zero merely because a finer packet tree has been allocated.
A data-dependent representation schedule also requires care. The pathwise bound applies to the realized schedule when all of its hypotheses hold along that schedule. Comparing two different schedules chosen by true and predicted histories introduces a selection discrepancy. No uniform statistical consistency of adaptive selection follows solely from the deterministic propagation law.
24. A growing operator grammar with coherent evidence
A large algebraic signal model should be able to add operator programs without changing the meaning of its existing evidence. This requires more than assigning a new softmax weight when a candidate is introduced.
Let \(\mathcal G\) be a countable declared program set, with one fixed probability mass function \(\pi(g)>0\) satisfying \(\sum_g\pi(g)=1\). Every program specifies a likelihood \(L_g(D)\) for the same observed data event or density relative to a common base measure. Assume
Let finite active sets \(\mathcal G_n\) increase to \(\mathcal G\). Their unnormalized evidence is \(Z_n=\sum_{g\in\mathcal G_n}\pi(g)L_g(D)\). Only stages with \(Z_n>0\) are considered.
Proposition 23. Exact accounting for an enumerated grammar. Define the full posterior \(P(g\mid D)=\pi(g)L_g(D)/Z\) and extend the active posterior \(P_n(g\mid D)=\pi(g)L_g(D)/Z_n\) by zero outside \(\mathcal G_n\). Then
For any common family of candidate predictive probability kernels, the corresponding predictive mixtures differ in total variation by at most \(\tau_n\). If the omitted unnormalized evidence has a bound
then
Proof. On the active set, \(P=(Z_n/Z)P_n\), so the total mass of the difference there is \(\tau_n\). The mass outside is also \(\tau_n\). Half the sum of absolute differences is therefore \(\tau_n\). Applying the same probability kernel to two mixing measures cannot increase total variation: integrate any event probability, a function taking values in \([0,1]\), against their difference. For the final bound, write \(Z=Z_n+W_n\) with \(0\leq W_n\leq b_n\) and use that \(W_n/(Z_n+W_n)\) increases with \(W_n\). \(\square\)
Since the active sets exhaust the grammar and the evidence sum is finite, \(Z_n\uparrow Z\) and \(\tau_n\downarrow0\). This is posterior approximation inside the declared grammar. It does not establish that the real signal-generating mechanism belongs to it.
For a discrete data sequence, a probability likelihood is at most one, so the omitted prior mass is a valid \(b_n\). A continuous-data density need not be at most one. It requires a density envelope or another justified bound on omitted evidence. A small unnormalized tail can also carry substantial posterior mass when the active evidence is very small.
The fixed prior need not be uniform. For an explicit enumeration \(g_1,g_2,\ldots\), \(\pi(g_k)=2^{-k}\) is a normalized example, with omitted prior mass \(2^{-n}\) after \(n\) programs. It expresses a particular enumeration bias, not a canonical prior on mathematical meaning.
A program outside the enumerated prefix can be accounted for only when its prior and likelihood are defined or bounded. A genuinely outside-grammar component needs its own model. Assigning an unexplained residual probability is not a Bayesian treatment of missing mechanisms.
24.1 Refining a program without double counting it
A representation-rich language can contain many syntactic versions of one mechanism. It can also split one broad mechanism into more specific children. These operations must conserve prior mass when they are intended as a re-expression of the same model.
Proposition 24. Prior-preserving program refinement. Suppose a parent program \(g\) is replaced by children \(g_j\) with weights \(\alpha_j\geq0\), \(\sum_j\alpha_j=1\), and suppose its likelihood is the mixture
Give child \(j\) prior mass \(\pi(g)\alpha_j\) and leave all other programs unchanged. Then total evidence is unchanged and the sum of the children's posterior masses equals the former parent's posterior mass. Posterior predictive distributions are also unchanged when the same mixture relation holds for the joint likelihood of \((D,y)\) for every predicted event \(y\).
Proof. The children's total unnormalized evidence is \(\sum_j\pi(g)\alpha_jL_{g_j}(D)=\pi(g)L_g(D)\). Thus both the total normalizer and every other program's posterior are unchanged. The same sum for joint evidence gives the predictive statement after normalization by the unchanged evidence of \(D\). \(\square\)
The mixture identity must hold for the whole data record. Averaging one-step likelihoods independently at each time generally describes a switching mechanism, not a fixed child selected once for the whole sequence.
If the children are only exact atlas recodings of the same operator and data model, their likelihoods are identical. Splitting the original prior mass leaves every prediction unchanged. Giving each duplicate a fresh full prior mass instead changes the model's prior preference. That may be an intentional new prior, but it is not invariance under representation.
The implementation's existing held-out predictive weights, described in Section 10, are not replaced by these identities merely by calling them posterior weights. The proposition supplies the probabilistic semantics that a future grammar-weighting implementation would need to enforce.
25. The finite learned-realization theorem
The preceding results combine into a concrete completion of one part of the learning problem. This is the integrated identification claim for the autonomous linear sector.
Proposition 25. Observationally grounded finite linear realization. Assume all of the following:
- The data come from an autonomous finite-dimensional linear process with declared outputs \(y_t=Cx_t\).
- A delay length \(h\) satisfies \(\operatorname{rank}O_{h+1}=\operatorname{rank}O_h\).
- Exact observed input windows in \(Z_-\) span \(V_h=\operatorname{im}O_h\), and \(Z_+\) records their true one-step successor windows.
Then an algorithm using only \((Z_-,Z_+)\) recovers the unique transition on \(V_h\). Initialized with a fully observed length-\(h\) window, it predicts every later output exactly. Any lossless atlas of \(V_h\), any complete packet decomposition of that space, and every coherently transported realization of the learned transition preserve the same predictions. Hidden variables outside the observable quotient need not be recovered.
If the fixed-space assumptions of Proposition 16 replace exact observations, that proposition bounds operator error. Whenever the resulting transition defects and Lipschitz constants satisfy Proposition 18, its product-sum bound controls the output prediction error. A common contraction margin yields the stated uniform bound.
Proof. Proposition 14 gives the unique true transition \(D_h\) on \(V_h\). Because the observed windows span \(V_h\), Proposition 15 identifies exactly its restriction, which is the entire required operator. Equality of the learned and true states follows by induction after initialization. The current output is the last block of the history state. Propositions 1--4 preserve state reconstruction, transport and complete decomposition, so they preserve that output sequence. The perturbation statement is the direct composition of Propositions 16 and 18 with the declared output decoder. \(\square\)
This establishes an actual route from observations to a composable algebraic signal model under explicit finite assumptions. It is more than a fit on a few points. It is also less than unrestricted nonlinear world-model discovery. In particular, the algorithm does not infer assumptions 1--3 merely by returning a matrix. The accompanying certificate keeps the data-verified identities separate from those model assumptions.
26. Connection to controlled prediction
Proposition 25 completes the autonomous observable linear sector under its declared coverage assumptions. It does not establish a common delay operator for arbitrary actions. Part III constructs controlled predictive states from interventional experiment probabilities. The two results are compatible routes, not interchangeable learning assumptions.
\newpage
Part III. Controlled stochastic realization
The next step is not to rename a deterministic predictor as a probabilistic one. A stochastic signal model must say what its state predicts, which actions intervene on the process, how observations update the state, and why every resulting number is a probability. These are separate obligations. A low dimensional linear representation can satisfy its reconstruction equations while failing probability positivity.
Part II established observable learning for an autonomous linear signal. Here the observable objects are probabilities of finite action and observation experiments. We derive a finite realization from them, show when it is minimal, give an explicit probability-valid state domain, and separate exact reconstruction from noisy estimation. The final construction connects back to Komorebi's packets through a normalized closure-defect identity.
The classical antecedents are predictive state representations, observable operator models and finite-rank realization [R5, R7, R9]. Those traditions already use probabilities of tests and products of operators. The purpose here is to make their assumptions and consequences explicit within the realization calculus. The twelve propositions below are a connected development, not twelve claims of historical novelty. No statement here replaces a proof with a numerical check.
27. Actions and observations are different inputs
Let \(\mathcal A\) and \(\mathcal O\) be finite nonempty action and observation sets. A symbol \(s=(a,o)\) records an action followed by an observation. A word \(w=s_1\cdots s_k\) has temporal order from left to right. The empty word is \(\epsilon\).
Write
The process is reset to one specified initial distribution before this experiment. The actions in this definition are fixed externally. This is not, in general, the distribution obtained by conditioning a passive record on an adaptively selected sequence of actions. We assume the consistency laws
These laws specify causal probabilities on finite observation trees. A fixed history-dependent policy can be evaluated by selecting the action on each branch and multiplying the appropriate conditional probabilities. This does not establish optimality of that policy or its identification from insufficient data.
27.1 A finite hidden-state construction
A concrete realization has a belief \(p\in\mathbb R^n_+\) with \(\mathbf1^\top p=1\). For each \(s=(a,o)\), let \(M_s\in\mathbb R^{n\times n}\) be entrywise nonnegative, with
The matrix entry \(M_{a,o}(j,i)\) is the probability of next hidden state \(j\) and observation \(o\), given current hidden state \(i\) and chosen action \(a\). This is a discrete probability instrument. It includes an action-specific transition followed by an observation, but does not require that factorization.
Put \(M_w=M_{s_k}\cdots M_{s_1}\). Then
on branches of positive probability. An observation with zero probability has no conditional posterior under that model. It must not be handled by silently dividing by an arbitrary numerical floor.
Proposition 26. Finite instrument closure. Starting with \(W_0=\operatorname{span}\{\mathbf1^\top\}\), define
This stabilizes after at most \(n-1\) strict enlargements at
Choose a full row rank matrix \(R\) with row space \(W_*\) and first row \(\mathbf1^\top\). There are unique matrices \(B_s\) with
For \(z=Rp\), both the next-observation likelihood and the posterior predictive state close:
Every linear summary determining all word probabilities by linear readouts for every initial belief has row space containing \(W_*\).
Proof. The nested-space argument from Proposition 8 applies to the finite family \(M_s\). Multiplication of a test row by \(M_s\) prepends the symbol \(s\) to that test, so the stable span is exactly the displayed word span. Invariance supplies the matrices \(B_s\), and full row rank makes them unique. Apply \(R\) to the unnormalized and normalized hidden-state updates. If all test probabilities are linear functions of another summary, then each row \(\mathbf1^\top M_w\) lies in its row space. Equality on all beliefs implies equality on their linear span, which is \(\mathbb R^n\). \(\square\)
This proposition uses a known hidden-state model. It is not yet an observable-data learner. It also concerns all initial beliefs, whereas a model initialized at one fixed \(p_0\) can have a smaller reachable predictive dimension. The next construction handles that distinction without requiring recovery of hidden states.
28. The state is a function of possible future experiments
For a history \(h\), define its unnormalized future-behavior column
where \(t\) ranges over every finite test word. Form the infinite history/test matrix
Its rank means the dimension of the span of its columns, equivalently the supremum of the ranks of finite submatrices. The columns are functions, not vectors obtained by physically observing every future. Their finite coordinate realization becomes learnable only when enough experiment probabilities are available.
Proposition 27. Minimal finite behavior realization. If \(\operatorname{rank}\mathscr H=r<\infty\), the vector space \(V=\operatorname{span}\{v_h\}\) has dimension \(r\). There are well-defined linear maps \(B_s:V\to V\), an initial vector \(v_\epsilon\), and a linear mass functional \(\lambda(v)=v(\epsilon)\) satisfying
Every linear realization of \(f\) has dimension at least \(r\). Any two reachable and observable realizations of dimension \(r\) differ by an invertible change of coordinates. For \(f(h)>0\), the normalized state \(v_h/f(h)\) determines every future test probability conditional on \(h\) and the declared future interventions.
Proof. If \(\sum_hc_hv_h=0\), then for every test \(t\),
because \(st\) is another test. Thus shifting columns preserves all linear relations and defines \(B_s\). The probability formula follows by induction. For any \(d\)-dimensional realization, the history/test matrix factors into test rows and reachable state columns, so its rank is at most \(d\).
For uniqueness, map the reachable vector after each history in one minimal realization to the corresponding vector in the other. If a linear combination vanishes in the first, every future test evaluates to zero on it. Those same evaluations vanish in the second, where observability makes the combination zero. The map is therefore well-defined and bijective, and intertwines every symbol operator and both boundary vectors. Finally \(f(ht)/f(h)\) is exactly the conditional experiment probability. \(\square\)
The coordinate dimension \(r\) includes probability mass. Normalized states lie in the affine hyperplane \(\lambda z=1\), so their affine dimension is at most \(r-1\). A finite hidden-state model of size \(n\) has behavior rank at most \(n\). The theorem itself requires finite behavior rank, not a preferred hidden-state ontology.
Two realized histories with equal normalized future columns are equivalent for every declared future test. This is a predictive equivalence, not equality of their physical pasts or lossless archives.
29. Recovering the operators from an observable probability table
Select \(r\) tests \(t_1,\ldots,t_r\) and \(r\) histories \(h_1,\ldots,h_r\), with \(t_1=\epsilon\), for which the square matrix
is invertible. Such a choice exists under the rank hypothesis of Proposition 27. The nonzero empty-test row can be extended to a row basis. Select independent history columns for that basis.
For each symbol \(s\), also obtain the shifted table and initial test vector
Proposition 28. Observable table reconstruction. Assume the full behavior rank is at most \(r\) and the displayed \(r\times r\) matrix \(H\) is invertible. Then the full rank is exactly \(r\), and
is a minimal realization of every word probability. Its state after history \(h\) is
The normalized state is \(b(h)/f(h)\) when \(f(h)>0\). Its update uses the formula in Proposition 26, with no hidden-state coordinates supplied to the reconstruction algorithm.
Proof. Invertibility gives rank at least \(r\), so the assumed upper bound is attained. Evaluation on the selected tests is an isomorphism \(V\to\mathbb R^r\), because it maps the chosen basis of history columns to the invertible matrix \(H\). The shifted images of those basis columns are precisely the columns of \(G_s\), so \(B_sH=G_s\). The first selected test evaluates the mass. Initial evaluation gives \(b_0\). Proposition 27 supplies correctness for every history, not only the selected columns. \(\square\)
There are two ways to justify coverage in this statement. One can know that the selected row and column families span the true behavior space. Alternatively, a justified upper bound \(r\) on the process rank, together with the invertible \(r\times r\) block, proves that coverage. The observed invertibility by itself only gives a lower bound. An unknown higher-rank process can agree on a finite table and differ later.
Exact experiment probabilities are population quantities. The reference example below supplies rational tables from a declared process and keeps that table generator separate from the learner. Section 36 supplies one finite sampling protocol for estimating those entries. Neither a hidden matrix nor its probability cone is secretly passed to the observable learner.
30. Probability validity is not a rank condition
A real matrix product can yield a negative number. Therefore a fitted factorization is not a stochastic model merely because it is low rank or exactly interpolates a probability table. Spectral estimators can encounter negative-probability problems in practice [R10]. We address the structural issue directly.
In the canonical behavior space \(V\), define
Proposition 29. Canonical probability cone. The set \(\mathcal C\) is a closed convex pointed cone that spans \(V\). Every \(B_s\) preserves it. The mass functional is strictly positive on its nonzero elements, and
Consequently the normalized base
is a compact convex probability-valid state space. Every positive-probability branch maps \(\mathcal K\) into itself by normalized updating. In the core-test coordinates of Proposition 28, \(\mathcal K\subseteq\{1\}\times[0,1]^{r-1}\).
Proof. Each test evaluation is a continuous linear functional on finite-dimensional \(V\), so the cone is closed and convex. If both \(v\) and \(-v\) belong to it, every test value is zero and \(v=0\). It spans \(V\) because it contains every column \(v_h\). Also \((B_sv)(t)=v(st)\), proving invariance.
The consistency law \(\sum_ov(t(a,o))=v(t)\) holds for each history column and hence throughout \(V\). Iterating it along any fixed action schedule expresses \(v(\epsilon)\) as a sum of nonnegative values for all observation words on that schedule. If \(v(\epsilon)=0\), each summand is zero, for every schedule and length, so \(v=0\). The same sums show \(0\leq v(t)\leq1\) when \(\lambda v=1\). They also give mass conservation. The base is closed and bounded in the finite core coordinates and thus compact. Invariance and positive normalization prove the posterior assertion. \(\square\)
This cone is defined by all finite tests. Its existence is not an algorithm for finding a finite list of inequalities. A finite list gives an outer approximation and can miss a longer-test violation. The cone also need not equal the convex hull of the states reachable from the fixed initial distribution.
A practical sufficient certificate can instead use a known closed cone \(C_0\) with nonempty normalized base, a strictly positive mass functional, \(b_0\in C_0\), \(B_sC_0\subseteq C_0\), and mass conservation. A nonnegative hidden-state realization supplies one after a coordinate change. Such a certificate may use information unavailable to an observable-only learner and must be labeled accordingly.
Negative entries of \(B_s\) do not, by themselves, prove invalidity. Coordinates can be signed while the represented state remains in the correct cone. Conversely, clipping negative coefficients can destroy the operator identities. Validity belongs to the pair consisting of the operators and their admissible cone, not to a coordinatewise heuristic.
31. Packets can reconstruct a belief without being beliefs
Let \(P_1,\ldots,P_m\) be a complete family of linear packet maps on a predictive space, with \(\sum_iP_i=I\). Algebraic reconstruction gives \(z=\sum_iP_iz\). It does not imply that \(P_iz\) is an unnormalized probability state.
Proposition 30. Positive packet decomposition and its boundary. Suppose a cone \(C\) has a strictly positive mass functional \(\lambda\). If every packet map preserves \(C\), then each normalized \(z\in C\) decomposes as a genuine mixture
A zero weight has \(P_iz=0\). In the standard cone \(\mathbb R^n_+\), however, a complete family of cone-preserving idempotents must consist of diagonal coordinate-selection matrices. In particular, a nontrivial Fourier or Haar packet family cannot generally be reinterpreted as a mixture of probability states merely because it is orthogonal and complete.
Proof. Cone preservation and strict positivity give the weights and the zero-weight claim. Apply \(\lambda\) to the complete decomposition to obtain their sum. Normalization of each nonzero term gives the mixture.
On \(\mathbb R^n_+\), cone preservation is entrywise nonnegativity. Since the matrices sum to \(I\), every off-diagonal entry in every matrix is zero. Idempotence forces each diagonal entry to be zero or one, and completeness makes the selections disjoint and exhaustive. \(\square\)
For a two-state belief \(p=(3/4,1/4)^\top\), the Haar detail projector
produces \((1/4,-1/4)^\top\). It is a valid signed detail packet, but not a subprobability distribution. Its total mass is zero despite being nonzero, which cannot happen inside a cone with a strictly positive mass functional.
This does not prohibit useful spectral packets. It fixes their type. Signed algebraic packets can be used for transfer, approximation and diagnostics. They must be synthesized before probability interpretation unless a separate positive decomposition has been established. That distinction prevents a reconstruction identity from being misused as a Bayesian mixture identity.
32. Action coverage is a mathematical requirement
Proposition 31. Missing-intervention nonidentifiability. There are finite controlled stochastic models that agree on every possible observation record generated using only action \(a\), but disagree under another action \(b\). Therefore even infinitely many passive records confined to \(a\) cannot identify the unrestricted controlled model.
Proof. Use one hidden state and two observations. In both models, action \(a\) gives probabilities \((1/2,1/2)\). Under \(b\), let the first model give \((1,0)\) and the second \((0,1)\). Every length-\(k\) record under \(a\) has probability \(2^{-k}\) in both models. Their first observation under \(b\) differs with probability one. \(\square\)
This is a coverage obstruction, not an optimization problem. It remains even when the hidden dimension is known exactly. Richer transforms of the same action-\(a\) observations do not create data about \(b\).
Action selection based on past observations requires a further distinction. Under a known policy \(\pi\), an action/observation word has joint probability
provided the action randomization has no access to hidden information beyond the declared history. The policy factors depend on the observed branch. Conditioning on a selected action sequence generally reweights branches rather than removing these factors. Recovering intervention probabilities from such data requires an appropriate coverage and reweighting argument. The reset experiments in Section 36 avoid that additional inference problem.
33. Learning the coordinate space under bounded noise
Part II bounded operator error after a state space and basis were fixed. Here we also bound one step of data-dependent space selection. This remains a finite-matrix statement. It does not certify that a chosen finite history/test table captures the full process.
Let \(X\in\mathbb R^{m\times n}\) have rank \(r\geq1\), let \(s=\sigma_r(X)>0\), and observe \(Y=X+E\) with \(\|E\|_2\leq\varepsilon<s/2\). Let \(P\) project onto \(\operatorname{im}X\), and let \(\widehat P\) project onto the leading \(r\) left singular vectors of \(Y\).
Proposition 32. Finite-table rank separation and subspace stability. One has
The \((r+1)\)st singular value is interpreted as zero when \(r\) is the maximum possible rank. If \(U\) and \(\widehat U\) are orthonormal bases of the two ranges, then
Thus \(x\mapsto\widehat U^\top x\) is injective on the true column space. Its inverse on that space is
with norm at most \([1-(2\varepsilon/s)^2]^{-1/2}\).
Proof. The min-max characterization of singular values and the triangle inequality imply that adding \(E\) changes each singular value by at most \(\varepsilon\). The leading-\(r\) singular subspace gives \(\|(I-\widehat P)Y\|_2=\sigma_{r+1}(Y)\leq\varepsilon\). Hence
For equal-rank orthogonal projections, the latter norm equals \(\|P-\widehat P\|_2\). To see the equality, diagonalize \(U^\top\widehat U\) by singular value decomposition. Each nontrivial paired principal plane contributes a sine to both norms, while common directions contribute zero. On the same planes the singular values of \(\widehat U^\top U\) are the corresponding cosines. This proves the last bound and the inverse formula. \(\square\)
A threshold \(\tau\) with \(\varepsilon<\tau<s-\varepsilon\) recovers the rank of this finite matrix. The unknown true gap must be justified or replaced by a valid confidence interval. Rank thresholds do not exclude modes absent from the selected tests or unexcited by the selected histories.
The inverse formula contains the unknown true basis \(U\). It proves identifiability and conditioning relative to that space, not that a learner already knows a perfect decoder. When using the estimated subspace as coordinates for a realization, all boundary vectors, symbol operators and probability constraints must be transformed together. Signed singular-vector coordinates are not new independent observations.
34. Noisy operator reconstruction in a fixed core
After choosing a core block, the square reconstruction formula yields a particularly direct error identity. It is useful even when the core was chosen on a separate discovery data set.
Let \(G_s=B_sH\) with \(s_H=\sigma_{\min}(H)>0\). Observe
Proposition 33. Noisy core-block realization. The estimate \(\widehat B_s=\widetilde G_s\widetilde H^{-1}\) satisfies
and therefore
A bound using only the fitted operator and the same perturbation budgets is
If \(\widetilde s=\sigma_{\min}(\widetilde H)\) is known and \(\|E\|_2\leq e<\widetilde s\), the last denominator may be replaced by \(\widetilde s-e\).
Proof. Since \(\widetilde H\) is invertible, subtract \(B_s\) from its fitted expression and collect the terms. The least singular value of \(\widetilde H\) is at least \(s_H-\|E\|_2\). For the second identity, use
and multiply by \(H^{-1}\). Finally \(s_H\geq\widetilde s-e\) by singular-value perturbation. \(\square\)
When the first test is empty, the mass row \(\lambda=e_1^\top\) is fixed rather than estimated. Set the first component of the initial vector to one as required by \(f(\epsilon)=1\). Errors in other initial test probabilities remain initial-state error. This convention removes an avoidable normalizer-estimation ambiguity but does not enforce cone preservation or mass conservation for the fitted symbol family.
Numerical implementation should solve the linear systems rather than explicitly form a floating-point inverse. The theorem describes the exact algebraic estimator. Finite-precision solve error is an additional residual, which the reference script measures separately.
35. Conditioning a prediction can amplify its error
Let \(C\) be a probability cone with compact normalized base \(K\) and mass row \(\lambda\). Work in a declared norm and its dual. Suppose true normalized states have norm at most \(R\). Let a true branch have
For a predicted state \(\widehat z\), form \(\widehat u=\widehat B_s\widehat z\), \(\widehat p=\lambda\widehat u\), and define \(d=\widehat u-u\). Write \(e=\|\widehat z-z\|\), \(\|\widehat B_s\|\leq L\), and \(\|\widehat B_s-B_s\|\leq\epsilon_B\). Then
Proposition 34. Normalization-aware filtering stability. If \(\|\lambda\|_*D<\eta\), then \(\widehat p>0\) and
so
Alternatively, if each realized fitted denominator is checked to satisfy \(\widehat p_t\geq\widehat\eta_t>0\), the pathwise recursion obeys
The product-sum bound of Proposition 18 applies. A uniform bound follows if these constants have a common contraction margin and bounded forcing. The algebraic bound alone does not certify that a fitted state or every unchosen branch is probability-valid.
Proof. The branch likelihood differs by \(\lambda d\), giving the lower bound on \(\widehat p\). Subtract \(u/p\) from \((u+d)/(p+\lambda d)\) and use \(u=pz^+\). The induced norm of \(I-z^+\lambda\) is at most \(1+R\|\lambda\|_*\). Applying that identity at each step gives the recursion. \(\square\)
Rare evidence can create a large posterior change from a small absolute error. For example, take a true unnormalized vector \((\eta,0)\) and a perturbed vector \((\eta,\eta)\). Their \(\ell^1\) difference is \(\eta\), while their normalized vectors differ by one in \(\ell^1\). Both arise as possible branch outputs of valid probability instruments when \(0<\eta\leq1/2\), after assigning the remaining mass to another branch. There is no uniform posterior Lipschitz constant independent of a probability margin.
When a known closed convex feasible base is available, Euclidean projection of a tentative posterior back onto that base cannot increase its distance from the true posterior. This is a declared repair to the estimator, not an identity for an unmodified update. It also does not repair negative branch probabilities of arbitrary fitted operators before normalization.
36. A finite experimental route to the probability table
The following deliberately simple protocol makes a finite-sample statement possible without confusing a deterministic identity with a statistical guarantee.
Fix positive integers \(M\) and \(N\) and choose \(M\) finite words before estimation. For each word, reset independently to the same initial distribution, execute its prescribed action sequence, and record the indicator that all observations match. Obtain \(N\) independent trials per word. Dependencies between different words' estimates are allowed, but the \(N\) trials used for each individual mean must be independent with the declared success probability.
Proposition 35. Reset-experiment probability certificate. Let \(\widehat f(w)\) be the sample mean for each selected word. For \(0<\delta<1\), with probability at least \(1-\delta\), all \(M\) estimates satisfy
If the selected words cover every entry of \(H\), every \(G_s\), and \(b_0\), then on the same event
These bounds feed Propositions 32 and 33 whenever their respective gap conditions hold. The empty-word probability may be set exactly to one instead of sampled.
Proof. For a centered Bernoulli variable \(X-p\), write \(g(t)=\log\mathbb E\exp(t(X-p))\). Then \(g(0)=g'(0)=0\), and \(g''(t)\) is the variance of a Bernoulli under exponential tilting, at most \(1/4\). Integrating twice gives \(g(t)\leq t^2/8\). Independence and the exponential Markov inequality give an upper-tail bound \(\exp(-2N\epsilon^2)\) after optimizing \(t=4\epsilon\). The lower tail has the same bound. A union bound over the \(M\) words proves the simultaneous result. The spectral norm is bounded by the Frobenius norm, yielding the matrix and vector estimates. \(\square\)
This is the Bernoulli form of Hoeffding's inequality [R11], proved here for the protocol used. It is a conservative entrywise-to-spectral bound, not a claim of optimal sample complexity. The number and lengths of chosen experiments are part of its cost.
If core tests and histories were selected from data, use a separate independent discovery sample and freeze the selected finite experiment list before these estimation trials. Conditional on that selection, the proposition applies. Alternatively, a uniform bound over a fixed larger candidate list can support data-dependent selection within that list. Reusing a post-selected table without either argument changes the claim.
A single dependent trajectory, an adaptive stopping rule, unknown resets, drifting sensors, or experiments chosen without coverage do not satisfy this proposition as stated. Those settings require their own concentration and identifiability arguments. Synthetic reset sampling in the reference example demonstrates the protocol, not its validity for every measurement source.
37. A filtering defect with a packet interpretation
The relation to the realization calculus is particularly direct before normalization. Let \(R\) be a linear summary of a hidden probability space, let \(\lambda R=\mathbf1^\top\), and let \(B_s\) be a proposed reduced branch operator. For a true belief \(p\), put
The defect \(\delta_s\) is a vector in the reduced representation, not a probability of model failure. It records the failure of the unnormalized instrument to commute with the summary.
Proposition 36. Exact normalized closure defect. Suppose the true and proposed branch likelihoods
are positive, and let \(\widehat z^+=B_sz/b\). Then
If \(M_s\neq0\) and both likelihoods and normalized updates agree on every strictly positive belief, then \(RM_s=B_sR\). Agreement of normalized updates alone need not force this equality.
Proof. The true unnormalized summary is \(B_sz+\delta_s\), with mass \(a=b+\lambda\delta_s\). Subtract the two normalized vectors and collect terms. Since a nonzero nonnegative branch has positive likelihood on every strictly positive belief, agreement permits multiplication by the common mass. This gives equality of unnormalized summaries on the interior of the probability simplex, whose linear span is the whole hidden vector space. For the last claim, a one-state model can assign different positive branch probabilities while its normalized state always remains the scalar one. \(\square\)
For an orthogonal retained projector \(P\) preserving the mass row, put \(Q=I-P\), \(R=P\) on its image and \(B_s=PM_sP\). Then
This is the unresolved-to-retained transfer block from Part I. Its posterior effect is not simply its norm: normalization applies the oblique map \(I-\widehat z^+\lambda\) and divides by the true branch probability. A component proportional to \(\widehat z^+\) changes only the branch mass and disappears from the normalized defect, while still changing likelihood-based inference over mechanisms.
A certificate should therefore retain both \(\lambda\delta_s\) and the posterior defect. A posterior that looks correct can hide a wrong evidence model. This is the stochastic counterpart of the distinction between reconstruction and predictive closure.
38. A complete two-action rational example
Use two hidden states only to generate a reference process. The learner does not receive these matrices. Let
The observation matrices are
The two transitions do not commute. Thus the process is genuinely action-dependent rather than one autonomous recurrence with an action label attached.
Choose tests \((\epsilon,(0,1))\) and histories \((\epsilon,(0,1))\). The observable tables give
The selected coordinate \(q\) is the probability of seeing observation one after choosing action zero. In the hidden realization,
The observable reconstruction \(B_s=G_sH^{-1}\) yields
The coefficients include negatives. Nevertheless, every operator preserves the cone \(R\mathbb R^2_+\), because \(B_sR=RM_s\). A sufficient normalized base is
This cone certificate uses the known generating model and is checked separately. It is not inferred merely from the finite probability table, and \(K_0\) is not claimed to be the maximal canonical base of Proposition 29.
The likelihoods and posterior maps can now be read from the matrices. For action zero and observation one,
Starting at \(q_0=13/24\), this gives probability \(13/24\) and \(q_1=15/26\). If action one then produces observation zero, its probability is
and the new predictive coordinate is
The joint probability of this two-step action/observation word is \(365/1152\). The reference checks compare these values and every word up to a declared finite test length against the hidden generator in exact rational arithmetic. The proposition, under its rank and coverage assumptions, is what extends the result to all word lengths.
The model's online state is not a sequence of newly measured future probabilities. It is a recursively maintained prediction initialized from learned experiment probabilities and updated after each actual action and observation. Uncertainty in those learned parameters is separate from the stochastic uncertainty in the next observation.
39. Integrated controlled stochastic realization
Proposition 37. Conditional observation-based stochastic completion. Suppose a causal finite-alphabet signal process has behavior rank at most \(r\), and a selected \(r\times r\) history/test block containing the empty test is invertible. Exact values of that block, every one-symbol shifted block, and the initial core-test vector determine a minimal realization of every future experiment probability. Its normalized state is recursively accessible after each positive-probability history, without reconstructing the hidden physical state.
The exact realization has the canonical probability cone of Proposition 29. Under any invertible atlas, transform the state, symbol operators, mass row and cone together:
Every likelihood and posterior prediction is then unchanged. A complete packet decomposition preserves these predictions after full synthesis. Packetwise probability interpretation requires the extra positivity condition of Proposition 30.
For estimates obtained under Proposition 35's reset protocol, the simultaneous event there bounds the finite table perturbations. A separated finite-table spectrum yields Proposition 32's subspace certificate. A well-conditioned fixed core yields Proposition 33's operator certificate. On branches satisfying Proposition 34's likelihood margins, the resulting normalized state errors obey its recursive bound. These perturbation statements alone do not certify a learned cone or global positivity of an estimated model.
Proof. Proposition 28 reconstructs the unique minimal behavior realization, and Proposition 27 gives its conditional interpretation. Proposition 29 supplies probability validity for the exact realization. Substituting the transformed quantities cancels every \(T^{-1}T\), including in the scalar normalization. Full packet synthesis is the identity, whereas Proposition 30 characterizes the additional probability-mixture condition. The perturbation conclusions are the stated chain of propositions with their separate assumptions retained. No missing gap, coverage or positivity assumption follows from matrix inversion alone. \(\square\)
This theorem closes a finite controlled and stochastic sector of the foundational model. It does not claim unrestricted discovery of a nonlinear environment from a single trajectory. It does identify what a learned predictive state represents, why its operations compose, where its evidence comes from, and which defects change likelihoods or posteriors.
40. From realization to admissible domains
Proposition 37 identifies a finite controlled realization from a covered population table. Existence of its canonical probability cone does not provide a finite certificate for a fitted model. Part IV supplies such certificates in finite positive-lift sectors, while keeping internal validity distinct from statistical coverage.
\newpage
Part IV. Certified probability domains and robust learning
A matrix learned from observations is not yet a probability model. Part III made that distinction exact: a finite predictive realization needs an admissible cone preserved by its action and observation operators. It supplied a canonical cone defined by all possible future tests, but did not supply a finite way to verify that cone for an estimated model.
This block develops a constructive finite sector of that missing step. It gives a certificate that can be checked by rational arithmetic, a search procedure that does not need hidden-state matrices, and a robust version that covers an entire declared uncertainty set. It also proves why a finite polyhedral certificate is not available for every valid finite-dimensional process. The resulting distinction is between an exact finite certificate, an uncertified candidate, and an explicitly modified approximation. Those are different mathematical outcomes.
The positive-realization and invariant-cone ingredients are classical. Jaeger's account includes the polyhedral-cone characterization and an irrational-rotation obstruction [R12]. Polyhedral invariance is also a standard subject of linear systems and optimization [R13]. The present contribution is a self-contained connection of these ingredients to Komorebi's observable learner, probability packets, confidence sets and approximation defects. The controlled rational rotation example is an adaptation of that classical obstruction, not a claim of discovering the phenomenon.
The next twelve propositions continue the numbering from the preceding volume. All have proofs below. The exact computations check stated finite instances and counterexamples. They are not a Lean formalization and do not change the production PR-LASM 2.3.0 runtime.
41. The object to be certified
Fix a real predictive space \(E=\mathbb R^r\), a nonzero row \(\lambda\), and a finite family \(B_s\), indexed by \(s=(a,o)\). The action and observation alphabets are finite and nonempty. Every declared observation branch is present, including branches represented by the zero operator. Temporal words act from right to left:
The mass identity is an exact constraint,
An initial state has \(\lambda z_0=1\). The proposed probability of a word is \(\lambda B_wz_0\). An invertible coordinate change must transform \(z_0\), \(B_s\), \(\lambda\) and the admissible cone together, as in Proposition 37.
A probability-domain certificate specifies a closed, pointed, generating cone \(C\subset E\), with \(\lambda\) strictly positive on \(C\setminus\{0\}\), such that
Its normalized base is \(K=\{z\in C:\lambda z=1\}\). These conditions imply nonnegative word probabilities and their per-action normalization at every depth. Positive-probability updates remain in \(K\). Zero-probability branches have no conditional posterior and must not be assigned one by denominator flooring.
There are three separate claims. Internal validity says the fitted model is a stochastic model. Data compatibility says it agrees with the observations within a specified error budget. Coverage of the generating process says the assumptions and uncertainty set actually include the process that generated those observations. The first can have a deterministic certificate even when the other two remain conditional.
42. Finite tests give an outer hierarchy, not automatic closure
Define the test cones
Each \(C_h\) is polyhedral because it is an intersection of finitely many homogeneous half-spaces. It need not be a pointed cone or have a bounded normalized base.
Proposition 38. Test-cone hierarchy and a finite stopping certificate. The cones satisfy
Their intersection \(C_\infty\) is the largest cone contained in \(C_0\) and invariant under all \(B_s\). If an exact equality \(C_{h+1}=C_h\) occurs, then \(C_h=C_\infty\) and the finite inequalities certify all future word probabilities. If the representation is observable, meaning that the rows \(\lambda B_w\) span \(E^*\), then \(C_\infty\) is pointed and \(\lambda\) is strictly positive on its nonzero elements. If reachable vectors \(B_wz_0\in C_\infty\) span \(E\), the cone is also generating.
Proof. For a first symbol \(s\) followed by a word \(v\), \(B_{sv}=B_vB_s\). The inequalities for all words of length at most \(h+1\) are therefore exactly the displayed recurrence. If \(z\in C_\infty\), every test on \(B_sz\) is another test on \(z\), so \(C_\infty\) is invariant. Conversely, invariance of a cone contained in \(C_0\) makes every successive test nonnegative, hence places it in \(C_\infty\). Equality at one finite step establishes invariance immediately and therefore gives the intersection.
If both \(z\) and \(-z\) lie in \(C_\infty\), every test row annihilates \(z\), so observability gives \(z=0\). For \(z\in C_\infty\) with \(\lambda z=0\), sum word probabilities along any fixed action sequence. Mass conservation makes the sum zero. Every term is nonnegative, so every word test vanishes and again \(z=0\). Spanning by reachable cone elements gives the final assertion. \(\square\)
A useful negative witness is a state admitted by the current inequalities and a specific longer word with negative proposed probability. It demonstrates that this state needs to be excluded. It does not prove that the initial state, or every reachable state, is invalid.
Finite stabilization is sufficient, not promised. Even three predictive coordinates can require a nonpolyhedral domain, as Section 51 shows. Checking a long list of trajectories is not the same as checking the invariance identity that makes every trajectory valid.
43. A finite certificate with a positive lift
Let \(V=[v_1\ \cdots\ v_m]\in\mathbb R^{r\times m}\) have full row rank and satisfy
Its columns generate \(C=V\mathbb R_+^m\) and its normalized base is \(K=\operatorname{conv}\{v_1,\ldots,v_m\}\). Finite generation makes the cone closed. Its normalized generators and full row rank make it pointed and generating, with a strictly positive mass functional.
Proposition 39. Finite positive-lift certificate. This cone is invariant under \(B_s\) if and only if there is an entrywise nonnegative matrix \(A_s\in\mathbb R^{m\times m}\) such that
If these identities hold for every symbol and \(z_0=V\alpha_0\) for \(\alpha_0\geq0\) with \(\mathbf1^\top\alpha_0=1\), then
Thus the finite witness defines a nonnegative probability instrument whose image under \(V\) is the predictive model. It certifies all word probabilities without enumerating them.
Proof. Invariance places each \(B_sv_j\) in the conical hull of the columns of \(V\). Choose its nonnegative coefficients as column \(j\) of \(A_s\). Conversely, \(B_sV\alpha=VA_s\alpha\in C\) for every \(\alpha\geq0\). Multiplying the intertwining identity by \(\lambda\) and summing over observations gives
Iterating \(B_sV=VA_s\) gives \(B_wV=VA_w\). Applying the mass row and the initial coefficients proves the word identity. Strict positivity of \(\lambda\) follows because a nonzero conical combination has positive total coefficient mass. \(\square\)
The \(m\) nonnegative coordinates need not be physical hidden states, and \(m\) need not equal the minimal predictive dimension \(r\). They are one finite realization witnessing validity. When \(V\) has dependent columns, the lift can be nonunique. Those redundant coordinates do not change the predictive state after synthesis.
There is a dual finite test. If a complete inequality description of the same cone is \(C=\{z:Gz\geq0\}\), then \(GB_sv_j\geq0\) for every generator \(j\) is equivalent to invariance. One may also certify \(GB_s=H_sG\) with \(H_s\geq0\). The latter certificate is sufficient by direct multiplication and necessary by the polyhedral cone form of Farkas' lemma [R13]. The reference calculation uses the direct generator test, not an unreported numerical optimization certificate.
44. Robustness is a quantified invariance statement
A certificate for one fitted operator does not automatically cover neighboring operators. For uncertainty, the quantifier must range over the whole proposed model set, including correlations required by mass conservation.
Let \(\mathcal U\) be the convex hull of finitely many complete operator families \(\{B_s^{(j)}\}_s\). Every family must satisfy the exact per-action mass identities. The convex coefficient is shared across that family's symbols. A product of independent uncertainty sets can instead be handled by its joint vertices.
Proposition 40. Robust finite-domain certificate. For a fixed normalized generator matrix \(V\), suppose every vertex family has witnesses
Then every family in \(\mathcal U\) preserves \(C=V\mathbb R_+^m\) and is mass conserving. If its initial state lies in \(K\), it defines valid finite observation-tree probabilities at every horizon. This remains true when a family in \(\mathcal U\) is selected separately at each history.
For a complete facet description \(C=\{z:Gz\geq0\}\), let \(g_i^\top\) be its rows. A sufficient norm-ball certificate around \(\widehat B_s\) is
It covers all perturbations with \(\|B_s-\widehat B_s\|_2\leq\epsilon_s\) that also satisfy the exact mass constraints.
Proof. For convex weights \(\theta_j\), take \(A_s=\sum_j\theta_jA_s^{(j)}\). The equality and nonnegativity survive convex combination, as does mass conservation. At each history, every selected family maps the current cone state back into the cone. Applying the mass identity at that node proves recursive normalization, regardless of how the family was selected. For the norm ball,
The margin makes every perturbed generator image satisfy all facets. Proposition 39, or the direct generator test, then supplies invariance. \(\square\)
The switching conclusion is structural validity, not an assertion that time-varying selection represents one fixed unknown parameter. A likelihood for a fixed mechanism and a likelihood for a mechanism selected anew at each history are different statistical models.
The norm-ball inequality is sufficient and can be conservative, especially at boundary images with zero slack. Failure of this margin is not a proof that no robust certificate exists. The exact vertex test retains uncertainty correlations that a separate bound for each entry may discard.
45. Fitting inside a certified domain
Once \(V\) is fixed, certificate-constrained fitting is a convex problem. Unknowns include the entries of \(\widetilde B_s\), \(A_s\), \(\widetilde z_0\) and \(\alpha_0\). The constraints
are linear. Entrywise confidence intervals and fixed-input regression tolerances add linear inequalities. Norm balls add convex constraints. Jointly optimizing \(V\) and \(A_s\) is different: the product \(VA_s\) is then bilinear. No general convexity is claimed for learning the domain and dynamics together.
Proposition 41. Certified correction and its error budget. For fixed \(V\), every feasible solution of the displayed constraints is an internally valid model. Suppose a center estimate \(\widehat B_s\) and a feasible corrected estimate \(\widetilde B_s\) satisfy
Then \(\|\widetilde B_s-B_s^*\|\leq2\epsilon_s\). The same conclusion holds for initial states in their specified norm. If the true model lies in these confidence sets and itself admits the fixed-domain certificate, the convex constraints are feasible. Feasibility alone does not establish that the true model preserves this domain.
Proof. Validity is Proposition 39. The error bound is the triangle inequality through the center estimate. Under the final assumptions, insert the true operators, true initial state and their nonnegative witnesses into the constraints. They provide a feasible point. The converse fails because existence of one valid point does not say that every point in a confidence set, or its true member, has the same invariant domain. \(\square\)
This is why existential correction and robust certification must be reported separately. A corrected model can be a well-defined approximation even when no common true-model cone has been proved. The trajectory bounds below require both models to be valid on a common specified cone.
45.1 An explicit positivity repair in a simplicial cone
When \(V\) is square and invertible, put \(A_s=V^{-1}\widehat B_sV\). Assume mass conservation has already been imposed. Let \(N_s\) be an entrywise positive reference instrument with the same per-action mass identities. For \(0\leq\eta\leq1\), the mixture
is nonnegative provided
An empty maximum is zero. If every negative entry has magnitude at most \(d\) and every reference entry is at least \(m_0>0\), the simpler sufficient choice is \(\eta=d/(m_0+d)\). These statements follow by solving each scalar inequality. Mass conservation is unchanged by convex mixing.
This is an explicit model modification, not coefficient clipping or an exact coordinate change. Relative to a true instrument on the same cone, the one-step error is at most \((1-\eta)\epsilon+2\eta\) in the summed base-norm metric defined next, when the original error was at most \(\epsilon\). The factor two follows because two valid instruments each distribute unit total mass. The reference example retains the repair amount as part of its approximation record.
46. A norm that measures probability mass before conditioning
Assume \(C\) is a closed, pointed, generating cone and \(\lambda\) is strictly positive on its nonzero elements. Its base \(K\) is compact. Define the base norm
The infimum is finite because the cone is generating. Compactness of \(K\) makes this a genuine norm. Its unit ball is \(\operatorname{conv}(K\cup(-K))\). In particular \(\|u\|_C=\lambda u\) on \(C\), and \(|\lambda x|\leq\|x\|_C\). On the ordinary probability orthant this is the \(\ell^1\) norm. These are standard ordered-space constructions [R15], used here to avoid mixing Euclidean coefficient magnitude with probability mass.
Proposition 42. Instrument contraction before normalization. If \(\{B_{a,o}\}_o\) is cone preserving and mass conserving for a fixed action, then
For a polyhedral base \(K=\operatorname{conv}\{v_j\}\), define the discrepancy between two instruments by
The supremum equals the maximum of the same expression at the finitely many generators \(v_j\). When both instruments are valid on \(C\), \(\epsilon_a\leq2\).
Proof. For any \(x=u-v\) with \(u,v\in C\), cone preservation gives
Take the infimum. For the generator statement, the discrepancy is convex in \(z\). At a convex combination its value is at most the same combination of generator values, hence at most their maximum. The generators themselves belong to \(K\), proving equality. Finally, for \(z\in K\), the triangle inequality bounds the sum by the total masses of both instruments, which are one each. \(\square\)
For a polyhedral certificate, the norm itself has a finite linear-program formulation:
If \(V\) is square and invertible, this reduces to \(\|V^{-1}x\|_1\). The exact worked interval certificate is of this type. With dependent generators, choosing one pseudoinverse does not in general compute the base norm.
47. Whole observation laws can be stable without rare-event margins
The posterior after one rare observation can be sensitive, but that does not force the entire distribution of observation records to be equally sensitive. Work before normalization and retain the probability of each branch.
Let two models share \((C,\lambda)\) and have initial states \(z_0,\widetilde z_0\in K\). At each time use the same policy, a specified function or randomization depending only on the observed history. Thus on any common history the action probabilities are identical in the two comparisons. This is not a comparison of two policies optimized separately for the two models.
Proposition 43. Observation-tree total-variation bound. If both instruments preserve \(C\) and conserve mass, and if every action at time \(t\) has the discrepancy in Proposition 42 bounded by \(\epsilon_t\), then the distributions \(P_T,\widetilde P_T\) of length-\(T\) observation records under the common policy satisfy
The same bound holds for joint action and observation records when the common policy is randomized. No positive lower bound on individual branch probabilities is required.
Proof. Let \(u_h,\widetilde u_h\in C\) be the unnormalized states on each history at depth \(t\), including the common policy factors. Their masses sum to one in each model. Define \(D_t=\sum_h\|u_h-\widetilde u_h\|_C\). At a history with selected action \(a\), write
Sum the first term over observations and use Proposition 42. When \(\lambda\widetilde u_h>0\), normalize only to apply the definition of \(\epsilon_t\) to the second term. Its sum is bounded by \(\epsilon_t\lambda\widetilde u_h\). If the mass is zero, strict positivity implies \(\widetilde u_h=0\). Summing histories yields \(D_{t+1}\leq D_t+\epsilon_t\). Averaging over common policy probabilities gives the same inequality for randomized actions. Finally each word-probability difference is bounded by its state difference because \(|\lambda x|\leq\|x\|_C\). Total variation is half the sum of those absolute differences. \(\square\)
This bound separates total law error from conditional-state error. It can grow linearly with the horizon and eventually become uninformative. It is not a uniform-in-time approximation theorem. It also does not control the regret of a changed planning policy without additional decision-theoretic assumptions.
48. Conditioning still has an unavoidable cost
Let \(P,Q\) be probability distributions on a finite record space. A record event \(H\) might be a particular observed prefix. Conditional continuation laws are obtained by restriction to \(H\) and normalization.
Proposition 44. Sharp event-probability dependence of conditioning. If \(\|P-Q\|_{\mathrm{TV}}\leq\delta\), \(p=P(H)>0\) and \(q=Q(H)>0\), then
When \(p>\delta\), positivity of \(q\) follows automatically. The factor \(1/p\) cannot be removed in general.
Proof. Write the subprobability restrictions on \(H\) as \(u,v\), of masses \(p,q\). On the complementary event, the sum of absolute probability differences is at least \(|p-q|\). Therefore
Adding and subtracting \(v/p\) gives
Also \(|p-q|\leq\delta\), establishing the positivity condition. For sharpness, put total mass \(p\) on an event with two outcomes. Let \(P\) assign all of it to the first, and let \(Q\) move mass \(\delta\leq p\) to the second, leaving the complement unchanged. The unconditional total variation is \(\delta\) and the conditional distance is \(\delta/p\). \(\square\)
Combining this proposition with Proposition 43 gives a conditional continuation bound when the evidence event has adequate probability. It does not give a coordinate-state error bound unless test probabilities control that coordinate norm. Validity at every branch and stable inference at every branch remain different properties.
49. Searching for a domain from observable operators
A cone certificate need not be imported from a hidden-state generator. With exact rational operators it can be searched for directly, then checked independently.
Proposition 45. Rational finite-certificate search. Suppose \(\lambda\), \(B_s\) and \(z_0\) are rational and satisfy the exact mass identities. There is an enumeration procedure that returns a valid finite certificate whenever a certificate of Proposition 39 exists with a rational normalized full-row-rank generator matrix \(V\). The procedure is sound whenever it returns. No termination guarantee is asserted when no such rational certificate exists.
Proof. Enumerate all finite rational matrices \(V\) in an order that eventually includes every size and every bounded numerator and denominator range. Retain those of full row rank with \(\lambda V=\mathbf1^\top\). For each retained matrix, the requirements \(B_sV=VA_s\), \(A_s\geq0\), \(V\alpha_0=z_0\), \(\alpha_0\geq0\) and \(\mathbf1^\top\alpha_0=1\) are rational linear feasibility problems. A nonempty rational polyhedron has a rational feasible point, as follows from elimination of rational linear inequalities. These problems can therefore be decided by exact rational linear programming, or witnesses can be enumerated in a dovetailed rational search. Recheck every returned equality and sign exactly. Proposition 39 proves soundness. If a rational generator certificate exists, its finite representation is eventually considered and its feasible witnesses are found. \(\square\)
The general enumeration is a mathematical existence algorithm, not a practical complexity claim. A bounded candidate family is often the useful reference experiment. Failure within that family means only that the chosen family did not produce a certificate.
49.1 A certificate learned without hidden matrices
The observable tables of Section 38 determine the two-dimensional operators \(B_s\) and initial vector \((1,13/24)^\top\). The new reference procedure receives only these fitted operators and that vector. It enumerates intervals with rational endpoints in \([0,1]\) of reduced denominator at most eight, ordered by increasing width and then by endpoints. The coordinate restriction \([0,1]\) comes from its declared test-probability meaning.
The first accepted interval is
It uses
Every entry of every \(A_s\) is at least \(1/64\). The full per-action column sums are exactly one. These rational identities certify all word probabilities of the fitted model, not only the words used in the fit. The certificate was not supplied to the search, and no hidden generator is imported by the standalone example.
The interval is the narrowest accepted interval in this specified finite grid, not a claim of the smallest invariant interval among all real endpoints. It also need not contain every predictive state associated with an arbitrary initialization of the old hidden model. It contains the declared initial state and preserves every subsequent branch of this fitted realization, which is the requirement being certified.
The positive slack gives a robust statement: any perturbation of the lifted entries by at most \(1/64\), coupled to retain the exact per-action column sums, remains nonnegative. Interior perturbations of size \(1/128\) are used in the worked comparison. Repeated composition is valid for all horizons, while Proposition 43 bounds the changed observation law for each fixed horizon.
50. Selection and noise do not invalidate an exact witness
A data-dependent cone proposal is not itself a statistical confidence statement. Once the matrices and witness are fixed, its exact invariance identities hold or fail independently of how it was proposed. What remains statistical is whether the chosen uncertainty set covers the generating model.
For a finite-sample interpretation, first obtain a valid uncertainty set for a fixed identifiable model class. Proposition 35 supplies one route under independent reset experiments and declared rank coverage. Proposition 33 translates those probability bounds into a fixed-core operator bound. A validity witness can then be sought for a fitted member or robustly for every member. Only the latter identifies a common true-model domain on the coverage event.
The admissible family must enforce mass conservation exactly. A free box around unconstrained matrices generally contains nonnormalized families. Discarding those families is legitimate only because the generating process is already assumed to obey the exact mass law. Other cuts on a confidence set need their own justification.
51. Why finite predictive rank need not mean finitely many positive states
The following controlled rational example makes the boundary explicit. It is a worked adaptation of the irrational-rotation obstruction in [R12, Section 7]. Related distinctions among finite linear, positive and quantum realizations are discussed in [R14]. Only the elementary cone argument below is used.
Let \(z=(t,x,y)^\top\), \(\lambda=(1,0,0)\), \(z_0=(1,1,0)^\top\), and
There are two actions and two observations. A rotate action has branches
A probe action has branches
The rotate outcome zero is certain. The probe has probabilities \((t+x)/2\) and \((t-x)/2\) and resets every nonzero branch to \(z_0\) after normalization. Both actions conserve the mass \(t\).
Proposition 46. A rational rank-three model with no polyhedral probability-domain certificate. The Lorentz cone
is probability valid for these operators. It is the canonical all-test cone. Every closed convex invariant cone containing \(z_0\) and admitting the nonnegative mass interpretation equals \(C_{\mathrm L}\). In particular no finite normalized generator certificate exists for this exact model. Its reachable and observable predictive dimension is three.
Proof. Rotation preserves \(t\) and \(x^2+y^2\). On \(C_{\mathrm L}\) both probe effects are nonnegative and each probe image is a nonnegative multiple of \(z_0\). Thus the displayed cone is invariant and valid.
Write the rotation angle as \(\theta\). It is not a rational multiple of \(2\pi\). Otherwise \(\zeta=e^{i\theta}\) would be a root of unity, so \(\zeta+\zeta^{-1}=6/5\) would be a rational algebraic integer. Every rational algebraic integer is an integer, by the rational-root argument for a monic integer polynomial, contradicting \(6/5\notin\mathbb Z\). The forward powers of an irrational circle rotation are dense. Compactness supplies powers approaching the identity, so their closure contains the inverse rotation and is a closed subgroup. An infinite closed subgroup has arbitrarily small positive angles. Integer multiples of such angles approximate every point on the circle, so that subgroup is the entire circle.
A state passing every test must pass a probe after any number \(k\) of rotations. Hence
Density gives \(t\geq\sqrt{x^2+y^2}\). The reverse inclusion was already proved, so the canonical cone is exactly \(C_{\mathrm L}\).
Any invariant cone with nonnegative mass lies inside this canonical cone, by Proposition 38. If it is closed, convex and contains \(z_0\), it contains the orbit \(T^kz_0\) and the closed convex hull of that orbit. On \(t=1\), this hull is the full unit disk. Conical scaling gives all of \(C_{\mathrm L}\). The unit disk has infinitely many extreme points, so it is not the convex hull of finitely many points.
Finally the reachable columns \(z_0,Tz_0,T^2z_0\) are independent. The mass row and probe rows before zero and one rotations span the dual space because the rotation has nonzero sine. Reachability and observability establish predictive dimension three. \(\square\)
This is not a failure of finite-dimensional prediction. The exact model has only three coordinates and a simple quadratic cone. It is a failure of the narrower requirement that every valid model admit a finite polyhedral domain. A certificate language limited to finitely many positive generators would exclude a legitimate low-rank signal model.
No general impossibility of all finite certificates follows. This very example has a finite quadratic certificate. The obstruction concerns polyhedral certificates, not invariant cones described by other mathematics.
52. A finite approximation with an explicit price
A curved probability domain can sometimes be replaced by a finite one after changing the dynamics. That change must be quantified, not described as exact refinement.
Let \(P_n\) be the regular \(n\)-gon inscribed in the unit disk, with a vertex at \((1,0)\) and \(n\geq3\). Replace the rotate branch by
while retaining the same probe branches. The spatial coordinates are attenuated, but total probability mass remains unchanged.
Proposition 47. Polyhedral approximation of the rotation model. If
then the cone over \(P_n\) is invariant under all modified operators and contains \(z_0\). Consequently the modified model admits a finite positive lift. The original and modified models are both valid on the larger cone \(C_{\mathrm L}\), and for every common policy their length-\(T\) laws satisfy
The angle-independent sufficient condition can be improved for a fixed polygon and rotation. For the rational diamond \(|x|+|y|\leq t\), the present rotation is invariant after attenuation whenever \(\alpha\leq5/7\). Thus \(\alpha=2/3\) gives an entirely rational four-generator certificate.
Proof. The regular polygon contains the disk of radius \(\cos(\pi/n)\) and lies in the unit disk. Every \(\alpha Rv\) with \(v\in P_n\) has Euclidean norm at most \(\alpha\), hence lies in \(P_n\) under the stated condition. The probe effects remain nonnegative because \(P_n\) lies in the disk, and their normalized images are the retained vertex \((1,0)\).
The base norm of \(C_{\mathrm L}\) is
Indeed the convex hull of the unit disks in the planes \(t=1\) and \(t=-1\) is the cylinder with this unit ball. The rotate-branch discrepancy on a normalized state is at most \(1-\alpha\) in this norm. Probe discrepancies are zero. Proposition 43 gives the law bound. For the diamond, the induced \(\ell^1\) norm of \(R\) is \(7/5\), so \(\alpha\|R\|_1\leq1\) proves invariance. \(\square\)
For a target attenuation loss \(0<\eta\leq1\), choosing \(\alpha=1-\eta\) and sufficiently many polygon vertices satisfies the certificate. The inequality \(1-\cos u\leq u^2/2\) shows that \(n\geq\pi/\sqrt{2\eta}\) is sufficient, with an integer ceiling and \(n\geq3\). This is a constructive upper bound on certificate size, not a lower bound or a proof of optimality.
The distinction is exact. At \(\alpha=1\) there is no finite polyhedral certificate for this model. At \(\alpha<1\) a sufficiently fine polygon works, at the cost of a declared change in the observation law. The reference figures show that tradeoff rather than labeling the approximation an identical model.
53. Repeated proposal and verification with a confidence budget
Model development may choose another core, propose another domain or collect another independent batch after an earlier certificate is withheld. Statistical coverage must survive that selection process.
Proposition 48. Sequential confidence accounting for certified proposals. At stage \(n\), let the proposal be measurable from previous data, and let a verification procedure return a set \(\mathcal U_n\) for the true model in the declared identifiable coordinate system. Assume conditional coverage
Then with probability at least \(1-\delta\), every returned uncertainty set covers its true model. On this event, any data-selected domain with an exactly verified robust certificate for all members of its corresponding \(\mathcal U_n\) is also valid for the true model, provided the true initial state is covered by the certified base. The conclusion remains valid at any stopping time chosen from the recorded results.
Proof. Take expectations of the conditional failure bounds, then apply the countable union bound. Off the union of failure events, all coverages hold simultaneously. Robust invariance over an entire set then applies to its true member by deterministic implication, independently of how its domain was chosen. A stopped stage is one of the already covered stages. \(\square\)
For example, \(\delta_n=\delta 2^{-n}\) spends a finite total confidence budget. Fresh reset trials for a proposal fixed by past data can use Proposition 35 conditionally, with this stage budget. An exact witness needs no statistical correction for how many candidate cones were tried. The coverage statement does need correction for repeatedly constructing confidence sets or inspecting post-selected data.
The proposition does not infer a global rank bound, validate an unknown reset operation, or transform dependent observations into independent trials. Those are assumptions of the coverage procedure. A known coordinate transformation can transport a confidence set exactly. An estimated transformation requires the subspace and conditioning controls already separated in Part III.
54. Integrated certificate-carrying realization
Proposition 49. Conditional learning with an all-horizon probability certificate. Suppose the observable identification assumptions of Proposition 37 hold for a declared finite-rank controlled process. Suppose a justified confidence procedure produces an operator and initial-state set containing the true realization on an event of probability at least \(1-\delta\). Assume the per-action mass constraints are exact.
If a fitted realization admits a checked finite witness of Proposition 39, it is an internally valid stochastic model at every horizon. If the witness is robust over the whole confidence set and covers its initial states, then on the coverage event the same probability domain is valid for both the fitted and true models. When their summed base-norm instrument discrepancies are at most \(\epsilon_t\) and their initial discrepancy is at most \(e_0\), every common policy has the bound
Conditioning on a positive-probability event then obeys Proposition 44. Exact coordinate transport preserves the certificate. A signed packet decomposition preserves it after complete synthesis but does not turn individual packets into probability states without Proposition 30's extra positivity requirement.
Proof. Observable identification supplies the specified realization and the confidence procedure supplies the stated event, not an inferred event of its own. Proposition 39 supplies internal validity for the fitted model. Proposition 40 extends the same domain to every member of the uncertainty set, hence to its true member on that event. Propositions 42 and 43 give the total-variation bound, and Proposition 44 gives the conditional bound. For a coordinate change \(T\), use \(V'=TV\), \(B_s'=TB_sT^{-1}\) and \(\lambda'=\lambda T^{-1}\). Then \(B_s'V'=V'A_s\) and \(\lambda'V'=\mathbf1^\top\) with the unchanged nonnegative witnesses. Full packet synthesis is the identity. \(\square\)
This closes a concrete foundational gap. A learner can now return not just a matrix family but a verifiable domain in which that family has probability meaning. It can attach a robust uncertainty statement when one is justified, or return an internally valid approximation with its correction budget explicitly recorded.
It does not turn a successful fit into unrestricted knowledge of the world. The model class, observation protocol, rank coverage, exact mass law, domain witness and approximation budget remain visible parts of the result. Failure of a finite search is a withheld certificate, not a proof of invalidity. The nonpolyhedral example explains why this distinction is necessary even for very small exact models.
55. From polyhedral to curved certificates
A finite positive lift is one sufficient certificate language. Proposition 46 shows that a low-dimensional model can require a nonpolyhedral domain. Part V represents that same model exactly through a curved matrix domain. A failed finite-generator search is therefore not a general invalidity theorem.
\newpage
Part V. Curved domains and continuous observations
Part IV separated a valid predictive model from a successful fit. Its finite positive lift gives a useful certificate, but its rank-three rotation example has no finite polyhedral probability domain. The same block left continuous observations outside the finite-alphabet construction. This part addresses both boundaries in a declared finite sector.
A curved domain can have a short matrix certificate even when it has infinitely many extreme states. A continuous observation can have a finite operator description even when its possible values are uncountable. Neither statement follows from finite predictive dimension alone. The constructions below require an explicit domain, a reference measure and a restricted observation-kernel family.
The block connects matrix positivity, integral observation semantics, population-moment identification, exact coarsening and approximation. It also identifies a defect that is easy to overlook: averaging a matrix transformation before applying its quadratic action does not equal averaging that action. The difference is a positive covariance contribution. Omitting it can destroy normalization even in a very small example.
There are established antecedents. Congruence representations and their positive semidefinite Gram matrices belong to the theory of completely positive maps [R16]. Polynomial Gram certificates are standard sum-of-squares machinery [R17]. Continuous-valued observable operator models already occur in Jaeger's work [R18]. The proofs here develop their specific connection to the realization and defect contracts of the large signal model. They do not claim to originate those ingredients.
This is a discrete-time theory with continuous-valued observations. It is not a continuous-time stochastic dynamics theorem. Production PR-LASM 2.3.0 is unchanged.
56. A matrix certificate for a curved domain
Let \(E=\mathbb R^r\) and let \(\mathbb S^d\) be the real symmetric \(d\times d\) matrices. Fix an injective linear map
Suppose some \(e\in E\) satisfies \(\mathcal J(e)\succ0\). Define
The map \(\mathcal J\) is an embedding, not a many-to-one decoder of a larger positive coordinate space. The cone is an inverse image of the positive semidefinite cone. This is different from the image cone \(V\mathbb R_+^m\) used in Proposition 39.
Proposition 50. Matrix-embedding probability certificate. The cone \(C_{\mathcal J}\) is closed, pointed and generating, and \(\lambda\) is strictly positive on its nonzero elements. Its normalized base is compact. Suppose a finite-alphabet predictive model has \(z_0\in K_{\mathcal J}\), exact mass conservation and identities
Then every \(B_s\) preserves \(C_{\mathcal J}\), and every finite observation word has a nonnegative, normalized probability. Zero maps may use an empty sum. An invertible change of predictive coordinates transports the certificate by
Proof. Closedness follows from continuity. If both \(z\) and \(-z\) are in the cone, \(\mathcal J(z)\) is both positive and negative semidefinite, hence zero. Injectivity gives \(z=0\). A sufficiently small neighborhood of \(e\) maps into positive definite matrices, so the cone has an interior point and is generating. A nonzero positive semidefinite matrix has positive trace. The trace-one positive semidefinite matrices form a compact set. The image of \(\mathcal J\) is a closed finite-dimensional subspace and its inverse on that image is continuous, so \(K_{\mathcal J}\) is compact.
Every term \(wKXK^{\mathsf T}\) is positive semidefinite when \(X\) is. The intertwining identity therefore proves invariance. Iteration gives invariance of every branch state, and the mass identity normalizes every finite observation tree. Substituting the transported quantities proves the final statement. \(\square\)
This is a sufficient certificate language. It is not asserted that every invariant cone has this form, or that every positive map on a matrix subspace extends through the displayed congruence representation. The extension requirement is stronger than positivity on the represented states alone.
The base norm of Section 46 is also not automatically the ambient trace norm. One always has
because every positive decomposition in \(E\) gives one in the matrix space. Equality holds when \(\mathcal J\) maps onto all of \(\mathbb S^d\), by the positive and negative spectral parts. For a proper matrix subspace, those spectral parts may leave the subspace. Norm transport therefore needs its own stated hypothesis.
57. The earlier disk obstruction has an exact finite matrix witness
For the rank-three model of Section 51, put
Then \(\operatorname{tr}\mathcal J=m\) and
Its eigenvalues are \((m\pm\sqrt{x^2+y^2})/2\), so positivity is exactly the Lorentz-cone condition. The normalized base is the unit disk.
Proposition 51. Exact rational matrix witness for the rotation model. Let \(T=\operatorname{diag}(1,R)\) with the unchanged rotation
The rotate branch and the two probe branches from Proposition 46 have the identities
where \(E_{ij}\) is the corresponding \(2\times2\) matrix unit. The other rotate branch is zero. Thus this exact model has a finite rational matrix certificate without attenuation or a polygonal replacement.
Proof. Direct multiplication gives the rotate identity. The two matrix-unit congruences retain, respectively, the upper and lower diagonal entries of \(\mathcal J(z)\) and place them in its upper diagonal entry. These are \((m+x)/2\) and \((m-x)/2\), exactly the two probe masses resetting to \((1,1,0)\). Also \(W^{\mathsf T}W=5I\) and
The actions conserve trace. The initial state embeds as \(\operatorname{diag}(1,0)\), and \((1,0,0)\) is a strictly positive interior witness. Proposition 50 applies. \(\square\)
Nothing has been changed in the observation law. Proposition 46 still rules out a finite positive-generator certificate for this model. What changes is the certificate language. A curved matrix domain and a finite set of positive coordinate states are not the same representation class.
This use of symmetric positive semidefinite matrices is mathematical. It neither requires quantum hardware nor claims a quantum description of the measured system. It borrows a well-defined operator representation rather than a physical interpretation.
58. An observation density is not a branch probability
Let \((Y,\mathcal Y)\) be a standard Borel observation space and let \(\mu\) be a fixed probability reference measure. Keep a finite action set. For each action \(a\), let \(B_a(y):E\to E\) be measurable and integrable in operator norm. Let \((C,\lambda)\) be a proper probability cone and mass functional as in Section 46. Assume that, outside one \(\mu\)-null set for each action,
The second identity is an equality of linear functionals. The first is a statement on the whole cone, not only on a finite selection of states.
Proposition 52. Continuous-observation instrument semantics. For \(z\in K=\{z\in C:\lambda z=1\}\),
define a probability kernel and, for \(\mu\)-almost every observation with positive denominator, a conditional state. They generate consistent finite record laws under any common measurable history-dependent policy. For prescribed actions, the joint observation density is
with respect to \(\mu^{\otimes T}\). Replacing the reference measure by \(\mu'=r\mu\), with \(r>0\) almost everywhere and \(\int r\,d\mu=1\), requires \(B'_a(y)=B_a(y)/r(y)\). Event probabilities and normalized states then remain unchanged.
Proof. Cone positivity makes the density nonnegative and the integrated mass identity makes its integral one. A positive-density update lies in \(K\). Assign any fixed state of \(K\) on the common exceptional null set and at zero-density points to obtain a measurable version of the update. This arbitrary choice has no effect on the probability kernel. Repeated integration, or conditioning one step at a time, proves normalization and consistency of the finite laws. Measurable policy factors are nonnegative and sum to one at each history. The reference-measure transformation cancels \(r\) in each integral and cancels \(1/r\) in the normalized update. \(\square\)
When \(\mu\) has no atoms, an individual point has probability zero even when its density is positive. A density may exceed one. Neither fact is a numerical pathology. Flooring a vanishing denominator changes the model and is not part of this theorem.
A continuous density at an exact point is a version-dependent object on null sets. The update formula specifies a convenient version wherever its denominator is positive. Event probabilities are the invariant statistical objects. Comparing model densities requires a common measure or a declared change-of-measure calculation.
A representation that lacks an operator-norm integrable density needs a more general operator-valued measure. Such representations are not silently covered by the density hypothesis above.
59. A finite certificate over an entire observation interval
Specialize to \(Y=[-1,1]\) and \(\mu(du)=du/2\). Suppose \(B_a(u)\) is a matrix polynomial. A finite congruence certificate can prove positivity for every \(u\), not merely for values on a grid.
For a real \(d^2\times d^2\) symmetric matrix \(F\), define
Vectorization orders the pairs with output index first inside each input column. If
then \(\Phi_F(X)=\sum_\nu w_\nu K_\nu X K_\nu^{\mathsf T}\). Every positive semidefinite \(F\) has a real factorization of this kind. The matrix is the usual Choi/Gram representation, with an explicit index convention [R16].
Proposition 53. Polynomial Gram certificate for interval observations. Fix the injective embedding \(\mathcal J\) and its strict interior witness from Proposition 50. Set
Suppose each action has constant positive semidefinite matrices \(G_{a,0},G_{a,1}\), with the second term omitted when \(q=0\), such that
Assume also \(z_0\in K_{\mathcal J}\) and the exact integrated mass identity. Then the model is internally probability-valid for every observation value in \([-1,1]\) and every finite horizon. The intertwining condition is finitely many coefficient equalities. The normalization condition uses only
Proof. Both Gram expressions are positive semidefinite on the interval because \(1-u^2\geq0\) there. Factoring \(F_a(u)\) at a fixed \(u\) gives a positive congruence map, whether or not its individual factors depend polynomially on \(u\). The intertwining identity therefore proves \(B_a(u)C_{\mathcal J}\subseteq C_{\mathcal J}\). Polynomials on a compact interval are integrable. The mass identity and Proposition 52 prove the assertion. Equality of two finite matrix polynomials is equivalent to equality of every coefficient. The moment values follow by elementary integration. \(\square\)
For fixed \(\mathcal J\) and fixed degree bounds, unknown entries of \(B_a\) and the Gram matrices obey linear equalities and positive semidefinite constraints. A convex data-fit objective therefore gives a semidefinite feasibility or optimization problem in this fixed representation. Jointly choosing \(\mathcal J\) changes the problem: the intertwining equation then contains products of unknown coefficients.
This is a sufficient certificate, not a completeness theorem for all positive kernels or all curved domains. A negative answer at one degree is not a proof of model invalidity. Approximate Gram feasibility is also not an exact witness. Rational factorizations with nonnegative rational weights, or rigorous matrix-positivity bounds, must survive a separate replay.
The reference checker accepts a factored version:
It checks the polynomial identities, weights, support, embedding, initial state and normalization exactly. It does not run a general semidefinite optimizer. The operator identity holds on the image of \(\mathcal J\); trace preservation of an arbitrary extension outside that image is not required.
60. Learning from integral observations rather than point probabilities
A finite polynomial family can be identified through moment operators. This gives a continuous-observation analogue of the core-table construction without pretending that a point has positive probability.
Assume the true operator density belongs to the declared class
Define the moment Gram matrix \(W_D\) by
A history is now a prescribed action sequence with observations constrained to measurable events of positive or zero probability. It produces an unnormalized predictive vector \(h_j\). A future event experiment produces a row functional \(t_i\). Choose \(r\) of each and form
The index \(j\) in \(G_{a,j}\) is a moment order. Odd-order entries can be signed expectations. They are not branch probabilities and are never normalized as if they were.
Proposition 54. Conditional finite-moment identification. Suppose the true predictive dimension is \(r\), the degree bound \(D\) is correct, and the selected \(r\times r\) core \(H\) is invertible. Exact population values of \(H\), every \(G_{a,j}\) and the initial test vector determine the complete polynomial operator family in test coordinates:
If every moment-operator error is at most \(\eta\) in one fixed induced norm, then coefficient reconstruction alone obeys
Proof. A nonzero polynomial has strictly positive squared integral, so its coefficient vector \(c\) satisfies \(c^{\mathsf T}W_Dc>0\). Thus \(W_D\) is invertible. Let \(T\) have rows \(t_i\) and \(S\) have columns \(h_j\). Since \(H=TS\) is an invertible square core in an \(r\)-dimensional space, both \(T\) and \(S\) are invertible. Therefore \(G_{a,j}H^{-1}=TM_{a,j}T^{-1}\). The relation \(M_{a,j}=\sum_k(W_D)_{jk}C_{a,k}\) is exact integration of the polynomial expansion. Invert it and conjugate by \(T\) to obtain the displayed reconstruction. The error bound is the triangle inequality for that same fixed linear inversion. \(\square\)
The rank and degree assumptions are not learned from the determinant of \(H\). The determinant establishes coverage only after the ambient predictive dimension has been justified. The moment Gram can become poorly conditioned at higher degree. An orthogonal polynomial basis can improve this representation without removing an incorrect degree assumption.
The entries have an operational population meaning. A reset experiment for an event word, with one observation weighted by \(u^j\), returns a bounded random variable in \([-1,1]\). Its expectation is the corresponding entry of \(G_{a,j}\). With \(N\) independent trials for each of \(M\) predeclared entries, the same exponential-moment argument as Proposition 35 gives the simultaneous bound
with failure probability at most \(\delta\). The factor differs from the Bernoulli bound because the value range is two. Dependent trajectories, unknown resets and post-selection are not covered by that observation protocol. This edition uses exact constructed population tables, not a new sampling experiment.
61. Finite moments cannot validate an unrestricted density class
A finite observation basis is a modeling assumption. Positivity certification does not turn that assumption into a consequence of data.
Proposition 55. A valid continuous counterexample to finite-moment identification. Let
Relative to \(du/2\), these are distinct strictly positive normalized densities with the same moments of orders zero, one and two. Their total-variation distance is
Taking independent observations gives two rank-one processes whose products of these low-order moments agree at every sequence length.
Proof. The endpoints and stationary points of \(L_3\) show \(|L_3(u)|\leq1\) on the interval. Direct integration gives \(\int u^jL_3(u)\,du/2=0\) for \(j=0,1,2\). This proves normalization and moment agreement, while strict positivity follows from \(\epsilon<1\). Its nonnegative root other than zero is \(\sqrt{3/5}\). Integrating \(|L_3|\) on the two subintervals of \([0,1]\) and using symmetry gives \(\int|L_3|\,du/2=13/40\). The density difference is \(2\epsilon L_3\), so its half-integral is the stated distance. Independence makes every product moment factor into the identical one-step moments. \(\square\)
This is not a defect in Proposition 54. These densities lie outside its quadratic class. It is a concrete reason to retain an explicit approximation or misspecification term when increasing the observation basis. Internal probability validity and correct density identification remain different claims even in one predictive dimension.
62. Exact observation bins preserve their own event law
Let \(A_1,\ldots,A_n\) be a measurable partition of \(Y\). Define the cell operators
Their role is to represent observed events. They are not point samples of the original density.
Proposition 56. Exact bin coarsening and its posterior. Every \(B_{a,j}\) preserves the probability cone, and \(\lambda\sum_jB_{a,j}=\lambda\). For a policy measurable with respect to bin histories, the finite model \(\{B_{a,j}\}\) reproduces exactly the binned observation law of the continuous model at every finite horizon. On an observed cell with \(p_j=\lambda B_{a,j}z>0\),
Proof. Integrals of vectors in a closed convex cone remain in the cone. Summing over a partition gives the original mass identity. For a fixed bin history, its policy factors are constant across the corresponding cells. Successive integration of the operator product therefore gives the product of the cell operators, in the original temporal order. This establishes every finite binned record probability. In the posterior identity, substitute \(z^+(y)\lambda B_a(y)z=B_a(y)z\) on the positive-density set. The remaining outputs are zero because they lie in the cone and have zero mass. \(\square\)
The coarsened posterior is a conditional average, not the posterior at a representative observation value. It preserves exactly the information retained by the event history. An arbitrary policy that reacts differently to two observations inside the same cell is not bin-measurable, so the preceding exact-policy conclusion does not apply to it.
Cell refinement is compatible with this semantics: a parent operator is the sum of its children. In this setting those cell operators are positive probability branches. This positivity is an extra property of event integration, not a property of every signed spectral packet in the signal archive.
63. Averaging transformations has a covariance defect
A common discretization replaces a family of matrices \(K(y)\) by its mean and then applies \(X\mapsto KXK^{\mathsf T}\). This is generally the wrong mean operation.
Proposition 57. Exact covariance correction for an averaged congruence. Let \(\nu\) be a probability measure and suppose \(K(y)\) has square-integrable entries. Put \(\bar K=\int K(y)\,\nu(dy)\). For every symmetric \(X\),
The correction is positive semidefinite when \(X\succeq0\). For an affine matrix function \(K(u)=K_0+uK_1\) with scalar \(u\) of finite second moment, it equals \(\operatorname{Var}_\nu(u)K_1XK_1^{\mathsf T}\). On a uniform interval of width \(h\), the variance is \(h^2/12\).
Proof. Substitute \(K=\bar K+(K-\bar K)\) and expand. The two cross terms integrate to zero. The remaining integrand is a positive congruence when \(X\) is positive. For an affine matrix function the centered term is \((u-\mathbb E_\nu u)K_1\), giving the variance formula. Integrating the square of a centered uniform variable gives \(h^2/12\). \(\square\)
This is a defect between two operations: averaging a transformation parameter and applying its quadratic action. It is not the same as the residual from dropping a signal packet, although both fit the framework's requirement to record failures of commutation.
For the example in Section 66, equal bins with midpoint evaluation lose total mass
when there are \(n\) bins. With one bin, midpoint evaluation retains only one quarter of the required mass. Exact integration or the covariance correction restores the missing contribution. Dividing an inaccurate approximation by a guessed normalizer would define another model rather than establish this identity.
64. A valid density approximation with a resolution bound
Exact coarsening gives a finite event model. To compare complete continuous observation laws, lift that model back to a density on the same observation space. Assume all cells have positive reference mass and put
This uses the normalized reference measure inside each cell. It is an explicit density approximation, not an assertion that the original conditional density was constant there.
Proposition 58. Positive cell-average approximation. The family \(B_a^\Pi\) is cone preserving and mass conserving. It has exactly the same cell operators as \(B_a\). Let \(Y\) have a metric and suppose
in the induced base norm, with every cell of diameter at most \(h\). Then
For equal intervals of width \(h\) under the uniform measure on \([-1,1]\), the sharper bound is \(\epsilon_\Pi\leq Lh/3\).
Proof. Each cell average is a positive scalar multiple of a positive cell operator. Its integral over that cell is \(B_{a,j}\), proving both normalization and exact cell agreement. For \(y\in A_j\), subtract the cell average and use
because \(\|z\|_C=1\). Integrating over \(y\) proves the diameter bound. Two independent uniform points in an interval of width \(h\) have mean absolute separation \(h/3\), obtained by integrating over its two triangles. Applying this inside every equal interval gives the sharper bound. \(\square\)
The Lipschitz constant is a declared operator bound. It is not a fit residual inferred from adjacent samples. The comparison needs a common domain and a common reference measure. If two kernels are only close in Euclidean coefficient norm, the relevant norm-conversion constant must be included.
In the full \(2\times2\) symmetric-matrix example, the pullback base norm is
There \(\|K(u)\|_2\leq2\) and \(\|K(u)-K(v)\|_2\leq(3/2)|u-v|\). The identity
and \(\|AXB\|_1\leq\|A\|_2\|X\|_1\|B\|_2\) give \(L=6\). The rotated action uses an orthogonal matrix before this map and has the same bound. Equal bins therefore give \(\epsilon_\Pi\leq4/n\), capped by two for two valid instruments. The inequality is sufficient, not asserted sharp.
65. Full continuous record laws inherit the stability argument
The base-norm proof from Part IV has an integral version. Its hypotheses must be supplied rather than inferred by replacing a finite sum with an integral in a program.
Proposition 59. Continuous-record total-variation bound. Let \(B_a(y)\) and \(\widetilde B_a(y)\) satisfy Proposition 52 on a common \((C,\lambda,\mu)\), and let the two models use the same measurable policy on every common observation history. Suppose their initial states satisfy \(\|z_0-\widetilde z_0\|_C\leq e_0\) and their time-\(t\) instruments satisfy
The length-\(T\) record laws, including actions when the policy is randomized, obey
No lower bound on each observation density is required. For identical initial states and the equal-bin approximation in Section 64, the example has the sufficient bound \(\min(1,2T/n)\).
Proof. For any positive decomposition \(x=u-v\), cone preservation and mass conservation give
Taking the infimum proves integrated instrument contraction. At each continuous history let \(v_h,\widetilde v_h\) be the unnormalized state densities, including the policy factors. Split the next difference as
Integrate the first term using contraction. The second is bounded by \(\epsilon_t\lambda\widetilde v_h\) after normalizing only its nonzero positive mass. Its integral over histories is at most \(\epsilon_t\). Tonelli's theorem applies to these nonnegative norm bounds, and operator integrability supplies well-defined finite-dimensional integrals. The integrated state discrepancy therefore satisfies \(D_{t+1}\leq D_t+\epsilon_t\). Finally \(|\lambda x|\leq\|x\|_C\) bounds the record-density difference, and total variation is half its absolute integral. \(\square\)
This result compares one common policy, not two policies optimized separately. It concerns the whole law and can become uninformative at long horizons. It is not a uniform posterior contraction theorem.
For a measurable event \(H\) of probability \(p=P_T(H)>0\), the restriction-and-normalization argument of Proposition 44 works for measures as well as finite vectors. When both conditional laws exist, their total-variation distance is at most \(\min(1,\delta/p)\) if the full-law distance is at most \(\delta\). The event probability is invariant under a reference-measure change. A density value at one point is not a substitute for that event probability.
The full-law comparison between \(B\) and \(B^\Pi\) permits any common measurable policy because both are continuous-output models on the same space. The exact equivalence to the finite bin model in Proposition 56 instead requires a bin-measurable policy. These are different conclusions with different observation information.
66. A complete continuous-observation example
Use the same three-coordinate matrix embedding as Section 57. Set
The direct action uses \(X\mapsto K(u)XK(u)^{\mathsf T}\). The rotated action first applies the unchanged rotation congruence from Section 57, giving
Both are exactly normalized because \(\int K(u)^{\mathsf T}K(u)\,du/2=I\) and \(W^{\mathsf T}W=5I\). In predictive coordinates, write
Although \(c(u)\) is negative for part of the interval, both actions preserve the disk cone. Entrywise sign is not the positivity criterion in these coordinates. The initial state is \(z_0=(1,0,0)^{\mathsf T}\), or the matrix \(I/2\).
66.1 Event tables and moment reconstruction
Let the three future tests be the empty test, a nonnegative observation after action zero, and a nonnegative observation after action one. Their rows form
Use histories consisting of the empty history, a nonnegative observation under action zero, and that history followed by a nonnegative observation under action one. Their unnormalized columns are
The learner is supplied only the resulting integral tables and the declared class. Its core is
The core is not required to be symmetric. Its rows are tests and its columns are histories, which are different experiments. The three moment orders use
Proposition 54 reconstructs \(QB_a(u)Q^{-1}\) from \(H\) and the six inserted-moment tables. The code that performs this reconstruction imports no generating matrices. The separate generator produces the tables and is used only in the verification suite. These are population values from an exact example, not learned sample estimates.
66.2 Two densities and a genuine event probability
For a direct observation value \(u=1/2\), the density and updated state are
If the rotated action next gives \(u=-1/2\), they are
Their product \(413/1280\) is a joint density with respect to \((du/2)^{\otimes2}\), not the probability of those two exact values. The two point values each have zero probability in this nonatomic model.
By contrast, the event \(u_1\geq0\) under action zero followed by \(u_2\geq0\) under action one has actual probability
and its conditional state is \((1,96/107,-24/107)^{\mathsf T}\). Exact bin operators produce this value without quadrature.
The final point-conditioned state remains strictly inside the disk:
A pure boundary state can instead have zero density at isolated observations. The checker rejects conditioning on such a zero-density branch. It does not manufacture a posterior by flooring its denominator.
66.3 The missing covariance is exactly calculable
On an interval with midpoint \(v\) and width \(h\), the average matrix instrument is
The correction is itself a positive congruence. The cell operator multiplies this mean by the cell mass \(h/2\). Equal bins have width \(2/n\), so omitting the covariance loses trace \(3/(4n^2)\) over the full partition. The reference calculation checks both the exact matrix identity and this mass defect for several rational partitions.
67. A declared template can be recognized without being rediscovered
General curved-domain learning remains difficult. There is a useful finite case between supplying a finished coordinate certificate and claiming unrestricted domain discovery.
Suppose a structural library supplies a certified template with polynomial coefficients \(C_{a,k}\), mass row \(\lambda\) and embedding \(\mathcal J\). Suppose observed population tables produce coefficients \(\widehat C_{a,k}\) and mass row \(\widehat\lambda\) of the same dimension.
Proposition 60. Exact recognition of a fixed certified template. The equations
are linear in \(T\). Every invertible solution transports the template certificate to the fitted coordinates by \(\widehat{\mathcal J}=\mathcal JT^{-1}\). If the fitted coefficient family is observable, meaning the rows \(\widehat\lambda\widehat C_w\) span the dual space, the solution is unique. The fitted initial state must still be tested in the transported normalized cone.
Proof. All coefficient matrices are fixed, so the displayed constraints are linear. An invertible solution yields \(\widehat B_a(u)=TB_a(u)T^{-1}\) coefficient by coefficient. Substitution in the template intertwining identity proves the transported certificate. To prove uniqueness, let \(D\) be the difference of two solutions. It satisfies \(\widehat C_{a,k}D=DC_{a,k}\) and \(\widehat\lambda D=0\). Therefore \(\widehat\lambda\widehat C_wD=0\) for every coefficient word. Observability makes every column of \(D\) zero. An initial state outside the transported cone is not repaired by this operator identity, which explains the separate check. \(\square\)
For the worked example, the linear system has the unique solution \(T=Q\) from Section 66. The learner receives the population tables and a fixed structural template. It reconstructs the operators, solves the linear match, transports the embedding and exports a certificate. Independent replay then needs only that certificate. It does not refit, search, or import the template or generator.
This is template recognition, not unrestricted geometry discovery. The supplied template includes the specific operator family. An incompatible model is rejected rather than changed to fit the template. The reference routine also withholds a result for an underdetermined or singular match. Noise would require a declared approximate match and a separate validity or correction argument. This edition's recognition example is exact.
For fixed embeddings, polynomial Gram optimization offers a broader construction route than fixed-template matching. That convex proposal problem is stated in Section 59 but is not implemented as a general solver here. The independent checker remains useful regardless of whether a proposal came from a template, symbolic derivation or an optimization routine.
68. The combined continuous polynomial sector
Proposition 61. Conditional continuous realization with a finite domain certificate. Assume a controlled continuous-observation process has the declared finite predictive dimension and polynomial operator-density degree, a fixed uniform reference measure on \([-1,1]\), and an invertible complete event-test core. Exact moment tables identify its operator density in test coordinates by Proposition 54.
If the fitted family and initial state admit a checked matrix/Gram certificate of Proposition 53, the supplied realization is internally valid at every observation value and finite horizon. A successful invertible fixed-template recognition is one sufficient way to construct such a certificate. Exact coordinate transport and complete packet synthesis preserve the represented event probabilities. Individual signed state packets still need the separate positivity condition before probabilistic interpretation.
Exact observation-cell integration gives the precise finite event model for bin-measurable policies. Positive cell-average density approximations have the resolution bound of Proposition 58. Two certified models on a common cone and common reference measure satisfy Proposition 59's full-law bound whenever the declared instrument discrepancy is controlled.
Proof. Moment identification uses the true dimension and degree hypotheses to apply Proposition 54. The independent Gram identity establishes positivity on the whole declared interval, while the exact integral mass identity supplies normalization through Proposition 52. Proposition 60 supplies a transported witness in the template case. Coordinate cancellation and full packet synthesis are exact identities. Proposition 56 proves the cell-model statement. Propositions 58 and 59 provide the two approximation conclusions with their common-domain and common-policy assumptions unchanged. \(\square\)
For a noisy fitted model, the certificate alone proves only internal validity. A guarantee about the generating process requires a valid uncertainty set covering it. A common-domain certificate over a convex uncertainty family can be checked at complete-family vertices with the same embedding, because the operator identities, Gram positivity and mass constraints survive convex combination. This does not justify an uncertainty set that excludes the true degree, rank or dynamics.
The construction leaves visible what is often collapsed into the phrase signal model: a represented state, an operator-valued observation density, a reference measure, a valid probability domain, an observation protocol and an approximation budget. Each is needed for a different reason.
69. From continuous semantics to discovery
Part V gives event-based probability semantics, exact binning and conditional finite-moment reconstruction. Fixed-template recognition in Proposition 60 is not unrestricted geometry discovery. Parts VI–VII add restricted discovery and robust verification with their own assumptions. The continuous reference measure and the observation basis remain explicit inputs.
\newpage
Part VI. Certificate discovery, approximation and sequential evidence
A signal model must distinguish a representation that preserves probability from a representation that predicts the observed process accurately. Part V gave an exact curved certificate and a way to recognize a supplied operator template. It also made the polynomial-degree assumption explicit. This block advances both boundaries without identifying them with each other.
The first development is a restricted discovery procedure. It recovers an ellipsoidal probability domain from fitted polynomial operators, without a supplied template, matrix embedding or hidden generator. The second is a positive approximation theorem for operator densities outside a finite polynomial class. The third is sequential evidence under dependent observations, including the cost of action coverage. Together they separate structural validity, approximation error and observable disagreement.
The invariant-form and quadratic-certificate ingredients belong to classical systems theory [R8, R13]. Positive polynomial averaging uses the classical Bernstein--Durrmeyer kernel [R19]. The sequential argument is an elementary martingale confidence sequence, not a new concentration principle [R20]. The contribution of this block is their connection to the large signal model's explicit representation contracts and to independently replayable examples. Propositions 62--70 have proofs below. No proof kernel or production training system is invoked.
70. An ellipsoidal certificate with an orientation test
Work in mass-first coordinates \(z=(m,x)\in\mathbb R\times\mathbb R^d\), with \(\lambda=(1,0,\ldots,0)\). For a center \(c\in\mathbb R^d\) and a symmetric positive definite matrix \(H\), define
Its normalized base is an ellipsoid. Set
Then \(C_{c,H}=\{z:m\geq0,\ z^{\mathsf T}Qz\geq0\}\). The cone is proper and generating, and \(m\) is strictly positive on every nonzero cone element. These facts follow by an invertible linear transformation to the Lorentz cone. In particular, an output with zero mass in the cone is the zero vector.
Let \(B_a(u)\) be polynomial operator densities on \([-1,1]\), relative to \(du/2\). Write
Proposition 62. Oriented quadratic probability certificate. Suppose that for each action and every \(u\in[-1,1]\) there are \(\tau_a(u)\geq0\) and \(R_a(u)\succeq0\) satisfying
and that
If \(z_0\in C_{c,H}\) has mass one and \(\lambda\int B_a(u)\,du/2=\lambda\), then every operator preserves \(C_{c,H}\). The supplied model has valid probabilities at every finite horizon. Polynomial sum-of-squares witnesses for the scalar inequalities and for \(R_a\) make this a finite certificate.
Proof. For \(z\in C_{c,H}\), the quadratic identity gives
For a nonzero \(z\), write \(x=m(c+w)\) with \(w^{\mathsf T}Hw\leq1\). Its output mass is \(m(\beta_a+b_a^{\mathsf T}w)\). Cauchy--Schwarz in the \(H\) inner product bounds this below by
Thus both inequalities defining the desired sheet of the quadratic cone hold. The zero input causes no difficulty. Integrated mass conservation gives normalization, and induction gives the finite-horizon conclusion. For polynomial data, identities between coefficients and nonnegative sums of polynomial squares establish the displayed hypotheses throughout the interval. Terms multiplied by \(1-u^2\) are also nonnegative on that interval. \(\square\)
The orientation test is essential. The matrix \(-I\) preserves every quadratic form, but maps positive mass to negative mass. A quadratic inequality without its mass orientation can certify the wrong sheet. Likewise, testing finitely many values of \(u\) is not a substitute for a polynomial identity and its positivity witness.
This is a sufficient language, not a completeness theorem for all cones. The reference checker accepts rational scalar and matrix sum-of-squares witnesses. Its discovery routine uses the narrower conformal case \(R_a=0\), but the checker does not confuse that search restriction with the theorem's scope.
71. Geometry can be proposed from the operators themselves
An operator identity can reveal a useful coordinate geometry without specifying that geometry in advance. Suppose two fitted operators evaluated at an anchor \(u_*\) give an invertible ratio
This ratio is an algebraic object. It need not be a physically available action or a positive map by itself. In mass-first coordinates write
If \(I-A\) is invertible, its fixed center is \(c=(I-A)^{-1}b\). The centered spatial action is \(A\). A candidate invariant metric obeys
With \(A\) fixed, this is a linear equation in the independent entries of symmetric \(H\). No coordinate template appears in that equation.
Proposition 63. Conditional operator-only ellipsoid discovery. Suppose the symmetric solution space of \(A^{\mathsf T}HA=H\) is one-dimensional and contains a positive definite matrix \(H_0\). Let a rank-one fitted branch at a declared observation have a nonzero-mass image, normalized to \(v=(1,x_v)\). If \(x_v\neq c\), define
This gives a unique positive definite metric in that solution space for which \(v\) is on the normalized boundary. If this proposed geometry passes every hypothesis of Proposition 62 for the complete fitted family and its initial state, it is an internally valid probability-domain certificate. The proposal requires no supplied template. Failure at any step withholds this certificate and does not establish model invalidity.
Proof. Invertibility of \(I-A\) gives the unique center. One-dimensionality makes every symmetric invariant form a scalar multiple of \(H_0\). The nonzero displacement has strictly positive squared \(H_0\) length, so the normalization is well-defined, positive and independent of the initially chosen positive scale. The final acceptance test is exactly Proposition 62. Neither invariance under the ratio nor the choice of a boundary image alone is asserted to certify the full family. \(\square\)
The rank-one image proposes a scale. It does not prove that every possible probability domain must place that image on its boundary. Some valid models have many admissible domains. The procedure can also fail because there is no suitable ratio, because invariant forms are not unique, because no rational rank-one anchor is found, or because the chosen ellipsoid fails another branch. These are explicit search outcomes.
71.1 The continuous example no longer needs its template
Use the same exact population moment tables as Section 66. The existing moment learner first reconstructs the two polynomial operators in test coordinates. The new discovery routine receives those operators, their mass row and initial vector. It does not read continuous_template.json or import the population generator.
At \(u_*=0\), the operator ratio gives the required invariant-form equation. A rational rank-one observation supplies a boundary image. The recovered center and metric are
The two metric eigenvalues are \(40/9\) and \(160/9\). For both actions the exact certificate is
Every scalar on the right has an exact rational sum-of-squares representation. Integrated mass is unchanged. The fitted initial state is the ellipsoid center, and the two densities \(13/16\) and \(413/1040\) from Section 66 are unchanged. This is a new certificate construction for the same model, not a change to its predictions.
The exported JSON contains the fitted operators, center, metric, initial state and all witnesses. Independent replay needs only that JSON. The discovery calculation is also checked on three additional rational coefficient families and affine coordinate changes. Those finite examples exercise the stated sector and do not establish unrestricted discovery.
The procedure currently uses three predictive coordinates, a two-dimensional ellipsoid, a finite rational anchor list, a one-dimensional invariant-form space and rational rank-one roots. It proposes only exactly conformal polynomial branches whose scalar nonnegativity can be resolved by its factor-based sum-of-squares routine. General noisy or jointly optimized curved-domain learning remains outside this implementation.
72. Approximate an unknown density without assuming its exact degree
Part V identified a polynomial operator density under an exact degree assumption. A different question is how to approximate a valid density that is not a polynomial at all.
In this section observations lie in \([0,1]\) with Lebesgue reference measure. This is an explicit change from the preceding \([-1,1]\) convention. Under \(v=(u+1)/2\), the probability measure \(du/2\) becomes \(dv\). A Lipschitz constant with respect to \(u\) becomes twice as large with respect to \(v\).
Write the Bernstein basis and the classical Bernstein--Durrmeyer kernel as
Both variables range over \([0,1]\). The kernel is nonnegative and integrates to one in either variable. Let \(B_a(v)\) be any operator-norm integrable density preserving a common probability cone \(C\), with mass functional \(\lambda\) and normalized base \(K\). Define
where
Proposition 64. Valid positive polynomial approximation. Every \(B_a^{[n]}\) preserves \(C\), has degree at most \(n\), and satisfies exact integrated mass conservation. If
for all actions, then
With identical initial states, the corresponding length-\(T\) laws under a common policy differ in total variation by at most
No finite true-degree hypothesis is needed.
Proof. Nonnegative integration preserves the closed convex cone. Since \(\int_0^1b_{n,k}=1/(n+1)\) and \(\sum_k b_{n,k}=1\), the kernel is doubly stochastic. Its integral in \(x\) gives \(\int B_a^{[n]}=\int B_a\), proving exact mass conservation.
For the error, subtract \(B_a(x)\) inside the kernel integral, apply the Lipschitz inequality and use \(\|z\|_C=1\). The integrated discrepancy is at most \(L\mathbb E|V-X|\) under the joint probability density \(D_n(x,v)\). This joint law can be constructed by first choosing \(J\) uniformly from \(\{0,\ldots,n\}\), and then drawing \(X,V\) independently, conditionally on \(J=k\), with density \((n+1)b_{n,k}\).
That conditional beta law has variance
The average of \((k+1)(n-k+1)\) over \(k=0,\ldots,n\) is \((n+2)(n+3)/6\), obtained by summing the first two integer powers. Hence
Cauchy--Schwarz gives the claimed integrated bound. The whole-law conclusion is Proposition 59 applied to the two valid instruments on the same domain. \(\square\)
This construction approximates the operator density itself. It does not average a congruence factor and omit its covariance. Its coefficients are positive operators even when their entries in predictive coordinates are signed.
Operationally, it replaces the raw observation \(v\) by a noisy observation \(x\) drawn from \(D_n(x,v)\). It preserves total mass and underlying unconditioned evolution, not the entire original observation law. That changed law is what the bound measures. A prediction policy must be compared as the same function of the reported history.
72.1 A nonpolynomial density and a failed ordinary truncation
Consider the valid scalar density
It integrates to one, takes values between \(1/2\) and \(3/2\), and has Lipschitz constant two. It is not a single polynomial on \([0,1]\). Proposition 64 nevertheless supplies normalized positive polynomial approximations at every degree, with an explicit approximation budget.
An ordinary orthogonal truncation need not preserve that validity. The valid density \(p(v)=3v^2\) has degree-one \(L^2\) projection
which is negative near zero. The positive construction at degree one instead gives
The two approximations solve different mathematical problems. Neither has permission to call its approximation error zero. The reference checks derive both expressions exactly.
73. Keep approximation bias separate from coefficient error
The positive approximation theorem assumes access to exact integral coefficients. A fitted model introduces another discrepancy. Its coefficient estimates need their own validity constraints.
Proposition 65. Certified coefficient error and approximation bias. Suppose every fitted coefficient \(\widetilde C_{a,n,k}\) preserves the same cone \(C\), and
Then \(\widetilde B_a^{[n]}(x)=\sum_k b_{n,k}(x)\widetilde C_{a,n,k}\) is internally valid. If
its integrated discrepancy from the original Lipschitz density is at most
A known additional model-correction discrepancy is added to this budget, not absorbed into the approximation term.
Proof. Nonnegative Bernstein weights preserve cone positivity, and their integrals prove normalization. For any \(z\in K\), the triangle inequality gives
Add the approximation error from Proposition 64. A separate correction is handled by another application of the triangle inequality. \(\square\)
This is a conditional approximation-and-estimation theorem, not a claim that more coefficients always improve a finite-data estimator. Increasing \(n\) reduces the stated bias bound while changing the number and conditioning of the quantities to be estimated. A cone witness proves that the resulting object is a probability model. It does not prove that the coefficient error budgets cover the generating process.
The moments needed for \(C_{a,n,k}\) are bounded Bernstein-weighted experiments. Their interpretation as full operator coefficients still requires suitable predictive-state coverage or a justified observable core. The next section provides sequential evidence for conditional experiments, not an automatic replacement for those identification assumptions.
74. Evidence from dependent observations
Let \(\mathcal F_{t-1}\) denote the available information immediately before observing the next measurement, including the chosen action, current model and any test selected using past data. A predictable test is a measurable function \(\phi_{t,j}\) of the upcoming measurement, chosen from this information, with values in \([0,1]\). Put
The model prediction is also predictable. It may change as the model learns. The true conditional mean \(m_{t,j}\) may depend on the entire past. No independence, stationarity or known reset state is assumed in this paragraph.
Choose budgets \(\alpha_j>0\) for a finite or countable prospectively specified family of test rules, with \(\sum_j\alpha_j\leq\delta\). Define
Proposition 66. Time-uniform predictable-residual certificate. With probability at least \(1-\delta\), for every time \(t\geq1\) and every registered rule \(j\),
For variables taking values in an interval of width \(c\), multiply the boundary by \(c\). The simultaneous guarantee remains valid when the observation time is chosen from the data.
Proof. For a random variable \(X\in[0,1]\) with conditional mean \(m\), its conditional centered logarithmic moment-generating function has value and first derivative zero at the origin. Its second derivative is the variance under exponential tilting. A variable supported on \([0,1]\) has variance at most \(1/4\), because its squared distance from \(1/2\) is at most \(1/4\). Integrating twice gives
Iterating this conditional inequality gives the same exponential bound for the martingale sum \(\sum_{s\leq t}(X_{s,j}-m_{s,j})\). Optimizing the exponential Markov inequality gives the two-sided bound \(2\exp(-2b^2/t)\). At the displayed \(b_{t,j}\) this is \(\alpha_j/[t(t+1)]\). Sum over \(t\), using \(\sum_{t\geq1}1/[t(t+1)]=1\), and then over \(j\). The resulting event holds at every time, so it also holds at a data-selected time. Scaling the variable proves the range version. \(\square\)
This elementary boundary pays a logarithmic time-union cost. It is not claimed to be optimal. Stronger confidence-sequence constructions are available in the literature [R20]. The proof here is supplied because the exact assumptions and the estimated object matter more than a sharp constant in this foundation.
The certificate estimates accumulated conditional predictive drift, not one fixed unconditional population probability. A prospectively chosen rule may use history, but an analyst cannot inspect an entire data sequence and then retroactively choose an unrestricted test with the same budget. New rules can be registered with fresh budgets and applied to future observations, using elapsed time after registration.
75. Actions need coverage, not only more observations
For intervention evidence use the information \(\mathcal H_{t-1}\) just before choosing the action. Let the action be randomized from a known history-dependent policy \(\pi_t\). The true next-observation kernel under action \(a\) is \(P_t(\cdot\mid\mathcal H_{t-1},a)\). Randomization uses no unrecorded future information. For a predictable test \(\phi_t\in[0,1]\), define
Proposition 67. Covered adaptive-action evidence. If the target action has \(\pi_t(a\mid\mathcal H_{t-1})\geq\pi_{\min}>0\) at every time, then
Consequently Proposition 66 gives a time-uniform certificate for \(\sum_{s\leq t}\mu_s^a\) from \(\sum_{s\leq t}W_s^a\), with its boundary multiplied by \(1/\pi_{\min}\). The target is the average action-conditioned quantity along the encountered contexts, not a reset-population table or the state distribution under a different long-term policy.
Proof. Condition first on the action and then on its observation. Only action \(a\) contributes, with probability \(\pi_t(a\mid\mathcal H_{t-1})\), cancelling the denominator. The range bound follows directly. Apply the scaled conditional moment-generating-function argument to the resulting martingale differences. \(\square\)
This simple bound is conservative. Variance-aware weighting can improve it, but cannot remove the need for information about the action of interest. When \(\pi_{\min}=0\), two models can give exactly the same law under the observed policy and disagree on a never-selected action. For example, both can emit a Bernoulli observation with parameter \(1/2\) under action zero, while their action-one parameters are \(1/10\) and \(9/10\). Repeating action zero cannot distinguish them.
The action randomization and its recorded probabilities are part of the experimental assumption. An observational association with a hidden selection mechanism is not automatically an intervention estimate. Long-horizon off-policy predictions require additional occupancy or trajectory-coverage arguments that this one-step statement does not provide.
76. What a residual can prove about misspecification
A certified operator family can still be the wrong model for the observed process. Predictable residuals offer an explicit way to find that mismatch without pretending that internal validity proves empirical adequacy.
Let \(P_{t,j}\) and \(Q_{t,j}\) denote the true and proposed next-observation laws on the same actual history and action. The proposed law must be defined there. Since \(0\leq\phi_{t,j}\leq1\),
Proposition 68. Sequential evidence gives a lower error certificate. On the simultaneous event of Proposition 66,
Thus crossing the boundary rejects zero conditional drift for that rule, with the registered false-rejection budget. Failing to cross supplies no general upper bound on total variation, no rank certificate and no guarantee for untested actions.
Proof. The triangle inequality gives \(|D_{t,j}|\geq(|R_{t,j}|-b_{t,j})_+\). Also \(|D_{t,j}|\leq\sum_s|m_{s,j}-p_{s,j}|\leq\sum_s\|P_{s,j}-Q_{s,j}\|_{\mathrm{TV}}\). Combine these inequalities. \(\square\)
76.1 Why one aggregate residual can miss conditional errors
The reference path has binary observations. Its true conditional probability of one is \(4/5\) after a one and \(1/5\) after a zero. A deliberately inadequate model predicts \(1/2\) at every step. Its long-run aggregate frequency can be near one half even though its history-conditioned predictions are wrong.
Three tests are registered before the path is generated. One tests the raw observation. A second tests it only after the preceding observation was one. The third applies the same gate but uses the correct conditional predictor. Gates are functions of past data, not selected after seeing the current measurement.
For the supplied 4,096-step seeded path, the aggregate mean residual is approximately \(0.00732\). The gated mean residual is approximately \(0.14917\), while its time-uniform mean boundary is approximately \(0.05114\). The lower certificate for average conditional total variation is approximately \(0.09803\). The correctly predicted gated residual is approximately \(-0.00303\).
These are illustrative path calculations, not estimated coverage rates. The theorem supplies the probability guarantee. The example supplies an explicit reason to test conditional structure rather than relying only on cancellation-prone aggregate residuals.
The exact reference suite also checks a dependent likelihood-ratio construction on every binary path through depth eight. Its expectation under the declared null is exactly one. This is an additional finite martingale calculation, not a replacement for the general argument or a claim that every possible test has been run.
77. Comparing a model with a process outside its representation
The earlier operator-norm law bounds compare two models on one common probability cone. A misspecified generating process may not have a representation in the fitted cone. It is useful to state the more basic comparison at the level of conditional observation laws.
Let \(P\) be any causal process on a standard Borel observation space with finite actions. Let \(Q\) be a proposed process. Both use the same measurable policy. Write \(P_t(\cdot\mid h,a)\) and \(Q_t(\cdot\mid h,a)\) for their next-observation kernels, and define
The proposed kernel must be specified on every history needed for the comparison, including histories impossible under \(Q\) but possible under \(P\). A declared fallback on a model-null history completes a causal kernel without changing \(Q\)'s own law. An implementation that rejects such a history must report that rejection rather than silently claim to have filtered it.
Proposition 69. Observation-law error without a shared state space. With the same initial history, the length-\(T\) joint action-and-observation laws satisfy
An initial-history discrepancy adds its total variation. In particular, uniform bounds \(d_t(h,a)\leq\epsilon_t\) give the upper bound \(\min(1,\sum_{t<T}\epsilon_t)\).
Proof. Introduce hybrid laws that follow \(P\) for the first \(t\) transitions and \(Q\) for the rest, with the common policy throughout. Two successive hybrids have the same \(P\) prefix distribution up to time \(t\) and the same future \(Q\) kernel after the next transition. Conditional on that prefix and its policy-chosen action, changing the next kernel costs \(d_t(h,a)\) in total variation. A common continuation kernel cannot increase total variation, because integration of a test bounded between zero and one produces another such test. Average over the shared prefix, sum the adjacent-hybrid bounds, and cap by one. The same contraction argument handles a different initial law. \(\square\)
This theorem does not make the right-hand side observable without further assumptions. Proposition 68 supplies lower evidence about visited conditional discrepancies, not an upper estimate for every future history. Turning a finite collection of residuals into an upper law guarantee requires additional regularity, coverage or a restricted class with a separating test family.
The distinction prevents a common logical jump: a model can pass several tests and still have a large untested discrepancy. Conversely, a statistically significant residual need not mean that its internal probability algebra is invalid. The model can be a valid distribution and an inadequate approximation at the same time.
78. An integrated theory of construction and correction
Proposition 70. Separate structural acceptance, approximation control and sequential rejection. Consider a fitted polynomial model in a declared predictive coordinate system, and prospectively registered observation tests with total budget \(\delta\).
A successfully replayed quadratic certificate from Proposition 62 establishes internal all-horizon validity. A successful operator-only proposal from Proposition 63 is one way to obtain it, not an assumption that all valid models have such a proposal.
For any generating density on a common probability cone with a justified Lipschitz constant, a positive degree-\(n\) approximation and certified coefficient budgets satisfy Propositions 64--65, hence a fixed-horizon upper law bound. If a generating process lies outside that representation, a separately justified bound on its conditional kernels instead invokes Proposition 69.
For predictable tests on actual dependent observations, Propositions 66 and 68 give simultaneous lower discrepancy certificates, including at data-selected stopping times. Covered randomized actions admit the contextual evidence of Proposition 67. These conclusions coexist without implying that structural acceptance or statistical nonrejection proves the generating process belongs to the model class.
Proof. Each structural statement is a deterministic consequence of its replayed equalities and inequalities. The approximation statements require their declared common-domain, regularity and coefficient conditions, or alternatively the conditional-kernel discrepancy assumptions. The statistical conclusions hold on the simultaneous probability event of Proposition 66. None of these implications uses nonrejection as an upper bound or infers coverage from a finite determinant. Their conjunction therefore retains the separate hypotheses and conclusions exactly as stated. \(\square\)
The central object is now more than a fitted operator. It includes an admissible domain with a replayable witness, an explicit approximation class, and an evidence rule specifying which observations could contradict its predictions. Geometry discovery proposes a representation. Verification checks its algebra. Approximation theory accounts for information lost by a chosen finite description. Sequential tests expose conditional discrepancies in observed signals.
Those roles are related but not interchangeable. A failed search means that this construction did not provide a witness. A failed witness means its submitted identities or signs do not establish the claimed validity. A detected residual means a registered predictive condition is inconsistent with the observations at the stated error budget. These three outcomes should be reported differently.
79. From exact geometry to robustness
Propositions 62–70 distinguish structural acceptance, controlled approximation and sequential evidence. Exact zero-margin preservation does not supply a positive noise radius. Part VII examines that obstruction and gives strict-margin witnesses for selected perturbed families. Statistical nonrejection is not used as an upper accuracy bound.
\newpage
Part VII. Robust curved certificates and noisy geometry
The preceding block discovered a curved probability domain from exact operators. Exact preservation is not the same as stability under uncertainty. In the worked continuous model, every pure boundary state can remain on the boundary. A perturbation too small to be conspicuous in a coefficient table can then move a predicted state outside the probability domain.
This block separates three operations. A numerical procedure proposes a geometry. An exact checker decides whether the submitted polynomial family has a declared margin on that geometry. A separately stated uncertainty set determines which other families share the guarantee. If a repair changes the operators, its prediction cost is retained rather than called an exact coordinate change.
The resulting construction is broader than the preceding exact conformal search. It accepts nonconformal polynomial families and certifies a full coefficient box around a repaired estimate. The reference example uses a deterministic perturbation, not an empirical confidence region. Statistical coverage is a conditional theorem, not a conclusion of that example.
Lorentz-positive maps, robust convex constraints and Bernstein positivity certificates have classical antecedents [R8, R13, R21, R22]. The following arguments are self-contained specializations connecting them to the signal model's representation and evidence contracts. The nine propositions continue the numbering at 71. Their count is not a novelty claim or a count of Lean proofs.
80. A margin that has probability meaning
Let \(\mathcal L_d=\{(m,v):m\geq\|v\|_2\}\), with mass row \(\lambda=(1,0,\ldots,0)\). A mass preserving invertible matrix \(T\) transports a proposed predictive domain to this cone. Its normalized base in the original coordinates is
Write \(F_a(u)=TB_a(u)T^{-1}\) for the transported operator density. Observations lie in \([-1,1]\) with reference measure \(du/2\). All operator densities in this block are polynomial and all equalities involving integrals use that reference measure. For each action define
The subscript \(\rm sp\) denotes spatial coordinates, not a separate measured signal. This margin is expressed before normalization. It simultaneously controls mass orientation and distance from the curved boundary. It is coordinate dependent, so its uncertainty budget must use the same \(T\).
Proposition 71. Robust oriented margin. Suppose \(\sigma_a(u)\geq\kappa_a(u)\) and an uncertain family has transported errors \(\Delta_a(u)\) satisfying, for every normalized input,
If \(\kappa_a(u)\geq e_{a,0}(u)+e_{a,1}(u)\), every uncertain branch preserves \(\mathcal L_d\). If the uncertain family has exact integrated mass and its initial state belongs to \(K_T\), it defines a valid probability model at every finite horizon. With a strict inequality, the residual slack is at least \(\kappa_a-e_{a,0}-e_{a,1}\).
Proof. For \(w=F_a(u)(1,v)\) and \(d=\Delta_a(u)(1,v)\),
The right side is nonnegative under the hypothesis. Homogeneity extends the assertion from normalized inputs to the entire cone. The only zero-mass cone element is zero. Integrated mass conservation and successive cone preservation give nonnegative normalized observation laws by the same induction as in Proposition 52. A positive-density branch may be normalized. A zero-density branch is not assigned a posterior. \(\square\)
This is a uniform deterministic guarantee. The inequalities must cover the whole observation interval and the whole normalized domain. A noise scale estimated from a few residuals does not automatically provide those inequalities. For continuous observations the displayed margin is a density margin relative to the declared measure, not a probability assigned to a single observation value.
81. Why exact boundary preservation can have zero noise radius
An exact certificate can be correct while having no tolerance to unrestricted perturbations. This is a geometrical fact, not a failure of the certificate checker.
Proposition 72. Boundary fragility under mass preserving perturbations. Assume \(d\geq1\). Let \(F\) preserve \(\mathcal L_d\). Suppose a normalized input \(z=(1,v)\) has a boundary image \(Fz=(m,w)\) with \(m=\|w\|_2\). For every \(\epsilon>0\), there is a linear perturbation \(\Delta\) with \(\lambda\Delta=0\) and
such that \((F+\Delta)z\notin\mathcal L_d\). Here the base norm is \(\|(m,w)\|_{\mathcal L_d}=\max(|m|,\|w\|_2)\).
Proof. If \(w\ne0\), put \(e=w/\|w\|_2\). If \(w=0\), choose any spatial unit vector \(e\). Define \(\Delta=\epsilon(0,e)\lambda\). Its mass row vanishes and its value on every normalized input is \((0,\epsilon e)\), of base norm \(\epsilon\). If \(w\ne0\), the new spatial norm is \(m+\epsilon\). If \(w=0\), the original output mass is zero and the new spatial norm is \(\epsilon\). Either image lies outside the cone. \(\square\)
The positive spatial dimension is essential. When \(d=0\), the state is mass alone and \(\lambda\) is the identity on \(\mathbb R\). The condition \(\lambda\Delta=0\) forces \(\Delta=0\), so a perturbation of prescribed norm \(\epsilon>0\) cannot exist. This explicit boundary was missing from the earlier wording. All previously worked examples have two spatial coordinates and are unaffected.
For a continuous operator family, a constant spatial perturbation of this form preserves every integrated mass row and produces the same obstruction at an observation with a boundary image. It need not preserve the algebraic restrictions of a narrower uncertainty class. Perturbing a positive congruence factor and rebuilding its congruence, for example, is different from perturbing an arbitrary predictive matrix.
In the preceding exact continuous example, the direct operator at \(u=0\) is \(I/4\). Every normalized boundary state has a boundary image. There is therefore no positive unrestricted matrix-error radius on that domain, even though the original family is exactly valid. A robust theorem that silently treats its zero slack as positive would be false.
82. An explicit repair and the price of its slack
In Lorentz coordinates let
This reference instrument reports a state-independent density and resets to the center. It is a mathematical model modification, not a new observation. Given an estimate \(\widehat F\), define
Proposition 73. Interior repair with an explicit error budget. Suppose a reference family \(F^*\) is valid on \(\mathcal L_d\), and the combined mass and spatial error of \(\widehat F-F^*\) is at most \(d_a(u)\) on normalized inputs. Then the margin of \(\widetilde F_a(u)\) is at least
An additional uncertainty family with combined output error at most \(\nu_a(u)\) is certified whenever
Assume exact integrated mass for the estimate. If its integrated base-norm discrepancy from \(F^*\) is at most \(\varepsilon\), the corrected discrepancy is at most
When the corrected model is valid on the same cone, with the same initial state, its length-\(T\) law under any common policy differs from the reference law by at most
Proof. The margin of \(\widehat F\) is at least \(-d_a\) by Proposition 71's inequality applied to \(F^*\). Adding the reference changes only output mass, by \(\eta g_a\), while scaling the remaining output by \(1-\eta\). This proves the first bound and the robust implication. Integrated mass is preserved by mixing.
For the error, write
On a normalized input, both \(N\) and \(F^*\) distribute total cone mass one. Their integrated difference has base norm at most two. Integrate the triangle inequality to obtain the stated budget. The final implication is the continuous observation-law bound of Proposition 59. It requires validity of both comparison models on the same cone. \(\square\)
For constant budgets, a sufficient repair weight is \((d+\nu)/(g_{\min}+d)\) when this number is at most one. This is not an optimality claim. The reference example instead fixes a repair weight and verifies the result independently. If an exact model already has a robust structured certificate, adding a reset may be unnecessary and would introduce avoidable bias.
The term \(2\eta\) is deliberately retained. A larger repair creates more room for uncertainty but can change future predictions more. The worst-case law bound can become uninformative at long horizons. Neither positivity nor mixing alone proves long-horizon predictive accuracy.
83. Noisy invariant forms are proposals, not certificates
The exact identity \(A^{\mathsf T}HA=H\) need not survive estimation error. A useful proposal can still be obtained from a least-squares invariant form, but the residual must not be interpreted as a positive-cone guarantee.
Let \(\mathcal M_A\) map symmetric matrices into a Euclidean product space by
Use the Frobenius norm on symmetric matrices and the induced product norm. Assume \(\mathcal M_{A^*}\) is injective and \(\mathcal M_{A^*}(H^*)=(0,1)=b\). This fixes the scale of a one-dimensional invariant-form space.
Proposition 74. A posteriori proposal error. Given \(\widehat H\), let \(r=\mathcal M_{\widehat A}(\widehat H)-b\). Suppose
Then
One may take \(e_M\) from the bound
For fixed centers \(c^*=(I-A^*)^{-1}b^*\) and \(\widehat c=(I-\widehat A)^{-1}\widehat b\), with \(\|\widehat A-A^*\|_2\leq e_A\) and \(\|\widehat b-b^*\|_2\leq e_b\), one has
when the denominator is positive.
Proof. The invariant-form difference satisfies
The least singular value of \(\mathcal M_{A^*}\) is at least that of \(\mathcal M_{\widehat A}\) minus \(e_M\). Apply this lower bound and the triangle inequality. To bound the map perturbation, insert and subtract \(\widehat A^{\mathsf T}HA^*\) and apply the Frobenius product inequality. The trace component is unchanged. Finally,
The same singular-value argument proves the center bound. \(\square\)
The theorem assumes an actual perturbation budget and a corresponding exact reference problem. The example reports numerical residuals but does not claim to have statistically certified these budgets. It uses a trace-normalized least-squares shape, a nearly rank-one image for scale, and rationalized center and factor. It accepts none of those computations as a proof of validity. The next checker uses only the submitted rational geometry and the complete polynomial operators.
A rationalized invertible factor \(L\) specifies \(H=L^{\mathsf T}L\) exactly. Positive definiteness is then an algebraic fact even when the factor was suggested by floating-point calculations. It still does not establish operator invariance.
84. Finite exact checking on a continuum
For three predictive coordinates, the spatial unit circle admits two rational charts. For \(s\in\{1,-1\}\) and \(t\in[-1,1]\), put
The charts cover the two semicircles, including their endpoints. For an operator \(F_a(u)\) and required margin \(\kappa\geq0\), form polynomial vectors and scalars
No square root occurs in these polynomials.
Proposition 75. Rational semicircle and Bernstein certificate. If \(f_{a,s}\) and \(h_{a,s}\) are nonnegative on \([-1,1]^2\) for both charts and every action, then \(\sigma_a(u)\geq\kappa\) for every observation. For rational operators and geometry, nonnegative tensor-Bernstein coefficients on a complete rational subdivision of that square give a finite exact witness. If the actual margin is uniformly strictly greater than \(\kappa\), sufficiently fine dyadic subdivision succeeds.
Proof. The first polynomial inequality establishes output orientation after subtracting \(\kappa d\). The second then gives
Divide by \(d>0\) to obtain the required margin on the circle. For fixed \(u\), the function consisting of affine output mass minus the norm of affine spatial output is concave in the input. Every point in the closed disk is a convex combination of circle points. Concavity therefore extends the lower bound to the disk.
On any rectangle mapped affinely to \([0,1]^2\), the tensor-Bernstein basis functions are nonnegative and sum to one. A polynomial with nonnegative coefficients in that basis is nonnegative throughout the rectangle. If the rectangles cover the square, their inequalities prove the global statement. The reference verifier checks coverage through a complete subdivision tree, not by accepting an arbitrary list of rectangles.
For termination under strict margin, both \(f\) and \(h\) are strictly positive on the compact square. Their minima are positive. For a fixed-degree polynomial on a small cell of side at most \(\delta\), its Taylor expansion at one corner has constant term equal to its value there and nonconstant coefficients bounded uniformly by \(O(\delta)\). The linear change from the fixed-degree power basis to the Bernstein basis is fixed and bounded. Thus every Bernstein coefficient differs from that constant value by \(O(\delta)\), uniformly over cells. At sufficiently small \(\delta\) all coefficients are positive. A finite dyadic subdivision attains that size. \(\square\)
This proof is a specialized use of classical Bernstein range bounds [R22]. It does not claim that a negative Bernstein coefficient proves the polynomial negative. Such a coefficient can merely mean that further subdivision is needed. A negative evaluated corner, by contrast, is an actual counterexample. At a zero-margin boundary, finite termination of this particular witness language is not promised. The exact conformal and matrix certificates of preceding blocks remain useful there.
The two tests are both necessary to this sufficient certificate. The squared inequality alone also admits reversed mass orientation. Both charts are required. A verifier that omits one chart has not covered the circle.
85. A finite uncertainty box and the role of convexity
A robust witness should identify its uncertainty family rather than supply an unexplained tolerance. Suppose the transformed error has polynomial coefficients
All bounds are in the submitted Lorentz coordinates. For three coordinates, normalized inputs have each coordinate of magnitude at most one.
Proposition 76. Coefficient boxes and convex uncertainty. The coefficient box satisfies the output bounds of Proposition 71 with
Only coefficient arrays obeying \(\lambda\int\Delta_a(u)\,du/2=0\) are included in the normalized uncertainty family. A checked margin at least \(e_0+e_1\) certifies every such array.
Alternatively, if finitely many complete polynomial operator families satisfy a common margin \(\kappa\) on the same domain, every convex combination of those families satisfies that margin. Initial states may also range over a certified convex hull.
For a fixed geometry, the constraints \(\sigma_a(u)\geq\kappa\), exact mass, and membership in a convex data-consistency set are convex constraints on the operator coefficients and initial state. Joint optimization of geometry and operators is not asserted to be convex.
Proof. Each output scalar is bounded by \(3r_k\) on a normalized input. Bound the spatial Euclidean norm by the sum of the absolute values of its two coordinates, obtaining \(6r_k\). Sum over coefficients using \(|u|\leq1\). Mass conservation is a separate linear constraint and follows for every family satisfying it.
For convex combinations, output mass is affine and the spatial norm is convex. Hence the margin of the combined output is at least the same combination of margins. A convex cone contains convex combinations of admissible initial states. For fixed geometry and input, the inequality affine mass minus a norm at least \(\kappa\) is convex in the coefficients. Intersect over all inputs and observations, then add the linear mass constraints and the convex data set. \(\square\)
The squared polynomial witness in Proposition 75 is a convenient verification language. Its coefficient constraints are not automatically a convex optimization formulation in unknown operators. Convexity of the underlying admissibility set must not be confused with convexity of every certificate parameterization.
Likewise, a free uncertainty box in all entries contains families that violate integrated mass. The affine restriction is essential. Removing physically valid candidates from a statistical confidence set for some other reason requires its own coverage argument.
86. What a search can guarantee
The numerical proposal is one practical route. There is also a precise conditional search theorem showing why a broad class of strict curved certificates can be found without an exact conformal identity.
Proposition 77. Conditional search for strict rational curved certificates. Consider finitely many rational polynomial operator families in three coordinates, each with exact mass, and finitely many rational normalized initial states. Suppose there exists a real mass preserving invertible \(T_*\) such that every transformed initial state lies strictly inside the Lorentz cone and every transformed operator family has a common positive margin on the observation interval. Then a rational mass preserving invertible \(T\), a positive rational margin, and finite rational Bernstein witnesses exist. A dovetailed enumeration of rational transforms, positive rational margins and subdivision depths eventually finds such a witness. Whenever it returns, the witness is sound for the convex hull of the supplied families.
Proof. Inversion and matrix multiplication are continuous near \(T_*\). Initial membership is strict, and the operator margin is continuous jointly in the transform, the observation, and the normalized input. The latter two range over compact sets and there are finitely many families. Thus the strict inequalities persist in a neighborhood of \(T_*\). Rational mass preserving matrices are dense in that affine matrix space, and invertibility persists in a neighborhood. Choose a rational transform there and a rational positive margin smaller than the remaining minimum. Proposition 75 supplies finite subdivisions because there is still a strict gap above the chosen margin. Enumerating all finite rational descriptions and depths eventually reaches this finite witness. Proposition 76 proves the convex-hull conclusion. \(\square\)
This is a conditional semidecision result, not a polynomial-time complexity theorem. It does not cover every valid model. Initial states on an essential boundary and families with zero margin may violate its assumptions. The implemented procedure uses a finite rationalized proposal and bounded subdivision, not this unbounded enumeration. Failure of that bounded procedure is a withheld certificate rather than a proof that no other domain works.
The theorem's uncertainty families are supplied in fixed predictive coordinates. It does not establish that those families cover an unknown generating process. It demonstrates a route from a strict robust model sector to finite machine-checkable witnesses.
87. Selection from data does not remove the uncertainty quantifier
A domain can be selected after examining fitted operators without a separate statistical penalty for each geometry tried. That statement holds only when the final verification ranges over an already justified confidence set.
Proposition 78. Data-selected geometry with honest coverage. Suppose a random set \(\mathcal U\) of operators and initial states, in a declared fixed predictive coordinate system, contains the true model with probability at least \(1-\delta\). A data-dependent procedure may propose any transform \(T\) and any domain. If its exact verifier proves domain membership and invariance for every member of the same set \(\mathcal U\), then on that coverage event the true model has the verified domain. No union bound over geometry proposals is needed for this deterministic implication.
If original-coordinate perturbations satisfy \(\|\Delta\|_2\leq e\), their combined mass and spatial output error after a known transform \(T\) is at most
The corresponding transformed uncertainty budget must be used even when \(T\) was chosen from the data.
Proof. On the coverage event the true model is a member of the set quantified by the verifier. Its domain properties follow by direct inclusion, regardless of how \(T\) was selected. For the error, a normalized Lorentz input has Euclidean norm at most \(\sqrt2\). Thus the Euclidean norm of \(T\Delta T^{-1}z\) is at most \(\sqrt2\|T\|_2\|T^{-1}\|_2e\). The sum of its mass magnitude and spatial norm is at most \(\sqrt2\) times its Euclidean norm. Multiply the two bounds. \(\square\)
The set must not be shrunk around a convenient candidate without preserving coverage. Moving the center by a reset mixture changes the family being certified. A tube around the repaired family need not contain the original true family. Repeatedly estimating new confidence sets still requires the stage accounting of Proposition 48. Statistical rank coverage and observation identifiability remain assumptions of the upstream procedure.
The final reference certificate explicitly says that no statistical coverage is claimed. Its uncertainty box is a deterministic, exactly described set around the repaired operator family. This supports a reproducible robustness statement without pretending that the example has estimated an unknown noise distribution.
88. A complete perturbed-operator construction
The input to the new proposer contains only two perturbed polynomial operators and the initial vector. A separate fixture generator creates those inputs from the preceding continuous model by adding \(\pm\delta E(u)\), with opposite signs for the two actions and
where
The perturbation preserves integrated mass. It is not a statistical sample. The proposer does not receive the predecessor metric or a hidden generator.
A noisy operator ratio supplies a center and a least-squares shape. A finite observation-anchor list supplies a scale proposal. Rationalizing the center and a triangular metric factor with denominator limit 48 gives
and hence
This is not the exact predecessor metric. The transformed family is not assumed conformal. The raw fitted family actually fails the proposed domain: for the direct action at \(u=-1\) and Lorentz input \((1,1,0)\), its output mass is \(106321/82000>0\), while the output quadratic form is
This is a genuine rejected input state, not a negative Bernstein coefficient of uncertain meaning. The normalized output lies outside the ellipse. No positivity certificate is claimed for the raw fit.
The construction next mixes each operator with the center-reset reference using \(\eta=1/8\). The exact checker proves a margin of at least \(1/16\) for the resulting family. It verifies eight polynomial inequalities across both actions and both semicircle charts. The witness contains 116 subdivision leaves and 1,780 exactly recomputed nonnegative Bernstein coefficients. These are arithmetic certificate sizes, not separate theorems.
The exported family admits mass and spatial perturbation budgets each equal to \(1/128\), leaving certified slack \(3/64\). In particular, it covers every degree-two transformed coefficient perturbation with each entry bounded by \(1/4096\), subject to the exact integrated mass constraints. Proposition 76 verifies that the whole coefficient box lies inside the declared output budgets. No uncertainty vertices or future trajectories are sampled to establish that inclusion.
The certificate is replayed independently from its JSON record. The record contains rational operators, initial state, transform, margin, uncertainty definition and complete subdivision trees. Floating-point proposal diagnostics are stored elsewhere and are not trusted by the replay routine.
89. The prediction cost needs a common comparison domain
A successful discovered domain does not imply that the known reference model has that same domain. To invoke a common-cone law bound, the example separately checks the corrected family on the predecessor's exact ellipse. This is a second certificate used only for comparison. It is not supplied to geometry discovery.
In that comparison domain, the coefficient-wise upper bounds on the raw fitting discrepancy are
Using \(\int_{-1}^1|u|^kdu/2=1/(k+1)\) gives the integrated upper budget \(241/40000\). The same center-reset reference is valid there. Both models have the same initial state.
Proposition 79. Certified repair with a retained comparison budget. For the exact construction of Section 88, independent replay establishes internal validity of the repaired model and its declared coefficient uncertainty family. Independent comparison-domain replay and the coefficient bounds above give
Consequently, for any common policy and any finite horizon \(T\),
This compares the unmodified reference process and the nominal repaired process. It is not a guarantee that the raw noisy model was valid or that the repaired uncertainty family contains the original reference.
Proof. The polynomial witnesses and exact mass identities imply the nominal and robust validity statements by Propositions 71, 75 and 76. The separate comparison witness establishes the required common domain for the nominal repaired and reference operators. The absolute coefficient bounds dominate the integrated raw error by the triangle inequality. Proposition 73 then gives the displayed repair budget and law bound. All hypotheses are either exact identities in the supplied witnesses or direct norm bounds from the declared rational coefficients. \(\square\)
The bound is conservative. It becomes one at horizon eight, so it makes no nontrivial long-horizon accuracy claim. The example's point is the separation of a verified all-horizon probability model from a finite-horizon approximation guarantee. A probability-valid model is not necessarily a close model.
The noisy proposal, the failed raw fit, the repair and the comparison bound all appear in the release records. None is replaced by a more favorable trajectory score. The coefficient box is around the repaired family in the specified coordinates, not an inferred distribution of sensor noise.
90. From deterministic robustness to data confidence
Part VII gives exact continuum certificates for deterministic uncertainty families. The worked repaired model is internally valid at every finite horizon, but its original comparison bound becomes trivial at horizon eight. Part VIII reduces that fixed repair and separately constructs a statistically covered domain for a supplied observable family. Its confidence statement is not transferred to the earlier repaired operators.
\newpage
Part VIII. Data-supported domains and smaller repairs
The preceding block certified a perturbed operator family after mixing it with a center reset. Its all-horizon validity guarantee was exact, but its error bound became trivial at horizon eight. It also kept deterministic coefficient boxes separate from statistical confidence sets. This block addresses both limitations in declared sectors rather than weakening either distinction.
First, the same fitted family and previously proposed geometry admit a smaller, independently checked repair. Second, a self-normalized martingale argument produces operator uncertainty sets from a dependent trajectory when the predictive state is an accessible, fixed observable. A nonlinear three-coordinate example connects those sets to a curved domain without a reset correction. Its geometry and feature dictionary are specified in advance. This is not a solution to discovering an unknown hidden state from the same data.
The sequential-regression result is a classical self-normalized concentration theorem [R23], proved here in the form used by the signal model. The contraction argument uses standard Markov-kernel and total-variation ideas [R24]. The new block's contribution is the explicit sequence of assumptions, finite witnesses, structured uncertainty, observation-measure distinctions and executable examples connecting these tools to the realization calculus. Proposition counts are not historical-novelty counts.
91. A smaller repair can be bracketed rather than guessed
Keep the fitted operators, proposed transform and initial state from Section 88 unchanged. In Lorentz coordinates use the center-reset operator \(N=e_0\lambda\), of constant density one under \(du/2\). Write
Fix a required margin \(0\leq\kappa\leq1\). Feasibility means that every action, observation and normalized input has margin at least \(\kappa\). Integrated mass and initial membership remain separate requirements.
Proposition 80. Convex repair bracketing. The feasible repair weights form a closed interval \([\eta_*,1]\). If a cone-positive functional \(\ell=(1,-v^{\mathsf T})\), with \(\|v\|_2\leq1\), and a normalized input \(z\) give \(c=\ell\widehat F_a(u)z<\kappa\), then every feasible weight satisfies
A complete certificate at a weight \(\eta_U\) gives an upper bound on \(\eta_*\). A negative Bernstein coefficient without a negative evaluated point does not give a lower bound.
Proof. For a fixed input and observation, output mass minus spatial norm is a concave continuous function of the operator, hence of \(\eta\). Its superlevel set is closed and convex. Intersecting these sets preserves those properties. The reset has margin one, so the intersection contains one and is an interval ending there. If \(F_\eta\) has margin at least \(\kappa\), every such functional satisfies \(\ell F_\eta z\geq\kappa\). Since \(\ell Nz=1\), this requires \((1-\eta)c+\eta\geq\kappa\). Here \(1-c>0\), and solving gives the bound. A complete positive witness proves feasibility. A negative coefficient only says that a particular sufficient polynomial bound failed. \(\square\)
91.1 The unchanged predecessor example
Use the previous transform with center \((1/2,1/2)\) and triangular factor
At the direct action, observation \(u=-1\), normalized input \((1,1,0)\) and functional \(\ell=(1,7/25,-24/25)\), exact evaluation gives
For the required margin \(\kappa=1/64\), this proves
The lower endpoint is approximately \(0.02299\). The upper endpoint is not a sampled estimate. A complete exact Bernstein certificate proves that \(\eta=1/32\) works. The weight \(1/64\) fails the requested \(1/64\) margin at the displayed separating functional. This does not by itself prove that its nominal model is probability-invalid.
The new certificate uses 194 leaves and 3,058 recomputed coefficients. A separate comparison-domain certificate uses 185 leaves and 2,833 coefficients. Both retain the same degree-two coefficient box of radius \(1/4096\) in their respective declared Lorentz coordinates, subject to exact integrated mass. The more general mass and spatial output budgets remain \(1/128\) each. Their sum equals the certified margin, so that larger output-error family has a nonnegative lower margin, not a claimed strictly positive residual margin.
The reset weight has fallen from \(1/8\) to \(1/32\), while the deterministic uncertainty box is unchanged. No global optimization over geometry or repair operators is claimed. The current bracket concerns this fixed geometry, this reset path and this required margin.
The comparison domain is checked again before invoking the law theorem. Keeping the same raw discrepancy budget \(241/40000\) from Section 89 gives
Consequently, under a common policy and identical initial states,
The first integer horizon where this upper bound reaches one is now 30 rather than eight. The bound remains conservative. The repaired uncertainty tube is still not a confidence set for the unmodified generating process. This improvement is deterministic certificate refinement, not a new statistical experiment.
92. Confidence from dependent observations
Let \(\mathcal F_{t-1}\) contain all information available before observation \(t\), including its already selected action. Let \(\phi_t\in\mathbb R^d\) be measurable with respect to that information. Assume
for every real \(h\), with \(\sigma>0\), and a declared bound \(\|\theta_*\|_2\leq S\). With ridge \(\gamma>0\), set
Proposition 81. Time-uniform regression confidence. For \(0<\delta<1\), with probability at least \(1-\delta\), simultaneously for every finite \(t\),
No independence of successive observations or of the regressors is required. The conditional mean, conditional noise and predictability hypotheses are required. The confidence set remains valid at a stopping time selected from the data.
Proof. Put \(S_t=\sum_{s\leq t}\phi_s\xi_s\). For every fixed \(q\), the conditional moment bound makes
a nonnegative supermartingale. Integrate against the fixed Gaussian law \(q\sim N(0,(\gamma I)^{-1})\). Conditional Tonelli and completing the square give the nonnegative supermartingale
For completeness, stop at the first crossing of \(1/\delta\), then at any deterministic finite horizon. The supermartingale expectation is at most one, so the crossing probability up to that horizon is at most \(\delta\). Increasing the horizon proves the same bound for any finite crossing. Off this event, taking logarithms bounds \(\|S_t\|_{V_t^{-1}}\) by the square-root term. Finally,
The triangle inequality and \(V_t\succeq\gamma I\) bound the ridge term by \(\sqrt\gamma S\). The simultaneous event already covers every selected finite stopping time. \(\square\)
This is the self-normalized mixture argument in [R23], not a new concentration principle. A bounded response in an interval of length \(b-a\) has conditionally sub-Gaussian centered noise with \(\sigma=(b-a)/2\). One direct proof differentiates its tilted log moment-generating function twice. The second derivative is a variance bounded by \((b-a)^2/4\). Integrating from zero gives the stated moment bound.
The determinant term accounts for the observed design, but does not create excitation. If \(\phi_t\) never probes a direction, the set need not shrink in that direction. Regressors extracted after inspecting the same response are not predictable. An estimated hidden state cannot be substituted for \(\phi_t\) while silently retaining the theorem's true conditional-mean equation.
93. From observable moment regression to operator sets
There is a direct connection to the unnormalized instrument. Assume that the normalized predictive state \(Z_t\) is an exact, fixed, measurable function of the available observation history. For action \(a\), suppose
on realized positive-density branches. Crucially, \(Z_{t+1}\) must also be available from the observations using the specified state map. This is an accessibility assumption, not a consequence of successful regression.
Let \(B_a(u)=\sum_{k=0}^{D}\psi_k(u)C_{a,k}\). Choose functions \(g_j\) so that the known matrix \(W_{jk}=\int g_j(u)\psi_k(u)\,\mu(du)\) is invertible.
Proposition 82. Sequential observable-moment operator sets.
Define \(M_{a,j}=\int g_j(u)B_a(u)\,\mu(du)\). Then
Under the rowwise noise and parameter-bound assumptions of Proposition 81, simultaneous confidence ellipsoids for the rows of \(M_{a,j}\) give a confidence set for every \(C_{a,k}\) by this known linear map. If row \(i\) of \(\widehat M_{a,j}-M_{a,j}\) has \(V_{a,t}\)-norm at most \(\beta_{a,j,i}\), then for any direction \(x\),
Proof. In the conditional integral, multiply the normalized state by the observation density. The denominator cancels, giving \(\int g_j(u)B_a(u)Z_t\,\mu(du)\). The polynomial expansion gives \(M_{a,j}=\sum_kW_{jk}C_{a,k}\) and hence the inverse formula. Apply Proposition 81 to each row, using the predictable regressor \(Z_t\) on visits to action \(a\), or its action-gated version. A union bound over the predeclared finite rows and moments gives simultaneous coverage. Independence between these response coordinates is unnecessary. Weighted Cauchy--Schwarz and the triangle inequality give the displayed directional bound. \(\square\)
Exact mass constraints can be imposed by intersecting the confidence set with the known mass affine subspace. This preserves coverage because the true model is assumed to obey that law. Arbitrarily shrinking a set around a convenient corrected operator does not have the same justification. Data-selected geometry can be checked against the whole resulting set using Proposition 78.
The theorem is conditional on a correct finite observation basis, accessible predictive state and parameter bounds. It does not obtain a hidden predictive state from a single dependent path without additional assumptions. Unknown dictionary discovery, state-estimation error and polynomial misspecification remain distinct problems. The following example makes the state accessible by construction rather than concealing that requirement.
94. A nonlinear observable state with a curved moment domain
For \(u\in[-1,1]\), put
The normalized states are moments of a scalar signal:
Their cone has the finite curved description
Assume that the declared probability measure \(\mu\) on \([-1,1]\) has moments zero at orders one and three, \(1/3\) at order two and \(1/5\) at order four. Define a two-parameter degree-four operator density
Here \(a,b\) are parameters, not action labels. This example is autonomous. Its dependence on the observed signal is nonlinear because the state includes its square.
Proposition 83. Structured curved-domain certificate. The set \(K\) is the convex hull of all \(v(u)\). If \(|a|+|b|\leq1-\rho\) for \(0\leq\rho\leq1\), then \(B_{a,b}(u)\) preserves \(C\) and has exact integrated mass. On a normalized state, its density is
Every positive-density posterior is exactly \(v(u)\). In particular, starting from \(v(s)\) produces the observable Markov law
Its next conditional moments satisfy
Proof. Every distribution on \([-1,1]\) has \(x^2\leq y\leq1\). Conversely, fix \(|x|<1\). A point mass at \(x\) gives \((x,x^2)\), while the endpoint mixture with mean \(x\) gives \((x,1)\). Their convex combinations give every intermediate \(y\). The cases \(x=\pm1\) are immediate. Homogenizing yields the matrix and scalar inequalities, which also force the only mass-zero cone element to be zero.
On \(K\), \(|x|\leq1\), \(0\leq y\leq1\), \(|u|\leq1\) and \(|L_2(u)|\leq1\). Thus the density is at least \(1-|a|-|b|\). Each operator output is that nonnegative scalar times \(v(u)\), so lies in \(C\). The first and second reference moments imply \(\lambda\int B_{a,b}(u)\,\mu(du)=\lambda\). Normalization cancels the scalar density, yielding \(v(u)\). Symmetry through order three gives \(\int uL_2(u)d\mu=0\). The second and fourth moments give \(\int L_2(u)^2d\mu=1/5\). Integrating \(u\) and \(L_2\) against the density proves the last two formulas. \(\square\)
The posterior ray lies on a curved boundary. Arbitrary perturbations of its output coordinates can destroy validity, as in Proposition 72. Parameter uncertainty inside the declared family behaves differently: it changes the likelihood but not the posterior ray. A certified rectangle in \((a,b)\) can therefore be robust even though a free matrix-entry ball has no positive validity radius.
This is an example of uncertainty that preserves the model's algebraic structure. It is not a claim that every real system obeys the two-parameter family. The feature dictionary \(1,u,u^2\), its geometry, the reference measure and the parameter model are declared inputs to learning.
95. A fully rational dependent-data reference experiment
Both the continuous uniform measure \(du/2\) and the following five-node measure satisfy the moment assumptions:
The reference experiment uses the five-node measure, not rounded observations from the continuous model. At each step, the conditional probability of a node is its reference weight times \(p_{a,b}(u\mid s)\). Exact rational probabilities are sampled by integer draws. The learner receives only the node sequence, the measure and the declared model. The true parameters are present only in the separate generator and verification code.
Matching four moments does not identify the two observation laws. The continuous model gives the set of five nodes probability zero, while the discrete model gives it probability one. Their single-observation total variation is therefore one. The continuous theorems and the exact discrete reference calculation are deliberately distinguished.
For a path of \(n\) transitions, use two scalar regressions,
The noise scales follow from response ranges \([-1,1]\) and \([-1/2,1]\). The prior parameter bounds are \(|a|,|b|\leq1\). Each regression uses ridge one and failure probability \(1/40\), giving a total failure budget \(1/20\) for both parameters at every time.
Proposition 84. Exact rational outer confidence intervals. Write \(V=1+\sum\phi_t^2\) and \(\widehat\theta=\sum\phi_tY_t/V\). Let \(k=\lceil\log_2 V\rceil\) and \(h=\lceil\log_2(1/\delta)\rceil\). A rational upper enclosure of
is an upper bound on Proposition 81's scalar radius. Rounding the lower interval endpoint downward and the upper endpoint upward, then intersecting with \([-1,1]\), preserves its conditional coverage. Both intervals cover their true parameters simultaneously at every time with probability at least \(19/20\).
Proof. The positive exponential series gives \(e^{7/10}>1+7/10+(7/10)^2/2+(7/10)^3/6+(7/10)^4/24>2\). Hence \(\log 2<7/10\), so \(\log V+2\log(1/\delta)\leq(7/10)(k+2h)\). Divide the scalar regression bound by \(\sqrt V\). For a nonnegative rational \(x\), the smallest integer \(j\) with \((j/10^{12})^2\geq x\) supplies a rational upper square-root enclosure. Integer square root and an exact final comparison compute it. Outward rounding only enlarges the interval. Intersecting with a true prior constraint preserves coverage. Apply the two-regression union bound on their time-uniform events. \(\square\)
All counts, sums, radius enclosures and interval endpoints in this experiment are rational. No floating-point logarithm is trusted in the confidence arithmetic. An exact replay verifies the arithmetic and the declared protocol. It does not authenticate the observations, prove that a physical process satisfies the model, or measure the procedure's coverage from one pseudorandom run.
95.1 Executed result
The seeded reference path uses \(a=2/5\), \(b=1/5\) and initial observation zero. It was first run to 65,536 transitions and then extended along the same seeded path to 262,144. Time-uniform coverage permits that extension without spending another time-selection budget. This is not a new independent coverage experiment.
At the final checkpoint, outward-rounded intervals are
Every member of this rectangle has density at least
The rounded fitted center is \((409/1000,203/1000)\). It lies inside the rectangle. No reset or other coefficient correction is applied.
At earlier recorded checkpoints the sufficient density-floor test withheld a whole-box certificate. Its first successful recorded checkpoint was 32,768 transitions, with floor \(1/100\). This is not a claim that 32,768 is the smallest possible sample size, or even the first passing individual step. Only the declared checkpoints were evaluated. A withheld whole-box certificate does not declare the nominal fit or true process invalid.
96. Structured uncertainty has an observable error bound
Consider two parameters in a certified family. For normalized \(z=(1,x,y)\), their density difference is
The state-update vector is the same \(v(u)\) in both models after observing \(u\). This does not make the models equivalent. Their likelihoods and hence their predicted observation laws can differ.
Proposition 85. Parameter-to-law error without a state-estimation term. Let
If \(|\Delta a|\leq d_a\) and \(|\Delta b|\leq d_b\), the conditional observation laws on any common normalized state obey
The same bound holds on every realized common observation history. With the same initial observation, complete future records of length \(T\) have total variation at most \(\min(1,T\Delta)\). For the five-node measure,
For \(du/2\), the respective constants are \(1/2\) and \(2/(3\sqrt3)\).
Proof. The first claim follows by integrating \(|\Delta p|/2\), using \(|x|\leq1\) and \(y\leq1\). After the same observed value both models have exactly \(v(u)\), so the conditional bound applies again without accumulated coordinate-state error. Telescoping conditional kernels gives the record bound, equivalently the common-history coupling argument of Proposition 69. Direct summation gives the discrete constants. For the continuous constant \(c_2\), split the integral at \(u=\pm1/\sqrt3\) and integrate the quadratic polynomial with its respective signs. \(\square\)
The use of \(\Delta\) here is an upper bound proved from a confidence set and a model class. It is not the lower discrepancy bound from a rejected residual test in Part VI. A test that does not reject cannot supply the displayed upper bound by itself.
For the final reference rectangle and fitted center,
These bounds compare the fitted model with every member of the rectangle. They apply to the true process on the conditional coverage event. Even when that event fails, the rectangle's internal probability certificate remains an exact statement about its members.
97. Long-horizon marginal control is not full-record control
The previous record bound grows with \(T\). A different question can have a uniform answer: how far apart are the distributions of the signal at a single future time? A lower density bound permits this distinction to be made exactly.
Let \(P,Q\) be Markov kernels on a common state space. Suppose \(P(s,\cdot)\geq\rho\mu(\cdot)\) for every \(s\), where \(\mu\) is one probability measure and \(0<\rho\leq1\). Suppose also \(\sup_s\|P(s,\cdot)-Q(s,\cdot)\|_{\mathrm{TV}}\leq\Delta\).
Proposition 86. Uniform marginal and rolling-window stability. For initial laws at distance \(e_0\), the one-time marginal laws satisfy
For \(k\) consecutive observations starting at time \(t\), their joint laws are at distance at most \(\min(1,e_t+(k-1)\Delta)\). A uniform small marginal bound does not in general bound the total variation of arbitrarily long complete records away from one.
Proof. For \(\rho<1\), write \(P(s,\cdot)=\rho\mu+(1-\rho)R(s,\cdot)\), where \(R\) is a Markov kernel. Applying a Markov kernel cannot increase total variation, by integration of bounded test functions. Thus \(\|\nu P-\widetilde\nu P\|_{\mathrm{TV}}\leq(1-\rho)\|\nu-\widetilde\nu\|_{\mathrm{TV}}\). Adding and subtracting \(\widetilde\nu P\) yields \(e_{t+1}\leq(1-\rho)e_t+\Delta\). Iterate this inequality. The case \(\rho=1\) follows directly. Applying the finite-record telescoping argument from the initial time \(t\) gives the rolling-window statement.
For the final assertion, take two distinct independent Bernoulli laws with parameters strictly between zero and one. Each is a Markov kernel with a positive common minorization. At any time their marginal distance is the fixed parameter difference. For a full record, test whether its empirical mean lies above the midpoint between the parameters. Under each law, Chebyshev's inequality makes the error probability tend to zero as the record length grows. The total variation of the record laws therefore tends to one. \(\square\)
This is the standard minorization/contraction mechanism [R24], with its error accounting made explicit. It is not a stability theorem for arbitrary feedback policies that remember the full history. The reference model is autonomous and its current scalar observation is the Markov state. It also does not bound long products of unnormalized likelihood ratios.
With \(e_0=0\), \(\rho=31/125\) and \(\Delta=2437/90000\), the reference confidence rectangle gives
This is a nontrivial all-time bound on one-time marginal observation distributions. The complete-record bound is still \(\min(1,2437T/90000)\), which becomes trivial at integer horizon 37. The model does not exchange those two different claims. In live filtering after a shared observed value, its state is exactly known and the next-observation law has the separate uniform upper bound \(2437/90000\).
98. A conditional data-to-domain theorem
Proposition 87. Observable structured learning with a covered domain. Suppose the scalar process obeys the conditional law of Proposition 83 with the declared reference measure and fixed parameters satisfying \(|a|,|b|\leq1\). Use the predictable regressions, fixed dictionary and confidence protocol of Proposition 84. At a data-selected finite time, form the full outward-rounded rectangle. If its density floor \(\rho\) is nonnegative, that rectangle has an exact all-horizon probability certificate on \(C\). If \(\rho>0\), every chosen center inside it has the marginal and record bounds of Propositions 85--86 relative to every member.
With probability at least \(19/20\) over the specified stochastic observation process, these implications cover its true parameter at every checked time, including a stopping time. Internal validity of a verified rectangle is deterministic. Statistical coverage remains conditional on the declared process, moment, parameter and observation assumptions.
Proof. Proposition 84 provides one simultaneous event on which both true parameters belong to all returned intervals. Proposition 83 supplies the domain certificate for the entire rectangle whenever its sufficient density floor is nonnegative. The selected center is one of its members, and the maximum distances to the endpoints bound each parameter difference. Proposition 85 gives conditional and record error. If the floor is positive, all true conditional rows have the same lower reference component, so Proposition 86 applies to the scalar Markov process. Selection of the stopping time does not alter the already simultaneous event. None of these implications verifies the assumptions from a single sample path. \(\square\)
No claim of empirical coverage follows from the saved run. It is one seeded finite realization of the exactly declared discrete analogue. The general confidence result is a theorem, with assumptions recorded in the certificate. Its independent checker rederives the sufficient statistics, rational radii, outer endpoints, density floor and error budgets from the saved path. It rejects a changed reference measure, a narrower claimed interval, an altered path checksum or an unsupported claim that the assumptions were independently verified.
This is the point where probability-domain certification and data uncertainty become connected in an executed dependent-signal example. The family is more restrictive than an arbitrary learned operator. That restriction is visible and does useful work: it preserves the nonlinear posterior curve while uncertainty changes its likelihoods.
99. Two distinct completed reference sectors
The smaller-repair construction retains the earlier operator family and geometry. Its coefficient box is deterministic, and its all-horizon validity is distinct from its finite-horizon comparison bound. The dependent-data construction learns two parameters of a different, supplied nonlinear observable family. It certifies a whole confidence rectangle without reset correction and distinguishes next-step, record and marginal bounds. Neither example establishes unrestricted representation discovery.
100. Scope and research boundary
The finite scope contains exact realization and transfer, closure and memory, conditional observable identification, controlled probability instruments, finite polyhedral and curved certificates, explicit approximation budgets, and a restricted dependent-data learning sector. The reference code and its executed examples are catalogued in the Proof and Evidence Register. General nonlinear dictionary discovery, unknown hidden-state learning, empirical sensor modeling and broad deployment remain research directions, not unstated conclusions or release prerequisites.
\newpage
References and source lineage
[K1] Komorebi P. Projector-Resolved Large Algebraic Signal Model. Historical 38-page capstone preserved in Komorebi_PR_LASM_2.3.0_Source.zip, paper/Projector_Resolved_Large_Algebraic_Signal_Model.tex. In particular, Sections 2–7 and the categorical-atlas distinction. Historical experimental claims are not used here.
[K2] Komorebi P. Projector-Resolved Seed Calculus VI, corrected presentation in Paper_VI_With_Visual_Appendix.pdf, Section 8, Theorem 8.4 and the adjacent coarsening criterion and one-packet boundary. This supplement restates and proves the finite result it uses. It does not claim a fresh Lean audit of that paper.
[K3] Komorebi P. PR-LASM 2.3.0 Release Audit, 9 October 2026, and the accompanying source. The source inspected includes exact_atlas.py, lasm_packetizer.py, lasm_model.py, docs/RELEASE_GUIDE.md and docs/claim-boundary.md.
[R1] Ronald R. Coifman and Mladen Victor Wickerhauser. “Entropy-based algorithms for best basis selection.” IEEE Transactions on Information Theory 38(2), 713–718, 1992. DOI: 10.1109/18.119732. Source checked: the authors' Washington University research record. The finite additive tree optimization is classical background.
[R2] MIT OpenCourseWare. 6.241J Dynamic Systems and Control, Chapter 24, “Observability,” Spring 2011, especially Section 24.5 on the largest invariant subspace contained in the observation kernel. Source checked: the MIT course notes, PDF page 6. Used to identify the classical control-theoretic antecedent of Proposition 8.
[R3] Michael Woodward, Yifeng Tian, Arvind Mohan, Yen Ting Lin, Christoph Hader, Hermann Fasel, Misha Chertkov and Daniel Livescu. Data-Driven Mori-Zwanzig: Approaching a Reduced Order Model for Hypersonic Boundary Layer Transition, arXiv:2301.07203, version 1, 2023, introduction and description of resolved, memory and unresolved contributions. Cited for the established projection-memory context, not as evidence for PR-LASM's performance or for any application discussed in that paper. Proposition 10 is proved directly by finite linear recurrence expansion here.
[R4] Peter Buchholz. “Exact and ordinary lumpability in finite Markov chains.” Journal of Applied Probability 31(1), 59–75, 1994. DOI: 10.2307/3215235. Source checked: Cambridge University Press publication record. Used for the classical finite Markov aggregation context.
[R5] Michael L. Littman, Richard S. Sutton and Satinder Singh. “Predictive Representations of State.” Advances in Neural Information Processing Systems 14, NIPS 2001. The paper connects state with predictions of future observations conditional on actions. Cited for conceptual lineage, not as a proof of the finite delay theorem or of this implementation. The original paper and proceedings record are the primary sources. The scanned paper carries all three author names.
[R6] Steven L. Brunton, Bingni W. Brunton, Joshua L. Proctor and J. Nathan Kutz. “Koopman Invariant Subspaces and Finite Linear Representations of Nonlinear Dynamical Systems for Control.” PLOS ONE 11(2), e0150171, 2016. DOI: 10.1371/journal.pone.0150171. Used for the established role and limitations of finite invariant observable families in nonlinear modeling. No general nonlinear feature-discovery guarantee is imported.
[R7] Daniel Hsu, Sham M. Kakade and Tong Zhang. “A Spectral Algorithm for Learning Hidden Markov Models.” Journal of Computer and System Sciences 78(5), 1460–1480, 2012. Primary author manuscript inspected: hmm-jcss.pdf, hosted by Daniel Hsu at Columbia University, especially its abstract and rank assumptions. The result illustrates observable-operator learning under explicit spectral hypotheses. Its statistical guarantees are not claimed for PR-LASM.
[R8] Stephen Boyd, Laurent El Ghaoui, Eric Feron and Venkataramanan Balakrishnan. Linear Matrix Inequalities in System and Control Theory. SIAM Studies in Applied Mathematics 15, 1994. The authors' Stanford publication page identifies the book and provides the primary text. Used to identify the classical common quadratic Lyapunov and matrix-inequality setting. The discrete contraction implication used in Section 19 is proved directly there.
[R9] Byron Boots, Sajid M. Siddiqi and Geoffrey J. Gordon. “Closing the Learning-Planning Loop with Predictive State Representations.” Robotics: Science and Systems VI, 2010. DOI: 10.15607/RSS.2010.VI.036. Primary eight-page proceedings paper inspected at the authors' Carnegie Mellon host, especially Section II on interventional tests and Section III on observable matrices. A journal version appeared in The International Journal of Robotics Research 30(7), 954–966, 2011. Used for the established PSR lineage, not as an assertion that this supplement implements that paper's complete algorithm or planning guarantees. Source: https://www.cs.cmu.edu/~ggordon/boots-siddiqi-gordon-closing-loop-psrs.pdf.
[R10] Han Zhao and Pascal Poupart. “A Sober Look at Spectral Learning.” ICML 2014 Workshop on Spectral Learning, 2014. Primary author-hosted manuscript inspected, including its discussion of negative probabilities in spectral estimates. Used only to identify that practical failure mode. The cone and packet statements in this supplement are proved independently. Source: https://hanzhaoml.github.io/papers/ICML2014/spectral-main.pdf.
[R11] Wassily Hoeffding. “Probability Inequalities for Sums of Bounded Random Variables.” Journal of the American Statistical Association 58(301), 13–30, 1963. DOI: 10.1080/01621459.1963.10500830. Publisher record checked. Proposition 35 supplies a self-contained Bernoulli moment-generating-function proof for the protocol it uses. Source: https://www.tandfonline.com/doi/abs/10.1080/01621459.1963.10500830.
[R7a] Daniel Hsu, Sham M. Kakade and Tong Zhang. “A Spectral Algorithm for Learning Hidden Markov Models.” COLT, 2009. Original conference manuscript inspected at the official Learning Theory proceedings site, particularly Sections 2.3–3.1 on rank, sampling and observable representations. The later journal version is [R7]. No statistical guarantee from either version is imported without its assumptions. Source: https://www.learningtheory.org/colt2009/papers/011.pdf.
[R12] Herbert Jaeger. “Observable operator models for discrete stochastic time series.” Neural Computation 12(6), 1371–1398, 2000. Author-hosted manuscript with the August 2017 erratum. Section 7 gives the polyhedral-cone connection to finite positive realizations and an irrational-rotation obstruction. The controlled rational example in Section 51 is an independently worked adaptation of that established obstruction, not a claim of discovering it. The corrected manuscript and its geometric argument were inspected. Source: https://www.ai.rug.nl/minds/uploads/oom_neco00.pdf.
[R13] Zoltán Horváth, Yunfei Song and Tamás Terlaky. A Novel Unified Approach to Invariance Conditions for a Linear Dynamical System. arXiv:1405.5167, version 2, 5 July 2016. Primary author manuscript, particularly the polyhedral-invariance discussion and use of Farkas' lemma. Used for the classical invariant-cone and optimization setting. The generator and robust-vertex statements used here are proved directly. Source: https://arxiv.org/html/1405.5167v2.
[R14] Marco Fanizza, Josep Lumbreras and Andreas Winter. Quantum theory in finite dimension cannot explain every general process with finite memory. arXiv:2209.11225, version 2, 5 May 2023. Primary manuscript, especially its discussion of finite linear realizations, positive realizations and polyhedral cones. Cited for the distinction between these representation classes. No quantum implementation or quantum advantage is claimed for PR-LASM. The present block uses a real three-dimensional cone construction proved in Section 51. Source: https://arxiv.org/html/2209.11225v2.
[R15] Anna Jenčová. Base norms and discrimination of generalized quantum channels. arXiv:1308.4030, version 1, 19 August 2013. Primary manuscript, Section 1.2, printed page 5, on finite-dimensional proper cones, normalized bases and base norms. That ordered-vector-space definition is used in Section 46. The classical instrument contraction and observation-tree bound needed here are derived in the text, without importing quantum-channel assumptions. Source: https://arxiv.org/pdf/1308.4030.
[R16] John Watrous. The Theory of Quantum Information. Cambridge University Press, 2018. Author-hosted prepublication text, Chapter 2, Theorem 2.22, printed page 82 (PDF page 90), and the trace-preserving characterization in Theorem 2.26. Used for the established Choi and congruence representations of completely positive maps, not to assign quantum physics to the signal model. The relevant theorem page was visually inspected. Source: https://cs.uwaterloo.ca/~watrous/TQI/TQI.pdf. This third-party text is not redistributed in the bundle.
[R17] Pablo A. Parrilo. Algebraic Techniques and Semidefinite Optimization, MIT 6.972, Lecture 10, 14 March 2006. MIT OpenCourseWare notes, pages 10-1 through 10-3, on polynomial Gram matrices, coefficient constraints and interval positivity. Used for the standard sum-of-squares certificate method. Proposition 53 proves only the sufficient matrix-valued certificate used here and does not import scalar completeness as a matrix theorem. The formula pages were visually inspected. Source: https://ocw.mit.edu/courses/6-972-algebraic-techniques-and-semidefinite-optimization-spring-2006/resources/lecture_10/.
[R18] Herbert Jaeger. Modeling and learning continuous-valued stochastic processes with OOMs. GMD technical report, 8 June 2000. Author-hosted report, title and abstract inspected, establishing that continuous-valued observable-operator modeling and learning are prior work. The present polynomial, certificate and binning arguments are proved within this volume. No guarantee from the report is assigned to the reference implementation. Source: https://www.ai.rug.nl/minds/uploads/jaeger.00.tr.contoom.pdf.
[R19] Ulrich Abel, Ana Maria Acu, Margareta Heilmann and Ioan Raşa. A representation for the integral kernel of the composition of multivariate Bernstein--Durrmeyer operators. arXiv:2512.17297v1, 19 December 2025. The authors' primary manuscript, Introduction, gives the classical univariate Bernstein--Durrmeyer kernel used in Section 72. Its composition theorems are not needed or claimed here. Proposition 64 derives the operator-valued positivity and error statements directly. Source: https://arxiv.org/html/2512.17297v1.
[R20] Steven R. Howard, Aaditya Ramdas, Jon McAuliffe and Jasjeet Sekhon. Time-uniform, nonparametric, nonasymptotic confidence sequences. Annals of Statistics 49(2), 1055--1080, 2021. DOI: 10.1214/20-AOS1991. The primary arXiv record and the conditional exponential-supermartingale discussion in the v4 HTML manuscript were inspected. Proposition 66 proves a simpler union-bound confidence sequence, not the paper's sharper rates. Source: https://arxiv.org/abs/1810.08240 and https://arxiv.org/html/1810.08240v4.
[R21] Roland Hildebrand. An LMI description for the cone of Lorentz-positive maps. Author manuscript, 17 January 2006. Primary manuscript inspected, especially the introduction and definitions on pages 1--2. It supplies classical context for Lorentz-positive maps and the existence of convex descriptions. The present margin and semicircle certificate arguments are proved directly and do not claim to implement the manuscript's full matrix-inequality construction. Source: https://optimization-online.org/wp-content/uploads/2005/12/1260.pdf.
[R22] Tareq Hamadneh, Hassan Al-Zoubi, Hamza Alzaareer and Rafael Wisniewski. Optimization and Positivity Certificates of Rational Functions using Bernstein Form. arXiv:1906.11037, 2019. Primary manuscript inspected, including its Bernstein range enclosure discussion on printed pages 3--4. It supplies established positivity-certificate context. Proposition 75 gives its own tensor-product rectangle argument and strict-subdivision proof. Source: https://arxiv.org/abs/1906.11037.
Inherited dependencies. Proposition 52 defines continuous instrument semantics. Proposition 59 supplies the common-domain continuous-law bound. Propositions 62--63 distinguish exact ellipsoid proposals from verification. Proposition 48 supplies repeated confidence-set accounting. The new block restates the specific repair and uncertainty inequalities it uses. The cited earlier results are present in this cumulative volume.
[R23] Yasin Abbasi-Yadkori, David Pál and Csaba Szepesvári. Online Least Squares Estimation with Self-Normalized Processes: An Application to Bandit Problems. arXiv:1102.2670, version 1, 14 February 2011. The primary manuscript's Section 2, Theorem 3, and Section 5 supply the classical martingale-mixture and adaptive-regression context. Proposition 81 includes its own proof, retains conditional noise and design assumptions, and does not claim a new concentration inequality. Source checked: https://arxiv.org/pdf/1102.2670.
[R24] David A. Levin and Yuval Peres, with contributions by Elizabeth L. Wilmer. Markov Chains and Mixing Times, second edition, 2017. Chapters 4 and 5 develop total variation and coupling. The present minorization argument is proved directly and does not import a finite-state assumption into its general kernel statement. The worked reference chain is finite. Source checked: the author-hosted second-edition text, https://pages.uoregon.edu/dlevin/MARKOV/mcmt2e.pdf.
Internal dependencies. Part VII, especially Propositions 71, 73, 75, 76 and 79, supplies the robust checker and predecessor repair comparison. Propositions 59 and 69 supply the continuous and conditional-kernel full-law comparisons. Their mathematical source passages are retained in this cumulative volume; prior archive identifiers remain in the provenance records.