A Perturbative-Based Generalized Series Expansion in Terms of Non-Orthogonal Component Functions ()
1. Introduction
We propose the following novel approximate generalized series expansion for a function
in terms of a finite number of component functions
where the component functions do not necessarily satisfy the orthogonality condition or form a basis for the function space:
(1)
Let
be a square-integrable function and let
(2)
be the Fourier transform of
.
In the expansion proposed in Equation (1), the component functions
may be orthogonal or non-orthogonal. In the event that there is an infinite number of component functions that are orthogonal, satisfying the condition that any
and
, on the interval
are such that
, the expansion is exactly yielding an equality in the expre- ssion shown in Equation (1). For a finite set of non-orthogonal component functions, the expansion is approximate.
The coefficients
may be determined using a perturbative approach where it is assumed that, at any given value of
, the contribution of each component function
to the function
is perturbed by contributions from the other component functions in the set. From the perspective of the Fourier domain, it may be said that each frequency component of
consists of contributions from the corresponding frequency component of each component function
where
(3)
is the Fourier transform of the
component function
.
The above description associated with a perturbed contribution of each of the component functions
to the corresponding frequency components of
can be written mathematically as shown in Equation (4). The coefficients
are the variables associated with the relative magnitude of the perturbative contribution of the other component functions to the frequency component of the component function
.
(4)
The expression in Equation (4) may be written more compactly as in Equation (5).
(5)
Evaluating Equation (4) or Equation (5) at a specific frequency
for
and re-arranging so as to isolate the terms with the variables
on the right hand side yields Equation (6).
(6)
The system of equations described by Equation (6) is such that
which constitutes an over-determined linear system. The value of
is the chosen number of frequency points at which (6) is evaluated. It is then possible to solve the linear equations, in the least squares sense, for the complex values of the variables
.
Using the decomposition in Equation (4), we define the following super- position of each component function
to be
:
(7)
The complex estimator
is defined as:
(8)
Since
is complex, its magnitude
will yield the multiplicative coefficient
by which each component function must be scaled as per Equation (1)
(9)
where
is the complex conjugate of
.
2. Examples
2.1. Example 1. Bessel Functions
In this first example, we utilize the proposed series expansion to estimate the coefficients associated with the sum of non-orthogonal Bessel functions
of the first kind, where
is a positive integer, as described in Equation (10) [1] :
(10)
A compound function was generated through the superposition of a series of Bessel functions of the first kind with a set of coefficients as per Equation (11):
(11)
Figure 1 illustrates the results of the proposed series expansion for determin- ing the coefficients of a series of superposed non-orthogonal Bessel functions. A comparison is shown in a) of the template compound function
with the estimated compound function
generated using the estimated coefficients
. The graph in b) is a comparison of the actual distribution coefficient values
, shown in the histogram, with the estimated distribution coefficient values
. A graph of the error magnitude is given in c) between the actual distribution coefficients and the estimated coefficient values
as well as the error magnitude between the actual distribution coefficients and the generalized Fourier series coefficients
. The generalized Fourier series coefficients
are found from Equation (12)
(12)
where
[2] .
2.2. Example 2. The Inverse Problem in Electrophysiology
The inverse problem in electrophysiology involves utilization of recorded biopotentials to determine the characteristics of the biological signal generator that gives rise to the recorded biopotential signals. A specific situation, in the context of the inverse problem, arises when it is necessary to determine the conduction velocity or size distribution of a population of peripheral nerve fibers from a recorded compound evoked potential. The compound evoked potential, to a first approximation, may be viewed as the superposition of the single fiber evoked potentials associated with the individual nerve fibers in the nerve trunk that have been electrically activated. Information related to the size distribution, which is linearly related to the conduction velocity distribution of nerve fibers in the evoked population, is of clinical use because different disease processes selectively impact different segments of the nerve population having different diameter or size [3] [4] .
When the separation of nerve fiber sizes is relatively large, the resultant single fiber evoked potentials are orthogonal to each other and a generalized Fourier series may be used to estimate the size distribution. In the event that the fiber size classes are such that there is significant temporal overlap between their associated single fiber evoked potentials, the generalized Fourier series can yield nonsensical results. In such cases, the perturbative decomposition outlined in this paper may be used to obtain an estimate of the coefficients associated with the series expansion shown in Equation (1) as first demonstrated by Szlavik [5] .
(13)
where
(14)
(a)(b)(c)
Figure 1. Graphs illustrating the results of the proposed series expansion for determining the coefficients of a series of superposed non-orthogonal Bessel functions. A comparison is shown in (a) of the template compound function
with the estimated compound function
generated using the estimated coefficients
. The graph in (b) is a comparison of the actual distribution coefficient values
, shown in the histogram, with the estimated distribution coefficient values
. A graph of the error magnitude is given in (c) between the actual distribution coefficients and the estimated coefficient values
as well as the error magnitude between the actual distribution coeffi- cients and the generalized Fourier series coefficients
. The generalized Fourier series coefficients
are found from Equation (12).
(15)
(16)
(17)
(18)
To illustrate the applicability of the perturbative expansion to the case where the single fiber evoked potentials are non-orthogonal, we generated a population of one hundred nerve fibers, using the distribution shown in Equation (19), with diameters ranging from
to
using a technique outlined by Szlavik et al. [6] . The population was generated with the parameters shown in Table 1 which were established from an empirical study [7] .
(19)
Each fiber
in the population of
fibers was then sorted into one of the fiber size classes
, for
, having a separation in diameter of
where the single fiber evoked potentials, shown in Equations (13) through (18), associated with the
size classes are not orthogonal and are generated using the function proposed by Fleisher et al. in Equation (6) of their paper [8] . All parameters in Equations (13) through (18) are as described in Fleisher and were assigned values
,
,
,
,
,
,
Table 1. Parameter values used in the fiber diameter distribution shown in Equation (19).
,
,
where
is the fiber diameter and
is the propagation distance. The function has been normalized to the current through the second pole such that
where
.
An activation function
was used to determine what stimulus current amplitude was sufficient to excite each fiber with diamater
as per Equation (20) where
and
.
(20)
The time domain representation of the compound evoked potential
may be written as a summation over
fibers in the population, as per Equation (21), where the contribution of a given single fiber evoked potential of a given class
associated with the
fiber is only added in when the stimulus current amplitude
is greater than the activation current for the given fiber diameter as quantified in Equation (20) [9] .
(21)
Figure 2 illustrates the results of the proposed series expansion for estimating the diameter distribution of the nerve fibers contributing to the compound evoked potential for the population generated as described earlier. A comparison is shown in a) of the template compound evoked potential
with the estimated compound evoked potential
generated using the estimated coefficients
. The graph in b) is a comparison of the actual fiber diameter distribution coefficient values
, shown in the histogram, with the estimated fiber diameter distribution coefficient values
. A graph of the error magni- tude is given in c) between the actual distribution coefficients and the estimated coefficient values
as well as the error magnitude between the actual distribution coefficients and the generalized Fourier series coefficients
where the coefficients
are found using Equation (12).
2.3. Example 3. Continous Time Signals
In the theory of continuous time signals, a periodic square waveform with a periodicity of
, a fundamental frequency of
or
in radians per second may be represented using a trigonometric Fourier series expansion as per (22) [10] :
(22)
For a square waveform of the type shown in Figure 3, the average value of the waveform is zero yielding a coefficient
. The waveform may be appro- ximated as a superposition of an infinite series of sinusoidal waveforms and thus the coefficients of the cosine term in the series expansion of (22) are
. Evaluation of the Fourier coefficients for the sinuosidal term in (22) yields
for odd
and
for even
.
(a)(b)(c)
Figure 2. Graphs illustrating the results of the proposed series expansion for determining the diameter distribution of the nerve fibers contributing to a compound evoked potential. A comparison is shown in (a) of the template compound evoked potential
with the estimated compound evoked potential
from the estimated coefficients
. The graph in (b) is a comparison of the actual fiber diameter distribution coefficient values
with the estimated fiber diameter distribution coefficient values
. A graph of the error is given in (c) between the actual distribution coefficients and the estimated coefficient values
as well as the error magnitude between the actual distribu- tion coefficients and the generalized Fourier series coefficients
where the coefficients
are found using Equation (12). The parameters for the neural simulation were
,
,
, where
is the number of diameter classes,
is the stimulus current amplitude,
is the propagation distance and
is the number of equations as per (6).
Figure 3. Square wave with fifty percent duty cycle and zero volt DC offset.
An approximation to (22) may be obtained if we limit the number of terms in the series expansion to a finite value as in (23):
(23)
If we apply the perturbative approach, the sinusoidal functions in (23) may be considered to form the component functions
. The coefficient estimates
may be determined using the technique proposed herein.
Figure 4 illustrates application of the perturbative approach to a square waveform of the type shown in Figure 3 for an approximate series consisting of
terms and a fundamental frequency
. A comparison is shown in a) of the template square waveform estimate
with the approxi- mated waveform
from the estimated Fourier coefficients
. The graph in b) is a comparison of the actual Fourier series coefficient values
with the estimated coefficient values
. A graph of the error is given in c) between the Fourier series coefficients and the estimated coefficient values
. The perturbative expansion provides a robust estimate of the Fourier series coeffi- cients. The accuracy of the estimated coefficients is not unexpected since the component functions are orthogonal.
3. Conclusions
The perturbative based series expansion proposed in this paper provides the mathematical framework that may be used to estimate coefficients of an appro- ximate series expansion of a compound function consisting of the superposition of a set of component functions that may or may not be orthogonal. It has been demonstrated, by way of the examples presented earlier, that, in the case of non-orthogonal component functions, the technique performs significantly better than the generalized Fourier series which can yield nonsensical results such as negative coefficient values. As may be seen from specific examples rela-
(a)(b)(c)
Figure 4. Graphs illustrating the results of the proposed perturbative expansion for determining estimates of the Fourier series coefficients. A comparison is shown in (a) of the template square waveform estimate
with the approximated waveform
from the estimated Fourier coefficients
. The graph in (b) is a comparison of the actual Fourier series coefficient values
with the estimated coefficient values
. A graph of the error is given in (c) between the Fourier series coefficients and the estimated coefficient values
. The parameters for the approximated square waveform were
and
.
ted to the inverse problem in electrophysiology presented in Szlavik [5] , the accuracy of the estimated series coefficients degrades as the degree of temporal overlap, or non-orthogonality, of the component functions increases.
The technique would appear to have broad applicability in electrophysiology and the neurosciences particularly with respect to determination of the charac- teristics of signal generators as related to the inverse problem. Currently investi- gation of the decomposition of postsynaptic potentials into the contributions of constituent receptor-ligand complexes is being undertaken. The technique is general in the sense that it may be applied in situations where a compound function is known to consist of the superposition of orthogonal or non-ortho- gonal component functions which are also known.