On Spatial Capacity of Wireless Ad Hoc Networks with Threshold Based Scheduling

Size: px
Start display at page:

Download "On Spatial Capacity of Wireless Ad Hoc Networks with Threshold Based Scheduling"

Transcription

1 On Spatal Capacty of Wreless Ad Hoc Networks wth Threshold Based Schedulng Yue Lng Che, Ru Zhang, Y Gong, and Lngje Duan Abstract arxv: v [cs.it] 9 Sep 24 Ths paper studes spatal capacty n a stochastc wreless ad hoc network, where mult-stage probng and data transmsson are sequentally performed. We propose a novel sgnal-to-nterference-rato (SIR) threshold based schedulng scheme: by startng wth ntal probng, each transmtter teratvely decdes to further probe or stay dle, dependng on whether the estmated SIR n the proceedng probng s larger or smaller than a predefned threshold. Snce only local SIR nformaton s requred for makng transmsson decson, the proposed scheme s approprate for dstrbuted mplementaton n practcal wreless ad hoc networks. Although one can assume that the transmtters are ntally deployed accordng to a homogeneous Posson pont process (PPP), the SIR based schedulng makes the PPP no longer applcable to model the locatons of retaned transmtters n the subsequent probng and data transmsson phases, due to the nterference nduced couplng n ther decsons. As the analyss becomes very complcated, we frst focus on sngle-stage probng and fnd that when the SIR threshold s set suffcently small to assure an acceptable nterference level n the network, the proposed scheme can greatly outperform the nonschedulng reference scheme n terms of spatal capacty. We clearly characterze the spatal capacty and obtan exact/approxmate closed-form expressons, by proposng a new approxmate approach to deal wth the correlated SIR dstrbutons over non-posson pont processes. Then we successfully extend to mult-stage probng by properly desgnng the multple SIR thresholds to assure gradual mprovement of the spatal capacty. Furthermore, we analyze the mpact of mult-stage probng overhead and present a probng-capacty tradeoff n schedulng desgn. Fnally, extensve numercal results are presented to demonstrate the performance of the proposed schedulng as compared to exstng schemes. Index Terms Wreless ad hoc network, threshold based schedulng, spatal capacty, stochastc geometry. I. INTRODUCTION Wreless ad hoc networks have emerged as a promsng technology that can provde seamless communcaton between wreless users (transmtter-recever pars) wthout relyng on any pre-exstng nfrastructure. In such Ths work has been presented n part at IEEE Internatonal Symposum on Informaton Theory (ISIT), Istanbul, Turkey, July 7-2, 23. Y. L. Che and Y. Gong are wth the School of Electrcal and Electronc Engneerng, Nanyang Technologcal Unversty, Sngapore (e-mal: chey4@ntu.edu.sg; eygong@ntu.edu.sg). R. Zhang s wth the Department of Electrcal and Computer Engneerng, Natonal Unversty of Sngapore (e-mal: elezhang@nus.edu.sg). He s also wth the Insttute for Infocomm Research, A*STAR, Sngapore. L. Duan s wth the Engneerng Systems and Desgn Pllar, Sngapore Unversty of Technology and Desgn (e-mal: lngje duan@sutd.edu.sg).

2 2 networks, the wreless users communcate wth each other n a dstrbuted manner. Due to the lack of centralzed coordnators to coordnate the transmssons among the users, the wreless ad hoc network s under compettve and nterference-domnant envronment n nature. Thereby, effcent transmsson schemes for transmtters to effectvely schedule/adapt ther transmssons are appealng for system performance mprovement, and thus have attracted wde research attentons n the past decade. Tradtonally, each transmtter s enabled to ndependently decde whether to transmt over a partcular channel based on ts own wllngness or channel strength []-[4], and the transmsson rate of each user can be maxmzed by fndng an optmal transmsson probablty or an optmal channel strength threshold, respectvely. Although easy to be mplemented, such ndependent transmsson schemes do not consder the resultng user nteractons n the wreless ad hoc networks due to the co-channel nterference, and thus do not acheve hgh system performance n general cases. Therefore, more complex transmsson schemes have been proposed to explot the user nteractons by explorng the nformaton of sgnal-to-nterference-rato (SIR). For example, by teratvely adaptng the transmt power level based on the estmated SIR, the Foschn-Mljanc algorthm [5] assures zero outage probablty and/or mnmum aggregate power consumpton for uplnk transmsson n a cellular network. In [6], Yates has studed power convergence condtons for such teratve power control algorthms. Moreover, there have been some recent studes (e.g. [7] and [8]) that extend the Foschn-Mljanc algorthm to the wreless ad hoc network through jont schedulng and power control transmsson schemes. In addton, by adaptng the transmsson probablty dependng on the receved SIR, [9] has studed varous random access schemes to mprove the system throughput and/or the user farness. However, [5]-[9] ether requre each transmtter to know at least the wreless envronment nformaton of ts neghbors, or are of hgh mplementaton complexty, and thus are not approprate for practcal large-scale wreless ad hoc networks. On the other hand, due to the randomzed locaton of each transmtter and the effects of channel fadng, the network-level performance analyss s fundamentally mportant for the study of wreless ad hoc networks. It s noted that Gupta and Kumar n [] studed scalng laws, whch quantfed the ncrease of the volume of capacty regon over the number of transmtters n ad hoc networks. Moreover, to determne the set of actve transmtters that can yeld maxmum aggregate Shannon capacty n the network, the authors n []-[3] addressed the capacty maxmzaton problem for an arbtrary wreless ad hoc network. However, []-[3] dd not consder the mpact of spatal confguraton of the ad hoc network, whch s a crtcal factor that determnes the ad hoc network capacty [4]. It came to our attenton that as a powerful tool to capture the mpact of wreless users spatal randomness on the network performance, stochastc geometry [5] s able to provde more comprehensve

3 3 characterzaton of the performance of wreless networks, and thus has attracted great attentons from both academy and ndustry [4], [6]. Among all the tools provded by stochastc geometry, homogeneous Posson pont process (PPP) [7] s the most wdely used one for network topology modelng and performance analyss. Under the assumpton that the transmtters are deployed accordng to a homogeneous PPP, the exact/approxmate capacty of a wreless ad hoc network under varous ndependent transmsson schemes, such as Aloha-based random transmsson [], channel-nverson based power control [2], and channel-threshold based schedulng [3], [4], can all be successfully characterzed by usng advanced tools from stochastc geometry. However, lmted work based on stochastc geometry has studed SIR-based transmsson schemes, where the user nteractons are nvolved. It s noted that [8] studed a probablty-based schedulng scheme, where each transmtter ndependently adjusts ts current transmsson probablty based on the receved SIR n the proceedng teraton. However, [8] only studed the convergence of the probablty-based schedulng, wthout addressng the network capacty wth spatalty dstrbuted users. To our best knowledge, there has been no exstng work on studyng the wreless ad hoc network capacty wth a SIR-based transmsson scheme. Hence, the mpact of SIR-based transmssons s lmtedly understood from the network-level pont of vew. A prncple goal of ths study s to use stochastc geometry to fll the vod of wreless network capacty characterzaton by an effcent SIR-based transmsson scheme. To ths end, we propose a novel SIR-threshold based schedulng scheme for a sngle-hop slotted wreless ad hoc network. We consder a probe-and-transmt protocol, where mult-stage probngs are sequentally performed to gradually determne the transmtters that are allowed to transmt data n each slot. Specfcally, we assume there are n total N probng phases and one data transmsson phase n each slot, N <. We sequentally label the N probng phases as P-Phase, P-Phase,..., and P-PhaseN, and label the data transmsson phase as D-Phase. As llustrated n Fg., f the feedback SIR from recever n P-Phase k, k N, s no smaller than a pre-defned threshold, transmtter decdes to transmt n P-Phase k; otherwse, to mprove the system throughput as well as save ts own energy, transmtter stays dle n the remanng tme of the slot as n [9], so as to let other transmtters that have hgher SIR levels re-contend the current transmsson opportunty. Snce each transmtter only requres drect-channel SIR feedback from ts ntended recever for lmted tmes, the proposed scheme s approprate for dstrbuted mplementaton n practcal wreless ad hoc networks. In ths paper, we characterze the wreless ad hoc network capacty wth a metrc called spatal capacty, whch has been used n [2] and gves the average number of successful transmtters per unt area for any gven ntal transmtter densty. We am at closed-form spatal capacty characterzaton and maxmzaton by explorng the SIR-threshold based transmsson.

4 4... Transmtter : P-Phase... Transmsson decson P-Phase k tme tme P-Phase P-Phase k-... Idle... tme Fg.. Illustraton of the SIR-threshold based transmsson n P-Phase k: If the SIR n P-Phase k s no smaller than the threshold γ k n P-Phase k, the transmtter decdes to transmt n P-Phase k; otherwse, t stays dle n the remanng tme of ths slot. The key contrbutons of ths paper are summarzed as follows. Novel SIR-threshold based schedulng scheme: In Secton II, we propose a novel SIR-threshold based transmsson scheme for a sngle-hop wreless ad hoc network, whch can be mplemented effcently n a dstrbuted manner. Though one can use a homogeneous PPP to model the stochastc locatons of the transmtters n the ntal probng phase, we fnd that due to the teratve SIR-based schedulng, the PPP model s no longer applcable to model the locatons of the retaned transmtters n all the subsequent probng or data transmsson phases. Furthermore, snce the SIR dstrbutons n all the probng and data transmsson phases are strongly correlated, t s challengng to analyze/characterze the spatal capacty of the proposed scheme. Sngle-stage probng for spacal capacty mprovement: In Secton III, we start up wth sngle-stage probng (N = ) to clearly decde the SIR threshold for the proposed scheme and characterze the spatal capacty. We show that a small SIR threshold can effcently reduce the retaned transmtter number and thus the nterference level n the data transmsson phase, whle a large SIR threshold wll overly reduce the retaned transmtter number and does not help mprove the spatal capacty. We also propose a new approxmate approach to characterze the spatal capacty n closed-form, whch s useful for analyzng performance of wreless networks wth nteracted transmtters. Mult-stage probng for spatal capacty mprovement: In Secton IV, we extend proposed schedulng scheme from the sngle-stage probng (N = ) to mult-stage probng (N > ) for greater spatal capacty mprovement. We show that once a sequence of ncreasng SIR thresholds are properly decded over probng phases, the spatal capacty s assured to gradually mprove. As mult-stage probng can ntroduce non-gnorable overhead n each tme slot, whch reduces the spatal capacty, we study an nterestng probng-capacty tradeoff over the probng-stage number N. Performance evaluatons for network desgn: In both Secton III and Secton IV, we also provde extensve

5 5 numercal results to further evaluate the mpact of key parameters of the proposed scheme. In partcular, we present a densty-capacty tradeoff n Secton III-C-), whch shows that a small ntal transmtter densty can help mprove the spatal capacty, whle a large one wll ntroduce hgh nterference level and thus reduce the spatal capacty. To hghlght the spatal capacty mprovement performance of the proposed scheme, we also compare the proposed scheme wth exstng dstrbuted schedulng schemes n Secton III-C-2). Moreover, we consder a practcal scenaro wth SIR estmaton and feedback errors and show that the proposed desgn s robust to the SIR errors n Secton III-C-3) by smulaton. In Secton IV, we study an example wth N = 2 and show the correspondng spatal capacty over both SIR thresholds n P-Phase and D-Phase. Interestngly, our numercal results show that the former SIR threshold plays a more crtcal role n determnng the spatal capacty than the latter one, snce the former SIR threshold decdes how many transmtters can have a second chance to contend the transmsson opportunty. It s noted that some of the exstng work has addressed the throughput/capacty analyss of a wreless communcaton system from the nformaton-theoretc pont of vew. For example, Tse and Hanly consdered a multpontto-pont system and characterzed the throughput capacty regon and delay-lmted capacty regon of the fadng multple-access channel n [2] and [22], respectvely, where the optmal power and/or rate allocaton that can acheve the boundary of the capacty regons was derved. Although appealng, both [2] and [22] have assumed multuser detecton at a centralzed recever and gnored the mpact of the random network topology drven by moble transmtters and recevers moblty, and thus cannot completely provde network-level system performance characterzaton wth dstrbuted sngle-user detecton (.e., treatng the multuser nterference as nose) recevers. Unlke Tse and Hanly s works n [2] and [22], we use stochastc geometry to model the large-scale random wreless ad hoc network topology, and novelly analyze the network-level performance of the teratve SIR-threshold based schedulng. In addton, t s also noted that some exstng work has adopted tools from stochastc geometry to study the non- PPP based wreless network. For example, by usng a PPP to approxmate the underlyng non-ppp based spatal dstrbuton of the transmtters locatons, [23]-[27] have successfully characterzed the non-ppp based wreless network capacty. Unlke [23]-[27], due to the teratve SIR-based schedulng of the proposed scheme, we need to address not only the non-ppp based spatal dstrbuton of the transmtters locatons, but also the resultng strongly-correlated SIR dstrbutons over all probng and data transmsson phases. To our best knowledge, such correlated SIR analyss/chracterzaton n non-ppp based wreless networks has not been addressed n the exstng work based on stochastc geometry.

6 6 II. SYSTEM MODEL AND PERFORMANCE METRIC In ths secton, we descrbe the consdered transmsson schemes n ths paper. We then develop the network model based on stochastc geometry. At last, we defne the spatal capacty as our performance metrc. A. Transmsson Schemes We focus on the proposed scheme wth SIR-threshold based schedulng. For comparson, we also consder a reference scheme wthout any transmsson schedulng. For both transmsson schemes, we assume that all transmtters transmt n a synchronzed tme-slotted manner. We also assume that all transmtters transmt at the same power level, whch s normalzed to be unty for convenence. ) SIR-Threshold Based Scheme: Based on the probe-and-transmt protocol, n each tme slot, N probng phases wth N < are sequentally mplemented before the data transmsson phase. We assume N s a pre-gven parameter and ts effects wll be studed later n Secton IV-B. Moreover, as shown n Fg., we denote the duraton of a tme slot and a probng phase as T and τ, respectvely, wth τ T, such that Nτ < T, as n [9]. By normalzng over T, the effectve data transmsson tme n a tme slot s obtaned as T Nτ T, whch reduces lnearly over N [28]. Furthermore, we assume f a transmtter transmts probng sgnals n a probng phase, ts ntended recever s able to measure the receved sgnal power over the total nterference power,.e., the SIR, and feeds t back to the transmtter at the end of the probng phase. The specfc algorthm desgn on SIR estmaton and feedback s out of the scope of ths paper and s not our focus. To obtan tractable analyss, we assume perfect SIR estmaton and feedback n ths paper, and thus the SIR value s exactly known at the transmtter; however, the mpact of fnte SIR estmaton and feedback errors on the network capacty s mportant to practcal desgn and thus wll also be evaluated by smulaton. Accordng to the feedback SIR level of ts own channel, each transmtter teratvely performs the threshold-based transmsson decson n each P-Phase or D-Phase, for whch the detals are gven as follows: In the ntal probng phase,.e., P-Phase, to ntalze the communcaton between each transmtter and recever par, all transmtters ndependently transmt probng sgnals to ther ntended recevers. Each recever then estmates the channel ampltude and phase (for possble coherent communcaton n the subsequent probng and data transmsson phases), and measures the receved SIR of the probng sgnal. Each transmtter receves the feedback SIR from ts ntended recever at the end of P-Phase. In general, each transmtters can transmt at dfferent power levels by teratvely adjustng ts transmt power based on the feedback SIR nformaton, as n [5] or [6]. However, n ths paper, we manly focus on SIR-based transmsson schedulng and thus restrct transmt power adaptaton to be bnary for smplcty.

7 7 In each of the remanng probng phases from P-Phase to P-Phase N, by explotng the feedback SIR n the proceedng probng phase, each transmtter decdes whether to transmt n the current probng phase wth a predefned SIR threshold. Specfcally, suppose a transmtter transmts n P-Phase k, k N. As shown n Fg., f the feedback SIR n P-Phase k s larger than or equal to the predefned SIRthreshold, denoted by γ k for P-Phase k, the transmtter contnues ts transmsson n P-Phase k and thus receves the feedback SIR n P-Phase k; otherwse, to mprove the system throughput as well as save ts energy, the transmtter decdes not to transmt any more n the remanng tme of ths slot and wll seek another transmsson opportunty n the next slot, so as to let other transmtters that have hgher SIR levels to re-contend the current transmsson opportunty. In the D-Phase, smlar to the SIR-threshold based schedulng from P-Phaseto P-PhaseN, f a transmtter transmts n P-Phase N and ts feedback SIR n P-Phase N s larger than or equal to the predefned threshold, denoted by γ N for the D-Phase, the transmtter sends data to ts ntended recever; otherwse, the transmtter remans slent n the rest tme of ths slot. The data transmsson s successful f the SIR at the recever s larger than or equal to the requred SIR level, denoted by β >. 2) Reference Scheme: There s no transmsson schedulng n the reference scheme. In each tme slot, we assume all transmtters transmt data drectly to ther ntended recevers n an ndependent manner. Thus, the effectve data transmsson tme for the reference scheme s. 2 The data transmsson s successful f the SIR at the recever s larger than or equal to the requred SIR level β as the proposed scheme. Note that by mplementng an ntal probng phase before the data transmsson, the reference scheme can be mproved to be a proposed scheme wth sngle-stage probng. B. Network Model In the next, we develop the network model based on stochastc geometry. For both consdered transmsson schemes, we focus on sngle-hop communcaton n one partcular tme slot. For both schemes, we assume that all transmtters are ndependently and unformly dstrbuted n the unbounded two-dmensonal plane R 2. We thus model the locatons of all the transmtters by a homogeneous PPP wth densty λ. Due to the lack of central nfrastructure for coordnaton n the wreless ad hoc network, we assume the transmtters have no knowledge about ther surroundng wreless envronment, and thus ntend to transmt ndependently n a tme slot wth probablty θ (,), as n []-[4]. Denote λ = λθ as the densty of the ntal 2 It s worth pontng out that for the reference scheme, an ntal tranng s needed pror to data transmsson for the recever to estmate the channel for coherent communcaton, smlar to the ntal probng of the proposed scheme wth N =, but wthout the SIR feedback to the transmtter. Here, we have assumed that such tranng ncurs a neglgble tme overhead as compared to each slot duraton.

8 8 transmtters that have the ntenton to transmt n a partcular tme slot. Accordng to the Colorng theory [5], the process of the ntal transmtters for both schemes s a homogeneous PPP wth densty λ, whch s denoted by Φ. Wthout loss of generalty, we assume λ and θ and hence λ are gven parameters, and wll dscuss the effects of λ later n Secton III-C. We assume each transmtter has one ntended recever, whch s unformly dstrbuted on a crcle of radus d meters (m) centered at the transmtter. We denote the locatons of the -th transmtter and ts ntended recever as x, wth x Φ, and r (not ncluded n Φ ), respectvely. The path loss between the -th transmtter and the j-th recever s gven by l j = x r j α, where α>2 s the path-loss exponent. We use h j to denote the dstance-ndependent channel fadng coeffcent from transmtter to recever j. We assume flat Raylegh fadng, where all h j s are ndependent and exponentally dstrbuted random varables wth unt mean. We also assume that h j s do not change wthn one tme-slot. We denote the SIR at the -th recever as SIR (), whch s gven by SIR () = h d α x j Φ,j h jl j. () Note that for the reference scheme wthout transmsson schedulng, SIR () gves the receved SIR level at the -th recever for the data transmsson of transmtter. As a result, n the reference scheme, the data transmsson of transmtter s successful f SIR () β s satsfed. Unlke the reference scheme, n the proposed scheme, SIR () only gves the receved SIR level at the -th recever n the ntal probng phase P-Phase. We then denote the pont process formed by the retaned transmtters n P-Phase k wth k N, or the D-Phase wth k = N, as Φ k. We also denote SIR (k) the -th recever n Φ k. Clearly, we have Φ k = {x Φ k : SIR (k ) as the receved SIR at γ k }, where the number of transmtters n Φ k s reduced as compared to that n Φ k. Thus, t s easy to verfy that SIR (k) SIR (k ) for any gven γ k, Φ k Φ k. Moreover, smlar to SIR (), for any Φ k, k {,...,N}, we can express SIR (k) as SIR (k) = h d α x j Φ k,j h jl j, k {,...,N}. (2) It s worth notng that due to the SIR-based schedulng, the transmtters are not retaned ndependently n Φ k. Thus, unlke SIR () n (), whch s determned by the homogeneous PPP Φ, SIR (k) n (2) s determned by the non-ppp Φ k n general [5]. For the proposed scheme, the data transmsson of transmtter s successful f SIR (k ) γ k, k {,...,N}, and SIR (N) β are all satsfed. C. Spatal Capacty Due to the statonarty of the homogeneous PPP Φ, t s easy to verfy that Φ k, k {,...N}, s also statonary [4]. We thus consder a typcal par of transmtter and recever n ths paper. Wthout loss of generalty, we assume

9 9 that the typcal recever s located at the orgn. The typcal par of transmtter and recever s named par,.e., =. Denote the successful transmsson probablty of the typcal par n the data transmsson phase of the proposed scheme wth N probng phases or the reference scheme as P p,n or P r, respectvely. We thus have P p,n = P(SIR () γ,..., SIR (N ) γ N, SIR (N) β). (3) P r = P(SIR () β). (4) We adopt spatal capacty as our performance metrc, whch s defned as the spatal densty of successful transmssons, or more specfcally the average number of transmtters wth successful data transmsson per unt area. Consderng the effectve data transmsson tme n a tme slot, we thus defne the spatal capacty by the proposed scheme wth N probng phases and the reference scheme as C p,n and C r, respectvely, gven by C p,n T Nτ λ P p,n, (5) T C r λ P r. (6) For the reference scheme, t s noted that P r, gven n (4), s the complementary cumulatve dstrbuton functon (CCDF) of SIR () taken at the value of β. We then have the followng proposton. Proposton 2.: The successful transmsson probablty n the reference scheme s where ρ = +v α/2 dv. When α = 4, we have ρ = π 2. P r = exp( πλ d 2 β 2 αρ), (7) The proof of Proposton 2. s smlar to that of [3, Theorem 2], whch s based on the probablty generatng functonal (PGFL) of the PPP, and thus s omtted here. Snce the network nterference level n the D-Phase ncreases over the ntal transmtter densty λ, we fnd that P r n (7) monotoncally decreases over λ as expected. Moreover, from (6) and (7), we can obtan the expresson of C r as C r = λ exp( πλ d 2 β 2 α ρ). (8) It s observed from (8) that unlke P r, the spatal capacty Cr does not vary monotoncally over λ, snce C r can be benefted by ncreasng λ f the resultng nterference s acceptable. Moreover, from (7) and (8), t s also expected that both P r and Cr monotoncally decrease over the dstance d between each transmtter and recever par, due to the reduced sgnal power receved at the recever, and decrease over the requred SIR level β. Unlke the reference scheme, whch s determned by the homogeneous PPP Φ, the proposed scheme s jontly determned by Φ and a sequence of non-ppps {Φ k }, k N, where the resultng SIR dstrbutons are

10 correlated. Therefore, t s very dffcult to analyze/characterze the spatal capacty of the proposed scheme wth N probng phases. To start up, n the next secton, we focus on a smple case wth sngle-stage probng (N = ) for some nsghtful results. III. SIR-THRESHOLD BASED SCHEME WITH SINGLE-STAGE PROBING In ths secton, we consder the proposed scheme wth sngle-stage probng,.e., N =. In ths case, there s only one round of SIR-based schedulng, whch s mplemented wth the threshold γ. For notatonal smplcty, for the case of N =, we omt the superscrpt N and use P p and Cp to represent the successful transmsson probablty and the spatal capacty of the typcal transmtter, respectvely. Based on (3), the successful transmsson probablty for the case of N = s reduced to P p = P(SIR() γ, SIR () β) (9) = P(SIR () γ )P(SIR () β SIR () γ ). () Moreover, when N =, the effectve data transmsson tme for the proposed scheme s T τ T. Snce τ T, we assume the sngle-stage probng overhead s neglgble; and thus, the effectve data transmsson tme becomes as the reference scheme. Consequently, based on (5), we can express the spatal capacty C p as Furthermore, by substtutng () to (), we can express C p alternatvely as C p =λ P p. () C p =λ P(SIR () γ )P(SIR () β SIR () γ ) =λ P(SIR () β SIR () γ ) (2) where λ =λ P(SIR () γ ) s the densty of Φ n the D-Phase, wth λ λ. Based on Proposton 2., by replacng β wth γ, t s easy to fnd that λ = λ exp ( πλ d 2 γ 2 α ρ ). (3) In the followng two subsectons, we compare the spatal capacty of the two consdered schemes, and characterze C p for the proposed scheme. A. Spatal Capacty Comparson and Closed-form Characterzaton wth γ = and γ β In ths subsecton, we compare the spatal capacty of the proposed scheme wth that of the reference scheme. We then characterze the spatal capacty C p for the proposed scheme and obtan closed-form expressons for the cases of γ = and γ β.

11 Frst, from (6) and (), to compare C p and C r, the key s to compare P p and Pr. In the reference scheme, denote the total nterference power receved at the typcal recever as I = x Φ, h l. In the proposed scheme, the receved total nterference power at the typcal recever n P-Phase s thus I, whle that n the D-Phase s gven by I = x Φ, h l. For any γ, we have I I snce Φ Φ, and thus SIR () SIR (). As a result, by changng over the value of γ [, ), we obtan the followng proposton. Proposton 3.: Gven the requred SIR level β>, for any γ [, ), we have C p > C r, f < γ < β ( conservatve transmsson regme ) C p = C r, f γ = or γ = β ( neutral transmsson regme ) C p < C r, f γ > β ( aggressve transmsson regme ). Proof: Please refer to Appendx A. (4) Remark 3.: Compared to the spatal capacty of the reference scheme, Proposton 3. shows that for the proposed scheme wth SIR-threshold based schedulng, due to the reduced nterference level n the D-Phase, the spatal capacty s mproved n the conservatve transmsson regme wth <γ < β. However, n the case of the aggressve transmsson regme wth γ >β, where the transmtters that are able to transmt successfully n the D-Phase may also be removed from transmsson, the retaned transmtters n the D-Phase are overly reduced. Consequently, the spatal capacty s reduced n the aggressve transmsson regme. It s also noted that n the neutral transmsson regme wth γ = or γ =β, the spatal capacty s dentcal for the two schemes. At last, t s worth notng that Proposton 3. holds regardless of the specfc channel fadng dstrbuton and/or transmtter locaton dstrbuton. Next, we characterze the spatal capacty C p for the proposed scheme wth N =. We focus on dervng the successful transmsson probablty P p n (9). Unlke Pr taken at value β, P p n (4), whch s gven by the margnal CCDF of SIR() s gven by the jont CCDF of SIR() and SIR () taken at values (γ,β). In the followng, we consder three cases γ =, γ β, and < γ < β, and fnd closed-form spatal capacty expressons for both cases of γ = and γ β. Specfcally, for the smple case wth γ =, we can nfer from Proposton 3. drectly that C p = C r, whch s gven n (8). For the case of γ β, snce SIR () > SIR (), we have P(SIR() β SIR () γ ) =. Accordng to (), we thus obtan P p = P(SIR() γ ) n ths case. By replacng β wth γ n Proposton 2., we further obtan that P p = exp( πλ d 2 γ 2 α ρ). As a result, based on (), we can express C p for the case of γ β as C p = λ exp( πλ d 2 γ 2 α ρ). (5) Smlar to C r gven n (8), t s observed that C p for both cases of γ = and γ β does not vary monotoncally over λ, but monotoncally decreases over the dstance d between each transmtter and recever par. Moreover,

12 2 unlke C r and C p for γ =, C p for γ β s not related to the requred SIR level β any more, snce all the retaned transmtters n the D-Phase meet the condton SIR () β n ths case. However, for the case of <γ <β, P p cannot be smply expressed by a margnal CCDF of SIR() as n the above two cases. Moreover, from (9), due to the correlaton between SIR () and SIR () as well as the underlyng non-ppp Φ that determnes SIR (), t s very dffcult, f not mpossble, to fnd an exact expresson of Pp and thus C p n ths case. As a result, n the next subsecton, we focus on fndng a tght approxmate to C p wth a tractable expresson for the case of <γ <β. B. Approxmate Approaches for Spatal Capacty Characterzaton wth < γ < β Ths subsecton focuses on approxmatng the spatal capacty of the proposed scheme for the case of <γ <β. We frst propose a new approxmate approach for C p and obtan an ntegral-based expresson. Next, to fnd a closed-form expresson for C p, we further approxmate the ntegral-based expresson obtaned by the proposed approach. At last, we apply the conventonal approxmate approach n the lterature and dscuss ts approxmate performance. The detals of the three approxmate approaches are gven as follows. ) Proposed Approxmaton: From (9), to fnd a good approxmate to P p and thus Cp, the key s to fnd a good approxmate to the jont SIR dstrbutons n Φ and Φ. Snce Φ Φ, we frst dvde the ntal PPP Φ nto two dsjont non-ppps: one s Φ, and the other s ts complementary set Φ c = Φ Φ, whch s the pont process formed by the non-retaned transmtters n the D-Phase. We denote the densty of Φ c as λ c =λ λ. Clearly, Φ and Φ c are mutually dependent. Denote the receved SIR level at the typcal recever n Φc as SIR (,c) = h d α / Φ h c l. SnceΦ Φ c =Φ andφ Φ c =, we have/sir() = /( SIR () ) +SIR(,c). As a result, (9) can be equally represented by usng the jont dstrbutons of SIR () and SIR (,c). Next, we state an assumpton, based on whch we can use a homogeneous PPP to approxmate Φ and Φ c, respectvely, such that the exstng results on PPP nterference dstrbuton n the lterature can be appled to approxmate the jont dstrbutons of SIR () and SIR (,c). Assumpton : In the proposed scheme wth N =, the transmtters are retaned ndependently n the D-Phase, wth probablty P(SIR () γ ). By applyng Assumpton, we denote the resultng pont processes formed by the retaned and non-retaned transmtters n the D-Phase as ˆΦ and ˆΦ c, respectvely. Clearly, both ˆΦ and ˆΦ c are homogeneous PPPs. Moreover, the densty of ˆΦ or ˆΦ c s the same as that of Φ or Φ c, respectvely. Snce the two homogeneous PPPs ˆΦ and ˆΦ c are dsjont, they are ndependent of each other [5]. Denote Î= ˆΦ h l and Îc = ˆΦ c h l as the receved nterference power at the typcal recever n ˆΦ and ˆΦ c, respectvely. We then use f Î (x ) and fîc (x 2 ) to

13 3 denote the probablty densty functons (pdfs) of Î and Îc, respectvely. The followng lemma gves the general nterference pdf n a homogeneous PPP-based network wth Raylegh fadng channels, whch s a well-known result n the lterature (e.g., [3]). Lemma 3.: For any homogeneous PPP of densty λ, f the channel fadng s Raylegh dstrbuted, the pdf of the receved nterference I at the typcal recever s gven by f I (x)= πx = ( ) + Γ(+2/α)sn(2π/α)! Moreover, when α = 4, (6) can be further expressed n a smpler closed-form as f I (x) = λ 4 ( π x) 3/2exp ( π4 λ 2 6x ( λπ 2 ) 2/α x 2/α. (6) sn(2π/α) ). (7) As a result, based on Lemma 3., by substtutng λ=λ to (6) and (7), we can obtan fî (x ) for the cases of general α and α=4, respectvely. Smlarly, wth λ=λ c, from (6) and (7) we can obtan f Î c (x 2 ) for general α and α = 4, respectvely. Therefore, by approxmatng Φ and Φ c by ˆΦ and ˆΦ c, respectvely, we can easly approxmate the jont dstrbuton of SIR () and SIR (,c) based on the nterference pdfs fî (x ) and fîc (x 2 ), and thereby obtan an ntegral-based approxmate to P p n the followng proposton. Proposton 3.2: The successful transmsson probablty by the proposed scheme for the case of < γ < β s approxmated as h P p e h βd α Proof: Please refer to Appendx B. fî (x ) h γ d α x fîc (x 2 )dx 2 dx dh. (8) Fnally, by multplyng λ wth the rght-hand sde of (8), we obtan an ntegral-based approxmate to C p for the case of < γ < β as h C p λ e h βd α fî (x ) h γ d α x fîc (x 2 )dx 2 dx dh. (9) Note that the proposed approxmate approach consders the correlaton between SIR () and SIR (), and only adopts PPP-based approxmaton to approxmate Φ and Φ c by ˆΦ and ˆΦ c, respectvely. Snce t has been shown n the lterature (e.g., [23]-[26]) that such PPP-based approxmaton can provde tght approxmate to the correspondng non-ppp, the proposed approxmate approach s able to provde tght spatal capacty approxmate to C p for the case of < γ < β. 2) Closed-form Approxmaton for (9): Although the spatal capacty expresson obtaned n (9) s easy to ntegrate, t s not of closed-form. Thus, based on (8), we focus on fndng a closed-form approxmate to P p and thus C p. We frst ncrease the upper lmt of fîc (x 2 ) n (8) from γ d α x to γ d α to obtan an upper bound

14 4 for the rght-hand sde of (8). Then by properly lower-boundng the obtaned upper bound based on Chebyshev s nequalty [32], we obtan a closed-form approxmate to P p, whch s shown n the followng proposton. Proposton 3.3: Based on the ntegral-based expresson gven n (8), a closed-form approxmate to P p for the case of < γ < β s obtaned as P p exp( πλ d 2 β 2 α ρ)exp( πλ c d 2 γ 2 α ρ). (2) Proof: Please refer to Appendx C. From (), (3) and (2), we obtan a closed-form approxmate to spatal capacty of the proposed scheme for the case of < γ < β as C p λ exp ( πλ d 2 γ 2 α ρ ) exp [ πλ exp( πλ d 2 γ 2 α ρ)d 2 ] β 2 α ρ exp [ πλ exp( πλ d 2 γ 2 α ρ)d 2 γ 2 α ρ ]. (2) 3) Conventonal Approxmaton: It s noted that the conventonal approxmate approach n the lterature (e.g., [23]-[26]), whch only focuses on dealng wth the non-ppp Φ, can often yeld a closed-form expresson. Thus, n the followng, we apply the conventonal approxmate approach and dscuss ts approxmate performance to C p. Frst, snce only the performance n Φ s concerned by the conventonal approxmate approach, t takes P ( SIR () β ) as the successful transmsson probablty of the typcal transmtter n the D-Phase. Next, the non-ppp Φ s approxmated by the homogeneous PPP ˆΦ under Assumpton. We denote the receved SIR at the typcal recever n ˆΦ as SIR (ˆ) = h d α / ˆΦ h l. Thus, P ( SIR () β ) s approxmated by P ( SIR (ˆ) β ). At last, by adoptng the product ofλ andp ( SIR (ˆ) β ) as an approxmate to the spatal capactyc p, a closed-form approxmate to C p for the case of < γ < β s obtaned as C p λ P ( SIR (ˆ) β ) (22) (a) =λ exp ( πλ d 2 γ 2 α ρ ) [ exp πλ exp ( πλ d 2 γ 2 α ρ ) ] d 2 β αρ 2 (23) where (a) follows by Proposton 2. and (3). Note that snce λ = λ P ( SIR () γ ), we can rewrte (22) as C p λ P ( SIR () γ ) P ( SIR (ˆ) β ) under the conventonal method. However, accordng to the defnton of C p forn =, whch s gven n (9) and (), we havec p = λ P ( SIR () γ, SIR () β ), where the dstrbuton of SIR () s strongly dependent on that of SIR () as Φ Φ. As a result, the conventonal approxmate approach only focuses on the PPP-based approxmate to Φ, but gnores the dependence between Φ and Φ. Therefore, (22) does not hold for representng, or reasonably approxmatng, the spatal capacty of the proposed scheme. In addton, by comparng (2) and (23), t s observed that for the case of < γ < β, gven any λ > and

15 5 d >, the closed-form spatal capacty obtaned based on the proposed approach s always outperformed that by the conventonal approach. C. Numercal Results Numercal results are presented n ths subsecton. Accordng to the method descrbed n [5], we generate a spatal Posson process, n whch the transmtters are placed unformly n a square of [m, 6m] [m, 6m]. To take care of the border effects, we focus on samplng the transmtters that locate n the nterm square of [2m, 4m] [2m, 4m]. We calculate the spatal capacty as the average of the network capacty over 2 ndependent network realzatons, where for each network realzaton, the network capacty s evaluated as the rato of the number of successful transmtters n the samplng square to the square area of 4 2 m 2. Unless otherwse specfed, n ths subsecton, we set α=4, β = 2.5, and d=m. We also observe by smulaton that smlar performance can be obtaned by usng other parameters. In the followng, we frst valdate our analytcal results on the spatal capacty of the proposed scheme and the reference scheme wthout schedulng. To hghlght the spatal capacty mprovement performance of the proposed scheme, we then compare the spatal capacty acheved by the propose scheme wth that by two exstng dstrbuted schedulng schemes: one s the probablty-based schedulng n [8], and the other s the channel-threshold based schedulng n [3] and [4]. At last, we consder a more practcal scenaro wth SIR estmaton and feedback errors, and show the effects of the SIR errors on the spatal capacty of the proposed scheme. ) Valdaton of the Spatal Capacty Analyss: We valdate our spatal capacty analyss n Secton III-A and Secton III-B for both proposed and reference schemes. Fg. 2 shows the spatal capacty versus the SIR threshold γ, for both the reference scheme wthout transmsson schedulng and the proposed scheme wth SIR-based schedulng. We set the ntal transmtter densty as λ =.25/m 2 n both schemes. The analytcal spatal capacty of the reference scheme s gven n (8). By comparng the smulaton results for the proposed scheme wth the analytcal results for the reference scheme, we observe that C r s constant overγ as expected. We also observe that ) when γ <β, C p >C r ; 2) when γ = or γ =β, C p =C r ; and 3) when γ >β, C p <C r. Ths s n accordance wth our analytcal results n Proposton 3.. Moreover, for the proposed scheme, we adopt (8) and (5) as the analytcal spatal capacty for the cases of γ = and γ β, respectvely, and observe that the analytcal results of the spatal capacty ft well to the smulaton counterparts. Furthermore, for the case of < γ < β of the proposed scheme, where only approxmate expressons for the spatal capacty are avalable, we compare the approxmate performance of the three approxmate approaches gven n Secton III-B. It s observed that the ntegral-based expresson by the proposed approxmate approach, gven n

16 6 8 x 4 7 (9) (2) 6 Spatal capacty (23) Reference scheme n (8) 2 <γ β γ β Proposed scheme: smulaton Proposed scheme: analytcal SIR threshold γ Fg. 2. Spatal capacty aganst γ wth λ =.25 and β = 2.5. (9), provdes a tght approxmate to C p for the case of < γ < β. In addton, as a cost of expressng n closedform, (2) s not as tght as (9), but (2) stll provdes a close approxmate to C p for the case of < γ < β. At last, t s observed that the closed-form expresson gven n (23) by the conventonal approxmate approach cannot properly approxmate C p for the case of < γ < β as expected. Fg. 3 shows the spatal capacty versus the ntal transmtter densty λ when γ < β. We set γ =.6. For the proposed scheme, smlar to the case n Fg. 2, we observe tght and close approxmates are provded by (9) and (2), respectvely, based on the proposed approxmate approach, whle mproper approxmate s provded by (23) based on the conventonal approxmate approach. Moreover, t s observed that the spatal capacty of the proposed scheme s always larger than that of the reference scheme, gven n (8), for all values of λ, whch s as expected from Proposton 3. snce γ <β n ths example. Furthermore, for both the proposed and reference schemes, we observe an nterestng densty-capacty tradeoff : by ncreasng λ, the spatal capacty frst ncreases due to more avalable transmtters, but as λ exceeds a certan threshold, t starts to decrease, due to the more domnant nterference effect. Thus, to maxmze the spatal capacty, under the system scenaro set n Fg. 3, the optmal λ should be set as.3/m 2. 2) Performance Comparson wth Exstng Dstrbuted Schemes: We consder two exstng dstrbuted schedulng schemes for performance comparson. The frst scheme s the teratve probablty-based schedulng as n [8].

17 7 8 x 4 7 (9) Spatal capacty (2) Reference scheme n (8) (23) Proposed scheme: smulaton Proposed scheme: analytcal Intal transmtter densty λ Fg. 3. Spatal capacty aganst λ. γ =.6. β = 2.5. Denote the transmsson probablty for transmtter n P-Phase k, k N, and the D-Phase as φ (k) or ( ) φ (N), respectvely. For any k {,..,N}, [8] sets φ (k) (k ) = mn SIR β,. Intutvely, [8] provdes a smple and proper way to teratvely adjust the transmsson probablty φ (k). The second scheme s the channel-threshold based schedulng wth sngle-stage probng as n [3] and [4], where the receved nterference power s not nvolved n the transmsson decson and each transmtter decdes to transmt n the D-Phase f ts drect channel strength n P-Phases no smaller than a predefned threshold γ,.e., h γ. For a far comparson, we consder sngle-stage probng wth N = for all the proposed SIR-threshold based scheme, the probablty-based schedulng n [8], and the channel-threshold based schedulng n [3] and [4]. Fg. 4 shows the spatal capactes acheved by the proposed scheme, the probablty-based schedulng, the channel-threshold based schedulng, and the reference scheme wthout schedulng. To clearly show the effects of nvolvng nterference n the transmsson decson for the proposed scheme, we set γ =γ =.4 for the channelthreshold based schedulng. We obtan the spatal capacty of the channel-threshold based schedulng by applyng ts exact expresson gven n [4]. Due to the lack of an exact spatal capacty expresson for the probablty-based schedulng, we obtan ts spatal capacty by smulaton. We lst our observatons from Fg. 4 as follows: SIR based schemes v.s. channel-threshold based scheme: It s observed that by adaptng the transmsson decson to the SIR, the acheved spatal capactes by both the proposed scheme and the probablty-based

18 Prob. based schedulng Proposed scheme Spatal capacty Reference scheme Channel threshold based schedulng smulaton analytcal Intal transmtter densty λ Fg. 4. Spatal capacty comparson wth exstng dstrbuted schedulng schemes. γ =.4. β = 2.5. schedulng are always hgher than that by the channel-threshold based schedulng, where the nterference nformaton s not exploted. Moreover, the spatal capacty of the channel-threshold based schedulng s smaller than that of the reference scheme when λ s small, and becomes larger when λ s suffcently large. Ths s n sharp contrast to the cases of the proposed scheme and the probablty-based schedulng, whch always guarantee capacty mprovement over the reference scheme wthout schedulng. SIR-threshold based schedulng v.s. probablty-based schedulng: It s nterestng to observe that although both the proposed scheme and the probablty-based schedulng adapt the transmsson decson to the SIR, the acheved spatal capacty by the former scheme s always hgher than that by the latter one n ths smulaton. Ths s because that the proposed scheme assures the mprovement of the successful transmsson probablty of each retaned transmtter n the D-Phase, whle the probablty-based schedulng only assures such mprovement wth some probablty. Moreover, t s observed that the optmal ntal transmtter densty that maxmzes the spatal capacty of the proposed scheme s λ =.36, whch s larger than that for the probablty-based schedulng locatng at λ =.26. Note that for the proposed scheme, a lower SIR threshold γ allows more transmtters to retan n the D-Phase, so as to have a second chance to transmt. Thus, by comparng the smulaton results of the proposed scheme n Fg. 4 wth that n Fg. 3, t s observed that the acheved optmal spatal capacty over λ wth γ =.4 n Fg. 4

19 Spatal capacty proposed scheme: w/o SIR error proposed scheme: w/ SIR error Intal transmtter densty λ Fg. 5. Effects of the SIR errors on the spatal capacty of the proposed scheme. γ =.4. β = 2.5. s larger than that wth γ =.6 n Fg. 3. In addton, for all the consdered schemes n Fg. 4, we observe a densty-capacty tradeoff, whch s smlar to that n Fg. 3. 3) Effects of the SIR Estmaton and Feedback Errors: We consder a more practcal scenaro, where SIR estmaton and feedback errors exst n the mplementaton of the proposed scheme, and show the effects of the SIR errors on the spatal capacty. Smlarly to [29], where the channel estmaton and feedback errors are assumed to be zero-mean Gaussan varables, respectvely, we assume the SIR estmaton and feedback errors follow zeromean Gaussan dstrbutons wth varance σest 2 and σ2 fed, respectvely. By further assumng that the two types of SIR errors are mutually ndependent, the sum of both SIR errors at transmtter, denoted by n, follows zero-mean Gaussan dstrbuton wth varance σ 2 = σest+σ 2 fed 2. Thus, n the presence of SIR errors, the feedback SIR level at transmtter n P-Phase s SIR () + n. Moreover, f the feedback SIR level SIR () + n γ for a gven SIR threshold γ n P-Phase, transmtter decdes to transmt n P-Phase ; otherwse, t decdes to be dle n the remanng tme of ths tme slot. Smlar to ts counterpart wthout SIR errors n Fg. 3 and Fg. 4, the spatal capacty wth SIR errors s calculated as an average value over all the transmtters random locatons, the random fadng channels, as well as the random SIR errors. Fg. 5 numercally shows the spatal capactes of the proposed scheme n both cases wth and wthout SIR errors. We set σ 2 = 2 and γ =.4 n ths example. It s observed from Fg. 5 that when λ s small, due to

20 2 the resultant small nterference n the network, each recever feeds back a suffcently hgh SIR level SIR () to ts assocated transmtter, such that n has a small probablty to affect the transmtter s decson. Thus, we observe that when λ s small, the spatal capacty wth SIR errors s tght to that wthout SIR errors. However, as λ ncreases, due to the decreased SIR () at each transmtter, the transmtters become more easly affected by the SIR errors n when decdng whether to transmt based on SIR () +n γ. It s noted that when λ s suffcently large, the average SIR level at each transmtter becomes very small; and even f SIR () γ for transmtter, SIR () s close to γ wth a large probablty. Thus, under the case wth zero-mean Gaussan dstrbuted error n, for the transmtters wth SIR () γ n the SIR error-free case, t s more lkely that these transmtters become SIR () +n < γ than SIR () +n γ n the SIR error-nvolved case. Smlarly, we can easly fnd that for the transmtters wth SIR () < γ n the SIR error-free case, t s also more lkely that these transmtters mantan SIR () + n < γ than SIR () + n γ n the SIR error-nvolved case. Thus, the number of transmtters wth SIR () +n < γ n the SIR error-nvolved case s larger than that wth SIR () < γ n the SIR error-free case n general. Hence, as compared to the case wthout SIR errors, more transmtters wll be refraned from transmttng n the D-Phase n the case wth SIR errors, whch mproves the successful transmsson probablty n the D-Phase due to the reduced nterference. As a result, t s nterestng to observe from Fg. 5 that when the ntal transmtter densty λ ncreases to some sgnfcant pont, the spatal capacty wth SIR errors becomes slghtly hgher than that wthout SIR errors; and ther gap slowly ncreases over λ after ths pont. Therefore, naccurate SIR may even help mprove the SIR-based schedulng performance n more nterference-lmted regme, whch makes the proposed desgn robust to SIR errors. IV. SIR-THRESHOLD BASED SCHEME WITH MULTI-STAGE PROBING In ths secton, we consder the proposed scheme wth mult-stage probng,.e., N >. In ths case, N probng phases are sequentally mplemented to gradually decde the transmtters that are allowed to transmt n the data transmsson phase. Accordng to (5), to fnd the spatal capacty C p,n wth N probng phases, we need to frst fnd the successful transmsson probablty P p,n gven n (3). However, due to the mutually coupled user transmssons over dfferent probng phases, the successful transmsson probablty n P-Phase k, < k N, s related to the SIR dstrbutons n all the proceedng probng phases (from P-Phase to P-Phase k ). Moreover, due to the dfferent pont process formed by the retaned transmtters n each probng phase, the SIR correlatons of any two probng phases are dfferent. Thus, t s challengng to express the successful transmsson probablty and thus the spatal capacty for the case wth N > n general. As a result, nstead of focusng on expressng the spatal capacty C p,n, we focus on studyng how the key system desgn parameters, such as the SIR thresholds and the

21 2 number of probng phases N, affect the spatal capacty of the proposed scheme wth N >. In partcular, unlke the case wth N =, where the sngle-stage overhead Nτ = τ T s neglgble, the mult-stage overhead Nτ wth N > may not be neglgble. In the followng, we frst study the mpact of multple SIR thresholds on the spatal capacty by extendng Proposton 3. for the case of N = to the case of N >. We then nvestgate the effects of the mult-stage probng overhead on the spatal capacty. A. Impact of SIR Thresholds From (3) and (5), the spatal capacty of the proposed scheme s determned by the values of SIR thresholds as well as the tme overhead Nτ for probng. To focus on the mpact of the SIR thresholds, n ths subsecton, we assume Nτ s neglgble and thus have C p,n = λ P ( SIR () γ,..., SIR (N ) γ N, SIR (N) β ) (24) where the dstrbutons of SIR (k) s, k N, are mutually dependent and all the Φ k s, k N, are non-ppps n general. It s also noted that for any k N, we have Φ k Φ k for γ k. Thus, the network nterference level n Φ k s reduced, as compared to that n Φ k. As a result, by extendng Proposton 3. for the case of N =, we obtan the followng proposton for the case of N >. Proposton 4.: Consder two proposed schemes wth arbtrary N and N probng phases, respectvely, N >. Suppose the two schemes adopt the same SIR threshold γ k n each Φ k, k {,...,N }. Then gven β >, by varyng the SIR threshold γ N [, ) n the data transmsson phase for the proposed scheme wth N probng phases, we have the followng relatonshp between C p,n and C p,n based on (24): C p,n > C p,n, f γ N < γ N < β ( conservatve transmsson regme ) C p,n = C p,n, f γ N γ N or γ N = β ( neutral transmsson regme ) C p,n < C p,n, f γ N > β ( aggressve transmsson regme ). Proof: Please refer to Appendx D. (25) Remark 4.: Smlar to the case of Proposton 3., n Proposton 4., n the conservatve transmsson regme wth γ N < γ N < β, we obtan mproved spatal capacty; n the aggressve transmsson regme wth γ N > β, we obtan reduced spatal capacty; and n the neutral transmsson regon wth γ N γ N or γ N = β, we obtan unchanged capacty. Moreover, based on the fact that the conservatve transmsson decson s benefcal for mprovng the spatal capacty of the proposed scheme, we obtan the followng corollary, whch gves a proper method to set the values of all the SIR-thresholds, such that the mprovement of spatal capacty over the number of probng phases s assured.

ECE559VV Project Report

ECE559VV Project Report ECE559VV Project Report (Supplementary Notes Loc Xuan Bu I. MAX SUM-RATE SCHEDULING: THE UPLINK CASE We have seen (n the presentaton that, for downlnk (broadcast channels, the strategy maxmzng the sum-rate

More information

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur Module 3 LOSSY IMAGE COMPRESSION SYSTEMS Verson ECE IIT, Kharagpur Lesson 6 Theory of Quantzaton Verson ECE IIT, Kharagpur Instructonal Objectves At the end of ths lesson, the students should be able to:

More information

Lecture Notes on Linear Regression

Lecture Notes on Linear Regression Lecture Notes on Lnear Regresson Feng L fl@sdueducn Shandong Unversty, Chna Lnear Regresson Problem In regresson problem, we am at predct a contnuous target value gven an nput feature vector We assume

More information

Error Probability for M Signals

Error Probability for M Signals Chapter 3 rror Probablty for M Sgnals In ths chapter we dscuss the error probablty n decdng whch of M sgnals was transmtted over an arbtrary channel. We assume the sgnals are represented by a set of orthonormal

More information

Composite Hypotheses testing

Composite Hypotheses testing Composte ypotheses testng In many hypothess testng problems there are many possble dstrbutons that can occur under each of the hypotheses. The output of the source s a set of parameters (ponts n a parameter

More information

A Robust Method for Calculating the Correlation Coefficient

A Robust Method for Calculating the Correlation Coefficient A Robust Method for Calculatng the Correlaton Coeffcent E.B. Nven and C. V. Deutsch Relatonshps between prmary and secondary data are frequently quantfed usng the correlaton coeffcent; however, the tradtonal

More information

Problem Set 9 Solutions

Problem Set 9 Solutions Desgn and Analyss of Algorthms May 4, 2015 Massachusetts Insttute of Technology 6.046J/18.410J Profs. Erk Demane, Srn Devadas, and Nancy Lynch Problem Set 9 Solutons Problem Set 9 Solutons Ths problem

More information

An Upper Bound on SINR Threshold for Call Admission Control in Multiple-Class CDMA Systems with Imperfect Power-Control

An Upper Bound on SINR Threshold for Call Admission Control in Multiple-Class CDMA Systems with Imperfect Power-Control An Upper Bound on SINR Threshold for Call Admsson Control n Multple-Class CDMA Systems wth Imperfect ower-control Mahmoud El-Sayes MacDonald, Dettwler and Assocates td. (MDA) Toronto, Canada melsayes@hotmal.com

More information

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification

2E Pattern Recognition Solutions to Introduction to Pattern Recognition, Chapter 2: Bayesian pattern classification E395 - Pattern Recognton Solutons to Introducton to Pattern Recognton, Chapter : Bayesan pattern classfcaton Preface Ths document s a soluton manual for selected exercses from Introducton to Pattern Recognton

More information

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X

3.1 Expectation of Functions of Several Random Variables. )' be a k-dimensional discrete or continuous random vector, with joint PMF p (, E X E X1 E X Statstcs 1: Probablty Theory II 37 3 EPECTATION OF SEVERAL RANDOM VARIABLES As n Probablty Theory I, the nterest n most stuatons les not on the actual dstrbuton of a random vector, but rather on a number

More information

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud

Resource Allocation with a Budget Constraint for Computing Independent Tasks in the Cloud Resource Allocaton wth a Budget Constrant for Computng Independent Tasks n the Cloud Wemng Sh and Bo Hong School of Electrcal and Computer Engneerng Georga Insttute of Technology, USA 2nd IEEE Internatonal

More information

Notes on Frequency Estimation in Data Streams

Notes on Frequency Estimation in Data Streams Notes on Frequency Estmaton n Data Streams In (one of) the data streamng model(s), the data s a sequence of arrvals a 1, a 2,..., a m of the form a j = (, v) where s the dentty of the tem and belongs to

More information

VQ widely used in coding speech, image, and video

VQ widely used in coding speech, image, and video at Scalar quantzers are specal cases of vector quantzers (VQ): they are constraned to look at one sample at a tme (memoryless) VQ does not have such constrant better RD perfomance expected Source codng

More information

Joint Statistical Meetings - Biopharmaceutical Section

Joint Statistical Meetings - Biopharmaceutical Section Iteratve Ch-Square Test for Equvalence of Multple Treatment Groups Te-Hua Ng*, U.S. Food and Drug Admnstraton 1401 Rockvlle Pke, #200S, HFM-217, Rockvlle, MD 20852-1448 Key Words: Equvalence Testng; Actve

More information

Lecture 3: Shannon s Theorem

Lecture 3: Shannon s Theorem CSE 533: Error-Correctng Codes (Autumn 006 Lecture 3: Shannon s Theorem October 9, 006 Lecturer: Venkatesan Guruswam Scrbe: Wdad Machmouch 1 Communcaton Model The communcaton model we are usng conssts

More information

A Lower Bound on SINR Threshold for Call Admission Control in Multiple-Class CDMA Systems with Imperfect Power-Control

A Lower Bound on SINR Threshold for Call Admission Control in Multiple-Class CDMA Systems with Imperfect Power-Control A ower Bound on SIR Threshold for Call Admsson Control n Multple-Class CDMA Systems w Imperfect ower-control Mohamed H. Ahmed Faculty of Engneerng and Appled Scence Memoral Unversty of ewfoundland St.

More information

Joint Scheduling and Power-Allocation for Interference Management in Wireless Networks

Joint Scheduling and Power-Allocation for Interference Management in Wireless Networks Jont Schedulng and Power-Allocaton for Interference Management n Wreless Networks Xn Lu *, Edwn K. P. Chong, and Ness B. Shroff * * School of Electrcal and Computer Engneerng Purdue Unversty West Lafayette,

More information

Lecture 12: Classification

Lecture 12: Classification Lecture : Classfcaton g Dscrmnant functons g The optmal Bayes classfer g Quadratc classfers g Eucldean and Mahalanobs metrcs g K Nearest Neghbor Classfers Intellgent Sensor Systems Rcardo Guterrez-Osuna

More information

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS Avalable onlne at http://sck.org J. Math. Comput. Sc. 3 (3), No., 6-3 ISSN: 97-537 COMPARISON OF SOME RELIABILITY CHARACTERISTICS BETWEEN REDUNDANT SYSTEMS REQUIRING SUPPORTING UNITS FOR THEIR OPERATIONS

More information

NUMERICAL DIFFERENTIATION

NUMERICAL DIFFERENTIATION NUMERICAL DIFFERENTIATION 1 Introducton Dfferentaton s a method to compute the rate at whch a dependent output y changes wth respect to the change n the ndependent nput x. Ths rate of change s called the

More information

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4)

Econ107 Applied Econometrics Topic 3: Classical Model (Studenmund, Chapter 4) I. Classcal Assumptons Econ7 Appled Econometrcs Topc 3: Classcal Model (Studenmund, Chapter 4) We have defned OLS and studed some algebrac propertes of OLS. In ths topc we wll study statstcal propertes

More information

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009

College of Computer & Information Science Fall 2009 Northeastern University 20 October 2009 College of Computer & Informaton Scence Fall 2009 Northeastern Unversty 20 October 2009 CS7880: Algorthmc Power Tools Scrbe: Jan Wen and Laura Poplawsk Lecture Outlne: Prmal-dual schema Network Desgn:

More information

Limited Dependent Variables

Limited Dependent Variables Lmted Dependent Varables. What f the left-hand sde varable s not a contnuous thng spread from mnus nfnty to plus nfnty? That s, gven a model = f (, β, ε, where a. s bounded below at zero, such as wages

More information

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011

Stanford University CS359G: Graph Partitioning and Expanders Handout 4 Luca Trevisan January 13, 2011 Stanford Unversty CS359G: Graph Parttonng and Expanders Handout 4 Luca Trevsan January 3, 0 Lecture 4 In whch we prove the dffcult drecton of Cheeger s nequalty. As n the past lectures, consder an undrected

More information

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests

Simulated Power of the Discrete Cramér-von Mises Goodness-of-Fit Tests Smulated of the Cramér-von Mses Goodness-of-Ft Tests Steele, M., Chaselng, J. and 3 Hurst, C. School of Mathematcal and Physcal Scences, James Cook Unversty, Australan School of Envronmental Studes, Grffth

More information

Inductance Calculation for Conductors of Arbitrary Shape

Inductance Calculation for Conductors of Arbitrary Shape CRYO/02/028 Aprl 5, 2002 Inductance Calculaton for Conductors of Arbtrary Shape L. Bottura Dstrbuton: Internal Summary In ths note we descrbe a method for the numercal calculaton of nductances among conductors

More information

Finite Mixture Models and Expectation Maximization. Most slides are from: Dr. Mario Figueiredo, Dr. Anil Jain and Dr. Rong Jin

Finite Mixture Models and Expectation Maximization. Most slides are from: Dr. Mario Figueiredo, Dr. Anil Jain and Dr. Rong Jin Fnte Mxture Models and Expectaton Maxmzaton Most sldes are from: Dr. Maro Fgueredo, Dr. Anl Jan and Dr. Rong Jn Recall: The Supervsed Learnng Problem Gven a set of n samples X {(x, y )},,,n Chapter 3 of

More information

Section 8.3 Polar Form of Complex Numbers

Section 8.3 Polar Form of Complex Numbers 80 Chapter 8 Secton 8 Polar Form of Complex Numbers From prevous classes, you may have encountered magnary numbers the square roots of negatve numbers and, more generally, complex numbers whch are the

More information

Kernel Methods and SVMs Extension

Kernel Methods and SVMs Extension Kernel Methods and SVMs Extenson The purpose of ths document s to revew materal covered n Machne Learnng 1 Supervsed Learnng regardng support vector machnes (SVMs). Ths document also provdes a general

More information

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems

Chapter 5. Solution of System of Linear Equations. Module No. 6. Solution of Inconsistent and Ill Conditioned Systems Numercal Analyss by Dr. Anta Pal Assstant Professor Department of Mathematcs Natonal Insttute of Technology Durgapur Durgapur-713209 emal: anta.bue@gmal.com 1 . Chapter 5 Soluton of System of Lnear Equatons

More information

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U)

ANSWERS. Problem 1. and the moment generating function (mgf) by. defined for any real t. Use this to show that E( U) var( U) Econ 413 Exam 13 H ANSWERS Settet er nndelt 9 deloppgaver, A,B,C, som alle anbefales å telle lkt for å gøre det ltt lettere å stå. Svar er gtt . Unfortunately, there s a prntng error n the hnt of

More information

Convergence of random processes

Convergence of random processes DS-GA 12 Lecture notes 6 Fall 216 Convergence of random processes 1 Introducton In these notes we study convergence of dscrete random processes. Ths allows to characterze phenomena such as the law of large

More information

Module 9. Lecture 6. Duality in Assignment Problems

Module 9. Lecture 6. Duality in Assignment Problems Module 9 1 Lecture 6 Dualty n Assgnment Problems In ths lecture we attempt to answer few other mportant questons posed n earler lecture for (AP) and see how some of them can be explaned through the concept

More information

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering /

P R. Lecture 4. Theory and Applications of Pattern Recognition. Dept. of Electrical and Computer Engineering / Theory and Applcatons of Pattern Recognton 003, Rob Polkar, Rowan Unversty, Glassboro, NJ Lecture 4 Bayes Classfcaton Rule Dept. of Electrcal and Computer Engneerng 0909.40.0 / 0909.504.04 Theory & Applcatons

More information

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1

On an Extension of Stochastic Approximation EM Algorithm for Incomplete Data Problems. Vahid Tadayon 1 On an Extenson of Stochastc Approxmaton EM Algorthm for Incomplete Data Problems Vahd Tadayon Abstract: The Stochastc Approxmaton EM (SAEM algorthm, a varant stochastc approxmaton of EM, s a versatle tool

More information

Lossy Compression. Compromise accuracy of reconstruction for increased compression.

Lossy Compression. Compromise accuracy of reconstruction for increased compression. Lossy Compresson Compromse accuracy of reconstructon for ncreased compresson. The reconstructon s usually vsbly ndstngushable from the orgnal mage. Typcally, one can get up to 0:1 compresson wth almost

More information

Difference Equations

Difference Equations Dfference Equatons c Jan Vrbk 1 Bascs Suppose a sequence of numbers, say a 0,a 1,a,a 3,... s defned by a certan general relatonshp between, say, three consecutve values of the sequence, e.g. a + +3a +1

More information

Negative Binomial Regression

Negative Binomial Regression STATGRAPHICS Rev. 9/16/2013 Negatve Bnomal Regresson Summary... 1 Data Input... 3 Statstcal Model... 3 Analyss Summary... 4 Analyss Optons... 7 Plot of Ftted Model... 8 Observed Versus Predcted... 10 Predctons...

More information

Optimal Pricing and Load Sharing for Energy Saving with Cooperative Communications

Optimal Pricing and Load Sharing for Energy Saving with Cooperative Communications Optmal Prcng and Load Sharng for Energy Savng wth Cooperatve Communcatons Ynghao Guo, Lnge Duan, and Ru Zhang arxv:1409.8402v4 [cs.it] 4 Aug 2015 Abstract Cooperatve communcatons has long been proposed

More information

More metrics on cartesian products

More metrics on cartesian products More metrcs on cartesan products If (X, d ) are metrc spaces for 1 n, then n Secton II4 of the lecture notes we defned three metrcs on X whose underlyng topologes are the product topology The purpose of

More information

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers

Psychology 282 Lecture #24 Outline Regression Diagnostics: Outliers Psychology 282 Lecture #24 Outlne Regresson Dagnostcs: Outlers In an earler lecture we studed the statstcal assumptons underlyng the regresson model, ncludng the followng ponts: Formal statement of assumptons.

More information

x = , so that calculated

x = , so that calculated Stat 4, secton Sngle Factor ANOVA notes by Tm Plachowsk n chapter 8 we conducted hypothess tests n whch we compared a sngle sample s mean or proporton to some hypotheszed value Chapter 9 expanded ths to

More information

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1]

Outline. Communication. Bellman Ford Algorithm. Bellman Ford Example. Bellman Ford Shortest Path [1] DYNAMIC SHORTEST PATH SEARCH AND SYNCHRONIZED TASK SWITCHING Jay Wagenpfel, Adran Trachte 2 Outlne Shortest Communcaton Path Searchng Bellmann Ford algorthm Algorthm for dynamc case Modfcatons to our algorthm

More information

DUE: WEDS FEB 21ST 2018

DUE: WEDS FEB 21ST 2018 HOMEWORK # 1: FINITE DIFFERENCES IN ONE DIMENSION DUE: WEDS FEB 21ST 2018 1. Theory Beam bendng s a classcal engneerng analyss. The tradtonal soluton technque makes smplfyng assumptons such as a constant

More information

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010

Parametric fractional imputation for missing data analysis. Jae Kwang Kim Survey Working Group Seminar March 29, 2010 Parametrc fractonal mputaton for mssng data analyss Jae Kwang Km Survey Workng Group Semnar March 29, 2010 1 Outlne Introducton Proposed method Fractonal mputaton Approxmaton Varance estmaton Multple mputaton

More information

Linear Approximation with Regularization and Moving Least Squares

Linear Approximation with Regularization and Moving Least Squares Lnear Approxmaton wth Regularzaton and Movng Least Squares Igor Grešovn May 007 Revson 4.6 (Revson : March 004). 5 4 3 0.5 3 3.5 4 Contents: Lnear Fttng...4. Weghted Least Squares n Functon Approxmaton...

More information

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE

CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE CHAPTER 5 NUMERICAL EVALUATION OF DYNAMIC RESPONSE Analytcal soluton s usually not possble when exctaton vares arbtrarly wth tme or f the system s nonlnear. Such problems can be solved by numercal tmesteppng

More information

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA

4 Analysis of Variance (ANOVA) 5 ANOVA. 5.1 Introduction. 5.2 Fixed Effects ANOVA 4 Analyss of Varance (ANOVA) 5 ANOVA 51 Introducton ANOVA ANOVA s a way to estmate and test the means of multple populatons We wll start wth one-way ANOVA If the populatons ncluded n the study are selected

More information

Cognitive Access Algorithms For Multiple Access Channels

Cognitive Access Algorithms For Multiple Access Channels 203 IEEE 4th Workshop on Sgnal Processng Advances n Wreless Communcatons SPAWC) Cogntve Access Algorthms For Multple Access Channels Ychuan Hu and Alejandro Rbero, Department of Electrcal and Systems Engneerng,

More information

Chapter 13: Multiple Regression

Chapter 13: Multiple Regression Chapter 13: Multple Regresson 13.1 Developng the multple-regresson Model The general model can be descrbed as: It smplfes for two ndependent varables: The sample ft parameter b 0, b 1, and b are used to

More information

Pricing and Resource Allocation Game Theoretic Models

Pricing and Resource Allocation Game Theoretic Models Prcng and Resource Allocaton Game Theoretc Models Zhy Huang Changbn Lu Q Zhang Computer and Informaton Scence December 8, 2009 Z. Huang, C. Lu, and Q. Zhang (CIS) Game Theoretc Models December 8, 2009

More information

Finding Dense Subgraphs in G(n, 1/2)

Finding Dense Subgraphs in G(n, 1/2) Fndng Dense Subgraphs n Gn, 1/ Atsh Das Sarma 1, Amt Deshpande, and Rav Kannan 1 Georga Insttute of Technology,atsh@cc.gatech.edu Mcrosoft Research-Bangalore,amtdesh,annan@mcrosoft.com Abstract. Fndng

More information

x = x 1 + :::+ x K and the nput covarance matrces are of the form ± = E[x x y ]. 3.2 Dualty Next, we ntroduce the concept of dualty wth the followng t

x = x 1 + :::+ x K and the nput covarance matrces are of the form ± = E[x x y ]. 3.2 Dualty Next, we ntroduce the concept of dualty wth the followng t Sum Power Iteratve Water-fllng for Mult-Antenna Gaussan Broadcast Channels N. Jndal, S. Jafar, S. Vshwanath and A. Goldsmth Dept. of Electrcal Engg. Stanford Unversty, CA, 94305 emal: njndal,syed,srram,andrea@wsl.stanford.edu

More information

Stat260: Bayesian Modeling and Inference Lecture Date: February 22, Reference Priors

Stat260: Bayesian Modeling and Inference Lecture Date: February 22, Reference Priors Stat60: Bayesan Modelng and Inference Lecture Date: February, 00 Reference Prors Lecturer: Mchael I. Jordan Scrbe: Steven Troxler and Wayne Lee In ths lecture, we assume that θ R; n hgher-dmensons, reference

More information

Chapter - 2. Distribution System Power Flow Analysis

Chapter - 2. Distribution System Power Flow Analysis Chapter - 2 Dstrbuton System Power Flow Analyss CHAPTER - 2 Radal Dstrbuton System Load Flow 2.1 Introducton Load flow s an mportant tool [66] for analyzng electrcal power system network performance. Load

More information

Global Sensitivity. Tuesday 20 th February, 2018

Global Sensitivity. Tuesday 20 th February, 2018 Global Senstvty Tuesday 2 th February, 28 ) Local Senstvty Most senstvty analyses [] are based on local estmates of senstvty, typcally by expandng the response n a Taylor seres about some specfc values

More information

Lecture 4. Instructor: Haipeng Luo

Lecture 4. Instructor: Haipeng Luo Lecture 4 Instructor: Hapeng Luo In the followng lectures, we focus on the expert problem and study more adaptve algorthms. Although Hedge s proven to be worst-case optmal, one may wonder how well t would

More information

Average Decision Threshold of CA CFAR and excision CFAR Detectors in the Presence of Strong Pulse Jamming 1

Average Decision Threshold of CA CFAR and excision CFAR Detectors in the Presence of Strong Pulse Jamming 1 Average Decson hreshold of CA CFAR and excson CFAR Detectors n the Presence of Strong Pulse Jammng Ivan G. Garvanov and Chrsto A. Kabachev Insttute of Informaton echnologes Bulgaran Academy of Scences

More information

STAT 511 FINAL EXAM NAME Spring 2001

STAT 511 FINAL EXAM NAME Spring 2001 STAT 5 FINAL EXAM NAME Sprng Instructons: Ths s a closed book exam. No notes or books are allowed. ou may use a calculator but you are not allowed to store notes or formulas n the calculator. Please wrte

More information

Real-Time Systems. Multiprocessor scheduling. Multiprocessor scheduling. Multiprocessor scheduling

Real-Time Systems. Multiprocessor scheduling. Multiprocessor scheduling. Multiprocessor scheduling Real-Tme Systems Multprocessor schedulng Specfcaton Implementaton Verfcaton Multprocessor schedulng -- -- Global schedulng How are tasks assgned to processors? Statc assgnment The processor(s) used for

More information

Comparison of Regression Lines

Comparison of Regression Lines STATGRAPHICS Rev. 9/13/2013 Comparson of Regresson Lnes Summary... 1 Data Input... 3 Analyss Summary... 4 Plot of Ftted Model... 6 Condtonal Sums of Squares... 6 Analyss Optons... 7 Forecasts... 8 Confdence

More information

The Geometry of Logit and Probit

The Geometry of Logit and Probit The Geometry of Logt and Probt Ths short note s meant as a supplement to Chapters and 3 of Spatal Models of Parlamentary Votng and the notaton and reference to fgures n the text below s to those two chapters.

More information

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan

Winter 2008 CS567 Stochastic Linear/Integer Programming Guest Lecturer: Xu, Huan Wnter 2008 CS567 Stochastc Lnear/Integer Programmng Guest Lecturer: Xu, Huan Class 2: More Modelng Examples 1 Capacty Expanson Capacty expanson models optmal choces of the tmng and levels of nvestments

More information

Secret Communication using Artificial Noise

Secret Communication using Artificial Noise Secret Communcaton usng Artfcal Nose Roht Neg, Satashu Goel C Department, Carnege Mellon Unversty, PA 151, USA {neg,satashug}@ece.cmu.edu Abstract The problem of secret communcaton between two nodes over

More information

Rethinking MIMO for Wireless Networks: Linear Throughput Increases with Multiple Receive Antennas

Rethinking MIMO for Wireless Networks: Linear Throughput Increases with Multiple Receive Antennas Retnng MIMO for Wreless etwors: Lnear Trougput Increases wt Multple Receve Antennas ar Jndal Unversty of Mnnesota Unverstat Pompeu Fabra Jont wor wt Jeff Andrews & Steven Weber MIMO n Pont-to-Pont Cannels

More information

Lecture 14 (03/27/18). Channels. Decoding. Preview of the Capacity Theorem.

Lecture 14 (03/27/18). Channels. Decoding. Preview of the Capacity Theorem. Lecture 14 (03/27/18). Channels. Decodng. Prevew of the Capacty Theorem. A. Barg The concept of a communcaton channel n nformaton theory s an abstracton for transmttng dgtal (and analog) nformaton from

More information

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity

Week3, Chapter 4. Position and Displacement. Motion in Two Dimensions. Instantaneous Velocity. Average Velocity Week3, Chapter 4 Moton n Two Dmensons Lecture Quz A partcle confned to moton along the x axs moves wth constant acceleraton from x =.0 m to x = 8.0 m durng a 1-s tme nterval. The velocty of the partcle

More information

Distributed Non-Autonomous Power Control through Distributed Convex Optimization

Distributed Non-Autonomous Power Control through Distributed Convex Optimization Dstrbuted Non-Autonomous Power Control through Dstrbuted Convex Optmzaton S. Sundhar Ram and V. V. Veeravall ECE Department and Coordnated Scence Lab Unversty of Illnos at Urbana-Champagn Emal: {ssrnv5,vvv}@llnos.edu

More information

Bayesian predictive Configural Frequency Analysis

Bayesian predictive Configural Frequency Analysis Psychologcal Test and Assessment Modelng, Volume 54, 2012 (3), 285-292 Bayesan predctve Confgural Frequency Analyss Eduardo Gutérrez-Peña 1 Abstract Confgural Frequency Analyss s a method for cell-wse

More information

Supplementary Notes for Chapter 9 Mixture Thermodynamics

Supplementary Notes for Chapter 9 Mixture Thermodynamics Supplementary Notes for Chapter 9 Mxture Thermodynamcs Key ponts Nne major topcs of Chapter 9 are revewed below: 1. Notaton and operatonal equatons for mxtures 2. PVTN EOSs for mxtures 3. General effects

More information

The optimal delay of the second test is therefore approximately 210 hours earlier than =2.

The optimal delay of the second test is therefore approximately 210 hours earlier than =2. THE IEC 61508 FORMULAS 223 The optmal delay of the second test s therefore approxmately 210 hours earler than =2. 8.4 The IEC 61508 Formulas IEC 61508-6 provdes approxmaton formulas for the PF for smple

More information

One-sided finite-difference approximations suitable for use with Richardson extrapolation

One-sided finite-difference approximations suitable for use with Richardson extrapolation Journal of Computatonal Physcs 219 (2006) 13 20 Short note One-sded fnte-dfference approxmatons sutable for use wth Rchardson extrapolaton Kumar Rahul, S.N. Bhattacharyya * Department of Mechancal Engneerng,

More information

A Hybrid Variational Iteration Method for Blasius Equation

A Hybrid Variational Iteration Method for Blasius Equation Avalable at http://pvamu.edu/aam Appl. Appl. Math. ISSN: 1932-9466 Vol. 10, Issue 1 (June 2015), pp. 223-229 Applcatons and Appled Mathematcs: An Internatonal Journal (AAM) A Hybrd Varatonal Iteraton Method

More information

Markov Chain Monte Carlo Lecture 6

Markov Chain Monte Carlo Lecture 6 where (x 1,..., x N ) X N, N s called the populaton sze, f(x) f (x) for at least one {1, 2,..., N}, and those dfferent from f(x) are called the tral dstrbutons n terms of mportance samplng. Dfferent ways

More information

Multigradient for Neural Networks for Equalizers 1

Multigradient for Neural Networks for Equalizers 1 Multgradent for Neural Netorks for Equalzers 1 Chulhee ee, Jnook Go and Heeyoung Km Department of Electrcal and Electronc Engneerng Yonse Unversty 134 Shnchon-Dong, Seodaemun-Ku, Seoul 1-749, Korea ABSTRACT

More information

Hidden Markov Models & The Multivariate Gaussian (10/26/04)

Hidden Markov Models & The Multivariate Gaussian (10/26/04) CS281A/Stat241A: Statstcal Learnng Theory Hdden Markov Models & The Multvarate Gaussan (10/26/04) Lecturer: Mchael I. Jordan Scrbes: Jonathan W. Hu 1 Hdden Markov Models As a bref revew, hdden Markov models

More information

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 30 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 2 Remedes for multcollnearty Varous technques have

More information

Tracking with Kalman Filter

Tracking with Kalman Filter Trackng wth Kalman Flter Scott T. Acton Vrgna Image and Vdeo Analyss (VIVA), Charles L. Brown Department of Electrcal and Computer Engneerng Department of Bomedcal Engneerng Unversty of Vrgna, Charlottesvlle,

More information

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper

Games of Threats. Elon Kohlberg Abraham Neyman. Working Paper Games of Threats Elon Kohlberg Abraham Neyman Workng Paper 18-023 Games of Threats Elon Kohlberg Harvard Busness School Abraham Neyman The Hebrew Unversty of Jerusalem Workng Paper 18-023 Copyrght 2017

More information

Rate-Memory Trade-off for the Two-User Broadcast Caching Network with Correlated Sources

Rate-Memory Trade-off for the Two-User Broadcast Caching Network with Correlated Sources Rate-Memory Trade-off for the Two-User Broadcast Cachng Network wth Correlated Sources Parsa Hassanzadeh, Antona M. Tulno, Jame Llorca, Elza Erkp arxv:705.0466v [cs.it] May 07 Abstract Ths paper studes

More information

Lecture 12: Discrete Laplacian

Lecture 12: Discrete Laplacian Lecture 12: Dscrete Laplacan Scrbe: Tanye Lu Our goal s to come up wth a dscrete verson of Laplacan operator for trangulated surfaces, so that we can use t n practce to solve related problems We are mostly

More information

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that

Online Appendix. t=1 (p t w)q t. Then the first order condition shows that Artcle forthcomng to ; manuscrpt no (Please, provde the manuscrpt number!) 1 Onlne Appendx Appendx E: Proofs Proof of Proposton 1 Frst we derve the equlbrum when the manufacturer does not vertcally ntegrate

More information

Power Allocation for Distributed BLUE Estimation with Full and Limited Feedback of CSI

Power Allocation for Distributed BLUE Estimation with Full and Limited Feedback of CSI Power Allocaton for Dstrbuted BLUE Estmaton wth Full and Lmted Feedback of CSI Mohammad Fanae, Matthew C. Valent, and Natala A. Schmd Lane Department of Computer Scence and Electrcal Engneerng West Vrgna

More information

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography CSc 6974 and ECSE 6966 Math. Tech. for Vson, Graphcs and Robotcs Lecture 21, Aprl 17, 2006 Estmatng A Plane Homography Overvew We contnue wth a dscusson of the major ssues, usng estmaton of plane projectve

More information

January Examinations 2015

January Examinations 2015 24/5 Canddates Only January Examnatons 25 DO NOT OPEN THE QUESTION PAPER UNTIL INSTRUCTED TO DO SO BY THE CHIEF INVIGILATOR STUDENT CANDIDATE NO.. Department Module Code Module Ttle Exam Duraton (n words)

More information

Assortment Optimization under MNL

Assortment Optimization under MNL Assortment Optmzaton under MNL Haotan Song Aprl 30, 2017 1 Introducton The assortment optmzaton problem ams to fnd the revenue-maxmzng assortment of products to offer when the prces of products are fxed.

More information

Time-Varying Systems and Computations Lecture 6

Time-Varying Systems and Computations Lecture 6 Tme-Varyng Systems and Computatons Lecture 6 Klaus Depold 14. Januar 2014 The Kalman Flter The Kalman estmaton flter attempts to estmate the actual state of an unknown dscrete dynamcal system, gven nosy

More information

A 2D Bounded Linear Program (H,c) 2D Linear Programming

A 2D Bounded Linear Program (H,c) 2D Linear Programming A 2D Bounded Lnear Program (H,c) h 3 v h 8 h 5 c h 4 h h 6 h 7 h 2 2D Lnear Programmng C s a polygonal regon, the ntersecton of n halfplanes. (H, c) s nfeasble, as C s empty. Feasble regon C s unbounded

More information

Appendix B: Resampling Algorithms

Appendix B: Resampling Algorithms 407 Appendx B: Resamplng Algorthms A common problem of all partcle flters s the degeneracy of weghts, whch conssts of the unbounded ncrease of the varance of the mportance weghts ω [ ] of the partcles

More information

Market structure and Innovation

Market structure and Innovation Market structure and Innovaton Ths presentaton s based on the paper Market structure and Innovaton authored by Glenn C. Loury, publshed n The Quarterly Journal of Economcs, Vol. 93, No.3 ( Aug 1979) I.

More information

Supplementary material: Margin based PU Learning. Matrix Concentration Inequalities

Supplementary material: Margin based PU Learning. Matrix Concentration Inequalities Supplementary materal: Margn based PU Learnng We gve the complete proofs of Theorem and n Secton We frst ntroduce the well-known concentraton nequalty, so the covarance estmator can be bounded Then we

More information

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications

MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.265/15.070J Fall 2013 Lecture 12 10/21/2013. Martingale Concentration Inequalities and Applications MASSACHUSETTS INSTITUTE OF TECHNOLOGY 6.65/15.070J Fall 013 Lecture 1 10/1/013 Martngale Concentraton Inequaltes and Applcatons Content. 1. Exponental concentraton for martngales wth bounded ncrements.

More information

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction

The Multiple Classical Linear Regression Model (CLRM): Specification and Assumptions. 1. Introduction ECONOMICS 5* -- NOTE (Summary) ECON 5* -- NOTE The Multple Classcal Lnear Regresson Model (CLRM): Specfcaton and Assumptons. Introducton CLRM stands for the Classcal Lnear Regresson Model. The CLRM s also

More information

Probability Theory. The nth coefficient of the Taylor series of f(k), expanded around k = 0, gives the nth moment of x as ( ik) n n!

Probability Theory. The nth coefficient of the Taylor series of f(k), expanded around k = 0, gives the nth moment of x as ( ik) n n! 8333: Statstcal Mechancs I Problem Set # 3 Solutons Fall 3 Characterstc Functons: Probablty Theory The characterstc functon s defned by fk ep k = ep kpd The nth coeffcent of the Taylor seres of fk epanded

More information

Using the estimated penetrances to determine the range of the underlying genetic model in casecontrol

Using the estimated penetrances to determine the range of the underlying genetic model in casecontrol Georgetown Unversty From the SelectedWorks of Mark J Meyer 8 Usng the estmated penetrances to determne the range of the underlyng genetc model n casecontrol desgn Mark J Meyer Neal Jeffres Gang Zheng Avalable

More information

Lecture 3. Ax x i a i. i i

Lecture 3. Ax x i a i. i i 18.409 The Behavor of Algorthms n Practce 2/14/2 Lecturer: Dan Spelman Lecture 3 Scrbe: Arvnd Sankar 1 Largest sngular value In order to bound the condton number, we need an upper bound on the largest

More information

Signal space Review on vector space Linear independence Metric space and norm Inner product

Signal space Review on vector space Linear independence Metric space and norm Inner product Sgnal space.... Revew on vector space.... Lnear ndependence... 3.3 Metrc space and norm... 4.4 Inner product... 5.5 Orthonormal bass... 7.6 Waveform communcaton system... 9.7 Some examples... 6 Sgnal space

More information

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS

A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS HCMC Unversty of Pedagogy Thong Nguyen Huu et al. A PROBABILITY-DRIVEN SEARCH ALGORITHM FOR SOLVING MULTI-OBJECTIVE OPTIMIZATION PROBLEMS Thong Nguyen Huu and Hao Tran Van Department of mathematcs-nformaton,

More information

Chapter 8 Indicator Variables

Chapter 8 Indicator Variables Chapter 8 Indcator Varables In general, e explanatory varables n any regresson analyss are assumed to be quanttatve n nature. For example, e varables lke temperature, dstance, age etc. are quanttatve n

More information

Consider the following passband digital communication system model. c t. modulator. t r a n s m i t t e r. signal decoder.

Consider the following passband digital communication system model. c t. modulator. t r a n s m i t t e r. signal decoder. PASSBAND DIGITAL MODULATION TECHNIQUES Consder the followng passband dgtal communcaton system model. cos( ω + φ ) c t message source m sgnal encoder s modulator s () t communcaton xt () channel t r a n

More information