The convolution computation for Perfectly Matched Boundary Layer algorithm in finite differences
|
|
- Alfred Osborne
- 5 years ago
- Views:
Transcription
1 The convoluton computaton for Perfectly Matched Boundary Layer algorthm n fnte dfferences Herman Jaramllo May 10, Introducton Ths s an exercse to help on the understandng on some mportant ssues on CPML) Convolutonal Perfect Matchng Layer algorthms to attenuate edge effects n fnte dfferences. From the many references, I cte Roden and Gedney s [4], and Pasalc and McGarry [2]. Also, some good notes 1 by Steven Johnson. 2 The Stretchng Factor The dea well explaned by Steven Johnson s as follows. There are two zones. The nsde zone where we want to preserve the numercs and a buffer zone the boundares and ts surroundngs) where we want to attenuate reflectons so that the numercal experments looks lke smulated n a free space medum wth boundares at nfnty. At the buffer zone some dstance away from the center), the wave equatons are analytcally contnued nto the upper or lower, dependng on the Fourer transform beng studed. For causal functons the sgns of the k vector and ω are opposte. That s, a plane wave s wrte n ts smplest form as e k x ωt.) half complex plane so that the oscllaton character changes to an 1 math.mt.edu/ stevenj/18.369/pml.pdf 1
2 exponental dampng character. The trck s done by the change of varable assumng 1D for smplcty, but the dea s general) x x + fx) 2.1) where 1. The functon f acts as the dampng factor. When f 0 then the problem s not changed. Ths occurs n the zone of nterest. Based on physcal ntuton f s chosen so that df dx σ xx), 2.2) ω and ths s explaned next. The reason for the ω n the denomnator s that ths electon wll make the problem ndependent of frequency for a gven reflecton angle. That s, t wll avod dspersve frequency dependent) behavor. To see ths, we observe from equatons 2.1 and 2.2 that xx) x + ω x σ x x )dx, wth the coordnate where the dampng layer starts. See for example that for f 0 n the new coordnate system a plane wave along the x wth k k x k 1 ) drecton, n the new coordnate x, could be represented as e k x e kx e k ω x σ xx )dx and snce c p k/ω s the phase velocty whch s constant n non dspersve meda) we see that the attenuaton represented by the factor e k ω x σ xx )dx. s ndependent of frequency for a gven propagaton drecton). To be sure that the attenuaton factor works correctly, we need to guarantee that x σ x x )dx s postve for postve k, assumng only postve frequences ω. Now, f the plane wave s defned as e k x+ωt the sgns for and k should be opposte). Note that f the wave s travelng along the negatve drecton k < 0 and the ntegral s negatve snce the upper ndex s smaller than the lower ndex and then stll the ntegral contrbuton s negatve, so stll the exponent s negatve) then also here the factor s a dampng factor. 2
3 The smplest electon of σ x s a constant n whch case we would get the mappng x 1 + σ ) x s x x ω where s x s called a stretchng factor. Roden and Gedney s [4] use a more general stretchng factor stll wth constant), s +, wth α > 0 and 1. Whle the factor can be seen as an overall scalng, the factor α s justfed as shft to the pole away from the real axs, whch wll mprove accuracy for grazng ncdence due to source proxmty to the boundary or large offsets. Martn et. al., [3] reaffrm Roden and Gedney s [4] statements. Here and n what follows the subndex means x, f 1, y f 2 and z f 3; for the three dmensonal space. Note that n ther work the sgn conventon s opposte, snce the magnary unt s n the denomnator. They descrbe the plane wave wth a mnus - sgn n front of the phase φ k x ωt. Martn et. al., [3] use the stretchng factor s x 1 + σ x ω where, as n Roden and Gedney, the sgn conventon s opposte to the one used here. Komattsch and Martn [1] also clam that wth the help of the new parameter α ntroduced by Roden and Gedney, better accuracy s obtaned at grazng ncdence angles. Snce a product n the frequency/wavenumber doman s convoluton n tme/space doman, and the FDTD Fnte Dfference Tme Doman) wave equaton s mplemented n tme doman, we need to fnd the nverse Fourer Laplace) transform of the stretchng factor. After fndng the stretchng factor we should apply convoluton. Next secton deals wth fndng the nverse Laplace transform of the stretchng factor, and the fnal secton shows a recurson formula that speeds up the convolutonal mplementaton. 3
4 3 Laplace Transform of the Stretchng factor The stretchng parameter s gven by so s + α + ωɛ ) 1 s α + ωɛ α + sɛ 0 α + sɛ 0 + ) 1 κ 1 α + sɛ 0 + So the nverse Laplace transform of 1/s s gven by ) L 1 1 δt) σ L 1 1 s κ 2 α + sɛ 0 + κ ) δt) σ κ 2 ɛ L s + α ɛ 0 + ɛ 0 δt) σ ) κ 2 ɛ e α ɛ + ɛ t 0 ut) 0 δt) + ζ t). wth γ ζ t) σ κ 2 ɛ e γ t ut) 0 α + σ ) 1ɛ0 α + σ ) ɛ 0 ɛ 0 Ths agrees wth Roden and Gedney s [4] result. 4 Implementaton by Convoluton At the end, we should convolve the functon ζ t) wth a dfferental operator. Let us refer to the tme varable t as n for the n th grd coordnate and 4
5 use a super ndex for t. Let us call the tme convoluton at some gven x, as ψ n x and at the grd pont x, n t). Then by the defnton of convoluton ψ n ζ x ) n n t 0 ) n t τ ζ x τ)dτ We now make use of the fact that we have a the data n a grd, and f a staggered grd method s used, the partal dervatve, s defned half a tme step between and m + 1) t, so ψ n n 1 m0 n 1 m+1) t ) n t τ ζ τ)dτ m+1) t ) n m+1/2) ζ τ)dτ m0 n 1 m0 Z m) ) n m+1/2), where Z m) m+1) t ɛ 0 κ 2 ɛ 0 κ 2 ζ τ)dτ m+1) t e γ t m+1) t e γ t ut) γ σ ɛ 0 κ 2 γ e mγ t e γ t 1 ) κ 2 α + e ) a e α + ɛ 0 α + ) [ ) ɛ 0 e α + t ɛ 0 1 ] wth a κ 2 α + [e α + ) t ɛ 0 1 ] 5 κ 2 α + b 1)
6 as n Roden and Gedney s [4]. Here as n Komattsch and Martn [1]. ) b e α + t ɛ The recurson formula We found ψ n wth so n 1 m0 Z m) ) n m+1/2) Z 0) ) n+1/2 + n 2 m0 Z m + 1) ) n m+1+1/2), Z m + 1) a e α + ) m+1) t ɛ 0 Z m)e α + ) t ɛ 0 Z m)b. ψ n Z 0) ) n+1/2 + n 2 m0 Z m + 1) ) n m+1+1/2) n 2 Z 0) ) n+1/2 + b Z m) ) n 1 m+1/2) m0 Z 0) ) n+1/2 + b ψ n 1 a ) n+1/2 + b ψ n 1 snce Z 0) a. Ths s equaton 26) n Komattsch and Martn [1]. Ths equaton, thanks to the recursve propertes of the exponental, provdes an effcent mplementaton of the boundary layer method. 4.2 A partcular mplementaton Komattsch and Martn [1] plemented the CPMBL, as follows: For the functon they used the symbol d and defned as d x ) d 0 x /L) N 6
7 where: x 1 x, x 2 y, x 3 z N 2 d 0 N + 1)v log R c 2L, v velocty R c L n d 100 thckness of absorbng layer Also they chose 1.0 The value of α s pcked lnear from the maxmum at the PML entrance wth a value of α max πf 0 wth f 0 the domnant frequency of the source wavelet) and 0 at the boundary. At the boundary they mposed Drechlet boundary condtons of v x v z 0. The varables: d, a, b, and α should be computed both at any nteger and nteger+half grd locaton accordng to the staggered grd recep. References [1] Komattsch M and R. Martn. An unsplt convolutonal perfectly matched layer mproved at grazng ncdence for sesmc wave equaton. Geophyscs, 725):SM155 SM167, [2] D. Pasalc and McGarry R. Convolutonal perfectly matched layer for sotropc and ansotropc acoustc wave equatons. SEG Techncal Program Expanded Abstracts, pages , [3] Martn R, D. Komattsch, and A. Ezzan. An unsplt convolutonal perfectly matched layer mproved at grazng ncdence for sesmc wave propagaton n poroelastc meda. Geophyscs, 734):T51 T61, [4] J. A. Roden and S. D. Gedney. Convoluton pml cpml): an effcent fdtd mplementaton of cfs-pml for arbtrary meda. Mcrowave and Optcal Technology Letters, 275): ,
Convolutional PML (CPML): An Efficient FDTD Implementation of the CFS-PML for Arbitrary Media
CPML: An Effcent FDTD Implementaton of the CFS-PML for Arbtrary Meda Not for Dstrbuton Convolutonal PML (CPML): An Effcent FDTD Implementaton of the CFS-PML for Arbtrary Meda J. Alan Roden Personal Systems
More information3: Central Limit Theorem, Systematic Errors
3: Central Lmt Theorem, Systematc Errors 1 Errors 1.1 Central Lmt Theorem Ths theorem s of prme mportance when measurng physcal quanttes because usually the mperfectons n the measurements are due to several
More informationII. Random Variables. Variable Types. Variables Map Outcomes to Numbers
II. Random Varables Random varables operate n much the same way as the outcomes or events n some arbtrary sample space the dstncton s that random varables are smply outcomes that are represented numercally.
More information15-451/651: Design & Analysis of Algorithms January 22, 2019 Lecture #3: Amortized Analysis last changed: January 18, 2019
5-45/65: Desgn & Analyss of Algorthms January, 09 Lecture #3: Amortzed Analyss last changed: January 8, 09 Introducton In ths lecture we dscuss a useful form of analyss, called amortzed analyss, for problems
More informationParallel Prefix addition
Marcelo Kryger Sudent ID 015629850 Parallel Prefx addton The parallel prefx adder presented next, performs the addton of two bnary numbers n tme of complexty O(log n) and lnear cost O(n). Lets notce the
More informationPrice and Quantity Competition Revisited. Abstract
rce and uantty Competton Revsted X. Henry Wang Unversty of Mssour - Columba Abstract By enlargng the parameter space orgnally consdered by Sngh and Vves (984 to allow for a wder range of cost asymmetry,
More information/ Computational Genomics. Normalization
0-80 /02-70 Computatonal Genomcs Normalzaton Gene Expresson Analyss Model Computatonal nformaton fuson Bologcal regulatory networks Pattern Recognton Data Analyss clusterng, classfcaton normalzaton, mss.
More informationRandom Variables. b 2.
Random Varables Generally the object of an nvestgators nterest s not necessarly the acton n the sample space but rather some functon of t. Techncally a real valued functon or mappng whose doman s the sample
More informationCS 286r: Matching and Market Design Lecture 2 Combinatorial Markets, Walrasian Equilibrium, Tâtonnement
CS 286r: Matchng and Market Desgn Lecture 2 Combnatoral Markets, Walrasan Equlbrum, Tâtonnement Matchng and Money Recall: Last tme we descrbed the Hungaran Method for computng a maxmumweght bpartte matchng.
More informationAppendix - Normally Distributed Admissible Choices are Optimal
Appendx - Normally Dstrbuted Admssble Choces are Optmal James N. Bodurtha, Jr. McDonough School of Busness Georgetown Unversty and Q Shen Stafford Partners Aprl 994 latest revson September 00 Abstract
More informationMgtOp 215 Chapter 13 Dr. Ahn
MgtOp 5 Chapter 3 Dr Ahn Consder two random varables X and Y wth,,, In order to study the relatonshp between the two random varables, we need a numercal measure that descrbes the relatonshp The covarance
More informationEvaluating Performance
5 Chapter Evaluatng Performance In Ths Chapter Dollar-Weghted Rate of Return Tme-Weghted Rate of Return Income Rate of Return Prncpal Rate of Return Daly Returns MPT Statstcs 5- Measurng Rates of Return
More informationMode is the value which occurs most frequency. The mode may not exist, and even if it does, it may not be unique.
1.7.4 Mode Mode s the value whch occurs most frequency. The mode may not exst, and even f t does, t may not be unque. For ungrouped data, we smply count the largest frequency of the gven value. If all
More informationAppendix for Solving Asset Pricing Models when the Price-Dividend Function is Analytic
Appendx for Solvng Asset Prcng Models when the Prce-Dvdend Functon s Analytc Ovdu L. Caln Yu Chen Thomas F. Cosmano and Alex A. Hmonas January 3, 5 Ths appendx provdes proofs of some results stated n our
More informationTests for Two Correlations
PASS Sample Sze Software Chapter 805 Tests for Two Correlatons Introducton The correlaton coeffcent (or correlaton), ρ, s a popular parameter for descrbng the strength of the assocaton between two varables.
More informationEDC Introduction
.0 Introducton EDC3 In the last set of notes (EDC), we saw how to use penalty factors n solvng the EDC problem wth losses. In ths set of notes, we want to address two closely related ssues. What are, exactly,
More information2) In the medium-run/long-run, a decrease in the budget deficit will produce:
4.02 Quz 2 Solutons Fall 2004 Multple-Choce Questons ) Consder the wage-settng and prce-settng equatons we studed n class. Suppose the markup, µ, equals 0.25, and F(u,z) = -u. What s the natural rate of
More informationCOMPARISON OF THE ANALYTICAL AND NUMERICAL SOLUTION OF A ONE-DIMENSIONAL NON-STATIONARY COOLING PROBLEM. László Könözsy 1, Mátyás Benke 2
COMPARISON OF THE ANALYTICAL AND NUMERICAL SOLUTION OF A ONE-DIMENSIONAL NON-STATIONARY COOLING PROBLEM László Könözsy 1, Mátyás Benke Ph.D. Student 1, Unversty Student Unversty of Mskolc, Department of
More informationNotes on experimental uncertainties and their propagation
Ed Eyler 003 otes on epermental uncertantes and ther propagaton These notes are not ntended as a complete set of lecture notes, but nstead as an enumeraton of some of the key statstcal deas needed to obtan
More informationUnderstanding Annuities. Some Algebraic Terminology.
Understandng Annutes Ma 162 Sprng 2010 Ma 162 Sprng 2010 March 22, 2010 Some Algebrac Termnology We recall some terms and calculatons from elementary algebra A fnte sequence of numbers s a functon of natural
More information2.1 Rademacher Calculus... 3
COS 598E: Unsupervsed Learnng Week 2 Lecturer: Elad Hazan Scrbe: Kran Vodrahall Contents 1 Introducton 1 2 Non-generatve pproach 1 2.1 Rademacher Calculus............................... 3 3 Spectral utoencoders
More informationMultifactor Term Structure Models
1 Multfactor Term Structure Models A. Lmtatons of One-Factor Models 1. Returns on bonds of all maturtes are perfectly correlated. 2. Term structure (and prces of every other dervatves) are unquely determned
More informationTime Domain Decomposition for European Options in Financial Modelling
Contemporary Mathematcs Volume 218, 1998 B 0-8218-0988-1-03047-8 Tme Doman Decomposton for European Optons n Fnancal Modellng Dane Crann, Alan J. Daves, Cho-Hong La, and Swee H. Leong 1. Introducton Fnance
More informationy\ 1 Target E-2 Extra Practice r i r Date: Name: 1. a) What is the approximate value of d when t = 3? Explain the method you used.
Target E-2 Extra Practce. a) What s the approxmate value of d when t = 3? Explan the method ou used. b) What s the approxmate value of t when d = 300? 4 c o d> taa - Ov 0 - ~ r H - \ V/ 00-06- - -
More informationConsumption Based Asset Pricing
Consumpton Based Asset Prcng Mchael Bar Aprl 25, 208 Contents Introducton 2 Model 2. Prcng rsk-free asset............................... 3 2.2 Prcng rsky assets................................ 4 2.3 Bubbles......................................
More informationElements of Economic Analysis II Lecture VI: Industry Supply
Elements of Economc Analyss II Lecture VI: Industry Supply Ka Hao Yang 10/12/2017 In the prevous lecture, we analyzed the frm s supply decson usng a set of smple graphcal analyses. In fact, the dscusson
More informationEconomic Design of Short-Run CSP-1 Plan Under Linear Inspection Cost
Tamkang Journal of Scence and Engneerng, Vol. 9, No 1, pp. 19 23 (2006) 19 Economc Desgn of Short-Run CSP-1 Plan Under Lnear Inspecton Cost Chung-Ho Chen 1 * and Chao-Yu Chou 2 1 Department of Industral
More informationLikelihood Fits. Craig Blocker Brandeis August 23, 2004
Lkelhood Fts Crag Blocker Brandes August 23, 2004 Outlne I. What s the queston? II. Lkelhood Bascs III. Mathematcal Propertes IV. Uncertantes on Parameters V. Mscellaneous VI. Goodness of Ft VII. Comparson
More informationECE 586GT: Problem Set 2: Problems and Solutions Uniqueness of Nash equilibria, zero sum games, evolutionary dynamics
Unversty of Illnos Fall 08 ECE 586GT: Problem Set : Problems and Solutons Unqueness of Nash equlbra, zero sum games, evolutonary dynamcs Due: Tuesday, Sept. 5, at begnnng of class Readng: Course notes,
More informationMonte Carlo Rendering
Last Tme? Monte Carlo Renderng Monte-Carlo Integraton Probabltes and Varance Analyss of Monte-Carlo Integraton Monte-Carlo n Graphcs Stratfed Samplng Importance Samplng Advanced Monte-Carlo Renderng Monte-Carlo
More informationTCOM501 Networking: Theory & Fundamentals Final Examination Professor Yannis A. Korilis April 26, 2002
TO5 Networng: Theory & undamentals nal xamnaton Professor Yanns. orls prl, Problem [ ponts]: onsder a rng networ wth nodes,,,. In ths networ, a customer that completes servce at node exts the networ wth
More informationSupplementary material for Non-conjugate Variational Message Passing for Multinomial and Binary Regression
Supplementary materal for Non-conjugate Varatonal Message Passng for Multnomal and Bnary Regresson October 9, 011 1 Alternatve dervaton We wll focus on a partcular factor f a and varable x, wth the am
More information3/3/2014. CDS M Phil Econometrics. Vijayamohanan Pillai N. Truncated standard normal distribution for a = 0.5, 0, and 0.5. CDS Mphil Econometrics
Lmted Dependent Varable Models: Tobt an Plla N 1 CDS Mphl Econometrcs Introducton Lmted Dependent Varable Models: Truncaton and Censorng Maddala, G. 1983. Lmted Dependent and Qualtatve Varables n Econometrcs.
More informationME 310 Numerical Methods. Differentiation
M 0 Numercal Metods fferentaton Tese presentatons are prepared by r. Cuneyt Sert Mecancal ngneerng epartment Mddle ast Tecncal Unversty Ankara, Turkey csert@metu.edu.tr Tey can not be used wtout te permsson
More informationCOS 511: Theoretical Machine Learning. Lecturer: Rob Schapire Lecture #21 Scribe: Lawrence Diao April 23, 2013
COS 511: Theoretcal Machne Learnng Lecturer: Rob Schapre Lecture #21 Scrbe: Lawrence Dao Aprl 23, 2013 1 On-Lne Log Loss To recap the end of the last lecture, we have the followng on-lne problem wth N
More informationQuiz on Deterministic part of course October 22, 2002
Engneerng ystems Analyss for Desgn Quz on Determnstc part of course October 22, 2002 Ths s a closed book exercse. You may use calculators Grade Tables There are 90 ponts possble for the regular test, or
More informationComparative analysis of CDO pricing models
Comparatve analyss of CDO prcng models ICBI Rsk Management 2005 Geneva 8 December 2005 Jean-Paul Laurent ISFA, Unversty of Lyon, Scentfc Consultant BNP Parbas laurent.jeanpaul@free.fr, http://laurent.jeanpaul.free.fr
More informationMicroeconomics: BSc Year One Extending Choice Theory
mcroeconomcs notes from http://www.economc-truth.co.uk by Tm Mller Mcroeconomcs: BSc Year One Extendng Choce Theory Consumers, obvously, mostly have a choce of more than two goods; and to fnd the favourable
More informationLecture 7. We now use Brouwer s fixed point theorem to prove Nash s theorem.
Topcs on the Border of Economcs and Computaton December 11, 2005 Lecturer: Noam Nsan Lecture 7 Scrbe: Yoram Bachrach 1 Nash s Theorem We begn by provng Nash s Theorem about the exstance of a mxed strategy
More informationASPECTS OF PRICING IRREGULAR SWAPTIONS WITH QUANTLIB Calibration and Pricing with the LGM Model
ASPECTS OF PRICING IRREGULAR SWAPTIONS WITH QUANTLIB Calbraton and Prcng wth the LGM Model HSH NORDBANK Dr. Werner Kürznger Düsseldorf, November 30th, 2017 HSH-NORDBANK.DE Dsclamer The content of ths presentaton
More informationEconomics 1410 Fall Section 7 Notes 1. Define the tax in a flexible way using T (z), where z is the income reported by the agent.
Economcs 1410 Fall 2017 Harvard Unversty Yaan Al-Karableh Secton 7 Notes 1 I. The ncome taxaton problem Defne the tax n a flexble way usng T (), where s the ncome reported by the agent. Retenton functon:
More informationFast Laplacian Solvers by Sparsification
Spectral Graph Theory Lecture 19 Fast Laplacan Solvers by Sparsfcaton Danel A. Spelman November 9, 2015 Dsclamer These notes are not necessarly an accurate representaton of what happened n class. The notes
More informationApplications of Myerson s Lemma
Applcatons of Myerson s Lemma Professor Greenwald 28-2-7 We apply Myerson s lemma to solve the sngle-good aucton, and the generalzaton n whch there are k dentcal copes of the good. Our objectve s welfare
More informationMULTIPLE CURVE CONSTRUCTION
MULTIPLE CURVE CONSTRUCTION RICHARD WHITE 1. Introducton In the post-credt-crunch world, swaps are generally collateralzed under a ISDA Master Agreement Andersen and Pterbarg p266, wth collateral rates
More informationHomework 9: due Monday, 27 October, 2008
PROBLEM ONE Homework 9: due Monday, 7 October, 008. (Exercses from the book, 6 th edton, 6.6, -3.) Determne the number of dstnct orderngs of the letters gven: (a) GUIDE (b) SCHOOL (c) SALESPERSONS. (Exercses
More informationGames and Decisions. Part I: Basic Theorems. Contents. 1 Introduction. Jane Yuxin Wang. 1 Introduction 1. 2 Two-player Games 2
Games and Decsons Part I: Basc Theorems Jane Yuxn Wang Contents 1 Introducton 1 2 Two-player Games 2 2.1 Zero-sum Games................................ 3 2.1.1 Pure Strateges.............................
More informationCorrelations and Copulas
Correlatons and Copulas Chapter 9 Rsk Management and Fnancal Insttutons, Chapter 6, Copyrght John C. Hull 2006 6. Coeffcent of Correlaton The coeffcent of correlaton between two varables V and V 2 s defned
More informationFoundations of Machine Learning II TP1: Entropy
Foundatons of Machne Learnng II TP1: Entropy Gullaume Charpat (Teacher) & Gaétan Marceau Caron (Scrbe) Problem 1 (Gbbs nequalty). Let p and q two probablty measures over a fnte alphabet X. Prove that KL(p
More informationComputational Finance
Department of Mathematcs at Unversty of Calforna, San Dego Computatonal Fnance Dfferental Equaton Technques [Lectures 8-10] Mchael Holst February 27, 2017 Contents 1 Modelng Fnancal Optons wth the Black-Scholes
More informationHewlett Packard 10BII Calculator
Hewlett Packard 0BII Calculator Keystrokes for the HP 0BII are shown n the tet. However, takng a mnute to revew the Quk Start secton, below, wll be very helpful n gettng started wth your calculator. Note:
More informationTaxation and Externalities. - Much recent discussion of policy towards externalities, e.g., global warming debate/kyoto
Taxaton and Externaltes - Much recent dscusson of polcy towards externaltes, e.g., global warmng debate/kyoto - Increasng share of tax revenue from envronmental taxaton 6 percent n OECD - Envronmental
More informationSingle-Item Auctions. CS 234r: Markets for Networks and Crowds Lecture 4 Auctions, Mechanisms, and Welfare Maximization
CS 234r: Markets for Networks and Crowds Lecture 4 Auctons, Mechansms, and Welfare Maxmzaton Sngle-Item Auctons Suppose we have one or more tems to sell and a pool of potental buyers. How should we decde
More information>1 indicates country i has a comparative advantage in production of j; the greater the index, the stronger the advantage. RCA 1 ij
69 APPENDIX 1 RCA Indces In the followng we present some maor RCA ndces reported n the lterature. For addtonal varants and other RCA ndces, Memedovc (1994) and Vollrath (1991) provde more thorough revews.
More informationParsing beyond context-free grammar: Tree Adjoining Grammar Parsing I
Parsng beyond context-free grammar: Tree donng Grammar Parsng I Laura Kallmeyer, Wolfgang Maer ommersemester 2009 duncton and substtuton (1) Tree donng Grammars (TG) Josh et al. (1975), Josh & chabes (1997):
More informationProblem Set 6 Finance 1,
Carnege Mellon Unversty Graduate School of Industral Admnstraton Chrs Telmer Wnter 2006 Problem Set 6 Fnance, 47-720. (representatve agent constructon) Consder the followng two-perod, two-agent economy.
More informationOPERATIONS RESEARCH. Game Theory
OPERATIONS RESEARCH Chapter 2 Game Theory Prof. Bbhas C. Gr Department of Mathematcs Jadavpur Unversty Kolkata, Inda Emal: bcgr.umath@gmal.com 1.0 Introducton Game theory was developed for decson makng
More informationCHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS
CHAPTER 9 FUNCTIONAL FORMS OF REGRESSION MODELS QUESTIONS 9.1. (a) In a log-log model the dependent and all explanatory varables are n the logarthmc form. (b) In the log-ln model the dependent varable
More informationNote on Cubic Spline Valuation Methodology
Note on Cubc Splne Valuaton Methodology Regd. Offce: The Internatonal, 2 nd Floor THE CUBIC SPLINE METHODOLOGY A model for yeld curve takes traded yelds for avalable tenors as nput and generates the curve
More informationData Mining Linear and Logistic Regression
07/02/207 Data Mnng Lnear and Logstc Regresson Mchael L of 26 Regresson In statstcal modellng, regresson analyss s a statstcal process for estmatng the relatonshps among varables. Regresson models are
More information- contrast so-called first-best outcome of Lindahl equilibrium with case of private provision through voluntary contributions of households
Prvate Provson - contrast so-called frst-best outcome of Lndahl equlbrum wth case of prvate provson through voluntary contrbutons of households - need to make an assumpton about how each household expects
More informationAn Application of Alternative Weighting Matrix Collapsing Approaches for Improving Sample Estimates
Secton on Survey Research Methods An Applcaton of Alternatve Weghtng Matrx Collapsng Approaches for Improvng Sample Estmates Lnda Tompkns 1, Jay J. Km 2 1 Centers for Dsease Control and Preventon, atonal
More informationMeasures of Spread IQR and Deviation. For exam X, calculate the mean, median and mode. For exam Y, calculate the mean, median and mode.
Part 4 Measures of Spread IQR and Devaton In Part we learned how the three measures of center offer dfferent ways of provdng us wth a sngle representatve value for a data set. However, consder the followng
More informationChapter 3 Descriptive Statistics: Numerical Measures Part B
Sldes Prepared by JOHN S. LOUCKS St. Edward s Unversty Slde 1 Chapter 3 Descrptve Statstcs: Numercal Measures Part B Measures of Dstrbuton Shape, Relatve Locaton, and Detectng Outlers Eploratory Data Analyss
More informationarxiv: v2 [math.co] 6 Apr 2016
On the number of equvalence classes of nvertble Boolean functons under acton of permutaton of varables on doman and range arxv:1603.04386v2 [math.co] 6 Apr 2016 Marko Carć and Modrag Žvkovć Abstract. Let
More informationCliquet Options and Volatility Models
Clquet Optons and olatlty Models Paul Wlmott paul@wlmott.com 1 Introducton Clquet optons are at present the heght of fashon n the world of equty dervatves. These contracts, llustrated by the term sheet
More informationCyclic Scheduling in a Job shop with Multiple Assembly Firms
Proceedngs of the 0 Internatonal Conference on Industral Engneerng and Operatons Management Kuala Lumpur, Malaysa, January 4, 0 Cyclc Schedulng n a Job shop wth Multple Assembly Frms Tetsuya Kana and Koch
More informationMacroeconomic Theory and Policy
ECO 209 Macroeconomc Theory and Polcy Lecture 7: The Open Economy wth Fxed Exchange Rates Gustavo Indart Slde 1 Open Economy under Fxed Exchange Rates Let s consder an open economy wth no captal moblty
More informationA DUAL EXTERIOR POINT SIMPLEX TYPE ALGORITHM FOR THE MINIMUM COST NETWORK FLOW PROBLEM
Yugoslav Journal of Operatons Research Vol 19 (2009), Number 1, 157-170 DOI:10.2298/YUJOR0901157G A DUAL EXTERIOR POINT SIMPLEX TYPE ALGORITHM FOR THE MINIMUM COST NETWORK FLOW PROBLEM George GERANIS Konstantnos
More information1 A Primer on Linear Models. 2 Chapter 1 corrections. 3 Chapter 2 corrections. 4 Chapter 3 corrections. 1.1 Corrections 23 May 2015
A Prmer on Lnear Models. Correctons 23 May 205 2 Chapter correctons Fx page: 9 lne -7 denomnator ( ρ 2 ) s mssng, should read Cov(e t, e s ) = V ar(a t )/( ρ 2 ) Fx page: lne -5 sgn needs changng: cos(a
More informationDependent jump processes with coupled Lévy measures
Dependent jump processes wth coupled Lévy measures Naoufel El-Bachr ICMA Centre, Unversty of Readng May 6, 2008 ICMA Centre Dscusson Papers n Fnance DP2008-3 Copyrght 2008 El-Bachr. All rghts reserved.
More informationCollective Motion from Consensus with Cartesian Coordinate Coupling - Part II: Double-integrator Dynamics
Proceedngs of the 47th IEEE Conference on Decson Control Cancun Mexco Dec. 9-8 TuB. Collectve Moton from Consensus wth Cartesan Coordnate Couplng - Part II: Double-ntegrator Dynamcs We Ren Abstract Ths
More informationFinite Math - Fall Section Future Value of an Annuity; Sinking Funds
Fnte Math - Fall 2016 Lecture Notes - 9/19/2016 Secton 3.3 - Future Value of an Annuty; Snkng Funds Snkng Funds. We can turn the annutes pcture around and ask how much we would need to depost nto an account
More informationChapter 5 Student Lecture Notes 5-1
Chapter 5 Student Lecture Notes 5-1 Basc Busness Statstcs (9 th Edton) Chapter 5 Some Important Dscrete Probablty Dstrbutons 004 Prentce-Hall, Inc. Chap 5-1 Chapter Topcs The Probablty Dstrbuton of a Dscrete
More informationarxiv: v1 [math-ph] 19 Oct 2007
Monomal ntegrals on the classcal groups T. Gorn and G. V. López Departamento de Físca, Unversdad de Guadalajara Blvd. Marcelno García Barragan y Calzada Olímpca 44840 Guadalajara, Jalsco, Méxco February,
More informationCh Rival Pure private goods (most retail goods) Non-Rival Impure public goods (internet service)
h 7 1 Publc Goods o Rval goods: a good s rval f ts consumpton by one person precludes ts consumpton by another o Excludable goods: a good s excludable f you can reasonably prevent a person from consumng
More informationSpatial Variations in Covariates on Marriage and Marital Fertility: Geographically Weighted Regression Analyses in Japan
Spatal Varatons n Covarates on Marrage and Martal Fertlty: Geographcally Weghted Regresson Analyses n Japan Kenj Kamata (Natonal Insttute of Populaton and Socal Securty Research) Abstract (134) To understand
More informationAnswers to exercises in Macroeconomics by Nils Gottfries 2013
. a) C C b C C s the ntercept o the consumpton uncton, how much consumpton wll be at zero ncome. We can thnk that, at zero ncome, the typcal consumer would consume out o hs assets. The slope b s the margnal
More informationThe Direct Control and Penalty Methods for American Put Options
The Drect Control and Penalty Methods for Amercan Put Optons by Ama Peprah Asare A thess presented to the Unversty of Waterloo n fulfllment of the thess requrement for the degree of Master of Mathematcs
More informationUnderstanding Predictability (JPE, 2004)
Understandng Predctablty (JPE, 2004) Lor Menzly, Tano Santos, and Petro Verones Presented by Peter Gross NYU October 27, 2009 Presented by Peter Gross (NYU) Understandng Predctablty October 27, 2009 1
More informationECO 209Y MACROECONOMIC THEORY AND POLICY LECTURE 8: THE OPEN ECONOMY WITH FIXED EXCHANGE RATES
ECO 209 MACROECONOMIC THEOR AND POLIC LECTURE 8: THE OPEN ECONOM WITH FIXED EXCHANGE RATES Gustavo Indart Slde 1 OPEN ECONOM UNDER FIXED EXCHANGE RATES Let s consder an open economy wth no captal moblty
More informationNotes on Debye-Hückel Theory
GEOS 0300/30300 Autumn 004 Thermodynamcs and Phase Change MS Ghorso Notes on Debye-Hückel Theory We seek: µ µ o + RT ln X + µ ex where µ ex RT ln λ or the devaton from deal behavor due to long range nteracton
More informationTHE IMPORTANCE OF THE NUMBER OF DIFFERENT AGENTS IN A HETEROGENEOUS ASSET-PRICING MODEL WOUTER J. DEN HAAN
THE IMPORTANCE OF THE NUMBER OF DIFFERENT AGENTS IN A HETEROGENEOUS ASSET-PRICING MODEL WOUTER J. DEN HAAN Department of Economcs, Unversty of Calforna at San Dego and Natonal Bureau of Economc Research
More informationTests for Two Ordered Categorical Variables
Chapter 253 Tests for Two Ordered Categorcal Varables Introducton Ths module computes power and sample sze for tests of ordered categorcal data such as Lkert scale data. Assumng proportonal odds, such
More informationCentre for International Capital Markets
Centre for Internatonal Captal Markets Dscusson Papers ISSN 1749-3412 Valung Amercan Style Dervatves by Least Squares Methods Maro Cerrato No 2007-13 Valung Amercan Style Dervatves by Least Squares Methods
More informationFinancial mathematics
Fnancal mathematcs Jean-Luc Bouchot jean-luc.bouchot@drexel.edu February 19, 2013 Warnng Ths s a work n progress. I can not ensure t to be mstake free at the moment. It s also lackng some nformaton. But
More informationPricing Variance Swaps with Cash Dividends
Prcng Varance Swaps wth Cash Dvdends Tmothy Klassen Abstract We derve a smple formula for the prce of a varance swap when the underlyng has cash dvdends. 1 Introducton The last years have seen renewed
More informationEfficient Project Portfolio as a Tool for Enterprise Risk Management
Effcent Proect Portfolo as a Tool for Enterprse Rsk Management Valentn O. Nkonov Ural State Techncal Unversty Growth Traectory Consultng Company Enterprse Rsk Management Symposum Socety of Actuares Chcago,
More informationAn annuity is a series of payments made at equal intervals. There are many practical examples of financial transactions involving annuities, such as
2 Annutes An annuty s a seres of payments made at equal ntervals. There are many practcal examples of fnancal transactons nvolvng annutes, such as a car loan beng repad wth equal monthly nstallments a
More informationGlobal sensitivity analysis of credit risk portfolios
Global senstvty analyss of credt rsk portfolos D. Baur, J. Carbon & F. Campolongo European Commsson, Jont Research Centre, Italy Abstract Ths paper proposes the use of global senstvty analyss to evaluate
More informationNumerical Optimisation Applied to Monte Carlo Algorithms for Finance. Phillip Luong
Numercal Optmsaton Appled to Monte Carlo Algorthms for Fnance Phllp Luong Supervsed by Professor Hans De Sterck, Professor Gregore Loeper, and Dr Ivan Guo Monash Unversty Vacaton Research Scholarshps are
More informationCreating a zero coupon curve by bootstrapping with cubic splines.
MMA 708 Analytcal Fnance II Creatng a zero coupon curve by bootstrappng wth cubc splnes. erg Gryshkevych Professor: Jan R. M. Röman 0.2.200 Dvson of Appled Mathematcs chool of Educaton, Culture and Communcaton
More informationMacroeconomic Theory and Policy
ECO 209 Macroeconomc Theory and Polcy Lecture 7: The Open Economy wth Fxed Exchange Rates Gustavo Indart Slde 1 Open Economy under Fxed Exchange Rates Let s consder an open economy wth no captal moblty
More informationA Single-Product Inventory Model for Multiple Demand Classes 1
A Sngle-Product Inventory Model for Multple Demand Classes Hasan Arslan, 2 Stephen C. Graves, 3 and Thomas Roemer 4 March 5, 2005 Abstract We consder a sngle-product nventory system that serves multple
More informationFinal Exam. 7. (10 points) Please state whether each of the following statements is true or false. No explanation needed.
Fnal Exam Fall 4 Econ 8-67 Closed Book. Formula Sheet Provded. Calculators OK. Tme Allowed: hours Please wrte your answers on the page below each queston. (5 ponts) Assume that the rsk-free nterest rate
More informationIntroduction to PGMs: Discrete Variables. Sargur Srihari
Introducton to : Dscrete Varables Sargur srhar@cedar.buffalo.edu Topcs. What are graphcal models (or ) 2. Use of Engneerng and AI 3. Drectonalty n graphs 4. Bayesan Networks 5. Generatve Models and Samplng
More informationDiscrete Dynamic Shortest Path Problems in Transportation Applications
17 Paper No. 98-115 TRANSPORTATION RESEARCH RECORD 1645 Dscrete Dynamc Shortest Path Problems n Transportaton Applcatons Complexty and Algorthms wth Optmal Run Tme ISMAIL CHABINI A soluton s provded for
More informationMaximum Likelihood Estimation of Isotonic Normal Means with Unknown Variances*
Journal of Multvarate Analyss 64, 183195 (1998) Artcle No. MV971717 Maxmum Lelhood Estmaton of Isotonc Normal Means wth Unnown Varances* Nng-Zhong Sh and Hua Jang Northeast Normal Unversty, Changchun,Chna
More informationMacroeconomic equilibrium in the short run: the Money market
Macroeconomc equlbrum n the short run: the Money market 2013 1. The bg pcture Overvew Prevous lecture How can we explan short run fluctuatons n GDP? Key assumpton: stcky prces Equlbrum of the goods market
More informationIntroduction to game theory
Introducton to game theory Lectures n game theory ECON5210, Sprng 2009, Part 1 17.12.2008 G.B. Ashem, ECON5210-1 1 Overvew over lectures 1. Introducton to game theory 2. Modelng nteractve knowledge; equlbrum
More informationOn the Moments of the Traces of Unitary and Orthogonal Random Matrices
Proceedngs of Insttute of Mathematcs of NAS of Ukrane 2004 Vol. 50 Part 3 1207 1213 On the Moments of the Traces of Untary and Orthogonal Random Matrces Vladmr VASILCHU B. Verkn Insttute for Low Temperature
More information