Stable block Toeplitz matrix for the processing of multichannel seismic data
|
|
- David Maxwell
- 5 years ago
- Views:
Transcription
1 Indian Journal of Marine Sciences Vol. 33(3), Sepember 2004, pp Sable block Toepliz marix for he processing of mulichannel seismic daa Kiri Srivasava* & V P Dimri Naional Geophysical Research Insiue, Uppal Road, Hyderabad , India *[ kiri_ngri@rediffmail.com] Received 29 Sepember 2003; revised 14 June 2004 Compuaion of deconvoluion operaors in he case of single channel secioned inpu/mulichannel seismic daa involves he inversion of a block Toepliz marix. The inversion of such a marix poses several problems. I is well esablished ha he error energy which measures he well posedness of he marix is seen o decrease wih an increase in he filer lengh. However, wih an increase in filer lengh he condiion number of he associaed marix increases. This means ha here is a rade off beween ill posedness and accuracy. The ill-posed problem has been made well posed by a process of (1) normalizaion of he block Toepliz marix and (2) by adding prewhiening parameer. The prewhiening parameer is aken as a few per cen of arihmeic or he geomeric mean of he main diagonal of he block Toepliz marix. Applicaion o a synheic as well as field seismic daa shows ha he condiion number of he associaed block Toepliz marix is reduced by a process of normalizaion and adding prewhiening parameer. Furher i is observed ha he condiion number is smaller when he prewhiening parameer is aken as a few per cenage of geomeric mean as compared o he arihmeic mean. Sabilizing he marix following he above procedure will help in obaining sable as well as accurae deconvoluion operaors. [Key words: Toepliz marix, mulichannel seismic, deconvoluion, filer lengh] Inroducion Mulichannel daa are acquired in many fields e.g. reflecion seismic, seismology, asronomy, sonar and radiomery. Various mulichannel filers have been designed for he suppression of coheren noise, sacking problems and velociy analysis. One algorihm used for he processing of mulichannel daa is known as he LWR (Levinson-Wiener- Robinson) algorihm which is fas and calculaes he required predicion filer for a ime-invarian sequence 1. However, for he ime-varying inpu sequence alernaive mehods o he LWR algorihm such as he gradien mehods have been proposed 2,3. The design of leas squares deconvoluion operaors for mulichannel seismic daa involves he inversion of he block Toepliz marix. The compuaion of he combined deconvoluion operaors for he imevarying sequence being secioned ino piece wise saionary secions also involves he block Toepliz marix 3,5. This block Toepliz marix is srucurally similar o he single channel auocorrelaion marix bu wih a difference ha is elemens are marices raher han scalars and he main propery is ha all he *Corresponding Auhor sub marices are idenical 6. In his invesigaion, a deailed sudy on he general properies of he block Toepliz marix regarding he sabiliy of he soluion and accuracy of he esimaed parameers, hrough he analysis of he condiion number has been carried ou. I has been well esablished ha when he condiion number is nearly one, he sysem is said o be well condiioned i.e. sable. However, unsable sysems are associaed wih large condiion number and are ermed ill condiioned. Hence, i becomes necessary o invesigae he behaviour of condiion number of a marix wih respec o various facors such as normalizaion of he block Toepliz marix and adding prewhiening parameer. The sabiliy of he block Toepliz marix has been discussed in ligh of he prewhiening parameer 7,8. Maerials and Mehods The compuaion of he deconvoluion operaors for a single channel secioned inpu and muli channel daa are designed using normal equaions ha show up he propery of he block Toepliz marix 2. The block Toepliz marix is mahemaically expressed as R P(0) P(1) = P ( 1) P (0) (1)
2 216 INDIAN J. MAR. SCI., VOL. 33, No.3, SEPTEMBER 2004 where P(0) = and P(1) = Rii (0) R (0) R (0) R (0) Rii (1) R (1) R (1) R (1) sonic daa and he covariance marix obained was normalized by: 1/ m ' m 1 τ = τ τ ii τ j= 1 /2 1/2 ( ) R ( ) R ( ) R ( )/( R ( ) R ( τ )) (3) The elemens of he marix are correlaion funcions where i/j indicae he number of channels or secions. Suppose we have wo secions/channels say X 1 () and X 2 (), hen he correlaion funcion is expressed as: R 11 (τ)=e[x 1 (+τ) X 1 ()] auo-correlaion funcion of he X 1 () R 12 (τ)=e X 1 (+τ) X 2 ()] cross-correlaion funcion of X 1 () and X 2 () R 21 (τ)=e[x 2 (+τ) X 1 ()] cross-correlaion funcion of X 2 () and X 1 () R 22 (τ)=e[x 2 (+τ) X 2 ()] auo-correlaion funcion of X 2 () where τ = lag, i.e τ = 0, 1, 2, 3. When τ = 0 and 1 we ge a wo-lengh filer which has a block-toepliz marix represened by Eq. (1) and each elemen of he marix R is a sub marix. In he design of accurae and sable deconvoluion operaors mahemaical operaions need o be performed on he associaed block Toepliz marix, R. The approach aken in his paper is o quanify he sabiliy of he block Toepliz marix hrough he analysis of is condiion number. The condiion number, CN, measures he well posedness of he sysem and is he raio of he maximum eigen value o he minimum eigen value 9. max CN = λ (2) λ min The sysem is said o be well condiioned if he condiion number is nearly one. However, in mos seismic deconvoluion sudies he block Toepliz marix associaed wih is design are associaed wih large condiion number. Hence, he behaviour of he condiion number is invesigaed wih respec o (1) normalizaion of he block Toepliz marix, (2) prewhiening parameer, and (3) lengh of he filer. Normalizaion of he block Toepliz marix In his sudy he esimaed covariance marix i.e. he block Toepliz marix is normalized 10, where he maximum likelihood mehod was applied o several where Rii () τ is he auocorrelaion of he i h channel and R () τ is he cross correlaion beween i h channel 1/ m and j h m channel and R ( τ ) is he geomeric j= 1 mean of he diagonal elemens of he block Toepliz marix. This normalizaion procedure forces he esimaed covariance marix o have equal diagonal elemens, a propery of he Toepliz marix. The prewhiening parameer The prewhiening parameer is applied o avoid he singulariy arising from he zero eigen value. In oher words, we add a small amoun of whie noise o sabilize he covariance marix i.e. prewhiening parameer is added o he main diagonal o make he inverse scheme sable. Prewhiening he Wiener filer improves he oupu signal o ambien noise raio, bu a he same ime i reduces he resoluion 11. The selecion of appropriae prewhiening parameer is a very edious procedure as seen in mos mehods like he ridge regression mehod 12, he singular value decomposiion 13,14 and he maximum likelihood mehod 15. In rouine seismic processing he prewhiening parameer is chosen as a few per cen of he arihmeic mean of he main diagonal. As he geomeric mean is known o resore he overall power of he covariance marix, he prewhiening parameer is aken as a few per cen of he geomeric mean of he main diagonal. Adding prewhiening parameer o he normalized block Toepliz marix resuls in '' ' R ( τ ) = R ( τ) + θi (4) where θ is prewhiening parameer and R (τ) is normalized block Toepliz marix given by Eq. (3) and I is he ideniy marix. Lengh of he filer The sabiliy of he block Toepliz marix associaed wih differen filer lenghs has been sudied by analyzing is condiion number. A well-
3 SRIVASTAVA & DIMRI: PROCESSING OF SEISMIC DATA 217 esablished fac is ha he error energy associaed wih he design of opimum filer decreases wih an increase in he filer lengh. In oher words, we can expec he soluion of he problem o be well posed by increasing he lengh of he filer. However, is accuracy needs o be esed and we shall demonsrae i hrough a numerical example. in boh cases. However, reducion in condiion number is more in he case of geomeric mean as compared o ha of he arihmeic mean. Finally, he effec of he filer lengh i.e. he size of he block Toepliz marix on he condiion number is analyzed. The condiion number for differen filer lenghs has been compued (Table 1), which shows Applicaion o Daa and Resuls Field daa A six-channel unprocessed marine seismic daa colleced from he wesern coninenal margin of India 16 has been aken for he sudy (Fig.1). The marix of auo and cross correlaion associaed wih he design of deconvoluion operaors for he wo inpu races has been obained. The condiion number for he above block Toepliz marix has been compued firs by normalizaion of he marix using Eq. (3). The condiion number has been obained for boh he cases i.e wih normalizaion and wihou normalizaion and i is observed ha he condiion number has reduced by a process of normalizaion (Figs 2 and 3). The effec of prewhiening parameer aken as a few per cen of arihmeic or geomeric mean on he condiion number has also been deal wih. The condiion number has been obained for he normalized block Toepliz marix adding prewhiening parameer (Eq. 4). The prewhiening parameer has been aken as 1 o 20% of he arihmeic or geomeric mean of he main diagonal of he marix and Figs 2 and 3 show he plo of condiion number versus prewhiening parameer for he wo cases. I is furher observed ha wih an increase of prewhiening parameer he condiion number reduces Fig. 2 Condiion number vs prewhiening parameer in he case of arihmeic mean for synheic daa (a) wih normalizaion and (b) wihou normalizaion Fig. 1 Six channels unprocessed marine seismic races Fig. 3 Condiion number vs prewhiening parameer in he case of geomeric mean for synheic daa (a) wih normalizaion and (b) wihou normalizaion.
4 218 INDIAN J. MAR. SCI., VOL. 33, No.3, SEPTEMBER 2004 Table 1 Condiion number for differen filer lenghs in case of marine seismic daa Filer lengh Wih prewhiening Condiion number Wihou prewhiening Fig. 4 Synheic seismic race divided ino wo secions ha he condiion number has increased wih an increase in he filer lengh. Synheic example A ime series can be represened by he convoluion of waveles and noise series i.e. he Woldian decomposiion heorem. Using his concep 17 synheic seismic race is generaed which is mahemaically represened by X ( ) = α A( ) + k β B( ) (5) where α is he source wavele, β is he noise wavele and 1/k is he signal o noise raio of he race respecively. A() is he refleciviy series generaed using random numbers beween [0,1] and B() is he random Gaussian noise series. The choice of he waveles depends on he naure of he race ha one would require. In several seismic exploraions he Ricker wavele is generally used 18. Synheic races have also been generaed using he minimum delay wavele for he source funcion. Hence, in his sudy he synheic seismic races have been obained for a given se of conrolling parameers, i.e., α = 5.,4.,3.,-1.,0.5. β = 4.,-2.,1.,-1.,0.5 k = 1/1.5 Figure 4 shows he plo of he race ha has been obained and has been divided ino wo secions as repored earlier 5. For hese wo secions he compuaion of he filer coefficiens requires he block Toepliz marix 4. The procedure 4 as discussed has been followed and for differen lenghs of he filer he error energy or he filer performance has been compued 3. Afer normalizing he block Toepliz marix and adding some prewhiening he condiion number has been obained. The resuls are ploed in Fig. 5. From he figure i is clear ha as he lengh of he filer increases, he error energy decreases bu he Fig. 5 Plo of (a) filer lengh versus normalized error energy (b) filer lengh versus condiion number for he synheic example. condiion number is increasing. We see ha for a filer lengh of 11 he wo curves are cuing each oher. The lengh of he filer for which he error energy becomes zero has been given earlier 4. For his example, if we were o look for accurae filers we would need o compue a filer of lengh 19. However, since we need o make a compromise he inersecion poin could be assumed o be an opimum lengh ha akes care of he sabiliy and accuracy of he filer design. Discussion In his sudy we have considered a six-channel marine seismic daa and synheic seismic races and for achieving he sabiliy of he deconvoluion operaors. Single channel secioned inpu and mulichannel seismic deconvoluion involve he inversion of he block Toepliz marix. To obain a sable soluion, he block Toepliz marix associaed wih he deconvoluion operaors has been normalized
5 SRIVASTAVA & DIMRI: PROCESSING OF SEISMIC DATA 219 o ge low values of he condiion number as his measures he well posedness of he problem. Furher, he sabiliy of he algorihm has been achieved by adding he prewhiening parameer. I is well known ha he prewhiening parameer prevens he inverse from blowing. However, i is observed ha wih an increase in he prewhiening parameer he condiion number decreases and also i has been found ha he prewhiening parameer compued from he geomeric mean as compared o he arihmeic mean of he elemens of he main diagonal of he block Toepliz marix reduces he condiion number furher. We have also seen ha wih an increase in he filer lengh he condiion number also increases. However, i is well known ha wih an increase in filer lengh he error energy decreases hereby improving he performance of he filer 2. A rade off beween well posedness and accuracy of he problem exiss and one needs o be careful in choosing he lengh of he filer in such a way ha sabiliy as well as accuracy is mainained. This seismic daa has been used o remove he muliples using he adapive deconvoluion 19. The combined convergence facor for he adapive algorihm is obained. The auo and cross correlaion marix has been sabilized using he same procedure o obain he maximum eigen value, ha is used o obain he convergen facor of he adapive algorihm. The procedure can be applied in various geophysical sudies bu care should be aken for appropriae compromise beween he sabiliy and accuracy of he esimaed parameer. This sudy will be very useful in designing deconvoluion operaors for a single channel secioned inpu/muli channel daa ha will accoun for is sabiliy, accuracy and well posedness. References 1 Robinson E A & Durrani T S, Geophysical signal processing (Prenice-Hall Inc, Englewood Cliffs, NJ), 1986, pp Dimri V P, Deconvoluion and inverse heory (Elsevier Science Publishers, Amserdam), 1992, pp Dimri V P, On he ime-varying Wiener filer, Geoph Pros, 34 (1986) Dimri V P & Srivasava K, Ideal performance crieria for deconvoluion operaors, Geoph Pros, 35 (1987) Dimri V P & Srivasava K, The opimum gae lengh for he ime-varying deconvoluion operaor, Geoph Pros, 38 (1991) Treiel S, Principles of digial mulichannel filering, Geophysics, 35 (1970) Wang R J & Treiel S, The deerminaion of digial Wiener filers by means of gradien mehods, Geophysics, 38 (1973) Treiel S & Lines L R, Linear inverse heory and deconvoluion, Geophysics, 47 (1982) Treiel S & Wang R J, The deerminaion of digial Wiener filers from an ill-condiioned sysem of normal equaions, Geoph Pros, 24 (1976) Hsu K & Baaggeroer A B, Applicaion of he maximumlikelihood mehod (MLM) for sonic velociy logging, Geophysics, 51 (1986) Bickel S H & Marinez D R, Resoluion performance of Wiener filers, Geophysics, 48 (1983) Inman, J R, Resisiviy inversion wih ridge regression, Geophysics, 40 (1975) Lines, L R & Treiel S, A review of leas-squares inversion and is applicaion o geophysical problems, Geoph Pros, 32 (1984) Ursin, B & Zheng Y, Idenificaion of seismic reflecions using singular value decomposiion, Geoph Pros, 33 (1985) Ursin, B, Mehods for esimaing he seismic reflecion response, Geophysics, 62 (1997) Reddy, S I, Roy Chowdhary K, Drollia R K, Ashalaha B, Mial G S, Subrahmanyam C S & Singh R N, On he srucure of he wesern coninenal margin off Mangalore coas, India, J Ass Expl Geoph, 9 (1988) Dash Bibhu P & Obaidullah K A, Deerminaion of signal and noise saisics using correlaion heory, Geophysics, 35 (1970) Ricker N, Wavele conracion, wavele expansion, and he conrol of seismic resoluion, Geophysics, 18 (1953) Srivasava K, An adapive scheme for processing mulirace marine seismic daa, Indian J Mar Sci, 29 (2000)
Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationGeorey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract
Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical
More informationRobust estimation based on the first- and third-moment restrictions of the power transformation model
h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,
More informationSTATE-SPACE MODELLING. A mass balance across the tank gives:
B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing
More informationt is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...
Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger
More informationReferences are appeared in the last slide. Last update: (1393/08/19)
SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be
More informationZürich. ETH Master Course: L Autonomous Mobile Robots Localization II
Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),
More informationCHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK
175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationRecursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems
8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear
More informationThe General Linear Test in the Ridge Regression
ommunicaions for Saisical Applicaions Mehods 2014, Vol. 21, No. 4, 297 307 DOI: hp://dx.doi.org/10.5351/sam.2014.21.4.297 Prin ISSN 2287-7843 / Online ISSN 2383-4757 The General Linear Tes in he Ridge
More informationModal identification of structures from roving input data by means of maximum likelihood estimation of the state space model
Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix
More informationSliding Mode Controller for Unstable Systems
S. SIVARAMAKRISHNAN e al., Sliding Mode Conroller for Unsable Sysems, Chem. Biochem. Eng. Q. 22 (1) 41 47 (28) 41 Sliding Mode Conroller for Unsable Sysems S. Sivaramakrishnan, A. K. Tangirala, and M.
More informationNature Neuroscience: doi: /nn Supplementary Figure 1. Spike-count autocorrelations in time.
Supplemenary Figure 1 Spike-coun auocorrelaions in ime. Normalized auocorrelaion marices are shown for each area in a daase. The marix shows he mean correlaion of he spike coun in each ime bin wih he spike
More informationStability and Bifurcation in a Neural Network Model with Two Delays
Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy
More informationPhysics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationNavneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi
Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec
More informationR t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t
Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,
More informationKriging Models Predicting Atrazine Concentrations in Surface Water Draining Agricultural Watersheds
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 Kriging Models Predicing Arazine Concenraions in Surface Waer Draining Agriculural Waersheds Paul L. Mosquin, Jeremy Aldworh, Wenlin Chen Supplemenal Maerial Number
More informationState-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter
Sae-Space Models Iniializaion, Esimaion and Smoohing of he Kalman Filer Iniializaion of he Kalman Filer The Kalman filer shows how o updae pas predicors and he corresponding predicion error variances when
More informationSub Module 2.6. Measurement of transient temperature
Mechanical Measuremens Prof. S.P.Venkaeshan Sub Module 2.6 Measuremen of ransien emperaure Many processes of engineering relevance involve variaions wih respec o ime. The sysem properies like emperaure,
More information10. State Space Methods
. Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he
More informationSUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL
HE PUBLISHING HOUSE PROCEEDINGS OF HE ROMANIAN ACADEMY, Series A, OF HE ROMANIAN ACADEMY Volume, Number 4/200, pp 287 293 SUFFICIEN CONDIIONS FOR EXISENCE SOLUION OF LINEAR WO-POIN BOUNDARY PROBLEM IN
More informationAn introduction to the theory of SDDP algorithm
An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking
More informationExponential Smoothing
Exponenial moohing Inroducion A simple mehod for forecasing. Does no require long series. Enables o decompose he series ino a rend and seasonal effecs. Paricularly useful mehod when here is a need o forecas
More informationAn recursive analytical technique to estimate time dependent physical parameters in the presence of noise processes
WHAT IS A KALMAN FILTER An recursive analyical echnique o esimae ime dependen physical parameers in he presence of noise processes Example of a ime and frequency applicaion: Offse beween wo clocks PREDICTORS,
More informationNumerical Dispersion
eview of Linear Numerical Sabiliy Numerical Dispersion n he previous lecure, we considered he linear numerical sabiliy of boh advecion and diffusion erms when approimaed wih several spaial and emporal
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More informationEcon107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)
I. Definiions and Problems A. Perfec Mulicollineariy Econ7 Applied Economerics Topic 7: Mulicollineariy (Sudenmund, Chaper 8) Definiion: Perfec mulicollineariy exiss in a following K-variable regression
More informationFor example, the comb filter generated from. ( ) has a transfer function. e ) has L notches at ω = (2k+1)π/L and L peaks at ω = 2π k/l,
Comb Filers The simple filers discussed so far are characeried eiher by a single passband and/or a single sopband There are applicaions where filers wih muliple passbands and sopbands are required The
More informationψ(t) = V x (0)V x (t)
.93 Home Work Se No. (Professor Sow-Hsin Chen Spring Term 5. Due March 7, 5. This problem concerns calculaions of analyical expressions for he self-inermediae scaering funcion (ISF of he es paricle in
More informationLecture 10 - Model Identification
Lecure - odel Idenificaion Wha is ssem idenificaion? Direc impulse response idenificaion Linear regression Regularizaion Parameric model ID nonlinear LS Conrol Engineering - Wha is Ssem Idenificaion? Experimen
More informationDEPARTMENT OF STATISTICS
A Tes for Mulivariae ARCH Effecs R. Sco Hacker and Abdulnasser Haemi-J 004: DEPARTMENT OF STATISTICS S-0 07 LUND SWEDEN A Tes for Mulivariae ARCH Effecs R. Sco Hacker Jönköping Inernaional Business School
More information14 Autoregressive Moving Average Models
14 Auoregressive Moving Average Models In his chaper an imporan parameric family of saionary ime series is inroduced, he family of he auoregressive moving average, or ARMA, processes. For a large class
More informationF This leads to an unstable mode which is not observable at the output thus cannot be controlled by feeding back.
Lecure 8 Las ime: Semi-free configuraion design This is equivalen o: Noe ns, ener he sysem a he same place. is fixed. We design C (and perhaps B. We mus sabilize if i is given as unsable. Cs ( H( s = +
More information4.1 Other Interpretations of Ridge Regression
CHAPTER 4 FURTHER RIDGE THEORY 4. Oher Inerpreaions of Ridge Regression In his secion we will presen hree inerpreaions for he use of ridge regression. The firs one is analogous o Hoerl and Kennard reasoning
More informationCombined Bending with Induced or Applied Torsion of FRP I-Section Beams
Combined Bending wih Induced or Applied Torsion of FRP I-Secion Beams MOJTABA B. SIRJANI School of Science and Technology Norfolk Sae Universiy Norfolk, Virginia 34504 USA sirjani@nsu.edu STEA B. BONDI
More information6.003 Homework #9 Solutions
6.00 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 0 a 0 5 a k sin πk 5 sin πk 5 πk for k 0 a k 0 πk j
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationExcel-Based Solution Method For The Optimal Policy Of The Hadley And Whittin s Exact Model With Arma Demand
Excel-Based Soluion Mehod For The Opimal Policy Of The Hadley And Whiin s Exac Model Wih Arma Demand Kal Nami School of Business and Economics Winson Salem Sae Universiy Winson Salem, NC 27110 Phone: (336)750-2338
More informationVectorautoregressive Model and Cointegration Analysis. Time Series Analysis Dr. Sevtap Kestel 1
Vecorauoregressive Model and Coinegraion Analysis Par V Time Series Analysis Dr. Sevap Kesel 1 Vecorauoregression Vecor auoregression (VAR) is an economeric model used o capure he evoluion and he inerdependencies
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More information3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon
3..3 INRODUCION O DYNAMIC OPIMIZAION: DISCREE IME PROBLEMS A. he Hamilonian and Firs-Order Condiions in a Finie ime Horizon Define a new funcion, he Hamilonian funcion, H. H he change in he oal value of
More informationdi Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.
di ernardo, M. (995). A purely adapive conroller o synchronize and conrol chaoic sysems. hps://doi.org/.6/375-96(96)8-x Early version, also known as pre-prin Link o published version (if available):.6/375-96(96)8-x
More informationOn Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature
On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check
More informationBias in Conditional and Unconditional Fixed Effects Logit Estimation: a Correction * Tom Coupé
Bias in Condiional and Uncondiional Fixed Effecs Logi Esimaion: a Correcion * Tom Coupé Economics Educaion and Research Consorium, Naional Universiy of Kyiv Mohyla Academy Address: Vul Voloska 10, 04070
More information2 int T. is the Fourier transform of f(t) which is the inverse Fourier transform of f. i t e
PHYS67 Class 3 ourier Transforms In he limi T, he ourier series becomes an inegral ( nt f in T ce f n f f e d, has been replaced by ) where i f e d is he ourier ransform of f() which is he inverse ourier
More information2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006
2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)
More informationON THE BEAT PHENOMENON IN COUPLED SYSTEMS
8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More information0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED
0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable
More informationL07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms
L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)
More informationAnalyze patterns and relationships. 3. Generate two numerical patterns using AC
envision ah 2.0 5h Grade ah Curriculum Quarer 1 Quarer 2 Quarer 3 Quarer 4 andards: =ajor =upporing =Addiional Firs 30 Day 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 andards: Operaions and Algebraic Thinking
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationKEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow
KEY Mah 334 Miderm III Winer 008 secion 00 Insrucor: Sco Glasgow Please do NOT wrie on his exam. No credi will be given for such work. Raher wrie in a blue book, or on your own paper, preferably engineering
More informationMATH 128A, SUMMER 2009, FINAL EXAM SOLUTION
MATH 28A, SUMME 2009, FINAL EXAM SOLUTION BENJAMIN JOHNSON () (8 poins) [Lagrange Inerpolaion] (a) (4 poins) Le f be a funcion defined a some real numbers x 0,..., x n. Give a defining equaion for he Lagrange
More informationDiebold, Chapter 7. Francis X. Diebold, Elements of Forecasting, 4th Edition (Mason, Ohio: Cengage Learning, 2006). Chapter 7. Characterizing Cycles
Diebold, Chaper 7 Francis X. Diebold, Elemens of Forecasing, 4h Ediion (Mason, Ohio: Cengage Learning, 006). Chaper 7. Characerizing Cycles Afer compleing his reading you should be able o: Define covariance
More informationHow to Deal with Structural Breaks in Practical Cointegration Analysis
How o Deal wih Srucural Breaks in Pracical Coinegraion Analysis Roselyne Joyeux * School of Economic and Financial Sudies Macquarie Universiy December 00 ABSTRACT In his noe we consider he reamen of srucural
More information6.003 Homework #9 Solutions
6.003 Homework #9 Soluions Problems. Fourier varieies a. Deermine he Fourier series coefficiens of he following signal, which is periodic in 0. x () 0 3 0 a 0 5 a k a k 0 πk j3 e 0 e j πk 0 jπk πk e 0
More informationOBJECTIVES OF TIME SERIES ANALYSIS
OBJECTIVES OF TIME SERIES ANALYSIS Undersanding he dynamic or imedependen srucure of he observaions of a single series (univariae analysis) Forecasing of fuure observaions Asceraining he leading, lagging
More informationTime series Decomposition method
Time series Decomposiion mehod A ime series is described using a mulifacor model such as = f (rend, cyclical, seasonal, error) = f (T, C, S, e) Long- Iner-mediaed Seasonal Irregular erm erm effec, effec,
More informationd 1 = c 1 b 2 - b 1 c 2 d 2 = c 1 b 3 - b 1 c 3
and d = c b - b c c d = c b - b c c This process is coninued unil he nh row has been compleed. The complee array of coefficiens is riangular. Noe ha in developing he array an enire row may be divided or
More informationUnderstanding the asymptotic behaviour of empirical Bayes methods
Undersanding he asympoic behaviour of empirical Bayes mehods Boond Szabo, Aad van der Vaar and Harry van Zanen EURANDOM, 11.10.2011. Conens 2/20 Moivaion Nonparameric Bayesian saisics Signal in Whie noise
More informationThe Potential Effectiveness of the Detection of Pulsed Signals in the Non-Uniform Sampling
The Poenial Effeciveness of he Deecion of Pulsed Signals in he Non-Uniform Sampling Arhur Smirnov, Sanislav Vorobiev and Ajih Abraham 3, 4 Deparmen of Compuer Science, Universiy of Illinois a Chicago,
More informationMean-square Stability Control for Networked Systems with Stochastic Time Delay
JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455
More informationOptimal Paired Choice Block Designs. Supplementary Material
Saisica Sinica: Supplemen Opimal Paired Choice Block Designs Rakhi Singh 1, Ashish Das 2 and Feng-Shun Chai 3 1 IITB-Monash Research Academy, Mumbai, India 2 Indian Insiue of Technology Bombay, Mumbai,
More informationTHE DISCRETE WAVELET TRANSFORM
. 4 THE DISCRETE WAVELET TRANSFORM 4 1 Chaper 4: THE DISCRETE WAVELET TRANSFORM 4 2 4.1 INTRODUCTION TO DISCRETE WAVELET THEORY The bes way o inroduce waveles is hrough heir comparison o Fourier ransforms,
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More information5 The fitting methods used in the normalization of DSD
The fiing mehods used in he normalizaion of DSD.1 Inroducion Sempere-Torres e al. 1994 presened a general formulaion for he DSD ha was able o reproduce and inerpre all previous sudies of DSD. The mehodology
More informationThe equation to any straight line can be expressed in the form:
Sring Graphs Par 1 Answers 1 TI-Nspire Invesigaion Suden min Aims Deermine a series of equaions of sraigh lines o form a paern similar o ha formed by he cables on he Jerusalem Chords Bridge. Deermine he
More informationLecture 20: Riccati Equations and Least Squares Feedback Control
34-5 LINEAR SYSTEMS Lecure : Riccai Equaions and Leas Squares Feedback Conrol 5.6.4 Sae Feedback via Riccai Equaions A recursive approach in generaing he marix-valued funcion W ( ) equaion for i for he
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More informationModule 4: Time Response of discrete time systems Lecture Note 2
Module 4: Time Response of discree ime sysems Lecure Noe 2 1 Prooype second order sysem The sudy of a second order sysem is imporan because many higher order sysem can be approimaed by a second order model
More informationFlow-Induced Vibration Analysis of Supported Pipes with a Crack
Flow-Induced Vibraion Analsis of Suppored Pipes wih a Crack Jin-Huk Lee, Samer Masoud Al-Said Deparmen of Mechanical Engineering American Universi of Sharjah, UAE Ouline Inroducion and Moivaion Aeroacousicall
More informationApplication of the radial basis function neural network to the prediction of log properties from seismic attributes
Applicaion of he RBF Applicaion of he radial basis funcion neural newor o he predicion of log properies from seismic aribues Brian H. Russell, Laurence R. Lines, and Daniel P. Hampson ABSTRACT In his paper,
More informationWavelet Variance, Covariance and Correlation Analysis of BSE and NSE Indexes Financial Time Series
Wavele Variance, Covariance and Correlaion Analysis of BSE and NSE Indexes Financial Time Series Anu Kumar 1*, Sangeea Pan 1, Lokesh Kumar Joshi 1 Deparmen of Mahemaics, Universiy of Peroleum & Energy
More informationLicenciatura de ADE y Licenciatura conjunta Derecho y ADE. Hoja de ejercicios 2 PARTE A
Licenciaura de ADE y Licenciaura conjuna Derecho y ADE Hoja de ejercicios PARTE A 1. Consider he following models Δy = 0.8 + ε (1 + 0.8L) Δ 1 y = ε where ε and ε are independen whie noise processes. In
More informationNCSS Statistical Software. , contains a periodic (cyclic) component. A natural model of the periodic component would be
NCSS Saisical Sofware Chaper 468 Specral Analysis Inroducion This program calculaes and displays he periodogram and specrum of a ime series. This is someimes nown as harmonic analysis or he frequency approach
More information1 Differential Equation Investigations using Customizable
Differenial Equaion Invesigaions using Cusomizable Mahles Rober Decker The Universiy of Harford Absrac. The auhor has developed some plaform independen, freely available, ineracive programs (mahles) for
More informationTime series model fitting via Kalman smoothing and EM estimation in TimeModels.jl
Time series model fiing via Kalman smoohing and EM esimaion in TimeModels.jl Gord Sephen Las updaed: January 206 Conens Inroducion 2. Moivaion and Acknowledgemens....................... 2.2 Noaion......................................
More informationBlock Diagram of a DCS in 411
Informaion source Forma A/D From oher sources Pulse modu. Muliplex Bandpass modu. X M h: channel impulse response m i g i s i Digial inpu Digial oupu iming and synchronizaion Digial baseband/ bandpass
More informationThe Rosenblatt s LMS algorithm for Perceptron (1958) is built around a linear neuron (a neuron with a linear
In The name of God Lecure4: Percepron and AALIE r. Majid MjidGhoshunih Inroducion The Rosenbla s LMS algorihm for Percepron 958 is buil around a linear neuron a neuron ih a linear acivaion funcion. Hoever,
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationRandom Walk with Anti-Correlated Steps
Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and
More informationTechnical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.
Technical Repor Doc ID: TR--203 06-March-203 (Las revision: 23-Februar-206) On formulaing quadraic funcions in opimizaion models. Auhor: Erling D. Andersen Convex quadraic consrains quie frequenl appear
More informationINTRODUCTION TO MACHINE LEARNING 3RD EDITION
ETHEM ALPAYDIN The MIT Press, 2014 Lecure Slides for INTRODUCTION TO MACHINE LEARNING 3RD EDITION alpaydin@boun.edu.r hp://www.cmpe.boun.edu.r/~ehem/i2ml3e CHAPTER 2: SUPERVISED LEARNING Learning a Class
More informationAnalysis of Microstrip Coupling Gap to Estimate Polymer Permittivity
Analysis of Microsrip Couplin Gap o Esimae Polymer Permiiviy Chanchal Yadav Deparmen of Physics & Elecronics Rajdhani Collee, Universiy of Delhi Delhi, India Absrac A ap in he microsrip line can be modeled
More information( ) ( ) if t = t. It must satisfy the identity. So, bulkiness of the unit impulse (hyper)function is equal to 1. The defining characteristic is
UNIT IMPULSE RESPONSE, UNIT STEP RESPONSE, STABILITY. Uni impulse funcion (Dirac dela funcion, dela funcion) rigorously defined is no sricly a funcion, bu disribuion (or measure), precise reamen requires
More informationLecture 2 October ε-approximation of 2-player zero-sum games
Opimizaion II Winer 009/10 Lecurer: Khaled Elbassioni Lecure Ocober 19 1 ε-approximaion of -player zero-sum games In his lecure we give a randomized ficiious play algorihm for obaining an approximae soluion
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationSTRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN
Inernaional Journal of Applied Economerics and Quaniaive Sudies. Vol.1-3(004) STRUCTURAL CHANGE IN TIME SERIES OF THE EXCHANGE RATES BETWEEN YEN-DOLLAR AND YEN-EURO IN 001-004 OBARA, Takashi * Absrac The
More informationVariational Iteration Method for Solving System of Fractional Order Ordinary Differential Equations
IOSR Journal of Mahemaics (IOSR-JM) e-issn: 2278-5728, p-issn: 2319-765X. Volume 1, Issue 6 Ver. II (Nov - Dec. 214), PP 48-54 Variaional Ieraion Mehod for Solving Sysem of Fracional Order Ordinary Differenial
More informationArticle from. Predictive Analytics and Futurism. July 2016 Issue 13
Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning
More informationTesting for a Single Factor Model in the Multivariate State Space Framework
esing for a Single Facor Model in he Mulivariae Sae Space Framework Chen C.-Y. M. Chiba and M. Kobayashi Inernaional Graduae School of Social Sciences Yokohama Naional Universiy Japan Faculy of Economics
More informationv A Since the axial rigidity k ij is defined by P/v A, we obtain Pa 3
The The rd rd Inernaional Conference on on Design Engineering and Science, ICDES 14 Pilsen, Czech Pilsen, Republic, Czech Augus Republic, 1 Sepember 1-, 14 In-plane and Ou-of-plane Deflecion of J-shaped
More informationA Bayesian Approach to Spectral Analysis
Chirped Signals A Bayesian Approach o Specral Analysis Chirped signals are oscillaing signals wih ime variable frequencies, usually wih a linear variaion of frequency wih ime. E.g. f() = A cos(ω + α 2
More information23.2. Representing Periodic Functions by Fourier Series. Introduction. Prerequisites. Learning Outcomes
Represening Periodic Funcions by Fourier Series 3. Inroducion In his Secion we show how a periodic funcion can be expressed as a series of sines and cosines. We begin by obaining some sandard inegrals
More informationChoice of Spectral Density Estimator in Ng-Perron Test: A Comparative Analysis
Inernaional Economeric Review (IER) Choice of Specral Densiy Esimaor in Ng-Perron Tes: A Comparaive Analysis Muhammad Irfan Malik and Aiq-ur-Rehman Inernaional Islamic Universiy Islamabad and Inernaional
More informationForecasting optimally
I) ile: Forecas Evaluaion II) Conens: Evaluaing forecass, properies of opimal forecass, esing properies of opimal forecass, saisical comparison of forecas accuracy III) Documenaion: - Diebold, Francis
More information