Fourier–Bessel series

In mathematics, Fourier–Bessel series is a particular kind of generalized Fourier series (an infinite series expansion on a finite interval) based on Bessel functions.

Fourier–Bessel series are used in the solution to partial differential equations, particularly in cylindrical coordinate systems.

Definition

The Fourier–Bessel series of a function f(x) with a domain of [0, b] satisfying f(b) = 0 $$f:[0,b]\to \mathbb {R}$$ is the representation of that function as a linear combination of many orthogonal versions of the same Bessel function of the first kind J_α, where the argument to each version n is differently scaled, according to $$(J_{\alpha })_{n}(x):=J_{\alpha }\left({\frac {u_{\alpha ,n}}{b}}x\right)$$ where u_α,n is a root, numbered n associated with the Bessel function J_α and c_n are the assigned coefficients: $$f(x)\sim \sum _{n=1}^{\infty }c_{n}J_{\alpha }\left({\frac {u_{\alpha ,n}}{b}}x\right).$$

Interpretation

The Fourier–Bessel series may be thought of as a Fourier expansion in the ρ coordinate of cylindrical coordinates. Just as the Fourier series is defined for a finite interval and has a counterpart, the continuous Fourier transform over an infinite interval, so the Fourier–Bessel series has a counterpart over an infinite interval, namely the Hankel transform.

Calculating the coefficients

As said, differently scaled Bessel Functions are orthogonal with respect to the inner product

$$\langle f,g\rangle =\int _{0}^{b}xf(x)g(x)\,dx$$

according to

$$\int _{0}^{1}xJ_{\alpha }(xu_{\alpha ,n})\,J_{\alpha }(xu_{\alpha ,m})\,dx={\frac {\delta _{mn}}{2}}[J_{\alpha +1}(u_{\alpha ,n})]^{2},$$

(where: $\delta _{mn}$ is the Kronecker delta). The coefficients can be obtained from projecting the function f(x) onto the respective Bessel functions:

$$c_{n}={\frac {\langle f,(J_{\alpha })_{n}\rangle }{\langle (J_{\alpha })_{n},(J_{\alpha })_{n}\rangle }}={\frac {\int _{0}^{b}xf(x)(J_{\alpha })_{n}(x)\,dx}{{\frac {1}{2}}(bJ_{\alpha \pm 1}(u_{\alpha ,n}))^{2}}}$$

where the plus or minus sign is equally valid.

One-to-one relation between order index (n) and continuous frequency ($F_{n}$)

Fourier–Bessel series coefficients are unique for a given signal, and there is one-to-one mapping between continuous frequency ($F_{n}$) and order index $(n)$ which can be expressed as follows:

$u_{n}={\frac {2\pi F_{n}L}{F_{s}}}$

Since, $u_{n}=u_{n-1}+\pi \approx n\pi$. So above equation can be rewritten as follows:

$F_{n}={\frac {F_{s}n}{2L}}$

where $L$ is the length of the signal and $F_{s}$ is the sampling frequency of the signal.

2-D- Fourier-Bessel series expansion

For an image $f(x,y)$ of size M×N, the synthesis equations for order-0 2D-Fourier–Bessel series expansion is as follows:

$f(x,y)=\sum _{m=1}^{M}\sum _{n=1}^{N}F(m,n)J_{0}{\bigg (}{\frac {u_{0,n}y}{N}}{\bigg )}J_{0}{\bigg (}{\frac {u_{0,n}x}{M}}{\bigg )}$

Where $F(m,n)$ is 2D-Fourier–Bessel series expansion coefficients whose mathematical expressions are as follows:

$F(m,n)={\frac {4}{\alpha _{1}}}\sum _{x=0}^{M-1}\sum _{y=0}^{N-1}xyf(x,y)J_{0}{\bigg (}{\frac {u_{0,n}y}{N}}{\bigg )}J_{0}{\bigg (}{\frac {u_{0,m}x}{M}}{\bigg )}$

where, $\alpha _{1}=(NM)^{2}(J_{1}(u_{0,m})J_{1}(u_{0,n}))^{2}$

Fourier-Bessel series expansion based entropies

For a signal of length $b$, Fourier-Bessel based spectral entropy such as Shannon spectral entropy ($H_{\text{SSE}}$), log energy entropy ($H_{\text{LLE}}$), and Wiener entropy ($H_{\text{WE}}$) are defined as follows:

$H_{\text{SSE}}=-\sum _{n=1}^{b}P(n)~{\text{log}}_{2}\left(P(n)\right)$

$H_{\text{WE}}=b{\frac {\sqrt {\displaystyle \prod _{n=1}^{b}E_{n}}}{\displaystyle \sum _{n=1}^{b}E_{n}}}$

$H_{\text{LE}}=-\sum _{n=1}^{b}~{\text{log}}_{2}\left(P(n)\right)$

where $P_{n}$ is is the normalized energy distribution which is mathematically defined as follows:

$P(n)={\frac {E_{n}}{\displaystyle \sum _{n=1}^{b}E_{n}}}$

$E_{n}$ is energy spectrum which is mathematically defined as follows:

$E_{n}={\frac {c_{n}^{2}b^{2}[J_{1}(u_{1,n})]^{2}}{2}}$

Fourier Bessel Series Expansion based Empirical Wavelet Transform

The Empirical wavelet transform (EWT) is a multi-scale signal processing approach for the decomposition of multi-component signal into intrinsic mode functions (IMFs) . The EWT is based on the design of empirical wavelet based filter bank based on the segregation of Fourier spectrum of the multi-component signals. The segregation of Fourier spectrum of multi-component signal is performed using the detection of peaks and then the evaluation of boundary points . For non-stationary signals, the Fourier Bessel Series Expansion (FBSE) is the natural choice as it uses Bessel function as basis for analysis and synthesis of the signal. The FBSE spectrum has produced the number of frequency bins same as the length of the signal in the frequency range [0, ${\frac {F_{s}}{2}}$]. Therefore, in FBSE-EWT, the boundary points are detected using the FBSE based spectrum of the non-stationary signal . Once, the boundary points are obtained, the empirical wavelet based filter-bank is designed in the Fourier domain of the multi-component signal to evaluate IMFs. The FBSE based method used in FBSE-EWT has produced higher number of boundary points as compared to FFT part in EWT based method . The features extracted from the IMFs of EEG and ECG signals obtained using FBSE-EWT based approach have shown better performance for the automated detection of Neurological and and cardiac ailments .

Fourier-Bessel Series Expansion Domain Discrete Stockwell Transform

For a discrete time signal, x(n), the FBSE domain discrete Stockwell transform (FBSE-DST) is evaluated as follows:$$T(n,l)=\sum _{m=1}^{L}Y{\Big (}m+l{\Big )}g(m,l)J_{0}{\Big (}{\frac {\lambda _{l}}{N}}n{\Big )}$$where Y(l) are the FBSE coefficients and these coefficients are calculated using the following expression as

$Y(l)={\frac {2}{N^{2}[J_{1}(\lambda _{l})]^{2}}}\sum _{n=0}^{N-1}nx(n)J_{0}{\Big (}{\frac {\lambda _{l}}{N}}n{\Big )}$

The $\lambda _{l}$is termed as the $l^{th}$ root of the Bessel function, and it is evaluated in an iterative manner based on the solution of $J_{0}(\lambda _{l})=0$using the Newton-Rapson method. Similarly, the g(m,l) is the FBSE domain Gaussian window and it is given as follows :

$g(m,l)={\text{e}}^{-{\frac {2\pi ^{2}\lambda _{m}^{2}}{\lambda _{l}^{2}}}},~{\{l,m=1,2,...L}\}$

Fourier–Bessel expansion-based discrete energy separation algorithm

For multicomponent amplitude and frequency modulated (AM-FM) signals, the discrete energy separation algorithm (DESA) together with the Gabor's filtering is a traditional approach to estimate the amplitude envelope (AE) and the instantaneous frequency (IF) functions<ref>Maragos, Petros; Kaiser, James F. (October 1993). "Energy Separation in Signal Modulations with Application to Speech Analysis". 1993 IEEE Transactions On Signal Processing: 3024–3051. doi:10.1109/78.277799. {{cite journal}}: line feed character in |title= at position 45 (help). It has been observed that the filtering operation distorts the amplitude and phase modulations in the separated monocomponent signals. In the work<ref>Pachori, Ram Bilas; Sircar, Pradip (January 2010). "Analysis of multicomponent AM-FM signals using FB-DESA method". 2010 Digital Signal Processing: 42–62. doi:10.1016/j.dsp.2009.04.013., the Fourier–Bessel expansion-based discrete energy separation algorithm (FB-DESA) for component separation and estimation of the AE and IF functions of a multicomponent AM-FM signal is proposed. The FB-DESA method gives accurate estimations of the AE and IF functions.

Advantages

The Fourier–Bessel series expansion does not require use of window function in order to obtain spectrum of the signal. It represents real signal in terms of real Bessel basis functions. It provides representation of real signals it terms of positive frequencies. The basis functions used are aperiodic in nature and converge. The basis functions include amplitude modulation in the representation. The Fourier–Bessel series expansion spectrum provides frequency points equal to the signal length.

Applications

The Fourier–Bessel series expansion employs aperiodic and decaying Bessel functions as the basis. The Fourier–Bessel series expansion has been successfully applied in diversified areas such as Gear fault diagnosis, discrimination of odorants in a turbulent ambient, postural stability analysis, detection of voice onset time, glottal closure instants (epoch) detection, separation of speech formants, EEG signal segmentation, speech enhancement, and speaker identification. The Fourier–Bessel series expansion has also been used to reduce cross terms in the Wigner–Ville distribution.

Dini series

A second Fourier–Bessel series, also known as Dini series, is associated with the Robin boundary condition $$bf'(b)+cf(b)=0,$$ where $c$ is an arbitrary constant. The Dini series can be defined by $$f(x)\sim \sum _{n=1}^{\infty }b_{n}J_{\alpha }(\gamma _{n}x/b),$$

where $\gamma _{n}$ is the n-th zero of $xJ'_{\alpha }(x)+cJ_{\alpha }(x)$.

The coefficients $b_{n}$ are given by $$b_{n}={\frac {2\gamma _{n}^{2}}{b^{2}(c^{2}+\gamma _{n}^{2}-\alpha ^{2})J_{\alpha }^{2}(\gamma _{n})}}\int _{0}^{b}J_{\alpha }(\gamma _{n}x/b)\,f(x)\,x\,dx.$$

See also

• Orthogonality
• Generalized Fourier series
• Hankel transform
• Kapteyn series
• Neumann polynomial
• Schlömilch's series

References

• Smythe, William R. (1968). Static and Dynamic Electricity (3rd ed.). New York: McGraw-Hill.
• Magnus, Wilhelm; Oberhettinger, Fritz; Soni, Raj Pal (1966). Formulas and Theorems for Special Functions of Mathematical Physics. Berlin: Springer.
• R.B. Pachori and P. Sircar, Non-stationary signal analysis: Methods based on Fourier-Bessel representation, LAP LAMBERT Academic Publishing, Saarbrucken, Germany, 2010, ISBN 978-3-8433-8807-8.
• J. Schroeder, Signal processing via Fourier–Bessel series expansion, Digital Signal Process. 3 (1993), 112–124.
• G. D'Elia, S. Delvecchio and G. Dalpiaz, On the use of Fourier–Bessel series expansion for gear diagnostics, Proc. of the Second Int. Conf. Condition Monitoring of Machinery in Non-Stationary Operations (2012), 267-275.
• A. Vergaraa, E. Martinelli, R. Huerta, A. D'Amico and C. Di Natale, Orthogonal decomposition of chemo-sensory signals: Discriminating odorants in a turbulent ambient, Procedia Engineering 25 (2011), 491–494.
• F.S. Gurgen and C. S. Chen, Speech enhancement by Fourier–Bessel coefficients of speech and noise, IEE Proc. Comm. Speech Vis. 137 (1990), 290–294.
• K. Gopalan, T. R. Anderson and E. J. Cupples, A comparison of speaker identification results using features based on cepstrum and Fourier–Bessel expansion, IEEE Trans. Speech Audio Process. 7 (1999), 289–294.
• R.B. Pachori and P. Sircar, Analysis of multicomponent AM-FM signals using FB-DESA method, Digital Signal Processing, vol. 20, pp. 42–62, January 2010.
• R.B. Pachori and P. Sircar, EEG signal analysis using FB expansion and second-order linear TVAR process, Signal Processing, vol. 88, no. 2, pp. 415–420, February 2008.
• R.B. Pachori and P. Sircar, A new technique to reduce cross terms in the Wigner distribution, Digital Signal Processing, vol. 17, no. 2, pp. 466–474, March 2007.
• A. Bhattacharyya, L. Singh, and R.B. Pachori, Fourier-Bessel series expansion based empirical wavelet transform for analysis of non-stationary signals, Digital Signal Processing, vol. 78, pp. 185–196, July 2018.
• P. Jain and R.B. Pachori, Event-based method for instantaneous fundamental frequency estimation from voiced speech based on eigenvalue decomposition of Hankel matrix, IEEE/ACM Transactions on Audio, Speech and Language Processing, vol. 22. issue 10, pp. 1467–1482, October 2014.
• R.K. Tripathy, A. Bhattacharyya, and R.B. Pachori, A novel approach for detection of myocardial infarction from ECG signals of multiple electrodes, IEEE Sensors Journal, DOI: 10.1109/JSEN.2019.2896308, 2019.
• V. Gupta and R.B. Pachori, Epileptic seizure identification using entropy of FBSE based EEG rhythms, Biomedical Signal Processing and Control, DOI: 10.1016/j.bspc.2019.101569, 2019.
• R. Katiyar, V. Gupta, and R.B. Pachori, FBSE-EWT based approach for the determination of respiratory rate from PPG signals, IEEE Sensors Letters, DOI: 10.1109/LSENS.2019.2926834, 2019.
• R.K. Tripathy, A. Bhattacharyya, and R.B. Pachori, Localization of myocardial infarction from multi lead electrocardiogram signals using multiscale convolution neural network, IEEE Sensors Journal, DOI: 10.1109/JSEN.2019.2935552, 2019.
• P. Gajbhiye, R.K. Tripathy, and R.B. Pachori, Elimination of ocular artifacts from single channel EEG signals using FBSE-EWT based rhythms, IEEE Sensors Journal, 2019.
• A.S. Hood, R.B. Pachori, V.K. Reddy, and P. Sircar, Parametric representation of speech employing multi-component AFM signal model, International Journal of Speech Technology, vol. 18, issue 03, pp. 287–303, September 2015.
• A. Anuragi, D. Sisodia, and RB Pachori, Automated alcoholism detection using Fourier-Bessel series expansion based empirical wavelet transform, IEEE Sensors Journal, DOI: 10.1109/JSEN.2020.2966766, 2020.
• T. Siddharth, P. Gajbhiye, R.K. Tripathy, and R.B. Pachori, EEG based detection of focal seizure area using FBSE-EWT rhythm and SAE-SVM network, IEEE Sensors Journal, DOI: 10.1109/JSEN.2020.2995749.
• V. Gupta and R.B. Pachori, Classification of focal EEG signals using FBSE based flexible time-frequency coverage wavelet transform, Biomedical Signal Processing and Control, DOI: https://doi.org/10.1016/j.bspc.2020.102124.
• A. Bhattacharyya, R.K. Tripathy, L. Garg, and R.B. Pachori, A novel multivariate-multiscale approach for computing EEG spectral and temporal complexity for human emotion recognition, IEEE Sensors Journal, DOI: 10.1109/JSEN.2020.3027181.
• P.K. Chaudhary and R.B. Pachori, Automatic diagnosis of glaucoma using two-dimensional Fourier-Bessel series expansion based empirical wavelet transform, Biomedical Signal Processing and Control, DOI: https://doi.org/10.1016/j.bspc.2020.102237.
• V. Gupta and R.B. Pachori, FBDM based time-frequency representation for sleep stages classification using EEG signals, Biomedical Signal Processing and Control, DOI: https://doi.org/10.1016/j.bspc.2020.102265.
• S.I. Khan and R.B. Pachori, Automated detection of posterior myocardial infarction from vectorcardiogram signals using Fourier Bessel series expansion based empirical wavelet transform, IEEE Sensors Letters, DOI: 10.1109/LSENS.2021.3070142.
• P.K. Chaudhary and R.B. Pachori, FBSED based automatic diagnosis of COVID-19 using X-ray and CT images, Computers in Biology and Medicine, DOI: https://doi.org/10.1016/j.compbiomed.2021.104454.
• A. Anuragi, D.S. Sisodia, and R.B. Pachori, Automated FBSE-EWT based learning framework for detection of epileptic seizures using time-segmented EEG signals, Computers in Biology and Medicine, DOI: https://doi.org/10.1016/j.compbiomed.2021.104708.
• A. Anuragi, D.S. Sisodia, and R.B. Pachori, Epileptic-seizure classification using phase-space representation of FBSE-EWT based EEG sub-band signals and ensemble learners, Biomedical Signal Processing and Control, DOI: https://doi.org/10.1016/j.bspc.2021.103138.
• P.K. Chaudhary and R.B. Pachori, Automatic diagnosis of different grades of diabetic retinopathy and diabetic macular edema using 2D-FBSE-FAWT, IEEE Transactions on Instrumentation & Measurement, DOI: 10.1109/TIM.2022.3140437.
• S.I. Khan, S.M. Qaisar, and R.B. Pachori, Automated classification of valvular heart diseases using FBSE-EWT and PSR based geometrical features, Biomedical Signal Processing and Control, DOI: https://doi.org/10.1016/j.bspc.2021.103445.
• S.I. Khan and R.B. Pachori, Derived vectorcardiogram based automated detection of posterior myocardial infarction using FBSE-EWT technique, Biomedical Signal Processing and Control, DOI: https://doi.org/10.1016/j.bspc.2021.103051.
• V. Gupta and R.B. Pachori, FB dictionary based SSBL-EM and its application for multi-class SSVEP classification using eight-channel EEG signals, IEEE Transactions on Instrumentation & Measurement, DOI: 10.1109/TIM.2022.3150848.
• A. Nalwaya, K. Das, and R.B. Pachori, Automated emotion identification using Fourier-Bessel domain-based entropies, Entropy, DOI: https://doi.org/10.3390/e24101322.
• A. Anuragi, D.S. Sisodia, and R.B. Pachori, EEG-based cross-subject emotion recognition using Fourier-Bessel series expansion based empirical wavelet transform and NCA feature selection method, Information Sciences, DOI: https://doi.org/10.1016/j.ins.2022.07.121.
• S. Dash, SK Ghosh, RK Tripathy, G. Panda, R.B. Pachori, Fourier-Bessel Domain based Discrete Stockwell Transform for the Analysis of Non-stationary Signals, IEEE India Council International Subsections conference (INDISCON), 2022, DOI: https://doi.org/10.1109/INDISCON54605.2022.9862863

External links

• "Fourier-Bessel series", Encyclopedia of Mathematics, EMS Press, 2001 [1994]
• Weisstein, Eric. W. "Fourier-Bessel Series". From MathWorld--A Wolfram Web Resource.
• Fourier–Bessel series applied to Acoustic Field analysis on Trinnov Audio's research page