Reading settings

Settings and reading history are stored only in this browser.

English中文0%

Select a passage or formula to leave a suggestion here. Saved suggestions remain private until reviewed.

Unit 3, continued: delaying the projection at a fixed depth of use

This part develops D Lemma 5.5, pp.35–37, equations (5.18)–(5.22). Fix one triple from the predictor family in the preceding part. The aim is to replace the actual function Fa,sF_{a,s} by the finer projection Pa,s+K0Fa,sP_{a,s+K_0}F_{a,s}. Every refinement acts on the same Fa,sF_{a,s}.

The delay error has two sources. Variation of the smooth test function on refined intervals contributes 2−K02^{-K_0}. The predictor's private coordinate contributes ϵ\sqrt\epsilon, and crossing K0K_0 projection levels gives K0ϵK_0\sqrt\epsilon. We derive these quantities from orthogonality and a finite telescoping identity.

1. The full assumptions and conclusion

Fix the two complementary positions a,ba,b and an index α\alpha from Corollary 5.4 in the preceding part. Abbreviate

ψa=ψa,α,ψb=ψb,α,λ=λα. \psi_a=\psi_{a,\alpha},\quad \psi_b=\psi_{b,\alpha},\quad \lambda=\lambda_\alpha.

The two functions are fixed and satisfy ∥ψa∥∞,∥ψb∥∞≤1\|\psi_a\|_\infty,\|\psi_b\|_\infty\le1. Inside a dyadic subroot K⊆J′K\subseteq J', choose a finite family U\mathcal U of intervals of use, retaining only those with

∣λ∣∣I∣2≤Dpred. |\lambda||I|^2\le D_{\mathrm{pred}}.

Position aa has inherited orthogonal projections Pa,hP_{a,h} on complete partitions, with a fixed column-function constant DhD_h. Here hh is the projection depth. The symbol DhD_h denotes one numerical bound for the column family in the complementary position, rather than a depth-indexed sequence of bounds.

Fix an integer K0≥0K_0\ge0 and a number 0<ϵ≤10<\epsilon\le1. Suppose that on every interval of use and all its further descendants, the projection family contains the same predictor column

Φv=(ψa‾,ϵ ev).(Predictor-column) \Phi_v=(\overline{\psi_a},\sqrt\epsilon\,e_v). \tag{Predictor-column}

The predictor column is included in the column-function bound, so ϵ≥Dh−1\sqrt\epsilon\ge D_h^{-1}. Further columns may be added, subject to the same prescribed bound DhD_h and with distinct private directions. The predictor column vv retains its own direction eve_v.

The original root is JJ, and every depth s(I)=log⁡2(∣J∣/∣I∣)s(I)=\log_2(|J|/|I|) keeps its original index. The sequence Fa,s=(fa,s,0)F_{a,s}=(f_{a,s},0) is the same complete normalized sequence restricted to KK. In position bb, the function zb,Iz_{b,I} is a sum of components from the three permitted opposing classes in the preceding part, with parameters Dh,LD_h,L.

The conclusion is

∑I∈U∣I∣ ∣TIα(Fa,s(I),zb,I)−TIα(Pa,s(I)+K0Fa,s(I),zb,I)∣≤C(2−K0+(1+K0)ϵ)poly⁡(Dpred,log⁡(2Dh),1+L)∣K∣.(Delayed-error) \begin{aligned} \sum_{I\in\mathcal U}|I|\, \bigl|\mathcal T_I^\alpha(F_{a,s(I)},z_{b,I}) -\mathcal T_I^\alpha(P_{a,s(I)+K_0}F_{a,s(I)},z_{b,I})\bigr| \\ \le C\bigl(2^{-K_0}+(1+K_0)\sqrt\epsilon\bigr) \operatorname{poly}\bigl(D_{\mathrm{pred}},\log(2D_h),1+L\bigr)|K|. \end{aligned} \tag{Delayed-error}

where the forms use only the physical scalar components of the functions.

2. The factor ϵ\epsilon in residual means and ϵ\sqrt\epsilon in private coordinates

Fix a depth of use ss. As the projection depth h≥sh\ge s varies, keep Fa,sF_{a,s} fixed. Write

rh[s]=(fa,s−(Pa,hFa,s)phys)ψa,μh[s]=Ehrh[s], r_h^{[s]} =\bigl(f_{a,s}-(P_{a,h}F_{a,s})_{\mathrm{phys}}\bigr)\psi_a, \qquad \mu_h^{[s]}=\mathbb E_h r_h^{[s]},

Here Eh\mathbb E_h is conditional expectation on the complete partition of KK at depth hh.

Take a small interval at depth hh, denoted by QQ, inside an interval of use. Let the coefficient of the predictor column in Pa,hFa,sP_{a,h}F_{a,s} be bh,v[s]b_{h,v}^{[s]}. On QQ, the residual is orthogonal to (Predictor-column). Since the inner product is linear in the first variable,

0=1∣Q∣⟨Fa,s−Pa,hFa,s,Φv⟩L2(Q)=1∣Q∣∫Qrh[s]−ϵbh,v[s]. 0=\frac1{|Q|}\langle F_{a,s}-P_{a,h}F_{a,s},\Phi_v\rangle_{L^2(Q)} =\frac1{|Q|}\int_Qr_h^{[s]}-\epsilon b_{h,v}^{[s]}.

Hence

μh[s]=ϵbh,v[s].(Mean-coefficient) \mu_h^{[s]}=\epsilon b_{h,v}^{[s]}. \tag{Mean-coefficient}

Set Δh=Pa,h+1−Pa,h\Delta_h=P_{a,h+1}-P_{a,h}. The projection's private component in direction vv is ϵ bh,v[s]\sqrt\epsilon\,b_{h,v}^{[s]}, so on the finer intervals

μh+1[s]−μh[s]=ϵ (ΔhFa,s)v,(Mean-increment) \mu_{h+1}^{[s]}-\mu_h^{[s]} =\sqrt\epsilon\,(\Delta_hF_{a,s})_v, \tag{Mean-increment}

and

rh[s]−rh+1[s]=ψa(ΔhFa,s)phys.(Residual-increment) r_h^{[s]}-r_{h+1}^{[s]} =\psi_a(\Delta_hF_{a,s})_{\mathrm{phys}}. \tag{Residual-increment}

These three identities remain valid after further columns are added, since the same predictor column and its private direction are retained. The physical component of the column must be ψa‾\overline{\psi_a}, while the residual is multiplied by ψa\psi_a. This conjugation can be checked directly on p.35 of the source PDF.

3. A finite identity for residual moments

Take an interval of use II at depth ss, and write h∗=s+K0h_*=s+K_0. Let ωI\omega_I be a smooth function satisfying

∥ωI∥∞+∣I∣∥ωI′∥∞≤Aω. \|\omega_I\|_\infty+|I|\|\omega_I'\|_\infty\le A_\omega.

For now, its integral need not be zero. For every finite n>h∗n>h_*, the source equation (5.20) reads

∫Irh∗[s]ωI=∫Iμh∗[s]ωI+∑h=h∗n−1∫I(rh[s]−rh+1[s])(ωI−EhωI)+∑h=h∗n−1∫I(μh+1[s]−μh[s])(Eh+1ωI−EhωI)+∫Irn[s](ωI−EnωI).(Telescope) \begin{aligned} \int_Ir_{h_*}^{[s]}\omega_I ={}&\int_I\mu_{h_*}^{[s]}\omega_I\\ &+\sum_{h=h_*}^{n-1}\int_I (r_h^{[s]}-r_{h+1}^{[s]})(\omega_I-\mathbb E_h\omega_I)\\ &+\sum_{h=h_*}^{n-1}\int_I (\mu_{h+1}^{[s]}-\mu_h^{[s]}) (\mathbb E_{h+1}\omega_I-\mathbb E_h\omega_I)\\ &+\int_Ir_n^{[s]}(\omega_I-\mathbb E_n\omega_I). \end{aligned} \tag{Telescope}

To verify it, define

Th=∫Irh[s](ωI−EhωI). T_h=\int_Ir_h^{[s]}(\omega_I-\mathbb E_h\omega_I).

Direct subtraction gives

Th−Th+1=∫I(rh[s]−rh+1[s])(ωI−EhωI)+∫Irh+1[s](Eh+1ωI−EhωI). \begin{aligned} T_h-T_{h+1} ={}&\int_I(r_h^{[s]}-r_{h+1}^{[s]}) (\omega_I-\mathbb E_h\omega_I)\\ &+\int_Ir_{h+1}^{[s]} (\mathbb E_{h+1}\omega_I-\mathbb E_h\omega_I). \end{aligned}

The test function in the second term is measurable for the depth-h+1h+1 partition, so rh+1[s]r_{h+1}^{[s]} may be replaced by μh+1[s]\mu_{h+1}^{[s]}. Moreover, μh[s]\mu_h^{[s]} is measurable for the depth-hh partition, so

∫Iμh[s](Eh+1ωI−EhωI)=0. \int_I\mu_h^{[s]} (\mathbb E_{h+1}\omega_I-\mathbb E_h\omega_I)=0.

These are the two terms in the telescoping sum. Finally,

∫Irh∗[s]ωI=∫Iμh∗[s]ωI+Th∗, \int_Ir_{h_*}^{[s]}\omega_I =\int_I\mu_{h_*}^{[s]}\omega_I+T_{h_*},

and summing Th−Th+1T_h-T_{h+1} over finitely many levels proves (Telescope). All integrals are complex bilinear pairings; these conditional-expectation identities require no additional conjugation.

4. The terminal term vanishes without convergence of projections to the original function

For h≥sh\ge s, the complete projection partition refines the interval of use II. Thus Pa,hP_{a,h} commutes with 1I\mathbf1_I. The contraction property of projection and the bound ∣ψa∣≤1|\psi_a|\le1 give

∥rh[s]∥L2(I)≤2∥Fa,s∥L2(I). \|r_h^{[s]}\|_{L^2(I)} \le2\|F_{a,s}\|_{L^2(I)}.

Smoothness gives

∥ωI−EnωI∥L2(I)≤CAω∣I∣1/22s−n. \|\omega_I-\mathbb E_n\omega_I\|_{L^2(I)} \le C A_\omega |I|^{1/2}2^{s-n}.

Cauchy–Schwarz therefore makes the terminal term in (Telescope) tend to zero.

If the original projections are defined only up to a finite depth, retain all existing columns and locally reproject on the finer intervals required. The column-function constant is unchanged. Apply the identity for each finite nn, then let n→∞n\to\infty. The estimate does not assume Pa,nFa,s→Fa,sP_{a,n}F_{a,s}\to F_{a,s}.

5. Two projection-energy estimates across depths of use

These are D equation (5.21), p.36, following from the previously developed Corollary 4.5.

For each fixed integer offset ℓ\ell, the intervals at depth of use ss are disjoint, so

∑s∑I∈Us(I)=s∥1IΔs+ℓFa,s∥22≤C(log⁡(2Dh))2∣K∣.(Offset-energy) \sum_s\sum_{\substack{I\in\mathcal U\\s(I)=s}} \|\mathbf1_I\Delta_{s+\ell}F_{a,s}\|_2^2 \le C\bigl(\log(2D_h)\bigr)^2|K|. \tag{Offset-energy}

Only valid projection depths are included. Similarly, for an integer Lproj≥0L_{\mathrm{proj}}\ge0,

∑s∑I∈Us(I)=s∥1I(Pa,s+Lproj−Pa,s)Fa,s∥22≤C(1+Lproj)2(log⁡(2Dh))2∣K∣.(Rolling-energy) \sum_s\sum_{\substack{I\in\mathcal U\\s(I)=s}} \|\mathbf1_I(P_{a,s+L_{\mathrm{proj}}}-P_{a,s})F_{a,s}\|_2^2 \le C(1+L_{\mathrm{proj}})^2 \bigl(\log(2D_h)\bigr)^2|K|. \tag{Rolling-energy}

Both estimates follow by reapplying the full row estimate on the subroot KK. However fine the projection partition becomes, its operand remains Fa,sF_{a,s} from the depth of use.

6. Error rows for the smooth mean-zero part

First replace e(λt2)\mathrm e(\lambda t^2) by one in the form. As in the preceding part, the resulting kernel has a tensor expansion in smooth, interior-supported functions with exactly zero mean in each variable. Fix a tensor index and apply (Telescope) to the test function in position aa, denoted by ωI\omega_I.

6.1 The two remainder sums contribute 2−K02^{-K_0}

For h=s+ℓh=s+\ell, smoothness gives

∥ωI−EhωI∥∞+∥Eh+1ωI−EhωI∥∞≤CAω2−ℓ. \|\omega_I-\mathbb E_h\omega_I\|_\infty +\|\mathbb E_{h+1}\omega_I-\mathbb E_h\omega_I\|_\infty \le C A_\omega2^{-\ell}.

The normalized test in the first remainder sum satisfies

∣I∣−1/2∣∫Iψa(Δs+ℓFa,s)phys(ωI−Es+ℓωI)∣≤CAω2−ℓ∥1IΔs+ℓFa,s∥2. |I|^{-1/2} \left|\int_I\psi_a(\Delta_{s+\ell}F_{a,s})_{\mathrm{phys}} (\omega_I-\mathbb E_{s+\ell}\omega_I)\right| \le C A_\omega2^{-\ell} \|\mathbf1_I\Delta_{s+\ell}F_{a,s}\|_2.

By (Mean-increment), the second remainder sum has an additional factor ϵ≤1\sqrt\epsilon\le1 and obeys the same bound.

For fixed ℓ\ell, apply (Offset-energy), then use Minkowski's inequality in the interval-indexed Hilbert direct sum to sum over ℓ≥K0\ell\ge K_0. Since

∑ℓ≥K02−ℓ=21−K0, \sum_{\ell\ge K_0}2^{-\ell}=2^{1-K_0},

the combined output norm of the two remainder sums is at most

CAω2−K0log⁡(2Dh)∣K∣1/2.(Smooth-remainders) C A_\omega2^{-K_0}\log(2D_h)|K|^{1/2}. \tag{Smooth-remainders}

6.2 The initial mean contributes K0ϵK_0\sqrt\epsilon

Here ωI\omega_I has integral zero, while μs[s]\mu_s^{[s]} is constant on II. Hence

∫Iμh∗[s]ωI=∫I(μh∗[s]−μs[s])ωI. \int_I\mu_{h_*}^{[s]}\omega_I =\int_I(\mu_{h_*}^{[s]}-\mu_s^{[s]})\omega_I.

Telescoping (Mean-increment) gives

μh∗[s]−μs[s]=ϵ ((Pa,s+K0−Pa,s)Fa,s)v. \mu_{h_*}^{[s]}-\mu_s^{[s]} =\sqrt\epsilon\,((P_{a,s+K_0}-P_{a,s})F_{a,s})_v.

Using the exact K0K_0 form of Corollary 4.5 and the same normalized test, the initial term has Hilbert direct-sum norm at most

CAωK0ϵlog⁡(2Dh)∣K∣1/2.(Zero-mean-initial) C A_\omega K_0\sqrt\epsilon\log(2D_h)|K|^{1/2}. \tag{Zero-mean-initial}

For K0=0K_0=0, this term is exactly zero.

7. Pairing with the row in the other position

At a fixed tensor index, ∣I∣|I| times the form equals the normalized residual moment in position aa times the normalized row value in position bb, multiplied by the fixed coordinate Jacobian factor and the tensor coefficient. Apply Cauchy–Schwarz over all intervals of use. Position aa is controlled by (Smooth-remainders) and (Zero-mean-initial). Position bb has three cases.

For an original function sequence, the full multiplier is the fixed function ψb\psi_b, with one smooth mean-zero row per interval of use. Its ordinary Bessel bound follows directly from the smooth-row proof in Unit 3: decompose the test function into martingale differences. The analysis operator for the hh-th part has norm at most 2Aω2−h2A_\omega2^{-h}; summing over h≥1h\ge1 gives at most 2Aω2A_\omega. Thus the single-multiplier case of D Lemma 4.6 applies and gives the square-sum row estimate for varying functions.

For a current projected function, first apply the H Compression operator statement recorded in the preceding part to this fixed single-multiplier row family, then apply D Lemma 4.4. Both the support and observation partitions are the complete partition at the current depth of use, with the same index ss. The cost is a fixed polynomial in log⁡(2Dh)\log(2D_h).

For an energy-type function Zb,IZ_{b,I}, Cauchy–Schwarz for one row gives

∣I∣−1/2∣∫IψbωIZb,I∣≤Aω∥Zb,I∥L2(I). |I|^{-1/2}\left|\int_I\psi_b\omega_I Z_{b,I}\right| \le A_\omega\|Z_{b,I}\|_{L^2(I)}.

The given total square-energy bound then bounds the output norm by AωL∣K∣1/2A_\omega L|K|^{1/2}.

The derivative bounds grow only polynomially in the tensor index, so the rapidly decaying tensor coefficients allow absolute summation of these estimates. The mean-zero part is therefore bounded by

C(2−K0+K0ϵ)poly⁡(log⁡(2Dh),1+L)∣K∣ C\bigl(2^{-K_0}+K_0\sqrt\epsilon\bigr) \operatorname{poly}\bigl(\log(2D_h),1+L\bigr)|K|

The H Compression result used here remains an explicit external assumption.

8. The quadratic-phase difference need not have zero mean

It remains to treat the kernel

wI(x,t)(e(λt2)−1). w_I(x,t)\bigl(\mathrm e(\lambda t^2)-1\bigr).

Fix a depth of use ss and set

xs=∣λ∣∣I∣2, x_s=|\lambda||I|^2,

This depends only on the depth, not on the interval of use at that depth. Only depths with xs≤Dpredx_s\le D_{\mathrm{pred}} are retained.

After rescaling to a unit interval, every required fixed-order derivative bound for this difference kernel contains the factor

xspoly⁡(Dpred). x_s\operatorname{poly}(D_{\mathrm{pred}}).

The one-variable functions in its ordinary smooth tensor expansion may now have nonzero mean. They still have polynomial derivative costs and rapidly decaying coefficients.

8.1 Ordinary projection orthogonality at fixed ss

The operand stays fixed at Fa,sF_{a,s} throughout refinement. Ordinary orthogonality of the nested projections therefore gives

∑h≥s+K0∥ΔhFa,s∥22≤∥1KFa,s∥22. \sum_{h\ge s+K_0}\|\Delta_hF_{a,s}\|_2^2 \le\|\mathbf1_KF_{a,s}\|_2^2.

This is a square sum with ss fixed and hh varying. It differs from the varying-function estimates across depths of use in Section 5.

The intervals of use at one depth are disjoint. Apply the smooth decay 2s−h2^{s-h} to the two remainder sums in (Telescope), then Cauchy–Schwarz in hh. The output norm is at most

CAω2−K0∥1KFa,s∥2. C A_\omega2^{-K_0}\|\mathbf1_KF_{a,s}\|_2.

8.2 The initial term with nonzero mean

We can no longer subtract μs[s]\mu_s^{[s]}. Instead, (Mean-coefficient) gives directly

μh∗[s]=ϵ (Pa,h∗Fa,s)v. \mu_{h_*}^{[s]} =\sqrt\epsilon\,(P_{a,h_*}F_{a,s})_v.

Disjointness at one depth and the contraction property of projection bound the initial test's output norm by

CAωϵ ∥1KFa,s∥2. C A_\omega\sqrt\epsilon\,\|\mathbf1_KF_{a,s}\|_2.

The output norm of the test in the other position at this fixed depth is at most CAω(1+L)∣K∣1/2C A_\omega(1+L)|K|^{1/2}, by the original function bound, local projection contraction, or the given energy bound. Normalization of the original function gives ∥1KFa,s∥2≤∣K∣1/2\|\mathbf1_KF_{a,s}\|_2\le|K|^{1/2}. The contribution at this depth is therefore at most

xspoly⁡(Dpred)(2−K0+ϵ)(1+L)∣K∣. x_s\operatorname{poly}(D_{\mathrm{pred}}) \bigl(2^{-K_0}+\sqrt\epsilon\bigr)(1+L)|K|.

8.3 A fixed unscaled λ\lambda makes the scale sum finite

For λ=0\lambda=0, the difference kernel vanishes. Otherwise, xsx_s decreases by the factor 1/41/4 as depth increases, so

∑s:xs≤Dpredxs≤CDpred. \sum_{s:x_s\le D_{\mathrm{pred}}}x_s\le C D_{\mathrm{pred}}.

The same bound holds for every selected finite subset of depths. Summing the estimate from Section 8.2 controls the entire quadratic-phase difference. Adding the bound from Section 7 proves (Delayed-error).

This requires one fixed unscaled phase parameter λ\lambda. Choosing unrelated values λs\lambda_s at different depths would invalidate the geometric-series argument.

9. Successive replacements stay within the permitted classes

This is the final paragraph of D Lemma 5.5, p.37.

If position bb has already undergone delayed projection, then on an interval of use II with s=s(I)s=s(I) we have the exact identity

Pb,s+K0Fb,s=Pb,sFb,s+(Pb,s+K0−Pb,s)Fb,s. P_{b,s+K_0}F_{b,s} =P_{b,s}F_{b,s} +(P_{b,s+K_0}-P_{b,s})F_{b,s}.

The first term is a current projected component. Restrict the entire second term to II. Equation (Rolling-energy) gives

∑I∈U∥1I(Pb,s(I)+K0−Pb,s(I))Fb,s(I)∥22≤C(1+K0)2(log⁡(2Dh))2∣K∣. \sum_{I\in\mathcal U} \left\|\mathbf1_I(P_{b,s(I)+K_0}-P_{b,s(I)})F_{b,s(I)}\right\|_2^2 \le C(1+K_0)^2\bigl(\log(2D_h)\bigr)^2|K|.

It is therefore a permitted energy-type component with parameter

L′=C(1+K0)log⁡(2Dh) L'=C(1+K_0)\log(2D_h)

Subtracting an old current projection leaves exactly three components with coefficients 1,1,−11,1,-1, as allowed by the preceding definition.

The entire delayed difference must be treated as one energy component. Splitting it into K0K_0 opposing components would no longer justify a bound of three on the number of components. The exact decomposition above allows the same error estimate to be used for the two successive replacements.

Source coverage and subsequent use

  • D p.35, equations (5.18)–(5.20): the full error statement, the conjugated predictor column, residual moments, and the finite telescoping identity.
  • D p.36, equations (5.21)–(5.22): varying-function energy at a fixed offset, the mean-zero initial term, and the quadratic difference-kernel estimate at a fixed depth of use.
  • D p.37: geometric summation with a fixed phase and the three-component decomposition for successive replacements.

The next step, D §5.6, must choose B,K0,ϵB,K_0,\epsilon, sum over all pairs of functions in the fixed-atom positions, and construct additional high-precision projections whose final column-family logarithmic size is compatible with these errors. The present estimate says that once the parameters and an actual inherited projection family satisfy the stated assumptions, delayed replacement has the uniform error displayed above.