A Distributed Spectrum Sharing Algorithm in Cognitive Radio Networks

Size: px
Start display at page:

Download "A Distributed Spectrum Sharing Algorithm in Cognitive Radio Networks"

Transcription

1 A Dstrbuted Spectrum Sharng Algorthm n Cogntve Rado Networks We Sun, Tong Lu, Yanmn Zhu, Bo L Shangha Jao Tong Unversty Hong Kong Unversty of Scence and Technology Emal: {sunwt,lutong ,yzhu}@sjtu.edu.cn, bl@cse.ust.hk Abstract In ths paper we study a socal welfare maxmzaton problem for spectrum sharng n cogntve rado networks. To fully use the spectrum resource, the lcensed spectrum owned by the prmary user PU) can be leased to secondary users SUs) for transmttng data. We frst formulate the socal welfare of a cogntve rado network, consderng the cost for the prmary user sharng spectrum and the utlty ganed for secondary users transmttng data. The socal welfare maxmzaton s a convex optmzaton, whch can be solved by standard methods n a centralzed manner. However, the utlty functon of each secondary user always contans the prvate nformaton, whch leads to the centralzed methods dsabled. To overcome ths challenge, we propose an teratve dstrbuted algorthm based on a prcng-based decomposton framework. It s theoretcally proved that our proposed algorthm converges to the optmal soluton. Numercal smulaton results are presented to show that our proposed algorthm acheves optmal socal welfare and fast convergence speed. Index Terms Decomposton, socal welfare maxmzaton, optmzaton, cogntve rado network. I. INTRODUCTION Wth the ncreasng development of wreless communcatons, spectrum becomes a more and more lmted resource. However, many studes [1] [2] [3] show that spectrum s under utlzed n realty. The technology of cogntve rado CR) [4] provdes a promsng mechansm for flexble usage of spectrum. The utlzaton of spectrum can be mproved by sharng vacant spectrum wth unlcensed users. We consder a cogntve rado network, whch conssts of a prmary user PU), the lcensed owner of spectrum, and several secondary users wthout lcences. Each secondary user conssts of a par of nodes, a sender and a recever, whch have the demand of transmttng data. To fully utlze spectrum, the prmary user can lease ts vacant spectrum to secondary users wth monetary rewards. We consder the problem of spectrum sharng among secondary users n a cogntve rado network. A cost s ncurred and pad by the prmary user for allocatng spectrum to secondary users. Each secondary user obtans an amount of utlty for transmttng data, whch contans ts prvate nformaton. We frst formulate the cost and the utltes, respectvely. Then, our objectve of such a cogntve rado network s to maxmze the socal welfare, whch s defned as the sum of the utlty obtaned by each secondary user mnus the cost pad by the prmary user. Fg. 1. PU Allocated spectrum SU1 SU SU 2 3 The llustraton of spectrum allocaton n a cogntve rado network. A number of exstng works [5] [6] have been done to solve the problem of spectrum sharng n cogntve rado networks. Some of these works [7] [8] [9] take advantage of the models n game theory to characterze cogntve rado networks, whch cannot acheve socal welfare maxmzaton. These works always assume that secondary users are ratonal and selfsh. Some other works [1] [11] [12] can really maxmze the network utlty for spectrum sharng. However, the algorthms desgned n these works are executed n a centralzed manner. Therefore, the prvacy of secondary users s not protected. Dfferent wth these works, we focus on desgnng a dstrbuted algorthm to maxmze the socal welfare. Such an optmal spectrum sharng problem s sgnfcantly challengng because of unque characterstcs of cogntve rado networks. Frst, the utlty functon of each secondary user contans the user-specfc nformaton, whch s ther prvacy. It means that the problem cannot be solved n a centralzed manner by the standard methods such as Newton method. A dstrbuted algorthm s necessary by decouplng the problem nto several subproblems, each of whch can be ndependently solved by a secondary user. Second, the objectve n the problem, maxmzng socal welfare, s coupled wth the cost functon of the prmary user and the utlty functon of each secondary user. In other words, the objectve functon cannot be drectly decoupled nto several ndependent subproblems. In response to the challenges mentoned above, we propose a fully dstrbuted algorthm for spectrum sharng n cogntve

2 2 rado networks. As the shared spectrum s prced by the prmary user, we fnd that the optmal socal welfare s obtaned when optmal soluton of an equvalent prcng problem s acheved. Then, we propose a prcng-based decomposton framework, whch decouples the problem nto several subproblems. The prmary user and each secondary user s responsble for a subproblem, whch only contans ts own prvacy. Based on the decomposton framework, a dstrbuted algorthm s provded, whch s executed teratvely. In each teraton, each secondary user nforms the prmary user ts demanded spectrum, accordng to the gven prce. The prmary user updates the prce based on the collected nformaton. We also propose a prce updatng rule, under whch our teratve algorthm s rgorously proved to converge to the optmal soluton n fnte steps. Extensve smulatons have been conducted to evaluate the socal welfare and convergence speed acheved by our dstrbuted algorthm. We have made three techncal contrbutons n our paper. We model the cost pad by the prmary user for spectrum sharng and the utlty obtaned by secondary users. The concept of socal welfare s formally defned as the objectve of a cogntve rado network to pursue. To maxmze the socal welfare n a dstrbuted manner, we propose a novel prcng-based decomposton framework. A dstrbuted algorthm s desgned based on the decomposton framework, where teratve operatons are executed by the prmary user and the secondary users n turn. Rgorous theoretcal analyss s conducted to demonstrate that our teratve algorthm converges to the optmal soluton n fnte steps. Extensve smulaton results show that our dstrbuted algorthm can acheve the optmal socal welfare, and converge quckly. The remander of ths paper s organzed as follows. We ntroduce the network model and problem formulaton n Secton II. The optmal dstrbuted approach s descrbed n Secton III. In Secton IV, smulaton results are reported to show the performance of our algorthm and compared algorthms. We revew related work n Secton V and conclude the paper n Secton VI. II. NETWORK MODEL AND PROBLEM FORMULATION In ths secton, we frst descrbe the model of cogntve rado networks, and then formulate the problem solved n ths paper. A. Network Model We consder a cogntve rado network whch conssts a prmary user and n secondary users. The set of secondary users s denoted by {SU 1, SU 2,, SU n }. Each secondary user actually conssts of a sender and a recever. The prmary user, the owner of the spectrum, allocates ts unused spectrum to secondary users. Each secondary user s allocated wth a porton of spectrum, namely bands. Each band s assocated wth a bandwdth representng a dvded frequency range. The bandwdth allocated to secondary user SU s noted by p, 1 n. We assume that the prmary user can communcate wth each secondary user drectly. Moreover, the frequency allocated to dfferent secondary user s orthogonal to avod nterference. The network model s llustrated n Fg. 1. B. Problem Formulaton In ths subsecton, we formulate the socal welfare of a cogntve rado network, ncludng two folds: the cost pad by the prmary user for spectrum allocaton and the utlty obtaned by the secondary users for data transmsson. Cost for spectrum allocaton. The spectrum sharng wth secondary users may mpact the data transmsson of the prmary user. Ths leads to a hghly ncreasng cost assocated wth the allocated bandwdth, whch s formulated as an exponental functon. We defne the cost functon as follows, Defnton 1. The cost functon of the prmary user s formulated as p C p) = βe 1 1) where p = p 1, p 2,, p m ) and β s a system parameter. The prmary user expects a payment n return from secondary users for usng spectrum. The prmary user provdes a prce q for the spectrum spent by the SU. Consderng ths payment, we formulate the proft of the prmary user n Defnton 2. Defnton 2. The prmary user obtans a proft for sharng spectrum as Φ p, q) = pq T C p). 2) Utlty for data transmsson. It s ntutve that each secondary user obtans a utlty for transmttng data. Accordng to the bandwdth p allocated to the secondary user, we defne the utlty functon n the followng, Defnton 3. The utlty functon of secondary user SU s defned as u p ) = a log 1 + p ), 3) where a s a user-specfc parameter. Accordng to the law of dmnshng margnal utlty n economcs [13], the log functon mples that the ncreasng rato of utlty gan obtaned by secondary users degrades as the allocated bandwdth ncreases. Note that the parameter a s the prvacy nformaton of secondary user SU. whch s unknown to the prmary user and other secondary users. Thus, the payoff of SU s formulated as follows, Defnton 4. The payoff of secondary user SU s equal to the utlty obtaned mnus the payment pad to the prmary user, P q, p ) = u p ) q p. 4) In ths paper, we consder secondary users are ratonal, tryng to maxmze ther own payoffs gven a prce q by the secondary user. Socal welfare. We ntroduce the concept of socal welfare, whch s equal to the utltes ganed by secondary users

3 3 mnus the cost pad by the prmary user, s formulated n the followng defnton. Defnton 5. The socal welfare of a cogntve rado network s defned as S p) = u p ) C p). 5) In ths paper, we am to maxmze the socal welfare of a cogntve rado network. The problem of socal welfare maxmzaton s defned n the followng. Defnton 6. The problem of maxmzng the socal welfare of a cogntve rado network s defned as max u p ) Cp) 6) s.t. p [, δ ], [1, n], 7) Varable SU n p q β a Cp) Φp, q) u p ) P q, p ) Sp) t q t) p t) ε TABLE I MAJOR ADOPTED NOTATIONS. Descrpton The -th secondary user The number of secondary users The bandwdth of spectrum allocated to the -th secondary user The prce gven to the -th secondary user The system parameter n the cost functon of the prmary user The user-specfc parameter n the utlty of the -th secondary user The cost functon of the prmary user The proft functon of the prmary user The utlty functon of the -th secondary user The payoff functon of the -th secondary user The socal welfare of a cogntve rado network The current number of teratons The prce updated for the -th secondary user n the t-th teraton The value of spectrum computed by the -th secondary user n the t-th teraton The system varable n the prce updatng rule where δ s the upper bound of p, determned by secondary user SU. Note that the problem defned n Defnton 6 s a convex optmzaton problem. Accordng to [14], the problem must have an optmal soluton p = [p 1, p 2,, p m]. The problem of maxmzng the socal welfare n a cogntve rado network should be solved n a dstrbuted manner. Frst, f the problem s centrally solved n one pont, t ncurs a huge computaton load on the sngle pont. Second, the utlty functon of each secondary user s prvacy, whch leads to unknowng the exact formulaton of socal welfare. A. Overvew III. OPTIMAL DISTRIBUTED APPROACH In response to the dffcultes descrbed above, we desgn a dstrbuted algorthm based on a prcng-based decomposton method. The algorthm can acheve the maxmal socal welfare and protect the prvate nformaton for each secondary user at the same tme. We frst convert the orgnal problem n 6) nto an e- quvalent prcng problem. The optmal spectrum allocaton can be obtaned when the optmal prces are acheved. Then, we propose a decomposton framework, whch decouples the optmzaton problem nto several subproblems. In the framework, each secondary user maxmzes ts own payoff based on the prce gven by the prmary user. The prmary user updates the prce after collectng the returned spectrum demand from secondary users. The above operatons are executed teratvely untl the spectrum allocaton converges to the optmal soluton. A prce updatng rule s also provded, gven the nsght of the condton that the optmal prces should satsfy. At last, we theoretcally prove that our teratve algorthm converges to the optmum n fnte steps. The notatons appeared n ths paper are summarzed n Table I. Algorthm 1 Optmal Dstrbuted Algorthm for Spectrum Sharng Input: A cogntve rado network consstng of many secondary users, denoted by {SU 1, SU 2,, SU n }. The cost functon of the prmary user s Cp) and the utlty functon of each secondary user s u p ), {1,, n}. Output: The optmal spectrum allocaton p for each secondary user. 1: t = // t counts the number of teratons. 2: The prmary user sets the same ntal value q ) = θ sent to all secondary users, where θ s a lttle postve real value. 3: Repeat for t =, 1,. 4: For each secondary user, the value of p t+1) computed accordng to p t+1) q ) = arg max u p ) q p. p [,δ ] 5: For the prmary user, accordng to the returned p t+1), the value of q t+1) s updated as q t+1) = 1 ε) q t) + ε C p) p p= p t+1). 6: Untl f p t+1) p t) ζ, where ζ s a tunable lttle real number. 7: p = p t+1), q = q t+1). B. Equvalent optmal prcng problem As mentoned above, the prmary user gves a prce q to each secondary user for leasng spectrum. Accordng to the prce, each secondary user decdes the bandwdth of the spectrum p t demands. Based on the prcng mechansm, the socal welfare can be rewrtten as follows, s

4 4 Fg. 2. SU 1 q 1 p 1 q 2 PU SU 2 p C q 1 q, p q 3 p 2 p 3 SU 3 p q arg maxu p qp p, p p The prcng-based decomposton framework context. The alternate operatons are executed teratvely by the prmary user and secondary users n turn untl the optmal soluton q and p s acheved. Intutvely, the prce updatng rule s sgnfcant to guarantee the teratve operatons eventually converge to the optmal soluton. Here, we propose a prcng updatng rule accordng to the the Karuth-Kuhn-Tucker KKT) condton. Accordng to the convex optmzaton theory [15], the optmal soluton q = [q1,, qn] of 1) must satsfy the KKT condton as shown n the followng proposton. [ Proposton 1. An optmal prce vector q = q 1, n], q should satsfy q = C p) p,, 11) p= p S p) = = = u p ) C p) u p ) p q ) + pq T C p) P q, p ) + Φ p, q) 8) Consderng secondary users are ratonal, they wll choose the bandwdth of allocated spectrum by maxmzng ther own payoffs gven a spectrum prce, p q ) = arg max u p ) q p. 9) p [,δ ] We denote p q ) as p for short n the next context. Therefore, the socal welfare maxmzaton problem 6) can be transformed nto an equvalent optmal prcng problem, whch determnes the optmal prce vector q = [ q 1,, q n] by maxmzng the socal welfare as follows, max q n Φ p, q) + P q, p ) 1) We llustrate why the prcng problem s equvalent to the orgnal problem brefly. Assume that the optmal soluton of problem 1) s q. The secondary user chooses p q ) dependng on the optmal prce q by maxmzng the payoff 9). After knowng p q ) and q, the socal welfare s maxmzed dependng on 8). Therefore, the optmal soluton of p can be obtaned when fndng the optmal prce vector q = [q 1,, q n]. C. Prcng-based decomposton framework In ths subsecton, we ntroduce a prcng-based decomposton framework, whch s shown n Fg.2. Each secondary user computes the value of p as shown n 9) locally and sends the result to the prmary user. Accordng to the collected value of p, the prmary user updates the prce q accordng to a updatng rule, whch s gven n the comng where p = [ p 1,, p n ] T wth p = p q ) gven n 9). Proof: Frst, accordng to the frst order condton, there exsts P p, q )/ p = for 2). Therefore, we obtan q = u p ). 12) Moreover, applyng the KKT condton S p)/ q = to ) 5), there s u p ) Cp) p p q p= p = snce p q <, we have u C p) p ) p =. 13) p = p The proposton s proved by combnng the 12) and 13). Gven the nsght of Proposton 1, we propose a prce updatng rule as follows, q = 1 ε) q +ε C p) p, 14) p= p, where ε, 1] s a tunable parameter. Note that n our decomposton framework, updatng prce s unnecessary to be executed after all p, {1,, n} are returned. The prmary user can update q for secondary user SU arbtrary tmes before q j j ) s updated. It means that our decomposton framework can be realzed n an asynchronous way so that the number of teratons can be largely degraded as dfferent secondary users have dfferent communcaton delay wth the prmary user. D. Dstrbuted Algorthm In ths subsecton, we propose a dstrbuted algorthm to maxmze the socal welfare of a cogntve rado network. In the dstrbuted algorthm, the prmary user and all secondary users should partcpate as follows, For each secondary user SU, t maxmzes the payoff P p, q ) = u p ) q p locally to obtan the value of p based on the prce q gven by the prmary user. Each secondary user only needs to know ts own utlty functon. For the prmary user, t updates the value of q based on p returned from each secondary user. The prce updatng rule s gven n 14).

5 5 The prmary user generates a new prce based on an orgnal prce and partal dervatve Cp) p, whch pushes the prce p= p q towards the optmal value q = Cp ) p. The detals of the dstrbuted algorthm are descrbed n algorthm 1. E. Convergence Analyss Contracton Mappng. Many teratve algorthms can be expressed as x l + 1) = φx l)), l =, 1,, where x ) X and l denotes the number of teratons. Mappng φ s called a contracton f φ x) φ y) κ x y, x, y X, 15) where s a norm and κ [, 1) s called a modulus of φ. Moreover, the mappng φ s called a pseudo-contracton f there exsts a fxed pont x X means x φ x )) and φ x) x x x, x X. 16) The convergence property of contracton or pseudocontractons s gven n Theorem 1. Theorem 1 Geometrc Convergence). Supposng the mappng or the pseudo-contracton and the modulus of φ s κ [, 1). { Then, φ has an} unque fxed pont x and a sequence xl), l =, 1, generated by x l + 1) = φ x l)) satsfes x l) x κ l x ) x, l, 17) for every choce of ntal x ) X. In partcular, x l) converges to x geometrcally. In ths secton, we analyze the convergence of our algorthm under the condton of ε = 1 and ε < 1 n Proposton 2 and Proposton 3, respectvely. For smplcty, the second order partal dervatves of Cp) s denoted by pjp C p) = Cp) p j p. Convergence of Algorthm 1 wth ε = 1. We frst defne a notaton [λ] + to denote the projecton of λ R onto the range [, τ ], [λ] + = arg max z=[, τ ] z λ. Brefly, a soluton to 9) s equvalent to p q ) = [arg max p u p ) q p ] + = [ u 1 q ) ] +. Therefore, when ε = 1, the prce updatng rule turns to q = c p). Substtutng t nto p q ) obtans p = φ p) = [ u 1 c p)) ] +. 18) Algorthm 1 appled wth ε = 1 eventually acheves convergence under a certan condton whch s gven n the followng proposton. Proposton 2. Supposng ε = 1, f we have pj p C p) < mn u p ), p p j=1 [, δ ], 19) )} { p q t) generated by Algorthm 1 converges geometrcally [16] to the optmal soluton p of 6), gvng an ntal value q ). Proof: For each, a functon g s defned as follows. g r) = u 1 c z r))) = u 1 c rx + 1 r) y)), where g r) s dfferentable and r [, 1]. We have φ x) φ y) = [ u 1 c x)) ] + [ u 1 c y)) ] + 1 g 1) g ) = dg r) dr 1 dg r) dr max dg r) r [,1] dr, where the frst nequalty s because [x ] + [y ] + x y for all x y R. Furthermore, applyng a chan rule, we have dg r) dr = n j=1 jc 1 c rx + 1 r) y) x j y j )) u 1 ) c z r))) n j=1 jc z r)) x j y j. If condton 19) holds, we have j=1 j c x) < mn u x )) u x u 1 c x)) ) for all x [,τ ]. Therefore, there exsts a number κ < 1 whch enables that dg r) dr κ max x j y j = κ x y j, r [, 1] Therefore, we have φ x) φ y) κ x y, x, y [, τ ], whch shows φ s a contracton wth modulus κ respect to a maxmum norm. Therefore, a proposton s proved. Convergence of Algorthm 1 wth ε < 1. We relax a value of ε by consderng a convergence condton wth ε < 1. In ths stuaton, we suppose that p q l)) generated wth each teraton l s always wthn a range [,τ ]. Therefore, we can rewrte a soluton to 9) as p q ) = u 1 q ), whch has a lttle dfference from ε = 1 stuaton. Combnng t wth q = 1 ε) q + εc p), we can obtan p = φ p) = u 1 1 ε) u p ) + εc p)) 2) Appled wth 2), the Algorthm 1 can acheve optmal under a convergence condton n the followng proposton. Proposton 3. Set p q ) ) to be for all or δ for all. If we have ) 1 ε 1 + p j p Cp) c p ) j p j p Cp), p [, δ ], 21) )} { p q t) generated by our algorthm 1 converges geometrcally [16] to the optmal soluton p of 6).

6 6 Proof: For each a functon g s defned as follows. g r) = u 1 1 ε) u z r))) + εc z r)),, where z r) = rx + 1 r) y. Supposng x < x, we wll show x < φ x) x. We have φ x) φ x ) = g 1) g ) = where dg r) dr s gven by 22). 1 dg r) dr, dr dg r) dr = u 1 ) 1 ε) u z r)) + εc z r)) ) 1 ε) u z r)) x x ) + ε j c z r)) x ) j x j j = u φ z r)))) 1 1 ε) u z r)) + ε c z r))) x x ) ) +ε j jc z r)) x ) j x j 22) Applyng h 1 ) to both sdes yelds x < φ x). Therefore, there exsts a number κ < 1 enables that φ x) x κ x x, whch s equvalent to φ x) x κ x x, x [, x ]. In other words, φ s a pseudo-contracton n [, x ]. Smlarly, [ f x ] > x, we can get φ, whch s a pseudo-contracton n x, τ. Therefore, the proposton s proved. IV. PERFORMANCE EVALUATION A. Methodology and Smulaton Setup We perform smulatons to evaluate the performance of our optmal dstrbuted algorthm, compared wth a baselne algorthm as follows. Dstrbuted weghted allocaton scheme DWAS). Frst, the prmary user provdes an ntal value q for all secondary users. Each secondary user SU returns the value of p q ) by maxmzng P q, p ). Next, the prmary user decdes the fnal value of q, {1,, n} by maxmzng Φq, p), assumng that u p ) s a lnear functon of p,.e., u p ) = b p. b represents the weght of secondary user SU, whch s computed as p q )/q. After recevng q, each secondary user computes ts optmal p q ). In addton to the DWAS, we also compare our algorthm wth the optmal value. The optmal value s obtaned by maxmzng the socal welfare drectly. The evaluaton of our algorthm s performed from two aspects, ncludng the network performance and the convergence speed. The metrcs of network performance are the socal welfare and the total allocated spectrum. We perform the smulatons wth varyng two factors: the number of secondary users and the value of β. The metrc of convergence speed s the number of teratons. The smulatons are performed under dfferent numbers of secondary users and values of ε. The ntal values of q, {1, 2,, n} n our proposed algorthm and the DWAS are set to 1. The default values of β and ε are.2 and.1, respectvely. The values of a are standard unformly dstrbuted n the open nterval.5, 1). Fg. 3. Socal welfare Algorthm 1ε =.1) Algorthm 1ε =.15) Number of teratons Socal welfare acheved wth the number of teratons B. Impacts of the number of secondary users We frst study the socal welfare and allocated spectrum of our algorthm and the DWAS under the condton of dfferent numbers of secondary users, compared wth the optmum. The results are shown n Fg. 4 and Fg. 7. Fg. 4 shows that the socal welfare obtaned by our algorthm performs as well as the optmum, and much better than the DWAS. The socal welfare obtaned by our algorthm and the DWAS ncreases wth the ncreasng number of secondary users. When there are 2 secondary users, the socal welfare obtaned by our algorthm s 5.34% hgher than the DWAS. In Fg. 7, the spectrum allocated by our algorthm s the same as the optmum, and more than the DWAS. The spectrum obtaned by our algorthm and the DWAS ncreases wth the ncreasng number of secondary users. However, the spectrum allocated by DWAS s always less than the spectrum allocated by our algorthm. When there are 2 secondary users, the allocated spectrum obtaned by our algorthm s 96.23% more than the DWAS. C. Effects of the value of β We next study the performance of our algorthm and the DWAS under the condton of dfferent values of β, compared wth the optmum. The results are shown n Fg. 5 and Fg. 8. Fg. 5 shows that the socal welfare acheved by our algorthm performs as well as the optmum, and more than the DWAS. The socal welfare obtaned by our algorthm and the DWAS decreases wth the ncreasng value of β. When there are 2 secondary users, the socal welfare obtaned by our algorthm s 4% more than the DWAS. In Fg. 8, the spectrum allocated by our algorthm s the same as the optmum, and more than the DWAS. For our algorthm and the DWAS, the allocated spectrum decreases, when the value of β ncreases. When there are 2 secondary users, the allocated spectrum obtaned by our algorthm s 97.43% more than the DWAS. D. Convergence Speed In ths subsecton, we study the convergence speed of our algorthm under dfferent numbers of secondary users and values of ε. The results are shown n Fg. 3, Fg. 6 and Fg. 9.

7 7 Socal welfare Optmum Algorthm 1 DWAS Socal welfare Optmum Algorthm 1 DWAS Number of teratons n = 2 n = 4 n = 6 n = 8 n = Number of secondary users Value of β varepslon =.1 varepslon =.15 Number of secondary users Fg. 4. Socal welfare comparson of dfferent Fg. 5. Socal welfare comparson of dfferent Fg. 6. Convergence of Algorthm 1 wth dfferent algorthms wth varyng number of secondary users algorthms wth varyng β number of secondary users Allocated spectrum Algorthm 1 Optmum DWAS Allocated spectrum Optmum Algorthm 1 DWAS Number of teratons Algorthm 1 wth 6 secondary users Algorthm 1 wth 1 secondary users Number of secondary users Value of β Value of ε Fg. 7. Allocated spectrum comparson of dfferent Fg. 8. Allocated spectrum comparson of dfferent Fg. 9. Convergence of Algorthm 1 wth dfferent algorthms wth varyng number of secondary users algorthms wth varyng β value of parameter ε Fg. 3 plots the socal welfare acheved by our algorthm n each teraton, when the value of system parameter ε s.1 and.15, respectvely. Under the two settngs, we can see that 99.9% of the optmal socal welfare s acheved n the frst ten teratons, whch mples our algorthm converges quckly. In Fg. 6, We study smulatons of our algorthm when the number of secondary users s 2, 4, 6, 8 and 1, respectvely. The value of parameter ε s.1 and.15. The numbers of teratons under dfferent numbers of secondary users and dfferent values of ε are close, whch llustrates our algorthm s adapted to large-scale cogntve rado networks. Fg. 9 shows that the convergence speed of our algorthm under the condton of the value of ε varyng from.1 to.15. The number of secondary users s 6 and 1. The number of teratons under dfferent values of ε s less than 32. The hgher value of system parameter ε, the faster the convergence speed. Our algorthm converges fast under varyng values of ε. V. RELATED WORK Increasng works focus on the problems n cogntve rado networks. Some works focus on estmaton of nterference temperature of the rado envronment and detecton of spectrum holes. such as [17] [18]. Works n [19] [2] [21] focus on estmaton of channel-state nformaton and predcton of channel capacty for usng by the transmtter. Some other works [5] [6] focus on transmt-power control and dynamc spectrum management. The spectrum allocaton n cogntve rado networks has been extensvely studed. By takng advantage of the models n game theory, some of exstng works solve the spectrum sharng problem by achevng Nash Equlbrum. Nyato et al. [7] formulate the problem of spectrum sharng as an olgopoly market competton and use a Cournot game to obtan the spectrum allocaton. However, the statc Cournot game only adapts to some specal envronment. Han et al. [9] propose a correlated equlbrum concept for users to have the dstrbutve opportunstc spectrum access. Wang et al. [8] propose a compettve spectrum-sharng scheme accordng to aucton theorem wthout achevng socal welfare maxmzaton. Some centralzed algorthms n some other works can actually acheve the socal welfare maxmzaton. However, these algorthms could not protect the prvacy of each user. Iler et al. [12] propose a framework under the regulaton of a spectrum polcy server for spectrum sharng. Peng et al. [1] propose a framework that defnes the spectrum sharng problem for some defntons of overall system utlty. The central sever s desgned for allocaton assgnments. Raman et al. [11] propose a centralzed spectrum server that coordnates users to sharng a common spectrum.

8 8 VI. CONCLUSION The problem of maxmzng socal welfare for spectrum sharng n cogntve rado networks s studed n ths paper. We frst defne the concept of socal welfare by formulatng the cost of the prmary user and the utlty of each secondary user. To protect the prvate nformaton of secondary users, we propose a dstrbuted algorthm based on a prcng-based decomposton framework. The dstrbuted algorthm s executed teratvely by the prmary user and each secondary user n turn. Theoretcal analyss s conducted to demonstrate that the algorthm converges to the optmal soluton. Extensve smulaton results show that our proposed acheves maxmal socal welfare and converges quckly. REFERENCES [1] P. Gupta and P. R. Kumar, The capacty of wreless networks, Informaton Theory, IEEE Transactons on, vol. 46, no. 2, pp , 2. [2] J. M. III and Q. G. J. Magure, Cogntve rado: makng software rados more personal, Personal Communcatons, IEEE, vol. 6, no. 4, pp , [3] S. Haykn, Cogntve rado: bran-empowered wreless communcatons, Selected Areas n Communcatons, IEEE Journal on, vol. 23, no. 2, pp , 25. [4] Q. Zhao and B. M. Sadler, A survey of dynamc spectrum access, Sgnal Processng Magazne, IEEE. [5] W. Yu, Competton and cooperaton n mult-user communcaton envronments, Ph.D. dssertaton, 22. [6] J. Zander, S. Km, M. Almgren, and O. Queseth, Rado resource management for wreless networks. Artech House, Inc., 21. [7] a. E. H. D. Nyato, A game-theoretc approach to compettve spectrum sharng n cogntve rado networks, n Wreless Communcatons and Networkng Conference, 27. [8] X. Wang, Z. L, P. c. Xu, Y. y. Xu, X. b. Gao, and H. H. Chen, Spectrum sharng n cogntve rado networksłan aucton-based approach, Systems, Man, and Cybernetcs, Part B: Cybernetcs, IEEE Transactons on, vol. 4, no. 3, pp , 21. [9] Z. Han and a. K. R. L. C. Pandana, Dstrbutve opportunstc spectrum access for cogntve rado usng correlated equlbrum and no-regret learnng, n Wreless Communcatons and Networkng Conference, 27. [1] C. Peng, H. t. Zheng, and B. Y. Zhao, Utlzaton and farness n spectrum assgnment for opportunstc spectrum access, Moble Networks and Applcatons, vol. 11, no. 4, pp , 26. [11] C. Raman, R. Yates, and N. B. Mandayam, Schedulng varable rate lnks va a spectrum server, n New Fronters n Dynamc Spectrum Access Networks, 25. [12] O. Iler, D. Samardzja, and N. B. Mandayam, Demand responsve prcng and compettve spectrum allocaton va a spectrum server, n New Fronters n Dynamc Spectrum Access Networks, 25. [13] A. Chateauneuf and M. Cohen, Rsk seekng wth dmnshng margnal utlty n a non-expected utlty model, Journal of Rsk and Uncertanty, vol. 9, no. 1, pp , [14] S. P. Boyd and L. Vandenberghe, Convex optmzaton. Cambrdge unversty press, 24. [15] D. P. Bertsekas, A. Ned, and A. E. Ozdaglar, Convex analyss and optmzaton. Athena Scentfc Belmont, 23. [16] K. H. Khall, Nonlnear systems. Prentce hall Upper Saddle Rver, 22, vol. 3. [17] S. Haykn, Communcaton systems, New York, 21. [18] J. D. Thomson, Spectrum estmaton and harmonc analyss, Proc. IEEE INFOCOM, vol. 7, no. 9, pp , [19] S. Haykn, K. Huber, and Z. Chen, Bayesan sequental state estmaton for mmo wreless communcatons, Proc. IEEE INFOCOM, vol. 92, no. 3, pp , 24. [2] B. M. Hochwald, T. L. Marzetta, and V. Tarokh, Multple-antenna channel hardenng and ts mplcatons for rate feedback and schedulng, Informaton Theory, IEEE Transactons on, vol. 5, no. 9, pp , 24. [21] G. J. Foschn and M. J. Gans, On lmts of wreless communcatons n a fadng envronment when usng multple antennas, Wreless personal communcatons, vol. 6, no. 3, pp , 1998.

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

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

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

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

Perfect Competition and the Nash Bargaining Solution

Perfect Competition and the Nash Bargaining Solution Perfect Competton and the Nash Barganng Soluton Renhard John Department of Economcs Unversty of Bonn Adenauerallee 24-42 53113 Bonn, Germany emal: rohn@un-bonn.de May 2005 Abstract For a lnear exchange

More information

Externalities in wireless communication: A public goods solution approach to power allocation. by Shrutivandana Sharma

Externalities in wireless communication: A public goods solution approach to power allocation. by Shrutivandana Sharma Externaltes n wreless communcaton: A publc goods soluton approach to power allocaton by Shrutvandana Sharma SI 786 Tuesday, Feb 2, 2006 Outlne Externaltes: Introducton Plannng wth externaltes Power allocaton:

More information

(1 ) (1 ) 0 (1 ) (1 ) 0

(1 ) (1 ) 0 (1 ) (1 ) 0 Appendx A Appendx A contans proofs for resubmsson "Contractng Informaton Securty n the Presence of Double oral Hazard" Proof of Lemma 1: Assume that, to the contrary, BS efforts are achevable under a blateral

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

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

Lecture 10 Support Vector Machines II

Lecture 10 Support Vector Machines II Lecture 10 Support Vector Machnes II 22 February 2016 Taylor B. Arnold Yale Statstcs STAT 365/665 1/28 Notes: Problem 3 s posted and due ths upcomng Frday There was an early bug n the fake-test data; fxed

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

Lecture 14: Bandits with Budget Constraints

Lecture 14: Bandits with Budget Constraints IEOR 8100-001: Learnng and Optmzaton for Sequental Decson Makng 03/07/16 Lecture 14: andts wth udget Constrants Instructor: Shpra Agrawal Scrbed by: Zhpeng Lu 1 Problem defnton In the regular Mult-armed

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

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

Infinitely Split Nash Equilibrium Problems in Repeated Games

Infinitely Split Nash Equilibrium Problems in Repeated Games Infntely Splt ash Equlbrum Problems n Repeated Games Jnlu L Department of Mathematcs Shawnee State Unversty Portsmouth, Oho 4566 USA Abstract In ths paper, we ntroduce the concept of nfntely splt ash equlbrum

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

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

An Admission Control Algorithm in Cloud Computing Systems

An Admission Control Algorithm in Cloud Computing Systems An Admsson Control Algorthm n Cloud Computng Systems Authors: Frank Yeong-Sung Ln Department of Informaton Management Natonal Tawan Unversty Tape, Tawan, R.O.C. ysln@m.ntu.edu.tw Yngje Lan Management Scence

More information

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium?

Welfare Properties of General Equilibrium. What can be said about optimality properties of resource allocation implied by general equilibrium? APPLIED WELFARE ECONOMICS AND POLICY ANALYSIS Welfare Propertes of General Equlbrum What can be sad about optmalty propertes of resource allocaton mpled by general equlbrum? Any crteron used to compare

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

Computing Correlated Equilibria in Multi-Player Games

Computing Correlated Equilibria in Multi-Player Games Computng Correlated Equlbra n Mult-Player Games Chrstos H. Papadmtrou Presented by Zhanxang Huang December 7th, 2005 1 The Author Dr. Chrstos H. Papadmtrou CS professor at UC Berkley (taught at Harvard,

More information

COS 521: Advanced Algorithms Game Theory and Linear Programming

COS 521: Advanced Algorithms Game Theory and Linear Programming COS 521: Advanced Algorthms Game Theory and Lnear Programmng Moses Charkar February 27, 2013 In these notes, we ntroduce some basc concepts n game theory and lnear programmng (LP). We show a connecton

More information

Mixed Taxation and Production Efficiency

Mixed Taxation and Production Efficiency Floran Scheuer 2/23/2016 Mxed Taxaton and Producton Effcency 1 Overvew 1. Unform commodty taxaton under non-lnear ncome taxaton Atknson-Stgltz (JPubE 1976) Theorem Applcaton to captal taxaton 2. Unform

More information

PROBLEM SET 7 GENERAL EQUILIBRIUM

PROBLEM SET 7 GENERAL EQUILIBRIUM PROBLEM SET 7 GENERAL EQUILIBRIUM Queston a Defnton: An Arrow-Debreu Compettve Equlbrum s a vector of prces {p t } and allocatons {c t, c 2 t } whch satsfes ( Gven {p t }, c t maxmzes βt ln c t subject

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

k t+1 + c t A t k t, t=0

k t+1 + c t A t k t, t=0 Macro II (UC3M, MA/PhD Econ) Professor: Matthas Kredler Fnal Exam 6 May 208 You have 50 mnutes to complete the exam There are 80 ponts n total The exam has 4 pages If somethng n the queston s unclear,

More information

The Second Anti-Mathima on Game Theory

The Second Anti-Mathima on Game Theory The Second Ant-Mathma on Game Theory Ath. Kehagas December 1 2006 1 Introducton In ths note we wll examne the noton of game equlbrum for three types of games 1. 2-player 2-acton zero-sum games 2. 2-player

More information

Economics 101. Lecture 4 - Equilibrium and Efficiency

Economics 101. Lecture 4 - Equilibrium and Efficiency Economcs 0 Lecture 4 - Equlbrum and Effcency Intro As dscussed n the prevous lecture, we wll now move from an envronment where we looed at consumers mang decsons n solaton to analyzng economes full of

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

The Order Relation and Trace Inequalities for. Hermitian Operators

The Order Relation and Trace Inequalities for. Hermitian Operators Internatonal Mathematcal Forum, Vol 3, 08, no, 507-57 HIKARI Ltd, wwwm-hkarcom https://doorg/0988/mf088055 The Order Relaton and Trace Inequaltes for Hermtan Operators Y Huang School of Informaton Scence

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

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem

Interactive Bi-Level Multi-Objective Integer. Non-linear Programming Problem Appled Mathematcal Scences Vol 5 0 no 65 3 33 Interactve B-Level Mult-Objectve Integer Non-lnear Programmng Problem O E Emam Department of Informaton Systems aculty of Computer Scence and nformaton Helwan

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

General viscosity iterative method for a sequence of quasi-nonexpansive mappings

General viscosity iterative method for a sequence of quasi-nonexpansive mappings Avalable onlne at www.tjnsa.com J. Nonlnear Sc. Appl. 9 (2016), 5672 5682 Research Artcle General vscosty teratve method for a sequence of quas-nonexpansve mappngs Cuje Zhang, Ynan Wang College of Scence,

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

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis

Appendix for Causal Interaction in Factorial Experiments: Application to Conjoint Analysis A Appendx for Causal Interacton n Factoral Experments: Applcaton to Conjont Analyss Mathematcal Appendx: Proofs of Theorems A. Lemmas Below, we descrbe all the lemmas, whch are used to prove the man theorems

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

Structure and Drive Paul A. Jensen Copyright July 20, 2003

Structure and Drive Paul A. Jensen Copyright July 20, 2003 Structure and Drve Paul A. Jensen Copyrght July 20, 2003 A system s made up of several operatons wth flow passng between them. The structure of the system descrbes the flow paths from nputs to outputs.

More information

Conjectures in Cournot Duopoly under Cost Uncertainty

Conjectures in Cournot Duopoly under Cost Uncertainty Conjectures n Cournot Duopoly under Cost Uncertanty Suyeol Ryu and Iltae Km * Ths paper presents a Cournot duopoly model based on a condton when frms are facng cost uncertanty under rsk neutralty and rsk

More information

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists *

How Strong Are Weak Patents? Joseph Farrell and Carl Shapiro. Supplementary Material Licensing Probabilistic Patents to Cournot Oligopolists * How Strong Are Weak Patents? Joseph Farrell and Carl Shapro Supplementary Materal Lcensng Probablstc Patents to Cournot Olgopolsts * September 007 We study here the specal case n whch downstream competton

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

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

Optimal Bandwidth and Power Allocation for

Optimal Bandwidth and Power Allocation for Optmal Bandwdth and Power Allocaton for 1 Sum Ergodc Capacty under Fadng Channels n Cogntve Rado Networks arxv:1006.5061v1 [cs.it] 25 Jun 2010 Xaowen Gong, Student Member, IEEE, Sergy A. Vorobyov, Senor

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

Supporting Information for: Two Monetary Models with Alternating Markets

Supporting Information for: Two Monetary Models with Alternating Markets Supportng Informaton for: Two Monetary Models wth Alternatng Markets Gabrele Camera Chapman Unversty & Unversty of Basel YL Chen St. Lous Fed November 2015 1 Optmal choces n the CIA model On date t, gven

More information

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming

EEL 6266 Power System Operation and Control. Chapter 3 Economic Dispatch Using Dynamic Programming EEL 6266 Power System Operaton and Control Chapter 3 Economc Dspatch Usng Dynamc Programmng Pecewse Lnear Cost Functons Common practce many utltes prefer to represent ther generator cost functons as sngle-

More information

MMA and GCMMA two methods for nonlinear optimization

MMA and GCMMA two methods for nonlinear optimization MMA and GCMMA two methods for nonlnear optmzaton Krster Svanberg Optmzaton and Systems Theory, KTH, Stockholm, Sweden. krlle@math.kth.se Ths note descrbes the algorthms used n the author s 2007 mplementatons

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

DUE to the global aspiration for energy saving and reduction. Relay Power Allocation in Auction-based Game Approach

DUE to the global aspiration for energy saving and reduction. Relay Power Allocation in Auction-based Game Approach Globecom 0 - Communcatons QoS, Relablty and Modellng Symposum Relay Power Allocaton n Aucton-based Game Approach Dan Wu, Yuemng Ca, Lang Zhou, and Joel Rodrgues Insttute of Communcatons Engneerng, PLA

More information

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION

ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EQUATION Advanced Mathematcal Models & Applcatons Vol.3, No.3, 2018, pp.215-222 ON A DETERMINATION OF THE INITIAL FUNCTIONS FROM THE OBSERVED VALUES OF THE BOUNDARY FUNCTIONS FOR THE SECOND-ORDER HYPERBOLIC EUATION

More information

Supporting Materials for: Two Monetary Models with Alternating Markets

Supporting Materials for: Two Monetary Models with Alternating Markets Supportng Materals for: Two Monetary Models wth Alternatng Markets Gabrele Camera Chapman Unversty Unversty of Basel YL Chen Federal Reserve Bank of St. Lous 1 Optmal choces n the CIA model On date t,

More information

APPENDIX A Some Linear Algebra

APPENDIX A Some Linear Algebra APPENDIX A Some Lnear Algebra The collecton of m, n matrces A.1 Matrces a 1,1,..., a 1,n A = a m,1,..., a m,n wth real elements a,j s denoted by R m,n. If n = 1 then A s called a column vector. Smlarly,

More information

Estimation: Part 2. Chapter GREG estimation

Estimation: Part 2. Chapter GREG estimation Chapter 9 Estmaton: Part 2 9. GREG estmaton In Chapter 8, we have seen that the regresson estmator s an effcent estmator when there s a lnear relatonshp between y and x. In ths chapter, we generalzed the

More information

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal Inner Product Defnton 1 () A Eucldean space s a fnte-dmensonal vector space over the reals R, wth an nner product,. Defnton 2 (Inner Product) An nner product, on a real vector space X s a symmetrc, blnear,

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

Lecture 21: Numerical methods for pricing American type derivatives

Lecture 21: Numerical methods for pricing American type derivatives Lecture 21: Numercal methods for prcng Amercan type dervatves Xaoguang Wang STAT 598W Aprl 10th, 2014 (STAT 598W) Lecture 21 1 / 26 Outlne 1 Fnte Dfference Method Explct Method Penalty Method (STAT 598W)

More information

On the Multicriteria Integer Network Flow Problem

On the Multicriteria Integer Network Flow Problem BULGARIAN ACADEMY OF SCIENCES CYBERNETICS AND INFORMATION TECHNOLOGIES Volume 5, No 2 Sofa 2005 On the Multcrtera Integer Network Flow Problem Vassl Vasslev, Marana Nkolova, Maryana Vassleva Insttute of

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

Credit Card Pricing and Impact of Adverse Selection

Credit Card Pricing and Impact of Adverse Selection Credt Card Prcng and Impact of Adverse Selecton Bo Huang and Lyn C. Thomas Unversty of Southampton Contents Background Aucton model of credt card solctaton - Errors n probablty of beng Good - Errors n

More information

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg

princeton univ. F 17 cos 521: Advanced Algorithm Design Lecture 7: LP Duality Lecturer: Matt Weinberg prnceton unv. F 17 cos 521: Advanced Algorthm Desgn Lecture 7: LP Dualty Lecturer: Matt Wenberg Scrbe: LP Dualty s an extremely useful tool for analyzng structural propertes of lnear programs. Whle there

More information

Solving Nonlinear Differential Equations by a Neural Network Method

Solving Nonlinear Differential Equations by a Neural Network Method Solvng Nonlnear Dfferental Equatons by a Neural Network Method Luce P. Aarts and Peter Van der Veer Delft Unversty of Technology, Faculty of Cvlengneerng and Geoscences, Secton of Cvlengneerng Informatcs,

More information

Analysis of Queuing Delay in Multimedia Gateway Call Routing

Analysis of Queuing Delay in Multimedia Gateway Call Routing Analyss of Queung Delay n Multmeda ateway Call Routng Qwe Huang UTtarcom Inc, 33 Wood Ave. outh Iseln, NJ 08830, U..A Errol Lloyd Computer Informaton cences Department, Unv. of Delaware, Newark, DE 976,

More information

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 )

Yong Joon Ryang. 1. Introduction Consider the multicommodity transportation problem with convex quadratic cost function. 1 2 (x x0 ) T Q(x x 0 ) Kangweon-Kyungk Math. Jour. 4 1996), No. 1, pp. 7 16 AN ITERATIVE ROW-ACTION METHOD FOR MULTICOMMODITY TRANSPORTATION PROBLEMS Yong Joon Ryang Abstract. The optmzaton problems wth quadratc constrants often

More information

Equilibrium with Complete Markets. Instructor: Dmytro Hryshko

Equilibrium with Complete Markets. Instructor: Dmytro Hryshko Equlbrum wth Complete Markets Instructor: Dmytro Hryshko 1 / 33 Readngs Ljungqvst and Sargent. Recursve Macroeconomc Theory. MIT Press. Chapter 8. 2 / 33 Equlbrum n pure exchange, nfnte horzon economes,

More information

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals

Simultaneous Optimization of Berth Allocation, Quay Crane Assignment and Quay Crane Scheduling Problems in Container Terminals Smultaneous Optmzaton of Berth Allocaton, Quay Crane Assgnment and Quay Crane Schedulng Problems n Contaner Termnals Necat Aras, Yavuz Türkoğulları, Z. Caner Taşkın, Kuban Altınel Abstract In ths work,

More information

Explicit and Implicit Temperature Constraints in Energy Harvesting Communications

Explicit and Implicit Temperature Constraints in Energy Harvesting Communications Explct and Implct Temperature Constrants n Energy Harvestng Communcatons Abdulrahman Baknna, Omur Ozel 2, and Sennur Ulukus Department of Electrcal and Computer Engneerng, Unversty of Maryland, College

More information

The Minimum Universal Cost Flow in an Infeasible Flow Network

The Minimum Universal Cost Flow in an Infeasible Flow Network Journal of Scences, Islamc Republc of Iran 17(2): 175-180 (2006) Unversty of Tehran, ISSN 1016-1104 http://jscencesutacr The Mnmum Unversal Cost Flow n an Infeasble Flow Network H Saleh Fathabad * M Bagheran

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

Dynamic Chinese Restaurant Game in Cognitive Radio Networks

Dynamic Chinese Restaurant Game in Cognitive Radio Networks 2013 Proceedngs IEEE INFOCO Dynamc Chnese Restaurant Game n Cogntve Rado Networks Chunxao Jang, Yan Chen, Yu-Han Yang, Chh-Yu Wang, and K. J. Ray Lu Department of Electrcal and Computer Engneerng, Unversty

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

Module 17: Mechanism Design & Optimal Auctions

Module 17: Mechanism Design & Optimal Auctions Module 7: Mechansm Desgn & Optmal Auctons Informaton Economcs (Ec 55) George Georgads Examples: Auctons Blateral trade Producton and dstrbuton n socety General Setup N agents Each agent has prvate nformaton

More information

Some modelling aspects for the Matlab implementation of MMA

Some modelling aspects for the Matlab implementation of MMA Some modellng aspects for the Matlab mplementaton of MMA Krster Svanberg krlle@math.kth.se Optmzaton and Systems Theory Department of Mathematcs KTH, SE 10044 Stockholm September 2004 1. Consdered optmzaton

More information

Hila Etzion. Min-Seok Pang

Hila Etzion. Min-Seok Pang RESERCH RTICLE COPLEENTRY ONLINE SERVICES IN COPETITIVE RKETS: INTINING PROFITILITY IN THE PRESENCE OF NETWORK EFFECTS Hla Etzon Department of Technology and Operatons, Stephen. Ross School of usness,

More information

Portfolios with Trading Constraints and Payout Restrictions

Portfolios with Trading Constraints and Payout Restrictions Portfolos wth Tradng Constrants and Payout Restrctons John R. Brge Northwestern Unversty (ont wor wth Chrs Donohue Xaodong Xu and Gongyun Zhao) 1 General Problem (Very) long-term nvestor (eample: unversty

More information

Minimisation of the Average Response Time in a Cluster of Servers

Minimisation of the Average Response Time in a Cluster of Servers Mnmsaton of the Average Response Tme n a Cluster of Servers Valery Naumov Abstract: In ths paper, we consder task assgnment problem n a cluster of servers. We show that optmal statc task assgnment s tantamount

More information

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan.

THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY. William A. Pearlman. References: S. Arimoto - IEEE Trans. Inform. Thy., Jan. THE ARIMOTO-BLAHUT ALGORITHM FOR COMPUTATION OF CHANNEL CAPACITY Wllam A. Pearlman 2002 References: S. Armoto - IEEE Trans. Inform. Thy., Jan. 1972 R. Blahut - IEEE Trans. Inform. Thy., July 1972 Recall

More information

MAXIMIZING THE THROUGHPUT OF CDMA DATA COMMUNICATIONS THROUGH JOINT ADMISSION AND POWER CONTROL

MAXIMIZING THE THROUGHPUT OF CDMA DATA COMMUNICATIONS THROUGH JOINT ADMISSION AND POWER CONTROL MAXIMIZI THE THROUHPUT OF CDMA DATA COMMUICATIOS THROUH JOIT ADMISSIO AD POWER COTROL Penna Orensten, Davd oodman and Zory Marantz Department of Electrcal and Computer Engneerng Polytechnc Unversty Brooklyn,

More information

Generalized Linear Methods

Generalized Linear Methods Generalzed Lnear Methods 1 Introducton In the Ensemble Methods the general dea s that usng a combnaton of several weak learner one could make a better learner. More formally, assume that we have a set

More information

COGNITIVE RADIO NETWORKS BASED ON OPPORTUNISTIC BEAMFORMING WITH QUANTIZED FEEDBACK

COGNITIVE RADIO NETWORKS BASED ON OPPORTUNISTIC BEAMFORMING WITH QUANTIZED FEEDBACK COGNITIVE RADIO NETWORKS BASED ON OPPORTUNISTIC BEAMFORMING WITH QUANTIZED FEEDBACK Ayman MASSAOUDI, Noura SELLAMI 2, Mohamed SIALA MEDIATRON Lab., Sup Com Unversty of Carthage 283 El Ghazala Arana, Tunsa

More information

Errors for Linear Systems

Errors for Linear Systems Errors for Lnear Systems When we solve a lnear system Ax b we often do not know A and b exactly, but have only approxmatons  and ˆb avalable. Then the best thng we can do s to solve ˆx ˆb exactly whch

More information

Auction-Based Spectrum Sharing

Auction-Based Spectrum Sharing Aucton-Based Spectrum Sharng Janwe Huang, Randall Berry and Mchael L. Hong Department of Electrcal and Computer Engneerng Northwestern Unversty, Evanston IL 60208, USA janwe@northwestern.edu, {rberry,mh}@ece.northwestern.edu

More information

Maximizing the number of nonnegative subsets

Maximizing the number of nonnegative subsets Maxmzng the number of nonnegatve subsets Noga Alon Hao Huang December 1, 213 Abstract Gven a set of n real numbers, f the sum of elements of every subset of sze larger than k s negatve, what s the maxmum

More information

Uncertainty in measurements of power and energy on power networks

Uncertainty in measurements of power and energy on power networks Uncertanty n measurements of power and energy on power networks E. Manov, N. Kolev Department of Measurement and Instrumentaton, Techncal Unversty Sofa, bul. Klment Ohrdsk No8, bl., 000 Sofa, Bulgara Tel./fax:

More information

The Study of Teaching-learning-based Optimization Algorithm

The Study of Teaching-learning-based Optimization Algorithm Advanced Scence and Technology Letters Vol. (AST 06), pp.05- http://dx.do.org/0.57/astl.06. The Study of Teachng-learnng-based Optmzaton Algorthm u Sun, Yan fu, Lele Kong, Haolang Q,, Helongang Insttute

More information

Power Allocation/Beamforming for DF MIMO Two-Way Relaying: Relay and Network Optimization

Power Allocation/Beamforming for DF MIMO Two-Way Relaying: Relay and Network Optimization Power Allocaton/Beamformng for DF MIMO Two-Way Relayng: Relay and Network Optmzaton Je Gao, Janshu Zhang, Sergy A. Vorobyov, Ha Jang, and Martn Haardt Dept. of Electrcal & Computer Engneerng, Unversty

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

Embedded Systems. 4. Aperiodic and Periodic Tasks

Embedded Systems. 4. Aperiodic and Periodic Tasks Embedded Systems 4. Aperodc and Perodc Tasks Lothar Thele 4-1 Contents of Course 1. Embedded Systems Introducton 2. Software Introducton 7. System Components 10. Models 3. Real-Tme Models 4. Perodc/Aperodc

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

Mixed Strategy Nash Equilibrium in Two-User Resource Allocation Games

Mixed Strategy Nash Equilibrium in Two-User Resource Allocation Games IEEE Internatonal Symposum on Informaton Theory Proceedngs Mxed Strategy Nash Equlbrum n Two-User Resource Allocaton Games Je Gao, Sergy A. Vorobyov, and Ha Jang Dept. of Electrcal and Computer Engneerng,

More information

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

On Spatial Capacity of Wireless Ad Hoc Networks with Threshold Based Scheduling On Spatal Capacty of Wreless Ad Hoc Networks wth Threshold Based Schedulng Yue Lng Che, Ru Zhang, Y Gong, and Lngje Duan Abstract arxv:49.2592v [cs.it] 9 Sep 24 Ths paper studes spatal capacty n a stochastc

More information

arxiv: v1 [math.ho] 18 May 2008

arxiv: v1 [math.ho] 18 May 2008 Recurrence Formulas for Fbonacc Sums Adlson J. V. Brandão, João L. Martns 2 arxv:0805.2707v [math.ho] 8 May 2008 Abstract. In ths artcle we present a new recurrence formula for a fnte sum nvolvng the Fbonacc

More information

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm

Design and Optimization of Fuzzy Controller for Inverse Pendulum System Using Genetic Algorithm Desgn and Optmzaton of Fuzzy Controller for Inverse Pendulum System Usng Genetc Algorthm H. Mehraban A. Ashoor Unversty of Tehran Unversty of Tehran h.mehraban@ece.ut.ac.r a.ashoor@ece.ut.ac.r Abstract:

More information

f(x,y) = (4(x 2 4)x,2y) = 0 H(x,y) =

f(x,y) = (4(x 2 4)x,2y) = 0 H(x,y) = Problem Set 3: Unconstraned mzaton n R N. () Fnd all crtcal ponts of f(x,y) (x 4) +y and show whch are ma and whch are mnma. () Fnd all crtcal ponts of f(x,y) (y x ) x and show whch are ma and whch are

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

Second Order Analysis

Second Order Analysis Second Order Analyss In the prevous classes we looked at a method that determnes the load correspondng to a state of bfurcaton equlbrum of a perfect frame by egenvalye analyss The system was assumed to

More information

Adaptive power control algorithm in cognitive radio based on game theory

Adaptive power control algorithm in cognitive radio based on game theory IET Communcatons Research Artcle Adaptve power control algorthm n cogntve rado based on game theory ISSN 1751-8628 Receved on 23rd November 2014 Revsed on 19th January 2015 Accepted on 24th February 2015

More information

A Game Theoretic Approach to Distributed Opportunistic Scheduling

A Game Theoretic Approach to Distributed Opportunistic Scheduling A Game Theoretc Approach to Dstrbuted Opportunstc Schedulng Albert Banchs, Senor Member, IEEE, Andres Garca-Saavedra, Pablo Serrano, Member, IEEE, and Joerg Wdmer, Senor Member, IEEE. Banchs et al.: A

More information

Arizona State University

Arizona State University SCHEDULING AND POWER ALLOCATION TO OPTIMIZE SERVICE AND QUEUE-WAITING TIMES IN COGNITIVE RADIO UPLINKS By arxv:1601.00608v1 [cs.it] 4 Jan 2016 Ahmed Emad Ewasha Commttee: Dr. Chan Tepedelenloğlu, Char

More information

An Interactive Optimisation Tool for Allocation Problems

An Interactive Optimisation Tool for Allocation Problems An Interactve Optmsaton ool for Allocaton Problems Fredr Bonäs, Joam Westerlund and apo Westerlund Process Desgn Laboratory, Faculty of echnology, Åbo Aadem Unversty, uru 20500, Fnland hs paper presents

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