Sequential RC-TGAN: Generating Relational Time Series with Spectral Envelope Loss
^††thanks: This work was supported by Mitacs through the Mitacs Accelerate program.

Let $\{S_{t},t\in\mathbb{Z}\}$ be a weakly stationary process with value on $\mathbb{R}$, with mean $\mu\in\mathbb{R}$ and autocovariance function $\gamma(h)=\text{cov}(S_{t+h},S_{t})$. The spectral density $f(\omega)$ describes how the variance of the process is distributed across frequencies $\omega\in[-1/2,1/2]$. The spectral density is defined as the Fourier transform of the autocovariance:

$$f(\omega)=\sum_{h=-\infty}^{\infty}\gamma(h)e^{-2\pi i\omega h}.$$

(1)

Inversely : $\gamma(h)=\int_{-1/2}^{1/2}{f(\omega)e^{2\pi i\omega h}d\omega}.$

In practice, for a finite time series $\{s_{1},\dots,s_{T}\}$, the spectral density is estimated using the periodogram $I(\omega_{k})$, calculated at Fourier frequencies $\omega_{k}=k/T$:

$$I(\omega_{k})=\left|d(\omega_{k})\right|^{2}=\frac{1}{T}\left|\sum_{t=1}^{T}s_{t}e^{-2\pi i\omega_{k}t}\right|^{2},$$

where $d(\omega_{k})=\frac{1}{\sqrt{T}}\sum_{t=1}^{\prime}s_{t}e^{-2\pi i\omega_{k}t}$, the Discrete Fourier Transform (DFT).

The concept of spectral density can be extended to the multivariate case. Let $\{\textbf{S}_{t}\}_{t\in\mathbb{Z}}$ be a weakly stationary process in $\mathbb{R}^{q}$ with mean $\mu$ and autocovariance matrix $\Gamma(h)=\mathbb{E}[(\textbf{S}_{t+h}-\mu)(\textbf{S}_{t}-\mu)^{\prime}]$. The $(j,p)$-th entry of this matrix is the cross-covariance function $\gamma_{jp}(h)=\text{cov}(\textbf{S}_{j,t+h},\textbf{S}_{p,t})$, which measures the covariance between component $j$ at time $t+h$ and component $p$ at time $t$.
The spectral density matrix $\textbf{f}(\omega)\in\mathbb{C}^{q\times q}$ is defined as the Fourier transform of $\Gamma(h)$:

$$\textbf{f}(\omega)=\sum_{h=-\infty}^{\infty}\Gamma(h)e^{-2\pi i\omega h},\quad\omega\in[-1/2,1/2].$$

The diagonal elements $\textbf{f}_{jj}(\omega)$ represent the univariate spectral densities, while the off-diagonal elements $\textbf{f}_{jp}(\omega)$ denote the cross-spectral densities.
For a finite observation $\{\textbf{s}_{1},\dots,\textbf{s}_{T}\}$, the spectral density is estimated using the periodogram. Let $\textbf{d}(\omega_{k})\in\mathbb{C}^{q}$ be the Discrete Fourier Transform (DFT) at frequency $\omega_{k}=k/T$:

$$\textbf{d}(\omega_{k})=\frac{1}{\sqrt{T}}\sum_{t=1}^{T}\textbf{s}_{t}e^{-2\pi i\omega_{k}t}.$$

The multivariate periodogram matrix $\textbf{I}(\omega_{k})$ is defined as the outer product:

$$\textbf{I}(\omega_{k})=\textbf{d}(\omega_{k})\textbf{d}(\omega_{k})^{*},$$

(2)

where ^∗ denotes the conjugate transpose. While $\textbf{I}(\omega_{k})$ is an asymptotically unbiased estimator of $\textbf{f}(\omega_{k})$, it is not consistent; its variance does not vanish as $T\to\infty$. Consequently, consistent estimation requires smoothing techniques, such as windowing or averaging over frequency bands.

Consider a categorical time series $X_{t}$ taking values in a finite set $\textbf{a}=\{a_{0},\dots,a_{K-1}\}$ that is stationary. Because standard frequency-domain tools cannot be directly applied to discrete qualitative symbols, we assign a vector of numerical scaling values $\beta=(\beta_{0},\dots,\beta_{K-1})^{\prime}\in\mathbb{R}^{K}$ to the categories in a. This transformation results in a real-valued numerical process, denoted $X_{t}(\beta)\in\mathbb{R}$, where $X_{t}(\beta)=\beta_{k}$ whenever the original series is in state $X_{t}=a_{k}$. By explicitly mapping the qualitative categories to quantitative scalars, we convert the discrete sequence into a standard univariate continuous-state time series. This mathematical conversion is a strict prerequisite, as it enables the calculation of autocovariance functions and the subsequent computation of the spectral density via the Fourier transform.

Instead of assigning arbitrary numbers to categories, the spectral envelope framework systematically derives optimal numerical values that expose hidden periodicities within a categorical time series. The primary objective is to find a scaling vector $\beta$ that maximizes the spectral density relative to the total variance at each specific frequency $\omega$. Formally, the Spectral Envelope $\lambda(\omega)$ is defined as:

$$\lambda(\omega)=\sup_{\beta\not\propto\mathbf{1}}\left\{\frac{f(\omega;\beta)}{\sigma^{2}(\beta)}\right\},\quad\forall\omega\in[-1/2,1/2],$$

(3)

where $f(\omega;\beta)$ and $\sigma^{2}(\beta)$ represent the spectral density and the variance of the transformed numerical process $X_{t}(\beta)$, respectively [18]. The condition $\beta\not\propto\mathbf{1}$ explicitly excludes trivial scalings where every category is assigned the exact same numerical value. If $\beta$ is proportional to a vector of all ones ($\beta\propto\mathbf{1}$), the transformed sequence $X_{t}(\beta)$ would merely become a flat, constant series. This would result in a variance of zero ($\sigma^{2}(\beta)=0$), thereby rendering the objective ratio undefined.

This optimization problem can be solved by representing the categorical process as a multivariate point process $Y_{t}\in\mathbb{R}^{K}$ (using one-hot vectors). Let $f_{Y}(\omega)$ be the spectral density matrix and $V_{Y}$ be the variance matrix of stationary process $Y_{t}$. The optimization problem in (3) can be re-written :

$$\lambda(\omega)=\sup_{\beta\not\propto\mathbf{1}}\left\{\frac{\beta^{\prime}f_{Y}(\omega)\beta}{\beta^{\prime}V_{Y}\beta}\right\},\quad\forall\omega\in[-1/2,1/2].$$

(4)

This expression is a generalized Rayleigh quotient. The solution $\lambda(\omega)$ is the largest eigenvalue of $f_{Y}(\omega)$ in the metric of $V_{Y}$. The corresponding eigenvector $\beta(\omega)$ is called the optimal scaling at frequency $\omega$.
The value $\lambda(\omega)$ is called the spectral envelope because it envelopes the normalized spectrum of any scaled process $X_{t}(\beta)$. In other words, for any normalized scaling $\beta$ (such that $\sigma^{2}(\beta)=1$), we have $f(\omega;\beta)\leq\lambda(\omega)$, with equality achieved if and only if $\beta$ is proportional to the optimal scaling $\beta(\omega)$.

While the spectral envelope provides a robust mechanism for uncovering the latent periodicities of a single categorical time series, leveraging this concept within a deep generative framework requires systematically comparing the structural properties of real and synthesized processes. To formulate a differentiable objective that minimizes the frequency-domain discrepancy between these temporal dynamics, we cannot merely view $\lambda(\omega)$ as a collection of point-wise maxima. Instead, we must formalize the spectral envelope as a distinct mathematical object residing within a well-defined functional space. This theoretical shift naturally motivates the construction of a metric space for spectral envelopes, providing the foundational distance metrics required to optimize our generative model via backpropagation.

Our formulation is grounded in the Probabilistic Relational Model (PRM) framework [4]. We consider a relational schema $\mathcal{S}=\{W,U\}$ containing two classes (tables) $W$ and $U$, where $W$ acts as the parent entity and $U$ as the child entity.

Let $\mathcal{A}(U)$ denote the set of attributes for table $U$. This set is partitioned into continuous attributes $\mathcal{A}_{cont}(U)=\{c_{1},\dots,c_{I}\}$ and categorical attributes $\mathcal{A}_{cat}(U)=\{d_{1},\dots,d_{J}\}$. The attribute space (or domain) for $U$ is defined as the Cartesian product of the domains of its individual attributes: $\mathcal{V}(U)=\bigotimes_{A\in\mathcal{A}(U)}\mathcal{V}(A)$. Similarly, we define $\mathcal{V}(W)$ as the attribute space for the parent table $W$.

In this relational structure, specific dependencies exist between instances of $W$ and $U$. Let $w\in W$ denote a specific row (instance) in the parent table, with feature values $w.\mathcal{A}\in\mathcal{V}(W)$. We define $\text{Children}(w)\subset U$ as the set of child rows in table $U$ that reference the parent $w$.

In the context of time series synthesis, the set $\text{Children}(w)$ is not merely a bag of rows but an ordered sequence associated with the parent entity. We denote this sequence as $\text{Children}(w).\mathcal{A}=(u_{1},\dots,u_{T})$, where each $u_{t}\in\mathcal{V}(U)$ represents the state of the child entity at time step $t$, and $T$ is the sequence length. Thus, the dataset consists of tuples $\left(w.\mathcal{A},\text{Children}(w).\mathcal{A}\right)$, pairing static parent features with dynamic child sequences.

To address the challenge of generating relational time series, we introduce the Sequential RC-TGAN (see Fig. 1), an extension of the Row Conditional-TGAN (RC-TGAN) [1] model enhanced by the temporal dimension modeling. The original RC-TGAN primarily focused on modeling inter-table relationships, employing a generator $\mathcal{G}$ to model the conditional distribution of a single child row given its parent: $\mathbb{P}(u|w.\mathcal{A})$.

The Sequential RC-TGAN adapts this paradigm to support inter-row relationships modeling inside a tabular data. Rather than mapping parent feature values $w.\mathcal{A}$ and a noise vector $z$ to a static point in the feature space $\mathcal{V}(U)$, our generator learns to map them to a temporal trajectory within $\mathcal{V}(U)^{T}$. Formally, the model approximates the conditional joint distribution of the child sequence given the parent attributes:

$$\mathbb{P}(\text{Children}(w).\mathcal{A}|w.\mathcal{A})=\mathbb{P}(u_{1},\dots,u_{T}|w.\mathcal{A}).$$

(6)

This formulation ensures that the generation process is explicitly conditioned on the static characteristics of the parent entity, thereby guaranteeing that the synthesized temporal dynamics remain consistent with their relational context.

The spectral envelope theory in [18] is inherently designed for categorical time series. However, relational datasets frequently contain continuous numerical attributes $\mathcal{A}_{cont}(U)=\{c_{1},\dots,c_{I}\}$ that exhibit significant periodic behavior (e.g., sales volume, temperature). To incorporate these attributes into our frequency-domain regularization, we first employ a discretization strategy based on Variational Gaussian Mixture Models (VGM) [21].

For each continuous attribute $c_{i}\in\mathcal{A}_{cont}(U)$, we fit a VGM to the training data to estimate the optimal number of modes $K_{c_{i}}$ and their parameters. The probability distribution of a value $u_{t,c_{i}}$ is modeled as a mixture of Gaussians:

$$\mathbb{P}(u_{t,c_{i}})=\sum_{k=1}^{K_{c_{i}}}\tau_{k}\mathcal{N}(u_{t,c_{i}};\mu_{k},\sigma_{k}).$$

(9)

To compute the spectral envelope for a continuous sequence $\mathbf{u}_{c_{i}}=(u_{1,c_{i}},\dots,u_{T,c_{i}})$, we transform it into a discrete sequence of mode indicators $\mathbf{m}_{c_{i}}=(m_{1,c_{i}},\dots,m_{T,c_{i}})$. At each time step $t$, the value $u_{t,c_{i}}$ is assigned to the mode $k$ that maximizes the posterior probability:

$$m_{t,c_{i}}=\arg\max_{k}\left(\tau_{k}\mathcal{N}(u_{t,c_{i}};\mu_{k},\sigma_{k})\right).$$

(10)

This process effectively maps the continuous domain $\mathcal{V}(c_{i})$ to a finite categorical set $\{1,\dots,K_{c_{i}}\}$. Consequently, we can calculate the spectral envelope $\lambda(\omega;c_{i})$ on this discretized sequence, allowing the spectral loss $\mathcal{L}_{spec}$ to enforce periodic consistency across both naturally categorical attributes $\mathcal{A}_{cat}(U)$ and discretized continuous attributes $\mathcal{A}_{cont}(U)$.

Beyond this discrete mode assignment, each continuous value is concurrently represented by a normalized scalar that captures its relative position within the assigned mode. Specifically, if the value $u_{t,c_{i}}$ is assigned to mode $k$, we compute an intra-mode scalar $v_{t,c_{i}}=\frac{u_{t,c_{i}}-\mu_{k}}{4\sigma_{k}}$. By concatenating the one-hot encoded discrete mode indicator $m_{t,c_{i}}$ with this normalized continuous scalar $v_{t,c_{i}}$, the model retains the complete information necessary to fully reconstruct the original continuous feature $u_{t,c_{i}}$. Therefore, while the categorical mode sequence $\mathbf{m}_{c_{i}}$ explicitly drives the frequency-domain regularization via the spectral envelope, the supplementary scalar sequence $v_{t,c_{i}}$ ensures no loss of localized continuous variance in the time domain.

Note that another way to incorporate the continuous attributes into frequency domain is to use the power spectrum of the signal. However as it will be shown in the ablation study section, we find that the discretization method is more effective.

The training of the generator is guided by a hybrid objective function designed to satisfy two complementary requirements: global statistical realism (via adversarial feedback) and frequency-domain fidelity (via spectral envelope matching).

The training procedure employs an alternating optimization strategy to balance the competing objectives. In each epoch, we execute the following three distinct phases:

1. 1.

   Discriminator Update: First, we optimize the discriminator $D$ to distinguish between real sequences and the current synthetic output. We perform $n_{critic}$ updates to the discriminator for every generator update to maintain an optimal gradient approximation for the WGAN objective.

2. 2.

   Adversarial Generator Update: Second, we update the generator $G$ by minimizing $\mathcal{L}_{adv}$. In this step, the generator weights are adjusted to fool the discriminator, ensuring global statistical coherence and adherence to the parent conditioning.

3. 3.

   Spectral Generator Update: Finally, we perform a specialized refinement step focused on frequency-domain fidelity. We update the generator by minimizing $\mathcal{L}_{spec}$. This update is repeated $n_{steps\_for\_spec}$ times per epoch.

Deriving the spectral envelope for a general transition matrix requires numerically solving the eigenvalue problem at every frequency. For the class of circulant transition matrices, we can derive an exact analytical form that links the stochastic parameters directly to the spectral shape.

A transition matrix $P$ is circulant if every row is a cyclic right shift of the preceding row. Consequently, the entire matrix is fully characterized by its first row vector $\textbf{b}=[b_{0},b_{1},\dots,b_{K-1}]$, where $b_{j}=\mathbb{P}(Y_{t+1}=e_{j}\mid Y_{t}=e_{0})$.

The general form of such a matrix is:

$$P=\begin{bmatrix}b_{0}&b_{1}&b_{2}&\dots&b_{K-1}\\ b_{K-1}&b_{0}&b_{1}&\dots&b_{K-2}\\ b_{K-2}&b_{K-1}&b_{0}&\dots&b_{K-3}\\ \vdots&\vdots&\vdots&\ddots&\vdots\\ b_{1}&b_{2}&b_{3}&\dots&b_{0}\end{bmatrix}.$$

(19)

This structural symmetry serves as a mathematical bridge between the time domain and the frequency domain: circulant matrices are diagonalized by the Inverse Discrete Fourier Transform matrix [6], a property we leverage to derive analytical spectral envelopes. Then, the eigenvalues $\gamma_{k}$ of the circulant matrix in (19) is given by:

$$\gamma_{k}=\sum_{j=0}^{K-1}b_{j}e^{i\frac{2\pi jk}{K}}\text{ for }k=0,\ldots,K-1.$$

(20)

The magnitude ($|\gamma_{k}|$) is determined by the concentration of the probability mass in b. If b is highly concentrated (low entropy), the magnitude approaches 1 ($|\gamma_{k}|\approx 1$), implying long memory, whereas a uniform b (high entropy) yields $|\gamma_{k}|\approx 0$, which is characteristic of a white noise process.

Using this lemma, we analyze two types of circulant chains representing distinct temporal dynamics: periodicity and inertia.

This process models periodic behavior with phase noise, serving as a robust benchmark for capturing seasonality and cyclic constraints.

A Noisy Cyclic Process (NCP) is defined by the transition matrix:

$$P_{ij}=\begin{cases}\alpha&\text{if }j\equiv(i+1)\pmod{K}\\ 1-\alpha&\text{if }j=i\end{cases}$$

where $\alpha\in(0.5,1)$ is the switching state parameter. An NCP is a circulant chain where $\textbf{b}=[1-\alpha,\alpha,0,\dots]$. Its spectral envelope is given by equation (21) where :

	$\displaystyle r_{k}$	$\displaystyle=\sqrt{1-2\alpha(1-\alpha)\left[1-\cos\left(\frac{2\pi k}{K}\right)\right]},$
	$\displaystyle\phi_{k}$	$\displaystyle=\arctan\left(\frac{\alpha\sin(2\pi k/K)}{(1-\alpha)+\alpha\cos(2\pi k/K)}\right).$

The NCP serves as a robust benchmark for modeling periodic behavior and cyclic constraints under varying degrees of phase noise. Its temporal dynamics are primarily controlled by the switching state parameter $\alpha\in(0.5,1)$, which dictates the strictness of the cycle. Figure 3 visualizes the spectral envelope of the NCP with a state space of $K=7$ for different values of $\alpha$. As shown, the process naturally exhibits distinct resonant peaks clustered around the fundamental frequency of $1/7\approx 0.14$ and its associated harmonics. When the cycle strength $\alpha$ approaches 1 (represented by the darker lines), the system mimics a deterministic cycle, concentrating the spectral energy into very sharp, Dirac-like peaks. Conversely, as $\alpha$ decreases toward its lower bound, the process introduces greater phase noise, which progressively broadens these sharp harmonic peaks into wide spectral hills, reflecting a more stochastic and relaxed periodic progression.

This process models systems with inertia, where the state tends to persist over time with no preferred direction of change.

A Symmetric Sticky Process (SSP) is defined by the transition matrix:

$$P_{ij}=\begin{cases}\alpha&\text{if }i=j\\ \frac{1-\alpha}{K-1}&\text{if }i\neq j\end{cases}$$

where $\alpha\in(1/K,1)$ is the switching state parameter. We can remark that the transition matrix of the SSP is circulant such that the first row $\textbf{b}=[\alpha,\frac{1-\alpha}{K-1},\dots]$. We can derive the spectral envelope of the SSP from equation (21):

$$\lambda(\omega)=\frac{1-\gamma^{2}}{1-2\gamma\cos(2\pi\omega)+\gamma^{2}},$$

where $\gamma=\frac{\alpha K-1}{K-1}$. The non-trivial eigenvalues are identical and real corresponding to $\gamma$.

The asymptotic behavior of the process is highly sensitive to the switching state parameter $\alpha$. As $\alpha\to 1$, the eigenvalue $\gamma\to 1$, causing the process to become extremely "sticky" and rarely switch states. In this regime, the spectral envelope forms a sharp peak at $\omega=0$, ultimately approaching a Dirac delta. Conversely, as $\alpha\to 1/K$, the eigenvalue $\gamma\to 0$, which reduces the process to pure random noise (see Figure 3). Consequently, the spectral energy becomes uniformly distributed, and the envelope flattens to a constant line where $\lambda(\omega)\approx 1$.

To rigorously evaluate the conditional generation capabilities of our model, we elaborate a method for building a synthetic relational database using a Bayesian hierarchical framework. Unlike using a single fixed parameter for the entire dataset, we model the switching parameter $\alpha$ as a random variable associated with each parent entity. This setup forces the generative model to learn the mapping $\alpha\mapsto\text{Spectral Envelope}(\alpha)$ rather than memorizing a static distribution.

Following the relational schema $\mathcal{S}$ introduced in Section IV, we construct a framework where the parent table $W$ governs the stochastic dynamics of the child time series in table $U$. Note that for any single simulated database, we select only one of these two types of circulant chains (either the Symmetric Sticky or the NCP) and fix the total number of categorical states $K$ to drive the temporal dynamics. Because both processes are fully parameterized by the switching parameter $\alpha$ for a given $K$, we apply a Bayesian hierarchical modeling to each choice: first, we sample and store $N$ values of $\alpha$ in the parent table; second, we generate a corresponding categorical time series for each of these $N$ values using the chosen Markov process.

For the NCP, where $\alpha$ represents the probability of advancing the cycle, we use a uniform prior over the valid range of directed cycles:

$$\alpha_{i}\sim\mathcal{U}(0.5,1.0),$$

where $\alpha_{i}$ is the random variable illustrating the value of $\alpha$ for the $i$-th row in the parent table. For the Symmetric Sticky, we simply replace $\mathcal{U}(0.5,1.0)$ distribution by $\mathcal{U}(1/K,1.0)$.

For each parent row $i$, we generate a categorical time series $\textbf{u}^{(i)}=\{u^{(i)}_{1},\dots,u^{(i)}_{T}\}$ of length $T$, which acts as part of the child table $X$. The dynamics of this series are conditioned strictly on the parent’s parameter $\alpha_{i}$.

The sequence is generated via the transition matrix $P_{\alpha_{i}}$ specific to the chosen type of circulant chain (Symmetric Sticky or Noisy Cyclic):

$$\mathbb{P}(u^{(i)}_{t+1}|u^{(i)}_{t},\alpha_{i})=[P_{\alpha_{i}}]_{u^{(i)}_{t},u^{(i)}_{t+1}}.$$

(22)

This hierarchical construction provides a "gold standard" dataset for conditional generative modeling. Since the true spectral envelope $\lambda(\omega;\alpha_{i})$ is analytically known for every parent $i$ (via Lemma 3), we can compute the exact expected spectral error. A successful generative model must produce synthetic children $\hat{\textbf{u}}^{(i)}$ such that their empirical spectral envelopes match the theoretical envelopes dictated by their sampled parent attributes $\hat{\alpha}_{i}$.

To evaluate the temporal fidelity of continuous features, we measure the divergence between the average power spectrums. We first normalize the spectral densities so that $\int_{-1/2}^{1/2}f(\omega)d\omega=1$, considering them as probability distributions over the frequency domain.

The individual divergence $div(\cdot,\cdot)$ for an attribute $c_{i}$ is defined as the divergence (e.g., Wasserstein or Kullback-Leibler) between the real and synthetic mean densities. The global Spectral Density Divergence (SDD) is the average over all continuous attributes:

$$\overline{\mathcal{D}}_{spec}=\frac{1}{I}\sum_{i=1}^{I}{\mathcal{D}_{spec}(c_{i})},$$

(25)

where $\mathcal{D}_{spec}(c_{i})=div\left(\bar{f}^{(real)}(\cdot;c_{i}),\bar{f}^{(synth)}(\cdot;c_{i})\right)$. It is worth noting that the aggregated quantities $\bar{f}^{(real)}(\cdot;c_{i})$ and $\bar{f}^{(synth)}(\cdot;c_{i})$ constitute well-defined normalized spectral densities. Mathematically, the set of valid normalized spectral densities is closed under convex combinations; since each individual $f(\omega;c_{i},w_{m})$ is a non-negative, real-valued, its integral under interval $[-1/2,1/2]$ is equal to one, and even function (for real-valued processes), their arithmetic mean preserves these fundamental properties. Consequently, $\bar{f}^{(real)}$ effectively represents the spectral density of a "representative" process for the class $U$, averaging out local idiosyncrasies to reveal the global frequency structure of the population. This validity allows us to treat these patterns as probability distributions and rigorously apply divergence metrics $div$ such as Kullback-Leibler (KL) or Wasserstein distance.

For categorical attributes, we assess the preservation of latent periodicities using the spectral envelope. The divergence for a single attribute $d_{j}$ is calculated as the $L^{2}$ distance between the mean envelopes:

$$\mathcal{D}_{env}(d_{j})=\frac{\big\|\bar{\lambda}^{(real)}(\cdot;d_{j})-\bar{\lambda}^{(synth)}(\cdot;d_{j})\big\|_{2}}{K-1}.$$

(26)

The global Spectral Envelope Divergence (SED) is the mean over all categorical attributes:

$$\overline{\mathcal{D}}_{env}=\frac{1}{J}\sum_{j=1}^{J}\mathcal{D}_{env}(d_{j}).$$

(27)

The aggregated spectral envelope $\bar{\lambda}^{(real)}(\omega;d_{j})$ inherits the fundamental algebraic properties of its constituents. Since $\bar{\lambda}^{(real)}$ is constructed as a finite linear combination of individual envelopes, and given that each component $\lambda(\cdot;d_{j},w_{m})$ is continuous on the compact interval $[-1/2,1/2]$ (Lemma 1), the mean envelope itself is necessarily a continuous function in $C^{0}([-1/2,1/2])$. Furthermore, the fundamental norm constraints established in Lemma 2 are preserved under this averaging operation. Specifically, by the convexity of the norm, the mean envelope satisfies the dimensionality constraint $1\leq\|\bar{\lambda}^{(real)}(\cdot;d_{j})\|_{1}\leq K-1$ (If $K$ is the number of categories of $d_{j}$) and maintains finite energy $1\leq\|\bar{\lambda}^{(real)}(\cdot;d_{j})\|_{2}<\infty$. These properties ensure that $\bar{\lambda}^{(real)}$ remains a well-defined term within the functional $L^{2}\left([-1/2,1/2]\right)$ and verify important properties of spectral envelopes, guaranteeing that the divergence metric $\overline{\mathcal{D}}_{env}$ is both bounded and mathematically stable.

The results are summarized in Tables II and I, alongside visual analyses in Figures 5 and 4. To assess model performance, we highlight relative improvements over the second-best baselines and evaluate statistical significance using a two-sample t-test ($p<0.05$).