EE6604 Personal & Mobile Communications. Week 9. Co-Channel Interference

Size: px
Start display at page:

Download "EE6604 Personal & Mobile Communications. Week 9. Co-Channel Interference"

Transcription

1 EE6604 Personal & Mobile Communications Week 9 Co-Channel Interference 1

2 Co-channel interference on the forward channel d 1 d 6 d 2 mobile subscriber d 0 d 5 d 3 d 4 The mobile station is being served by the center base station. 2

3 At a particular location, let d = (d 0, d 1,, d NI ) be the vector of distances between a mobile station and the serving base station BS 0 and N I co-channel base stations BS k,k = 1,...,N I. The received signal power power at distance d, Ω p (dbm) (d), is a Gaussian random variable that depends on the distance d through the path loss model, i.e., µ Ωp (dbm) (d) = E[Ω p (dbm) (d)] = µ Ωp (dbm) (d o ) 10βlog 10 (d/d o ) Experiments have verified that co-channel interferers add noncoherently (power addition) rather than coherently (amplitude addition). The C/I a function of the vector d is or in decibel units Λ(d) = Ω p(d 0 ) N I k=1ω p (d k ) Λ(d) (db) = Ω p (dbm) (d 0 ) 10log 10 The outage probability given vector d is N I k=1 O(d) = P r ( Λ(d)(dB) < Λ th(db) ) Ω p (d k ) Although the Ω p (d k ) are log-normal random variables, the sum N I k=1ω p (d k ) is not a log-normal random variable. 3

4 Multiple Log-normal Interferers Consider the sum of N I log-normal random variables I = N I k=1 Ω k = N I 10 Ω k(dbm) /10 k=1 where the Ω k (dbm) are independent Gaussian random variables with mean µ Ωk (dbm) and variance σ 2 Ω k. The sum I is commonly approximated by another log-normal random variable with appropriately chosen parameters, i.e., I = N I k=1 10 Ω k(dbm) /10 10 Z (dbm) /10 = Î where Z (dbm) is a Gaussian random variable with mean µ Z (dbm) and variance σ 2 Z. The task is to find µ Z (dbm) and σ 2 Z. 4

5 Fenton-Wilkinson Method The mean µ Z (dbm) and variance σ 2 Z of Z (dbm) are obtained by matching the first two moments of I and Î. Switching from base 10 to base e: Ω k = 10 Ω k (dbm) /10 = e ξω k (dbm) = eˆω k where ˆΩ k = ξω k (dbm) and ξ = (ln10)/10 = Note that µˆωk = ξµ Ωk (dbm) σ 2ˆΩk = ξ 2 σ 2 Ω k The nth moment of the log-normal random variable Ω k can be obtained from the moment generating function of the Gaussian random variable ˆΩ k as E[Ω n k] = E[e nˆω k ] = e nµˆωk +(1/2)n 2 σ 2ˆΩ Here we have assumed identical shadow variances, σ 2ˆΩk = σ 2ˆΩ, which is a reasonable assumption. 5

6 Suppose that ˆΩ 1,..., ˆΩNI are independent with means µˆω1,..., µˆωni and identical variances σ 2ˆΩ. The appropriate moments of the log-normal approximation are obtained by equating the means on both sides of where Ẑ = ξz (dbm). This gives µ I = E[I] = N I N I k=1 k=1 e µˆωk Also equate the variances on both sides of E[eˆΩ k ] E[eẐ] = E[Î] = µ Î e (1/2)σ2ˆΩ = e µ Ẑ +(1/2)σ2 Ẑ (1) σ 2 I = E[I 2 ] µ 2 I E[Î2 ] µ 2 Î = σ2 Î This gives N I k=1 e 2µˆΩk e σ2ˆω(e σ2ˆω 1) = e 2µ Ẑe σ2 Ẑ(e σ2 Ẑ 1) (2) 6

7 To obtain µẑ and σ 2 Ẑ 1. Square Eq. (1) and divide by Eq. (2) to obtain σ 2 Ẑ. 2. Obtain µẑ from Eq. (1) The above procedure yields σ 2 Ẑ = ln µẑ = σ2ˆω σ 2 Ẑ 2 (eσ2ˆω 1) ( NI +ln N I k=1e 2µˆΩk k=1e µˆωk ) 2 +1 N I k=1 e µˆωk Given the means µˆω1,..., µˆωni and variance σ 2ˆΩ, µẑ and σ 2 Ẑ are easily obtained. Finally, we convert back to base 10 by scaling, such that µ Z (dbm) = ξ 1 µẑ σz 2 = ξ 2 σ 2 Ẑ where ξ =

8 Fenton s method breaks down in the prediction of the first and second moments for σ Ω > 4 db. Schwartz and Yeh s method yields the exact first and second moments. However, Fenton s method accurately predicts the tails of the complementary distribution function cdfc F c I (x) = P r(i x) and the cdf F I (x) = 1 F c I (x) = P r (I < x). We are interested in the accuracy of the approximations FI(x) c Q lnx µẑ F I (x) 1 Q when x is large and small, respectively. σẑ lnx µẑ σẑ The cdfc is more important than the cdf for outage calculations and predictions, since outages typically occur when the interference is large. 8

9 Wilkinson Schwartz and Yeh Farley Simulation P(I<x) N I =2 N I = x Comparison of the cdf for the sum of two and six log-normal random variables with various approximations; σ Ω = 6 db. 9

10 P(I>x) Wilkinson Schwartz and Yeh Farley Simulation x Comparison of the cdfc for the sum of two log-normal random variables with various approximations; σ Ω = 6 db. 10

11 P(I>x) Wilkinson Schwartz and Yeh Farley Simulation x Comparison of the cdfc for the sum of six log-normal random variables with various approximations; σ Ω = 6 db. 11

12 P(I>x) Wilkinson Schwartz and Yeh Farley Simulation x Comparison of the cdfc for the sum of six log-normal random variables with various approximations; σ Ω = 12 db. 12

13 Outage with Multiple Interferers 1. First obtain the mean and variance µ Z = µẑ/ξ σ 2 Z = σ 2 Ẑ /ξ2 ξ = Treat the average CIR as Gaussian distributed with mean and variance µ Λ(d) = µ Ω(d0 ) µ Z (dbm) σ 2 Λ(d) = σ 2 Ω +σ2 Z. 3. Compute the outage for a given location, described by d O(d) = Q µ Ω(d0 ) µ Z Λ th(db) σ 2 Ω +σ 2 Z 4. Average over all locations d by Monte Carlo integration O = R N O(d)p d (d)dd 13

14 Single Co-channel Interferer For a single co-channel interferer where p Λ(d)(dB) (x) = The outage for a given d is 1 exp 4πσΩ (x µ Λ(d) (db) ) 2 4σ 2 Ω µ Λ(d)(dB) = µ Ω(d0 ) (db) µ Ω(d1 ) (db) O(d) = Pr(Λ(d) (db) < Λ th(db) ) = Λ th(db) = Q 1 4πσΩ exp µ Λ(d)(dB) Λ th(db) 2σΩ (x µ Λ(d) (db) ) 2 4σ 2 Ω dx 14

15 R D-R serving MS co-channel BS BS Worst case interference from a single co-channel base-station. In this case d = (R,D R). The worst case outage due to a single co-channel interferer is O(R) = Q µ Ω(R)(dB) µ Ω(D R)(dB) Λ th (db) 2σΩ 15

16 Using a simple inverse-β path loss characteristic µ Ω(dB) = Ω (db) (d o ) 10βlog 10 (d/d o ) gives O(R) = Q 10log 10 ( D R 1) β Λth (db) 2σΩ The minimum CIR margin on the cell fringe is M Λ = 10log 10 D R 1 For an ideal hexagonal layout D R = 3N, so that N = β Λ th (db) M Λ +Λ 2 th (db) 10β +1 A small cluster size is achieved by making the margin M Λ and receiver threshold Λ th small. 16

17 Rician/Multiple Rayleigh Interferers Sometimes propagation conditions exist such that the received signals experience fading, but not shadowing. In this section, we calculate the outage probability for the case of fading only. The received signal may consist of a direct line of sight (LoS) component, or perhaps a specular component, accompanied by a diffuse component. The envelope of the received desired signal experiences Ricean fading. The interfering signals are often assumed to be Rayleigh faded, because a direct LoS is unlikely to exist due to the larger physical distances between the co-channel interferers and the receiver. Let the instantaneous power in the desired signal and the N I interfering signals be denoted by s 0 and s k, k = 1,, N I, respectively. Note that s i = α 2 i, where α2 i is the squared-envelope. The carrier-to-interference ratio is defined as λ = s 0 / N I k=1s k, and for a specified receiver threshold λ th, the outage probability is O I = P(λ < λ th ). 17

18 Single Interferer For the case of a single interferer, the outage probability reduces to the simple closed form O I = λ th exp λ th +A KA 1 1 λ th +A 1, where K is the Rice factor of the desired signal, A 1 = Ω 0 /(K + 1)Ω 1, and Ω k = E[s k ]. If the desired signal is Rayleigh faded, then the outage probability can be obtained by setting K = 0. 18

19 Multiple Interferers For the case of multiple interferers, each with mean power Ω k, the outage probability has the closed form O I = 1 N I k=1 1 λ th exp λ th +A KA k k λ th +A k N I j=1 j k A j A j A k, where A k = Ω 0 /(K +1)Ω k. This expression is only valid if Ω i Ω j when i j, i.e., the different interferers have different mean power. If all the interferers have the same mean power, then the outage probability can be derived as λ th O I = exp λ th +A KA 1 1 λ th +A 1 N I 1 k=0 A 1 (λ th +A 1 k k m=0 k m 1 m! Kλ th λ th +A 1 If the desired signal is Rayleigh faded, then the probability of outage with multiple Rayleigh faded interferers can be obtained by setting K = 0. m. 19

20 Probability of Outage, O I K = 0, N I = 1 K = 7, N I = 1 K = 7, N I = 2 K = 7, N I = Λ (db) Probability of outage with multiple interferers. The desired signal is Ricean faded with various Rice factors, while the interfering signals are Rayleigh faded and of equal power; λ th = 10.0 db. Λ = Ω 0 N I Ω 1 20

EE6604 Personal & Mobile Communications. Week 7. Path Loss Models. Shadowing

EE6604 Personal & Mobile Communications. Week 7. Path Loss Models. Shadowing EE6604 Personal & Mobile Communications Week 7 Path Loss Models Shadowing 1 Okumura-Hata Model L p = A+Blog 10 (d) A+Blog 10 (d) C A+Blog 10 (d) D for urban area for suburban area for open area where A

More information

ECE6604 PERSONAL & MOBILE COMMUNICATIONS. Lecture 3. Interference and Shadow Margins, Handoff Gain, Coverage

ECE6604 PERSONAL & MOBILE COMMUNICATIONS. Lecture 3. Interference and Shadow Margins, Handoff Gain, Coverage ECE6604 PERSONAL & MOBILE COMMUNICATIONS Lecture 3 Interference and Shadow Margins, Handoff Gain, Coverage 1 Interference Margin As the subscriber load increases, additional interference is generated from

More information

Indoor Propagation Models

Indoor Propagation Models Indoor Propagation Models Outdoor models are not accurate for indoor scenarios. Examples of indoor scenario: home, shopping mall, office building, factory. Ceiling structure, walls, furniture and people

More information

White Paper: Comparison of Narrowband and Ultra Wideband Channels. January 2008

White Paper: Comparison of Narrowband and Ultra Wideband Channels. January 2008 White Paper: Comparison of Narrowband and Ultra Wideband Channels January 28 DOCUMENT APPROVAL: Author signature: Satisfied that this document is fit for purpose, contains sufficient and correct detail

More information

A Model of Coverage Probability under Shadow Fading

A Model of Coverage Probability under Shadow Fading A Model of Coverage Probability under Shadow Fading Kenneth L. Clarkson John D. Hobby August 25, 23 Abstract We give a simple analytic model of coverage probability for CDMA cellular phone systems under

More information

EE Large Scale Path Loss Log Normal Shadowing. The Flat Fading Channel

EE Large Scale Path Loss Log Normal Shadowing. The Flat Fading Channel EE447- Large Scale Path Loss Log Normal Shadowing The Flat Fading Channel The channel functions are random processes and hard to characterize We therefore use the channel correlation functions Now assume:

More information

EE 577: Wireless and Personal Communications

EE 577: Wireless and Personal Communications EE 577: Wireless and Personal Communications Large-Scale Signal Propagation Models 1 Propagation Models Basic Model is to determine the major path loss effects This can be refined to take into account

More information

Probability distributions relevant to radiowave propagation modelling

Probability distributions relevant to radiowave propagation modelling Rec. ITU-R P.57 RECOMMENDATION ITU-R P.57 PROBABILITY DISTRIBUTIONS RELEVANT TO RADIOWAVE PROPAGATION MODELLING (994) Rec. ITU-R P.57 The ITU Radiocommunication Assembly, considering a) that the propagation

More information

Pathloss and Link Budget From Physical Propagation to Multi-Path Fading Statistical Characterization of Channels. P r = P t Gr G t L P

Pathloss and Link Budget From Physical Propagation to Multi-Path Fading Statistical Characterization of Channels. P r = P t Gr G t L P Path Loss I Path loss L P relates the received signal power P r to the transmitted signal power P t : P r = P t Gr G t L P, where G t and G r are antenna gains. I Path loss is very important for cell and

More information

The Impact of Fading on the Outage Probability in Cognitive Radio Networks

The Impact of Fading on the Outage Probability in Cognitive Radio Networks 1 The Impact of Fading on the Outage obability in Cognitive Radio Networks Yaobin Wen, Sergey Loyka and Abbas Yongacoglu Abstract This paper analyzes the outage probability in cognitive radio networks,

More information

Table of Contents. Kocaeli University Computer Engineering Department 2011 Spring Mustafa KIYAR Optimization Theory

Table of Contents. Kocaeli University Computer Engineering Department 2011 Spring Mustafa KIYAR Optimization Theory 1 Table of Contents Estimating Path Loss Exponent and Application with Log Normal Shadowing...2 Abstract...3 1Path Loss Models...4 1.1Free Space Path Loss Model...4 1.1.1Free Space Path Loss Equation:...4

More information

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies

Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies Limit Theorems for the Empirical Distribution Function of Scaled Increments of Itô Semimartingales at high frequencies George Tauchen Duke University Viktor Todorov Northwestern University 2013 Motivation

More information

Equity correlations implied by index options: estimation and model uncertainty analysis

Equity correlations implied by index options: estimation and model uncertainty analysis 1/18 : estimation and model analysis, EDHEC Business School (joint work with Rama COT) Modeling and managing financial risks Paris, 10 13 January 2011 2/18 Outline 1 2 of multi-asset models Solution to

More information

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions.

UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. UQ, STAT2201, 2017, Lectures 3 and 4 Unit 3 Probability Distributions. Random Variables 2 A random variable X is a numerical (integer, real, complex, vector etc.) summary of the outcome of the random experiment.

More information

Bivariate Birnbaum-Saunders Distribution

Bivariate Birnbaum-Saunders Distribution Department of Mathematics & Statistics Indian Institute of Technology Kanpur January 2nd. 2013 Outline 1 Collaborators 2 3 Birnbaum-Saunders Distribution: Introduction & Properties 4 5 Outline 1 Collaborators

More information

2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises

2 Control variates. λe λti λe e λt i where R(t) = t Y 1 Y N(t) is the time from the last event to t. L t = e λr(t) e e λt(t) Exercises 96 ChapterVI. Variance Reduction Methods stochastic volatility ISExSoren5.9 Example.5 (compound poisson processes) Let X(t) = Y + + Y N(t) where {N(t)},Y, Y,... are independent, {N(t)} is Poisson(λ) with

More information

EENG473 Mobile Communications Module 3 : Week # (11) Mobile Radio Propagation: Large-Scale Path Loss

EENG473 Mobile Communications Module 3 : Week # (11) Mobile Radio Propagation: Large-Scale Path Loss EENG473 Mobile Communications Module 3 : Week # (11) Mobile Radio Propagation: Large-Scale Path Loss Practical Link Budget Design using Path Loss Models Most radio propagation models are derived using

More information

Research Article BEP/SEP and Outage Performance Analysis of L-Branch Maximal-Ratio Combiner for κ-μ Fading

Research Article BEP/SEP and Outage Performance Analysis of L-Branch Maximal-Ratio Combiner for κ-μ Fading Hindawi Publishing Corporation International Journal of Digital Multimedia Broadcasting Volume 29, Article ID 57344, 8 pages doi:.55/29/57344 Research Article BEP/ and Outage Performance Analysis of L-Branch

More information

IEOR E4602: Quantitative Risk Management

IEOR E4602: Quantitative Risk Management IEOR E4602: Quantitative Risk Management Basic Concepts and Techniques of Risk Management Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Edinburgh Research Explorer

Edinburgh Research Explorer Edinburgh Research Explorer The Distribution of Path Losses for Uniformly Distributed Nodes in a Circle Citation for published version: Bharucha, Z & Haas, H 2008, 'The Distribution of Path Losses for

More information

FINANCIAL OPTION ANALYSIS HANDOUTS

FINANCIAL OPTION ANALYSIS HANDOUTS FINANCIAL OPTION ANALYSIS HANDOUTS 1 2 FAIR PRICING There is a market for an object called S. The prevailing price today is S 0 = 100. At this price the object S can be bought or sold by anyone for any

More information

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER

Two hours. To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER Two hours MATH20802 To be supplied by the Examinations Office: Mathematical Formula Tables and Statistical Tables THE UNIVERSITY OF MANCHESTER STATISTICAL METHODS Answer any FOUR of the SIX questions.

More information

Statistics Class 15 3/21/2012

Statistics Class 15 3/21/2012 Statistics Class 15 3/21/2012 Quiz 1. Cans of regular Pepsi are labeled to indicate that they contain 12 oz. Data Set 17 in Appendix B lists measured amounts for a sample of Pepsi cans. The same statistics

More information

EELE 6333: Wireless Commuications

EELE 6333: Wireless Commuications EELE 6333: Wireless Commuications Chapter # 2 : Path Loss and Shadowing (Part Two) Spring, 2012/2013 EELE 6333: Wireless Commuications - Ch.2 Dr. Musbah Shaat 1 / 23 Outline 1 Empirical Path Loss Models

More information

Binomial model: numerical algorithm

Binomial model: numerical algorithm Binomial model: numerical algorithm S / 0 C \ 0 S0 u / C \ 1,1 S0 d / S u 0 /, S u 3 0 / 3,3 C \ S0 u d /,1 S u 5 0 4 0 / C 5 5,5 max X S0 u,0 S u C \ 4 4,4 C \ 3 S u d / 0 3, C \ S u d 0 S u d 0 / C 4

More information

Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014

Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014 NTNU Page 1 of 5 Institutt for fysikk Contact during the exam: Professor Ingve Simonsen Exam in TFY4275/FY8907 CLASSICAL TRANSPORT THEORY Feb 14, 2014 Allowed help: Alternativ D All written material This

More information

Financial Risk Management

Financial Risk Management Financial Risk Management Professor: Thierry Roncalli Evry University Assistant: Enareta Kurtbegu Evry University Tutorial exercices #4 1 Correlation and copulas 1. The bivariate Gaussian copula is given

More information

M5MF6. Advanced Methods in Derivatives Pricing

M5MF6. Advanced Methods in Derivatives Pricing Course: Setter: M5MF6 Dr Antoine Jacquier MSc EXAMINATIONS IN MATHEMATICS AND FINANCE DEPARTMENT OF MATHEMATICS April 2016 M5MF6 Advanced Methods in Derivatives Pricing Setter s signature...........................................

More information

SAMPLE STANDARD DEVIATION(s) CHART UNDER THE ASSUMPTION OF MODERATENESS AND ITS PERFORMANCE ANALYSIS

SAMPLE STANDARD DEVIATION(s) CHART UNDER THE ASSUMPTION OF MODERATENESS AND ITS PERFORMANCE ANALYSIS Science SAMPLE STANDARD DEVIATION(s) CHART UNDER THE ASSUMPTION OF MODERATENESS AND ITS PERFORMANCE ANALYSIS Kalpesh S Tailor * * Assistant Professor, Department of Statistics, M K Bhavnagar University,

More information

Estimation Errors and SCR Calculation

Estimation Errors and SCR Calculation Estimation Errors and Calculation E. KARAM & F. PLANCHET February, 2013 Université de Lyon, Université Lyon 1, ISFA, laboratoire SAF EA2429, 69366 Lyon France Abstract The risk of not accurately estimating

More information

The Wireless Communication Channel. Objectives

The Wireless Communication Channel. Objectives The Wireless Communication Channel muse Objectives Understand fundamentals associated with free space propagation. Define key sources of propagation effects both at the large and small scales Understand

More information

Model Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16

Model Estimation. Liuren Wu. Fall, Zicklin School of Business, Baruch College. Liuren Wu Model Estimation Option Pricing, Fall, / 16 Model Estimation Liuren Wu Zicklin School of Business, Baruch College Fall, 2007 Liuren Wu Model Estimation Option Pricing, Fall, 2007 1 / 16 Outline 1 Statistical dynamics 2 Risk-neutral dynamics 3 Joint

More information

Point Estimators. STATISTICS Lecture no. 10. Department of Econometrics FEM UO Brno office 69a, tel

Point Estimators. STATISTICS Lecture no. 10. Department of Econometrics FEM UO Brno office 69a, tel STATISTICS Lecture no. 10 Department of Econometrics FEM UO Brno office 69a, tel. 973 442029 email:jiri.neubauer@unob.cz 8. 12. 2009 Introduction Suppose that we manufacture lightbulbs and we want to state

More information

Which GARCH Model for Option Valuation? By Peter Christoffersen and Kris Jacobs

Which GARCH Model for Option Valuation? By Peter Christoffersen and Kris Jacobs Online Appendix Sample Index Returns Which GARCH Model for Option Valuation? By Peter Christoffersen and Kris Jacobs In order to give an idea of the differences in returns over the sample, Figure A.1 plots

More information

PROPAGATION PATH LOSS IN URBAN AND SUBURBAN AREA

PROPAGATION PATH LOSS IN URBAN AND SUBURBAN AREA PROPAGATION PATH LOSS IN URBAN AND SUBURBAN AREA Divyanshi Singh 1, Dimple 2 UG Student 1,2, Department of Electronics &Communication Engineering Raj Kumar Goel Institute of Technology for Women, Ghaziabad

More information

Asian options and meropmorphic Lévy processes

Asian options and meropmorphic Lévy processes Asian options and meropmorphic Lévy processes July 9 13, 2013 Outline Introduction Existing pricing approaches and the exponential functional I q Methodology overview Meromorphic Lévy processes Theoretical

More information

SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD

SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1. By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The Annals of Applied Probability 1999, Vol. 9, No. 2, 493 53 SADDLEPOINT APPROXIMATIONS TO OPTION PRICES 1 By L. C. G. Rogers and O. Zane University of Bath and First Chicago NBD The use of saddlepoint

More information

Bias Reduction Using the Bootstrap

Bias Reduction Using the Bootstrap Bias Reduction Using the Bootstrap Find f t (i.e., t) so that or E(f t (P, P n ) P) = 0 E(T(P n ) θ(p) + t P) = 0. Change the problem to the sample: whose solution is so the bias-reduced estimate is E(T(P

More information

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model

Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility model 1(23) Valuing volatility and variance swaps for a non-gaussian Ornstein-Uhlenbeck stochastic volatility

More information

Multilevel quasi-monte Carlo path simulation

Multilevel quasi-monte Carlo path simulation Multilevel quasi-monte Carlo path simulation Michael B. Giles and Ben J. Waterhouse Lluís Antoni Jiménez Rugama January 22, 2014 Index 1 Introduction to MLMC Stochastic model Multilevel Monte Carlo Milstein

More information

EE266 Homework 5 Solutions

EE266 Homework 5 Solutions EE, Spring 15-1 Professor S. Lall EE Homework 5 Solutions 1. A refined inventory model. In this problem we consider an inventory model that is more refined than the one you ve seen in the lectures. The

More information

Machine Learning for Quantitative Finance

Machine Learning for Quantitative Finance Machine Learning for Quantitative Finance Fast derivative pricing Sofie Reyners Joint work with Jan De Spiegeleer, Dilip Madan and Wim Schoutens Derivative pricing is time-consuming... Vanilla option pricing

More information

Continous time models and realized variance: Simulations

Continous time models and realized variance: Simulations Continous time models and realized variance: Simulations Asger Lunde Professor Department of Economics and Business Aarhus University September 26, 2016 Continuous-time Stochastic Process: SDEs Building

More information

Strategies for Improving the Efficiency of Monte-Carlo Methods

Strategies for Improving the Efficiency of Monte-Carlo Methods Strategies for Improving the Efficiency of Monte-Carlo Methods Paul J. Atzberger General comments or corrections should be sent to: paulatz@cims.nyu.edu Introduction The Monte-Carlo method is a useful

More information

ECE 295: Lecture 03 Estimation and Confidence Interval

ECE 295: Lecture 03 Estimation and Confidence Interval ECE 295: Lecture 03 Estimation and Confidence Interval Spring 2018 Prof Stanley Chan School of Electrical and Computer Engineering Purdue University 1 / 23 Theme of this Lecture What is Estimation? You

More information

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay. Solutions to Final Exam The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2012, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Consider

More information

Statistical Methods in Financial Risk Management

Statistical Methods in Financial Risk Management Statistical Methods in Financial Risk Management Lecture 1: Mapping Risks to Risk Factors Alexander J. McNeil Maxwell Institute of Mathematical Sciences Heriot-Watt University Edinburgh 2nd Workshop on

More information

SOCIETY OF ACTUARIES Quantitative Finance and Investment Advanced Exam Exam QFIADV AFTERNOON SESSION

SOCIETY OF ACTUARIES Quantitative Finance and Investment Advanced Exam Exam QFIADV AFTERNOON SESSION SOCIETY OF ACTUARIES Exam QFIADV AFTERNOON SESSION Date: Friday, May 2, 2014 Time: 1:30 p.m. 3:45 p.m. INSTRUCTIONS TO CANDIDATES General Instructions 1. This afternoon session consists of 6 questions

More information

Assicurazioni Generali: An Option Pricing Case with NAGARCH

Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: An Option Pricing Case with NAGARCH Assicurazioni Generali: Business Snapshot Find our latest analyses and trade ideas on bsic.it Assicurazioni Generali SpA is an Italy-based insurance

More information

Valuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005

Valuation of Volatility Derivatives. Jim Gatheral Global Derivatives & Risk Management 2005 Paris May 24, 2005 Valuation of Volatility Derivatives Jim Gatheral Global Derivatives & Risk Management 005 Paris May 4, 005 he opinions expressed in this presentation are those of the author alone, and do not necessarily

More information

EELE 5414 Wireless Communications. Chapter 4: Mobile Radio Propagation: Large-Scale Path Loss

EELE 5414 Wireless Communications. Chapter 4: Mobile Radio Propagation: Large-Scale Path Loss EELE 5414 Wireless Communications Chapter 4: Mobile Radio Propagation: Large-Scale Path Loss In the last lecture Outline Diffraction. Scattering. Practical link budget design. Log-distance model Log-normal

More information

Lecture Stat 302 Introduction to Probability - Slides 15

Lecture Stat 302 Introduction to Probability - Slides 15 Lecture Stat 30 Introduction to Probability - Slides 15 AD March 010 AD () March 010 1 / 18 Continuous Random Variable Let X a (real-valued) continuous r.v.. It is characterized by its pdf f : R! [0, )

More information

Information Processing and Limited Liability

Information Processing and Limited Liability Information Processing and Limited Liability Bartosz Maćkowiak European Central Bank and CEPR Mirko Wiederholt Northwestern University January 2012 Abstract Decision-makers often face limited liability

More information

King s College London

King s College London King s College London University Of London This paper is part of an examination of the College counting towards the award of a degree. Examinations are governed by the College Regulations under the authority

More information

Chapter 2 Uncertainty Analysis and Sampling Techniques

Chapter 2 Uncertainty Analysis and Sampling Techniques Chapter 2 Uncertainty Analysis and Sampling Techniques The probabilistic or stochastic modeling (Fig. 2.) iterative loop in the stochastic optimization procedure (Fig..4 in Chap. ) involves:. Specifying

More information

STRESS-STRENGTH RELIABILITY ESTIMATION

STRESS-STRENGTH RELIABILITY ESTIMATION CHAPTER 5 STRESS-STRENGTH RELIABILITY ESTIMATION 5. Introduction There are appliances (every physical component possess an inherent strength) which survive due to their strength. These appliances receive

More information

The Normal Distribution

The Normal Distribution The Normal Distribution The normal distribution plays a central role in probability theory and in statistics. It is often used as a model for the distribution of continuous random variables. Like all models,

More information

Counterparty Credit Risk Simulation

Counterparty Credit Risk Simulation Counterparty Credit Risk Simulation Alex Yang FinPricing http://www.finpricing.com Summary Counterparty Credit Risk Definition Counterparty Credit Risk Measures Monte Carlo Simulation Interest Rate Curve

More information

Calculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the

Calculating VaR. There are several approaches for calculating the Value at Risk figure. The most popular are the VaR Pro and Contra Pro: Easy to calculate and to understand. It is a common language of communication within the organizations as well as outside (e.g. regulators, auditors, shareholders). It is not really

More information

MTH6154 Financial Mathematics I Stochastic Interest Rates

MTH6154 Financial Mathematics I Stochastic Interest Rates MTH6154 Financial Mathematics I Stochastic Interest Rates Contents 4 Stochastic Interest Rates 45 4.1 Fixed Interest Rate Model............................ 45 4.2 Varying Interest Rate Model...........................

More information

Stochastic Processes and Brownian Motion

Stochastic Processes and Brownian Motion A stochastic process Stochastic Processes X = { X(t) } Stochastic Processes and Brownian Motion is a time series of random variables. X(t) (or X t ) is a random variable for each time t and is usually

More information

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations Stan Stilger June 6, 1 Fouque and Tullie use importance sampling for variance reduction in stochastic volatility simulations.

More information

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam

The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay. Solutions to Final Exam The University of Chicago, Booth School of Business Business 41202, Spring Quarter 2017, Mr. Ruey S. Tsay Solutions to Final Exam Problem A: (40 points) Answer briefly the following questions. 1. Describe

More information

Simple Robust Hedging with Nearby Contracts

Simple Robust Hedging with Nearby Contracts Simple Robust Hedging with Nearby Contracts Liuren Wu and Jingyi Zhu Baruch College and University of Utah October 22, 2 at Worcester Polytechnic Institute Wu & Zhu (Baruch & Utah) Robust Hedging with

More information

1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and

1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and CHAPTER 13 Solutions Exercise 1 1. In this exercise, we can easily employ the equations (13.66) (13.70), (13.79) (13.80) and (13.82) (13.86). Also, remember that BDT model will yield a recombining binomial

More information

Lecture 10: Point Estimation

Lecture 10: Point Estimation Lecture 10: Point Estimation MSU-STT-351-Sum-17B (P. Vellaisamy: MSU-STT-351-Sum-17B) Probability & Statistics for Engineers 1 / 31 Basic Concepts of Point Estimation A point estimate of a parameter θ,

More information

Path Loss Models and Link Budget

Path Loss Models and Link Budget Path Loss Models and Link Budget A universal path loss model P r dbm = P t dbm + db Gains db Losses Gains: the antenna gains compared to isotropic antennas Transmitter antenna gain Receiver antenna gain

More information

High-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5]

High-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5] 1 High-Frequency Data Analysis and Market Microstructure [Tsay (2005), chapter 5] High-frequency data have some unique characteristics that do not appear in lower frequencies. At this class we have: Nonsynchronous

More information

Conditional Density Method in the Computation of the Delta with Application to Power Market

Conditional Density Method in the Computation of the Delta with Application to Power Market Conditional Density Method in the Computation of the Delta with Application to Power Market Asma Khedher Centre of Mathematics for Applications Department of Mathematics University of Oslo A joint work

More information

3 ˆθ B = X 1 + X 2 + X 3. 7 a) Find the Bias, Variance and MSE of each estimator. Which estimator is the best according

3 ˆθ B = X 1 + X 2 + X 3. 7 a) Find the Bias, Variance and MSE of each estimator. Which estimator is the best according STAT 345 Spring 2018 Homework 9 - Point Estimation Name: Please adhere to the homework rules as given in the Syllabus. 1. Mean Squared Error. Suppose that X 1, X 2 and X 3 are independent random variables

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Generating Random Variables and Stochastic Processes Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Stochastic Volatility (SV) Models

Stochastic Volatility (SV) Models 1 Motivations Stochastic Volatility (SV) Models Jun Yu Some stylised facts about financial asset return distributions: 1. Distribution is leptokurtic 2. Volatility clustering 3. Volatility responds to

More information

IEOR E4703: Monte-Carlo Simulation

IEOR E4703: Monte-Carlo Simulation IEOR E4703: Monte-Carlo Simulation Simulating Stochastic Differential Equations Martin Haugh Department of Industrial Engineering and Operations Research Columbia University Email: martin.b.haugh@gmail.com

More information

Operational Risk Aggregation

Operational Risk Aggregation Operational Risk Aggregation Professor Carol Alexander Chair of Risk Management and Director of Research, ISMA Centre, University of Reading, UK. Loss model approaches are currently a focus of operational

More information

Contagion models with interacting default intensity processes

Contagion models with interacting default intensity processes Contagion models with interacting default intensity processes Yue Kuen KWOK Hong Kong University of Science and Technology This is a joint work with Kwai Sun Leung. 1 Empirical facts Default of one firm

More information

No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate

No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Fuzzy Optim Decis Making 217 16:221 234 DOI 117/s17-16-9246-8 No-arbitrage theorem for multi-factor uncertain stock model with floating interest rate Xiaoyu Ji 1 Hua Ke 2 Published online: 17 May 216 Springer

More information

The Vasicek Distribution

The Vasicek Distribution The Vasicek Distribution Dirk Tasche Lloyds TSB Bank Corporate Markets Rating Systems dirk.tasche@gmx.net Bristol / London, August 2008 The opinions expressed in this presentation are those of the author

More information

Chapter 2. Random variables. 2.3 Expectation

Chapter 2. Random variables. 2.3 Expectation Random processes - Chapter 2. Random variables 1 Random processes Chapter 2. Random variables 2.3 Expectation 2.3 Expectation Random processes - Chapter 2. Random variables 2 Among the parameters representing

More information

Dependence Modeling and Credit Risk

Dependence Modeling and Credit Risk Dependence Modeling and Credit Risk Paola Mosconi Banca IMI Bocconi University, 20/04/2015 Paola Mosconi Lecture 6 1 / 53 Disclaimer The opinion expressed here are solely those of the author and do not

More information

induced by the Solvency II project

induced by the Solvency II project Asset Les normes allocation IFRS : new en constraints assurance induced by the Solvency II project 36 th International ASTIN Colloquium Zürich September 005 Frédéric PLANCHET Pierre THÉROND ISFA Université

More information

Value at Risk and Self Similarity

Value at Risk and Self Similarity Value at Risk and Self Similarity by Olaf Menkens School of Mathematical Sciences Dublin City University (DCU) St. Andrews, March 17 th, 2009 Value at Risk and Self Similarity 1 1 Introduction The concept

More information

Course information FN3142 Quantitative finance

Course information FN3142 Quantitative finance Course information 015 16 FN314 Quantitative finance This course is aimed at students interested in obtaining a thorough grounding in market finance and related empirical methods. Prerequisite If taken

More information

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets

Chapter 15: Jump Processes and Incomplete Markets. 1 Jumps as One Explanation of Incomplete Markets Chapter 5: Jump Processes and Incomplete Markets Jumps as One Explanation of Incomplete Markets It is easy to argue that Brownian motion paths cannot model actual stock price movements properly in reality,

More information

ELEMENTS OF MONTE CARLO SIMULATION

ELEMENTS OF MONTE CARLO SIMULATION APPENDIX B ELEMENTS OF MONTE CARLO SIMULATION B. GENERAL CONCEPT The basic idea of Monte Carlo simulation is to create a series of experimental samples using a random number sequence. According to the

More information

FINA 695 Assignment 1 Simon Foucher

FINA 695 Assignment 1 Simon Foucher Answer the following questions. Show your work. Due in the class on March 29. (postponed 1 week) You are expected to do the assignment on your own. Please do not take help from others. 1. (a) The current

More information

THE MARTINGALE METHOD DEMYSTIFIED

THE MARTINGALE METHOD DEMYSTIFIED THE MARTINGALE METHOD DEMYSTIFIED SIMON ELLERSGAARD NIELSEN Abstract. We consider the nitty gritty of the martingale approach to option pricing. These notes are largely based upon Björk s Arbitrage Theory

More information

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours

NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 MAS3904. Stochastic Financial Modelling. Time allowed: 2 hours NEWCASTLE UNIVERSITY SCHOOL OF MATHEMATICS, STATISTICS & PHYSICS SEMESTER 1 SPECIMEN 2 Stochastic Financial Modelling Time allowed: 2 hours Candidates should attempt all questions. Marks for each question

More information

Lattice Model of System Evolution. Outline

Lattice Model of System Evolution. Outline Lattice Model of System Evolution Richard de Neufville Professor of Engineering Systems and of Civil and Environmental Engineering MIT Massachusetts Institute of Technology Lattice Model Slide 1 of 48

More information

Numerical Simulation of Stochastic Differential Equations: Lecture 2, Part 2

Numerical Simulation of Stochastic Differential Equations: Lecture 2, Part 2 Numerical Simulation of Stochastic Differential Equations: Lecture 2, Part 2 Des Higham Department of Mathematics University of Strathclyde Montreal, Feb. 2006 p.1/17 Lecture 2, Part 2: Mean Exit Times

More information

Modeling Credit Exposure for Collateralized Counterparties

Modeling Credit Exposure for Collateralized Counterparties Modeling Credit Exposure for Collateralized Counterparties Michael Pykhtin Credit Analytics & Methodology Bank of America Fields Institute Quantitative Finance Seminar Toronto; February 25, 2009 Disclaimer

More information

Analysis of kurtosis-based LOS/NLOS Identification based on indoor MIMO Channel Measurements

Analysis of kurtosis-based LOS/NLOS Identification based on indoor MIMO Channel Measurements Post-print of: Zhang, J., Salmi, J. and Lohan, E-S. Analysis of kurtosis-based LOS/NLOS identification using indoor MIMMO channel measurement in IEEE transactions on vehicular technology, vol. 62, no.

More information

1 Geometric Brownian motion

1 Geometric Brownian motion Copyright c 05 by Karl Sigman Geometric Brownian motion Note that since BM can take on negative values, using it directly for modeling stock prices is questionable. There are other reasons too why BM is

More information

On modelling of electricity spot price

On modelling of electricity spot price , Rüdiger Kiesel and Fred Espen Benth Institute of Energy Trading and Financial Services University of Duisburg-Essen Centre of Mathematics for Applications, University of Oslo 25. August 2010 Introduction

More information

The Binomial Lattice Model for Stocks: Introduction to Option Pricing

The Binomial Lattice Model for Stocks: Introduction to Option Pricing 1/27 The Binomial Lattice Model for Stocks: Introduction to Option Pricing Professor Karl Sigman Columbia University Dept. IEOR New York City USA 2/27 Outline The Binomial Lattice Model (BLM) as a Model

More information

Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days

Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days 1. Introduction Richard D. Christie Department of Electrical Engineering Box 35500 University of Washington Seattle, WA 98195-500 christie@ee.washington.edu

More information

Robust X control chart for monitoring the skewed and contaminated process

Robust X control chart for monitoring the skewed and contaminated process Hacettepe Journal of Mathematics and Statistics Volume 47 (1) (2018), 223 242 Robust X control chart for monitoring the skewed and contaminated process Derya Karagöz Abstract In this paper, we propose

More information

Online Appendix (Not intended for Publication): Federal Reserve Credibility and the Term Structure of Interest Rates

Online Appendix (Not intended for Publication): Federal Reserve Credibility and the Term Structure of Interest Rates Online Appendix Not intended for Publication): Federal Reserve Credibility and the Term Structure of Interest Rates Aeimit Lakdawala Michigan State University Shu Wu University of Kansas August 2017 1

More information

continuous rv Note for a legitimate pdf, we have f (x) 0 and f (x)dx = 1. For a continuous rv, P(X = c) = c f (x)dx = 0, hence

continuous rv Note for a legitimate pdf, we have f (x) 0 and f (x)dx = 1. For a continuous rv, P(X = c) = c f (x)dx = 0, hence continuous rv Let X be a continuous rv. Then a probability distribution or probability density function (pdf) of X is a function f(x) such that for any two numbers a and b with a b, P(a X b) = b a f (x)dx.

More information

Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach

Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach Lei Jiang Tsinghua University Ke Wu Renmin University of China Guofu Zhou Washington University in St. Louis August 2017 Jiang,

More information

2.1 Mathematical Basis: Risk-Neutral Pricing

2.1 Mathematical Basis: Risk-Neutral Pricing Chapter Monte-Carlo Simulation.1 Mathematical Basis: Risk-Neutral Pricing Suppose that F T is the payoff at T for a European-type derivative f. Then the price at times t before T is given by f t = e r(t

More information