---
title: Spectral Subband Centroid Frequencies (SSCFs)
url: https://www.emergentmind.com/topics/spectral-subband-centroid-frequencies-sscfs
type: topic
---

# Spectral Subband Centroid Frequencies (SSCFs)

Searching arXiv for recent papers on SSCFs and related ASR usage to ground the article.
arxiv.search({"query":"\"Spectral Subband Centroid Frequencies\" OR SSCF ASR polar coordinate centroid frequency", "max_results": 10, "sort_by": "submittedDate"})
Fetching the most relevant arXiv records for SSCFs and centroid-based speech features.
Spectral Subband Centroid Frequencies (SSCFs) are short-time speech features that estimate the “centre of gravity” of the power spectrum within predefined subbands, thereby encoding where spectral mass is concentrated rather than only how much energy is present. In the cited literature, SSCFs are treated as frequency-centroid descriptors computed over fixed or Mel-spaced filter-bank bands, and are positioned as an alternative or complement to Mel-Frequency Cepstral Coefficients (MFCCs) for automatic speech recognition (ASR) and word recognition. The recent literature emphasizes three recurring roles for SSCFs: as static subband-centroid features, as the basis for trajectory descriptors in a two-dimensional SSCF plane, and as ingredients in ratio- or polar-coordinate parameterizations intended to capture dynamic structure while reducing discontinuities and inter-speaker variability [2212.01245], [2206.07176], [2507.22964].

## 1. Formal definition and signal-theoretic interpretation

SSCFs are defined as subband-wise frequency centroids of a short-time spectrum. In one formulation, the speech spectrum \(P(f)\) is divided into \(M\) overlapping or non-overlapping subbands with lower and upper limits \(l_m\) and \(h_m\), and weighting function \(w_m(f)\). The SSCF in band \(m\) is then

\[
\mathrm{SSCF}_m
=
\frac{
\int_{l_m}^{h_m} f \cdot w_m(f) \cdot P^\gamma(f)\, df
}{
\int_{l_m}^{h_m} w_m(f) \cdot P^\gamma(f)\, df
}.
\]

In the studies considered here, \(\gamma\) is a dynamic-range coefficient and is often set to \(1\) [2212.01245]. An equivalent time-indexed form uses the short-time power spectrum \(P(t,f)\):

\[
\mathrm{SSCF}_{m}(t)
=
\frac{
\displaystyle \int_{l_m}^{h_m} f\,w_m(f)\,\bigl[P(t,f)\bigr]^\gamma\,df
}{
\displaystyle \int_{l_m}^{h_m} w_m(f)\,\bigl[P(t,f)\bigr]^\gamma\,df
},
\quad (\gamma=1 \text{ in the reported experiments}).
\]

A discrete version is also reported, replacing the integrals by sums over FFT bins [2507.22964].

The defining interpretation is consistent across the literature: SSCFs estimate the “centre of gravity” of spectral power in each subband, which makes them akin to formant tracking in spirit, but computed over fixed filter-bank bands rather than by explicit formant estimation [2212.01245]. In the frequency-centroid formulation for non-native English word recognition, the same object is termed a Frequency-Centroid (FC) feature and explicitly identified as equivalent to SSCFs; there, each Mel subband yields a single real-valued centroid in hertz [2206.07176].

This representation differs fundamentally from MFCCs. MFCCs summarize the spectral energy distribution after Mel filtering, logarithmic compression, and a DCT, whereas SSCFs preserve the frequency-weighted location of energy within each subband. The papers therefore frame SSCFs as encoding spectral shape and local spectral shift rather than only energy magnitude [2206.07176].

## 2. Subband construction, extraction, and preprocessing

The cited works use filter-bank-defined subbands, but with two distinct operating regimes. In the non-native English word-recognition study, SSCFs are extracted with a 24-filter triangular Mel filterbank spanning from \(0\) Hz to the Nyquist frequency, with subbands defined by the overlapping triangular windows \(W_k(f)\). The centroid in subband \(k\) is computed as

\[
\mathrm{FC}_k
=
\frac{\displaystyle\int_{0}^{f_{\mathrm{Nyq}}} f\;\lvert S(f)\rvert^2\,W_k(f)\,df}
{\displaystyle\int_{0}^{f_{\mathrm{Nyq}}} \lvert S(f)\rvert^2\,W_k(f)\,df}
=
\frac{1}{E_k}
\int_{f_{k-1}}^{f_{k+1}} f\;\lvert S(f)\rvert^2\,W_k(f)\,df,
\]

with \(K=24\) filters and triangular windows anchored at Mel-warped frequencies \(\{f_0,\dots,f_{K+1}\}\) [2206.07176].

In the ASR studies using SSCF trajectories, the frontend is more compact. One system uses \(M=6\) subbands, yielding SSCF\(_1\) through SSCF\(_6\), and constructs five adjacent-band planes from these values [2212.01245]. A later Vietnamese system likewise uses six triangular Mel-spaced bandpass filters indexed \(m=0,\ldots,5\), identical to the 6-filter MFCC bank employed in that work [2507.22964].

The extraction pipeline reported for SSCF-derived polar features is explicit: 25 ms frames, 10 ms hop, Hamming window, pre-emphasis \(0.97\), and 3-frame moving average smoothing [2212.01245]. The Vietnamese study reports 25 ms windows, 10 ms shift, 0.97 pre-emphasis, 512-point FFT, Hamming window, and liftering \(=22\) for the accompanying MFCCs [2507.22964]. The frequency-centroid study uses a similar but not identical frontend: 20 ms frame length, 10 ms frame shift, pre-emphasis \(0.98\), Hamming window, FFT-based power spectrum computation, and weighted sums over the filterbank [2206.07176].

A high-level extraction procedure common to these reports consists of pre-emphasis, framing, Hamming windowing, FFT, power-spectrum estimation, filter-bank weighting, and then, for each band, a centroid numerator divided by the corresponding band energy denominator [2206.07176]. The additional computation relative to MFCC extraction is described as one extra weighted sum per filter, and the reported cost is negligible relative to FFT and moderate-sized CNN inference [2206.07176].

## 3. Trajectory modeling: angles, polar coordinates, and ratio planes

A central line of work uses SSCFs not merely as framewise static features but as coordinates in low-dimensional acoustic trajectories. Earlier work described in the 2022 ASR paper observed that vowel-to-vowel transitions in the SSCF\(_1\)–SSCF\(_2\) plane are quasi-straight, motivating the definition of a transition angle

\[
\theta = \arctan\!\left(\frac{\Delta \mathrm{SSCF}_2}{\Delta \mathrm{SSCF}_1}\right),
\]

where \(\Delta \mathrm{SSCF}_i\) denotes the difference between the endpoint and start of a short sliding window on axis \(i\) [2212.01245]. The stated difficulty is that \(\arctan(\cdot)\) has a branch cut at \(\pm \pi\), which produces abrupt steps in the measured angle despite smooth underlying trajectories.

To address this, the paper proposes converting adjacent SSCF pairs into polar coordinates. For a point on the \(\mathrm{SSCF}_i\)–\(\mathrm{SSCF}_{i+1}\) plane at time \(t\),

\[
x_t = \mathrm{SSCF}_i(t), \qquad y_t = \mathrm{SSCF}_{i+1}(t),
\]

the polar coordinates are

\[
r_t = \sqrt{x_t^2 + y_t^2}, \qquad
\phi_t = \arctan2(y_t,x_t).
\]

The radius measures distance from the origin in the selected SSCF plane, while the angle measures direction using a continuous two-argument arctangent [2212.01245]. The paper’s motivation is explicitly twofold: to avoid discontinuities and to capture both direction and magnitude of the spectral trajectory.

The 2025 Vietnamese study adopts a related but not identical construction by first forming a ratio plane:

\[
x(t)=\frac{\mathrm{SSCF}_1(t)}{\mathrm{SSCF}_3(t)},\qquad
y(t)=\frac{\mathrm{SSCF}_2(t)}{\mathrm{SSCF}_3(t)}.
\]

These ratio features are then mapped to polar coordinates,

\[
r(t)=\sqrt{x(t)^2+y(t)^2}, \qquad
\theta(t)=\arctan\!\Bigl(\frac{\mathrm{SSCF}_2(t)}{\mathrm{SSCF}_1(t)}\Bigr),
\]

with the interpretation that the radius captures the overall magnitude of spectral shift in the ratio plane and the angle captures the direction of transition [2507.22964]. The same study also notes an alternative formulation in which a two-dimensional point is built from current and previous SSCF values in one band, followed by polar conversion, although the reproducible implementation details specify that the reported system computes polar ratios only on the \(\mathrm{SSCF}_1/\mathrm{SSCF}_3\) versus \(\mathrm{SSCF}_2/\mathrm{SSCF}_3\) plane [2507.22964].

The claim that these trajectory parameters are relatively speaker- and rate-independent appears in both the 2022 and 2025 ASR-oriented papers, though their empirical conclusions about gender independence differ. This suggests that the trajectory abstraction is attractive as a dynamical descriptor, but that its invariance properties depend materially on the exact embedding, normalization, and task [2212.01245], [2507.22964].

## 4. Dynamic parameters, derivatives, and pseudo-\(F_0\)

Dynamic augmentation is a recurring theme in SSCF-based ASR. In the 2022 polar-coordinate study, first- and second-order finite-difference derivatives are defined in the usual way over a frame shift \(\Delta\), for example

\[
\Delta r_t = r_t-r_{t-\Delta}, \qquad
\Delta \phi_t = \phi_t-\phi_{t-\Delta},
\]

with \(\Delta \phi_t\) wrapped into \([-\pi,+\pi]\) if necessary, and similarly

\[
\Delta\Delta r_t = \Delta r_t-\Delta r_{t-\Delta}, \qquad
\Delta\Delta \phi_t = \Delta \phi_t-\Delta \phi_{t-\Delta}.
\]

Using six SSCF bands yields five adjacent-band planes; each plane contributes \(\{r,\phi\}\), producing a 10-dimensional static polar representation and a 30-dimensional representation when \(\Delta\) and \(\Delta\Delta\) are included [2212.01245].

The Vietnamese work instead combines 6 static MFCCs with their \(\Delta\) and \(\Delta\Delta\), producing 18 dimensions, and then appends the ratio-plane \(r(t)\) and \(\theta(t)\), with optional inclusion of SSCF\(_0\) as a pseudo-\(F_0\). This yields either 20 dimensions without SSCF\(_0\) or 21 dimensions with SSCF\(_0\) [2507.22964].

The treatment of SSCF\(_0\) is notable. The 2025 study defines SSCF\(_0\) as the centroid over the lowest-frequency band and uses it as a robust stand-in for the true fundamental frequency. To compensate gender-dependent pitch scaling, utterance-level mean-variance normalization is applied:

\[
\widehat{\mathrm{SSCF}_0}(t)
=
\frac{\mathrm{SSCF}_0(t)-\mu_0}{\sigma_0},
\]

where \(\mu_0\) and \(\sigma_0\) are the utterance-level mean and standard deviation [2507.22964].

At the same time, the 2022 paper lists SSCF\(_0\) as a limitation and open question: SSCF\(_0\), described there as relative \(F_0\), moves very differently from higher SSCF bands, and mixing it into trajectories may be physically meaningless unless careful normalization or exclusion is used [2212.01245]. These statements are not contradictory so much as method-dependent. One paper warns against unnormalized inclusion of SSCF\(_0\) in generic trajectory planes; the other uses SSCF\(_0\) as a separate pseudo-pitch channel with explicit normalization. A plausible implication is that SSCF\(_0\) is useful when isolated and normalized, but problematic when indiscriminately fused into higher-band geometric trajectories.

## 5. Integration into recognition pipelines

SSCFs have been integrated into both hybrid DNN–HMM ASR systems and CNN-based closed-set word recognizers.

In the French ASR study, the acoustic model is a hybrid DNN–HMM implemented in Kaldi. The input uses a \(\pm 4\)-frame context splice, producing a \(9 \times D\) input where \(D\) is the feature dimension. The network contains 3 hidden layers of 512 units with tanh activation, followed by a fully connected softmax over approximately 2000 triphone states. Training proceeds through monophone and triphone GMM stages and then DNN training over 20 epochs with learning rate \(0.01 \rightarrow 0.001\). Decoding uses a WFST decoder with a 3-gram LM built by SRILM and a lexicon generated by Phonetisaurus [2212.01245].

The Vietnamese gender-independent ASR study reports essentially the same DNN–HMM backbone: Kaldi nnet2 recipe, 9-frame context, 3 hidden layers of 512 tanh units, 20 epochs, and learning rate \(0.01 \rightarrow 0.001\) [2507.22964]. In that system, SSCF-derived parameters are concatenated with MFCC-based features rather than replacing them outright.

By contrast, the non-native English word-recognition study uses a two-stage CNN. Each utterance is transformed into a time–filterbank matrix for MFCCs and another for FC/SSCFs; the matrices are zero-padded or interpolated to \(256 \times 256\). A single feature type yields a \(256 \times 256\) single-channel input, while fused MFCC+FC input forms a \(256 \times 256 \times 2\) tensor. The CNN contains a first convolutional layer with 32 filters and ReLU activation, a second convolutional layer with 64 filters and ReLU activation, and a fully connected softmax output with 5 units corresponding to the five target words [2206.07176].

The following table summarizes the principal configurations reported in the cited literature.

| Study | SSCF role | Model |
|---|---|---|
| [2212.01245] | Adjacent-band polar coordinates, with optional \(\Delta,\Delta\Delta\) | Hybrid DNN–HMM (Kaldi) |
| [2507.22964] | Ratio-plane polar parameters plus optional SSCF\(_0\) | DNN–HMM backbone, Kaldi nnet2 |
| [2206.07176] | 24-band SSCF/FC maps fused with MFCC maps | Two-stage CNN |

Across these systems, normalization is applied either per speaker or per utterance, following standard Kaldi-style mean and variance normalization in the ASR papers, while the CNN study reports mean normalization and does not report additional whitening or cepstral mean variance normalization beyond zero-padding [2212.01245], [2507.22964], [2206.07176].

## 6. Empirical performance and comparative findings

The three papers support distinct but related empirical conclusions.

For French read-speech ASR on BRAF100, the 2022 study reports the following DNN–HMM word error rates (WER). In the normal train/test split, 13-MFCC+\(\Delta+\Delta\Delta\) achieves \(7.14\%\), while raw angle features perform extremely poorly at \(94.19\%\). Raw polar coordinates reduce WER to \(16.44\%\), and Polar+\(\Delta+\Delta\Delta\) further reduces it to \(13.01\%\). The paper states that Polar+\(\Delta+\Delta\Delta\) significantly outperforms raw polar and angle features, with \(p<0.05\) via paired \(t\)-test over multiple test sentences [2212.01245].

In cross-gender recognition on the same corpus, 13-MFCC+\(\Delta+\Delta\Delta\) yields \(16.74\%\) WER in the cross-male condition and \(13.45\%\) in the cross-female condition, whereas Polar+\(\Delta+\Delta\Delta\) yields \(21.51\%\) and \(29.20\%\), respectively. The conclusion drawn in that study is that, although polar coordinates are superior to transition angles for defining acoustic trajectories, they are not much more gender-independent than conventional MFCCs [2212.01245].

The non-native English word-recognition study reaches a different conclusion about SSCFs’ role: SSCFs alone are not competitive in quiet conditions, but fusion with MFCCs consistently improves performance. Clean-speech accuracy for MFCC alone is \(0.78\), \(0.80\), and \(0.75\) for Arabic, French, and Spanish accents, respectively; FC alone yields \(0.48\), \(0.49\), and \(0.52\); MFCC+FC yields \(0.80\), \(0.82\), and \(0.80\). Under NOISEX-92 additive noise, MFCC+FC outperforms MFCC alone by 3–8 percentage points at every SNR level and for each of White, Babble, and HF-channel noise. At \(0\) dB White noise, the reported improvement is from approximately \(45\%\) to approximately \(55\%\) [2206.07176].

The Vietnamese study reports that adding only the SSCF-derived polar-ratio parameters and SSCF\(_0\) to a 6-MFCC+\(\Delta+\Delta\Delta\) baseline reduces WER from \(12.77\%\) to \(9.50\%\) when SSCF\(_0\) is mean-variance normalized, and reduces the male/female WER gap from \(3.59\%\) to \(1.55\%\). The paper explicitly states that the proposed parameters exhibit greater gender independence than the baseline MFCCs and match or exceed the 13-MFCC+\(\Delta+\Delta\Delta\) baseline, which reports \(9.46\%\) all-speaker WER and a \(1.98\%\) male/female gap [2507.22964].

These findings can be organized succinctly as follows.

| Setting | Main comparison | Reported outcome |
|---|---|---|
| French BRAF100 ASR [2212.01245] | Angle vs Polar vs Polar+\(\Delta+\Delta\Delta\) | Polar features improve strongly over angles; derivatives further improve WER |
| French cross-gender ASR [2212.01245] | Polar+\(\Delta+\Delta\Delta\) vs 13-MFCC+\(\Delta+\Delta\Delta\) | No significant relative advantage in cross-gender robustness |
| Non-native English word recognition [2206.07176] | MFCC vs FC vs MFCC+FC | SSCFs alone weak; fusion improves clean and noisy accuracy |
| Vietnamese ASR [2507.22964] | 6-MFCC baseline vs +Polar-ratio vs +Polar-ratio+SSCF\(_0\) | Lower WER and smaller male/female gap with SSCF-derived additions |

A common empirical pattern is therefore that SSCFs are most effective when used as complementary or structured descriptors rather than as stand-alone replacements for MFCCs.

## 7. Limitations, interpretive issues, and research directions

The available literature presents SSCFs as technically promising but methodologically sensitive. Several limitations are stated explicitly in the 2022 polar-coordinate study. First, SSCF\(_0\) behaves very differently from higher SSCF bands, so mixing it into the same trajectory representation may be physically meaningless without careful normalization or exclusion. Second, beyond the first SSCF planes, trajectories are not quasi-linear, suggesting that the straight-line intuition motivating angle or polar parameterization may weaken for higher bands. Third, \(\Delta\) and \(\Delta\Delta\) are described as crude approximations of dynamic gestures, implying that more sophisticated temporal representations may be preferable [2212.01245].

The same paper identifies future directions including per-speaker or per-gender normalization of SSCFs, especially for the \(F_0\)-related band; end-to-end embeddings of SSCF trajectories via RNNs or CNNs on two-dimensional “SSCF maps”; and unsupervised spatial normalization of polar trajectories to reduce inter-speaker variability [2212.01245]. The Vietnamese work can be read as partially instantiating one of these directions by applying utterance-level normalization to SSCF\(_0\) and using ratio-based geometry to reduce spectral variation [2507.22964].

The differing conclusions about gender independence deserve careful interpretation. One paper concludes that SSCF-derived polar coordinates are not inherently more speaker- or gender-independent than MFCCs under the BRAF100 cross-gender setup [2212.01245], whereas another reports greater gender independence than baseline MFCCs on Vietnamese speech when using polar-ratio parameters and normalized SSCF\(_0\) [2507.22964]. The concrete experimental conditions differ in corpus, language, feature composition, and normalization. This suggests that gender independence is not an intrinsic property of SSCFs alone, but depends on how SSCFs are parameterized and normalized within the larger frontend.

Another recurring interpretive issue is the relationship between static SSCFs and dynamic SSCF-derived geometry. The frequency-centroid study supports SSCFs primarily as complementary static features, especially under additive noise, and explains their utility by their encoding of the shape and shift of spectral subbands [2206.07176]. The ASR studies instead emphasize geometric trajectories in SSCF planes, with polar coordinates intended to characterize direction and magnitude of movement while avoiding angle discontinuities [2212.01245], [2507.22964]. A plausible implication is that SSCFs support two distinct methodological families: centroid augmentation of conventional spectral frontends, and low-dimensional dynamical modeling of articulatory-acoustic transitions.

Taken together, the literature characterizes SSCFs as a compact spectral representation with multiple operational forms: first-moment descriptors of Mel or fixed subbands, geometric coordinates in SSCF planes, and normalized low-frequency proxies for \(F_0\). Their practical value lies less in displacing MFCCs universally than in exposing spectral-center information, transition direction, and trajectory magnitude that MFCC-only pipelines do not directly encode [2206.07176], [2212.01245], [2507.22964].

Source: https://www.emergentmind.com/topics/spectral-subband-centroid-frequencies-sscfs