( x) SUMMARY (2) α z. (3) Now approximating z derivative by forward difference
|
|
- Joy Morris
- 6 years ago
- Views:
Transcription
1 Depth etrapolatio of seismic wavefields usig cubic splie approimatio Suhas Phadke* ad Dheera Bhardwa Cetre for Developmet of Advaced Computig, Pue Uiversity Campus, GaeshKhid, Pue 4007, Idia SUMMARY Depth etrapolatio equatio used for seismic migratio is ofte solved by fiite-differece techique. The most commoly used migratio method is based upo the Crak- Nicolso implicit scheme. Some modificatios to this scheme are i practice which improve the impulse respose of the migratio operator. I this paper we propose a differet approach where the wavefield is approimated by cubic splie fuctio. The approimatio to secod derivative resultig from this approach is similar to that of the /6th trick of Claerbout. There is aother splie parameter θ which cotrols the ature of impulse respose. The scheme is ucoditioally stable for 05. θ 0.. As the value of θ icreases the dispersed evaescet eergy gets more ad more atteuated. The method is demostrated by calculatig the impulse respose ad by applyig it to a sythetic data set. INTRODUCTION Seismic migratio is a key step i imagig of the udergroud geological structures. For laterally varyig velocity structures, the methods based upo the parabolic approimatio (Claerbout 985) of the wave equatio are most commo i use. Several migratio methods have bee proposed which use fiite-differece approimatio to the oe way wave equatio for dowward etrapolatio of the wavefield. Both implicit (Claerbout, 985) ad eplicit (Hale, 99) methods are i practice. Eplicit methods eed special care i their implemetatio because of the stability criteria, whereas implicit methods are ucoditioally stable. I this paper we have developed a implicit scheme for the depth etrapolatio of seismic wavefield usig cubic splie approimatio. The first part of the paper gives a mathematical descriptio of the method. Net we demostrate the usefuless of this method by showig the impulse resposes. The method is fially tested by applyig it to a sythetic data set, followed by discussio ad coclusios. MATHEMATICAL FORMULATION Wave etrapolatio equatio, which is accurate for migratig dips upto 45 degree is derived from the dispersio relatio kv α kv z ω ω kv β () ω where ω is the agular frequecy, v is the wave velocity ad k ad k z are horizotal ad vertical waveumbers respectively. For 45 degree accuracy α ad β are ad respectively. Accuracy for larger dips ca be obtaied by choosig α ad β i some optimal way (Yilmaz, 987). I terms of the retarded wavefield (Claerbout 985) the differetial equatio resultig from the diffractio term is give by im Q Q + Q 0 m z () α z α where m (ω/v) ad i is the square root of egative oe. The other term called thi les term is solved aalytically. Equatio () is usually solved usig some fiite differece approimatio. Here we preset a differet approach for solvig equatio (). Rewrite equatio () as β z i Q αm im Q + α Q 0 (3) Now approimatig z derivative by forward differece approimatio ad ( Q ) by secod order derivatives ( M ) of cubic splie fuctio S ( ) (where is mesh iterval) iterpolatig Q ( 0,,,... J), we get αm + + ( M M ) im ( Q Q ) z α z + + θm + ( θ) M 0 (4) where, θ [0,] is a cubic splie parameter ad z is icremet i z-directio. Now the cotiuity coditio for secod order derivative of cubic splie fuctio S ( ) (Ahlberg et. al., 967) gives 6 ( ) ( Q Q Q ) + + M + 4M + M + (5) or ( ) 6 ( ) ( Q Q Q ) + + M M + M + 6M + δ M 6 + 6M δ Q ( ) δ M Q δ + 6 where δ Q Q Q + Q (7) +. (6)
2 Depth Etrapolatio usig Cubic Splies Approimatio of the secod derivative by the epressio give i (6) is the same as that of the /6 trick of Claerbout (985). He calls it a less obvious epressio that offers more accuracy at less cost, without givig a eplaatio for the derivatio of the epressio used. Here we have derived it by approimatig the wavefield by cubic splie fuctio approimatio. This approimatio gives better results because cubic splies provide better approimatio to a give fuctio (Ahlberg et. al., 967). Now substitutig M from equatio (6) i equatio (4), ad multiplyig by z + δ 6, we obtai ( )( ) θ( z) δ + im δ + + δ Q + Q + Q α 6 im δ θ( z) δ Q + Q δ Q α 6 (8) 0. Distace (km) 0.. θ 0. Usig (7) i (8) we obtai after some algebraic maipulatios Distace (km) AQ+ + BQ + AQ AQ + BQ + AQ (9) θ Figure : Impulse respose of migratio algorithm for differet values of θ. N as the value of θ icreases. Steep dips are also atteuated with the icrease i where 06 im θ( z) A + 6α im θ ( z) B 3 α im ( θ)( z) A 6α im ( θ)( z) B + 3 α
3 Depth Etrapolatio usig Cubic Splies i Figure for differet values of θ. For θ oe ca observe a large amout of dispersed evaescet eergy. Also at steep dips there is dispersio of low ad high frequecies. For θ 055. the dispersed evaescet eergy has bee atteuated to a great etet. Also all the dips are correctly positioed. For θ 06. ad θ 07. the dispersed evaescet eergy is further reduced. Steep dips are also atteuated, but the remaiig dips are correctly positioed alog cocetric semicircles. Oe has to decide o a optimal choice for θ, which reduces the udesired evaescet eergy as well as correctly positios all the required dips. θ 06. seems to be the best choice from the practical poit of view. SYNTHETIC DATA EXAMPLE Net we applied the migratio method based upo the cubic splie approimatio to a sythetic data set. The velocity model used for geeratig the sythetic data is show i Figure. The model comprises of a weathered layer o top of a dippig layer ad a graba like structure. Sythetic seismograms were calculated for this model for 46 source positios with a source iterval of 50m ad a receiver iterval of 5m. A higher order fiite differece modellig algorithm based upo acoustic wave equatio was used with a grid spacig of 5m. A Ricker wavelet with a domiat frequecy of 30Hz was used as the source fuctio. Distace (km) m/s 500 m/s 850 m/s 3400 m/s Equatio (9) is a tridiagoal system of equatio. The solutio of (9) gives the wavefield at depth (+) i terms of the wavefield at depth (). This implicit scheme is ucoditioally stable for 05. θ 0.. Absorbig boudary coditios are used o the sides of the model. IMPULSE RESPONSE I order to test the etrapolatio method described above, a 45 degree migratio program was developed. The iput to the migratio program was a sectio cotaiig three bad limited Ricker wavelets i the cetre of the sectio. The domiat frequecy of the wavelets was 30 Hz. Spatial samplig of 8m ad a samplig rate of ms was used. The migratio of this data set yields the impulse resposes show Figure : Velocity model used for the geeratio of sythetic data set. The acoustic wave modellig code was implemeted i a distributed computig eviromet usig PVM (Parallel Virtual Machie) message passig calls (Geist et. al. 994). For both, modellig ad migratio eamples preseted i this paper we have used 8 UltraSparc Workstatios etworked usig a fast etheret switch. This is a part of the facility called PARAM OpeFrame. Each UltraSparc Workstatio has a 00 MHz CPU with 8 MB of RAM ad Solaris operatig system.
4 Depth Etrapolatio usig Cubic Splies Distace (km) Distace (km).5.0 Two-way time (sec) Figure 3: A CDP stacked sectio of the sythetic data set. A CDP stacked sectio of the sythetic data set is show i Figure 3. Net this stacked data set was migrated usig the above algorithm. The migrated output sectios for two differet values of θ are show i Figure 4. O both the migrated sectios all the evets are properly imaged. Oe ca otice that for θ 05. there is more evaescet eergy o the migrated sectio as compared to the migrated sectio with θ 06.. A user ca decide o the appropriate value for θ. DISCUSSION AND CONCLUSIONS I this paper we have developed ad demostrated a implicit scheme for the depth etrapolatio of the seismic wavefield usig cubic splie approimatio. The approimatio to the secod derivative is derived from the cotiuity coditio for the secod order derivative of cubic splie fuctio. The differece approimatio to the oe-way wave equatio resultig from this approach has aother splie parameter θ, which cotrols the dispersed evaescet eergy. As the value of θ icreases more ad more evaescet eergy gets atteuated. Steep dips also get atteuated with icrease i the value of θ. A proper choice of this parameter helps i reducig the udesired oise o the migrated sectios. This approach ca be easily icorporated i D prestack migratio algorithm as well as i 3D prestack ad poststack migratio algorithms. ACKNOWLEDGEMENTS Authors wish to epress their gratitude to the Departmet of Sciece ad Techology (DST), Govermet of Idia, for fudig the seismic data processig proect. Authors also wish to thak the Cetre for Developmet of Advaced Computig (C-DAC), Pue for providig the computatioal facilities o PARAM OpeFrame ad permissio to publish this work..5 θ Distace (km).5.0 θ Figure 4: The migrated seismic sectios of the stacked data show i Figure 3 for θ 05. ad θ 06.. Notice the reductio i evaescet eergy for θ 06.. REFERENCES Ahlberg, J. H., Nilso, E. N., ad Walsh, J. L., 967, The theory of splies ad their applicatios, Academic Press. Claerbout, J. F., 985, Imagig the Earth s iterior, Blackwell Scietific Publicatios. Geist, A., Begueli, A., Dogarra, J., Jiag, W., Machek, R., ad Suderam, V., 994, PVM: Parallel Virtual Machie, A users guide ad tutorial for etworked parallel computig, MIT Press, Cambridge, Massachusetts. Hale, D., 99, Stable eplicit depth etrapolatio of seismic wavefields, Geophysics, 56, p
5 Yilmez, O., 987, Seismic Data Processig, Society of Eploratio Geophysicists. Depth Etrapolatio usig Cubic Splies
DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS. Park Road, Islamabad, Pakistan
Mathematical ad Computatioal Applicatios, Vol. 9, No. 3, pp. 30-40, 04 DECOMPOSITION METHOD FOR SOLVING A SYSTEM OF THIRD-ORDER BOUNDARY VALUE PROBLEMS Muhammad Aslam Noor, Khalida Iayat Noor ad Asif Waheed
More informationx x x 2x x N ( ) p NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS By Newton-Raphson formula
NUMERICAL METHODS UNIT-I-SOLUTION OF EQUATIONS AND EIGENVALUE PROBLEMS. If g( is cotiuous i [a,b], te uder wat coditio te iterative (or iteratio metod = g( as a uique solutio i [a,b]? '( i [a,b].. Wat
More informationA numerical Technique Finite Volume Method for Solving Diffusion 2D Problem
The Iteratioal Joural Of Egieerig d Sciece (IJES) Volume 4 Issue 10 Pages PP -35-41 2015 ISSN (e): 2319 1813 ISSN (p): 2319 1805 umerical Techique Fiite Volume Method for Solvig Diffusio 2D Problem 1 Mohammed
More informationComparison Study of Series Approximation. and Convergence between Chebyshev. and Legendre Series
Applied Mathematical Scieces, Vol. 7, 03, o. 6, 3-337 HIKARI Ltd, www.m-hikari.com http://d.doi.org/0.988/ams.03.3430 Compariso Study of Series Approimatio ad Covergece betwee Chebyshev ad Legedre Series
More informationNumerical Solutions of Second Order Boundary Value Problems by Galerkin Residual Method on Using Legendre Polynomials
IOSR Joural of Mathematics (IOSR-JM) e-issn: 78-578, p-issn: 319-765X. Volume 11, Issue 6 Ver. IV (Nov. - Dec. 15), PP 1-11 www.iosrjourals.org Numerical Solutios of Secod Order Boudary Value Problems
More informationWarped, Chirp Z-Transform: Radar Signal Processing
arped, Chirp Z-Trasform: Radar Sigal Processig by Garimella Ramamurthy Report o: IIIT/TR// Cetre for Commuicatios Iteratioal Istitute of Iformatio Techology Hyderabad - 5 3, IDIA Jauary ARPED, CHIRP Z
More informationBIBECHANA A Multidisciplinary Journal of Science, Technology and Mathematics
Mohd Yusuf Yasi / BIBECHANA 8 (2012) 31-36 : BMHSS, p. 31 BIBECHANA A Multidiscipliary Joural of Sciece, Techology ad Mathematics ISSN 2091-0762 (olie) Joural homepage: http://epjol.ifo/ide.php/bibechana
More informationFinite Difference Approximation for Transport Equation with Shifts Arising in Neuronal Variability
Iteratioal Joural of Sciece ad Research (IJSR) ISSN (Olie): 39-764 Ide Copericus Value (3): 64 Impact Factor (3): 4438 Fiite Differece Approimatio for Trasport Equatio with Shifts Arisig i Neuroal Variability
More informationFinite Difference Approximation for First- Order Hyperbolic Partial Differential Equation Arising in Neuronal Variability with Shifts
Iteratioal Joural of Scietific Egieerig ad Research (IJSER) wwwiseri ISSN (Olie): 347-3878, Impact Factor (4): 35 Fiite Differece Approimatio for First- Order Hyperbolic Partial Differetial Equatio Arisig
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationSolution of Differential Equation from the Transform Technique
Iteratioal Joural of Computatioal Sciece ad Mathematics ISSN 0974-3189 Volume 3, Number 1 (2011), pp 121-125 Iteratioal Research Publicatio House http://wwwirphousecom Solutio of Differetial Equatio from
More informationμ are complex parameters. Other
A New Numerical Itegrator for the Solutio of Iitial Value Problems i Ordiary Differetial Equatios. J. Suday * ad M.R. Odekule Departmet of Mathematical Scieces, Adamawa State Uiversity, Mubi, Nigeria.
More informationPOWER SERIES SOLUTION OF FIRST ORDER MATRIX DIFFERENTIAL EQUATIONS
Joural of Applied Mathematics ad Computatioal Mechaics 4 3(3) 3-8 POWER SERIES SOLUION OF FIRS ORDER MARIX DIFFERENIAL EQUAIONS Staisław Kukla Izabela Zamorska Istitute of Mathematics Czestochowa Uiversity
More informationLecture 8: Solving the Heat, Laplace and Wave equations using finite difference methods
Itroductory lecture otes o Partial Differetial Equatios - c Athoy Peirce. Not to be copied, used, or revised without explicit writte permissio from the copyright ower. 1 Lecture 8: Solvig the Heat, Laplace
More informationNumerical Solution of the Two Point Boundary Value Problems By Using Wavelet Bases of Hermite Cubic Spline Wavelets
Australia Joural of Basic ad Applied Scieces, 5(): 98-5, ISSN 99-878 Numerical Solutio of the Two Poit Boudary Value Problems By Usig Wavelet Bases of Hermite Cubic Splie Wavelets Mehdi Yousefi, Hesam-Aldie
More informationECE-S352 Introduction to Digital Signal Processing Lecture 3A Direct Solution of Difference Equations
ECE-S352 Itroductio to Digital Sigal Processig Lecture 3A Direct Solutio of Differece Equatios Discrete Time Systems Described by Differece Equatios Uit impulse (sample) respose h() of a DT system allows
More informationOrthogonal Gaussian Filters for Signal Processing
Orthogoal Gaussia Filters for Sigal Processig Mark Mackezie ad Kiet Tieu Mechaical Egieerig Uiversity of Wollogog.S.W. Australia Abstract A Gaussia filter usig the Hermite orthoormal series of fuctios
More informationStopping oscillations of a simple harmonic oscillator using an impulse force
It. J. Adv. Appl. Math. ad Mech. 5() (207) 6 (ISSN: 2347-2529) IJAAMM Joural homepage: www.ijaamm.com Iteratioal Joural of Advaces i Applied Mathematics ad Mechaics Stoppig oscillatios of a simple harmoic
More informationFREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING
Mechaical Vibratios FREE VIBRATION RESPONSE OF A SYSTEM WITH COULOMB DAMPING A commo dampig mechaism occurrig i machies is caused by slidig frictio or dry frictio ad is called Coulomb dampig. Coulomb dampig
More informationAlgebra II Notes Unit Seven: Powers, Roots, and Radicals
Syllabus Objectives: 7. The studets will use properties of ratioal epoets to simplify ad evaluate epressios. 7.8 The studet will solve equatios cotaiig radicals or ratioal epoets. b a, the b is the radical.
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordiary Differetial Equatios Braislav K. Nikolić Departmet of Physics ad Astroomy, Uiversity of Delaware, U.S.A. PHYS 460/660: Computatioal Methods of Physics http://www.physics.udel.edu/~bikolic/teachig/phys660/phys660.html
More informationStreamfunction-Vorticity Formulation
Streamfuctio-Vorticity Formulatio A. Salih Departmet of Aerospace Egieerig Idia Istitute of Space Sciece ad Techology, Thiruvaathapuram March 2013 The streamfuctio-vorticity formulatio was amog the first
More informationNumerical Solution of the First-Order Hyperbolic Partial Differential Equation with Point-Wise Advance
Iteratioal oural of Sciece ad Research (ISR) ISSN (Olie): 39-74 Ide Copericus Value (3): 4 Impact Factor (3): 4438 Numerical Solutio of the First-Order Hyperbolic Partial Differetial Equatio with Poit-Wise
More informationarxiv: v1 [cond-mat.stat-mech] 25 Jan 2017
Iteractig Multi-particle Classical Szilard Egie P. S. Pal a,b ad A. M. Jayaavar a,b a Istitute of Physics, Sachivalaya Marg, Bhubaeswar-75005, Idia b Homi Bhabha atioal Istitute, Traiig School Comple,
More informationPeriod Function of a Lienard Equation
Joural of Mathematical Scieces (4) -5 Betty Joes & Sisters Publishig Period Fuctio of a Lieard Equatio Khalil T Al-Dosary Departmet of Mathematics, Uiversity of Sharjah, Sharjah 77, Uited Arab Emirates
More informationPAijpam.eu ON TENSOR PRODUCT DECOMPOSITION
Iteratioal Joural of Pure ad Applied Mathematics Volume 103 No 3 2015, 537-545 ISSN: 1311-8080 (prited versio); ISSN: 1314-3395 (o-lie versio) url: http://wwwijpameu doi: http://dxdoiorg/1012732/ijpamv103i314
More informationNumerical Conformal Mapping via a Fredholm Integral Equation using Fourier Method ABSTRACT INTRODUCTION
alaysia Joural of athematical Scieces 3(1): 83-93 (9) umerical Coformal appig via a Fredholm Itegral Equatio usig Fourier ethod 1 Ali Hassa ohamed urid ad Teh Yua Yig 1, Departmet of athematics, Faculty
More informationMath 257: Finite difference methods
Math 257: Fiite differece methods 1 Fiite Differeces Remember the defiitio of a derivative f f(x + ) f(x) (x) = lim 0 Also recall Taylor s formula: (1) f(x + ) = f(x) + f (x) + 2 f (x) + 3 f (3) (x) +...
More informationChapter 9: Numerical Differentiation
178 Chapter 9: Numerical Differetiatio Numerical Differetiatio Formulatio of equatios for physical problems ofte ivolve derivatives (rate-of-chage quatities, such as velocity ad acceleratio). Numerical
More informationL 5 & 6: RelHydro/Basel. f(x)= ( ) f( ) ( ) ( ) ( ) n! 1! 2! 3! If the TE of f(x)= sin(x) around x 0 is: sin(x) = x - 3! 5!
aylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. At ay poit i the eighbourhood of =0, the fuctio ca be represeted as a power series of the followig form: X 0 f(a) f() ƒ() f()= ( ) f( ) (
More informationHALF-SWEEP GAUSS-SEIDEL ITERATION APPLIED TO TIME-FRACTIONAL DIFFUSION EQUATIONS
Jural Karya Asli Loreka Ahli Matematik Vol. 8 No. () Page 06-0 Jural Karya Asli Loreka Ahli Matematik HALF-SWEEP GAUSS-SEIDEL ITERATION APPLIED TO TIME-FRACTIONAL DIFFUSION EQUATIONS A. Suarto J. Sulaima
More informationExact Solutions for a Class of Nonlinear Singular Two-Point Boundary Value Problems: The Decomposition Method
Exact Solutios for a Class of Noliear Sigular Two-Poit Boudary Value Problems: The Decompositio Method Abd Elhalim Ebaid Departmet of Mathematics, Faculty of Sciece, Tabuk Uiversity, P O Box 741, Tabuki
More informationResearch Article Laplace Transform Collocation Method for Solving Hyperbolic Telegraph Equation
Hidawi Iteratioal Joural of Egieerig Mathematics Volume 7, Article ID 596, 9 pages https://doi.org/.55/7/596 Research Article Laplace Trasform Collocatio Method for Solvig Hyperbolic Telegraph Equatio
More informationExplicit Group Methods in the Solution of the 2-D Convection-Diffusion Equations
Proceedigs of the World Cogress o Egieerig 00 Vol III WCE 00 Jue 0 - July 00 Lodo U.K. Explicit Group Methods i the Solutio of the -D Covectio-Diffusio Equatios a Kah Bee orhashidah Hj. M. Ali ad Choi-Hog
More information. (24) If we consider the geometry of Figure 13 the signal returned from the n th scatterer located at x, y is
.5 SAR SIGNA CHARACTERIZATION I order to formulate a SAR processor we first eed to characterize the sigal that the SAR processor will operate upo. Although our previous discussios treated SAR cross-rage
More informationGenerating Functions for Laguerre Type Polynomials. Group Theoretic method
It. Joural of Math. Aalysis, Vol. 4, 2010, o. 48, 257-266 Geeratig Fuctios for Laguerre Type Polyomials α of Two Variables L ( xy, ) by Usig Group Theoretic method Ajay K. Shula* ad Sriata K. Meher** *Departmet
More informationFrequency Response of FIR Filters
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we itroduce the idea of the frequecy respose of LTI systems, ad focus specifically o the frequecy respose of FIR filters.. Steady-state
More informationPractical Spectral Anaysis (continue) (from Boaz Porat s book) Frequency Measurement
Practical Spectral Aaysis (cotiue) (from Boaz Porat s book) Frequecy Measuremet Oe of the most importat applicatios of the DFT is the measuremet of frequecies of periodic sigals (eg., siusoidal sigals),
More informationFREE VIBRATIONS OF SIMPLY SUPPORTED BEAMS USING FOURIER SERIES
Abdullah : FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS FREE VIBRATIONS OF SIMPY SUPPORTED BEAMS USING FOURIER SERIES SAWA MUBARAK ABDUAH Assistat ecturer Uiversity of Mosul Abstract Fourier series will be
More informationTaylor polynomial solution of difference equation with constant coefficients via time scales calculus
TMSCI 3, o 3, 129-135 (2015) 129 ew Treds i Mathematical Scieces http://wwwtmscicom Taylor polyomial solutio of differece equatio with costat coefficiets via time scales calculus Veysel Fuat Hatipoglu
More informationFinite Difference Derivations for Spreadsheet Modeling John C. Walton Modified: November 15, 2007 jcw
Fiite Differece Derivatios for Spreadsheet Modelig Joh C. Walto Modified: November 15, 2007 jcw Figure 1. Suset with 11 swas o Little Platte Lake, Michiga. Page 1 Modificatio Date: November 15, 2007 Review
More informationFive Steps Block Predictor-Block Corrector Method for the Solution of ( )
Applied Mathematics, 4,, -66 Published Olie May 4 i SciRes. http://www.scirp.org/oural/am http://dx.doi.org/.46/am.4.87 Five Steps Block Predictor-Block Corrector y = f x, y, y Method for the Solutio of
More informationSection 1 of Unit 03 (Pure Mathematics 3) Algebra
Sectio 1 of Uit 0 (Pure Mathematics ) Algebra Recommeded Prior Kowledge Studets should have studied the algebraic techiques i Pure Mathematics 1. Cotet This Sectio should be studied early i the course
More informationCourse Outline. Designing Control Systems. Proportional Controller. Amme 3500 : System Dynamics and Control. Root Locus. Dr. Stefan B.
Amme 3500 : System Dyamics ad Cotrol Root Locus Course Outlie Week Date Cotet Assigmet Notes Mar Itroductio 8 Mar Frequecy Domai Modellig 3 5 Mar Trasiet Performace ad the s-plae 4 Mar Block Diagrams Assig
More informationLainiotis filter implementation. via Chandrasekhar type algorithm
Joural of Computatios & Modellig, vol.1, o.1, 2011, 115-130 ISSN: 1792-7625 prit, 1792-8850 olie Iteratioal Scietific Press, 2011 Laiiotis filter implemetatio via Chadrasehar type algorithm Nicholas Assimais
More informationOn the convergence, consistence and stability of a standard finite difference scheme
AMERICAN JOURNAL OF SCIENTIFIC AND INDUSTRIAL RESEARCH 2, Sciece Huβ, ttp://www.sciub.org/ajsir ISSN: 253-649X, doi:.525/ajsir.2.2.2.74.78 O te covergece, cosistece ad stabilit of a stadard fiite differece
More informationSignal Processing in Mechatronics. Lecture 3, Convolution, Fourier Series and Fourier Transform
Sigal Processig i Mechatroics Summer semester, 1 Lecture 3, Covolutio, Fourier Series ad Fourier rasform Dr. Zhu K.P. AIS, UM 1 1. Covolutio Covolutio Descriptio of LI Systems he mai premise is that the
More informationRun-length & Entropy Coding. Redundancy Removal. Sampling. Quantization. Perform inverse operations at the receiver EEE
Geeral e Image Coder Structure Motio Video (s 1,s 2,t) or (s 1,s 2 ) Natural Image Samplig A form of data compressio; usually lossless, but ca be lossy Redudacy Removal Lossless compressio: predictive
More informationA Lattice Green Function Introduction. Abstract
August 5, 25 A Lattice Gree Fuctio Itroductio Stefa Hollos Exstrom Laboratories LLC, 662 Nelso Park Dr, Logmot, Colorado 853, USA Abstract We preset a itroductio to lattice Gree fuctios. Electroic address:
More informationStability Analysis of the Euler Discretization for SIR Epidemic Model
Stability Aalysis of the Euler Discretizatio for SIR Epidemic Model Agus Suryato Departmet of Mathematics, Faculty of Scieces, Brawijaya Uiversity, Jl Vetera Malag 6545 Idoesia Abstract I this paper we
More informationUNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA
UNDERWATER OBJECT CLASSIFICATION BY MEANS OF AN ACOUSTIC METHOD EUGENIUSZ KOZACZKA Naval Uiversity of Gdyia 81-13 Gdyia, Śmidowicza 69, Polad Gdańsk Uiversity of Techology 8-95 Gdańsk, Narutowicza 11/1,
More informationChapter 2 The Solution of Numerical Algebraic and Transcendental Equations
Chapter The Solutio of Numerical Algebraic ad Trascedetal Equatios Itroductio I this chapter we shall discuss some umerical methods for solvig algebraic ad trascedetal equatios. The equatio f( is said
More informationA STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD
IRET: Iteratioal oural of Research i Egieerig ad Techology eissn: 39-63 pissn: 3-7308 A STUDY ON MHD BOUNDARY LAYER FLOW OVER A NONLINEAR STRETCHING SHEET USING IMPLICIT FINITE DIFFERENCE METHOD Satish
More informationEXPERIMENT OF SIMPLE VIBRATION
EXPERIMENT OF SIMPLE VIBRATION. PURPOSE The purpose of the experimet is to show free vibratio ad damped vibratio o a system havig oe degree of freedom ad to ivestigate the relatioship betwee the basic
More informationAverage Number of Real Zeros of Random Fractional Polynomial-II
Average Number of Real Zeros of Radom Fractioal Polyomial-II K Kadambavaam, PG ad Research Departmet of Mathematics, Sri Vasavi College, Erode, Tamiladu, Idia M Sudharai, Departmet of Mathematics, Velalar
More informationTeaching Mathematics Concepts via Computer Algebra Systems
Iteratioal Joural of Mathematics ad Statistics Ivetio (IJMSI) E-ISSN: 4767 P-ISSN: - 4759 Volume 4 Issue 7 September. 6 PP-- Teachig Mathematics Cocepts via Computer Algebra Systems Osama Ajami Rashaw,
More informationKinetics of Complex Reactions
Kietics of Complex Reactios by Flick Colema Departmet of Chemistry Wellesley College Wellesley MA 28 wcolema@wellesley.edu Copyright Flick Colema 996. All rights reserved. You are welcome to use this documet
More informationCO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS
CO-LOCATED DIFFUSE APPROXIMATION METHOD FOR TWO DIMENSIONAL INCOMPRESSIBLE CHANNEL FLOWS C.PRAX ad H.SADAT Laboratoire d'etudes Thermiques,URA CNRS 403 40, Aveue du Recteur Pieau 86022 Poitiers Cedex,
More informationNewton Homotopy Solution for Nonlinear Equations Using Maple14. Abstract
Joural of Sciece ad Techology ISSN 9-860 Vol. No. December 0 Newto Homotopy Solutio for Noliear Equatios Usig Maple Nor Haim Abd. Rahma, Arsmah Ibrahim, Mohd Idris Jayes Faculty of Computer ad Mathematical
More informationApproximating the ruin probability of finite-time surplus process with Adaptive Moving Total Exponential Least Square
WSEAS TRANSACTONS o BUSNESS ad ECONOMCS S. Khotama, S. Boothiem, W. Klogdee Approimatig the rui probability of fiite-time surplus process with Adaptive Movig Total Epoetial Least Square S. KHOTAMA, S.
More informationNUMERICAL METHODS FOR SOLVING EQUATIONS
Mathematics Revisio Guides Numerical Methods for Solvig Equatios Page 1 of 11 M.K. HOME TUITION Mathematics Revisio Guides Level: GCSE Higher Tier NUMERICAL METHODS FOR SOLVING EQUATIONS Versio:. Date:
More informationTime-Domain Representations of LTI Systems
2.1 Itroductio Objectives: 1. Impulse resposes of LTI systems 2. Liear costat-coefficiets differetial or differece equatios of LTI systems 3. Bloc diagram represetatios of LTI systems 4. State-variable
More informationLinear Differential Equations of Higher Order Basic Theory: Initial-Value Problems d y d y dy
Liear Differetial Equatios of Higher Order Basic Theory: Iitial-Value Problems d y d y dy Solve: a( ) + a ( )... a ( ) a0( ) y g( ) + + + = d d d ( ) Subject to: y( 0) = y0, y ( 0) = y,..., y ( 0) = y
More informationTaylor expansion: Show that the TE of f(x)= sin(x) around. sin(x) = x - + 3! 5! L 7 & 8: MHD/ZAH
Taylor epasio: Let ƒ() be a ifiitely differetiable real fuctio. A ay poit i the eighbourhood of 0, the fuctio ƒ() ca be represeted by a power series of the followig form: X 0 f(a) f() f() ( ) f( ) ( )
More informationSolving a Nonlinear Equation Using a New Two-Step Derivative Free Iterative Methods
Applied ad Computatioal Mathematics 07; 6(6): 38-4 http://www.sciecepublishiggroup.com/j/acm doi: 0.648/j.acm.070606. ISSN: 38-5605 (Prit); ISSN: 38-563 (Olie) Solvig a Noliear Equatio Usig a New Two-Step
More information( ) (( ) ) ANSWERS TO EXERCISES IN APPENDIX B. Section B.1 VECTORS AND SETS. Exercise B.1-1: Convex sets. are convex, , hence. and. (a) Let.
Joh Riley 8 Jue 03 ANSWERS TO EXERCISES IN APPENDIX B Sectio B VECTORS AND SETS Exercise B-: Covex sets (a) Let 0 x, x X, X, hece 0 x, x X ad 0 x, x X Sice X ad X are covex, x X ad x X The x X X, which
More informationFind a formula for the exponential function whose graph is given , 1 2,16 1, 6
Math 4 Activity (Due by EOC Apr. ) Graph the followig epoetial fuctios by modifyig the graph of f. Fid the rage of each fuctio.. g. g. g 4. g. g 6. g Fid a formula for the epoetial fuctio whose graph is
More informationU8L1: Sec Equations of Lines in R 2
MCVU U8L: Sec. 8.9. Equatios of Lies i R Review of Equatios of a Straight Lie (-D) Cosider the lie passig through A (-,) with slope, as show i the diagram below. I poit slope form, the equatio of the lie
More informationVisco-acoustic Reverse-time Migration Using the Pseudo-analytical Method
Visco-acoustic Reverse-time Migratio Usig the Pseudo-aalytical Method J. Ramos-Martiez* (Petroleum Geo-Services), A.A. Valeciao (Petroleum Geo-Services) & M. Orlovich (Petroleum Geo-Services) SUMMARY We
More informationPARALEL PREDICTOR-CORRECTOR SCHEMES FOR PARABOLIC PROBLEMS ON GRAPHS
COMPUTATIONAL METHODS IN APPLIED MATHEMATICS, Vol. 10(2010, No. 3, pp. 275 282 c 2010 Istitute of Matematics of te Natioal Academy of Scieces of Belarus PARALEL PREDICTOR-CORRECTOR SCHEMES FOR PARABOLIC
More informationIN many scientific and engineering applications, one often
INTERNATIONAL JOURNAL OF COMPUTING SCIENCE AND APPLIED MATHEMATICS, VOL 3, NO, FEBRUARY 07 5 Secod Degree Refiemet Jacobi Iteratio Method for Solvig System of Liear Equatio Tesfaye Kebede Abstract Several
More informationAn Optimal Nearly Analytic Discrete Method for 2D Acoustic and Elastic Wave Equations
Bulleti of the Seismological Society of America, Vol. 94, No. 5, pp. 198 1991, October 004 A Optimal Nearly Aalytic Discrete Method for D Acoustic ad Elastic Wave Equatios by Dighui Yag, Mig Lu, Rusha
More informationComputational Tutorial of Steepest Descent Method and Its Implementation in Digital Image Processing
Computatioal utorial of Steepest Descet Method ad Its Implemetatio i Digital Image Processig Vorapoj Pataavijit Departmet of Electrical ad Electroic Egieerig, Faculty of Egieerig Assumptio Uiversity, Bagkok,
More informationSome New Iterative Methods for Solving Nonlinear Equations
World Applied Scieces Joural 0 (6): 870-874, 01 ISSN 1818-495 IDOSI Publicatios, 01 DOI: 10.589/idosi.wasj.01.0.06.830 Some New Iterative Methods for Solvig Noliear Equatios Muhammad Aslam Noor, Khalida
More informationFree Surface Hydrodynamics
Water Sciece ad Egieerig Free Surface Hydrodyamics y A part of Module : Hydraulics ad Hydrology Water Sciece ad Egieerig Dr. Shreedhar Maskey Seior Lecturer UNESCO-IHE Istitute for Water Educatio S. Maskey
More informationInversion of Earthquake Rupture Process:Theory and Applications
Iversio of Earthquake Rupture Process:Theory ad Applicatios Yu-tai CHEN 12 * Yog ZHANG 12 Li-sheg XU 2 1School of the Earth ad Space Scieces, Pekig Uiversity, Beijig 1871 2Istitute of Geophysics, Chia
More informationNew Version of the Rayleigh Schrödinger Perturbation Theory: Examples
New Versio of the Rayleigh Schrödiger Perturbatio Theory: Examples MILOŠ KALHOUS, 1 L. SKÁLA, 1 J. ZAMASTIL, 1 J. ČÍŽEK 2 1 Charles Uiversity, Faculty of Mathematics Physics, Ke Karlovu 3, 12116 Prague
More informationGUIDELINES ON REPRESENTATIVE SAMPLING
DRUGS WORKING GROUP VALIDATION OF THE GUIDELINES ON REPRESENTATIVE SAMPLING DOCUMENT TYPE : REF. CODE: ISSUE NO: ISSUE DATE: VALIDATION REPORT DWG-SGL-001 002 08 DECEMBER 2012 Ref code: DWG-SGL-001 Issue
More informationUsing Spreadsheets as a Computational Tool in Teaching Mechanical. Engineering
Proceedigs of the th WSEAS Iteratioal Coferece o COMPUTERS, Vouliagmei, Athes, Greece, July 3-5, 6 (pp35-3) Usig Spreadsheets as a Computatioal Tool i Teachig Mechaical Egieerig AHMADI-BROOGHANI, ZAHRA
More informationSignals & Systems Chapter3
Sigals & Systems Chapter3 1.2 Discrete-Time (D-T) Sigals Electroic systems do most of the processig of a sigal usig a computer. A computer ca t directly process a C-T sigal but istead eeds a stream of
More informationl -State Solutions of a New Four-Parameter 1/r^2 Singular Radial Non-Conventional Potential via Asymptotic Iteration Method
America Joural of Computatioal ad Applied Mathematics 8, 8(): 7-3 DOI:.593/j.ajcam.88. l -State Solutios of a New Four-Parameter /r^ Sigular Radial No-Covetioal Potetial via Asymptotic Iteratio Method
More informationSolution of Quantum Anharmonic Oscillator with Quartic Perturbation
ISS -79X (Paper) ISS 5-0638 (Olie) Vol.7, 0 Abstract Solutio of Quatum Aharmoic Oscillator with Quartic Perturbatio Adelaku A.O. Departmet of Physics, Wesley Uiversity of Sciece ad Techology, Odo, Odo
More information( a) ( ) 1 ( ) 2 ( ) ( ) 3 3 ( ) =!
.8,.9: Taylor ad Maclauri Series.8. Although we were able to fid power series represetatios for a limited group of fuctios i the previous sectio, it is ot immediately obvious whether ay give fuctio has
More informationFind quadratic function which pass through the following points (0,1),(1,1),(2, 3)... 11
Adrew Powuk - http://www.powuk.com- Math 49 (Numerical Aalysis) Iterpolatio... 4. Polyomial iterpolatio (system of equatio)... 4.. Lier iterpolatio... 5... Fid a lie which pass through (,) (,)... 8...
More informationSEQUENCES AND SERIES
Chapter 9 SEQUENCES AND SERIES Natural umbers are the product of huma spirit. DEDEKIND 9. Itroductio I mathematics, the word, sequece is used i much the same way as it is i ordiary Eglish. Whe we say that
More informationSNAP Centre Workshop. Basic Algebraic Manipulation
SNAP Cetre Workshop Basic Algebraic Maipulatio 8 Simplifyig Algebraic Expressios Whe a expressio is writte i the most compact maer possible, it is cosidered to be simplified. Not Simplified: x(x + 4x)
More informationNumerical Study on MHD Flow And Heat Transfer With The Effect Of Microrotational Parameter In The Porous Medium
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 4 (April. 5), V PP 8-7 www.iosrje.org Numerical Study o MHD Flow Ad Heat rasfer With he Effect Of Microrotatioal Parameter
More informationA New Method to Order Functions by Asymptotic Growth Rates Charlie Obimbo Dept. of Computing and Information Science University of Guelph
A New Method to Order Fuctios by Asymptotic Growth Rates Charlie Obimbo Dept. of Computig ad Iformatio Sciece Uiversity of Guelph ABSTRACT A ew method is described to determie the complexity classes of
More informationwavelet collocation method for solving integro-differential equation.
IOSR Joural of Egieerig (IOSRJEN) ISSN (e): 5-3, ISSN (p): 78-879 Vol. 5, Issue 3 (arch. 5), V3 PP -7 www.iosrje.org wavelet collocatio method for solvig itegro-differetial equatio. Asmaa Abdalelah Abdalrehma
More informationDiscrete-Time Signals and Systems. Discrete-Time Signals and Systems. Signal Symmetry. Elementary Discrete-Time Signals.
Discrete-ime Sigals ad Systems Discrete-ime Sigals ad Systems Dr. Deepa Kudur Uiversity of oroto Referece: Sectios. -.5 of Joh G. Proakis ad Dimitris G. Maolakis, Digital Sigal Processig: Priciples, Algorithms,
More informationA.1 Algebra Review: Polynomials/Rationals. Definitions:
MATH 040 Notes: Uit 0 Page 1 A.1 Algera Review: Polyomials/Ratioals Defiitios: A polyomial is a sum of polyomial terms. Polyomial terms are epressios formed y products of costats ad variales with whole
More informationModified Ratio Estimators Using Known Median and Co-Efficent of Kurtosis
America Joural of Mathematics ad Statistics 01, (4): 95-100 DOI: 10.593/j.ajms.01004.05 Modified Ratio s Usig Kow Media ad Co-Efficet of Kurtosis J.Subramai *, G.Kumarapadiya Departmet of Statistics, Podicherry
More informationSubstantiation of the Water Filling Theorem Using Lagrange Multipliers
J.A. Crawford U048 Water Fillig o Multiple Gaussia Chaels.doc Substatiatio of the Water Fillig heorem Usig Lagrage Multipliers Shao s capacity theorem states that i the case of parallel statistically idepedet
More informationLESSON 2: SIMPLIFYING RADICALS
High School: Workig with Epressios LESSON : SIMPLIFYING RADICALS N.RN.. C N.RN.. B 5 5 C t t t t t E a b a a b N.RN.. 4 6 N.RN. 4. N.RN. 5. N.RN. 6. 7 8 N.RN. 7. A 7 N.RN. 8. 6 80 448 4 5 6 48 00 6 6 6
More informationAP Calculus BC Review Applications of Derivatives (Chapter 4) and f,
AP alculus B Review Applicatios of Derivatives (hapter ) Thigs to Kow ad Be Able to Do Defiitios of the followig i terms of derivatives, ad how to fid them: critical poit, global miima/maima, local (relative)
More informationNumerical Solution of Non-linear Equations
Numerical Solutio of Noliear Equatios. INTRODUCTION The most commo reallife problems are oliear ad are ot ameable to be hadled by aalytical methods to obtai solutios of a variety of mathematical problems.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Electrical Engineering and Computer Science. BACKGROUND EXAM September 30, 2004.
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departmet of Electrical Egieerig ad Computer Sciece 6.34 Discrete Time Sigal Processig Fall 24 BACKGROUND EXAM September 3, 24. Full Name: Note: This exam is closed
More informationMathematical Series (You Should Know)
Mathematical Series You Should Kow Mathematical series represetatios are very useful tools for describig images or for solvig/approimatig the solutios to imagig problems. The may be used to epad a fuctio
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationNumber of fatalities X Sunday 4 Monday 6 Tuesday 2 Wednesday 0 Thursday 3 Friday 5 Saturday 8 Total 28. Day
LECTURE # 8 Mea Deviatio, Stadard Deviatio ad Variace & Coefficiet of variatio Mea Deviatio Stadard Deviatio ad Variace Coefficiet of variatio First, we will discuss it for the case of raw data, ad the
More information