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

Size: px
Start display at page:

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

Transcription

1 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 Unversty, Morgantown, WV, U.S.A. E-mal: and arxv: v1 [cs.it] 14 Sep 013 Abstract Ths paper nvestgates the problem of adaptve power allocaton for dstrbuted best lnear unbased estmaton BLUE of a random parameter at the fuson center FC of a wreless sensor network WSN. An optmal power-allocaton scheme s proposed that mnmzes the L -norm of the vector of local transmt powers, gven a maxmum varance for the BLUE estmator. Ths scheme results n the ncreased lfetme of the WSN compared to smlar approaches that are based on the mnmzaton of the sum of the local transmt powers. The lmtaton of the proposed optmal power-allocaton scheme s that t requres the feedback of the nstantaneous channel state nformaton CSI from the FC to local sensors, whch s not practcal n most applcatons of large-scale WSNs. In ths paper, a lmted-feedback strategy s proposed that elmnates ths requrement by desgnng an optmal codebook for the FC usng the generalzed Lloyd algorthm wth modfed dstorton metrcs. Each sensor amplfes ts analog nosy observaton usng a quantzed verson of ts optmal amplfcaton gan, whch s receved by the FC and used to estmate the unknown parameter. Index Terms Lmted feedback, best lnear unbased estmator BLUE, generalzed Lloyd algorthm, L -norm, power allocaton, dstrbuted estmaton, parameter estmaton, fuson center, wreless sensor networks. I. INTRODUCTION Dstrbuted estmaton s a technology that enables a wde range of wreless sensor network WSN applcatons, such as event detecton, classfcaton, and object trackng [1 6]. In a WSN performng dstrbuted estmaton, the frst step s for the spatally dstrbuted sensors to locally process ther nosy observatons that are correlated wth an unknown parameter to be estmated. Each sensor ether transmts ts analog local observatons usng an amplfy-and-forward strategy [1 4] or sends a quantzed verson of ts local observatons to the fuson center FC [4 6]. In ths paper, we wll consder the former approach due to ts smplcty and practcal feasblty and wll concentrate on the best lnear unbased estmaton BLUE of an unknown random parameter at the FC. In order to fnd the BLUE estmator of the unknown parameter, the FC combnes lnearly processed, nosy observatons of local sensors receved through orthogonal channels corrupted by fadng and addtve Ths work was supported n part by the Offce of Naval Research under Award No. N The contrbutons of M. Fanae and M.C. Valent were sponsored n part by the Natonal Scence Foundaton under Award No. CNS The work of M. Fanae s sponsored n part by the Natonal Scence Foundaton under Award No. IIA Gaussan nose. Ths paper wll address one of the man ssues n the case of analog amplfy-and-forward local processng, whch s fndng the optmal local amplfcaton gans [1 4]. The values of these gans set the nstantaneous transmt power of sensors; therefore, we refer to ther determnaton as the optmal power allocaton to sensors. Cu et al. [] have proposed an optmal power-allocaton scheme to mnmze the sum of the local transmt powers, gven a maxmum estmaton dstorton defned as the varance of the BLUE estmator of a random scalar parameter at the FC of a WSN. Although optmal wth respect to the total transmt power n the network, ths strategy could result n assgnng very hgh transmt powers to sensors wth hgh qualty observatons and less nosy channels, whle assgnng zero power to other sensors. The drect consequence of such power allocaton s that some sensors wll de quckly, whch could n turn result n a network partton, whle the remanng sensors have ether low observaton qualty or too nosy communcaton channels. In order to allevate ths drawback, we propose an adaptve power-allocaton strategy that mnmzes thel -norm of the local transmt power vector, gven a maxmum estmaton dstorton as defned above. Ths approach prevents the assgnment of hgh transmt powers to sensors by puttng a hgher penalty on them, whch n tself reduces the chances of those sensors dyng and the network becomng parttoned. Furthermore, the total transmt power used n the entre network stll stays bounded. As t wll be seen n the next sectons, the optmal local amplfcaton gans found based on the proposed power-allocaton scheme depend on the nstantaneous fadng coeffcents of the channels between the sensors and FC, as s the case n []. Therefore, the FC must feed the exact channel fadng gans back to sensors through nfnte-rate, error-free lnks. Ths requrement s not practcal n most WSN applcatons, especally when the number of sensors n the network s large. In the remander of ths paper, we propose a lmted-feedback strategy to allevate ths requrement. The proposed approach s based on desgnng an optmal codebook usng the generalzed Lloyd algorthm wth modfed dstorton functons, whch s used to quantze the space of the optmal power-allocaton vectors used by the sensors to set ther local amplfcaton gans. In our prevous work [7], we have addressed the same drawback of the power-allocaton scheme proposed n [].

2 Fg. 1: System model of a WSN n whch the FC fnds an estmate of θ. In summary, the man contrbutons of ths paper are as follows: An adaptve power-allocaton scheme s proposed to mnmze the L -norm of the local transmt power vector, gven a maxmum estmaton dstorton at the FC. Ths scheme allevates the problem of assgnng very hgh transmt powers to some sensors, whle turnng off the other ones. Furthermore, a lmted-feedback strategy s proposed to quantze the vector space of the optmal local amplfcaton gans. Approprate dstorton functons are defned for the applcaton of the generalzed Lloyd algorthm n the doman of adaptve power allocaton for dstrbuted estmaton. The rest of ths paper s organzed as follows: In Secton II, the system model of the WSN under study s descrbed. The proposed adaptve power-allocaton strategy s derved n Secton III. A bref dscusson on the motvaton for and mplementaton of the lmted feedback for the proposed power-allocaton scheme s presented n Secton IV. Detals of the mplementaton of the proposed lmted-feedback scheme are dscussed n Secton V. Secton VI provdes the numercal results to show the applcablty of the proposed schemes. Fnally, the paper s concluded n Secton VII. II. SYSTEM MODEL Consder a WSN composed of K spatally dstrbuted sensors, as depcted n Fg. 1. The goal of the WSN s to relably estmate an unknown random parameter θ at ts fuson center FC usng lnearly amplfed versons of local nosy observatons receved through parallel orthogonal coherent channels corrupted by fadng and addtve Gaussan nose. An example of the unknown parameter to be estmated can be the ntensty of the sgnal broadcast by an energy-emttng source and sensed by a set of locally dstrbuted sgnal detectors. Ths estmated varable along wth the propagaton model of the gven sgnal n the observaton envronment could then be used to estmate the locaton of the source. It s assumed that θ has zero mean and unt power, and s otherwse unknown. Suppose that the local nosy observaton at each sensor s a lnear functon of the unknown random parameter as x h θ +n, 1,,...,K, 1 where h s the fxed local observaton gan of sensor, known at the sensor and FC, and n s the spatally ndependent and dentcally dstrbuted..d. addtve observaton nose wth zero mean and known varance σo. Note that no further assumpton s made on the dstrbuton of the random parameter to be estmated and that of the observaton nose. We defne the observaton sgnal-to-nose rato SNR at sensor as β h σ, where denotes the absolute-value operaton. o We assume that there s no nter-sensor communcaton and/or collaboraton among spatally dstrbuted sensors. Each sensor uses an amplfy-and-forward scheme to amplfy ts local nosy observaton before sendng t to the FC as z a x a h θ +a n, 1,,...,K, where z s the sgnal transmtted from sensor to the FC and a s the local amplfcaton gan at sensor. Note that the nstantaneous transmt power of sensor can be found as P a h +σo a σo 1+β. 3 As t can be seen n 3, the value of the local amplfcaton gan at each sensor determnes the nstantaneous transmt power allocated to that sensor. Therefore, we wll call any strategy that assgns a set of local amplfcaton gans to sensors a power-allocaton scheme. All locally processed observatons are transmtted to the FC through orthogonal fadng channels. The receved sgnal from sensor at the FC can be descrbed as y g z +w, 1,,...,K, 4 where g s the multplcatve fadng coeffcent of the channel between sensor and the FC, and w s the spatally ndependent and dentcally dstrbuted addtve Gaussan nose wth zero mean and varanceσc. We assume that the FC can relably estmate the fadng coeffcent of the channel between each sensor and tself. Note that n the above model, we have also assumed that each sensor s synchronzed wth the FC. We defne the channel sgnal-to-nose rato of the sgnal receved from sensor as γ g σ. c III. OPTIMAL POWER ALLOCATION WITH MINIMAL L -NORM OF TRANSMIT-POWER VECTOR Gven a power-allocaton scheme and a realzaton of the fadng gans, the FC combnes the set of receved sgnals from dfferent sensors to fnd the best lnear unbased estmator BLUE for the unknown parameter θ as [8, Chapter 6] K 1 h θ a g K h a g y a g σ o +σc a g σ o +σc, 5 where the correspondng estmator varance can be found as K 1 Var θ a,g h a g a g σ o +σ c K 1 β γ a σ o 1+γ a, 6 σ o

3 n whch a [a 1,a,...,a K ] T and g [g 1,g,...,g K ] T are column vectors contanng the set of local amplfcaton gans a and fadng coeffcents of the channels g, respectvely. A goal of ths paper s to fnd the optmal local amplfcaton gans or equvalently, the optmal power-allocaton scheme that mnmzes the L -norm of the vector of local transmt powers defned as P [P 1,P,...,P K ] T, gven a constrant on the varance of the estmate as defned n 6. Ths objectve can be formulated as the followng convex optmzaton problem: mnmze {P } K K 7 subject to Var θ a,g By replacng P and Var θ a,g from Equatons 3 and 6, respectvely, Equaton 7 s converted to the followng form, whose optmzaton varables are the local amplfcaton gans: [ mnmze a σ {a } K o 1+β ] 8 β γ a subject to σ o 1+γ a 1 σ o Let b be defned as b βγa σ o 1+γ a. The above constraned σ o optmzaton problem could be re-wrtten as b 1+β mnmze {b } K β b γ σo 9 subject to b 1 AN b < β whch s a convex optmzaton problem n terms of b. The Lagrangan functon for ths optmzaton problem s Lb,λ 0,µ P 1 b 1+β β b γ σo 1 +λ 0 b µ b, 10 where b [b 1,b,...,b K ] T s the column vector of target optmzed varables and µ [µ 1,µ,...,µ K ] T s the Lagrangan multpler vector. The Karush-Kuhn-Tucker KKT condtons for ths optmzaton problem can be wrtten as Lb,λ 0,µ β b 1+β b β b 3 λ 0 µ 0, 11a γ σ4 o b 1, 11b µ b 0, 1,,...,K, 11c µ 0 and b 0, 1,,...,K. 11d It can be shown that the cubc equaton defned n 11a only ALGORITHM I: The water-fllng-based teratve process to fnd the unque values for the number of actve sensors K 1 and the constant λ 0. Requre: K, {β } K, and {γ } K. 1. Intalzaton. for 1,,...,K do 3. δ 1+β β γ 4. end for 5. Sort the sensors based on the ascendng values of δ so that δ 1 δ δ K. 6. K 1 K 7. EndIntalzaton 8. repeat 9. Usng the gven value for K 1, fnd the value of λ 0 by solvng Replace the value of λ 0 n Eq. 1 and fnd the new values of b, 1,,...,K. 11. K 1 K untl The values of b do not change from the prevous teraton. In partcular, b > 0 for all K 1, and b 0 for all > K return K 1 and λ 0. has a unque real root, whch s n the nterval 0 < b < β as β b β 1 3 δ T 1 3 β δ + λ 0 3 λ 0 T, 1 where δ 1+β β γ, 1,,...,K, T s defned as T β δ 7λ 0, 13 and the operator [ ] + s defned such that [x] + x f x > 0, and [x] + 0 f x 0. Note that n dervng 1, the complementary slackness requrement 11c s used based on whch µ 0 when b > 0, and b 0 when µ > 0. It should also be noted that as the observaton SNR β or channel SNR γ decreases, the value of δ ncreases, whch n turn ncreases the value of T and decreases the value of b. Therefore, f the sensors are sorted so that δ 1 δ δ K, only the frst K 1 sensors wth the least values of δ wll have a postve value for b, and b 0 for all > K 1. The values of the number of actve sensors K 1 for whch b > 0, and the equalty-constrant Lagrangan multpler λ 0 are unque and can be found by replacng b from 1 nto 11b to derve the followng relatonshp between them: K 1 3 β β δ T λ β δ λ 0 T K 1 β The values of K 1 and λ 0 can be found through the waterfllng-based teratve process summarzed n Algorthm I. It can be shown that the soluton of the above teratve algorthm n terms of K 1 and λ 0 always exsts and s unque. Havng found b through the above process, the local amplfcaton gan a can be found as follows: λ β δ a γ σ T o 1, K β δ λ 0 T 0, > K 1 The above power-allocaton strategy assgns a zero amplfca-

4 ton gan or equvalently, zero transmt power to the sensors for whch δ s large, because ether the sensor s observaton SNR or ts channel SNR s too low. The assgned nstantaneous transmt power to other sensors s non-zero and based on the value of δ for each sensor. Note that based on the above power-allocaton scheme, there s a unque one-to-one mappng between g and a that could be denoted as a f g. IV. LIMITED FEEDBACK FOR POWER ALLOCATION The optmal power-allocaton scheme proposed n the prevous secton s based on the assumpton that the complete forward channel state nformaton CSI s avalable at local sensors. In other words, Equaton 15 shows that the optmal value of the local amplfcaton gan at sensor s a functon of ts channel SNR γ, whch n tself s a functon of the nstantaneous fadng coeffcent of the channel between sensor and the FC. Therefore, n order to acheve the mnmum L -norm of the vector of local transmt powers, the FC must feed the nstantaneous amplfcaton gan a back to each sensor. 1 Ths requrement s not practcal n most applcatons, especally n large-scale WSNs, snce the feedback nformaton s typcally transmtted through fnte-rate dgtal feedback lnks. In the rest of ths paper, we propose a lmted-feedback strategy to allevate the above-mentoned requrement for nfnte-rate dgtal feedback lnks from the FC to the local sensors. For each channel realzaton, the FC frst fnds the optmal power-allocaton scheme usng the approach proposed n the prevous secton. Note that the FC has access to the perfect backward CSI;.e., the nstantaneous fadng gan of the channel between each sensor and tself. Therefore, t can fnd the exact power-allocaton strategy of the entre network based on 15, gven any channel realzaton. In the next step, the FC sends back the ndex of the quantzed verson of the optmzed power-allocaton vector to all sensors. In the lmted-feedback strategy summarzed above, the FC and local sensors must agree on a codebook of the local amplfcaton gans or equvalently, a codebook of possble power-allocaton schemes. The optmal codebook can be desgned offlne by quantzng the space of the optmzed powerallocaton vectors usng the generalzed Lloyd algorthm [9] wth modfed dstorton metrcs. Let L be the number of feedback bts that the FC uses to quantze the space of the optmal local power-allocaton vectors nto L dsjont regons. Note that L s the total number of feedback bts broadcast by the FC, and not the number of bts fed back to each sensor. A codeword s chosen n each quantzaton regon. The length of each codeword s K, and ts th entry s a real-valued number representng a quantzed verson of the 1 Note that nstead of feedng a back to each sensor, the FC could send back the fadng coeffcent of the channel between each sensor and the FC. However, the knowledge of g alone s not enough for sensor to compute the optmal value of ts local amplfcaton gan a. The sensor must also know whether t needs to transmt or stay slent. There are two ways that the extra data can be fed back to the sensors: Ths nformaton could be encoded n an extra one-bt command nstructng each sensor to transmt or stay slent, or the sensor could lsten for the entre vector of g sent by the FC over a broadcast channel. Sendng back each value of a avods ths extra communcaton. optmal local amplfcaton gan for sensor. The proposed quantzaton scheme could then be thought of as a mappng from the space of channel state nformaton to a dscrete set of L length-k real-valued power-allocaton vectors. Detals of ths quantzaton method are descrbed n the next secton. V. CODEBOOK DESIGN USING LLOYD ALGORITHM Let C [a 1 a a L] T be a L K codebook matrx of the optmal local amplfcaton gans, where [C] l, denotes ts element n row l and column as the optmal gan of sensor n codeword l. Note that each a l, l 1,,..., L s assocated wth a realzaton of the fadng coeffcents of the channels between local sensors and the FC. We apply the generalzed Lloyd algorthm wth modfed dstorton metrcs to solve the problem of vector quantzaton n the space of the optmal local amplfcaton gans. Ths algorthm desgns the optmal codebook C n an teratve process, as explaned n the followng dscussons. In order to mplement the generalzed Lloyd algorthm, a dstorton metrc must be defned for the codebook and for each codeword. Let D B C denote the average dstorton for codebook C defned as D B C E a [ mn l {1,,..., L } ] D W a l,a, 16 where E a [ ] denotes the expectaton operaton wth respect to the optmal vector of local amplfcaton gans and D W a l,a represents the dstance between codeword a l and the optmal power-allocaton vector a, defned as D W a l,a J a l J a, 17 where J s the optmzaton cost of the power-allocaton vector. Let P l and P be the vectors of local transmt powers, when the vector of local amplfcaton gans s a l and a, respectvely. The cost functonja s defned as the L -norm of the correspondng vector of transmt powers P,.e., K 1 K 1 J a P [ a σ o 1+β ]. 18 Let A R K+ be the K-dmensonal vector space of the optmal local amplfcaton gans, whose entres are chosen from the set of real-valued non-negatve numbers. Gven the dstorton functon for the codebook C and that for each one of ts codewords defned n Equatons 16 and 17, respectvely, the two man condtons of the generalzed Lloyd algorthm could be reformulated for our vector-quantzaton problem as follows [9, Chapter 11]: Nearest Neghbor Condton: Ths condton fnds the optmal Vorono cells of the vector space to be quantzed, gven a fxed codebook. Based on ths condton, gven a codebook C, the space A of optmzed power-allocaton vectors s dvded nto L dsjont quantzaton regons or Vorono cells wth the lth regon represented by codeword a l C and defned as A l {a A : D W a l,a D W a k,a, k l}. 19

5 ALGORITHM II: The process of optmal codebook desgn based on the generalzed Lloyd algorthm wth modfed dstorton functons. Requre: K and L. Requre: Fadng model of the channel between local sensors and the FC. Requre: M. M s the number of tranng vectors n space A. Requre: ǫ. ǫ s the dstorton threshold to stop the teratons. 1. Intalzaton. G s A set of M length-k vectors of channel-fadng realzatons based on the gven fadng model of the channels between local sensors and the FC. M L. 3. A s The set of optmal local power-allocaton vectors assocated wth the channel fadng vectors n G s, found by applyng Eq. 15. A s s the set of tranng vectors, and A s A. { } 4. a 0 L l l1 Randomly select L optmal power-allocaton vectors from the set A s as the ntal set of codewords. [ ] T 5. C 0 a 0 1 a0 a0 C 0 s the ntal codebook. L 6. NewCost D B C 0 and j 0. The average dstorton of codebook s found usng Eq EndIntalzaton 8. repeat 9. j j +1 and OldCost NewCost. 10. Gven codebook C j 1, optmally partton the set A s nto L dsjont subsets based on the Nearest-Neghbor Condton usng Eq. 19. Denote the resulted optmal parttons by A j 1 l, l 1,,..., L. 11. for all A j 1 l, l 1,,..., L do 1. a j l Optmal codeword assocated wth partton Aj 1 l 13. end for 14. C j found based on the Centrod Condton usng Eq. 0. [ ] a j T 1 aj aj L 15. NewCost D B C j 16. untl OldCost NewCost ǫ 17. return C OPT C j. C j s the new codebook. Centrod Condton: Ths condton fnds the optmal codebook, gven a specfc parttonng of the vector space to be quantzed. Based on ths condton, gven a specfc parttonng of the space of the optmzed powerallocaton vectors {A 1,A,...,A L}, the optmal codeword assocated wth each Vorono cell A l A s the centrod of that cell wth respect to the dstance functon defned n 17 as a l arg mn a l A l E a Al [D W a l,a], 0 where the expectaton operaton s performed over the set of members of partton A l The optmal codebook s desgned offlne by the FC usng the above two condtons. It can be shown that the average codebook dstorton defned n 16 wll monotoncally decrease through the teratve usage of the Centrod Condton and the Nearest-Neghbor Condton [9, Chapter 11]. Detals of the codebook-desgn process are summarzed n Algorthm II. The optmal codebook s stored n the FC and all sensors. Upon observng a realzaton of the channel fadng vector g, the FC fnds ts assocated optmal power-allocaton vector a OPT, usng 15 to calculate each one of ts elements. It then dentfes the closest codeword n the optmal codebook C to a OPT wth respect to the dstance metrc defned n 17. Fnally, the FC broadcasts the L-bt ndex of that codeword over an error-free dgtal feedback channel to all sensors as l arg mn D W ak,a OPT. 1 k {1,,..., L },a k C Upon recepton of the ndex l, each sensor knows ts quantzed local amplfcaton gan or equvalently, ts powerallocaton weght as [C] l,, where l and are the row and column ndexes of the codebook C, respectvely. VI. NUMERICAL RESULTS In ths secton, numercal results are provded to assess the performance of the optmal power-allocaton scheme proposed n Secton III and to verfy the effectveness of the lmtedfeedback strategy proposed n Secton V n achevng the energy effcency close to that of a WSN wth full CSI feedback. In ths paper, the energy effcency of a powerallocaton scheme s defned as the L -norm of the vector of local transmt powers formulated n 18. In our smulatons, the local observaton gans are randomly chosen from a Gaussan dstrbuton wth unt mean and varance In all smulatons, the average power of h across all sensors s set to be 1.. The observaton and channel nose varances are set to σo 10 dbm and σc 90 dbm, respectvely. The followng fadng model s consdered for the channels between local sensors and the FC: g η 0 d d 0 α f, 1,,...,K, where η 0 30 db s the nomnal fadng gan at the reference dstance set to be d 0 1 meter, d s the dstance between sensor and the FC n meters, α s the path-loss exponent, and f s the ndependent and dentcally dstrbuted..d. Raylegh-fadng random varable wth unt varance. The dstance between sensors and the FC s unformly dstrbuted between 50 and 150 meters. The sze of the tranng set n the optmal codebook-desgn process descrbed n Algorthm II s set to M 5,000, and the codebookdstorton threshold for stoppng the teratve algorthm s assumed to beǫ The results are obtaned by averagng over 10,000 Monte Carlo smulatons. Fgure llustrates the energy effcency of the adaptve power-allocaton scheme proposed n Secton III. The fgure depcts the average L -norm of the vector of local transmt powers versus the maxmum dstorton target for dfferent values of the number of sensors n the network K. The energy effcency for the case of equal power allocaton,.e., the mnmum transmt power requred to acheve the gven target dstorton at the FC, s also shown wth dotted lne as a benchmark. As t can be seen n ths fgure, the energy effcency of the network mproves as the number of sensors ncreases. Ths s due to the fact that when there are fewer sensors n the network, each one of them needs to transmt wth a hgher power n order for the FC to acheve the same estmaton dstorton. Note that n our analyss, there s no constrant on the total transmt power consumed n the entre network. Another observaton from Fg. s that the proposed adaptve power allocaton scheme acheves a hgher energy

6 Average L Norm of Transmt Power Vector Equal Power Allocaton, K10 Optmal Power Allocaton, K10 Equal Power Allocaton, K50 Optmal Power Allocaton, K50 Average L Norm of Transmt Power Vector L L 3 L 4 Full Feedback Maxmum Dstorton Target Fg. : Average energy effcency versus the target estmaton dstorton for the proposed adaptve power-allocaton scheme and the equal power-allocaton strategy. effcency than the equal power-allocaton strategy. As the maxmum dstorton constrant at the FC s relaxed,.e., the value of s ncreased, the gan n the energy effcency decreases slghtly. Fgure 3 llustrates the effect of L as the number of feedback bts from the FC to local sensors on the energy effcency of the proposed power-allocaton scheme. It should be emphaszed that L s the total number of feedback bts broadcast by the FC, and not the number of bts fed back to each sensor. Ths fgure depcts the average L -norm of the vector of local transmt powers versus the maxmum dstorton target for dfferent values of the number of feedback bts L, when there are K 50 sensors n the network. As t can be seen n ths fgure, the energy effcency of the proposed adaptve power allocaton wth lmted feedback s close to that wth full feedback, and gets closer to t as the number of feedback bts s ncreased. VII. CONCLUSIONS In ths paper, an adaptve power-allocaton scheme was proposed that mnmzes the L -norm of the vector of local transmt powers n a WSN, gven a maxmum varance for the BLUE estmator of a random scalar parameter at the FC. Ths approach results n an ncrease n the lfetme of the network at the expense of a potental slght ncrease n the sum total transmt power of all sensors. The next contrbuton of ths paper was to propose a lmted-feedback strategy to elmnate the requrement of nfnte-rate feedback of the nstantaneous forward CSI from the FC to local sensors. Ths scheme desgns an optmal codebook by quantzng the vector space of the optmal local amplfcaton gans usng the generalzed Lloyd algorthm wth modfed dstorton functons. Numercal results showed that the proposed adaptve power-allocaton scheme acheves a hgh energy effcency, and that even wth a lmted number of feedback bts small codebook, ts average energy effcency based on the proposed lmted-feedback strategy s close to that of a WSN wth full CSI feedback Maxmum Dstorton Target Fg. 3: Average energy effcency of the proposed power allocaton scheme versus the target estmaton dstorton for dfferent values of the number of feedback bts L, when there are K 50 sensors n the network. REFERENCES [1] J.-J. Xao, S. Cu, Z.-Q. Luo, and A. J. Goldsmth, Lnear coherent decentralzed estmaton, IEEE Transactons on Sgnal Processng, vol. 56, no., pp , February 008. [] S. Cu, J.-J. Xao, A. J. Goldsmth, Z.-Q. Luo, and H. V. Poor, Estmaton dversty and energy effcency n dstrbuted sensng, IEEE Transactons on Sgnal Processng, vol. 55, no. 9, pp , Sep [3] M. K. Banavar, C. Tepedelenloğlu, and A. Spanas, Estmaton over fadng channels wth lmted feedback usng dstrbuted sensng, IEEE Transactons on Sgnal Processng, vol. 58, no. 1, pp , January 010. [4] M. Fanae, M. C. Valent, N. A. Schmd, and M. M. Alkhweld, Dstrbuted parameter estmaton n wreless sensor networks usng fused local observatons, n Proceedngs of SPIE Wreless Sensng, Localzaton, and Processng VII, vol. 8404, Baltmore, MD, May 01. [5] J.-J. Xao, S. Cu, Z.-Q. Luo, and A. J. Goldsmth, Power schedulng of unversal decentralzed estmaton n sensor networks, IEEE Transactons on Sgnal Processng, vol. 54, no., pp , February 006. [6] A. Rbero and G. B. Gannaks, Bandwdth-constraned dstrbuted estmaton for wreless sensor networks part I: Gaussan case, IEEE Transactons on Sgnal Processng, vol. 54, no. 3, pp , March 006. [7] M. Fanae, M. C. Valent, and N. A. Schmd, Lmtedfeedback-based channel-aware power allocaton for lnear dstrbuted estmaton, n Proceedngs of Aslomar Conference on Sgnals, Systems, and Computers, Pacfc Grove, CA, November 013. [8] S. M. Kay, Fundamentals of Statstcal Sgnal Processng: Estmaton Theory, 1st ed. NJ: Prentce Hall, [9] A. Gersho and R. M. Gray, Vector Quantzaton and Sgnal Compresson. Boston: Kluwer Academc Publshers, 199.

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

Distributed parameter estimation in wireless sensor networks using fused local observations

Distributed parameter estimation in wireless sensor networks using fused local observations Dstrbuted parameter estmaton n wreless sensor networks usng fused local observatons Mohammad Fanae, Matthew C. Valent, Natala A. Schmd, and Marwan M. Alkhweld Lane Department of Computer Scence and Electrcal

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

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

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

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

Asymptotic Quantization: A Method for Determining Zador s Constant

Asymptotic Quantization: A Method for Determining Zador s Constant Asymptotc Quantzaton: A Method for Determnng Zador s Constant Joyce Shh Because of the fnte capacty of modern communcaton systems better methods of encodng data are requred. Quantzaton refers to the methods

More information

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur

Dr. Shalabh Department of Mathematics and Statistics Indian Institute of Technology Kanpur Analyss of Varance and Desgn of Exerments-I MODULE III LECTURE - 2 EXPERIMENTAL DESIGN MODELS Dr. Shalabh Deartment of Mathematcs and Statstcs Indan Insttute of Technology Kanur 2 We consder the models

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

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

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

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

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

Feature Selection: Part 1

Feature Selection: Part 1 CSE 546: Machne Learnng Lecture 5 Feature Selecton: Part 1 Instructor: Sham Kakade 1 Regresson n the hgh dmensonal settng How do we learn when the number of features d s greater than the sample sze n?

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

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

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity

LINEAR REGRESSION ANALYSIS. MODULE IX Lecture Multicollinearity LINEAR REGRESSION ANALYSIS MODULE IX Lecture - 31 Multcollnearty Dr. Shalabh Department of Mathematcs and Statstcs Indan Insttute of Technology Kanpur 6. Rdge regresson The OLSE s the best lnear unbased

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

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

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

Lecture 4: November 17, Part 1 Single Buffer Management

Lecture 4: November 17, Part 1 Single Buffer Management Lecturer: Ad Rosén Algorthms for the anagement of Networs Fall 2003-2004 Lecture 4: November 7, 2003 Scrbe: Guy Grebla Part Sngle Buffer anagement In the prevous lecture we taled about the Combned Input

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

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

The Gaussian classifier. Nuno Vasconcelos ECE Department, UCSD

The Gaussian classifier. Nuno Vasconcelos ECE Department, UCSD he Gaussan classfer Nuno Vasconcelos ECE Department, UCSD Bayesan decson theory recall that we have state of the world X observatons g decson functon L[g,y] loss of predctng y wth g Bayes decson rule s

More information

I + HH H N 0 M T H = UΣV H = [U 1 U 2 ] 0 0 E S. X if X 0 0 if X < 0 (X) + = = M T 1 + N 0. r p + 1

I + HH H N 0 M T H = UΣV H = [U 1 U 2 ] 0 0 E S. X if X 0 0 if X < 0 (X) + = = M T 1 + N 0. r p + 1 Homework 4 Problem Capacty wth CSI only at Recever: C = log det I + E )) s HH H N M T R SS = I) SVD of the Channel Matrx: H = UΣV H = [U 1 U ] [ Σr ] [V 1 V ] H Capacty wth CSI at both transmtter and

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

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

Which Separator? Spring 1

Which Separator? Spring 1 Whch Separator? 6.034 - Sprng 1 Whch Separator? Mamze the margn to closest ponts 6.034 - Sprng Whch Separator? Mamze the margn to closest ponts 6.034 - Sprng 3 Margn of a pont " # y (w $ + b) proportonal

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

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

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

The Expectation-Maximization Algorithm

The Expectation-Maximization Algorithm The Expectaton-Maxmaton Algorthm Charles Elan elan@cs.ucsd.edu November 16, 2007 Ths chapter explans the EM algorthm at multple levels of generalty. Secton 1 gves the standard hgh-level verson of the algorthm.

More information

EEE 241: Linear Systems

EEE 241: Linear Systems EEE : Lnear Systems Summary #: Backpropagaton BACKPROPAGATION The perceptron rule as well as the Wdrow Hoff learnng were desgned to tran sngle layer networks. They suffer from the same dsadvantage: they

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

What would be a reasonable choice of the quantization step Δ?

What would be a reasonable choice of the quantization step Δ? CE 108 HOMEWORK 4 EXERCISE 1. Suppose you are samplng the output of a sensor at 10 KHz and quantze t wth a unform quantzer at 10 ts per sample. Assume that the margnal pdf of the sgnal s Gaussan wth mean

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

Chapter Newton s Method

Chapter Newton s Method Chapter 9. Newton s Method After readng ths chapter, you should be able to:. Understand how Newton s method s dfferent from the Golden Secton Search method. Understand how Newton s method works 3. Solve

More information

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013

ISSN: ISO 9001:2008 Certified International Journal of Engineering and Innovative Technology (IJEIT) Volume 3, Issue 1, July 2013 ISSN: 2277-375 Constructon of Trend Free Run Orders for Orthogonal rrays Usng Codes bstract: Sometmes when the expermental runs are carred out n a tme order sequence, the response can depend on the run

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

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

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

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

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

Chapter 7 Channel Capacity and Coding

Chapter 7 Channel Capacity and Coding Wreless Informaton Transmsson System Lab. Chapter 7 Channel Capacty and Codng Insttute of Communcatons Engneerng atonal Sun Yat-sen Unversty Contents 7. Channel models and channel capacty 7.. Channel models

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

Classification as a Regression Problem

Classification as a Regression Problem Target varable y C C, C,, ; Classfcaton as a Regresson Problem { }, 3 L C K To treat classfcaton as a regresson problem we should transform the target y nto numercal values; The choce of numercal class

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 4: Universal Hash Functions/Streaming Cont d

Lecture 4: Universal Hash Functions/Streaming Cont d CSE 5: Desgn and Analyss of Algorthms I Sprng 06 Lecture 4: Unversal Hash Functons/Streamng Cont d Lecturer: Shayan Oves Gharan Aprl 6th Scrbe: Jacob Schreber Dsclamer: These notes have not been subjected

More information

Two-Way and Multiple-Access Energy Harvesting Systems with Energy Cooperation

Two-Way and Multiple-Access Energy Harvesting Systems with Energy Cooperation Two-Way and Multple-Access Energy Harvestng Systems wth Energy Cooperaton Berk Gurakan, Omur Ozel, Jng Yang 2, and Sennur Ulukus Department of Electrcal and Computer Engneerng, Unversty of Maryland, College

More information

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1

j) = 1 (note sigma notation) ii. Continuous random variable (e.g. Normal distribution) 1. density function: f ( x) 0 and f ( x) dx = 1 Random varables Measure of central tendences and varablty (means and varances) Jont densty functons and ndependence Measures of assocaton (covarance and correlaton) Interestng result Condtonal dstrbutons

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

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

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD

CHALMERS, GÖTEBORGS UNIVERSITET. SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS. COURSE CODES: FFR 135, FIM 720 GU, PhD CHALMERS, GÖTEBORGS UNIVERSITET SOLUTIONS to RE-EXAM for ARTIFICIAL NEURAL NETWORKS COURSE CODES: FFR 35, FIM 72 GU, PhD Tme: Place: Teachers: Allowed materal: Not allowed: January 2, 28, at 8 3 2 3 SB

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

p 1 c 2 + p 2 c 2 + p 3 c p m c 2

p 1 c 2 + p 2 c 2 + p 3 c p m c 2 Where to put a faclty? Gven locatons p 1,..., p m n R n of m houses, want to choose a locaton c n R n for the fre staton. Want c to be as close as possble to all the house. We know how to measure dstance

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

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

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

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

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

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin

LOW BIAS INTEGRATED PATH ESTIMATORS. James M. Calvin Proceedngs of the 007 Wnter Smulaton Conference S G Henderson, B Bller, M-H Hseh, J Shortle, J D Tew, and R R Barton, eds LOW BIAS INTEGRATED PATH ESTIMATORS James M Calvn Department of Computer Scence

More information

Identification of Linear Partial Difference Equations with Constant Coefficients

Identification of Linear Partial Difference Equations with Constant Coefficients J. Basc. Appl. Sc. Res., 3(1)6-66, 213 213, TextRoad Publcaton ISSN 29-434 Journal of Basc and Appled Scentfc Research www.textroad.com Identfcaton of Lnear Partal Dfference Equatons wth Constant Coeffcents

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

Entropy Coding. A complete entropy codec, which is an encoder/decoder. pair, consists of the process of encoding or

Entropy Coding. A complete entropy codec, which is an encoder/decoder. pair, consists of the process of encoding or Sgnal Compresson Sgnal Compresson Entropy Codng Entropy codng s also known as zero-error codng, data compresson or lossless compresson. Entropy codng s wdely used n vrtually all popular nternatonal multmeda

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

Solutions to exam in SF1811 Optimization, Jan 14, 2015

Solutions to exam in SF1811 Optimization, Jan 14, 2015 Solutons to exam n SF8 Optmzaton, Jan 4, 25 3 3 O------O -4 \ / \ / The network: \/ where all lnks go from left to rght. /\ / \ / \ 6 O------O -5 2 4.(a) Let x = ( x 3, x 4, x 23, x 24 ) T, where the varable

More information

Chapter 8 SCALAR QUANTIZATION

Chapter 8 SCALAR QUANTIZATION Outlne Chapter 8 SCALAR QUANTIZATION Yeuan-Kuen Lee [ CU, CSIE ] 8.1 Overvew 8. Introducton 8.4 Unform Quantzer 8.5 Adaptve Quantzaton 8.6 Nonunform Quantzaton 8.7 Entropy-Coded Quantzaton Ch 8 Scalar

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

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl

Suppose that there s a measured wndow of data fff k () ; :::; ff k g of a sze w, measured dscretely wth varable dscretzaton step. It s convenent to pl RECURSIVE SPLINE INTERPOLATION METHOD FOR REAL TIME ENGINE CONTROL APPLICATIONS A. Stotsky Volvo Car Corporaton Engne Desgn and Development Dept. 97542, HA1N, SE- 405 31 Gothenburg Sweden. Emal: astotsky@volvocars.com

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

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

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

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

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

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

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

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS

8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS SECTION 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS 493 8.4 COMPLEX VECTOR SPACES AND INNER PRODUCTS All the vector spaces you have studed thus far n the text are real vector spaces because the scalars

More information

Flexible Quantization

Flexible Quantization wb 06/02/21 1 Flexble Quantzaton Bastaan Klejn KTH School of Electrcal Engneerng Stocholm wb 06/02/21 2 Overvew Motvaton for codng technologes Basc quantzaton and codng Hgh-rate quantzaton theory wb 06/02/21

More information

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M

CIS526: Machine Learning Lecture 3 (Sept 16, 2003) Linear Regression. Preparation help: Xiaoying Huang. x 1 θ 1 output... θ M x M CIS56: achne Learnng Lecture 3 (Sept 6, 003) Preparaton help: Xaoyng Huang Lnear Regresson Lnear regresson can be represented by a functonal form: f(; θ) = θ 0 0 +θ + + θ = θ = 0 ote: 0 s a dummy attrbute

More information

A NEW DISCRETE WAVELET TRANSFORM

A NEW DISCRETE WAVELET TRANSFORM A NEW DISCRETE WAVELET TRANSFORM ALEXANDRU ISAR, DORINA ISAR Keywords: Dscrete wavelet, Best energy concentraton, Low SNR sgnals The Dscrete Wavelet Transform (DWT) has two parameters: the mother of wavelets

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

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

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION

CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING INTRODUCTION CONTRAST ENHANCEMENT FOR MIMIMUM MEAN BRIGHTNESS ERROR FROM HISTOGRAM PARTITIONING N. Phanthuna 1,2, F. Cheevasuvt 2 and S. Chtwong 2 1 Department of Electrcal Engneerng, Faculty of Engneerng Rajamangala

More information

The L(2, 1)-Labeling on -Product of Graphs

The L(2, 1)-Labeling on -Product of Graphs Annals of Pure and Appled Mathematcs Vol 0, No, 05, 9-39 ISSN: 79-087X (P, 79-0888(onlne Publshed on 7 Aprl 05 wwwresearchmathscorg Annals of The L(, -Labelng on -Product of Graphs P Pradhan and Kamesh

More information

Convexity preserving interpolation by splines of arbitrary degree

Convexity preserving interpolation by splines of arbitrary degree Computer Scence Journal of Moldova, vol.18, no.1(52), 2010 Convexty preservng nterpolaton by splnes of arbtrary degree Igor Verlan Abstract In the present paper an algorthm of C 2 nterpolaton of dscrete

More information

Equal-Optimal Power Allocation and Relay Selection Algorithm Based on Symbol Error Probability in Cooperative Communication

Equal-Optimal Power Allocation and Relay Selection Algorithm Based on Symbol Error Probability in Cooperative Communication INTERNATIONAL JOURNAL OF COUNICATIONS Volume 1, 18 Equal-Optmal Power Allocaton and Relay Selecton Algorthm Based on Symbol Error Probablty n Cooperatve Communcaton Xn Song, Syang Xu and ngle Zhang Abstract

More information

PHYS 705: Classical Mechanics. Calculus of Variations II

PHYS 705: Classical Mechanics. Calculus of Variations II 1 PHYS 705: Classcal Mechancs Calculus of Varatons II 2 Calculus of Varatons: Generalzaton (no constrant yet) Suppose now that F depends on several dependent varables : We need to fnd such that has a statonary

More information

Computing MLE Bias Empirically

Computing MLE Bias Empirically Computng MLE Bas Emprcally Kar Wa Lm Australan atonal Unversty January 3, 27 Abstract Ths note studes the bas arses from the MLE estmate of the rate parameter and the mean parameter of an exponental dstrbuton.

More information

Linear Regression Analysis: Terminology and Notation

Linear Regression Analysis: Terminology and Notation ECON 35* -- Secton : Basc Concepts of Regresson Analyss (Page ) Lnear Regresson Analyss: Termnology and Notaton Consder the generc verson of the smple (two-varable) lnear regresson model. It s represented

More information

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016

U.C. Berkeley CS294: Spectral Methods and Expanders Handout 8 Luca Trevisan February 17, 2016 U.C. Berkeley CS94: Spectral Methods and Expanders Handout 8 Luca Trevsan February 7, 06 Lecture 8: Spectral Algorthms Wrap-up In whch we talk about even more generalzatons of Cheeger s nequaltes, and

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

Antenna Combining for the MIMO Downlink Channel

Antenna Combining for the MIMO Downlink Channel Antenna Combnng for the IO Downlnk Channel arxv:0704.308v [cs.it] 0 Apr 2007 Nhar Jndal Department of Electrcal and Computer Engneerng Unversty of nnesota nneapols, N 55455, USA Emal: nhar@umn.edu Abstract

More information

the absence of noise and fading and also show the validity of the proposed algorithm with channel constraints through simulations.

the absence of noise and fading and also show the validity of the proposed algorithm with channel constraints through simulations. 1 Channel Aware Decentralzed Detecton n Wreless Sensor Networks Mlad Kharratzadeh mlad.kharratzadeh@mal.mcgll.ca Department of Electrcal & Computer Engneerng McGll Unversty Montreal, Canada Abstract We

More information

Communication-efficient Distributed Solutions to a System of Linear Equations with Laplacian Sparse Structure

Communication-efficient Distributed Solutions to a System of Linear Equations with Laplacian Sparse Structure Communcaton-effcent Dstrbuted Solutons to a System of Lnear Equatons wth Laplacan Sparse Structure Peng Wang, Yuanq Gao, Nanpeng Yu, We Ren, Janmng Lan, and D Wu Abstract Two communcaton-effcent dstrbuted

More information

4DVAR, according to the name, is a four-dimensional variational method.

4DVAR, according to the name, is a four-dimensional variational method. 4D-Varatonal Data Assmlaton (4D-Var) 4DVAR, accordng to the name, s a four-dmensonal varatonal method. 4D-Var s actually a drect generalzaton of 3D-Var to handle observatons that are dstrbuted n tme. The

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

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

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