Implicit Splitting Finite Difference Scheme for Multi-dimensional Wave Simulation
|
|
- Edward Flowers
- 5 years ago
- Views:
Transcription
1 Iplicit Splittig Fiite Differece Schee for Multi-diesioal Wave Siulatio Houhu (Jaes Zhag, Yu Zhag, CGGVeritas, Housto Jaes Su, CGGVeritas, Sigapore Suar I this abstract, we propose a ew fiite-differece schee for solvig wave equatios. This schee splits the ultidiesioal sste ito differet directios ad solves each directio iplicitl. Ulie ost splittig ethods i the literature which produce uerical aisotrop i diagoal directios, this ethod gives perfect circular ipulse resposes ad allows lateral velocit variatios. I this paper, we prove that the proposed schee is ucoditioall stable. I the uerical eaples, we show soe ipulse respose tests ad copare the with the results fro soe high-order eplicit fiite-differece ethods. The ew ethod allows larger tie step ad requires less eor storage durig the reverse tie igratio. Itroductio Seisic odelig usig fiite-differece (FD ethods is a ivaluable tool for seisic iterpretatio ad plas the e role i reverse tie igratio. Geerall speaig, published fiite-differece schees ca be divided ito two classes: iplicit FD schees (Mufti, 1985 ad eplicit FD schees. For 1D wave odelig, the iplicit schees are superior for two reasos. First the ca be ade ucoditioall stable so the tie step sie ca be arbitraril chose to save coputatio tie ad, secod, the lead to tridiagoal sstes which ca be solved efficietl b forward/bacward substitutio (Trefethe, However, i higher diesios, the iplicit FD schees require solvig ulti-diesioal atrices which greatl icreases eor usage ad coputatioal cost. Therefore, i ulti-diesios, seisic wave equatio odelig algoriths rel alost eclusivel o eplicit FD ethods (Kosloff, et.al., 198; Forberg Though eplicit FD schees are eas to solve, theoretical aalsis shows that the ca ol be coditioall stable ad the archig tie step sie is liited b the stabilit coditio. Such a restrictio put a theoretical upper boud to the coputatioal speed of this tpe of FD ethods. I this abstract, we propose a ew iplicit FD schee for solvig ulti-diesioal wave equatios. This schee splits the sste ito several oe-diesioal FD schees at differet spatial directios ad solves each directio iplicitl. Ulie ost of splittig ethods i the literature which produce uerical aisotrop i diagoal directios, this ethod gives perfect circles for ipulse respose testig. The ucoditioal stabilit has bee proved ad the dispersio has bee aaled i this paper. We also tested the ew ethods with differet ipulse resposes ad copared the results with those fro the high-order eplicit FD ethods. Siilar schee was proposed b Mufti (1985. But ours ca achieve higher order accurac i space ad is sipler to ipleet i higher diesios. Iplicit Splittig Fiite Differece Schee Let us start with the three-diesioal wave equatio: 1 p p p p = + +, (1 c (,, t where c (,, is the edia velocit. Deotig p = p( l, j,, we propose the followig Iplicit Splittig Fiite Differece (ISFD schee to solve equatio (1 uericall U δ (1 3 t V δ (1 3 t W δ (1 3 t where γ ad a are chose costats, + au + av + aw, ( δ is the cetral fiite differece operator δ p( = p( + p( + p(, ad the auiliar wavefields U, V ad W are defied as U = u p + u (3 V = v p + v W = w p + w +1 The pressure wavefield p ca be coputed as u + v + w p =. (4 3 The proposed ISFD ( cosists of three decoupled iplicit fiite-differece schees alog, ad directios. Each oe of the leads to a tri-diagoal atri ad ca be efficietl solved b a forward ad bacward substitutio ethod. Therefore, uericall FD schee ( is ore cost-effective tha the covetio iplicit FD schees. Most of the eistig directioal splittig ethods for SEG/Sa Atoio 007 Aual Meetig 011
2 Iplicit Splittig Fiite Differece to guaratee stabilit. Whe a is big copared with γ, (8 allows a larger tie step tha the covetioal eplicit fiite-differece schees. Figure 1: Dispersio error of the eplicit fiite-differece schees with d, 4 th, 6 th ad 8 th order i space ad of ISFD for differet a ad γ values. Tie i d order. solvig wave equatios itroduce uerical aisotrop, which causes big ieatic errors i odelig ad igratio. But our proposed ethod does ot suffer fro this proble. Furtherore, D ad 3D ipleetatios follows the sae uerical frae ad higher order accurac i space ca be achieved b adjustig a ad γ. These ae our schee differet fro that proposed b Mufti (1985. Stabilit ad Dispersio Aalsis To aale the stabilit ad dispersio, we first trasfor the fiite differece schee ( to tie ad waveuber doai as P P + P ~ ~ ~ = c ( + + P. (5 t Where is defied as ad ~ ~ θ + = θ 1 γ 3ac t θ, (6 4si ( / = θ. (7 Siilar epressios ca be writte for Λ ad Λ. It is eas to chec that equatio (5 is secod order i tie ad is a approiatio to the origial wave equatio (1. Therefore, o splittig error is itroduced i this ethod. We ca also prove that (5 is ucoditioall stable if a 1/ 4. For arbitraril give γ ad a, t ca be chose so that c t γ ( a To aale the spatial dispersio, igorig the errors caused b tie discretiatio i (5, we have the followig equatio ~ ~ ~ ω = c + +. (9 ( The above dispersio relatio is siilar to those of the covetioal high-order space eplicit FD schees ~ ~ ~ ω = c (, +, +,, (10 ~ ~ ~ where,, ad, are the spatial dispersio ters, of stadard th ~ -order FD schees. For eaples, is, the sae as θ defied b equatio (7, ad, ~ ~ 1, 4 =,(1 + si ( /, (11 3 etc.. I Figure (1, we copare the dispersio error for the covetioal spatial FD schees fro the d order to the 8 th order ad for ISFD for various a ad γ values. The relative dispersio error is defied as ~ / 1 ~, ad / 1. Ulie the covetioal FD schees, the dispersio relatio of ISFD ca be adjusted b choosig differet a ad γ values ad it also depeds o the odelig grid sie,, tie step t ad local velocit c. B properl choosig these paraeters, ISFD ca achieve a sae dispersio as higher order eplicit fiite-differece schees ad still use larger tie step if a is selected to be larger tha γ. Nuerical Eaples The proposed algorith was first tested o a ipulse respose. Figure (a shows the sap shot at 1s coputed usig the D eplicit fiite-differece with 8 th -order i space ad d -order i tie. Figure (b shows the sap shot at the sae tie step but coputed usig ISFD. I both cases, velocit is 000/s, t=0.005sec, =0 ad =0. ISFD has less dispersio o the frot edge tha the eplicit schee. A D dispersio error is plotted i Figure (3 for 8 th order eplicit schee ad i Figure (4 for ISFD usig the sae paraeters as for producig Figure (b. As i Figure (, the dispersio give b ISFD is ver close to 8 th -order fiite-differece but soother at high waveubers. Siilar results are show i Figure (5a ad Figure (5b ecept i this case = =40. Obviousl, ISFD has less dispersio. SEG/Sa Atoio 007 Aual Meetig 01
3 Iplicit Splittig Fiite Differece a Figure 4: Dispersio error of ISFD. The D dispersio error was plotted based the followig equatios: ~ ~ ~ ~,8 +,8 + E FD = 1, E = 1. (1 ISFD + + To test how well the algorith ca hadle coplicated geolog eviroets, we desiged a salt odel with cople geoetr ad high velocit. The odel is show i Figure (6. The velocit of the right salt is 4500/s. = =0. The source is located at (,=(4000,3000. Figure : The ipulse repose of the 8 th order eplicit fiite-differece schees (a ad the ISFD ethod (b. Tie is 1 secod. Spatial grid icreet is 0. b Eighth-order eplicit fiite differece becoes stable util t is reduced to 0.005sec. Figure (7 shows the sap shot at 1sec. Figure (8 shows the correspodig result with t=0.005sec usig ISFD. The eplicit schee ad ISFD have alost idetical ieatics. But the tie step for ISFD is twice as large as that for the eplicit ethod (ad ca be bigger thereb reducig coputatio cost. Coclusios Figure 3: Dispersio error of 8 th order eplicit fiitedifferece schee. We have proposed a iplicit splittig fiite differece schee for uerical solutio of ulti-diesioal wave equatio. Ulie ost splittig ethods, this ew schee provides good approiatio i all spatial directios ad this has bee validated b the ipulse resposes. B splittig a ulti-diesioal sste ito separate directios, each of the decoupled sste leads to a tridiagoal sste which ca be solved efficietl. We have show that the dispersio error of this schee ca be iiied b properl selectig odelig paraeters. Nuerical results show that the ew ethod, b usig larger tie step, ca give results of the sae qualit as those usig eplicit high-order fiite-differece schees. SEG/Sa Atoio 007 Aual Meetig 013
4 Iplicit Splittig Fiite Differece Acowledgets We tha CGGVeritas for the perissio to publish this paper ad David Yigst for a helpful discussios. Figure 6: A cople velocit odel. a Figure 7: The sap shot at 1 secod coputed usig the eplicit fiite-differece with 8 th -order i space. t=0.005sec. b Figure 5: The ipulse repose of the 8 th order eplicit fiitedifferece schees (a ad the ISFD ethod (b. Tie is 1 secod. Spatial grid icreet is 40. Figure 8: The sap shot at 1 secod coputed usig ISFD. t=0.005sec. SEG/Sa Atoio 007 Aual Meetig 014
5 EDITED REFERENCES Note: This referece list is a cop-edited versio of the referece list subitted b the author. Referece lists for the 007 SEG Techical Progra Epaded Abstracts have bee cop edited so that refereces provided with the olie etadata for each paper will achieve a high degree of liig to cited sources that appear o the Web. REFERENCES Dablai, M. A., 1986, The applicatio of high-order differecig to the scalar wave equatio: Geophsics, 51, Etge, J. T., 1988, Evaluatig fiite-differece operators applied to wave siulatio: Staford Eploratio Project 57. Forberg, B., 1987, The pseudospectral ethod: Coparisos with fiite-differece for the elastic wave equatio: Geophsics, 5, Kidela, M., A. Kael, ad P. Sguaero, 1990, O the costructio ad efficiec of staggered uerical differetiators for the wave equatio: Geophsics, 55, Kosloff, D. D., ad E. Basal, 198, Forward odelig b a Fourier ethod: Geophsics, 47, Mufti, I. R., 1985, Seisic odelig i the iplicit odel: Geophsical Prospectig, 33, Trefethe, L. N., 1996, Fiite-differece ad spectral ethods for ordiar ad partial differetial equatios, SEG/Sa Atoio 007 Aual Meetig 015
Lecture 11. Solution of Nonlinear Equations - III
Eiciecy o a ethod Lecture Solutio o Noliear Equatios - III The eiciecy ide o a iterative ethod is deied by / E r r: rate o covergece o the ethod : total uber o uctios ad derivative evaluatios at each step
More informationEnergy Identities of ADI-FDTD Method with Periodic Structure
Applied Matheatics 5 6 65-73 Published Olie Februar 5 i SciRes. http://www.scirp.org/oural/a http://d.doi.org/.436/a.5.65 erg Idetities of ADI-FDTD Method with Periodic Structure Regag Shi aitia Yag 3
More informationRAYLEIGH'S METHOD Revision D
RAYGH'S METHOD Revisio D B To Irvie Eail: toirvie@aol.co Noveber 5, Itroductio Daic sstes ca be characterized i ters of oe or ore atural frequecies. The atural frequec is the frequec at which the sste
More informationThe Differential Transform Method for Solving Volterra s Population Model
AASCIT Couicatios Volue, Issue 6 Septeber, 15 olie ISSN: 375-383 The Differetial Trasfor Method for Solvig Volterra s Populatio Model Khatereh Tabatabaei Departet of Matheatics, Faculty of Sciece, Kafas
More informationLebesgue Constant Minimizing Bivariate Barycentric Rational Interpolation
Appl. Math. If. Sci. 8, No. 1, 187-192 (2014) 187 Applied Matheatics & Iforatio Scieces A Iteratioal Joural http://dx.doi.org/10.12785/ais/080123 Lebesgue Costat Miiizig Bivariate Barycetric Ratioal Iterpolatio
More informationThe heat equation. But how to evaluate: F F? T T. We first transform this analytical equation into the corresponding form in finite-space:
L9_0/Cop. Astro.-, HS 0, Uiversität Basel/AAH The heat equatio T t T x. We first trasfor this aalytical equatio ito the correspodig for i fiite-space:.. t T δt T T F F F F F x x x x T T T F x But how to
More informationAdv. Studies Theor. Phys., Vol. 5, 2011, no. 4, A Numerical Computation of Non-Dimensional. Form of a Mathematical Model of Soil Salinity
Adv. Studies Theor. Phs. Vol. 5 11 o. 4 185-191 A Nuerical Coputatio of No-Diesioal For of a Matheatical Model of Soil Saliit Profile i a Rice Field ear Marie Shrip Aquaculture Far Nattawoot Pogoo 1 ad
More informationDiscrete population models
Discrete populatio odels D. Gurarie Ratioal: cclic (seasoal) tiig of reproductio ad developet, schroizatio Topics:. Reewal odels (Fiboacci). Discrete logistic odels (Verhulst vs. Ricker); cobwebs; equilibria,
More informationEvaluation of Bessel Functions Using a Computer Program
Evaluatio of Bessel Fuctios Usig a Coputer Progra P. S. Yeh, Ph.D. Abstract I cylidrical coordiate, there are two types of Bessel fuctios. These fuctios are the Bessel fuctio ad the odified Bessel fuctio.
More informationThe Binomial Multi-Section Transformer
4/15/2010 The Bioial Multisectio Matchig Trasforer preset.doc 1/24 The Bioial Multi-Sectio Trasforer Recall that a ulti-sectio atchig etwork ca be described usig the theory of sall reflectios as: where:
More informationAutomated Proofs for Some Stirling Number Identities
Autoated Proofs for Soe Stirlig Nuber Idetities Mauel Kauers ad Carste Scheider Research Istitute for Sybolic Coputatio Johaes Kepler Uiversity Altebergerstraße 69 A4040 Liz, Austria Subitted: Sep 1, 2007;
More informationCrank-Nicolson Implicit Method For The Nonlinear Schrodinger Equation With Variable Coefficient
Crak-Nicolso Iplicit Method For The Noliear Schrodiger Equatio With Variable Coefficiet Yaa Yee Choy a Wooi Nee Ta b Ki Gaik Tay c ad Chee Tiog Og d a Faculty of Sciece Techology & Hua Developet Uiversiti
More informationOrthogonal Function Solution of Differential Equations
Royal Holloway Uiversity of Loo Departet of Physics Orthogoal Fuctio Solutio of Differetial Equatios trouctio A give oriary ifferetial equatio will have solutios i ters of its ow fuctios Thus, for eaple,
More informationLecture 20 - Wave Propagation Response
.09/.093 Fiite Eleet Aalysis of Solids & Fluids I Fall 09 Lecture 0 - Wave Propagatio Respose Prof. K. J. Bathe MIT OpeCourseWare Quiz #: Closed book, 6 pages of otes, o calculators. Covers all aterials
More informationChapter 2. Asymptotic Notation
Asyptotic Notatio 3 Chapter Asyptotic Notatio Goal : To siplify the aalysis of ruig tie by gettig rid of details which ay be affected by specific ipleetatio ad hardware. [1] The Big Oh (O-Notatio) : It
More information5.6 Binomial Multi-section Matching Transformer
4/14/21 5_6 Bioial Multisectio Matchig Trasforers 1/1 5.6 Bioial Multi-sectio Matchig Trasforer Readig Assiget: pp. 246-25 Oe way to axiize badwidth is to costruct a ultisectio Γ f that is axially flat.
More informationBinomial transform of products
Jauary 02 207 Bioial trasfor of products Khristo N Boyadzhiev Departet of Matheatics ad Statistics Ohio Norther Uiversity Ada OH 4580 USA -boyadzhiev@ouedu Abstract Give the bioial trasfors { b } ad {
More information(s)h(s) = K( s + 8 ) = 5 and one finite zero is located at z 1
ROOT LOCUS TECHNIQUE 93 should be desiged differetly to eet differet specificatios depedig o its area of applicatio. We have observed i Sectio 6.4 of Chapter 6, how the variatio of a sigle paraeter like
More informationFUZZY RELIABILITY ANALYSIS OF COMPOUND SYSTEM BASED ON WEIBULL DISTRIBUTION
IJAMML 3:1 (2015) 31-39 Septeber 2015 ISSN: 2394-2258 Available at http://scietificadvaces.co.i DOI: http://dx.doi.org/10.18642/ijal_7100121530 FUZZY RELIABILITY ANALYSIS OF COMPOUND SYSTEM BASED ON WEIBULL
More informationX. Perturbation Theory
X. Perturbatio Theory I perturbatio theory, oe deals with a ailtoia that is coposed Ĥ that is typically exactly solvable of two pieces: a referece part ad a perturbatio ( Ĥ ) that is assued to be sall.
More informationStatistics and Data Analysis in MATLAB Kendrick Kay, February 28, Lecture 4: Model fitting
Statistics ad Data Aalysis i MATLAB Kedrick Kay, kedrick.kay@wustl.edu February 28, 2014 Lecture 4: Model fittig 1. The basics - Suppose that we have a set of data ad suppose that we have selected the
More informationLecture 10: Bounded Linear Operators and Orthogonality in Hilbert Spaces
Lecture : Bouded Liear Operators ad Orthogoality i Hilbert Spaces 34 Bouded Liear Operator Let ( X, ), ( Y, ) i i be ored liear vector spaces ad { } X Y The, T is said to be bouded if a real uber c such
More informationUncertainty Principle of Mathematics
Septeber 27 Ucertaity Priciple of Matheatics Shachter Mourici Israel, Holo ourici@walla.co.il Preface This short paper prove that atheatically, Reality is ot real. This short paper is ot about Heiseberg's
More informationThe Binomial Multi- Section Transformer
4/4/26 The Bioial Multisectio Matchig Trasforer /2 The Bioial Multi- Sectio Trasforer Recall that a ulti-sectio atchig etwork ca be described usig the theory of sall reflectios as: where: ( ω ) = + e +
More information5.6 Binomial Multi-section Matching Transformer
4/14/2010 5_6 Bioial Multisectio Matchig Trasforers 1/1 5.6 Bioial Multi-sectio Matchig Trasforer Readig Assiget: pp. 246-250 Oe way to axiize badwidth is to costruct a ultisectio Γ f that is axially flat.
More informationFourth Order Positively Smoothed Padé Schemes. for Parabolic Partial Differential Equations with. Nonlocal Boundary Conditions
Applied Matheatical Scieces, Vol. 4, 200, o. 42, 2065-2080 Fourth Order Positively Soothed Padé Schees for Parabolic Partial Differetial Equatios with Nolocal Boudary Coditios Mohaad Siddique Departet
More information2D DSP Basics: Systems Stability, 2D Sampling
- Digital Iage Processig ad Copressio D DSP Basics: Systes Stability D Saplig Stability ty Syste is stable if a bouded iput always results i a bouded output BIBO For LSI syste a sufficiet coditio for stability:
More informationVerification of continuous predictands
barbara.casati@ec.gc.ca Verificatio of cotiuous predictads Barbara Casati 9 Ja 007 Exploratory ethods: joit distributio Scatter-plot: plot of observatio versus forecast values Perfect forecast obs, poits
More information1 xn, and wavenumber index m as follows. ( m =0,,M-1, M denotes the total sample points of the interval) space axis, g (-)
Drealet prestack depth igratio usig Local Cosie Basis ad Local Expoetial Fraes Bagyu Wu *, Ru-Sha Wu, ad Jighuai Gao Modelig ad Iagig Laboratory, IGPP, Uiversity of Califoria, Sata Cruz, CA 95064 Suary
More informationx+ 2 + c p () x c p () x is an arbitrary function. ( sin x ) dx p f() d d f() dx = x dx p cosx = cos x+ 2 d p () x + x-a r (1.
Super Derivative (No-iteger ties Derivative). Super Derivative ad Super Differetiatio Defitio.. p () obtaied by cotiuig aalytically the ide of the differetiatio operator of Higher Derivative of a fuctio
More informationA New Family of Multipoint Iterative Methods for Finding Multiple Roots of Nonlinear Equations
Aerica Joural o Coputatioal ad Applied Matheatics 3, 3(3: 68-73 DOI:.593/j.ajca.333.3 A New Faily o Multipoit Iterative Methods or Fidig Multiple Roots o Noliear Equatios R. Thukral Padé Research Cetre,
More informationFinite element analysis of nonlinear structures with Newmark method
Iteratioal Joural of the Physical Scieces Vol. 6(6), 95-40, 8 March, 0 Available olie at http://www.acadeicjourals.org/ijps ISSN 99-950 0 Acadeic Jourals Full Legth Research Paper Fiite eleet aalysis of
More informationAPPLICATION OF UNCERTAIN NONLINEAR SYSTEMS PARTIAL STATE VARIABLES CONTROL TO A CLASS OF PENDULUM SYSTEMS
3 st Deceber 0 Vol 46 No 005-0 JATIT & LLS All rights reserved ISSN: 99-8645 wwwjatitorg E-ISSN: 87-395 APPLICATION OF UNCERTAIN NONLINEAR SYSTEMS PARTIAL STATE VARIABLES CONTROL TO A CLASS OF PENDULUM
More informationExercise 8 CRITICAL SPEEDS OF THE ROTATING SHAFT
Exercise 8 CRITICA SEEDS OF TE ROTATING SAFT. Ai of the exercise Observatio ad easureet of three cosecutive critical speeds ad correspodig odes of the actual rotatig shaft. Copariso of aalytically coputed
More informationAfter the completion of this section the student
Chapter II CALCULUS II. Liits ad Cotiuity 55 II. LIMITS AND CONTINUITY Objectives: After the copletio of this sectio the studet - should recall the defiitios of the it of fuctio; - should be able to apply
More informationReliability Equivalence Analysis of a Parallel-Series System Subject to Degradation Facility
Sciece Joural of Applied Mateatics ad Statistics 5; 3(3): 6-64 Publised olie Jue 6 5 (ttp://www.sciecepublisiggroup.co/j/sjas) doi:.648/j.sjas.533.9 ISSN: 376-949 (Prit); ISSN: 376-953 (Olie) Reliability
More informationMath 4707 Spring 2018 (Darij Grinberg): midterm 2 page 1. Math 4707 Spring 2018 (Darij Grinberg): midterm 2 with solutions [preliminary version]
Math 4707 Sprig 08 Darij Griberg: idter page Math 4707 Sprig 08 Darij Griberg: idter with solutios [preliiary versio] Cotets 0.. Coutig first-eve tuples......................... 3 0.. Coutig legal paths
More informationApplication 10.5D Spherical Harmonic Waves
Applicatio 10.5D Spherical Haroic Waves I probles ivolvig regios that ejoy spherical syetry about the origi i space, it is appropriate to use spherical coordiates. The 3-diesioal Laplacia for a fuctio
More informationx !1! + 1!2!
4 Euler-Maclauri Suatio Forula 4. Beroulli Nuber & Beroulli Polyoial 4.. Defiitio of Beroulli Nuber Beroulli ubers B (,,3,) are defied as coefficiets of the followig equatio. x e x - B x! 4.. Expreesio
More informationChapter 9 Computation of the Discrete. Fourier Transform
Chapter 9 Coputatio of the Discrete Fourier Trasfor Itroductio Efficiet Coputatio of the Discrete Fourier Trasfor Goertzel Algorith Deciatio-I-Tie FFT Algoriths Deciatio-I-Frequecy FFT Algoriths Ipleetatio
More informationNew Simple Methods of Tuning Three-Term Controllers for Dead-Time Processes
New Siple Methods of Tuig Three-Ter Cotrollers for Dead-Tie Processes.G.ARVANTS Departet of Natural Resources Maageet ad Agricultural Egieerig Agricultural Uiversity of Athes era Odos 75, Botaikos 8 55,
More informationNumerical Solution of Nonlinear Integro-Differential Equations with Initial Conditions by Bernstein Operational Matrix of Derivative
Iteratioal Joural of Moder Noliear heory ad Applicatio 3-9 http://ddoiorg/36/ijta38 Published Olie Jue 3 (http://wwwscirporg/joural/ijta) Nuerical Solutio of Noliear Itegro-Differetial Equatios with Iitial
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 informationAVERAGE MARKS SCALING
TERTIARY INSTITUTIONS SERVICE CENTRE Level 1, 100 Royal Street East Perth, Wester Australia 6004 Telephoe (08) 9318 8000 Facsiile (08) 95 7050 http://wwwtisceduau/ 1 Itroductio AVERAGE MARKS SCALING I
More informationOn Exact Finite-Difference Scheme for Numerical Solution of Initial Value Problems in Ordinary Differential Equations.
O Exact Fiite-Differece Sceme for Numerical Solutio of Iitial Value Problems i Ordiar Differetial Equatios. Josua Suda, M.Sc. Departmet of Matematical Scieces, Adamawa State Uiversit, Mubi, Nigeria. E-mail:
More informationA Tabu Search Method for Finding Minimal Multi-Homogeneous Bézout Number
Joural of Matheatics ad Statistics 6 (): 105-109, 010 ISSN 1549-3644 010 Sciece Publicatios A Tabu Search Method for Fidig Miial Multi-Hoogeeous Bézout Nuber Hassa M.S. Bawazir ad Ali Abd Raha Departet
More information) is a square matrix with the property that for any m n matrix A, the product AI equals A. The identity matrix has a ii
square atrix is oe that has the sae uber of rows as colus; that is, a atrix. he idetity atrix (deoted by I, I, or [] I ) is a square atrix with the property that for ay atrix, the product I equals. he
More informationA New Type of q-szász-mirakjan Operators
Filoat 3:8 07, 567 568 https://doi.org/0.98/fil7867c Published by Faculty of Scieces ad Matheatics, Uiversity of Niš, Serbia Available at: http://www.pf.i.ac.rs/filoat A New Type of -Szász-Miraka Operators
More informationSummer MA Lesson 13 Section 1.6, Section 1.7 (part 1)
Suer MA 1500 Lesso 1 Sectio 1.6, Sectio 1.7 (part 1) I Solvig Polyoial Equatios Liear equatio ad quadratic equatios of 1 variable are specific types of polyoial equatios. Soe polyoial equatios of a higher
More informationPerturbation Theory, Zeeman Effect, Stark Effect
Chapter 8 Perturbatio Theory, Zeea Effect, Stark Effect Ufortuately, apart fro a few siple exaples, the Schrödiger equatio is geerally ot exactly solvable ad we therefore have to rely upo approxiative
More informationASYMPTOTIC STABILITY OF POSITIVE FRACTIONAL 2D LINEAR SYSTEMS WITH DELAYS
ASYMPTOTIC STABILITY OF POSITIVE FRACTIONAL D LINEAR SYSTEMS WITH DELAYS Tadeusz Kaczorek Faculty of Electrical Egieerig Bialystok Techical Uiversity Wiejska 45D 5-35 Bialystok Polad kaczorek@isep.pw.edu.pl
More informationDefine a Markov chain on {1,..., 6} with transition probability matrix P =
Pla Group Work 0. The title says it all Next Tie: MCMC ad Geeral-state Markov Chais Midter Exa: Tuesday 8 March i class Hoework 4 due Thursday Uless otherwise oted, let X be a irreducible, aperiodic Markov
More informationECE Spring Prof. David R. Jackson ECE Dept. Notes 20
ECE 6341 Sprig 016 Prof. David R. Jackso ECE Dept. Notes 0 1 Spherical Wave Fuctios Cosider solvig ψ + k ψ = 0 i spherical coordiates z φ θ r y x Spherical Wave Fuctios (cot.) I spherical coordiates we
More informationWe have also learned that, thanks to the Central Limit Theorem and the Law of Large Numbers,
Cofidece Itervals III What we kow so far: We have see how to set cofidece itervals for the ea, or expected value, of a oral probability distributio, both whe the variace is kow (usig the stadard oral,
More informationAppendix: The Laplace Transform
Appedix: The Laplace Trasform The Laplace trasform is a powerful method that ca be used to solve differetial equatio, ad other mathematical problems. Its stregth lies i the fact that it allows the trasformatio
More informationME 501A Seminar in Engineering Analysis Page 1
Accurac, Stabilit ad Sstems of Equatios November 0, 07 Numerical Solutios of Ordiar Differetial Equatios Accurac, Stabilit ad Sstems of Equatios Larr Caretto Mecaical Egieerig 0AB Semiar i Egieerig Aalsis
More informationApril 1980 TR/96. Extrapolation techniques for first order hyperbolic partial differential equations. E.H. Twizell
TR/96 Apri 980 Extrapoatio techiques for first order hyperboic partia differetia equatios. E.H. Twize W96086 (0) 0. Abstract A uifor grid of step size h is superiposed o the space variabe x i the first
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 informationCOMP 2804 Solutions Assignment 1
COMP 2804 Solutios Assiget 1 Questio 1: O the first page of your assiget, write your ae ad studet uber Solutio: Nae: Jaes Bod Studet uber: 007 Questio 2: I Tic-Tac-Toe, we are give a 3 3 grid, cosistig
More informationA Class of Blended Block Second Derivative Multistep Methods for Stiff Systems
Iteratioal Joural of Iovative Mathematics, Statistics & Eerg Policies ():-6, Ja.-Mar. 7 SEAHI PUBLICATIONS, 7 www.seahipa.org ISSN: 67-8X A Class of Bleded Bloc Secod Derivative Multistep Methods for Stiff
More informationFIR Filter Design: Part II
EEL335: Discrete-Time Sigals ad Systems. Itroductio I this set of otes, we cosider how we might go about desigig FIR filters with arbitrary frequecy resposes, through compositio of multiple sigle-peak
More informationThe Hypergeometric Coupon Collection Problem and its Dual
Joural of Idustrial ad Systes Egieerig Vol., o., pp -7 Sprig 7 The Hypergeoetric Coupo Collectio Proble ad its Dual Sheldo M. Ross Epstei Departet of Idustrial ad Systes Egieerig, Uiversity of Souther
More informationand then substitute this into the second equation to get 5(11 4 y) 3y
Math E-b Lecture # Notes The priary focus of this week s lecture is a systeatic way of solvig ad uderstadig systes of liear equatios algebraically, geoetrically, ad logically. Eaple #: Solve the syste
More informationAcoustic Field inside a Rigid Cylinder with a Point Source
Acoustic Field iside a Rigid Cylider with a Poit Source 1 Itroductio The ai objectives of this Deo Model are to Deostrate the ability of Coustyx to odel a rigid cylider with a poit source usig Coustyx
More informationBertrand s postulate Chapter 2
Bertrad s postulate Chapter We have see that the sequece of prie ubers, 3, 5, 7,... is ifiite. To see that the size of its gaps is ot bouded, let N := 3 5 p deote the product of all prie ubers that are
More informationData Analysis and Statistical Methods Statistics 651
Data Aalysis ad Statistical Methods Statistics 651 http://www.stat.tau.edu/~suhasii/teachig.htl Suhasii Subba Rao Exaple The itroge cotet of three differet clover plats is give below. 3DOK1 3DOK5 3DOK7
More informationPARTIAL DIFFERENTIAL EQUATIONS SEPARATION OF VARIABLES
Diola Bagayoko (0 PARTAL DFFERENTAL EQUATONS SEPARATON OF ARABLES. troductio As discussed i previous lectures, partial differetial equatios arise whe the depedet variale, i.e., the fuctio, varies with
More informationContents Two Sample t Tests Two Sample t Tests
Cotets 3.5.3 Two Saple t Tests................................... 3.5.3 Two Saple t Tests Setup: Two Saples We ow focus o a sceario where we have two idepedet saples fro possibly differet populatios. Our
More informationApplication of Homotopy Analysis Method for Solving various types of Problems of Ordinary Differential Equations
Iteratioal Joural o Recet ad Iovatio Treds i Coputig ad Couicatio IN: 31-8169 Volue: 5 Issue: 5 16 Applicatio of Hootopy Aalysis Meod for olvig various types of Probles of Ordiary Differetial Equatios
More informationCEE 522 Autumn Uncertainty Concepts for Geotechnical Engineering
CEE 5 Autum 005 Ucertaity Cocepts for Geotechical Egieerig Basic Termiology Set A set is a collectio of (mutually exclusive) objects or evets. The sample space is the (collectively exhaustive) collectio
More informationGoal. Adaptive Finite Element Methods for Non-Stationary Convection-Diffusion Problems. Outline. Differential Equation
Goal Adaptive Fiite Elemet Methods for No-Statioary Covectio-Diffusio Problems R. Verfürth Ruhr-Uiversität Bochum www.ruhr-ui-bochum.de/um1 Tübige / July 0th, 017 Preset space-time adaptive fiite elemet
More informationRecursive computation of spherical harmonic rotation coefficients of large degree
Recursive coputatio of spherical haroic rotatio coefficiets of large degree Nail A. Guerov 1 ad Raai Duraiswai 1, 1 Istitute for Advaced Coputer Studies Departet of Coputer Sciece, Uiversity of Marylad,
More informationThe Jordan Normal Form: A General Approach to Solving Homogeneous Linear Systems. Mike Raugh. March 20, 2005
The Jorda Normal Form: A Geeral Approach to Solvig Homogeeous Liear Sstems Mike Raugh March 2, 25 What are we doig here? I this ote, we describe the Jorda ormal form of a matrix ad show how it ma be used
More informationS. A. ALIEV, Y. I. YELEYKO, Y. V. ZHERNOVYI. STEADY-STATE DISTRIBUTIONS FOR CERTAIN MODIFICATIONS OF THE M/M/1/m QUEUEING SYSTEM
Trasactios of Azerbaija Natioal Acadey of Scieces, Series of Physical-Techical ad Matheatical Scieces: Iforatics ad Cotrol Probles 009 Vol XXIX, 6 P 50-58 S A ALIEV, Y I YELEYKO, Y V ZHERNOVYI STEADY-STATE
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 informationGENERALIZED KERNEL AND MIXED INTEGRAL EQUATION OF FREDHOLM - VOLTERRA TYPE R. T. Matoog
GENERALIZED KERNEL AND MIXED INTEGRAL EQUATION OF FREDHOLM - VOLTERRA TYPE R. T. Matoog Assistat Professor. Deartet of Matheatics, Faculty of Alied Scieces,U Al-Qura Uiversity, Makkah, Saudi Arabia Abstract:
More informationarxiv: v1 [cs.ds] 23 Nov 2012
Aalysis of a radoized approxiatio schee for atrix ultiplicatio Daiel Hsu 1, Sha M. Kakade 1, ad Tog Zhag 2 arxiv:1211.5414v1 cs.ds] 23 Nov 2012 1 Microsoft Research, New Eglad 2 Departet of Statistics,
More informationOn the Derivation and Implementation of a Four Stage Harmonic Explicit Runge-Kutta Method *
Applied Mathematics, 05, 6, 694-699 Published Olie April 05 i SciRes. http://www.scirp.org/joural/am http://dx.doi.org/0.46/am.05.64064 O the Derivatio ad Implemetatio of a Four Stage Harmoic Explicit
More informationWavelet Transform Theory. Prof. Mark Fowler Department of Electrical Engineering State University of New York at Binghamton
Wavelet Trasfor Theory Prof. Mark Fowler Departet of Electrical Egieerig State Uiversity of New York at Bighato What is a Wavelet Trasfor? Decopositio of a sigal ito costituet parts Note that there are
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 informationAl Lehnen Madison Area Technical College 10/5/2014
The Correlatio of Two Rado Variables Page Preliiary: The Cauchy-Schwarz-Buyakovsky Iequality For ay two sequeces of real ubers { a } ad = { b } =, the followig iequality is always true. Furtherore, equality
More informationSequences, Mathematical Induction, and Recursion. CSE 2353 Discrete Computational Structures Spring 2018
CSE 353 Discrete Computatioal Structures Sprig 08 Sequeces, Mathematical Iductio, ad Recursio (Chapter 5, Epp) Note: some course slides adopted from publisher-provided material Overview May mathematical
More informationMost text will write ordinary derivatives using either Leibniz notation 2 3. y + 5y= e and y y. xx tt t
Itroductio to Differetial Equatios Defiitios ad Termiolog Differetial Equatio: A equatio cotaiig the derivatives of oe or more depedet variables, with respect to oe or more idepedet variables, is said
More informationObservations on Derived K-Fibonacci and Derived K- Lucas Sequences
ISSN(Olie): 9-875 ISSN (Prit): 7-670 Iteratioal Joural of Iovative Research i Sciece Egieerig ad Techology (A ISO 97: 007 Certified Orgaizatio) Vol. 5 Issue 8 August 06 Observatios o Derived K-iboacci
More informationis also known as the general term of the sequence
Lesso : Sequeces ad Series Outlie Objectives: I ca determie whether a sequece has a patter. I ca determie whether a sequece ca be geeralized to fid a formula for the geeral term i the sequece. I ca determie
More informationPartial Differential Equations
Partial Differetial Equatios Part 2 Massimo Ricotti ricotti@astro.umd.edu Uiversity of Marylad Partial Differetial Equatios p.1/15 Upwid differecig I additio to amplitude errors (istability or dampig),
More informationChapter 2: Numerical Methods
Chapter : Numerical Methods. Some Numerical Methods for st Order ODEs I this sectio, a summar of essetial features of umerical methods related to solutios of ordiar differetial equatios is give. I geeral,
More informationTHE CLOSED FORMS OF CONVERGENT INFINITE SERIES ESTIMATION OF THE SERIES SUM OF NON-CLOSED FORM ALTERNATING SERIES TO A HIGH DEGREE OF PRECISION.
THE CLSED FRMS F CNERGENT INFINITE SERIES ESTIMATIN F THE SERIES SUM F NN-CLSED FRM ALTERNATING SERIES T A HIGH DEGREE F PRECISIN. Peter G.Bass. PGBass M er..0.0. www.relativityoais.co May 0 Abstract This
More informationECE 901 Lecture 4: Estimation of Lipschitz smooth functions
ECE 9 Lecture 4: Estiatio of Lipschitz sooth fuctios R. Nowak 5/7/29 Cosider the followig settig. Let Y f (X) + W, where X is a rado variable (r.v.) o X [, ], W is a r.v. o Y R, idepedet of X ad satisfyig
More information2. F ; =(,1)F,1; +F,1;,1 is satised by thestirlig ubers of the rst kid ([1], p. 824). 3. F ; = F,1; + F,1;,1 is satised by the Stirlig ubers of the se
O First-Order Two-Diesioal Liear Hoogeeous Partial Dierece Equatios G. Neil Have y Ditri A. Gusev z Abstract Aalysis of algoriths occasioally requires solvig of rst-order two-diesioal liear hoogeeous partial
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 informationA New Hybrid in the Nonlinear Part of Adomian Decomposition Method for Initial Value Problem of Ordinary Differential Equation
Joural of Matematics Researc; Vol No ; ISSN - E-ISSN - Publised b Caadia Ceter of Sciece ad Educatio A New Hbrid i te Noliear Part of Adomia Decompositio Metod for Iitial Value Problem of Ordiar Differetial
More informationNull Spaces, Column. Transformations. Remarks. Definition ( ) ( ) Null A is the set of all x's in R that go to 0 in R 3
Null Spaces Colu Spaces ad Liear rasforatios earks Whe we study fuctios f ( ) we fid that a uderstadig of the doai D f itervals of - values ad the rage f itervals of y - values are ecessary to descrie
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 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 informationLecture 19. Curve fitting I. 1 Introduction. 2 Fitting a constant to measured data
Lecture 9 Curve fittig I Itroductio Suppose we are preseted with eight poits of easured data (x i, y j ). As show i Fig. o the left, we could represet the uderlyig fuctio of which these data are saples
More information8.3 Perturbation theory
8.3 Perturbatio theory Slides: Video 8.3.1 Costructig erturbatio theory Text referece: Quatu Mechaics for Scietists ad gieers Sectio 6.3 (u to First order erturbatio theory ) Perturbatio theory Costructig
More information3D acoustic wave modeling with a time-space domain dispersion-relation-based Finite-difference scheme
P-8 3D acoustic wave odeling with a tie-space doain dispersion-relation-based Finite-difference schee Yang Liu * and rinal K. Sen State Key Laboratory of Petroleu Resource and Prospecting (China University
More information6.867 Machine learning, lecture 13 (Jaakkola)
Lecture topics: Boostig, argi, ad gradiet descet copleity of classifiers, geeralizatio Boostig Last tie we arrived at a boostig algorith for sequetially creatig a eseble of base classifiers. Our base classifiers
More informationLinear Programming and the Simplex Method
Liear Programmig ad the Simplex ethod Abstract This article is a itroductio to Liear Programmig ad usig Simplex method for solvig LP problems i primal form. What is Liear Programmig? Liear Programmig is
More information