EE6604 Personal & Mobile Communications. Week 9. Co-Channel Interference
|
|
- Julia Hawkins
- 5 years ago
- Views:
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 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 informationECE6604 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 informationIndoor 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 informationWhite 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 informationA 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 informationEE 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 informationEE 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 informationProbability 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 informationPathloss 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 informationThe 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 informationTable 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 informationLimit 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 informationEquity 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 informationUQ, 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 informationBivariate 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 information2 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 informationEENG473 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 informationResearch 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 informationIEOR 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 informationEdinburgh 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 informationFINANCIAL 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 informationTwo 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 informationStatistics 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 informationEELE 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 informationBinomial 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 informationExam 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 informationFinancial 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 informationM5MF6. 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 informationSAMPLE 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 informationEstimation 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 informationThe 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 informationModel 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 informationPoint 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 informationWhich 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 informationPROPAGATION 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 informationAsian 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 informationSADDLEPOINT 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 informationBias 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 informationValuing 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 informationMultilevel 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 informationEE266 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 informationMachine 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 informationContinous 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 informationStrategies 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 informationECE 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 informationThe 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 informationStatistical 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 informationSOCIETY 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 informationAssicurazioni 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 informationValuation 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 informationEELE 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 informationLecture 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 informationInformation 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 informationKing 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 informationChapter 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 informationSTRESS-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 informationThe 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 informationCounterparty 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 informationCalculating 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 informationMTH6154 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 informationStochastic 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 informationThe 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 informationThe 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 informationSimple 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 information1. 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 informationLecture 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 informationPath 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 informationHigh-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 informationConditional 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 information3 ˆθ 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 informationIEOR 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 informationStochastic 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 informationIEOR 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 informationOperational 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 informationContagion 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 informationNo-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 informationThe 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 informationChapter 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 informationDependence 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 informationinduced 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 informationValue 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 informationCourse 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 informationChapter 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 informationELEMENTS 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 informationFINA 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 informationTHE 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 informationNEWCASTLE 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 informationLattice 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 informationNumerical 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 informationModeling 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 informationAnalysis 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 information1 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 informationOn 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 informationThe 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 informationMaximum 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 informationRobust 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 informationOnline 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 informationcontinuous 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 informationInternet 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 information2.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