NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL. Georgeta Budura
|
|
- Aldous Spencer
- 5 years ago
- Views:
Transcription
1 NONLINEAR SYSTEMS IDENTIFICATION USING THE VOLTERRA MODEL Georgeta Budura Politenica University of Timisoara, Faculty of Electronics and Telecommunications, Comm. Dep., Abstract: Nonlinear adaptive filtering tecniques are widely used for te nonlinearities identification in manny applications. Tis paper investigates te performances of te Volterra estimator by considering a nonlinear system identification application. Te Volterra estimator parameters are compared wit tose of a linear estimator. For te nonlinear estimator, based on a second order RLS Volterra filter, a new implementation is proposed. Te experimental results sow tat te proposed Volterra identifier as better performances tan te linear one. Different degrees of nonlinearity for te nonlinear system are considered. Keywords: Recursive least squares algoritm, nonlinear estimator, linear estimator, Volterra filter.. INTRODUCTION Detection, representation and identification of nonlinearities in different telecommunication systems represent important tasks in many applications and ad a major contribution to te development of te main nonlinear modeling tecniques. Te current trend in te telecommunication systems design is te identification and compensation of unwanted nonlinearities. It was demonstrated tat unwanted nonlinearities in te system will ave a determinant effect on is performance (Tsimbinos, 995). Tere are various ways of reducing te effects of undesired nonlinearities (Stenger, 000 a,b; Küc, 00 ). Te use of nonlinear models considered in tis paper to caracterize and compensate armful nonlinearities offer a possible solution. Te Volterra series ave been widely applied as nonlinear system modeling tecnique wit considerable success. However, at present, none general metod exists to calculate te Volterra kernels for nonlinear systems, altoug tey can be calculated for systems wose order is known and finite. Wen te nonlinear system order is unknown, adaptive metods and algoritms are widely used for te Volterra kernel estimation. Te accuracy of te Volterra kernels will determine te accuracy of te system model and te accuracy of te inverse system used for compensation. Te speed of kernel estimation process is also important. A fast kernel estimation metod may allow te user to construct a iger order model tat give an even better system representation. Tis paper investigates te performances of te Volterra estimator by considering a nonlinear system identification application. Te Volterra estimator parameters are compared wit tose of a linear estimator. Te nonlinear estimator is based on a second order RLS Volterra filter. A new implementation of te second order RLS Volterra filter based on te extended input vector and on te extended filter coefficients vector is also proposed. Due to te linearity of te input-output relation of te Volterra model wit respect to filter coefficients, te implementation of te RLS algoritm was realized as an extension of te RLS algoritm for linear filters.. THE VOLTERRA MODEL Tis section will discuss some important aspects of te second order Volterra model. For a discrete-time and causal nonlinear system wit memory, wit input x[n] and output y[n], te Volterra series expansion is given by (Scetzen, 980) :
2 N M y [ n] =... r[ k,..., kr ] x[ n k]... x[ n kr ] () r = k = 0 M kr = 0 were N represent te model nonlinearity degree. Coosing N =, te input-output relationsip of te second order Volterra filter (FV ) can be expressed as: y M [ n] = 0 [ k] x[ n k] k 0 M M k = 0 k = 0 = () [ k, k ] x[ n k ] x[ n k ] Te input-output relation can also be written in terms of nonlinear operators as indicated in relation (3). [ n] H [ x[ n ] H [ x[ n ] H [ x[ n ] y 0 = (3) In te above representations, te functions i, i = 0, represent te kernels associated to te nonlinear operators [ x[ n ] H i. Te nonlinear model described by te relations () and (3) is called a second order Volterra model. Note tat te above representations as te same memory for all nonlinearity orders. In te most general case te relation () may used different memory for eac nonlinearity order. A furter simplification can be made to relation () by considering symmetric Volterra kernels. Te kernel i [ k, Λ k i ] is symmetric if te indices can be intercanged witout affecting its value. Te second order Volterra kernel is a (M M) matrix. If we consider symmetric kernels of memory M, te second order Volterra kernel requires te determination of M(M)/ coefficients as indicated in relation (4). M M = k = 0 k = k [ k, k ] x[ n k ] x[ n ] y [ n] k (4) Te kernel estimation accuracy becomes te major problem in te practical applications. It was sown tat te Volterra operators are omogeneous and generally not ortogonal. As a consequence of tis last caracteristic te Volterra kernels can not be measured using te cross correlation tecniques and te values of te Volterra kernels will depend on te order of te Volterra representation being used. If te order of te Volterra model is canged te Volterra kernels will cange and tey must be recalculated (Scetzen, 980; Budura and Nafornita 00; Budura and Botoca, 00a). However, for an input aving a symmetric amplitude density function, suc as te Gaussian noise, te odd order Volterra functionals are ortogonal to te even order Volterra functionals. It follows tat for tis type of input, a nd order Volterra model, wit zero DC component, is an ortogonal model (Budura and Botoca, 005). Tis leads to direct Volterra kernel measurement by te cross-correlation metods. In practical situations wen confronted wit system wic parameters are time varying, adaptive metods for kernels estimation are widely used. Due to te linearity of te input-output relation according to te kernels, or filter coefficients, te application of adaptive algoritms for te Volterra filters implementation is quite simple. Te nonlinearity is reflected only by multiple products between te delayed versions of te input signal. Te nonlinear estimator proposed in te article is based on a second order Volterra filter (FV ) descriebed by relation ()-(4) and considering 0 = 0. Te filter structure is indicated in Fig.. x[n] Fig.. Te second order Volterra filter structure Next we will introduce te input vectors corresponding to different orders kernels. Te first order input vector is defined as follows: [ x[ n] x[ n... ] x[ n ] X (5) = M Te second order input vector can be expressed by: X * X X = (6) Wen we deal wit symmetric Volterra kernels only M(M)/ elements of te MxM matrix indicated in relation (6) are used in te input output relation (). Based on tis observation we introduced te extended input vector for te second order Volterra filter as indicated in relation (7): [ n]... x[ n M ] x [ n] x[ n] x[ n M ] X = x... x [ n M ] [k ] [k, k ] (7) and te extended filter coefficients vector according to (8). [ ] H (8) = 0... M M M M M 3. RLS VOLTERRA ESTIMATOR Adaptive metods and algoritms are widely used for te purpose of kernel estimation. A Volterra filter of FV y V [n]
3 fixed order and fixed memory adapts to te unknown nonlinear system using one of te various adaptive algoritms. Te use of adaptive tecniques for Volterra kernel estimation as been well publised. Most of te previous work considers nd order Volterra filters and some consider te 3 rd order case. A simple and commonly used algoritm uses an LMS adaptation criterion (Matews, 99; Budura and Botoca, 00b; Budura and Botoca 004). Adaptive Volterra filters based on te LMS adaptation algoritm are computational simple but suffer from slow and input signal dependent convergence beavior and ence are not useful in many applications. Te aim of tis section is to discuss te efficient implementation of te RLS adaptive algoritm on a second order Volterra filter. As in te linear case te adaptive nonlinear system minimizes te following cost function at eac time: J n n k T ( n) = λ y( k) H ( n) X ( k) k = 0 were H ( n) and ( n) (9) X are te coefficients and te input signal vectors, respectively, as defined in (8) and (7), is a factor tat controls te memory span of te adaptive filter and y( k) represents te desired output. Te solution of equation (8) can be obtained recursively using te RLS algoritm. Based on te second order Volterra filter presented in Section we ave implemented te RLS nonlinear estimator. Te filter coefficients are adapted according to te following steps: I. Initialization: -define te filter memory(lengt for H(n) and X(n)); ( 0 ) = [ ] H ; ( 0) δ I C = * ; II. Operations: for C n =, nr. of iterations. Create te extended input vector: X ( n) ;. Compute te error: e ( n / n ) = d( n) H ( n ) * X ( n) ; 3. Compute te scalar: ( n) = X ( n) * C ( n ) * X ( n) µ ; 4. Compute te matrix: ( n) = C ( n ) * H ( n ) /( λ µ ); G 5. Updates te filter vector: H ( n) H ( n ) e( n / n ) * G ( n) = ; 6. Updates te matrix C : ( n) = λ *( C ( n ) G( n) * X ( n) * C ( n ) ) ; In te relations above C denotes te inverse autocorrelation matrix of te extended input signal. Inversion was done according to te matrix inversion lemma (Mortensen, 987). 4. EXPERIMENTS AND RESULTS Te RLS Volterra estimator performances were studied and compared wit te linear estimator performances in a typical nonlinear system identification application presented in figure. were is a small positive constant; FV y V [n] _ Noise generator z[n] FIR x [n] Nonlinear system y[n] Adaptive linear filter y l [n] _ Fig.. Te nonlinear system identification
4 Te nonlinear system wit memory being identified is represented in figure 3 and consists of a linear filter(fir) wit impulse response given by: ( ) [ n] = 0 n M n (0) followed by a nonlinear system witout memory wic input-output relation is: y [ n] = u[ n] bu [ n] () Te linear filter memory in relation (0) is M=0. Te coeficient b permits to cange te nonlinearity degree. Te input signal x[n] is generated by colouring te wite Gaussian sequence z[n], wit te autoregresive filter (FIR) described by: [ n] x[ n ] 0,9x[ n] 0, z[ n] x = 5 () Fig.4. Te error of te RLS Volterra identifier b = 0.0 Te error using te linear adaptive filter as identifier is presented in Fig.5. Experiments ave been done regarding te identification of a second order nonlinear system wit different degrees of nonlinearity. x[n] FIR u[n] Nonlinear system witout memory y[n] Nonlinear system Fig. 3. Te nonlinear system structure Te performance of te RLS adaptive Volterra filter was appreciated by comparing te error of te nonlinear identifier wit te error of a linear identifier. Te same memory M=0 as been cosen for te linear and for te second order Volterra kernel according to te relation (). For te linear identifier we fixed te same memory. Te factor was cosen equal to 0,995. Te simulations ave been done in te MATLAB software. Te following degrees of nonlinearity are considered: -Nonlinear system identified: b=0,0; Te error using te RLS Volterra identifier is indicated in Fig.4. Fig.5. Te error of te linear identifier b = 0, 0 -Nonlinear system identified: b=0,; Te error using te RLS Volterra identifier is indicated in Fig.6 and te error using te linear identifier is depicted in Fig.7. Fig.6. Te error of te RLS Volterra identifier b = 0.
5 Table Te identifiers parameters b m el σ el m ev σ ev 0, , , 0,03 0,068 6, ,353 0,675 7,50-5 0, CONCLUSIONS Fig.7. Te error of te RLS linear identifier b = 0. -Nonlinear system identified: b=; Te error signal using te RLS Volterra identifier is indicated in Fig.8 and te error using te linear identifier is depicted in Fig.9. Fig.8. Te error of te RLS Volterra identifier b = Te identification and compensation of unwanted nonlinearities are required in manny practical application in order to improve te system performances. Usualy tis are week nonlinearities and can be well represented wit te Volterra model. Tis paper investigates te performances of te Volterra estimator by considering a nonlinear system identification application. Te nonlinear identifier is based on te RLS Volterra filter. A new implementation of te second order RLS Volterra filter based on te extended input vector and on te extended filter coefficients vector is also proposed. Due to te linearity of te input-output relation of te Volterra model wit respect to filter coefficients, te implementation of te RLS algoritm was realized as an extension of te RLS algoritm for linear filters. For simplicity, tere were considered only second order nonlinearities, but te proposed tecnique can be extented to iger order nonlinearities. Te nonlinear adaptive filter performances were evaluated in a typical system identification application and compared wit te performances of a linear identifier. Te experimental results sowed tat te RLS Volterra filter performed better tan te linear filter wic performances were inacceptable wen te nonlinearity degree was increased. Te costs of tese performances are paid by te computational complexity required by te nonlinear adaptive Volterra filter implementation. REFERENCES Fig. 9. Te error of te RLS linear identifier b = Te mean and te dispersion of te error signals were calculated too, in order to caracterize te performances of bot identifiers. Te corresponding values are indicated in Table. As it can be seen te nonlinear identifier performances ( mev, σ ev ) are very good for different nonlinearity degrees, wile te linear identifier performances ( mel, σ el, ) are unsatisfactory wen te nonlinearity increases. Budura Georgeta., I. Nafornia, (00),"Kernels Measurement Tecniques for Constructing Nonlinear Models", Scientific Bulletin of te Politenica" University of Timioara, Romania, Vol.I.,Timioara, pag Budura Georgeta, C.Botoca (00a), Analyses and Modeling of Systems wit Nonlinearities, Proceedings of te Symposium on Electronics and Telecomm. `ETC 00`, Vol.I, Timioara, pp Budura Georgeta, C.Botoca (00b), "Applications of te Volterra Models in Nonlinear Systems Identification", Proceedings of te Symposium on Electronics and Telecomm. `ETC 00`, Vol.I, Timioara, pag Budura Georgeta, C. Botoca (004), "Efficient Implementation of Te Second Order Volterra Filter", Scientific Bulletin of te Politenica" University of Timioara, Romania, Transactions
6 on Electronics and Communications, 49 tom, Fascicle, pag Budura Georgeta., C.Botoca, (005) "La construction dun modele non lineaire a laide de series Volterra et Wiener", accepted to be publised in Revue de lacademie roumaine, Bucarest. Küc, F.,W. Kellermann, (00), Nonlinear Eco Cancellation using a Second Order Volterra Filter, ICCASSP 00, {kuec,wk}@lnt.de. Matews, V. J. (99), "Adaptive Polynomial Filters", IEEE, Signal Processing Magazine, July pp Mortensen R.E., (987), "Random Signals and Systems", New York, Jon Wiley & Sons, 05 pages. Scetzen M. (980), "Te Volterra and Wiener teories of nonlinear systems", New York, Jon Wiley and Sons,. Stenger A., W. Kellermann (000), RLS-Adapted Polynomial for Nonlinear Acoustic Eco Cancelling, Proceeedings of te X EUSIPCO 000, Tampere, Finland, pp Stenger A., W. Kellermann (000), Adaptation of a Memoryless Preprocesor for Nonlinear Acoustic Eco Cancelling, Signal Processing, No. 9, 000, pp Tsimbinos J. (995), "Identification And Compensation Of Nonlinear Distortion", Tesis, Institute for Telecommunications Researc, University of Sout Australia, ttp:// 37 pages.
Numerical Differentiation
Numerical Differentiation Finite Difference Formulas for te first derivative (Using Taylor Expansion tecnique) (section 8.3.) Suppose tat f() = g() is a function of te variable, and tat as 0 te function
More informationDigital Filter Structures
Digital Filter Structures Te convolution sum description of an LTI discrete-time system can, in principle, be used to implement te system For an IIR finite-dimensional system tis approac is not practical
More informationThese errors are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More informationHOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS
HOW TO DEAL WITH FFT SAMPLING INFLUENCES ON ADEV CALCULATIONS Po-Ceng Cang National Standard Time & Frequency Lab., TL, Taiwan 1, Lane 551, Min-Tsu Road, Sec. 5, Yang-Mei, Taoyuan, Taiwan 36 Tel: 886 3
More informationch (for some fixed positive number c) reaching c
GSTF Journal of Matematics Statistics and Operations Researc (JMSOR) Vol. No. September 05 DOI 0.60/s4086-05-000-z Nonlinear Piecewise-defined Difference Equations wit Reciprocal and Cubic Terms Ramadan
More informationCopyright IEEE, 13th Stat. Signal Proc. Workshop (SSP), July 2005, Bordeaux (F)
Copyrigt IEEE, 13t Stat. Signal Proc. Worksop SSP), 17-0 July 005, Bordeaux F) ROBUSTNESS ANALYSIS OF A GRADIENT IDENTIFICATION METHOD FOR A NONLINEAR WIENER SYSTEM Ernst Ascbacer, Markus Rupp Institute
More informationChapter 5 FINITE DIFFERENCE METHOD (FDM)
MEE7 Computer Modeling Tecniques in Engineering Capter 5 FINITE DIFFERENCE METHOD (FDM) 5. Introduction to FDM Te finite difference tecniques are based upon approximations wic permit replacing differential
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximatinga function fx, wose values at a set of distinct points x, x, x,, x n are known, by a polynomial P x suc
More informationA MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES
A MONTE CARLO ANALYSIS OF THE EFFECTS OF COVARIANCE ON PROPAGATED UNCERTAINTIES Ronald Ainswort Hart Scientific, American Fork UT, USA ABSTRACT Reports of calibration typically provide total combined uncertainties
More informationLIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION
LIMITATIONS OF EULER S METHOD FOR NUMERICAL INTEGRATION LAURA EVANS.. Introduction Not all differential equations can be explicitly solved for y. Tis can be problematic if we need to know te value of y
More informationA = h w (1) Error Analysis Physics 141
Introduction In all brances of pysical science and engineering one deals constantly wit numbers wic results more or less directly from experimental observations. Experimental observations always ave inaccuracies.
More informationThe Verlet Algorithm for Molecular Dynamics Simulations
Cemistry 380.37 Fall 2015 Dr. Jean M. Standard November 9, 2015 Te Verlet Algoritm for Molecular Dynamics Simulations Equations of motion For a many-body system consisting of N particles, Newton's classical
More informationLecture 15. Interpolation II. 2 Piecewise polynomial interpolation Hermite splines
Lecture 5 Interpolation II Introduction In te previous lecture we focused primarily on polynomial interpolation of a set of n points. A difficulty we observed is tat wen n is large, our polynomial as to
More information1 The concept of limits (p.217 p.229, p.242 p.249, p.255 p.256) 1.1 Limits Consider the function determined by the formula 3. x since at this point
MA00 Capter 6 Calculus and Basic Linear Algebra I Limits, Continuity and Differentiability Te concept of its (p.7 p.9, p.4 p.49, p.55 p.56). Limits Consider te function determined by te formula f Note
More informationHARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS
HARMONIC ALLOCATION TO MV CUSTOMERS IN RURAL DISTRIBUTION SYSTEMS V Gosbell University of Wollongong Department of Electrical, Computer & Telecommunications Engineering, Wollongong, NSW 2522, Australia
More informationlecture 26: Richardson extrapolation
43 lecture 26: Ricardson extrapolation 35 Ricardson extrapolation, Romberg integration Trougout numerical analysis, one encounters procedures tat apply some simple approximation (eg, linear interpolation)
More informationCORRELATION TEST OF RESIDUAL ERRORS IN FREQUENCY DOMAIN SYSTEM IDENTIFICATION
IAC Symposium on System Identification, SYSID 006 Marc 9-3 006, Newcastle, Australia CORRELATION TEST O RESIDUAL ERRORS IN REQUENCY DOMAIN SYSTEM IDENTIICATION István Kollár *, Ri Pintelon **, Joan Scouens
More informationCombining functions: algebraic methods
Combining functions: algebraic metods Functions can be added, subtracted, multiplied, divided, and raised to a power, just like numbers or algebra expressions. If f(x) = x 2 and g(x) = x + 2, clearly f(x)
More informationSymmetry Labeling of Molecular Energies
Capter 7. Symmetry Labeling of Molecular Energies Notes: Most of te material presented in tis capter is taken from Bunker and Jensen 1998, Cap. 6, and Bunker and Jensen 2005, Cap. 7. 7.1 Hamiltonian Symmetry
More informationThese error are made from replacing an infinite process by finite one.
Introduction :- Tis course examines problems tat can be solved by metods of approximation, tecniques we call numerical metods. We begin by considering some of te matematical and computational topics tat
More information5 Ordinary Differential Equations: Finite Difference Methods for Boundary Problems
5 Ordinary Differential Equations: Finite Difference Metods for Boundary Problems Read sections 10.1, 10.2, 10.4 Review questions 10.1 10.4, 10.8 10.9, 10.13 5.1 Introduction In te previous capters we
More informationPolynomial Interpolation
Capter 4 Polynomial Interpolation In tis capter, we consider te important problem of approximating a function f(x, wose values at a set of distinct points x, x, x 2,,x n are known, by a polynomial P (x
More information5.1 We will begin this section with the definition of a rational expression. We
Basic Properties and Reducing to Lowest Terms 5.1 We will begin tis section wit te definition of a rational epression. We will ten state te two basic properties associated wit rational epressions and go
More informationComputer Derivations of Numerical Differentiation Formulae. Int. J. of Math. Education in Sci. and Tech., V 34, No 2 (March-April 2003), pp
Computer Derivations o Numerical Dierentiation Formulae By Jon H. Matews Department o Matematics Caliornia State University Fullerton USA Int. J. o Mat. Education in Sci. and Tec. V No (Marc-April ) pp.8-87.
More informationBayesian ML Sequence Detection for ISI Channels
Bayesian ML Sequence Detection for ISI Cannels Jill K. Nelson Department of Electrical and Computer Engineering George Mason University Fairfax, VA 030 Email: jnelson@gmu.edu Andrew C. Singer Department
More informationConsider a function f we ll specify which assumptions we need to make about it in a minute. Let us reformulate the integral. 1 f(x) dx.
Capter 2 Integrals as sums and derivatives as differences We now switc to te simplest metods for integrating or differentiating a function from its function samples. A careful study of Taylor expansions
More informationLTI Discrete-Time Systems in Transform Domain Ideal Filters Zero Phase Transfer Functions Linear Phase Transfer Functions
LTI Discrete-Time Systems in Transform Domain Ideal Filters Zero Pase Transfer Functions Linear Pase Transfer Functions Tania Stataki 811b t.stataki@imperial.ac.uk Types of Transfer Functions Te time-domain
More informationMulti-User Communication: Capacity, Duality, and Cooperation. Nihar Jindal Stanford University Dept. of Electrical Engineering
Multi-User Communication: Capacity, Duality, and Cooperation Niar Jindal Stanford University Dept. of Electrical Engineering February 3, 004 Wireless Communication Vision Cellular Networks Wireless LAN
More informationDifferential Calculus (The basics) Prepared by Mr. C. Hull
Differential Calculus Te basics) A : Limits In tis work on limits, we will deal only wit functions i.e. tose relationsips in wic an input variable ) defines a unique output variable y). Wen we work wit
More informationConvergence Analysis of Two-Dimensional LMS FIR Filters
Convergence Analysis of Two-Dimensional LMS FIR Filters Maa Sadayde Masayuki Kawamata Graduate Scool of Engineering, Tooku University Aoba, Aramaki, Aoba-ku, Sendai 980 77, APAN maa@mk.ecei.tooku.ac.jp
More informationarxiv: v1 [math.oc] 18 May 2018
Derivative-Free Optimization Algoritms based on Non-Commutative Maps * Jan Feiling,, Amelie Zeller, and Cristian Ebenbauer arxiv:805.0748v [mat.oc] 8 May 08 Institute for Systems Teory and Automatic Control,
More informationTe comparison of dierent models M i is based on teir relative probabilities, wic can be expressed, again using Bayes' teorem, in terms of prior probab
To appear in: Advances in Neural Information Processing Systems 9, eds. M. C. Mozer, M. I. Jordan and T. Petsce. MIT Press, 997 Bayesian Model Comparison by Monte Carlo Caining David Barber D.Barber@aston.ac.uk
More informationMATH1151 Calculus Test S1 v2a
MATH5 Calculus Test 8 S va January 8, 5 Tese solutions were written and typed up by Brendan Trin Please be etical wit tis resource It is for te use of MatSOC members, so do not repost it on oter forums
More informationThe total error in numerical differentiation
AMS 147 Computational Metods and Applications Lecture 08 Copyrigt by Hongyun Wang, UCSC Recap: Loss of accuracy due to numerical cancellation A B 3, 3 ~10 16 In calculating te difference between A and
More informationSECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY
(Section 3.2: Derivative Functions and Differentiability) 3.2.1 SECTION 3.2: DERIVATIVE FUNCTIONS and DIFFERENTIABILITY LEARNING OBJECTIVES Know, understand, and apply te Limit Definition of te Derivative
More informationChapter Seven The Quantum Mechanical Simple Harmonic Oscillator
Capter Seven Te Quantum Mecanical Simple Harmonic Oscillator Introduction Te potential energy function for a classical, simple armonic oscillator is given by ZÐBÑ œ 5B were 5 is te spring constant. Suc
More informationFast Exact Univariate Kernel Density Estimation
Fast Exact Univariate Kernel Density Estimation David P. Hofmeyr Department of Statistics and Actuarial Science, Stellenbosc University arxiv:1806.00690v2 [stat.co] 12 Jul 2018 July 13, 2018 Abstract Tis
More informationRobust Parametric Modeling of Speech in Additive White Gaussian Noise
Journal of Signal and Information Processing, 05, 6, 99-08 Publised Online May 05 in SciRes. ttp://www.scirp.org/journal/jsip ttp://d.doi.org/0.436/jsip.05.600 Robust Parametric Modeling of Speec in Additive
More informationarxiv: v3 [math.na] 15 Dec 2009
THE NAVIER-STOKES-VOIGHT MODEL FOR IMAGE INPAINTING M.A. EBRAHIMI, MICHAEL HOLST, AND EVELYN LUNASIN arxiv:91.4548v3 [mat.na] 15 Dec 9 ABSTRACT. In tis paper we investigate te use of te D Navier-Stokes-Voigt
More informationEfficient algorithms for for clone items detection
Efficient algoritms for for clone items detection Raoul Medina, Caroline Noyer, and Olivier Raynaud Raoul Medina, Caroline Noyer and Olivier Raynaud LIMOS - Université Blaise Pascal, Campus universitaire
More information1 Calculus. 1.1 Gradients and the Derivative. Q f(x+h) f(x)
Calculus. Gradients and te Derivative Q f(x+) δy P T δx R f(x) 0 x x+ Let P (x, f(x)) and Q(x+, f(x+)) denote two points on te curve of te function y = f(x) and let R denote te point of intersection of
More informationMATH745 Fall MATH745 Fall
MATH745 Fall 5 MATH745 Fall 5 INTRODUCTION WELCOME TO MATH 745 TOPICS IN NUMERICAL ANALYSIS Instructor: Dr Bartosz Protas Department of Matematics & Statistics Email: bprotas@mcmasterca Office HH 36, Ext
More information1 1. Rationalize the denominator and fully simplify the radical expression 3 3. Solution: = 1 = 3 3 = 2
MTH - Spring 04 Exam Review (Solutions) Exam : February 5t 6:00-7:0 Tis exam review contains questions similar to tose you sould expect to see on Exam. Te questions included in tis review, owever, are
More informationSolving Continuous Linear Least-Squares Problems by Iterated Projection
Solving Continuous Linear Least-Squares Problems by Iterated Projection by Ral Juengling Department o Computer Science, Portland State University PO Box 75 Portland, OR 977 USA Email: juenglin@cs.pdx.edu
More informationPreconditioning in H(div) and Applications
1 Preconditioning in H(div) and Applications Douglas N. Arnold 1, Ricard S. Falk 2 and Ragnar Winter 3 4 Abstract. Summarizing te work of [AFW97], we sow ow to construct preconditioners using domain decomposition
More informationNumerical Analysis MTH603. dy dt = = (0) , y n+1. We obtain yn. Therefore. and. Copyright Virtual University of Pakistan 1
Numerical Analysis MTH60 PREDICTOR CORRECTOR METHOD Te metods presented so far are called single-step metods, were we ave seen tat te computation of y at t n+ tat is y n+ requires te knowledge of y n only.
More informationMaterial for Difference Quotient
Material for Difference Quotient Prepared by Stepanie Quintal, graduate student and Marvin Stick, professor Dept. of Matematical Sciences, UMass Lowell Summer 05 Preface Te following difference quotient
More informationMathematics 5 Worksheet 11 Geometry, Tangency, and the Derivative
Matematics 5 Workseet 11 Geometry, Tangency, and te Derivative Problem 1. Find te equation of a line wit slope m tat intersects te point (3, 9). Solution. Te equation for a line passing troug a point (x
More information158 Calculus and Structures
58 Calculus and Structures CHAPTER PROPERTIES OF DERIVATIVES AND DIFFERENTIATION BY THE EASY WAY. Calculus and Structures 59 Copyrigt Capter PROPERTIES OF DERIVATIVES. INTRODUCTION In te last capter you
More informationDefinition of the Derivative
Te Limit Definition of te Derivative Tis Handout will: Define te limit grapically and algebraically Discuss, in detail, specific features of te definition of te derivative Provide a general strategy of
More informationTechnology-Independent Design of Neurocomputers: The Universal Field Computer 1
Tecnology-Independent Design of Neurocomputers: Te Universal Field Computer 1 Abstract Bruce J. MacLennan Computer Science Department Naval Postgraduate Scool Monterey, CA 9393 We argue tat AI is moving
More informationThe derivative function
Roberto s Notes on Differential Calculus Capter : Definition of derivative Section Te derivative function Wat you need to know already: f is at a point on its grap and ow to compute it. Wat te derivative
More informationEffect of the Dependent Paths in Linear Hull
1 Effect of te Dependent Pats in Linear Hull Zenli Dai, Meiqin Wang, Yue Sun Scool of Matematics, Sandong University, Jinan, 250100, Cina Key Laboratory of Cryptologic Tecnology and Information Security,
More informationESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT
ESTIMATION OF TIME-DOMAIN FREQUENCY STABILITY BASED ON PHASE NOISE MEASUREMENT P. C. Cang, H. M. Peng, and S. Y. Lin National Standard Time & Frequenc Laborator, TL, Taiwan, Lane 55, Min-Tsu Road, Sec.
More informationService Outage Based Power and Rate Allocation
1 Service Outage Based Power and Rate Allocation Jiangong Luo, Lang Lin, Roy Yates, and Predrag Spasojević Abstract Tis paper combines te concepts of ergodic capacity and capacity versus outage for fading
More informationAN IMPROVED WEIGHTED TOTAL HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES
AN IMPROVED WEIGHTED TOTA HARMONIC DISTORTION INDEX FOR INDUCTION MOTOR DRIVES Tomas A. IPO University of Wisconsin, 45 Engineering Drive, Madison WI, USA P: -(608)-6-087, Fax: -(608)-6-5559, lipo@engr.wisc.edu
More informationProfit Based Unit Commitment in Deregulated Electricity Markets Using A Hybrid Lagrangian Relaxation - Particle Swarm Optimization Approach
Profit Based Unit Commitment in Deregulated Electricity Markets Using A Hybrid Lagrangian Relaxation - Particle Swarm Optimization Approac Adline K. Bikeri, Cristoper M. Maina and Peter K. Kiato Proceedings
More informationTaylor Series and the Mean Value Theorem of Derivatives
1 - Taylor Series and te Mean Value Teorem o Derivatives Te numerical solution o engineering and scientiic problems described by matematical models oten requires solving dierential equations. Dierential
More informationHarmonic allocation to MV customers in rural distribution systems
University of Wollongong Researc Online Faculty of Engineering - Papers (Arcive) Faculty of Engineering and Information Sciences 2007 Harmonic allocation to MV customers in rural distribution systems Victor
More informationDedicated to the 70th birthday of Professor Lin Qun
Journal of Computational Matematics, Vol.4, No.3, 6, 4 44. ACCELERATION METHODS OF NONLINEAR ITERATION FOR NONLINEAR PARABOLIC EQUATIONS Guang-wei Yuan Xu-deng Hang Laboratory of Computational Pysics,
More informationOn the Identifiability of the Post-Nonlinear Causal Model
UAI 9 ZHANG & HYVARINEN 647 On te Identifiability of te Post-Nonlinear Causal Model Kun Zang Dept. of Computer Science and HIIT University of Helsinki Finland Aapo Hyvärinen Dept. of Computer Science,
More informationUnilateral substitutability implies substitutable completability in many-to-one matching with contracts
Unilateral substitutability implies substitutable completability in many-to-one matcing wit contracts Sangram Vilasrao Kadam April 6, 2015 Abstract We prove tat te unilateral substitutability property
More informationSimulation and verification of a plate heat exchanger with a built-in tap water accumulator
Simulation and verification of a plate eat excanger wit a built-in tap water accumulator Anders Eriksson Abstract In order to test and verify a compact brazed eat excanger (CBE wit a built-in accumulation
More informationarxiv: v1 [cs.sd] 28 Feb 2017
Nonlinear Volterra Model of a Loudspeaker Behavior Based on Laser Doppler Vibrometry Alessandro Loriga, Parvin Moyassari, and Daniele Bernardini Intranet Standard GmbH, Ottostrasse 3, 80333 Munich, Germany
More informationIntroduction to Derivatives
Introduction to Derivatives 5-Minute Review: Instantaneous Rates and Tangent Slope Recall te analogy tat we developed earlier First we saw tat te secant slope of te line troug te two points (a, f (a))
More informationSolutions to the Multivariable Calculus and Linear Algebra problems on the Comprehensive Examination of January 31, 2014
Solutions to te Multivariable Calculus and Linear Algebra problems on te Compreensive Examination of January 3, 24 Tere are 9 problems ( points eac, totaling 9 points) on tis portion of te examination.
More informationFast Explicit and Unconditionally Stable FDTD Method for Electromagnetic Analysis Jin Yan, Graduate Student Member, IEEE, and Dan Jiao, Fellow, IEEE
Tis article as been accepted for inclusion in a future issue of tis journal. Content is final as presented, wit te exception of pagination. IEEE TRANSACTIONS ON MICROWAVE THEORY AND TECHNIQUES 1 Fast Explicit
More informationContinuous Stochastic Processes
Continuous Stocastic Processes Te term stocastic is often applied to penomena tat vary in time, wile te word random is reserved for penomena tat vary in space. Apart from tis distinction, te modelling
More informationModel Selection in Functional Networks via Genetic Algorithms
Model Selection in Functional Networs via Genetic Algoritms R E Pruneda and B Lacruz Abstract Several statistical tools and most recentl Functional Networs (FN) ave been used to solve nonlinear regression
More information4. The slope of the line 2x 7y = 8 is (a) 2/7 (b) 7/2 (c) 2 (d) 2/7 (e) None of these.
Mat 11. Test Form N Fall 016 Name. Instructions. Te first eleven problems are wort points eac. Te last six problems are wort 5 points eac. For te last six problems, you must use relevant metods of algebra
More informationERROR BOUNDS FOR THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BRADLEY J. LUCIER*
EO BOUNDS FO THE METHODS OF GLIMM, GODUNOV AND LEVEQUE BADLEY J. LUCIE* Abstract. Te expected error in L ) attimet for Glimm s sceme wen applied to a scalar conservation law is bounded by + 2 ) ) /2 T
More informationEDML: A Method for Learning Parameters in Bayesian Networks
: A Metod for Learning Parameters in Bayesian Networks Artur Coi, Kaled S. Refaat and Adnan Darwice Computer Science Department University of California, Los Angeles {aycoi, krefaat, darwice}@cs.ucla.edu
More informationBootstrap prediction intervals for Markov processes
arxiv: arxiv:0000.0000 Bootstrap prediction intervals for Markov processes Li Pan and Dimitris N. Politis Li Pan Department of Matematics University of California San Diego La Jolla, CA 92093-0112, USA
More informationCopyright c 2008 Kevin Long
Lecture 4 Numerical solution of initial value problems Te metods you ve learned so far ave obtained closed-form solutions to initial value problems. A closedform solution is an explicit algebriac formula
More informationAdaptive Neural Filters with Fixed Weights
Adaptive Neural Filters wit Fixed Weigts James T. Lo and Justin Nave Department of Matematics and Statistics University of Maryland Baltimore County Baltimore, MD 150, U.S.A. e-mail: jameslo@umbc.edu Abstract
More information232 Calculus and Structures
3 Calculus and Structures CHAPTER 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS FOR EVALUATING BEAMS Calculus and Structures 33 Copyrigt Capter 17 JUSTIFICATION OF THE AREA AND SLOPE METHODS 17.1 THE
More informationarxiv: v1 [math.na] 17 Jul 2014
Div First-Order System LL* FOSLL* for Second-Order Elliptic Partial Differential Equations Ziqiang Cai Rob Falgout Sun Zang arxiv:1407.4558v1 [mat.na] 17 Jul 2014 February 13, 2018 Abstract. Te first-order
More informationDe-Coupler Design for an Interacting Tanks System
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 2278-1676,p-ISSN: 2320-3331, Volume 7, Issue 3 (Sep. - Oct. 2013), PP 77-81 De-Coupler Design for an Interacting Tanks System
More informationDepartment of Statistics & Operations Research, Aligarh Muslim University, Aligarh, India
Open Journal of Optimization, 04, 3, 68-78 Publised Online December 04 in SciRes. ttp://www.scirp.org/ournal/oop ttp://dx.doi.org/0.436/oop.04.34007 Compromise Allocation for Combined Ratio Estimates of
More informationLECTURE 14 NUMERICAL INTEGRATION. Find
LECTURE 14 NUMERCAL NTEGRATON Find b a fxdx or b a vx ux fx ydy dx Often integration is required. However te form of fx may be suc tat analytical integration would be very difficult or impossible. Use
More informationResearch Article New Results on Multiple Solutions for Nth-Order Fuzzy Differential Equations under Generalized Differentiability
Hindawi Publising Corporation Boundary Value Problems Volume 009, Article ID 395714, 13 pages doi:10.1155/009/395714 Researc Article New Results on Multiple Solutions for Nt-Order Fuzzy Differential Equations
More informationDiscontinuous Galerkin Methods for Relativistic Vlasov-Maxwell System
Discontinuous Galerkin Metods for Relativistic Vlasov-Maxwell System He Yang and Fengyan Li December 1, 16 Abstract e relativistic Vlasov-Maxwell (RVM) system is a kinetic model tat describes te dynamics
More informationQuasiperiodic phenomena in the Van der Pol - Mathieu equation
Quasiperiodic penomena in te Van der Pol - Matieu equation F. Veerman and F. Verulst Department of Matematics, Utrect University P.O. Box 80.010, 3508 TA Utrect Te Neterlands April 8, 009 Abstract Te Van
More informationNumerical Solution of One Dimensional Nonlinear Longitudinal Oscillations in a Class of Generalized Functions
Proc. of te 8t WSEAS Int. Conf. on Matematical Metods and Computational Tecniques in Electrical Engineering, Bucarest, October 16-17, 2006 219 Numerical Solution of One Dimensional Nonlinear Longitudinal
More informationACCURATE SYNTHESIS FORMULAS OBTAINED BY USING A DIFFERENTIAL EVOLUTION ALGORITHM FOR CONDUCTOR-BACKED COPLANAR WAVEG- UIDES
Progress In Electromagnetics Researc M, Vol. 10, 71 81, 2009 ACCURATE SYNTHESIS FORMULAS OBTAINED BY USING A DIFFERENTIAL EVOLUTION ALGORITHM FOR CONDUCTOR-BACKED COPLANAR WAVEG- UIDES S. Kaya, K. Guney,
More informationVolume 29, Issue 3. Existence of competitive equilibrium in economies with multi-member households
Volume 29, Issue 3 Existence of competitive equilibrium in economies wit multi-member ouseolds Noriisa Sato Graduate Scool of Economics, Waseda University Abstract Tis paper focuses on te existence of
More informationQuantum Mechanics Chapter 1.5: An illustration using measurements of particle spin.
I Introduction. Quantum Mecanics Capter.5: An illustration using measurements of particle spin. Quantum mecanics is a teory of pysics tat as been very successful in explaining and predicting many pysical
More informationSOLUTIONS OF FOURTH-ORDER PARTIAL DIFFERENTIAL EQUATIONS IN A NOISE REMOVAL MODEL
Electronic Journal of Differential Equations, Vol. 7(7, No., pp.. ISSN: 7-669. URL: ttp://ejde.mat.txstate.edu or ttp://ejde.mat.unt.edu ftp ejde.mat.txstate.edu (login: ftp SOLUTIONS OF FOURTH-ORDER PARTIAL
More informationDifferential equations. Differential equations
Differential equations A differential equation (DE) describes ow a quantity canges (as a function of time, position, ) d - A ball dropped from a building: t gt () dt d S qx - Uniformly loaded beam: wx
More informationCOMPARISON OF FUZZY LOGIC CONTROLLERS FOR A MULTIVARIABLE PROCESS
COMPARISON OF FUZZY LOGIC CONTROLLERS FOR A MULTIVARIABLE PROCESS KARTHICK S, LAKSHMI P, DEEPA T 3 PG Student, DEEE, College of Engineering, Guindy, Anna University, Cennai Associate Professor, DEEE, College
More informationADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS
ADCP MEASUREMENTS OF VERTICAL FLOW STRUCTURE AND COEFFICIENTS OF FLOAT IN FLOOD FLOWS Yasuo NIHEI (1) and Takeiro SAKAI (2) (1) Department of Civil Engineering, Tokyo University of Science, 2641 Yamazaki,
More informationChemical Engineering & Process Techniques
emical Engineering & Process Tecniques eview Article eedback ontrol for Liquid Level in a Gravity-Drained Multi-Tank System Larry K Jang* Department of emical Engineering, alifornia State University, USA
More informationMath 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006
Mat 102 TEST CHAPTERS 3 & 4 Solutions & Comments Fall 2006 f(x+) f(x) 10 1. For f(x) = x 2 + 2x 5, find ))))))))) and simplify completely. NOTE: **f(x+) is NOT f(x)+! f(x+) f(x) (x+) 2 + 2(x+) 5 ( x 2
More informationBoosting Kernel Density Estimates: a Bias Reduction. Technique?
Boosting Kernel Density Estimates: a Bias Reduction Tecnique? Marco Di Marzio Dipartimento di Metodi Quantitativi e Teoria Economica, Università di Cieti-Pescara, Viale Pindaro 42, 65127 Pescara, Italy
More informationDerivation Of The Schwarzschild Radius Without General Relativity
Derivation Of Te Scwarzscild Radius Witout General Relativity In tis paper I present an alternative metod of deriving te Scwarzscild radius of a black ole. Te metod uses tree of te Planck units formulas:
More informationChemical Engineering & Process Techniques
emical Engineering & Process Tecniques eview Article eedback ontrol for Liquid Level in a Gravity-Drained Multi-Tank System Larry K. Jang* Department of emical Engineering, alifornia State University,
More informationBounds on the Moments for an Ensemble of Random Decision Trees
Noname manuscript No. (will be inserted by te editor) Bounds on te Moments for an Ensemble of Random Decision Trees Amit Durandar Received: Sep. 17, 2013 / Revised: Mar. 04, 2014 / Accepted: Jun. 30, 2014
More informationKernel Density Based Linear Regression Estimate
Kernel Density Based Linear Regression Estimate Weixin Yao and Zibiao Zao Abstract For linear regression models wit non-normally distributed errors, te least squares estimate (LSE will lose some efficiency
More informationTHE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS. L. Trautmann, R. Rabenstein
Worksop on Transforms and Filter Banks (WTFB),Brandenburg, Germany, Marc 999 THE STURM-LIOUVILLE-TRANSFORMATION FOR THE SOLUTION OF VECTOR PARTIAL DIFFERENTIAL EQUATIONS L. Trautmann, R. Rabenstein Lerstul
More informationCS522 - Partial Di erential Equations
CS5 - Partial Di erential Equations Tibor Jánosi April 5, 5 Numerical Di erentiation In principle, di erentiation is a simple operation. Indeed, given a function speci ed as a closed-form formula, its
More information