High performance digital signal processing: Theory, design, and applications in finance

HIGH PERFORMANCE DIGITAL SIGNAL PROCESSING: THEORY, DESIGN, AND APPLICATIONS IN FINANCE by Mustafa Ugur Torun The way scientific research and business is conducted has drastically changed over the last decade. Big data and data-intensive scientific discovery are two terms that have been coined recently. They describe the tremendous amounts of noisy data, created extremely rapidly by various sensing devices and methods that need to be explored for information inference. Researchers and practitioners who can obtain meaningful information out of big data in the shortest time gain a competitive advantage. Hence, there is more need than ever for a variety of high performance computational tools for scientific and business analytics. Interest in developing efficient data processing methods like compression and noise filtering tools enabling real-time analytics of big data is increasing. A common concern in digital signal processing applications has been the lack of fast handling of observed data. This problem has been an active research topic being addressed by the progress in analytical tools allowing fast processing of big data. One particular tool is the Karhunen-Loève transform (KLT) (also known as the principal component analysis) where covariance matrix of a stochastic process is decomposed into its eigenvectors and eigenvalues as the optimal orthonormal transform. Specifically, eigenanalysis is utilized to determine the KLT basis functions. KLT is a widely employed signal analysis method used in applications including noise filtering of measured data and compression. However, defining KLT basis for a given signal covariance matrix demands prohibitive computational resources in many real-world scenarios. In this dissertation, engineering implementation of KLT as well as the theory of eigenanalysis for auto-regressive order one, AR(1), discrete stochastic processes are investigated and novel improvements are proposed. The new findings are applied to well-known problems in quantitative finance (QF). First, an efficient method to derive the explicit KLT kernel for AR(1) processes that utilizes a simple root finding method for the transcendental equations is introduced. Performance improvement over a popular numerical eigenanalysis algorithm, called divide and conquer, is shown. Second, implementation of parallel Jacobi algorithm for eigenanalysis on graphics processing units is improved such that the access to the dynamic random access memory is entirely coalesced. The speed is improved by a factor of 68.5 by the proposed method compared to a CPU implementation for a square matrix of size 1,024. Third, several tools developed and implemented in the dissertation are applied to QF problems such as risk analysis and portfolio risk management. In addition, several topics in QF, such as price models, Epps effect, and jump processes are investigated and new insights are suggested from a multi-resolution (multi-rate) signal processing perspective. It is expected to see this dissertation to make contributions in better understanding and bridging the analytical methods in digital signal processing and applied mathematics, and their wider utilization in the finance sector. The emerging joint research and technology development efforts in QF and financial engineering will benefit the investors, bankers, and regulators to build and maintain more robust and fair financial markets in the future. HIGH PERFORMANCE DIGITAL SIGNAL PROCESSING: THEORY, DESIGN, AND APPLICATIONS IN FINANCE by Mustafa Ugur Torun A Dissertation Submitted to the Faculty of New Jersey Institute of Technology in Partial Fulfillment of the Requirements for the Degree of Doctor of Philosophy in Electrical Engineering Department of Electrical and Computer Engineering May 2013 Copyright © 2013 by Mustafa Ugur Torun ALL RIGHTS RESERVED APPROVAL PAGE HIGH PERFORMANCE DIGITAL SIGNAL PROCESSING: THEORY, DESIGN, AND APPLICATIONS IN FINANCE Mustafa Ugur Torun Dr. Ali N. Akansu, Dissertation Advisor Date Professor of Electrical and Computer Engineering, NJIT Dr. Michael A. Ehrlich, Committee Member Date Assistant Professor of Finance, NJIT Dr. Richard A. Haddad, Committee Member Date Professor of Electrical and Computer Engineering, NJIT Dr. Sanjeev R. Kulkarni, Committee Member Date Professor of Electrical Engineering, Princeton University Dr. Osvaldo Simeone, Committee Member Date Associate Professor of Electrical and Computer Engineering, NJIT

[1]  W. Root,et al.  An introduction to the theory of random signals and noise , 1958 .

[2]  N. Ahmed,et al.  Discrete Cosine Transform , 1996 .

[3]  V. S Pugachev A method for the determination of the eigenvalues and eigenfunctions of a certain class of linear integral equations , 1959 .

[4]  V. S. Pugachev,et al.  A method of solving the basic integral equation of statistical theory of optimum systems in finite form , 1959 .

[5]  S. Heston A Closed-Form Solution for Options with Stochastic Volatility with Applications to Bond and Currency Options , 1993 .

[6]  Ali N. Akansu,et al.  A fast derivation of Karhunen-Loeve transform kernel for first-order autoregressive discrete process , 2014, PERV.

[7]  E. E. Burniston,et al.  An exact analytical solution of Kepler's equation , 1972 .

[8]  R. Clarke,et al.  Relation between the Karhunen Loève and cosine transforms , 1981 .

[9]  David R. Butenhof Programming with POSIX threads , 1993 .

[10]  G. G. Stokes "J." , 1890, The New Yale Book of Quotations.

[11]  Ali N. Akansu,et al.  Risk management for trading in multiple frequencies , 2011, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[12]  J. H. Wilkinson The algebraic eigenvalue problem , 1966 .

[13]  Ali N. Akansu,et al.  Optimal design of phase function in Generalized DFT , 2010, Phys. Commun..

[14]  Franklin T. Luk,et al.  A Proof of Convergence for Two Parallel Jacobi SVD Algorithms , 1989, IEEE Trans. Computers.

[15]  M. Melamed Detection , 2021, SETI: Astronomy as a Contact Sport.

[16]  Ali N. Akansu,et al.  FPGA based eigenfiltering for real-time portfolio risk analysis , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[17]  R. Cont,et al.  Financial Modelling with Jump Processes , 2003 .

[18]  E. E. Burniston,et al.  The use of Riemann problems in solving a class of transcendental equations , 1973, Mathematical Proceedings of the Cambridge Philosophical Society.

[19]  Gilbert Strang,et al.  Introduction to applied mathematics , 1988 .

[20]  Vedran Novakovic,et al.  A GPU-based hyperbolic SVD algorithm , 2010, ArXiv.

[21]  C. Siewert,et al.  On Solutions of a Transcendental Equation Basic to the Theory of Vibrating Plates , 1979 .

[22]  Ali N. Akansu,et al.  A novel method to derive explicit KLT kernel for AR(1) process , 2013, 2013 IEEE International Conference on Acoustics, Speech and Signal Processing.

[23]  F. Black,et al.  The Pricing of Options and Corporate Liabilities , 1973, Journal of Political Economy.

[24]  M. Avellaneda,et al.  Statistical Arbitrage in the U.S. Equities Market , 2008 .

[25]  Gordon Bell,et al.  Beyond the Data Deluge , 2009, Science.

[26]  Ali N. Akansu,et al.  On Epps effect and rebalancing of hedged portfolio in multiple frequencies , 2011, 2011 4th IEEE International Workshop on Computational Advances in Multi-Sensor Adaptive Processing (CAMSAP).

[27]  E. Isaacson,et al.  Numerical Analysis for Applied Science , 1997 .

[28]  L. Bachelier,et al.  Théorie de la spéculation , 1900 .

[29]  Ali N. Akansu,et al.  On toeplitz approximation to empirical correlation matrix of financial asset returns , 2012, 2012 46th Annual Conference on Information Sciences and Systems (CISS).

[30]  E. G. Anastasselou,et al.  A Formal Comparison of the Delves—Lyness and Burniston—Siewert Methods for Locating the Zeros of Analytic Functions , 1986 .

[31]  Kiyosi Itô On a stochastic integral equation , 1946 .

[32]  T. Alderweireld,et al.  A Theory for the Term Structure of Interest Rates , 2004, cond-mat/0405293.

[33]  Anirvan M. Sengupta,et al.  Distributions of singular values for some random matrices. , 1997, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[34]  T. W. Epps Comovements in Stock Prices in the Very Short Run , 1979 .

[35]  J. W. Brown,et al.  Complex Variables and Applications , 1985 .

[36]  Jason Sanders,et al.  CUDA by example: an introduction to general purpose GPU programming , 2010 .

[37]  Gene H. Golub,et al.  Matrix computations , 1983 .

[38]  James Demmel,et al.  Jacobi's Method is More Accurate than QR , 1989, SIAM J. Matrix Anal. Appl..

[39]  V. Plerou,et al.  Random matrix approach to cross correlations in financial data. , 2001, Physical review. E, Statistical, nonlinear, and soft matter physics.

[40]  James W. Stevens,et al.  Explicit Solutions for Transcendental Equations , 2002, SIAM Rev..

[41]  Alan V. Oppenheim,et al.  Discrete-Time Signal Pro-cessing , 1989 .

[42]  Milton Abramowitz,et al.  Handbook of Mathematical Functions with Formulas, Graphs, and Mathematical Tables , 1964 .

[43]  A. Kyle,et al.  The Flash Crash: The Impact of High Frequency Trading on an Electronic Market , 2011 .

[44]  Bruno Dupire Pricing with a Smile , 1994 .

[45]  Tony Hey,et al.  The Fourth Paradigm: Data-Intensive Scientific Discovery , 2009 .

[46]  Ali N. Akansu,et al.  On basic price model and volatility in multiple frequencies , 2011, 2011 IEEE Statistical Signal Processing Workshop (SSP).

[47]  Ian T. Jolliffe,et al.  Principal Component Analysis , 2002, International Encyclopedia of Statistical Science.

[48]  R. Haddad,et al.  Multiresolution Signal Decomposition: Transforms, Subbands, and Wavelets , 1992 .

[49]  W. Hackbusch Singular Integral Equations , 1995 .

[50]  Aarnout Brombacher,et al.  Probability... , 2009, Qual. Reliab. Eng. Int..

[51]  J. Bouchaud,et al.  RANDOM MATRIX THEORY AND FINANCIAL CORRELATIONS , 2000 .

[52]  Mary Hall,et al.  Takagi Factorization on GPU using CUDA , 2011 .

[53]  W. D. Ray,et al.  Further decomposition of the Karhunen-Loève series representation of a stationary random process , 1970, IEEE Trans. Inf. Theory.

[54]  Lan Zhang Estimating Covariation: Epps Effect, Microstructure Noise , 2006 .

[55]  Andrew Rennie,et al.  Financial Calculus: An Introduction to Derivative Pricing , 1996 .

[56]  Ali N. Akansu,et al.  Novel GPU implementation of Jacobi algorithm for Karhunen-Loève transform of dense matrices , 2012, 2012 46th Annual Conference on Information Sciences and Systems (CISS).

[57]  A. Stuart,et al.  Portfolio Selection: Efficient Diversification of Investments , 1959 .

[58]  C. Jacobi,et al.  C. G. J. Jacobi's Gesammelte Werke: Über ein leichtes Verfahren, die in der Theorie der Sacularstorungen vorkommenden Gleichungen numerisch aufzulosen , 1846 .

[59]  D. Owen Handbook of Mathematical Functions with Formulas , 1965 .

[60]  Ali N. Akansu,et al.  Implementation of generalized DFT on Field Programmable Gate Array , 2012, 2012 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[61]  Ali N. Akansu,et al.  A novel GPU implementation of eigenanalysis for risk management , 2012, 2012 IEEE 13th International Workshop on Signal Processing Advances in Wireless Communications (SPAWC).

[62]  Kalyanmoy Deb,et al.  Multi-objective optimization using evolutionary algorithms , 2001, Wiley-Interscience series in systems and optimization.

[63]  Alan G. White,et al.  The Pricing of Options on Assets with Stochastic Volatilities , 1987 .

[64]  Robert M. Farber,et al.  CUDA Application Design and Development , 2011 .

[65]  J. Miller Numerical Analysis , 1966, Nature.

[66]  Athanasios Papoulis,et al.  Probability, Random Variables and Stochastic Processes , 1965 .

[67]  Ali N. Akansu,et al.  Portfolio Risk in Multiple Frequencies , 2011, IEEE Signal Processing Magazine.

[68]  Steven Kay,et al.  Modern Spectral Estimation: Theory and Application , 1988 .

[69]  Ali N. Akansu,et al.  Toeplitz Approximation to Empirical Correlation Matrix of Asset Returns: A Signal Processing Perspective , 2012, IEEE Journal of Selected Topics in Signal Processing.

[70]  Robert A. Leathers,et al.  Closed-Form Solutions for Transcendental Equations of Heat Transfer , 1996 .

[71]  Emmanuel Bacry,et al.  Modeling microstructure noise using Hawkes processes , 2011, 2011 IEEE International Conference on Acoustics, Speech and Signal Processing (ICASSP).

[72]  Ulf Schlichtmann,et al.  A Successive Approach to Compute the Bounded Pareto Front of Practical Multiobjective Optimization Problems , 2009, SIAM J. Optim..