---
title: 'Linear-Signal Analysis: Foundations and Applications'
url: https://www.emergentmind.com/topics/linear-signal-analysis
type: topic
---

# Linear-Signal Analysis: Foundations and Applications

Linear-signal analysis denotes the study of signals and systems through linear operators, linear time-invariant models, linear descriptor systems, and linearized or multilinear local approximations of more complex dynamics. Its central objects are impulse responses, transfer functions, convolution laws, Fourier- and Laplace-domain representations, and covariance- or Jacobian-based surrogates that preserve superposition or local differential structure. In this broad sense, the term covers continuous and discrete SISO LTI systems, algebraic signal models \((\mathcal{A},M,\varphi)\), linear response to arbitrary periodic excitation, Laplace-based analog-circuit analysis, high-dimensional linear inverse problems, and small-signal stability models for networked power systems [2208.11891], [0612077], [1208.3995], [2001.06702], [2510.16534].

## 1. Foundational formalism

The classical core of linear-signal analysis is the LTI input–output map. For continuous time,
\[
y(t)=(h*u)(t)=\int_{-\infty}^{\infty} h(\tau)\,u(t-\tau)\,d\tau,
\]
and for discrete time,
\[
y[k]=(h*u)[k]=\sum_{l=-\infty}^{\infty} h[l]\,u[k-l].
\]
Causality imposes \(h(t)=0\) for \(t<0\) and \(h[k]=0\) for \(k<0\), while BIBO stability requires absolute integrability or summability of the impulse response [2208.11891].

The same framework yields the transform-domain description. Complex exponentials are eigenfunctions of LTI systems, so the Laplace- and \(z\)-domain relations reduce to
\[
Y(s)=H(s)U(s),\qquad Y(z)=H(z)U(z),
\]
with \(H(s)\) and \(H(z)\) the transfer functions. This is the operational basis for frequency response, pole–zero analysis, and the equivalence between convolution in time and multiplication in transform space [2208.11891].

An algebraic formulation makes this structure explicit by defining a linear signal model as a triple \((\mathcal{A},M,\varphi)\), where the filter space is cast as an algebra \(\mathcal{A}\), the signal space as a module \(M\), and \(\varphi\) generalizes the \(z\)-transform to bijective linear mappings from a vector space of signal samples into \(M\). In that view, the shift operator is the generator of the algebra of filters, and once the shift is chosen, notions of filtering, spectrum, and Fourier transform follow systematically [0612077]. This suggests that linear-signal analysis is not a single technique but a unifying operator-theoretic framework.

## 2. Periodic excitation, transfer functions, and spectral shaping

A central extension of the LTI viewpoint is linear response to arbitrary periodic forcing. In Kubo’s LRT, a weak periodic input \(H(t)=H_0\eta(t)\) drives an observable \(M(t)\) through the causal convolution
\[
M(t)=\int_{-\infty}^{t}\chi(t-t')\,H(t')\,dt',
\]
with response kernel \(\chi(t)\) and susceptibility \(\chi(\omega)\). Expanding \(\eta(t)\) in a Fourier series yields
\[
M(t)=H_0\sum_{n=-\infty}^{\infty} c_n\,\chi(n\omega)\,e^{in\omega t},
\]
so each harmonic is weighted only by its Fourier coefficient and the susceptibility at that harmonic [1208.3995].

For dissipative systems with a single relaxation time \(C(t)=e^{-t/\tau}\), the paper introduces the dimensionless efficiency
\[
\Omega(x)=\frac{A}{\pi H_0^2}
\]
with \(x=\omega\tau\), and derives
\[
\Omega(x)=\chi_0\sum_{n=1}^{\infty}(a_n^2+b_n^2)\frac{n^2x}{1+(nx)^2}.
\]
This separates waveform information from system dynamics and makes waveform design a spectral-allocation problem. Closed forms are obtained for several signals, including the harmonic wave \(\Omega(x)=x/(1+x^2)\), the square wave \(\Omega(x)=\frac{2}{\pi}\tanh\!\left(\frac{\pi}{2x}\right)\), the sawtooth wave \(\Omega(x)=\frac{3}{\pi}\left[\coth\!\left(\frac{\pi}{x}\right)-\frac{x}{\pi}\right]\), and the \(\delta\)-pulse limit \(\Omega_\infty(x)=1/x\) [1208.3995]. This suggests that, under fixed input-energy normalization, changing only the signal shape can change dissipation substantially.

Linear-signal analysis also clarifies ostensibly nonlinear demodulation schemes. The Carré 4-step phase-shifting algorithm can be decomposed into a tunable linear 4-step PSA and a nonlinear phase-step estimator. The linear part has impulse response
\[
\begin{aligned}
h(t)=&\left[-\delta\!\left(t+\frac{3}{2}\right)+\delta\!\left(t+\frac{1}{2}\right)+\delta\!\left(t-\frac{1}{2}\right)-\delta\!\left(t-\frac{3}{2}\right)\right]\cos\!\left(\frac{\omega_0}{2}\right)\\
&+i\left[\delta\!\left(t+\frac{3}{2}\right)+\delta\!\left(t+\frac{1}{2}\right)-\delta\!\left(t-\frac{1}{2}\right)-\delta\!\left(t-\frac{3}{2}\right)\right]\sin\!\left(\frac{\omega_0}{2}\right),
\end{aligned}
\]
and magnitude-relevant transfer function
\[
H(\omega)=\left[1-e^{-i\omega}\right]\left[1-e^{-i(\omega+\omega_0)}\right]\left[1-e^{-i(\omega+\pi)}\right].
\]
Its zeros are at \(\omega=0\), \(\omega=-\omega_0\), and \(\omega=\pi\), and the SNR gain
\[
G(\omega_0)=\frac{|H(\omega_0)|^2}{\frac{1}{2\pi}\int_{-\pi}^{\pi}|H(\omega)|^2\,d\omega}
\]
attains its maximum value \(G_{\max}=4\) at \(\omega_0=\pi/2\) [1203.1947].

A related frequency-domain perspective appears in double-integrator design. The linear double integrator
\[
\dot{x}_1=x_2,\qquad \dot{x}_2=x_3,\qquad \varepsilon^4\dot{x}_3=-k_1\varepsilon x_1-k_2\varepsilon^2x_2-k_3(x_3-a(t))
\]
has transfer functions
\[
G_j(s)=\frac{X_j(s)}{A(s)}=\frac{k_3\,s^{j-3}}{\varepsilon^4s^3+k_3s^2+k_2\varepsilon^2s+k_1\varepsilon},\qquad j\in\{1,2,3\},
\]
and tends to \(1/s^{3-j}\) as \(\varepsilon\to 0\), so \(x_3\) tracks the signal, \(x_2\) approximates \(1/s\), and \(x_1\) approximates \(1/s^2\). The nonlinear counterpart is analyzed by a frequency-sweep method, showing that both integrators estimate low-frequency content accurately and suppress high-frequency noise, while the nonlinear integrator exhibits better estimation performance and stronger robustness [1506.02102].

## 3. Formal and computational realizations

In analog circuit analysis, linear-signal methods typically begin with Kirchhoff laws, proceed to an ODE, apply the Laplace transform, and solve for a rational transfer function
\[
H(s)=\frac{Y(s)}{X(s)}=\frac{b_ms^m+b_{m-1}s^{m-1}+\dots+b_0}{a_ns^n+a_{n-1}s^{n-1}+\dots+a_0}.
\]
The FASiM framework automates this pipeline for Simulink models of linear analog circuits by extracting transfer-function coefficients from MATLAB, translating them into HOL Light lists, constructing the corresponding ODE with `diff_eq_n_order`, and proving ODE \(\leftrightarrow\) transfer-function equivalence in higher-order logic [2001.06702].

Its workflow is explicitly
\[
\text{Simulink block diagram} \rightarrow \text{transfer function }H(s) \rightarrow \text{coefficient lists} \rightarrow \text{rational function and ODE in HOL Light},
\]
and the automation uses the tactics `DIFF_EQ_2_TRANS_FUN_TAC` and `TRANS_FUN_2_DIFF_EQ_TAC` [2001.06702]. For the Sallen–Key low-pass example, the formally verified transfer function is
\[
\frac{V_o(s)}{V_i(s)}=\frac{7812500000.48828}{s^2+62500.0000039063\,s+3906249999.75586},
\]
under differentiability, zero-initial-condition, and Laplace-existence assumptions [2001.06702].

This formal branch of linear-signal analysis makes explicit assumptions that are often left implicit in engineering calculations: linearity, time invariance, exponential order, zero initial conditions, and nonzero denominators. A plausible implication is that, in high-assurance settings, the traditional distinction between “analysis” and “verification” becomes part of the subject itself rather than a post-processing step.

## 4. High-dimensional inference, quantization, and statistical linear models

Modern linear-signal analysis extends beyond deterministic transforms to high-dimensional inference in random linear systems. In generalized linear mixing,
\[
\mathbf{w}=\mathbf{\Phi}\mathbf{x},\qquad
f_{\mathbf{Y|W}}(\mathbf{y}\mid\mathbf{w})=\prod_{i=1}^M f_{Y|W}(y_i\mid w_i),
\]
relaxed BP yields scalar sufficient statistics
\[
q_j=x_j+v_j,\qquad v_j\sim\mathcal{N}(0,\mu),
\]
and reconstruction under any additive error metric
\[
D(\widehat{\mathbf{x}},\mathbf{x})=\sum_{j=1}^N d(\widehat{x}_j,x_j)
\]
reduces to coordinatewise Bayes decisions
\[
\widehat{x}_{\mathrm{opt},j}=\arg\min_{\widehat{x}_j}\mathbb{E}[d(\widehat{x}_j,X_j)\mid q_j].
\]
This decoupling supports MMSE, MMAE, support recovery error, and large-\(p\) \(\ell_p\)-type losses within the same linear mixing framework [1302.0321].

Quantized linear channels preserve the same structure but add a nonlinear observation layer. For the linear vector channel
\[
\mathbf{y}=H\mathbf{x}_0+\sigma_0\boldsymbol{\nu},
\]
followed by scalar quantization, replica analysis in the large-system limit yields a bit-error formula
\[
P_b=Q(\sqrt{E}),
\]
where \(E\) is determined by coupled saddle-point equations. In the noiseless case \((\sigma_0=0)\), the optimal detector admits a perfect-detection solution \(m=1\), \(E=\infty\), hence \(P_b=0\), for any \(d>0\) and \(\beta\). In the noisy case, detection ability decreases monotonically as the quantization step size increases [0805.3406].

A complementary line studies recovery from \(b\)-bit quantized Gaussian measurements while deliberately ignoring the quantizer nonlinearity at reconstruction time. The canonical linear estimator uses
\[
\eta=\frac{1}{m}A^\top y,\qquad 
\widehat{x}\in \arg\min_{x\in\mathcal{C}} -\eta^\top x,
\]
and satisfies
\[
\|\widehat{x}-x_u^*\|_2\le 24\,\frac{\Psi_{b,\sigma}(\mathbf{t},\boldsymbol{\mu})}{\lambda_{b,\sigma}(\mathbf{t},\boldsymbol{\mu})}\,\frac{w(\overline{\Delta}(\mathcal{C}))}{\sqrt{m}}.
\]
Under a fixed bit budget \(B=mb\), the paper concludes that \(b=1\) is optimal for estimating the unit vector corresponding to the signal, while \(b=2\) is optimal for estimating the direction and the norm; Lloyd–Max quantization is optimal with respect to \(\ell_2\)-estimation error [1607.02649].

In 1-bit MU-MISO downlink, Bussgang decomposition restores an effective linear model:
\[
\bar{\mathbf{t}}=\mathcal{Q}_t(\bar{\mathbf{x}})=\mathbf{A}\bar{\mathbf{x}}+\mathbf{q},
\]
so the received signal becomes a linear effective channel plus colored quantization distortion. The achievable sum-rate lower bound is then expressed through an SQINR matrix, and the analysis indicates the merit of reconsidering traditional signal processing techniques, specifically transmit covariance matrices whose rank exceeds the channel rank. The resulting linear precoder design supports higher-order modulation schemes under 1-bit DACs [1802.10329].

Linear diagnostics remain meaningful even when the underlying parameter vector is not. In the high-dimensional model
\[
y=X\beta+\varepsilon,
\]
with \(d/n\to\rho\in[0,\infty)\), the paper proposes consistent, asymptotically normal estimators of residual variance \(\sigma^2\), \(\ell_2\)-signal strength \(\tau^2=\|\beta\|_2^2\), and signal-to-noise ratio, without sparsity assumptions and even when \(d>n\) and the signal itself is non-estimable [1209.0012]. This broadens linear-signal analysis from system characterization to high-dimensional statistical identifiability.

## 5. Adaptive time–frequency analysis and direct component separation

For non-stationary multicomponent signals,
\[
x(t)=\sum_{k=1}^{K}A_k(t)e^{i2\pi\phi_k(t)},
\]
adaptive time–frequency methods retain linear analysis operators while making the window or wavelet parameters time-varying. The adaptive STFT uses
\[
V_x(t,\eta)=\int_{-\infty}^{\infty}x(t+\tau)\,\frac{1}{\sigma(t)}\,g\!\left(\frac{\tau}{\sigma(t)}\right)e^{-i2\pi\eta\tau}\,d\tau,
\]
which remains linear in \(x\) for fixed \(\sigma(t)\) [2010.10760].

Under a sinusoidal local approximation,
\[
x_k(t+\tau)\approx x_k(t)e^{i2\pi\phi_k'(t)\tau},
\]
the adaptive STFT expands as
\[
V_x(t,\eta)=\sum_{k=0}^{K}x_k(t)\,\widehat g\big(\sigma(t)(\eta-\phi_k'(t))\big)+R_0(t,\eta).
\]
Under a linear chirp local approximation,
\[
x_k(t+\tau)\approx x_k(t)e^{i2\pi(\phi_k'(t)\tau+\frac12\phi_k''(t)\tau^2)},
\]
the analogous expansion is
\[
V_x(t,\eta)=\sum_{k=0}^{K}x_k(t)\,G_k\big(\sigma(t)(\eta-\phi_k'(t))\big)+R_0(t,\eta),
\]
with \(G_k\) the Fourier transform of a chirped window [2010.10760].

This leads to direct signal separation operation: extract a ridge \(\eta_k(t)\) by maximizing \(|V_x(t,\eta)|\) in the component’s TF zone, then reconstruct by evaluating the transform on that ridge. For the sinusoidal model,
\[
x_k(t)\approx V_x(t,\eta_k(t)),
\]
whereas the linear-chirp correction gives
\[
x_k(t)\approx \frac{1}{G_k(0)}V_x(t,\eta_k(t)).
\]
For Gaussian windows, the practical real-signal formula becomes
\[
x_k(t)\approx 2\,\Re\left(\sqrt{1-i\,2\pi\sigma(t)^2\phi_k''(t)}\,V_x(t,\eta_k(t))\right),
\]
and explicit error bounds are obtained for both IF estimation and component recovery [2010.10760].

An adaptive continuous wavelet-like transform develops the same idea in a time–scale setting:
\[
W_x(a,b)=\int_{-\infty}^{\infty}x(b+at)\,\frac{1}{\sigma(b)}\,g\!\left(\frac{t}{\sigma(b)}\right)e^{-i2\pi\mu t}\,dt.
\]
Its sinusoidal and linear-chirp local models lead to analogous ridge-based IF estimators and direct reconstruction formulas, including the chirp-corrected form
\[
x_k(b)\approx \frac{1}{G_k(0,a_k(b),b)}\,W_x(a_k(b),b).
\]
The method is presented as an alternative to SST: SST estimates IF and then reconstructs by a definite integral along the estimated IF curve, whereas SSO and the adaptive CWLT avoid the second step and recover the component directly from the transform value on the ridge [2010.13448].

## 6. Applications, scope, and limits

Linear-signal analysis is applied across markedly different physical domains. In magnetic hyperthermia, the LRT framework predicts waveform-dependent heating efficiencies and shows that discontinuous waveforms and narrow pulses can outperform harmonic driving under equal input-energy normalization; the same formalism is then connected to magnetic Langevin or Fokker–Planck modeling of nanoparticles used as local heating centers in oncological treatments [1208.3995].

In analog electronics, the common collector amplifier is analyzed with an Early equivalent transistor model that preserves enough nonlinearity to estimate average voltage gain, total harmonic distortion, and average input and output resistances. The results confirm the importance of the trade-off of current gain and show that sub-optimal performance can occur when base and emitter resistances are not properly chosen; they also report that THD tends to be larger for PNP devices than NPN counterparts with the same average current gain [1805.02705].

In power systems, small-signal stability is formulated directly from implicit multilinear models. After suitable transformations of variables, including representations of trigonometric terms, the multilinear DAE is linearized into a descriptor model, and equilibrium stability is determined from generalized eigenvalues of the pair \((A,E)\). On a 3-bus network example, time-domain simulations of the implicit multilinear model agree with those of the nonlinear model, and the generalized eigenvalues agree with those of the linearized nonlinear model; the decomposed tensor representation also allows faster linearization compared to conventional methods in MATLAB Simulink [2510.16534].

The literature also places clear validity conditions on these methods. LRT requires linearity, small driving amplitude, and near-equilibrium conditions; ideal discontinuous waveforms require infinite harmonic content, and their nonzero \(\Omega(0)\) is explicitly identified as an artifact of ideal discontinuity [1208.3995]. Formal Laplace-domain verification assumes linearity, time invariance, ideal components, and Laplace existence conditions [2001.06702]. Linear recovery from quantized measurements is competitive in moderate- and high-noise settings but can fall short of more sophisticated competitors in a low-noise setting [1607.02649]. This suggests that linear-signal analysis is most powerful when its approximation regime is made explicit: it is neither restricted to exactly linear physics nor universally exact once nonlinearity, finite bandwidth, quantization, or model mismatch dominate.

Source: https://www.emergentmind.com/topics/linear-signal-analysis