Research / Mathematics / Information Science · Number Geometry
Arithmetic Sensing III — Optimized Quadrature and Unavoidable Aliases
Positive cosine windows, fixed-grid aliases, and a 117-fold complete-tail certificate
Originating question and research environment: TGN's human founder
Technical continuation · theorem, computation and hypothesis separated
Plain-text Markdown · Structured JSON · Number Geometry index · Arithmetic Sensing I–II
Why this continuation matters
Arithmetic Sensing III — Optimized Quadrature and Unavoidable Aliases
Arithmetic Sensing II used a Hann-weighted midpoint grid. This continuation asks whether the window can be designed for logarithmic integer frequencies themselves. A positive cosine-series family gives a globally solved finite linear program within its declared family and truncation, followed by an independent exact-kernel and million-term tail audit. The same analysis proves a hard limitation: weighting a fixed uniform midpoint grid can suppress pre-alias leakage but cannot remove exact aliases.
1. A positive arithmetic window family
Exact statement or rational verificationFor midpoint times t_j=(j+1/2)T/m, the order-H cosine family uses w_j(c)=m^-1[1+2 sum(c_r cos(2 pi r t_j/T))]. Midpoint cosine sums preserve total mass one. Requiring the bracketed density to be nonnegative makes the weights a genuine positive quadrature rather than a signed cancellation device.
The centered response is real and affine in the coefficients. The finite design therefore exposes conditioning, density and arithmetic leakage as explicit constraints rather than as informal window folklore.
w_j(c)=\frac1m\left[1+2\sum_{r=1}^{H}c_r\cos\left(\frac{2\pi r t_j}{T}\right)\right]R_c(\omega)=e^{-iT\omega/2}K_c(\omega)2. The unavoidable midpoint-grid alias theorem
Exact statement or rational verificationLet arbitrary complex weights satisfy sum_j w_j=1. At every integer grid alias omega_q=2 pi q m/T, every midpoint has the same phase e^(i t_j omega_q)=(-1)^q. The common phase factors out of the weighted sum, so K_w(omega_q)=(-1)^q and |K_w(omega_q)|=1.
This result does not require positivity or symmetry. Weight optimization can suppress leakage before an alias, but it cannot impose a strict alias-height constraint on a fixed uniform grid. Alias location must be controlled jointly through m/T, and the remote tail must still be bounded analytically.
K_w(\omega_q)=(-1)^q,\qquad |K_w(\omega_q)|=1\omega_q=\frac{2\pi q m}{T}- Exact aliases are a structural property of the sampling schedule, not an optimization failure.
- A finite optimizer that ignores alias location can certify the wrong problem.
Sources: [1]
3. Finite linear programming and conditioning
Exact statement or rational verificationFor target frequencies Delta_nl=log(n/l), auxiliary variables bound the absolute centered responses. A Gershgorin constraint sum_{l != n} z_nl <= 1-gamma gives lambda_min(G_c) >= gamma. Sampled positivity and an upper density cap C keep the measurement positive without concentrating all weight on a few readings.
The divisor-weighted design tail N<k<=M_d is minimized by a linear program. When feasible, its optimum is the global minimum of that declared finite leakage proxy over the declared order-H cosine family, sampled positivity constraints, density cap and Gershgorin floor. It is not a claim of global optimality over all windows or schedules.
\sum_{\ell\ne n}z_{n\ell}\le 1-\gamman^\sigma\sum_{N<k\le M_d}\tau(k)k^{-\sigma}u_{nk}\le q4. Independent complete-tail certification
Exact statement or rational verificationAfter solving the finite program, the implementation discards its auxiliary variables and recomputes the exact realized kernel, exact target Gram spectrum, divisor-envelope correlations through a separate truncation, an analytic remainder beyond that truncation, the final coefficient bound and the weighted Gaussian sensor-noise covariance.
For the order-eight cosine window, the safe pre-alias endpoint is L_H=pi m/T-H(2 pi/T). The million-term certificate lies inside that endpoint for the reference parameters; the analytic remainder then controls the infinite tail, including the exact remote aliases that weighting cannot remove.
L_H=\frac{\pi m}{T}-H\frac{2\pi}{T}\log(1,000,001)=13.8155<15.65775. Stage III reference result
Numerical experimentAt N=50, sigma=2, T=1,000 and m=5,000, the eight-harmonic tail-optimized point has minimum Gram eigenvalue 0.98182107, complete quadratic-tail bound 0.000232014 and effective sample count 2759.61. It is 117.17 times smaller than the re-audited Hann certificate at the same resources.
The balanced cap-2.0 point has complete bound 0.000495807 and is 54.83 times smaller than Hann. Reporting both points makes the conditioning/leakage/noise Pareto tradeoff visible.
These are discrete reported certificate values, not a continuous response curve. The density cap exposes the leakage/noise tradeoff rather than hiding it.
| design | maximum density | exact lambda_min | effective samples | complete-tail bound | improvement over Hann |
|---|---|---|---|---|---|
| Hann | 2.000 | 0.984390 | 3333.33 | 0.0271845 | 1.0x |
| balanced optimized | 2.000 | 0.982433 | 3253.48 | 0.000495807 | 54.83x |
| tail-optimized | 2.500 | 0.981821 | 2759.61 | 0.000232014 | 117.17x |
6. Noise and held-out arithmetic control
Numerical experimentAt 5,000 readings the optimized design certifies noiseless tail recovery, but the stated sensor noise sigma=0.01 is too large for the Gaussian union bound at that density. At 100,000 readings the computed failure bound is 6.09*10^-10. Tail leakage and sensor noise are separate resources.
A held-out quadratic-field experiment through norm 2,000 had maximum complex coefficient error 0.000002105 and integer-rounding success. This is an empirical falsification control, not part of the universal proof.
The tail bias and sensor-noise term are separate resources; these are reported markers, not a fitted law.
| samples | complete-tail bound | Gaussian failure bound at noise sigma=0.01 |
|---|---|---|
| 5,000 | 0.000232014 | 1.0 |
| 20,000 | 0.0000478558 | 0.06778 |
| 50,000 | 0.0000478554 | 5.13*10^-5 |
| 100,000 | 0.0000478553 | 6.09*10^-10 |
7. Interpretation: arithmetic-aware measurement geometry
Interpretive synthesisThe optimized weights do not reveal a new axiom of numbers. They are a measurement geometry adapted to a known arithmetic spectrum: logarithmic ratios of integers with a divisor-weighted risk. A small number of positive trigonometric corrections can exploit that structure more effectively than a generic window while displaying the conditioning and noise cost.
Arithmetic Sensing IV supersedes and sharpens the remote-alias part of this numerical certificate. It does not erase the exact-alias theorem; it supplies a stronger continuum-positive and alias-band audit.
8. What is proved, computed and open
Exact statement or rational verificationProved here: exact unit-magnitude aliases for every weighting of the midpoint grid; affine centered responses; the finite-family linear program; Gershgorin conditioning from its constraints; and the declared alias-safe complete-tail certificate.
Computed and independently checked: global finite LP solutions, positive realized weights, exact Gram spectra, million-term certificates, density-cap Pareto points and held-out quadratic recovery.
Open: global optimality among all positive schedules, a closed-form extremal window, optimal joint scaling in N, T and m, and sharper all-alias bounds beyond the declared family.
- Floating-point optimization and spectral calculations are computations/certifications, not symbolic formal proofs.
- The active continuation is technical work, not a claim that the central Hann result was invalid.
Reproduction and evidence boundary
From a clone of the public repository, run python optimized_arithmetic_quadrature.py with the reference parameters described in the README, then run python -m unittest discover -v. The full optimization is slower than evaluation of the published coefficient vector; the repository and source files are the reproduction resources rather than pasted code.
Repository CI: unittest passed on Python 3.11 and Python 3.12 at dbec038b.
optimized_arithmetic_quadrature.pyARITHMETIC_SENSING_III.mdREADME.md
Public-safe source manifest
Filenames and hashes identify the reviewed research inputs without exposing local filesystem paths or private caches.
ARITHMETIC_SENSING_III.md— canonical technical manuscript;2536c91dc78ad1563f6f11410167a91725dff529f4448d7973720ba2870c925eoptimized_arithmetic_quadrature.py— public optimizer and certifier;b980e84d8501a6aca794c8265ff679b31136d6e185d3f25a23bcf10b8af69349README.md— public reproduction instructions;b24a76233c15a7fd59c350b7c06c5d7b72441fc6a0929c86bdb318470cf3ede4
Theorem / computation boundary
The alias identity and finite-family LP statement are mathematical results for their declared assumptions. The reported coefficients, spectra, million-term sums and noise bounds are independently checked floating-point computations. They do not establish global window optimality, degree-five recovery, a Riemann result or a physical law.
Sources consulted
- Arithmetic Sensing III: Optimized Quadrature and Unavoidable Aliases — Adelic Arithmetic Research; software repository; retrieved 2026-08-13.
- Arithmetic Sensing optimizer and complete-tail certifier — Adelic Arithmetic Research; software repository; retrieved 2026-08-13.
- Adelic Arithmetic Research reproduction instructions — Adelic Arithmetic Research; software repository; retrieved 2026-08-13.
- Arithmetic Sensing IV: Continuum Positivity, Alias-Band Tails, and Local Information — Adelic Arithmetic Research; software repository; retrieved 2026-08-13.
These classical sources and the public reproduction repository provide context and verification routes; they do not establish project novelty or literature priority.