The Note-Chord-Voice Framework: Structured Source Separation and Causal Inference for EV Charging Data
Authors: Jiajie Chen, Jinfeng Li
Organizations: Xuteli School, Beijing Institute of Technology, No 8 and 9 Yards, Liangxiang East Road, Fangshan District, Beijing, 102488, China · School of Integrated Circuits and Electronics, Beijing Institute of Technology, No. 5 South Zhongguancun Street, Haidian District, Beijing, 100081, China
Real-world EV charging data exhibit three interlocking pathologies: hardware fragmentation (network timeouts and billing resets split sessions), physical violations (independent energy/duration models produce impossible states like 50 kWh in 10 min on a 7 kW charger), and collider bias (clustering on post-treatment outcomes opens backdoor paths for price elasticity). We propose the Note-Chord-Voice framework, a music-inspired, axiom-driven pipeline that separates data cleaning (Repair Chords), structural pattern discovery (Harmonic Chords), descriptive source separation (NMF Voices), and causal inference into distinct, falsifiable stages. Key innovations: (i) falsification gates (A1-A5, G3, G10) that test data suitability before modeling; (ii) Gamma-initialized NMF with input rescaling for convergence stability from STL decomposition; (iii) tag-based coupon grading (A/B/C/D) to isolate quasi-random treatment from night-time confounders and targeted promotions; (iv) separate per-voice OLS to avoid simplex collinearity; (v) Foote novelty curves for structural regime detection. Applied to the Jiangmen dataset (495,707 sessions, 20 stations, from July 2024 to March 2025), all core axioms pass except G3 (no strong 168 h cycle). NMF achieves R^2=0.9921; the physically constrained duration model yields aggregate R^2=0.5409. Two voices are price-sensitive (beta = -11 to -14 min, p<0.001), of which one is stable (Voice 3, beta=-14.16) and one treatment-driven (Voice 1, beta=-11.10); only the stable voice supports causal claims. Counterfactual simulation shows targeting discounts to price-sensitive voices recovers 52.8% of discount expenditures (~0.85M CNY/year); restricting to the single stable price-sensitive voice yields a more conservative estimate.
Figures & tables
Concept
Symbol
Definition
Role
Note
ni
(ti,Ei,Δti,pi,ci) subject to physical manifold M
Atomic validated event
Repair Chord
Cr
Information-theoretic reconstruction operator merging fragmented notes when ΔMDL>0
Data reconstruction
Harmonic Chord
Ch
Recurrent local tensor motif discovered via Matrix Profile on residuals
Residual pattern discovery
Voice
Vk
Horizontal latent stream operationalized as non-negative basis vector in NMF factorization
Source separation
Measure
Γ(t)
Hierarchical metrical field (e.g. 24 h , 168 h , and holidays) via STL decomposition
Baseline intensity
Movement
Mj
Structural regime change in station load, detected by Foote novelty method
Regime segmentation
Table 1: Topological definitions for the Note–Chord–Voice framework.
Axiom
Statement
Test
Pass Cond.
Result
A1 Physical Manifold
A non-trivial fraction of sessions obey E/Δt≤1.2Pmax
Binomial proportion: H0 : > 95% violate
p<0.05
PASS (0.03% violations, p≈0 )
A2 Topological Contractibility
Short-gap sessions show MDL compression when merged
Wilcoxon on ΔMDL for gap ∈(0,2]
med >0 , p<0.01
PASS (median log-likelihood ratio =+0.43 , p=2.3×10−25 )
A3 Motif Recurrence
Residuals after Γ(t) contain recurring patterns
Permutation test on MP distance vs. phase-randomized surrogate
p<0.05
PASS (best p=0.000 over 9 tests)
G3 Multi-periodicity
Load exhibits both 24 h and 168 h periodicity
FFT: rank of 24 h and 168 h frequency bins ≤3
both
FAIL (only 24 h strong; 168 h rank ≥7 )
A4 Source Separability
NMF outperforms station-mean baseline
Bootstrap error + ARI stability
err < base, ARI >0.7
PASS ( R2=0.9921 , stable K=5 )
A5 Interaction Structure
(a) Station–hour interaction >5% (b) NMF stable to treatment inclusion
Table 2: Axiom falsification gates: statement, test, pass condition, and result on the Jiangmen dataset.
Figure 1: Note–Chord–Voice pipeline overview. The Axiom Gate (Step 0.5) tests data suitability before any complex model is built, ensuring graceful degradation if assumptions fail.
Grade
Description
Identification rationale
Role
Grade A
Generic platform-wide coupons, no targeting, available to all users
Exogenous: orthogonal to session timing and user characteristics
Primary instrument
Grade B
Night-time specific coupons (e.g., “off-peak discount”)
Endogenous: severely confounded by circadian timing
Correlated with user behavior and station characteristics
Excluded
Grade D
Other/unknown promotional tags
Ambiguous; excluded to maintain instrument purity
Excluded
None
No discount applied
Natural control group
Control
Table 3: Coupon grading taxonomy for Jiangmen promotions.
Figure 2: Physical manifold: energy vs. duration on log–log axes, colored by charger type. Diagonal lines show power limits ( E/Δt=Pmax ) for 7 kW, 30 kW, 120 kW, and 480 kW chargers. The Γ -initialized NMF is fitted to data within these boundaries, while unconstrained models might produce impossible states.
Figure 3: A2 axiom test results by gap bracket. Left: Median log-likelihood ratio (negative for back-to-back sessions, positive for short gaps). Right: Fraction of session pairs where merging is statistically beneficial. The fragmentation signal concentrates in the (0,2] minute gap bracket, contradicting the assumption that back-to-back sessions are the primary fragmentation source.
Figure 4: NMF voice selection: reconstruction error (left axis) and stability (ARI, right axis) vs. number of voices K . K=5 is selected as the smallest K passing both A4 and A5(b). K=3 achieves higher ARI but fails A5(b) (only 1/3 voices stable).
Figure 5: NMF voice temporal profiles (rows of H ). Each panel shows one voice’s normalized energy intensity across 168 hours of the week. Solid lines = stable voices (support causal claims); dashed lines = treatment-driven voices (descriptive only).
Figure 6: Station voice composition (normalized NMF responsibilities). Each bar shows how a station’s load is distributed across the 5 voices. Stations dominated by a single voice have clear behavioural signatures; mixed stations reflect overlapping user populations.
Figure 7: Movement detection for a representative station. Top: cosine self-similarity matrix of daily load profiles. Middle: Foote novelty curve (globally normalized) with detected changepoints. Bottom: daily energy load with Major (solid lines) and Minor (dashed lines) Movement boundaries.
Figure 8: Confounder check by coupon grade. Left axis: average price per kWh (blue); right axis: percentage of night-time sessions (orange). Grade B has an 86% night-time session concentration and average price 0.575 CNY/kWh, confirming that night-time electricity pricing—not coupons—drives its low price.
Figure 9: Voice-specific heterogeneous treatment effects (Grade A vs. None). Circles = stable voices (support causal claims); diamonds = treatment-driven voices (descriptive only). Red = price-sensitive ( β<−5 minutes, p<0.1 ); blue = price-insensitive. Voice 3 (stable, β=−14.16 minutes) is the clean causal finding.
Figure 10: Revenue counterfactuals under alternative discount strategies. Comparison of baseline expenditure, targeted allocation to price-sensitive voices (Voices 1 and 3), and the new-targeted expansion scenario, illustrating net expenditure recovery and revenue impacts.
School of Engineering and Architecture, University College Cork, Cork, Ireland · MaREI Centre, Sustainability Institute, University College Cork, Cork, Ireland