Collision and Packet Loss Analysis in a LoRaWAN Network
|
|
- Charles Caldwell
- 5 years ago
- Views:
Transcription
1 7 5th Europea Sigal Processig Cferece EUSIPCO Collisi ad Packet Loss Aalysis i a LoRaWAN Network Guillaume Ferré IMS Laboratory - Uiversity of Bordeaux - Bordeaux INP 35 cours de la libérati 334 Talece FRANCE guillaume.ferre@ims-bordeaux.fr Abstract Iteret of thigs IoT is csidered as the ext techological revoluti. Therefore, may solutis are developed either i free, i.e. ISM bads or i free bads with the ultimate aim of affordig cectivity over several kilometers. Based this feature, i urba evirmet the desity of IoT devices will be extremely high. I this paper we propose to aalyze the collisi ad packet loss whe LoRaWAN is csidered. Based the LoRaWAN features, we develop closed-form expressis of collisi ad packet loss probabilities. Simulati results cfirm our theoretical developmets. We also show that our theoretical expressis are more accurate tha the Poiss distributed process to describe the collisis. I. INTRODUCTION The Iteret of thigs IoT deploymet i free bads is maily based two commuicati techologies which are Sigfox [ ad LoRa [. Such techologies are called low power wide area etwork LPWAN ad they share the same objectives, i.e the establishmet of lg rage, low power ad low data rate commuicatis. Despite these similarities, they are techically ad ecomically opposites. Ideed the Sigfox physical layer is based ultra arrow baommuicati, as LoRa uses a spreadig spectrum techique to exchage iformati. I this paper we focus our atteti LoRa i order to aalyze the collisi ad packet loss whe LoRaWAN is csidered. Ideed, if we believe to the IoT success i our every day life, it is predictable to have several millis of IoT devices usig such modulati techiques. Therefore, the coexistece of all these devices, whe the commuicati is performed usig the same rage, will threat IoT commuicati with iterfereces. The chael access i LoRaWAN class A is based a ALOHA priciple. For etworks usig a ALOHA type chael access, may packet collisi studies have bee de before. The csidered system i this paper is differet tha the e used for classical aalysis of packet collisi i a u-slotted ALOHA based protocol [3. Ideed, packet time air ad time betwee two successive trasmissis deped the IoT applicati. Thus, the model allowig to describe the collisi effect based Poiss distributed process PDP is ot accurate eough. However, i our case this modelig is useful to defie lower ad upper bouds of the probability of success for a give spreadig factor. I this paper we propose a more accurate ad specific approach to predict the collisi ad packet loss. Thus, based the LoRaWAN MAC mechaisms we develop closed-form expressis of the probability of collisi ad packet loss ad our theoretical aalyses are cfirmed by the simulatis. The paper is orgaized as follows. I secti II the LoRa modulati ad the LoRaWAN MAC layer are itroduced. I secti III we develop the received sigal model that we csider. Based this model, secti IV is devoted to the collisi ad packet loss theoritical aalyse. Simulati results are give i secti V. II. LORA AND LORAWAN The followig sectis are dedicated to itroduce the LoRa physical ad MAC layers which is based the patet [4. For more details the reader is referred to [5. It should be oted that the LoRa physical Layer is ot published. A. LoRa physical layer priciple LoRa is based Chirp Spread Spectrum CSS modulati. CSS was proposed for the first time for commuicati systems by Wikler [6 ad applicati to digital commuicati by Beri [7. CSS is csidered as a subcategory of Direct-Sequece Spread Spectrum. CSS is compliat with IoT etwork eeds because it permits to come over the receiver s sesitivity issue ad icrease the commuicati rage at the cost of a reduced spectral efficiecy. The spectrum spreadig i LoRa is achieved usig a chirp sigal that ca be described by its istataeous phase φt or a specific time fucti f c t. f c t is called the raw chirp that: icreases liearly, for a up raw chirp, from a iitial value B to a fial value B, decreases liearly, for a dow raw chirp, from a iitial value B to a fial value B, B stads for the ISM sigal badwidth used for the commuicati. The raw chirp time durati is equal to the symbol period T s. f c t is defied as follows: f ct = ± B T s t The relatiship betwee the badwidth ad the symbol period is give by T s = SF B SF stads for the spreadig factor expet SF [7,...,. Let D s be the symbol rate of the B depeds the used ISM bad aa be chose equal to 5, 5 or 5 khz ISBN EURASIP 7 655
2 7 5th Europea Sigal Processig Cferece EUSIPCO trasmitted sigal ad D b the bit rate, the D s = D b /SF. Lger rage is achieved by varyig the spreadig factor, however to meet highly robust commuicati it is possible to vary the codig rate. With LoRa, symbols are obtaied from a biary combiati of SF bits. Each symbol is associated to a uique chirp. The differet chirps are orthogal to each other i order to retrieve at the receiver the symbols without iter-symbol iterferece IES. If we ote M the set of symbols, the chirp associated to the symbol m, m [, M, is obtaied by delayig the raw chirp f c t by τ m = m B. The chirp outside [ Ts, Ts Ts is cyclically shifted i the iterval [, Ts +τ m. Thus, the chirp associated to the trasmissi of the m th symbol is decomposed of parts: from t [ Ts, Ts + τ m[, raw chirp up ou dow advaced of T s τ m, from t [ Ts + τ m, Ts, raw chirp up ou dow delayed of τ m. For a up chirp, we obtai: f m c t = B T s t τ m + B f m c t = B T s t τ m for t [ T s, T s + m B [ for t [ T s + m B, T s Thus, the expressi of the basebad trasmitted sigal by the ode is give as follows: r t = k S l e jπf c,kt kt st kt s+jφ f c,k t represets the trasmittehirp at time kt s, S l the set of trasmitted symbols iside the packet p ad φ a iitial phase. If we ote K the S l size, thus the silet time durati of the LoRa ode will be at least equal to KTs KT s. Thus, from t + KT s to t + KTs the ode will be silet, is the duty cycle. B. LoRaWAN: a LoRa Mac layer LoRaWAN is a ope stadard developed by the LoRa Alliace. It s e of the possible MAC layer for the LoRa modulati ad obviously the well kow. The LoRaWAN specificati defies 3 categories of odes: Class A: a basic class of LoRa that is implemeted i all LoRa chips. It allows bi-directial commuicatis which is usually origially started by the ode i a asychrous way. The uplik trasmissi triggers two short dowlik receive widows. The trasmissi slot is scheduled whe eeded by the ode i a radom time basis. Accordig to LoRaWAN specificatis, class A is a ALOHA based-protocol. Class B: this class is cceived to guaratee uplik ad dowlik separati. Nodes are sychrized usig a beac trasmitted by the gateway. Thus, they ca receive iformati from Iteret without sedig requests. Class C: the ode has ctiuously ope receive widows that are closed ly while trasmittig. Compared to A ad B classes, C class csumes more eergy to operate but it offers the lowest latecy. The packet time durati is called the time air. This value depeds several parameters, such SF, B, the size of the payload, the codig rate, etc. After each uplik packet trasmissi, the ode waits for a gateway ackowledgmet ACK of the correct packet recepti. LoRaWAN allows two possible ackowledgmets two differet chaels. The first is trasmitted by the gateway with a cstat delay of T after the ed of the uplik packet recepti. The gateway uses the same chael tha the precedig uplik for this ACK. The secd ACK is trasmitted a differet chael tha the uplik after a time T > T + T g, T g is the time air of the packet trasmitted by the gateway. For more details the LoRaWAN MAC structure see [6. III. RECEIVED SIGNAL MODEL Based the complex evelope defiiti of, we express i this secti the model of the received sigal that we use. We csider a system which is composed of e LoRa gateway ad N LoRa odes. The th ode trasmit P packets, ad p represets the packet umber p trasmitted by the ode with a time air T. The badwidth is arrow eough to make the assumpti of flat fadig propagati chael. We ote h the chael coefficiet associated to the ode. If xt is the received sigal at the LoRa gateway, we have: N P xt = h r t Tp 3 = p= Let Tp be the begiig of the LoRa packet p. Based this, we defie: p Tp = p T + Toff + T r u = p T + u= p T r u 4 u= T r u is a statiary iid radom process with etries uiformly distributed. T r u is used for modelig the differet ode cstraits i terms of commuicatis with the gateways. I the followig sectis we csider that T r u is uiformly distributed i the set [T mi, T. This assumpti is justified due to the huge umber of IoT applicatis. T ad Toff are the time air ad time off air, respectively. IV. COLLISIONS AND PACKET LOSS AT THE LORAWAN GATEWAY A. Probability of Collisi at the LoRaWAN gateway I the followig we express the probability of collisi betwee p ad, give that p represets the ode of iterest. If C p deotes this collisi, based 4, the collisi evet is defied as follows: C p = {Tp Ω C = [Tp T, Tp + T} Based this we focus our atteti ly T. ISBN EURASIP 7 656
3 7 5th Europea Sigal Processig Cferece EUSIPCO A collisi betwee betwee p ad will happe if Tp Ω C. P C p represets the probability of this collisi ad we have: PC p = Tp Ω C 5 The radom variable Tp is the summati of multiple uiform radom variables see 4. For p >> we suppose, usig the cetral limit theorem, that T = p u= T ru is Gaussia distributed with a mea ad a variace µ = p T + T mi 6 = p T T mi 7 Uder this assumpti, Tp is ow a summati betwee two radom variables T ad T ad T T N µ, 8 UT mi, T 9 Ideed, due to the varity of IoT applicatis, the time air T of the ode ca be csidered uiformly distributed betwee the miimal time air T mi ad the imal time air T. Based this, PC p ca be rewritte give T as PC p = T T = T T mi Tp Ω C T dt T mi. With 8 we have T Tp Ω C T p T = Q µ T p + T Q µ µ = µ + p T give Qx = Q x, ca be rewritte as follows: PC p = T T T mi T T mi Q αt + β dt Q γt + δ dt α = p 3 β = [ µ T T + p 4 γ = p + = α + 5 δ = [ µ T T p = β + T 6 with the two chage of variables y = αt + β ad z = γt + δ, is rewritte as PC p = y Q y dy α T y y = αt mi z = γt mi γ T + β, y = αt + δ, z = γt + β + δ z z Q z dz usig: [ x Qxdx = x Qx x Qx + e x e x x π we get the probability of collisi: 7 PC p = fy, y fz, z 8 α T γ T PC p represets the probability of collisi betwee p ad. As we csider p as the ode of iterest, a fiite set Ω p of packets p of the ode have a ull probability to be i collisi with p. We ote Np the umber of packet of that ca be i collisi with p. Thus, the more we icrease p, the higher Np will be. For p, the evets C p with p Ω p are disjoit. If we ote C p the evet associated to the collisi of p with, we have C p = p Ω p C p 9 Thus the probability of collisi of p with is: PC p = PC p p Ω p Numerically, we observed that is idepedet of p. This result is iterestig because it demstrates that is ot depedet the packet umber p, csequetly it does t deped the time. Logically PC p is a fucti of T, ad the T spreadig i time i.e the support of the distributi. Based we express the probability of at least e collisi i a evirmet composed of N idepedet odes. Thus the probability of havig at least e collisi betwee ad the rest of the odes usig the same SF is give by: N [ P c N = PC p = B. Packet loss at the LoRaWAN gateway If o collisi occurs p we csider that the packet will be ackowledge by the gateway. I this case, after a fixed delay T the LoRaWAN gateway will use the same badwidth to trasmit the ACK. We ote T g the time air used by the gateway to aswer ad we suppose that it ca be csidered as a cstat value. Durig T g, the gateway is uable to receive a packet. Thus, a packet loss will occurs if a ode sed a ISBN EURASIP 7 657
4 7 5th Europea Sigal Processig Cferece EUSIPCO iformati durig the gateway ACK of p. This meas that the packet loss of e ode depeds the success of aother ode. Our goal is to express the probability of packet loss of based the success of p. This will occur whe: T p Ω L = [T p + T, τ T, τ + T g τ = Tp + T + T, ad x, y refers to x if x > y ad y otherwise. Based 3 cases must be csidered: case whe T < T mi, case whe T > T ad case 3 whe T [T mi, T. We defie the probability of packet loss of as P L p, L p is the evet associated to the loss give the success of p. Usig the same approach as the e developed i the previous secti, we have : I case : P L p p = T p T Ω L T dt I case : P L p p = T I case 3: P L p T mi T mi =.5 T +.5 T T mi T Tp Ω L T dt Tp Ω L T Tp Ω L 3 4 dt 5 T dt 6 Ω L = [T p + T, τ + T, g Ω L = [τ T, τ + T, g T = T T mi ad T = T T. Based 8 we have: Tp Ω L T = Q αt + µ τ T g Q αt + β 7 ad Tp Ω L T = Q αt T g Q γt 8 Followig the same maipulatis as the es from to 8, we ca express P for the 3 previous cases. L p Regardless the value of T, we ote L p the evet associated to the packet loss of give the success of p. Thus we have: P L p = P L p 9 p Ω l p Ω l p is the set of p that ca be lost give the success of p. Thus the probability of packet loss i.e a collisi with the ACK trasmitted by the gateway to p give the success of p is: P l N = P s N N = [ P L p 3 P s N = P c N is the probability of success of p i.e o collisi. V. SIMULATION RESULTS We csider a system which is composed of e gateway that ca establish cectis with N odes usig the same spreadig factor. We also csider that the sigals trasmitted by odes usig differet SF are orthogal. The odes access the chael radomly i time. The badwidth is fixed to B = 5kHz ad = %. T is uiformly distributed. The support of his distributi begis from the miimum time air to the imum time air. These values are reported i table I ad obtaied usig the lie LoRaWAN calculator. The miimum ad imum umber of packet per day are fixed by the miimum ad imum chael time durati used, respectively. For our simulatis we csider a miimum chael time use of 3s ad a imum of 864s. Time air Payload size bytes SF mi ms ms mi Table I: LoRaWAN imum ad miimum time air ad associated umber of bytes i the payload. Time air values obtaied from the LoRaWAN calculator. First of all, we compare our theoretical results with simulatis. Fig. shows the evoluti of PC p whe p =, SF = for T = s ad s. As it ca be see, our theoretical results are csistet with the simulatis. Ideed, we verify the accuracy of the theoretical aalysis, both theoretical expressi ad simulati curves coicide. I fig. we show the evoluti of P c N. O the e had we compare our theoretical results with the simulati whe T = s ad s. I both cases, simulati ad theory are superimposed. O the other had, we compare our results with a Poiss distributed process PDP geerally use to describe collisis i such a etwork ad defied as follows: D p = T P P DP N = e T + T+Tmi. Dp N As we ca see whe T is closed to the mea value of T, our approach ad the e based PDP give the same results. However, whe it s ot the case the PDP based approach is ot accurate to predict the collisi. Ideed, whe T = s we ca see a differece with our results. For example, for P c N =.3 our theoretical results ad also the simulati show that 58 odes ca be maaged by the gateway, as the PDP predict 8 odes. I fig. 3 we compare the probability of at least e collisi for a payload size betwee ad 59 bytes, whe odes are csidered. These sizes are available for each SF i LoRaWAN. The results ca be grouped i three parts providig the same performace. The first is composed of SF = 7 ad ISBN EURASIP 7 658
5 7 5th Europea Sigal Processig Cferece EUSIPCO 6 5 to a importat umber of collisi at the gateway. C p Probability of collisi of p - P T = s, SF = - Theory T = s, SF = - Simulati T = s, SF = - Theory T = s, SF = - Simulati Probability of at least e packet loss T = T mi, T = s T = T mi, T = T mi T = T mi, T T T = T, T = s T = T, T = T mi T = T, T T Packet umber p Figure : Probability of collisi P C p whe p =, SF = ad for T = s ad s 4 6 8,,,4,6,8,,,4,6,8 3, Number of odes Figure 4: Probability of at least e packet loss, SF = Probability of at least e collisi T = s, SF = - Theory T = s, SF = - Simulati T = s, Poiss distributed process T = s, SF = - Theory T = s, SF = - Simulati T = s, Poiss distributed process 5,,5,,5 3, Number of odes Figure : Probability of at least e collisi for SF =, whe T = s ad s 8, the secd of SF = 9 ad the third of SF =,,. From the probability of collisi poit of view it s better to use SF = 8, SF = 9 or SF =. Ideed, SF = ad give the same probability of collisi tha SF = but they offer a lower sesitivity. Probability of at least e collisi SF=7 SF=8 SF=9 SF= SF= SF= Payload size i Bytes Figure 3: Probability of at least e collisi for the differet SF whe odes are csidered I fig. 4 ad 5 we show the evoluti of P l N whe SF = ad SF = 7, respectively. This allows to show the worst ad the best case. For each SF, T g = T mi which correspds to a ACK composed of byte. We ca see that the P l N is bouded ad pass through a imum value. From the etwork poit of view ly the icreasig part of the curves are iterestig. Ideed the decreasig part correspds Probability of at least e packet loss ,,,4,6,8,,,4,6,8 3, Number of odes T = T mi, T = s T = T mi, T = T mi T = T mi, T T T = T, T = s T = T, T = T mi T = T, T T Figure 5: Probability of at least e packet loss, SF = 7 VI. CONCLUSION I this paper we proposed a aalysis of packet collisi ad loss whe LoRaWAN is csidered. Based the LoRaWAN features, we developed theoretical expressis for both the collisi ad the packet loss. These developmets have bee cfirmed by simulatis results. We have also showed that our approach allows to more accurately describe the collisi tha the classical PDP approach. The perspectives of this work are umerous. We are curretly workig the theoretical demstrati of the idepedece of PC p i p ad the developmet of a software that ca predict the iterferece level of each odes based their locatis, duty-cycle, time air, etc. Csider that the odes with differet SF ca be orthogal is also a iterestig perspective of this work. REFERENCES [ SigFox, 6, [accessed [Olie. Available: [ Lora modem desig guide : Sx7/3/6/7/8. [Olie. Available: [3 R. Rom ad M. Sidi, Multiple Access Protocols: Performace ad Aalysis. New York, NY, USA: Spriger-Verlag New York, Ic., 99. [4 O. Seller ad N. Sori, Low power lg rage trasmitter, Aug. 7 4, us Patet App. 4/7,7. [Olie. Available: [5 A. Augusti, J. Yi, T. Clause, ad W. M. Towsley, A study of lora: Lg rage ad low power etworks for the iteret of thigs, Sesors, vol. 6, o. 9, p. 466, 6. [Olie. Available: [6 M. Wikler, Chirp sigals for commuicatis, IEEE WESCON Cveti Record, p. 7, 96. [7 A. Beri ad W. Gregg, O the utility of chirp modulati for digital sigalig, IEEE Trasactis Commuicatis, vol., o. 6, pp , Jue 973. ISBN EURASIP 7 659
1 of 7 7/16/2009 6:06 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 6. Order Statistics Defiitios Suppose agai that we have a basic radom experimet, ad that X is a real-valued radom variable
More informationAn Introduction to Randomized Algorithms
A Itroductio to Radomized Algorithms The focus of this lecture is to study a radomized algorithm for quick sort, aalyze it usig probabilistic recurrece relatios, ad also provide more geeral tools for aalysis
More informationLecture 27. Capacity of additive Gaussian noise channel and the sphere packing bound
Lecture 7 Ageda for the lecture Gaussia chael with average power costraits Capacity of additive Gaussia oise chael ad the sphere packig boud 7. Additive Gaussia oise chael Up to this poit, we have bee
More informationApproximations and more PMFs and PDFs
Approximatios ad more PMFs ad PDFs Saad Meimeh 1 Approximatio of biomial with Poisso Cosider the biomial distributio ( b(k,,p = p k (1 p k, k λ: k Assume that is large, ad p is small, but p λ at the limit.
More informationECE 564/645 - Digital Communication Systems (Spring 2014) Final Exam Friday, May 2nd, 8:00-10:00am, Marston 220
ECE 564/645 - Digital Commuicatio Systems (Sprig 014) Fial Exam Friday, May d, 8:00-10:00am, Marsto 0 Overview The exam cosists of four (or five) problems for 100 (or 10) poits. The poits for each part
More informationThe Maximum-Likelihood Decoding Performance of Error-Correcting Codes
The Maximum-Lielihood Decodig Performace of Error-Correctig Codes Hery D. Pfister ECE Departmet Texas A&M Uiversity August 27th, 2007 (rev. 0) November 2st, 203 (rev. ) Performace of Codes. Notatio X,
More informationProblem Set 4 Due Oct, 12
EE226: Radom Processes i Systems Lecturer: Jea C. Walrad Problem Set 4 Due Oct, 12 Fall 06 GSI: Assae Gueye This problem set essetially reviews detectio theory ad hypothesis testig ad some basic otios
More informationIntegrable Functions. { f n } is called a determining sequence for f. If f is integrable with respect to, then f d does exist as a finite real number
MATH 532 Itegrable Fuctios Dr. Neal, WKU We ow shall defie what it meas for a measurable fuctio to be itegrable, show that all itegral properties of simple fuctios still hold, ad the give some coditios
More informationSignals & Systems Chapter3
Sigals & Systems Chapter3 1.2 Discrete-Time (D-T) Sigals Electroic systems do most of the processig of a sigal usig a computer. A computer ca t directly process a C-T sigal but istead eeds a stream of
More information1 Approximating Integrals using Taylor Polynomials
Seughee Ye Ma 8: Week 7 Nov Week 7 Summary This week, we will lear how we ca approximate itegrals usig Taylor series ad umerical methods. Topics Page Approximatig Itegrals usig Taylor Polyomials. Defiitios................................................
More informationLecture 14: Graph Entropy
15-859: Iformatio Theory ad Applicatios i TCS Sprig 2013 Lecture 14: Graph Etropy March 19, 2013 Lecturer: Mahdi Cheraghchi Scribe: Euiwoog Lee 1 Recap Bergma s boud o the permaet Shearer s Lemma Number
More informationStochastic Simulation
Stochastic Simulatio 1 Itroductio Readig Assigmet: Read Chapter 1 of text. We shall itroduce may of the key issues to be discussed i this course via a couple of model problems. Model Problem 1 (Jackso
More informationEECS564 Estimation, Filtering, and Detection Hwk 2 Solns. Winter p θ (z) = (2θz + 1 θ), 0 z 1
EECS564 Estimatio, Filterig, ad Detectio Hwk 2 Sols. Witer 25 4. Let Z be a sigle observatio havig desity fuctio where. p (z) = (2z + ), z (a) Assumig that is a oradom parameter, fid ad plot the maximum
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2016 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationBER results for a narrowband multiuser receiver based on successive subtraction for M-PSK modulated signals
results for a arrowbad multiuser receiver based o successive subtractio for M-PSK modulated sigals Gerard J.M. Jasse Telecomm. ad Traffic-Cotrol Systems Group Dept. of Iformatio Techology ad Systems Delft
More informationConvergence of random variables. (telegram style notes) P.J.C. Spreij
Covergece of radom variables (telegram style otes).j.c. Spreij this versio: September 6, 2005 Itroductio As we kow, radom variables are by defiitio measurable fuctios o some uderlyig measurable space
More informationFig. 2. Block Diagram of a DCS
Iformatio source Optioal Essetial From other sources Spread code ge. Format A/D Source ecode Ecrypt Auth. Chael ecode Pulse modu. Multiplex Badpass modu. Spread spectrum modu. X M m i Digital iput Digital
More informationDistribution of Random Samples & Limit theorems
STAT/MATH 395 A - PROBABILITY II UW Witer Quarter 2017 Néhémy Lim Distributio of Radom Samples & Limit theorems 1 Distributio of i.i.d. Samples Motivatig example. Assume that the goal of a study is to
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More informationStatistical Noise Models and Diagnostics
L. Yaroslavsky: Advaced Image Processig Lab: A Tutorial, EUSIPCO2 LECTURE 2 Statistical oise Models ad Diagostics 2. Statistical models of radom iterfereces: (i) Additive sigal idepedet oise model: r =
More informationt distribution [34] : used to test a mean against an hypothesized value (H 0 : µ = µ 0 ) or the difference
EXST30 Backgroud material Page From the textbook The Statistical Sleuth Mea [0]: I your text the word mea deotes a populatio mea (µ) while the work average deotes a sample average ( ). Variace [0]: The
More informationIt is always the case that unions, intersections, complements, and set differences are preserved by the inverse image of a function.
MATH 532 Measurable Fuctios Dr. Neal, WKU Throughout, let ( X, F, µ) be a measure space ad let (!, F, P ) deote the special case of a probability space. We shall ow begi to study real-valued fuctios defied
More informationCS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 5
CS434a/54a: Patter Recogitio Prof. Olga Veksler Lecture 5 Today Itroductio to parameter estimatio Two methods for parameter estimatio Maimum Likelihood Estimatio Bayesia Estimatio Itroducto Bayesia Decisio
More informationLecture 6 Chi Square Distribution (χ 2 ) and Least Squares Fitting
Lecture 6 Chi Square Distributio (χ ) ad Least Squares Fittig Chi Square Distributio (χ ) Suppose: We have a set of measuremets {x 1, x, x }. We kow the true value of each x i (x t1, x t, x t ). We would
More information6.3 Testing Series With Positive Terms
6.3. TESTING SERIES WITH POSITIVE TERMS 307 6.3 Testig Series With Positive Terms 6.3. Review of what is kow up to ow I theory, testig a series a i for covergece amouts to fidig the i= sequece of partial
More informationEnterprise interoperability measurement - Basic concepts
Eterprise iteroperability measuremet - Basic cocepts Nicolas Dacli, David Che, Bruo Vallespir LAPS/GRAI, Uiversity Bordeaux 1, ENSEIRB, UMR 5131 CNRS, 351 cours de la Libératio, 33405 Talece Cedex, Frace
More informationLecture 2: Monte Carlo Simulation
STAT/Q SCI 43: Itroductio to Resamplig ethods Sprig 27 Istructor: Ye-Chi Che Lecture 2: ote Carlo Simulatio 2 ote Carlo Itegratio Assume we wat to evaluate the followig itegratio: e x3 dx What ca we do?
More informationAxioms of Measure Theory
MATH 532 Axioms of Measure Theory Dr. Neal, WKU I. The Space Throughout the course, we shall let X deote a geeric o-empty set. I geeral, we shall ot assume that ay algebraic structure exists o X so that
More informationLecture 5: April 17, 2013
TTIC/CMSC 350 Mathematical Toolkit Sprig 203 Madhur Tulsiai Lecture 5: April 7, 203 Scribe: Somaye Hashemifar Cheroff bouds recap We recall the Cheroff/Hoeffdig bouds we derived i the last lecture idepedet
More informationLecture 3: August 31
36-705: Itermediate Statistics Fall 018 Lecturer: Siva Balakrisha Lecture 3: August 31 This lecture will be mostly a summary of other useful expoetial tail bouds We will ot prove ay of these i lecture,
More informationOn Random Line Segments in the Unit Square
O Radom Lie Segmets i the Uit Square Thomas A. Courtade Departmet of Electrical Egieerig Uiversity of Califoria Los Ageles, Califoria 90095 Email: tacourta@ee.ucla.edu I. INTRODUCTION Let Q = [0, 1] [0,
More informationSECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES
SECTION 1.5 : SUMMATION NOTATION + WORK WITH SEQUENCES Read Sectio 1.5 (pages 5 9) Overview I Sectio 1.5 we lear to work with summatio otatio ad formulas. We will also itroduce a brief overview of sequeces,
More informationELEC1200: A System View of Communications: from Signals to Packets Lecture 3
ELEC2: A System View of Commuicatios: from Sigals to Packets Lecture 3 Commuicatio chaels Discrete time Chael Modelig the chael Liear Time Ivariat Systems Step Respose Respose to sigle bit Respose to geeral
More informationInfinite Sequences and Series
Chapter 6 Ifiite Sequeces ad Series 6.1 Ifiite Sequeces 6.1.1 Elemetary Cocepts Simply speakig, a sequece is a ordered list of umbers writte: {a 1, a 2, a 3,...a, a +1,...} where the elemets a i represet
More information32 estimating the cumulative distribution function
32 estimatig the cumulative distributio fuctio 4.6 types of cofidece itervals/bads Let F be a class of distributio fuctios F ad let θ be some quatity of iterest, such as the mea of F or the whole fuctio
More informationStat 421-SP2012 Interval Estimation Section
Stat 41-SP01 Iterval Estimatio Sectio 11.1-11. We ow uderstad (Chapter 10) how to fid poit estimators of a ukow parameter. o However, a poit estimate does ot provide ay iformatio about the ucertaity (possible
More informationOutput Analysis (2, Chapters 10 &11 Law)
B. Maddah ENMG 6 Simulatio Output Aalysis (, Chapters 10 &11 Law) Comparig alterative system cofiguratio Sice the output of a simulatio is radom, the comparig differet systems via simulatio should be doe
More informationENGI 4421 Confidence Intervals (Two Samples) Page 12-01
ENGI 44 Cofidece Itervals (Two Samples) Page -0 Two Sample Cofidece Iterval for a Differece i Populatio Meas [Navidi sectios 5.4-5.7; Devore chapter 9] From the cetral limit theorem, we kow that, for sufficietly
More informationHomework 5 Solutions
Homework 5 Solutios p329 # 12 No. To estimate the chace you eed the expected value ad stadard error. To do get the expected value you eed the average of the box ad to get the stadard error you eed the
More informationRandom Variables, Sampling and Estimation
Chapter 1 Radom Variables, Samplig ad Estimatio 1.1 Itroductio This chapter will cover the most importat basic statistical theory you eed i order to uderstad the ecoometric material that will be comig
More informationSAMPLING LIPSCHITZ CONTINUOUS DENSITIES. 1. Introduction
SAMPLING LIPSCHITZ CONTINUOUS DENSITIES OLIVIER BINETTE Abstract. A simple ad efficiet algorithm for geeratig radom variates from the class of Lipschitz cotiuous desities is described. A MatLab implemetatio
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationOPTIMAL ALGORITHMS -- SUPPLEMENTAL NOTES
OPTIMAL ALGORITHMS -- SUPPLEMENTAL NOTES Peter M. Maurer Why Hashig is θ(). As i biary search, hashig assumes that keys are stored i a array which is idexed by a iteger. However, hashig attempts to bypass
More informationMath 155 (Lecture 3)
Math 55 (Lecture 3) September 8, I this lecture, we ll cosider the aswer to oe of the most basic coutig problems i combiatorics Questio How may ways are there to choose a -elemet subset of the set {,,,
More informationKurskod: TAMS11 Provkod: TENB 21 March 2015, 14:00-18:00. English Version (no Swedish Version)
Kurskod: TAMS Provkod: TENB 2 March 205, 4:00-8:00 Examier: Xiagfeg Yag (Tel: 070 2234765). Please aswer i ENGLISH if you ca. a. You are allowed to use: a calculator; formel -och tabellsamlig i matematisk
More informationSequences A sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece 1, 1, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet
More informationThere is no straightforward approach for choosing the warmup period l.
B. Maddah INDE 504 Discrete-Evet Simulatio Output Aalysis () Statistical Aalysis for Steady-State Parameters I a otermiatig simulatio, the iterest is i estimatig the log ru steady state measures of performace.
More informationIntroduction to Signals and Systems, Part V: Lecture Summary
EEL33: Discrete-Time Sigals ad Systems Itroductio to Sigals ad Systems, Part V: Lecture Summary Itroductio to Sigals ad Systems, Part V: Lecture Summary So far we have oly looked at examples of o-recursive
More informationInformation Theory and Coding
Sol. Iformatio Theory ad Codig. The capacity of a bad-limited additive white Gaussia (AWGN) chael is give by C = Wlog 2 ( + σ 2 W ) bits per secod(bps), where W is the chael badwidth, is the average power
More informationAdvanced Stochastic Processes.
Advaced Stochastic Processes. David Gamarik LECTURE 2 Radom variables ad measurable fuctios. Strog Law of Large Numbers (SLLN). Scary stuff cotiued... Outlie of Lecture Radom variables ad measurable fuctios.
More informationIntroduction to Probability. Ariel Yadin
Itroductio to robability Ariel Yadi Lecture 2 *** Ja. 7 ***. Covergece of Radom Variables As i the case of sequeces of umbers, we would like to talk about covergece of radom variables. There are may ways
More informationResampling Methods. X (1/2), i.e., Pr (X i m) = 1/2. We order the data: X (1) X (2) X (n). Define the sample median: ( n.
Jauary 1, 2019 Resamplig Methods Motivatio We have so may estimators with the property θ θ d N 0, σ 2 We ca also write θ a N θ, σ 2 /, where a meas approximately distributed as Oce we have a cosistet estimator
More informationFormation of A Supergain Array and Its Application in Radar
Formatio of A Supergai Array ad ts Applicatio i Radar Tra Cao Quye, Do Trug Kie ad Bach Gia Duog. Research Ceter for Electroic ad Telecommuicatios, College of Techology (Coltech, Vietam atioal Uiversity,
More information( ) = p and P( i = b) = q.
MATH 540 Radom Walks Part 1 A radom walk X is special stochastic process that measures the height (or value) of a particle that radomly moves upward or dowward certai fixed amouts o each uit icremet of
More informationNotes 5 : More on the a.s. convergence of sums
Notes 5 : More o the a.s. covergece of sums Math 733-734: Theory of Probability Lecturer: Sebastie Roch Refereces: Dur0, Sectios.5; Wil9, Sectio 4.7, Shi96, Sectio IV.4, Dur0, Sectio.. Radom series. Three-series
More information4. Partial Sums and the Central Limit Theorem
1 of 10 7/16/2009 6:05 AM Virtual Laboratories > 6. Radom Samples > 1 2 3 4 5 6 7 4. Partial Sums ad the Cetral Limit Theorem The cetral limit theorem ad the law of large umbers are the two fudametal theorems
More informationNUMERICAL METHODS FOR SOLVING EQUATIONS
Mathematics Revisio Guides Numerical Methods for Solvig Equatios Page 1 of 11 M.K. HOME TUITION Mathematics Revisio Guides Level: GCSE Higher Tier NUMERICAL METHODS FOR SOLVING EQUATIONS Versio:. Date:
More informationReliability and Queueing
Copyright 999 Uiversity of Califoria Reliability ad Queueig by David G. Messerschmitt Supplemetary sectio for Uderstadig Networked Applicatios: A First Course, Morga Kaufma, 999. Copyright otice: Permissio
More informationChapter 6 Sampling Distributions
Chapter 6 Samplig Distributios 1 I most experimets, we have more tha oe measuremet for ay give variable, each measuremet beig associated with oe radomly selected a member of a populatio. Hece we eed to
More information5.1 A mutual information bound based on metric entropy
Chapter 5 Global Fao Method I this chapter, we exted the techiques of Chapter 2.4 o Fao s method the local Fao method) to a more global costructio. I particular, we show that, rather tha costructig a local
More information4.3 Growth Rates of Solutions to Recurrences
4.3. GROWTH RATES OF SOLUTIONS TO RECURRENCES 81 4.3 Growth Rates of Solutios to Recurreces 4.3.1 Divide ad Coquer Algorithms Oe of the most basic ad powerful algorithmic techiques is divide ad coquer.
More informationLecture 7: MIMO Architectures Theoretical Foundations of Wireless Communications 1. Overview. CommTh/EES/KTH
: Theoretical Foudatios of Wireless Commuicatios 1 Thursday, May 19, 2016 12:30-15:30, Coferece Room SIP 1 Textbook: D. Tse ad P. Viswaath, Fudametals of Wireless Commuicatio 1 / 1 Overview Lecture 6:
More informationTopic 9: Sampling Distributions of Estimators
Topic 9: Samplig Distributios of Estimators Course 003, 2018 Page 0 Samplig distributios of estimators Sice our estimators are statistics (particular fuctios of radom variables), their distributio ca be
More informationOutput Analysis and Run-Length Control
IEOR E4703: Mote Carlo Simulatio Columbia Uiversity c 2017 by Marti Haugh Output Aalysis ad Ru-Legth Cotrol I these otes we describe how the Cetral Limit Theorem ca be used to costruct approximate (1 α%
More information1 Introduction to reducing variance in Monte Carlo simulations
Copyright c 010 by Karl Sigma 1 Itroductio to reducig variace i Mote Carlo simulatios 11 Review of cofidece itervals for estimatig a mea I statistics, we estimate a ukow mea µ = E(X) of a distributio by
More informationChapter 2 The Monte Carlo Method
Chapter 2 The Mote Carlo Method The Mote Carlo Method stads for a broad class of computatioal algorithms that rely o radom sampligs. It is ofte used i physical ad mathematical problems ad is most useful
More informationLecture Chapter 6: Convergence of Random Sequences
ECE5: Aalysis of Radom Sigals Fall 6 Lecture Chapter 6: Covergece of Radom Sequeces Dr Salim El Rouayheb Scribe: Abhay Ashutosh Doel, Qibo Zhag, Peiwe Tia, Pegzhe Wag, Lu Liu Radom sequece Defiitio A ifiite
More informationMATH 320: Probability and Statistics 9. Estimation and Testing of Parameters. Readings: Pruim, Chapter 4
MATH 30: Probability ad Statistics 9. Estimatio ad Testig of Parameters Estimatio ad Testig of Parameters We have bee dealig situatios i which we have full kowledge of the distributio of a radom variable.
More informationMA131 - Analysis 1. Workbook 3 Sequences II
MA3 - Aalysis Workbook 3 Sequeces II Autum 2004 Cotets 2.8 Coverget Sequeces........................ 2.9 Algebra of Limits......................... 2 2.0 Further Useful Results........................
More informationLecture 4: April 10, 2013
TTIC/CMSC 1150 Mathematical Toolkit Sprig 01 Madhur Tulsiai Lecture 4: April 10, 01 Scribe: Haris Agelidakis 1 Chebyshev s Iequality recap I the previous lecture, we used Chebyshev s iequality to get a
More informationA note on log-concave random graphs
A ote o log-cocave radom graphs Ala Frieze ad Tomasz Tocz Departmet of Mathematical Scieces, Caregie Mello Uiversity, Pittsburgh PA53, USA Jue, 08 Abstract We establish a threshold for the coectivity of
More informationChannel coding, linear block codes, Hamming and cyclic codes Lecture - 8
Digital Commuicatio Chael codig, liear block codes, Hammig ad cyclic codes Lecture - 8 Ir. Muhamad Asial, MSc., PhD Ceter for Iformatio ad Commuicatio Egieerig Research (CICER) Electrical Egieerig Departmet
More informationCooperative Communication Fundamentals & Coding Techniques
3 th ICACT Tutorial Cooperative commuicatio fudametals & codig techiques Cooperative Commuicatio Fudametals & Codig Techiques 0..4 Electroics ad Telecommuicatio Research Istitute Kiug Jug 3 th ICACT Tutorial
More informationFall 2013 MTH431/531 Real analysis Section Notes
Fall 013 MTH431/531 Real aalysis Sectio 8.1-8. Notes Yi Su 013.11.1 1. Defiitio of uiform covergece. We look at a sequece of fuctios f (x) ad study the coverget property. Notice we have two parameters
More informationA proposed discrete distribution for the statistical modeling of
It. Statistical Ist.: Proc. 58th World Statistical Cogress, 0, Dubli (Sessio CPS047) p.5059 A proposed discrete distributio for the statistical modelig of Likert data Kidd, Marti Cetre for Statistical
More informationDefinition 4.2. (a) A sequence {x n } in a Banach space X is a basis for X if. unique scalars a n (x) such that x = n. a n (x) x n. (4.
4. BASES I BAACH SPACES 39 4. BASES I BAACH SPACES Sice a Baach space X is a vector space, it must possess a Hamel, or vector space, basis, i.e., a subset {x γ } γ Γ whose fiite liear spa is all of X ad
More informationA sequence of numbers is a function whose domain is the positive integers. We can see that the sequence
Sequeces A sequece of umbers is a fuctio whose domai is the positive itegers. We ca see that the sequece,, 2, 2, 3, 3,... is a fuctio from the positive itegers whe we write the first sequece elemet as
More informationGeneralized Semi- Markov Processes (GSMP)
Geeralized Semi- Markov Processes (GSMP) Summary Some Defiitios Markov ad Semi-Markov Processes The Poisso Process Properties of the Poisso Process Iterarrival times Memoryless property ad the residual
More informationB. Maddah ENMG 622 ENMG /27/07
B. Maddah ENMG 622 ENMG 5 3/27/7 Queueig Theory () What is a queueig system? A queueig system cosists of servers (resources) that provide service to customers (etities). A Customer requestig service will
More informationCHAPTER 5. Theory and Solution Using Matrix Techniques
A SERIES OF CLASS NOTES FOR 2005-2006 TO INTRODUCE LINEAR AND NONLINEAR PROBLEMS TO ENGINEERS, SCIENTISTS, AND APPLIED MATHEMATICIANS DE CLASS NOTES 3 A COLLECTION OF HANDOUTS ON SYSTEMS OF ORDINARY DIFFERENTIAL
More informationRandom Models. Tusheng Zhang. February 14, 2013
Radom Models Tusheg Zhag February 14, 013 1 Radom Walks Let me describe the model. Radom walks are used to describe the motio of a movig particle (object). Suppose that a particle (object) moves alog the
More informationLecture 2. The Lovász Local Lemma
Staford Uiversity Sprig 208 Math 233A: No-costructive methods i combiatorics Istructor: Ja Vodrák Lecture date: Jauary 0, 208 Origial scribe: Apoorva Khare Lecture 2. The Lovász Local Lemma 2. Itroductio
More informationLecture 9: Hierarchy Theorems
IAS/PCMI Summer Sessio 2000 Clay Mathematics Udergraduate Program Basic Course o Computatioal Complexity Lecture 9: Hierarchy Theorems David Mix Barrigto ad Alexis Maciel July 27, 2000 Most of this lecture
More informationFIR Filters. Lecture #7 Chapter 5. BME 310 Biomedical Computing - J.Schesser
FIR Filters Lecture #7 Chapter 5 8 What Is this Course All About? To Gai a Appreciatio of the Various Types of Sigals ad Systems To Aalyze The Various Types of Systems To Lear the Skills ad Tools eeded
More informationFACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING. Lectures
FACULTY OF MATHEMATICAL STUDIES MATHEMATICS FOR PART I ENGINEERING Lectures MODULE 5 STATISTICS II. Mea ad stadard error of sample data. Biomial distributio. Normal distributio 4. Samplig 5. Cofidece itervals
More informationLecture 2: Concentration Bounds
CSE 52: Desig ad Aalysis of Algorithms I Sprig 206 Lecture 2: Cocetratio Bouds Lecturer: Shaya Oveis Ghara March 30th Scribe: Syuzaa Sargsya Disclaimer: These otes have ot bee subjected to the usual scrutiy
More information7-1. Chapter 4. Part I. Sampling Distributions and Confidence Intervals
7-1 Chapter 4 Part I. Samplig Distributios ad Cofidece Itervals 1 7- Sectio 1. Samplig Distributio 7-3 Usig Statistics Statistical Iferece: Predict ad forecast values of populatio parameters... Test hypotheses
More informationLecture 01: the Central Limit Theorem. 1 Central Limit Theorem for i.i.d. random variables
CSCI-B609: A Theorist s Toolkit, Fall 06 Aug 3 Lecture 0: the Cetral Limit Theorem Lecturer: Yua Zhou Scribe: Yua Xie & Yua Zhou Cetral Limit Theorem for iid radom variables Let us say that we wat to aalyze
More information7. Modern Techniques. Data Encryption Standard (DES)
7. Moder Techiques. Data Ecryptio Stadard (DES) The objective of this chapter is to illustrate the priciples of moder covetioal ecryptio. For this purpose, we focus o the most widely used covetioal ecryptio
More informationCastiel, Supernatural, Season 6, Episode 18
13 Differetial Equatios the aswer to your questio ca best be epressed as a series of partial differetial equatios... Castiel, Superatural, Seaso 6, Episode 18 A differetial equatio is a mathematical equatio
More informationSimulation. Two Rule For Inverting A Distribution Function
Simulatio Two Rule For Ivertig A Distributio Fuctio Rule 1. If F(x) = u is costat o a iterval [x 1, x 2 ), the the uiform value u is mapped oto x 2 through the iversio process. Rule 2. If there is a jump
More informationThe Boolean Ring of Intervals
MATH 532 Lebesgue Measure Dr. Neal, WKU We ow shall apply the results obtaied about outer measure to the legth measure o the real lie. Throughout, our space X will be the set of real umbers R. Whe ecessary,
More informationA Question. Output Analysis. Example. What Are We Doing Wrong? Result from throwing a die. Let X be the random variable
A Questio Output Aalysis Let X be the radom variable Result from throwig a die 5.. Questio: What is E (X? Would you throw just oce ad take the result as your aswer? Itroductio to Simulatio WS/ - L 7 /
More informationComparing your lab results with the others by one-way ANOVA
Comparig your lab results with the others by oe-way ANOVA You may have developed a ew test method ad i your method validatio process you would like to check the method s ruggedess by coductig a simple
More informationlim za n n = z lim a n n.
Lecture 6 Sequeces ad Series Defiitio 1 By a sequece i a set A, we mea a mappig f : N A. It is customary to deote a sequece f by {s } where, s := f(). A sequece {z } of (complex) umbers is said to be coverget
More informationJoint Probability Distributions and Random Samples. Jointly Distributed Random Variables. Chapter { }
UCLA STAT A Applied Probability & Statistics for Egieers Istructor: Ivo Diov, Asst. Prof. I Statistics ad Neurology Teachig Assistat: Neda Farziia, UCLA Statistics Uiversity of Califoria, Los Ageles, Sprig
More informationChapter 4. Fourier Series
Chapter 4. Fourier Series At this poit we are ready to ow cosider the caoical equatios. Cosider, for eample the heat equatio u t = u, < (4.) subject to u(, ) = si, u(, t) = u(, t) =. (4.) Here,
More informationChapter 11 Output Analysis for a Single Model. Banks, Carson, Nelson & Nicol Discrete-Event System Simulation
Chapter Output Aalysis for a Sigle Model Baks, Carso, Nelso & Nicol Discrete-Evet System Simulatio Error Estimatio If {,, } are ot statistically idepedet, the S / is a biased estimator of the true variace.
More informationInformation Theory Tutorial Communication over Channels with memory. Chi Zhang Department of Electrical Engineering University of Notre Dame
Iformatio Theory Tutorial Commuicatio over Chaels with memory Chi Zhag Departmet of Electrical Egieerig Uiversity of Notre Dame Abstract A geeral capacity formula C = sup I(; Y ), which is correct for
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.436J/15.085J Fall 2008 Lecture 19 11/17/2008 LAWS OF LARGE NUMBERS II THE STRONG LAW OF LARGE NUMBERS
MASSACHUSTTS INSTITUT OF TCHNOLOGY 6.436J/5.085J Fall 2008 Lecture 9 /7/2008 LAWS OF LARG NUMBRS II Cotets. The strog law of large umbers 2. The Cheroff boud TH STRONG LAW OF LARG NUMBRS While the weak
More information