Sideband-Aware Stability Monitor for Periodic Training
Implementation & benchmark of arXiv:2607.09275 — A Multi-Frequency Input-Admittance Model of Locomotive Rectifier Considering PWM Sideband Harmonic Coupling in Electrical Railways
Source paper: A Multi-Frequency Input-Admittance Model of Locomotive Rectifier Considering PWM Sideband Harmonic Coupling in Electrical Railways arXiv:2607.09275 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Replace the usual averaged Jacobian test for a periodically modulated neural update with a finite harmonic-transfer model that explicitly couples perturbation frequencies separated by the modulation frequency. Use the resulting lifted spectral radius to cap the learning rate or reduce modulation amplitude when sideband interactions create an instability that is invisible in the averaged model.
Formulas
Mathematical statement
Linearize one period of a periodically modulated training or recurrent update as \(\dot{\delta x}=A(t)\delta x\), where \(\delta x\) is a parameter or state perturbation and \(A(t)=\sum_{m=-M}^{M}A_m e^{im\omega_s t}\) has period \(T_s=2\pi/\omega_s\). A perturbation at frequency \(\omega\) has Fourier components \(\delta x_k\) at sidebands \(\omega+k\omega_s\), and harmonic balance gives \(j(\omega+k\omega_s)\delta x_k=\sum_m A_m\delta x_{k-m}\). Truncating to \(k\in[-K,K]\) produces a lifted matrix whose spectral radius or resolvent gain determines amplification. For a discrete update, stability requires \(\rho(\mathcal{M})<1\), where \(\mathcal{M}=\prod_{r=0}^{S-1}(I+\eta A_r)\) is the one-period monodromy matrix, \(\eta\) is the base step size, and \(A_r\) is the Jacobian at phase \(r\). Retaining off-diagonal sideband blocks can reveal instability near or above half the modulation frequency that is absent from the averaged model.
Implementation notes
1. Integration point: apply this to an optimizer with a periodic learning-rate, momentum, gradient-clipping, or data-order schedule, or to an SSM/RNN whose recurrent matrix is periodically gated. Every \(S\) optimizer steps, estimate the local parameter-space Jacobian at each schedule phase. For a cheap MVP, use a two-layer MLP on MNIST with cosine learning-rate modulation of periods \(S=8,16,32\), and compare ordinary SGD or momentum SGD with the sideband-aware controller. 2. Pseudocode: for each phase \(r=0,\ldots,S-1\), compute Jacobian-vector products \(A_rv\) with autograd at the current minibatch or a fixed probe batch; apply the phase operators successively to Lanczos or power-iteration vectors using \(v\leftarrow v+\eta A_rv\); estimate \(\rho(\mathcal M)\). If the estimate exceeds \(1-\epsilon\), shrink \(\eta\leftarrow0.8\eta\) or reduce modulation amplitude. Optionally estimate Fourier coefficients \(A_m=S^{-1}\sum_r A_re^{-2\pi imr/S}\) and construct the block sideband matrix with \(K=1\) or \(2\). 3. Computed versus empirical: monodromy and lifted operators are computed from Jacobian-vector products; minibatch noise and nonlinear basin effects are measured empirically. The averaged baseline uses only \(\bar A=S^{-1}\sum_rA_r\), deleting off-diagonal sideband blocks. 4. First experiment: train the same MLP on MNIST across learning-rate sweeps and modulation periods \(S=8,16,32\), recording loss divergence, gradient norm, \(\rho(\mathcal M)\), and the averaged-Jacobian prediction. The falsifiable prediction is that the actual divergence threshold follows \(\rho(\mathcal M)=1\) within 20%, while the averaged model increasingly overestimates the safe learning rate when modulation amplitude is large or perturbation content is near the schedule's Nyquist frequency. A second signature is a sharp rise in sideband gain near \(\omega\approx\omega_s/2\) or its aliases.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Not reproduced ✗
Stage 2 · Practical benchmark result: Significantly worse ✗
Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works
Stage 1 — Mechanism check agent confidence 8/10
Built a reproducible periodic two-phase shear MVP with exact monodromy, averaged-Jacobian prediction, empirical perturbation growth, Fourier coefficients, and a finite lifted-model prototype. The mechanism manifested: at g=0 exact and averaged models agree; for g=0.05–0.5 the rho gap at eta=0.3 scales quadratically, with gap/g^2=0.0227–0.0249; at g=1.5 the exact stability boundary is eta≈0.50 versus eta≈1.1429 for the averaged model, and eta=0.6 is exact-unstable but averaged-stable. This supports the sideband/noncommutativity claim in the toy setting, but is not evidence of a neural-network training win.
- Agent confidence
- 8/10
- Baseline
- Averaged-Jacobian threshold: eta_c=1.142857 at g=1.5; rho_avg(eta=0.6)=0.8281 and remains below 1.
- Idea
- Exact monodromy threshold: eta_c≈0.500000 at g=1.5; rho_exact(eta=0.6)=1.0498, with empirical perturbation norm ratio after 30 periods ≈15.77. At eta=0.3, exact-minus-averaged rho gap is 0.06650. Small-g quadratic coefficients gap/g^2 range 0.02274–0.02490.
Limitations: Only a 2D synthetic linear periodic system was tested; no MNIST/MLP, optimizer controller, JVP estimation, CUDA run, FLOP/speed comparison, noise robustness, or period/frequency sweep was performed. The prototype lifted matrix is poorly normalized in this discrete formulation and returned rho≈1 across modulation values, so the validated result is the exact monodromy versus averaged model rather than the lifted spectral-radius implementation.
How to run: python3 sideband_monitor.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2607.09275)
Significantly WORSE than the tuned baseline (+1.0%, p=0.0324, wins 2 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.3012
- Idea mean
- 0.3044
- Effect (Δ)
- 0.0031 (+1.0%; negative = idea better)
- Wins
- 2 / 8 paired seeds
- p-value
- 0.0324 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±0.7%
- Mechanism
- Not confirmed ✗
- Practical effect
- Hurts
- Baseline tuning
- swept over 6 configs
The benchmark used the built-in dynamics track, 12 epochs, 800 training samples, and a small three-setting idea sweep. The monitor used a one-vector finite-difference phase-ordered estimate rather than a full harmonic-transfer/JVP Lanczos implementation; periods 16/32, momentum, and speed/FLOP comparisons were not tested.
How to run: python3 bench_sideband.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Stage 2 — Benchmark latest Failed ✗
Agent confidence: 9/10
Implemented and evaluated the sideband-aware periodic learning-rate controller on the registered dynamics track with the shared rnn_small model. The idea reduced modulation amplitude, but test MSE was worse than tuned periodic SGD: 0.304377 versus 0.301230, with paired delta +0.003148 and permutation p=0.03245. Therefore the idea did not transfer as a training win, and its mechanism signature was not quantitatively confirmed.
- Baseline
- Mean test MSE 0.3012298942 over 8 seeds; tuned config lr=0.004, amp=0.8, period=8.
- Idea
- Mean test MSE 0.3043774869 over 8 seeds; config lr=0.004, amp=0.8, period=8, threshold=0.98; paired delta +0.0031475928, 2/8 wins, p=0.03245.
Benchmark result
Significantly WORSE than the tuned baseline (+1.0%, p=0.0324, wins 2 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.3012
- Idea mean
- 0.3044
- Effect (Δ)
- 0.0031 (+1.0%; negative = idea better)
- Wins
- 2 / 8 paired seeds
- p-value
- 0.0324 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±0.7%
- Mechanism
- Not confirmed ✗
- Practical effect
- Hurts
- Baseline tuning
- swept over 6 configs
The benchmark used the built-in dynamics track, 12 epochs, 800 training samples, and a small three-setting idea sweep. The monitor used a one-vector finite-difference phase-ordered estimate rather than a full harmonic-transfer/JVP Lanczos implementation; periods 16/32, momentum, and speed/FLOP comparisons were not tested.
How to run: python3 bench_sideband.py
Artifacts
- bench_report.json 7.8 KB View Raw
- bench_sideband.py 6.6 KB View Raw
- report.md 1.8 KB View
- report_bench_2026-08-31T000703.md 4.8 KB View
- results.json 4.2 KB View Raw
- sideband_monitor.py 3.8 KB View Raw