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

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

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

Indoor Propagation Models

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

A Model of Coverage Probability under Shadow Fading

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

EE 577: Wireless and Personal Communications

Probability distributions relevant to radiowave propagation modelling

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

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

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

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

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

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

Bivariate Birnbaum-Saunders Distribution

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

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

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

IEOR E4602: Quantitative Risk Management

Edinburgh Research Explorer

FINANCIAL OPTION ANALYSIS HANDOUTS

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

Statistics Class 15 3/21/2012

EELE 6333: Wireless Commuications

Binomial model: numerical algorithm

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

Financial Risk Management

M5MF6. Advanced Methods in Derivatives Pricing

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

Estimation Errors and SCR Calculation

The Wireless Communication Channel. Objectives

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

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

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

PROPAGATION PATH LOSS IN URBAN AND SUBURBAN AREA

Asian options and meropmorphic Lévy processes

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

Bias Reduction Using the Bootstrap

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

Multilevel quasi-monte Carlo path simulation

EE266 Homework 5 Solutions

Machine Learning for Quantitative Finance

Continous time models and realized variance: Simulations

Strategies for Improving the Efficiency of Monte-Carlo Methods

ECE 295: Lecture 03 Estimation and Confidence Interval

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

Statistical Methods in Financial Risk Management

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

Assicurazioni Generali: An Option Pricing Case with NAGARCH

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

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

Lecture Stat 302 Introduction to Probability - Slides 15

Information Processing and Limited Liability

King s College London

Chapter 2 Uncertainty Analysis and Sampling Techniques

STRESS-STRENGTH RELIABILITY ESTIMATION

The Normal Distribution

Counterparty Credit Risk Simulation

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

MTH6154 Financial Mathematics I Stochastic Interest Rates

Stochastic Processes and Brownian Motion

The Use of Importance Sampling to Speed Up Stochastic Volatility Simulations

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

Simple Robust Hedging with Nearby Contracts

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

Lecture 10: Point Estimation

Path Loss Models and Link Budget

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

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

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

IEOR E4703: Monte-Carlo Simulation

Stochastic Volatility (SV) Models

IEOR E4703: Monte-Carlo Simulation

Operational Risk Aggregation

Contagion models with interacting default intensity processes

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

The Vasicek Distribution

Chapter 2. Random variables. 2.3 Expectation

Dependence Modeling and Credit Risk

induced by the Solvency II project

Value at Risk and Self Similarity

Course information FN3142 Quantitative finance

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

ELEMENTS OF MONTE CARLO SIMULATION

FINA 695 Assignment 1 Simon Foucher

THE MARTINGALE METHOD DEMYSTIFIED

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

Lattice Model of System Evolution. Outline

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

Modeling Credit Exposure for Collateralized Counterparties

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

1 Geometric Brownian motion

On modelling of electricity spot price

The Binomial Lattice Model for Stocks: Introduction to Option Pricing

Maximum Likelihood Estimates for Alpha and Beta With Zero SAIDI Days

Robust X control chart for monitoring the skewed and contaminated process

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

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

Internet Appendix for Asymmetry in Stock Comovements: An Entropy Approach

2.1 Mathematical Basis: Risk-Neutral Pricing

Transcription:

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

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

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

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

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 = 0.23026. 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

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

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 ξ = 0.23026. 7

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

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

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

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

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

Outage with Multiple Interferers 1. First obtain the mean and variance µ Z = µẑ/ξ σ 2 Z = σ 2 Ẑ /ξ2 ξ = 0.23026 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

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

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

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 = 1 3 10 β Λ 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

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

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

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

Probability of Outage, O I 10 0 10 1 10 2 10 3 K = 0, N I = 1 K = 7, N I = 1 K = 7, N I = 2 K = 7, N I = 6 10 4 10 15 20 25 30 35 Λ (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