MULTIPLE input multiple output (MIMO) systems offer

Size: px
Start display at page:

Download "MULTIPLE input multiple output (MIMO) systems offer"

Transcription

1 IEEE TRANSACTIONS ON COMMUNICATIONS, ACCEPTED FOR PUBLICATION 1 Generalised Sphere Decoding for Spaial Modulaion Abdelhamid Younis, Sinan Sinanović, Marco Di Renzo, Raed Mesleh and Harald Haas arxiv: v3 [cs.it] 3 May 13 Absrac In his paper, Sphere Decoding (SD) algorihms for Spaial Modulaion (SM) are developed o reduce he compuaional complexiy of Maximum Likelihood (ML) deecors. Two SDs specifically designed for SM are proposed and analysed in erms of Bi Error Raio (BER) and compuaional complexiy. Using Mone Carlo simulaions and mahemaical analysis, i is shown ha by carefully choosing he iniial radius he proposed sphere decoder algorihms offer he same BER as ML deecion, wih a significan reducion in he compuaional complexiy. A igh closed form expression for he BER performance of SM SD is derived in he paper, along wih an algorihm for choosing he iniial radius which provides near o opimum performance. Also, i is shown ha none of he proposed SDs are always superior o he ohers, bu he bes SD o use depends on he arge specral efficiency. The compuaional complexiy rade off offered by he proposed soluions is sudied via analysis and simulaion, and is shown o validae our findings. Finally, he performance of SM SDs are compared o Spaial Muliplexing (SMX) applying ML decoder and applying SD. I is shown ha for he same specral efficiency, SM SD offers up o 84% reducion in complexiy compared o SMX SD, wih up o 1 db beer BER performance han SMX ML decoder. Index Terms Muliple inpu muliple oupu (MIMO) sysems, spaial modulaion (SM), spaial muliplexing (SMX), sphere decoding (SD), large scale MIMO. I. INTRODUCTION MULTIPLE inpu muliple oupu (MIMO) sysems offer a significan increase in specral efficiency in comparison o single anenna sysems [1]. An example is Spaial Muliplexing (SMX) [], which ransmis simulaneously over all he ransmi anennas. This mehod achieves a specral efficiency ha increases linearly wih he number of ransmi anennas. However, hese sysems canno cope wih he exponenial increase of wireless daa raffic, and a larger number The associae edior coordinaing he review of his paper and approving i for publicaion was Prof. G. Bauch. Manuscrip received July 9, 1; revised December 4, 1. This work was presened in par a he IEEE GLOBECOM 1, Miami, USA, and IEEE ICC 11, Kyoo, Japan. A. Younis and H. Haas are wih The Universiy of Edinburgh, College of Science and Engineering, Insiue for Digial Communicaions, Join Research Insiue for Signal and Image Processing, King s Buildings, Mayfield Road, Edinburgh, EH9 3JL, UK (e mail: a.younis, h.haas}@ed.ac.uk). S. Sinanović is wih Glasgow Caledonian Universiy, School of Engineering and Buil Environmen, George Moore Building, Cowcaddens Road, Glasgow, G4 BA, UK (e mail: sinan.sinanovic@gcu.ac.uk). M. Di Renzo is wih he Laboraory of Signals and Sysems (LS), French Naional Cener for Scienific Research (CNRS), École Supérieure d Élecricié (SUPÉLEC), Universiy of Paris Sud XI (UPS), 3 rue Jolio Curie, 9119 Gif sur Yvee (Paris), France (e mail: marco.direnzo@lss.supelec.fr). R. Mesleh is wih he Elecrical Engineering Deparmen and SNCS research cener, Universiy of Tabuk, P.O.Box: Tabuk, Saudi Arabia (e mail: rmesleh.sncs@u.edu.sa). Digial Objec Idenifier 1.119/TCOMM /1$5. c 13 IEEE of ransmi anennas (large scale MIMO) should be used [3]. Large scale MIMO sudied in [4] [6], offers higher daa raes and beer bi error rae (BER) performance.however, his comes a he expense of an increase in: 1) Compuaional complexiy: A SMX maximum likelihood (ML) opimum receiver searches across all possible combinaions, and ries o resolve he iner channel inerference (ICI) caused by ransmiing from all anennas simulaneously, on he same frequency. Sphere decoder (SD) was proposed o reduce he complexiy of he SMX ML algorihm while reaining a near opimum performance [7], [8]. The SD reduces he complexiy of he ML decoder by limiing he number of possible combinaions. Only hose combinaions ha lie wihin a sphere cenred a he received signal are considered. However, even hough SMX SD offers a large reducion in complexiy compared o SMX ML, i sill has a high complexiy which increases wih he increase of he number of ransmi anennas. ) Hardware complexiy: In SMX he number of radio frequency (RF) chains is equal o he number of ransmi anennas. From [9], RF chains are circuis ha do no follow Moore s law in progressive improvemen. Therefore, increasing he number of ransmi anennas and consequenly he number of RF chains increases significanly he cos of real sysem implemenaion [1]. 3) Energy consumpion: RF chains conain Power Amplifiers (PAs) which are responsible for 5 8% of he oal power consumpion in he ransmier [11]. Therefore, increasing he number of RF chains resuls in a decrease in he energy efficiency [1]. Thus, SMX may no always be feasible and a more energy efficien and low complexiy soluion should be considered. Spaial Modulaion (SM) is a ransmission echnology proposed for MIMO wireless sysems. I aims o increase he specral efficiency, (m), of single anenna sysems while avoiding ICI [1]. This is aained hrough he adopion of a new modulaion and coding scheme, which foresees: i) he acivaion, a each ime insance, of a single anenna ha ransmis a given daa symbol (consellaion symbol), and ii) he exploiaion of he spaial posiion (index) of he acive anenna as an addiional dimension for daa ransmission (spaial symbol) [13]. Boh he consellaion symbol and he spaial symbol depend on he incoming daa bis. An overall increase by he base wo logarihm of he number of ransmi anennas of he specral efficiency is achieved. This limis he number of ransmi anennas o be a power of wo unless fracional bi encoding SM (FBE SM) [14], or generalised SM (GSM) [15]

2 IEEE TRANSACTIONS ON COMMUNICATIONS, ACCEPTED FOR PUBLICATION are used. Acivaing only one anenna a a ime means ha only one RF chain is needed, which significanly reduces he hardware complexiy of he sysem [16]. Moreover, as only one RF chain is needed, SM offers a reducion in he energy consumpion which scales linearly wih he number of ransmi anennas [1], [17]. Furhermore, he compuaional complexiy of SM ML is equal o he complexiy of single inpu muliple oupu (SIMO) sysems [18], i.e. he complexiy of SM ML depends only on he specral efficiency and he number of receive anennas, and does no depend on he number of ransmi anennas. Recenly he poenial benefis of SM have been validaed no only by simulaions bu also via experimens [19], []. Moreover, in [1] for he firs ime he performance of SM is analysed using real world channel measuremens. Accordingly, SM appears o be a good candidae for large scale MIMO [] [5]. In spie of is low complexiy implemenaion, here is sill poenial for furher reducions, by limiing he number of possible combinaions using he SD principle. However, exising SD algorihms in lieraure do no consider he basic and fundamenal principle of SM, ha only one anenna is acive a any given ime insance. Therefore, wo modified SD algorihms based on he ree search srucure ha are ailored o SM are proposed. The firs SD will be called receiver cenric SD (SM Rx), which was firs presened in [6]. The algorihm in [6] combines he received signal from he muliple receive anennas, as long as he Euclidean disance from he received poin is less han a given radius. This SD based deecor is especially suiable when he number of receive anennas is very large. This echnique reduces he size of he search space relaed o he muliple anennas a he receiver (we denoe his search space as receive search space ). I will be shown laer ha here is no loss in eiher he diversiy order or he coding gain, i.e. he BER is very close o ha of he ML deecor. However, he main limiaion is ha i does no reduce he search space relaed o he number of possible ransmied poins (we denoe his as ransmi search space ). This limiaion prevens he deecor from achieving a significan reducion in compuaional complexiy when high daa raes are required. The second SD, which is called Transmi cenric (SM Tx) was firs presened in [7]. I aims a reducing he ransmi search space by limiing he number of possible spaial and consellaion poins. The SM Tx algorihm avoids an exhausive search by examining only hose poins ha lie inside a sphere wih a given radius. However, SM Tx is limied o he overdeermined MIMO seup (N N r ), where N and N r are he number of ransi and receiver anennas respecively. In [8], [9], i is shown ha SM Tx in [7] can sill be used for he case of (N r 1) N > N r, where SM Tx is referred o as enhanced Tx SD (E Tx SD). Moreover, in [8], [9] a deecor for he case of N > N r referred o as generalised Tx SD (G Tx SD) is proposed. By using he division algorihm he G Tx SD echnique: 1) Divides he se of possible anennas o a number of subses. ) Performs E Tx SD over each subse. 3) Takes he minimum soluion of all he ses. However, in his paper we propose a simple soluion, in which all ha is needed is o se a consan ϕ o for N N r and ϕ = σ n for N > N r, where σ n is he noise variance. In [8], [9], he normalised expeced number of nodes visied by he SM Tx algorihm is used o compare is complexiy wih he complexiy of he SM ML algorihm. This does no ake ino accoun pre compuaions needed by SM Tx. In his paper, when comparing he complexiy of SM Tx wih he complexiy of SM ML and SM Rx, we ake ino accoun he pre compuaions needed by he SM Tx. We show ha because of hose pre compuaions, he SM Tx echnique is no always he bes soluion, where in some cases i is even more complex han SM ML. The performance of boh SDs is closely ied o he choice of he iniial radius. The chosen radius should be large enough for he sphere o conain he soluion. On he one hand, he larger he radius is, he larger he search space, which increases he complexiy. On he oher hand, a small radius may cause he algorihm o fail in finding a poin inside he sphere. In his paper, a careful sudy of he performance of hese wo deecors is provided, along wih an accurae comparison of heir compuaional complexiy. Numerical resuls show ha wih no loss in he BER performance, he proposed soluions provide a subsanial reducion in compuaional complexiy wih respec o he SM ML decoder. We also derive a closed form expression for he BER performance of SM SD and show ha he iniial radius can be chosen such ha SM SD gives an opimum performance. Furhermore, i is shown ha SM Rx is less complex han SM Tx for lower specral efficiencies, while SM Tx is he bes soluion for higher specral efficiencies. Finally, using numerical resuls we show ha SM SD offers a significan reducion and nearly he same performance when compared o SMX wih ML decoder or SD. The remainder of his paper is organised as follows: In secion II, he sysem model along wih he ML opimum deecor is summarised. In secion III, he new SM Rx and SM Tx receivers are described. In secion IV, an accurae analysis of he compuaional complexiy of boh SM Rx and SM Tx is performed. In secion V, he analyical BER performance for SM SDs is derived along wih he iniial radius selecion mehod. Numerical resuls are presened in secion VI, and he paper is concluded in secion VII. A. SM Modulaor II. SYSTEM MODEL In SM, he bi sream emied by a binary source is divided ino blocks conaining m = log (N ) +log (M) bis each, where M is he consellaion size. Then he following mapping rule is used [1]: The firs log (N ) bis are used o selec he anenna which is swiched on for daa ransmission, while he oher ransmi anennas are kep silen. In his paper, he acual ransmi anenna which is acive for ransmission is denoed by l wih l 1,,...,N }. The second log (M) bis are used o choose a symbol in he signal consellaion diagram. Wihou loss of generaliy, Quadraure Ampliude Modulaion (QAM) is considered. In his paper, he acual complex symbol emied by he ransmi anenna l is denoed by s, wih s s 1,s,...,s M }.

3 YOUNIS e al.: GENERALISED SPHERE DECODING FOR SPATIAL MODULATION 3 Accordingly, he N 1 ransmied vecor is: x l,s = [ 1 (l 1),s, 1 (N l )] T, (1) where [ ] T denoes ranspose operaion, and p q is a p q marix wih all zero enries. Noe, a power consrain on he average energy per ransmission of uniy is assumed (i.e.e s = E[x H x] = 1), where E } is he expecaion operaor. From above, he maximum achievable specral efficiency by SM is, However, for SMX, m SM = log (N )+log (M) () m SMX = N log (M) (3) From () and (3), we can see ha he specral efficiency of SM does no increase linearly wih he number of ransmi anennas as SMX does. Therefore, SM needs a larger number of ransmi anennas/ larger consellaion size o arrive a he same specral efficiency as SMX. However, because in SM only one anenna is acive: The compuaional complexiy of SM does no depend on he number of ransmi anennas. Unlike SMX where he compuaional complexiy increases linearly wih he number of ransmi anennas, he compuaional complexiy of SM is he same as he compuaional complexiy of SIMO sysems. The number of RF chains needed by SM is significanly less han he number of RF chains needed by SMX. In fac, only one RF chain is required for SM. For hese reasons we believe ha SM is a srong candidae for large scale MIMO sysems, which srongly moivaes his work. B. Channel Model The modulaed vecor, x l,s, in (1) is ransmied hrough a frequency fla N r N MIMO fading channel wih ransfer funcion H, where N r is he number of receive anennas. In his paper, a Rayleigh fading channel model is assumed. Thus, he enries of H are modelled as complex independen and idenically disribued (i.i.d.) enries according o CN(, 1). Moreover, a perfec channel sae informaion (CSI) a he receiver is assumed, wih no CSI a he ransmier. Thus, he N r 1 received vecor can be wrien as follows: y = Hx l,s +n = h l s +n (4) where n is he N r dimensional Addiive Whie Gaussian Noise (AWGN) wih zero mean and variance σ per dimension a he receiver inpu, and h l is he l h column of H. Noe, he signal-o-noise-raio is SNR = E s /N o = 1/σ n. C. ML Opimum Deecor The ML opimum receiver for MIMO sysems can be wrien as, ˆx (ML) } = argmin y Hx x Q m F where Q m is a m space conaining all possible (N 1) ransmied vecors, F is he Frobenius norm, andˆ denoes he esimaed spaial and consellaion symbols. Noe, in SM only one ransmi anenna is acive a a ime. Therefore, he opimal receiver in (5) can be simplified o, [ˆl(ML),ŝ (ML) ] = arg min l 1,,...N } s s 1,s,...s M} = arg min l 1,,...N } s s 1,s,...s M} y h l s F Nr } } y r h l,r s r=1 where y r andh l,r are he r h enries of y andh l respecively. III. SPHERE DECODERS FOR SM In his secion we inroduce wo SDs ailored for SM, SM Rx and SM Tx. SM Rx aims a reducing he number of summaions over N r required by he ML receiver in (6). SM Tx aims a reducing he number of poins (l,s) he ML receiver searches over. Firs, for ease of derivaion, we inroduce he real valued equivalen of he complex valued model in (4) following [3], where, ȳ = H x l,s + n (5) (6) = h l s + n (7) ȳ = [ Re y T},Im y T}] T [ ] ReH} ImH} H = ImH} ReH} (8) (9) x l,s = [ Re x T l,s},im x T l,s }] T (1) n = [ Re n T},Im n T}] T h l = [ Hl, H ] l+n [ ] Res} s = Ims} (11) (1) (13) where Re } and Im } denoe real and imaginary pars respecively, and H l is he l h column of H. A. SM Rx Deecor The SM Rx is a reduced complexiy and close o opimum BER achieving decoder, which aims a reducing he receive search space. The deecor can formally be wrien as follows: [ˆl(Rx),ŝ (Rx) ] = arg min l 1,,...N } r=1 s s 1,s,...s M} Ñ r(l,s) ȳr h l,r s (14)

4 4 IEEE TRANSACTIONS ON COMMUNICATIONS, ACCEPTED FOR PUBLICATION Ri+1 = z D xl,s F = N (R i R ( )+ zν D (ν,l) Res} D ) (ν,l+n)ims} ν=1 () where h l,r is he r h row of h l, and, } n Ñ r (l,s) = max n y r h l,r s R n 1,,...N r} r=1 (15) The idea behind SM Rx is ha i keeps combining he received signals as long as he Euclidean disance in (14) is less or equal o he radius R. Whenever a poin is found o be inside he sphere, he radius, R, is updaed wih he Euclidean disance of ha poin. The poin wih he minimum Euclidean disance and Ñr(, ) = N r is considered o be he soluion. B. SM Tx Deecor The convenional SD is designed for SMX, where all anennas are acive a each ime insance [7], [31] [33]. However, in SM only one anenna is acive a a ime. Therefore, a modified SD algorihm ailored for SM, named SM Tx, is presened. More specifically, similar o convenional SDs, he SM Tx scheme reduces he number of poins (l,s) for l 1,,...N } and s s 1,s,...s M } o be searched hrough in (6), i.e., he ransmi search space, by compuing he Euclidean disances only for hose poins ha lie inside a sphere wih radius R and are cenred around he received signal. However, unlike convenional SDs, in our scheme he se of poins inside he sphere are much simpler o compue, as here is only a single acive anenna in SM. In his secion, we show how o compue he se of poins. The Cholesky facorisaion of he (N N ) posiive definie marix Ḡ = H T H+ϕĪ N is Ḡ = D D T, where σ ϕ = n N > N r (16) N N r Then he SM Tx can be formally wrien as follow, [ˆl(Tx SD),ŝ (Tx SD) ] } = argmin z D x l,s F (l,s) Θ R (17) where Θ R is he subse of poins (l,s) for l 1,,...N } and s s 1,s,...s M } in he ransmi search space ha lie inside a sphere wih radius R and cenred around he received signal z, z = D ρ and ρ = Ḡ 1 H T ȳ. Unlike SM Rx, SM Tx reduces he compuaional complexiy of he ML receiver by reducing he ransmi search space, which is done by he efficien compuaion of he subse Θ R. Afer some algebraic manipulaions as shown in Appendix A, he subse of poins Θ R lie in he inervals: R i + z l+n Ims} R i + z l+n D (l+n,l+n ) D (l+n,l+n ) R + z l l+n D l,l (18) Res} R + z l l+n D l,l (19) where z a b = z a D (a,b) Ims} () R = R N ν=n +1 z ν l+n (1) Noe, every ime a poin is found inside he sphere, he radius R is updaed as shown in (), wih he Euclidean disance of ha poin. Moreover, (19) needs o be compued only for hose poins ha lie inside he inerval in (18), for he reason ha (19) depends implicily on (18). Because of he unique properies of SM he inervals in (18) and (19) needs o be calculaed only once for each possible ransmi poin, unlike convenional SDs where he inervals have o be calculae N imes for each ransmi poin. Furhermore, we noe ha SM Tx works for boh underdeermined MIMO seup wihn > N r, and overdeermined MIMO seup wih N N r. As opposed o he SM Rx scheme, he SM Tx scheme uses some pre compuaions o esimae he poins ha lie inside he sphere of radius R. These addiional compuaions are carefully aken ino accoun in he analysis of he compuaional complexiy of he SM Tx scheme and is comparison wih he ML opimum deecor in secion IV. IV. COMPUTATIONAL COMPLEXITY OF SM RX AND SM TX In his secion, we analyse he compuaional complexiy of SM ML, SM Rx and SM Tx. The complexiy is compued as he number of real muliplicaion and division operaions needed by each algorihm [34]. A. SM ML The compuaional complexiy of SM ML receiver in (6), yields, C SM ML = 8N r m, (3) as he ML deecor searches hrough he whole ransmi and receive ( search spaces. Noe, evaluaing he Euclidean disance y r h l,r s ) requires 8 real muliplicaions. The compuaional complexiy of SMX ML receiver in (5) is equal o C SMX ML = 4(N +1)N r m. (4) ( Noe, y Hx ) in (5) requires (N +1) complex muliplicaions. From (3) and (4), he complexiy of SM does no depend on he number of ransmi anennas, and i is equal o he complexiy of SIMO sysems. However, he complexiy of SMX increases linearly wih he number of ransmi anennas.

5 YOUNIS e al.: GENERALISED SPHERE DECODING FOR SPATIAL MODULATION 5 Thus, he reducion of SM ML receiver complexiy relaive o he complexiy of he SMX ML decoder for he same specral efficiency is given by, C ML rel = 1 ( 1 N +1 ). (5) From (5), he reducion in complexiy offered by SM increases wih he increase in he number of ransmi anennas. For example for N = 4 SM offers a 6% reducion in complexiy compared o SMX, and as he number of ransmi anennas increases he reducion increases. B. SM Rx The complexiy of he SM Rx receiver is given by: N M C Rx SD = 3 Ñ r (l,s) (6) l=1 s=1 I is easy o show ha C Rx SD lies in he inerval 3 m C Rx SD 6N r m, where he lower bound corresponds o he scenario where Ñ r (l,s) = 1, and he upper bound corresponds o he scenario where Ñ r (l,s) = N r for l 1,,...N } and s s 1,s,...s M }. An ineresing observaion is ha SM Rx offers a reducion in complexiy even for he case of N r = 1, where he complexiy lies in he inerval 3 m C Rx SD 6 m. We noe ha he SM Rx soluion requires no pre compuaions wih respec o he ML opimum deecor. In fac, Ñ r (l,s) for l 1,,...N } and s s 1,s,...s M } in (15) are implicily compued when solving he deecion problem in (14). C. SM Tx The compuaional complexiy of SM Tx can be upper bounded by, C Tx SD C ΘR +3N cardθ R } (7) where card } denoes he cardinaliy of a se, and C ΘR is he complexiy of finding he poins in he subse Θ R, C ΘR = C Pre-Comp +C Inerval (8) where, 1) C Pre-Comp is he number of operaions needed o compue he Cholesky decomposiion. Calculaing he upper riangular marix D using Cholesky decomposiion has he complexiy [34], C CH = 4N 3 /3 (9) I can be easily shown ha he calculaion of Ḡ, ρ and z requires N r N (N +1), N (N +N r +1) and N (N +1) real operaions respecively, where back subsiuion algorihm was used for calculaing ρ [34]. Hence, C Pre-Comp = C CH +N (4N r N +6N r +6N +3) (3) ) C Inerval is he number of operaions needed o compue he inervals in (18),(19), where, C inerval = N +(N +3)N (19) (31) For (18): N real divisions are needed. For (19): (N +3)N (19) real muliplicaions are needed, where (N +3) is he number of real compuaions needed o compue (19), and N (19) is he number of imes (19) is compued, which is calculaed by simulaions. Noe, i) he inerval in (19) depends on he anenna index l and only he imaginary par of he symbol s, ii) some symbols share he same imaginary par. Therefore, (19) is only calculaed for hose poins (l,s) which lie in he inerval in (18) and does no have he same l and Ims} as a previously calculaed poin. V. ERROR PROBABILITY OF SM SDS AND INITIAL RADIUS SELECTION METHOD In his secion, we derive an analyical expresion for he BER performance of SM SD, and we show ha SM SD offers a near opimum performance. The BER for SM SD is esimaed using he union bound [35], which can be expressed as follows, BER SM SD l,s l,s N ( x l,s, x l,s ) E H Pr e,sm SD } m m (3) where N ( x l,s, x l,s ) is he number of bis in error beween x l,s and x l,s, and, ) ) Pr (( lsm SD = Pr, s SM SD (l,s ) (33) e,sm SD is he pairwise) error probabiliy of deciding on he poin ( lsm SD, s SM SD given ha he poin (l,s ) is ransmied. The probabiliy of error Pr e,sm SD can be hough of as wo muually exclusive evens depending on wheher he ransmied poin (l,s ) is inside he sphere. In oher words, he probabiliy of error for SM SD can be separaed in wo pars, as shown in (34) [36]: ( ) Pr ( l SM ML, s SM ML ) (l,s ) : The probabiliy of deciding on he incorrec ransmied symbol and/or used anenna combinaion, given ha he ransmied poin (l,s ) is inside he sphere. Pr((l,s ) / Θ R ): The probabiliy ha he ransmied poin (l,s ) is ouside he se of poins Θ R considered by he SD algorihm. ( ( ) ) Pr Pr ( l SM ML, s SM ML ) (l,s ) +Pr((l,s ) / Θ R ) e,sm SD (34)

6 6 IEEE TRANSACTIONS ON COMMUNICATIONS, ACCEPTED FOR PUBLICATION E H Pr e,sm SD } = [ ( σ ζ s 4σn )] Nr N r 1 r= ( Nr 1+r r )[ 1 ζ ( σ s 4σ n )] r (46) However, he probabiliy of error for he ML decoder is, ( ) Pr Pr ( l, s) (l,s ) (35) e,sm ML Thus, SM SD will have a near opimum performance when, ( ) Pr((l,s ) / Θ R ) << Pr ( l, s) (l,s ) (36) The probabiliy of no having he ransmied poin (l,s ) inside Θ R can be wrien as, ( Nr ) Pr((l,s ) / Θ R ) = Pr ȳr h l,r s > R where, r=1 ( ( ) ) R = Pr κ > σ n / ( ( ) ) R γ N r, σ n = 1 Γ(N r ) N r κ = n r σ n / r=1 (37) (38) is a cenral chi-squared random variable wih N r degree of freedom having a cumulaive disribuion funcion (CDF) equal o [35], F κ (a,b) = γ(b/,a/) (39) Γ(b/) where γ(c, d) is he lower incomplee gamma funcion given by, γ(c,d) = d and Γ(c) is he gamma funcion given by, Γ(c) = c 1 e d (4) c 1 e d (41) The iniial sphere radius considered in SM SD is a funcion of he noise variance as given in [37], R = αn r σ n (4) where α is a consan chosen o saisfy (36). This can be done by seing Pr((l,s ) / Θ R ) = 1 6 and back solving (37). For N r = 1,,4, α = 13.8,8.3,5.3 respecively. Finally, Pr e,sm SD can be formulaed as, ( Pr ȳ hl s = Pr > ȳ hl s ) e,sm SD where, ξ = Re ( hl s h ) } T l s n = Pr ( ξ > h l s h l s ) (43) N ( (,σn hl s h l s )) (44) Thus, hl s Pr = Q h l s (45) e,sm SD σ n where Q(x) = (1/ π) + e / d. x In he case of Rayleigh fading, we can derive he closed form soluion for E H Pr e,sm SD } in (3) by employing he soluion from [38, eq. (6)]. Noe ha he argumen of he Q-funcion in (45) can be represened as he summaion of N r squared Gaussian random variables, wih zero mean and variance equal o 1. This means ha he argumen in he Q funcion can be described by a cenral chi squared disribuion wih N r degrees of freedom. The resul for E H Pr e,sm SD } is as given in (46), where σs = x l,s x l,s F and, ζ(c) = 1 ( ) c 1 (47) 1+c Plugging (46) ino (3) gives a closed form expression for he BER of SM SD. In he nex secion, we show ha (3) gives a igh approximaion of he BER of SM SD, and ha SM SD offers a near opimum performance. VI. RESULTS In he following, Mone Carlo simulaion resuls for a leas 1 6 Rayleigh fading channel realisaions are shown o compare he performance and compuaional complexiy of large scale MIMO, SM ML, SM SD and SMX SD. A. Analyical performance of SM SD Figs. 1- show he BER simulaion resuls for SM ML, SM Rx and SM Tx compared wih he analyical bound Bi Error Raio SM ML, SM ML, SM ML, SM Rx, SM Rx, SM Rx, SM Tx, SM Tx, SM Tx, N =, 3 QAM N = 8, 8 QAM N = 3, BPSK N =, 3 QAM N = 8, 8 QAM N = 3, BPSK N =, 3 QAM N = 8, 8 QAM N = 3, BPSK Analyical, N =, 3 QAM Analyical, N = 8, 8 QAM Analyical, N = 3, BPSK Fig. 1. BER agains SNR. m = 6, and N r = 4.

7 YOUNIS e al.: GENERALISED SPHERE DECODING FOR SPATIAL MODULATION 7 Bi Error Raio SM ML, N =, 18 QAM SM ML, N = 8, 3 QAM SM ML, SM ML, N = 3, 8 QAM N = 18, BPSK SM Rx, N =, 18 QAM SM Rx, N = 8, 3 QAM SM Rx, SM Rx, N = 3, 8 QAM N = 18, BPSK SM Tx, N =, 18 QAM SM Tx, N = 8, 3 QAM SM Tx, SM Tx, N = 3, 8 QAM N = 18, BPSK Analyical, N =, 18 QAM Analyical, N = 8, 3 QAM Analyical, N = 3, 8 QAM Analyical, N = 18, BPSK Fig.. BER agains SNR. m = 8, and N r = 4. derived in secion V, where m = 6,8 and N r = 4. From he figures we can see ha boh SM Tx and SM Rx offer a near opimum performance, where he resuls overlap wih SM ML. Furhermore, Figs. 1- validae our analyical bound as for BER < 1 all graphs closely mach he analyical resuls. Noe, i is will known ha he union bound is loose for low SNR [35]. B. Comparison of he BER performance of SM and SMX Figs. 3 and 4 show a BER comparison beween all possible combinaions of SM and SMX for m = 6 and N r =,4. In Fig. 3, we can observe ha he BER performance depends on he he number of ransmi anennas used and, consequenly, he consellaion size. The smaller he consellaion size, he beer he performance. Anoher observaion ha can be made is ha SM and SMX offer nearly he same performance when using he same consellaion size. In Fig. 4, where he number of receive anennas is increased, we noice ha SM performs Bi Error Raio 1 SM SD, N =, 3 QAM SM SD, N = 8, 8 QAM SM SD, N = 3, BPSK 1 1 SMX SD, N =, 8 QAM SMX SD, N = 6, BPSK Fig. 4. BER agains SNR. m = 6, and N r = 4. 1 Bi Error Raio SM SD, N =, 18 QAM SM SD, N = 8, 3 QAM SM SD, N = 3, 8 QAM SM SD, N = 18, BPSK SMX SD, N =, 16 QAM SMX SD, N = 8, BPSK Fig. 5. BER agains SNR. m = 8, and N r = SM SD, N =, 3 QAM SM SD, N = 8, 8 QAM SM SD, N = 3, BPSK SMX SD, N =, 8 QAM SMX SD, N = 6, BPSK SM SD, N =, 18 QAM SM SD, N = 8, 3 QAM SM SD, N = 3, 8 QAM SM SD, N = 18, BPSK SMX SD, N =, 16 QAM SMX SD, N = 8, BPSK Bi Error Raio 1 Bi Error Raio Fig. 3. BER agains SNR. m = 6, and N r = Fig. 6. BER agains SNR. m = 8, and N r = 4.

8 8 IEEE TRANSACTIONS ON COMMUNICATIONS, ACCEPTED FOR PUBLICATION C rel (%) SM Tx, N=, 3 QAM SM Rx, N= 4, 16 QAM SM Rx, N= 3, BPSK SMX SD, N=, 8 QAM SMX SD, N= 3= QPSK C rel (%) SM Tx, N=, 18 QAM SM Rx, N= 4, 64 QAM SM Rx, N= 18, BPSK SMX SD, N=, 16 QAM SMX SD, N= 4, QPSK Fig. 7. Compuaional complexiy agains SNR. m = 6, and N r =. Fig. 9. Compuaional complexiy agains SNR. m = 8, and N r = SM Tx, N=, 3 QAM SM Rx, N= 4, 16 QAM SM Rx, N= 3, BPSK SMX SD, N=, 8 QAM SMX SD, N= 3, QPSK C rel (%) C rel (%) SM Tx, N=,18 QAM SM Tx, N= 4, 64 QAM SM Rx, N= 8, 3 QAM SM Rx, N= 18, BPSK SMX SD, N=, 16 QAM SMX SD, N= 4, QPSK SMX SD, N= 8, BPSK Fig. 8. Compuaional complexiy agains SNR. m = 6, and N r = 4. Fig. 1. Compuaional complexiy agains SNR. m = 8, and N r = 4. beer han SMX. In paricular, BPSK SM provides a 1 db beer performance han BPSK SMX. Also 8 QAM SM offers a slighly beer performance (.5 db) han 8-QAM SMX. In Figs. 5 and 6, he BER comparisons for m = 8 and N r =,4 are shown. In Fig. 5, SM and SMX offer similar performance for he same consellaion size. However, SM offers a beer performance when he number of receive anennas increases as shown in Fig. 6. In summary, SM offers a similar or beer performance han SMX, where he performance of boh sysems depends on he size of he consellaion diagram and he number of receive anennas. We also noe ha he BER performance of SM can be improved by increasing he number of receive anennas. C. Complexiy Analysis In Figs. 7-1, he compuaional complexiy of SM Rx and SM Tx provided in (6) and (7) respecively is compared wih he compuaional complexiy of SMX SD, where he iniial radius is chosen according o (4). In paricular, he figures show he relaive compuaional complexiy of he SDs wih respec o he SM ML deecor, i.e C rel (%) = 1 (C SD /C ML ). Noe, for SM he SD wih he lowes complexiy is chosen. In Figs. 7 and 8, he relaive compuaional complexiies for m = 6 and N r =,4 are shown. Fig. 7, shows ha for large consellaion sizes he lowes relaive compuaional complexiy is offered by SM Tx N =. The relaive compuaional complexiy ranges beween 4% for low SNR and 16% for high SNR. However, for lower consellaion sizes SM Rx provides he lowes relaive compuaional complexiy, which is beween 56% for low SNR and 6% for high SNR. As SM Rx reduces he receive search space, he reducion in he compuaional complexiy offered by SM Rx does no depend on he number of ransmi anennas. Therefore, only SM Rx wihn = 4,3 are shown, where boh scenarios offer nearly he same relaive compuaional complexiy. Finally, from Fig. 7 we can see ha SMX SD N = and N = 3

9 YOUNIS e al.: GENERALISED SPHERE DECODING FOR SPATIAL MODULATION 9 are less complex han SM ML wih a relaive compuaional complexiy 48% and 79% 8% respecively. However, comparing SM SD o SMX SD N =, for 3 QAM SM SD is 3% less complex han SMX SD, and for BPSK SM SD is % less complex han SMX SD. In Fig. 8, i can be seen ha for large consellaion sizes SM Tx is sill he bes choice wih a relaive complexiy ha ranges beween % for low SNR and 1% for high SNR, which is 15% less han SMX SD N =. For smaller consellaion sizes SM Rx is he bes choice wih relaive complexiy ha ranges beween 55% for low SNR and 14% for high SNR, offering a 3% exra reducion in complexiy when compared o SMX SD N =. Noe, i) SMX SD N = 6 is no shown in he figure, because his scenario has a complexiy higher han he complexiy of SM ML, ii) he complexiy of SMX SD N = 3 increased wih he increase of SNR, for he reason ha, in he underdeermined case ϕ depends on he SNR (5). The relaive complexiy for m = 8 and N r =,4 is shown in Fig. 9 and 1. Since SM Tx reduces he ransmi search space, he reducion in complexiy increased by more han 1% wih he increase in he wordsize and consequenly he consellaion size. In Fig. 9 for high consellaion sizes SM Tx N = is he bes choice wih a relaive complexiy ha reaches 4% for high SNR,. In Fig. 1 for high consellaion sizes SM Tx N = and N = 4 are he bes choice wih a relaive complexiy ha reaches 3% and 1% respecively. On he oher hand, SM Rx reduces he receive search space, herefore, i sill offers nearly he same relaive complexiy. However, he complexiy reduces wih he increase of N r, where SM Rx N r = 4 is ( 1%) less complex han SM Rx N r =. Finally, from boh figures i can be seen ha alhough SM ML is much less complex han SMX ML, SMX SD is less complex han SM ML. For ha reason, SM SD has o be developed, where SM SD is ( %) less complex han SMX SD for N r =, and ( 1%) less complex han SMX SD for N r = 4. Noe, he complexiy of boh SM Tx and SMX SD decreases wih he increase of N r, because for he case of N r < N, he less under-deermined he sysem, he fewer pre compuaions are needed. To summarize, wo SDs for SM are inroduced: SM Tx which reduces he ransmi search space, and SM Rx which reduces he receive search space. Boh deecion algorihms are shown o offer a significan reducion in compuaional complexiy while mainaining a near opimum BER performance. For sysems wih few ransmi anennas, SM Tx is shown o be he beer choice. For sysems wih wih a larger number of receive anennas, SM Rx is shown o be he beer candidae in erms of complexiy reducion. The decision for he mos appropriae SD depends on he paricular deploymen scenario. VII. CONCLUSION In his paper we have inroduced and analysed he performance/complexiy rade off of wo SDs designed specifically for SM. The proposed SDs provide a subsanial reducion in he compuaional complexiy while reaining he same BER as he ML opimum deecor. The closed form analyical performance of SM in i.i.d. Rayleigh fla fading channels has been derived, and analyical and simulaion resuls were shown o closely agree. Furhermore, numerical resuls have highlighed ha no SD is superior o he ohers, and ha he bes soluion o use depends on he MIMO seup, and he SNR a he receiver. In general, SM Rx is he bes choice for lower specral efficiencies, and SM Tx is he bes opion for higher specral efficiencies. Finally, simulaion resuls showed ha SM using SD offers a significan reducion in compuaional complexiy and nearly he same BER performance as SMX using ML decoder or SD. Overall, SM SD offers i) hardware complexiy and power consumpion ha does no depend on he number of ransmi anennas, ii) BER performance ha increases wih he increase of he number of ransmi anennas, and iii) a large reducion in compuaional complexiy compared o SMX. Thus, we believe ha SM SD is an ideal candidae for large scale MIMO sysems. APPENDIX PROOF OF THE INTERVALS (18), (19) Proof: 1) Firs (17) can be hough of as an inequaliy, R ȳ H x l,s } F (48) Then add ϕ x H l,s x l,s o boh sides of (48) o ge (49), where Ḡ = H H H + ϕī N is a (N N ) posiive definie marix, wih a Cholesky facorisaion defined as Ḡ = D H D, where D is a (N N ) upper riangular marix. Now by defining ρ = Ḡ 1 H H ȳ, and adding ρ D H D ρ o boh sides of (49), i can be re wrien as, R ϕ z D x l,s } F N N z i D i,j x l,s (j) i=1 where, z = D ρ and, j=i (5) Rϕ = R +ϕ x T l,s x l,s +ȳ T H ρ ȳ T ȳ (51) σ ϕ = n N > N r (5) N N r For simpliciy, in his paper we assume R ϕ = R. R +ϕ x H l,s x l,s ȳ H x l,s F +ϕ xh l,s x } l,s ȳ H ȳ ȳ H H xl,s + x H l,s H H ȳ+ x H l,sḡ x } l,s (49)

10 1 IEEE TRANSACTIONS ON COMMUNICATIONS, ACCEPTED FOR PUBLICATION ) Second, we noe a necessary condiion ha he poins of he ransmi search space need o saisfy o belong o he subse Θ R is (for all i = 1,,...,N ): N R z i D i,j x l,s (j) (53) j=i which is a condiion similar o convenional SD algorihms [31]. 3) We need o ake ino accoun ha in SM only a single anenna is acive a any ime insance. In he equivalen real valued signal model in (1), his is equivalen o having only wo, ou of N, non zero enries in he signal vecors x l,s and x l,s, respecively. By aking his remark ino accoun, i follows ha: a) if i = N + 1,N +,...,N, hen only he imaginary par of x l,s plays a role in (53), and, hus, only one enry x l,s (l+n ) can be non zero; and b) if i = 1,,...,N, hen boh real and imaginary pars of x l,s play a role in (53), and, hus, only wo enries x l,s (l), x l,s (l+n ) can be non zero. The consideraions in a) and b) lead o he inervals in (18) and (19), respecively, which are direcly obained by solving he inequaliy in (53). REFERENCES [1] E. Telaar, Capaciy of muli-anenna Gaussian channels, European Trans. Telecommun., vol. 1, no. 6, pp , Nov. / Dec [] G. J. Foschini, Layered space-ime archiecure for wireless communicaion in a fading environmen when using muli-elemen anennas, Bell Labs Tech. J., vol. 1, no., pp , [3] J. Hoydis, S. en Brink, and M. Debbah, Massive MIMO: How many anennas do we need? in 11 49h Annual Alleron Conf. Commun., Conrol, Compu. (Alleron), Sep. 11, pp [4] S. Mohammed, A. Zaki, A. Chockalingam, and B. Rajan, High-rae space-ime coded large-mimo sysems: Low-complexiy deecion and channel esimaion, IEEE J. Sel. Topics Signal Process., vol. 3, no. 6, pp , Dec. 9. [5] B. Cerao and E. Vierbo, Hardware implemenaion of a lowcomplexiy deecor for large MIMO, in IEEE In. Symp. Circuis Sys. (ISCAS 9), May 9, pp [6] S. Mohammed, A. Chockalingam, and B. Sundar Rajan, A Lowcomplexiy near-ml performance achieving algorihm for large MIMO deecion, in IEEE In. Symp. Inf. Theory (ISIT 8), July 8, pp [7] E. Vierbo and J. Bouros, A universal laice code decoder for fading channels, IEEE Trans. Inf. Theory, vol. 45, no. 5, pp , July [8] O. Damen, A. Chkeif, and J.-C. Belfiore, Laice code decoder for space-ime codes, IEEE Commun. Le., vol. 4, no. 5, pp , May. [9] A. Molisch and M. Win, MIMO sysems wih anenna selecion, IEEE Microwave Mag., vol. 5, no. 1, pp , Mar. 4. [1] A. Savridis, S. Sinanović, M. D. Renzo., H. Haas, and P. Gran, An energy saving base saion employing spaial modulaion, in IEEE 17h In. Workshop Compu. Aided Modeling Design Commun. Links New. (CAMAD), Sep , pp [11] L. Correia, D. Zeller, O. Blume, D. Ferling, A. Kangas, I. Godor, G. Auer, and L. Van Der Perre, Challenges and enabling echnologies for energy aware mobile radio neworks, IEEE Commun. Mag., vol. 48, no. 11, pp. 66 7, Nov. 1. [1] R. Mesleh, H. Haas, S. Sinanović, C. W. Ahn, and S. Yun, Spaial modulaion, IEEE Trans. Veh. Technol., vol. 57, no. 4, pp. 8 41, July 8. [13] M. Di Renzo, H. Haas, and P. M. Gran, Spaial modulaion for muliple-anenna wireless sysems: A survey, IEEE Commun. Mag,, vol. 49, no. 11, pp , Nov. 11. [14] N. Serafimovski, M. Di Renzo, S. Sinanović, R. Y. Mesleh, and H. Haas, Fracional bi encoded spaial modulaion (FBE SM), IEEE Commun. Le., vol. 14, no. 5, pp , May 1. [15] A. Younis, N. Serafimovski, R. Mesleh, and H. Haas, Generalised Spaial modulaion, in Asilomar Conf. Signals, Sys., Compu., Pacific Grove, CA, USA, Nov. 1. [16] J. Jeganahan, A. Ghrayeb, L. Szczecinski, and A. Ceron, Space shif keying modulaion for MIMO channels, IEEE Trans. Wireless Commun., vol. 8, no. 7, pp , July 9. [17] A. Savridis, S. Sinanović, M. D. Renzo, and H. Haas, A power saving dual-hop archiecure based on hybrid spaial modulaion, in 1 Conf. Record Fory Sixh Asilomar Conf. Signals, Sys. Compu. (ASILOMAR), Nov , pp [18] J. Jeganahan, A. Ghrayeb, and L. Szczecinski, Spaial modulaion: Opimal deecion and performance analysis, IEEE Commun. Le., vol. 1, no. 8, pp , 8. [19] N. Serafimovski, A. Younis, R. Mesleh, P. Chambers, M. D. Renzo, C.-X. Wang, P. M. Gran, M. A. Beach, and H. Haas, Pracical implemenaion of spaial modulaion, in IEEE Trans. Veh. Technol., o appear. [Online]. Available: hp://arxiv.org/abs/ [] N. Serafimovski, A. Younis, and H. Haas, Spaial modulaion, Oc. 1. [Online]. Available: hp://youu.be/cpkibxredho [1] A. Younis, W. Thompson, M. D. Renzo, C.-X. Wang, M. A. Beach, H. Haas, and P. M. Gran, Performance of Spaial Modulaion using Measured Real World Channels, in Proc. of he IEEE Veh. Tech. Conf. (VTC Fall 13), Las Vegas, USA, Sep. 5, 13. [Online]. Available: hp://arxiv.org/abs/ [] M. Di Renzo and H. Haas, Bi Error probabiliy of spaial modulaion (SM) MIMO over generalized fading channels, IEEE Trans. Veh. Technol., vol. 61, no. 3, pp , Mar. 1. [3] M. Di Renzo and H. Haas, On ransmi-diversiy for spaial modulaion MIMO: Impac of spaial-consellaion diagram and shaping filers a he ransmier, in IEEE Trans. Veh. Technol., 13, o appear. Available: IEEE Xplore Early Access. [4] M. D. Renzo, D. D. Leonardis, F. Graziosi, and H. Haas, Space shif keying (SSK) MIMO wih pracical channel esimaes, IEEE Trans. Commun., vol. 6, no. 4, pp , Apr. 1. [5] M. Di Renzo and H. Haas, Bi error probabiliy of space-shif keying MIMO over muliple-access independen fading channels, IEEE Trans. Veh. Technol., vol. 6, no. 8, pp , Oc. 11. [6] A. Younis, R. Mesleh, H. Haas, and P. M. Gran, Reduced complexiy sphere decoder for spaial modulaion deecion receivers, in 1 IEEE Global Telecommun. Conf. (GLOBECOM), Miami, USA, Dec. 1, pp [7] A. Younis, M. Di Renzo, R. Mesleh, and H. Haas, Sphere decoding for spaial modulaion, in Proc. IEEE In. Conf. Commun. (ICC), Kyoo, Japan, 5 9 Jun. 11, pp [8] R. Rajashekar and K. V. S. Hari, Sphere decoding for spaial modulaion sysems wih arbirary N, CoRR, vol. abs/1.5187, 1. [9], Low complexiy maximum likelihood deecion in spaial modulaion sysems, CoRR, vol. abs/16.619, 1. [3] H. V. T. Kailah and B. Hassibi, Space-Time Wireless Sysems: From Array Processing o MIMO Communicaions, C. P. H. Bolcskei, D. Gesber and A. J. van der Veen, Eds. Cambridge Universiy Press, 6. [31] B. Hassibi and H. Vikalo, On he sphere-decoding algorihm I. Expeced complexiy, IEEE Trans. Signal Process., vol. 53, no. 8, pp , Aug. 5. [3] T. Cui and C. Tellambura, An efficien generalized sphere decoder for rank-deficien MIMO sysems, IEEE Commun. Le., vol. 9, no. 5, pp , May 5. [33] P. Wang and T. Le-Ngoc, A low-complexiy generalized sphere decoding approach for underdeermined linear communicaion sysems: Performance and complexiy evaluaion, IEEE Trans. Commun., vol. 57, no. 11, pp , Nov. 9. [34] G. H. Golub and C. F. van Loan, Marix Compuaions. The John Hopkins Universiy Press, [35] J. G. Proakis, Digial Communicaions, 4h ed. McGraw Hill,. [36] J. Jalden, L. Barbero, B. Oersen, and J. Thompson, The error probabiliy of he fixed-complexiy sphere decoder, IEEE Trans. Signal Process., vol. 57, no. 7, pp , July 9. [37] X. Xia, H. Hu, and H. Wang, Reduced iniial searching radius for sphere decoder, in Proc. IEEE 18h In. Symp. Personal, Indoor, Mobile Radio Commun. (PIMRC), Ahens, Greece, Sep. 3 7, 7, pp [38] M.-S. Alouini and A. Goldsmih, A unified approach for calculaing error raes of linearly modulaed signals over generalized fading channels, IEEE Trans. Commun., vol. 47, no. 9, pp , 1999.

11 YOUNIS e al.: GENERALISED SPHERE DECODING FOR SPATIAL MODULATION 11 Abdelhamid Younis received a BSc in Elecrical and Elecronic Engineering wih honours in 7 from he Universiy of Benghazi, Libya and an MSc wih disincion in Signal Processing and Communicaion Engineering in 9 from he Universiy of Edinburgh, UK. He is currenly compleing a Ph.D. in Communicaion Engineering a he Insiue of Digial Communicaions (IDCOM) a he Universiy of Edinburgh where in 1 he was awarded he Overseas Research Suden Award (ORS) in recogniion of his work. His main research ineress lie in he area of wireless communicaion and digial signal processing wih a paricular focus on spaial modulaion, MIMO wireless communicaions, reduced complexiy MIMO design and opical wireless communicaions. Raed Mesleh (S -M 8-SM 13) Dr Mesleh holds a Ph.D. in Elecrical Engineering from Jacobs Universiy in Bremen, Germany and several years of pos-docoral wireless communicaion and opical wireless communicaion research experience in Germany. In Ocober 1, he joined Universiy of Tabuk in Saudi Arabia where he is now an assisan professor and he direcor of research excellence uni. His main research ineress are in spaial modulaion, MIMO cooperaive wireless communicaion echniques and opical wireless communicaion. Dr Mesleh publicaions received more han 8 ciaions since 7. He has published more han 5 publicaions in op-ier journals and conferences, and he holds 7 graned paens. He also serves as on he TPC for academic conferences and is a regular reviewer for mos of IEEE/OSA Communicaion Sociey s journals and IEEE/OSA Phoonics Sociey s journals. Sinan Sinanović (S 98-M 7) is a lecurer a Glasgow Caledonian Universiy. He has obained his Ph.D. in elecrical and compuer engineering from Rice Universiy, Houson, Texas, in 6. In he same year, he joined Jacobs Universiy Bremen in Germany as a pos docoral fellow. In 7, he joined he Universiy of Edinburgh in he UK where he has worked as a research fellow in he Insiue for Digial Communicaions. While working wih Halliburon Energy Services, he has developed acousic elemery receiver which was paened. He has also worked for Texas Insrumens on developmen of ASIC esing. He is a member of he Tau Bea Pi engineering honor sociey and a member of Ea Kappa Nu elecrical engineering honor sociey. He won an honorable menion a he Inernaional Mah Olympiad in Marco Di Renzo (S 5 AM 7 M 9) was born in L Aquila, Ialy, in He received he Laurea (cum laude) and he Ph.D. degrees in Elecrical and Informaion Engineering from he Deparmen of Elecrical and Informaion Engineering, Universiy of L Aquila, Ialy, in April 3 and in January 7, respecively. From Augus o January 8, he was wih he Cener of Excellence for Research DEWS, Universiy of L Aquila, Ialy. From February 8 o April 9, he was a Research Associae wih he Telecommunicaions Technological Cener of Caalonia (CTTC), Barcelona, Spain. From May 9 o December 9, he was an EPSRC Research Fellow wih he Insiue for Digial Communicaions (IDCOM), The Universiy of Edinburgh, Edinburgh, Unied Kingdom (UK). Since January 1, he has been a Tenured Researcher ( Chargé de Recherche Tiulaire ) wih he French Naional Cener for Scienific Research (CNRS), as well as a faculy member of he Laboraory of Signals and Sysems (LS), a join research laboraory of he CNRS, he École Supérieure d Élecricié (SUPÉLEC), and he Universiy of Paris Sud XI, Paris, France. His main research ineress are in he area of wireless communicaions heory. He is a Principal Invesigaor of hree European funded research projecs (Marie Curie ITN GREENET, Marie Curie IAPP WSN4QoL, and Marie Curie ITN CROSSFIRE). Dr. Di Renzo is he recipien of he special menion for he ousanding five year (1997 3) academic career, Universiy of L Aquila, Ialy; he THALES Communicaions fellowship for docoral sudies (3 6), Universiy of L Aquila, Ialy; he Bes Spin Off Company Award (4), Abruzzo Region, Ialy; he Torres Quevedo award for research on ulra wide band sysems and cooperaive localizaion for wireless neworks (8 9), Minisry of Science and Innovaion, Spain; he Dérogaion pour l Encadremen de Thèse (1), Universiy of Paris Sud XI, France; he 1 IEEE CAMAD Bes Paper Award from he IEEE Communicaions Sociey; and he 1 Exemplary Reviewer Award from he IEEE WIRELESS COMMUNICATIONS LETTERS of he IEEE Communicaions Sociey. He currenly serves as an Edior of he IEEE COMMUNICATIONS LETTERS. Professor Harald Haas (S98-A-M3) holds he Chair of Mobile Communicaions in he Insiue for Digial Communicaions (IDCOM) a he Universiy of Edinburgh, and he currenly is he CTO of a universiy spin-ou company purevlc Ld. His main research ineress are in inerference coordinaion in wireless neworks, spaial modulaion and opical wireless communicaion. Prof. Haas holds 3 paens. He has published more han 6 journal papers including a Science Aricle and more han 16 peer-reviewed conference papers. Nine of his papers are invied papers. Prof. Haas has co-auhored a book eniled Nex Generaion Mobile Access Technologies: Implemening TDD wih Cambridge Universiy Press. Since 7 Prof. Haas has been a Regular High Level Visiing Scienis suppored by he Chinese 111-program a Beijing Universiy of Poss and Telecommunicaions (BUPT). He was an invied speaker a he TED Global conference 11, and his work on opical wireless communicaion was lised among he 5 bes invenions in 11 in Time Magazine. He recenly has been awarded a presigious Fellowship of he Engineering and Physical Sciences Research Council (EPSRC) in he UK.

Solutions to the Exam Digital Communications I given on the 11th of June = 111 and g 2. c 2

Solutions to the Exam Digital Communications I given on the 11th of June = 111 and g 2. c 2 Soluions o he Exam Digial Communicaions I given on he 11h of June 2007 Quesion 1 (14p) a) (2p) If X and Y are independen Gaussian variables, hen E [ XY ]=0 always. (Answer wih RUE or FALSE) ANSWER: False.

More information

Block Diagram of a DCS in 411

Block Diagram of a DCS in 411 Informaion source Forma A/D From oher sources Pulse modu. Muliplex Bandpass modu. X M h: channel impulse response m i g i s i Digial inpu Digial oupu iming and synchronizaion Digial baseband/ bandpass

More information

Suboptimal MIMO Detector based on Viterbi Algorithm

Suboptimal MIMO Detector based on Viterbi Algorithm Proceedings of he 7h WSEAS Inernaional Conference on ulimedia Sysems & Signal Processing, Hangzhou, China, April 5-7, 007 9 Subopimal IO Deecor based on Vierbi Algorihm Jin Lee and Sin-Chong Park Sysem

More information

UNIVERSITY OF CALIFORNIA College of Engineering Department of Electrical Engineering and Computer Sciences EECS 121 FINAL EXAM

UNIVERSITY OF CALIFORNIA College of Engineering Department of Electrical Engineering and Computer Sciences EECS 121 FINAL EXAM Name: UNIVERSIY OF CALIFORNIA College of Engineering Deparmen of Elecrical Engineering and Compuer Sciences Professor David se EECS 121 FINAL EXAM 21 May 1997, 5:00-8:00 p.m. Please wrie answers on blank

More information

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB Elecronic Companion EC.1. Proofs of Technical Lemmas and Theorems LEMMA 1. Le C(RB) be he oal cos incurred by he RB policy. Then we have, T L E[C(RB)] 3 E[Z RB ]. (EC.1) Proof of Lemma 1. Using he marginal

More information

EE3723 : Digital Communications

EE3723 : Digital Communications EE373 : Digial Communicaions Week 6-7: Deecion Error Probabiliy Signal Space Orhogonal Signal Space MAJU-Digial Comm.-Week-6-7 Deecion Mached filer reduces he received signal o a single variable zt, afer

More information

st semester. Kei Sakaguchi

st semester. Kei Sakaguchi 0 s semeser MIMO Communicaion Sysems #5: MIMO Channel Capaciy Kei Sakaguchi ee ac May 7, 0 Schedule ( s half Dae Tex Conens # Apr. A-, B- Inroducion # Apr. 9 B-5, B-6 Fundamenals

More information

Lecture 2 October ε-approximation of 2-player zero-sum games

Lecture 2 October ε-approximation of 2-player zero-sum games Opimizaion II Winer 009/10 Lecurer: Khaled Elbassioni Lecure Ocober 19 1 ε-approximaion of -player zero-sum games In his lecure we give a randomized ficiious play algorihm for obaining an approximae soluion

More information

Joint Optimization of Rate Allocation and BLAST Ordering to Minimize Outage Probability

Joint Optimization of Rate Allocation and BLAST Ordering to Minimize Outage Probability Join Opimizaion of Rae Allocaion and BLAST Ordering o Minimize Ouage Probabiliy Arumugam Kannan, Badri Varadarajan and John R. Barry School of Elecrical and Compuer Engineering Georgia Insiue of Technology,

More information

Resource Allocation in Visible Light Communication Networks NOMA vs. OFDMA Transmission Techniques

Resource Allocation in Visible Light Communication Networks NOMA vs. OFDMA Transmission Techniques Resource Allocaion in Visible Ligh Communicaion Neworks NOMA vs. OFDMA Transmission Techniques Eirini Eleni Tsiropoulou, Iakovos Gialagkolidis, Panagiois Vamvakas, and Symeon Papavassiliou Insiue of Communicaions

More information

An introduction to the theory of SDDP algorithm

An introduction to the theory of SDDP algorithm An inroducion o he heory of SDDP algorihm V. Leclère (ENPC) Augus 1, 2014 V. Leclère Inroducion o SDDP Augus 1, 2014 1 / 21 Inroducion Large scale sochasic problem are hard o solve. Two ways of aacking

More information

1 Review of Zero-Sum Games

1 Review of Zero-Sum Games COS 5: heoreical Machine Learning Lecurer: Rob Schapire Lecure #23 Scribe: Eugene Brevdo April 30, 2008 Review of Zero-Sum Games Las ime we inroduced a mahemaical model for wo player zero-sum games. Any

More information

Vehicle Arrival Models : Headway

Vehicle Arrival Models : Headway Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where

More information

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model

Modal identification of structures from roving input data by means of maximum likelihood estimation of the state space model Modal idenificaion of srucures from roving inpu daa by means of maximum likelihood esimaion of he sae space model J. Cara, J. Juan, E. Alarcón Absrac The usual way o perform a forced vibraion es is o fix

More information

Problem Formulation in Communication Systems

Problem Formulation in Communication Systems Problem Formulaion in Communicaion Sysems Sooyong Choi School of Elecrical and Elecronic Engineering Yonsei Universiy Inroducion Problem formulaion in communicaion sysems Simple daa ransmission sysem :

More information

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS NA568 Mobile Roboics: Mehods & Algorihms Today s Topic Quick review on (Linear) Kalman Filer Kalman Filering for Non-Linear Sysems Exended Kalman Filer (EKF)

More information

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing

Application of a Stochastic-Fuzzy Approach to Modeling Optimal Discrete Time Dynamical Systems by Using Large Scale Data Processing Applicaion of a Sochasic-Fuzzy Approach o Modeling Opimal Discree Time Dynamical Sysems by Using Large Scale Daa Processing AA WALASZE-BABISZEWSA Deparmen of Compuer Engineering Opole Universiy of Technology

More information

Robust estimation based on the first- and third-moment restrictions of the power transformation model

Robust estimation based on the first- and third-moment restrictions of the power transformation model h Inernaional Congress on Modelling and Simulaion, Adelaide, Ausralia, 6 December 3 www.mssanz.org.au/modsim3 Robus esimaion based on he firs- and hird-momen resricions of he power ransformaion Nawaa,

More information

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II Roland Siegwar Margaria Chli Paul Furgale Marco Huer Marin Rufli Davide Scaramuzza ETH Maser Course: 151-0854-00L Auonomous Mobile Robos Localizaion II ACT and SEE For all do, (predicion updae / ACT),

More information

CHAPTER 2 Signals And Spectra

CHAPTER 2 Signals And Spectra CHAPER Signals And Specra Properies of Signals and Noise In communicaion sysems he received waveform is usually caegorized ino he desired par conaining he informaion, and he undesired par. he desired par

More information

control properties under both Gaussian and burst noise conditions. In the ~isappointing in comparison with convolutional code systems designed

control properties under both Gaussian and burst noise conditions. In the ~isappointing in comparison with convolutional code systems designed 535 SOFT-DECSON THRESHOLD DECODNG OF CONVOLUTONAL CODES R.M.F. Goodman*, B.Sc., Ph.D. W.H. Ng*, M.S.E.E. Sunnnary Exising majoriy-decision hreshold decoders have so far been limied o his paper a new mehod

More information

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

Mean-square Stability Control for Networked Systems with Stochastic Time Delay JOURNAL OF SIMULAION VOL. 5 NO. May 7 Mean-square Sabiliy Conrol for Newored Sysems wih Sochasic ime Delay YAO Hejun YUAN Fushun School of Mahemaics and Saisics Anyang Normal Universiy Anyang Henan. 455

More information

Západočeská Univerzita v Plzni, Czech Republic and Groupe ESIEE Paris, France

Západočeská Univerzita v Plzni, Czech Republic and Groupe ESIEE Paris, France ADAPTIVE SIGNAL PROCESSING USING MAXIMUM ENTROPY ON THE MEAN METHOD AND MONTE CARLO ANALYSIS Pavla Holejšovsá, Ing. *), Z. Peroua, Ing. **), J.-F. Bercher, Prof. Assis. ***) Západočesá Univerzia v Plzni,

More information

Linear Response Theory: The connection between QFT and experiments

Linear Response Theory: The connection between QFT and experiments Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and

More information

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK

CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 175 CHAPTER 10 VALIDATION OF TEST WITH ARTIFICAL NEURAL NETWORK 10.1 INTRODUCTION Amongs he research work performed, he bes resuls of experimenal work are validaed wih Arificial Neural Nework. From he

More information

Solutions for Assignment 2

Solutions for Assignment 2 Faculy of rs and Science Universiy of Torono CSC 358 - Inroducion o Compuer Neworks, Winer 218 Soluions for ssignmen 2 Quesion 1 (2 Poins): Go-ack n RQ In his quesion, we review how Go-ack n RQ can be

More information

STATE-SPACE MODELLING. A mass balance across the tank gives:

STATE-SPACE MODELLING. A mass balance across the tank gives: B. Lennox and N.F. Thornhill, 9, Sae Space Modelling, IChemE Process Managemen and Conrol Subjec Group Newsleer STE-SPACE MODELLING Inroducion: Over he pas decade or so here has been an ever increasing

More information

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems

Recursive Least-Squares Fixed-Interval Smoother Using Covariance Information based on Innovation Approach in Linear Continuous Stochastic Systems 8 Froniers in Signal Processing, Vol. 1, No. 1, July 217 hps://dx.doi.org/1.2266/fsp.217.112 Recursive Leas-Squares Fixed-Inerval Smooher Using Covariance Informaion based on Innovaion Approach in Linear

More information

INTRODUCTION TO MACHINE LEARNING 3RD EDITION

INTRODUCTION TO MACHINE LEARNING 3RD EDITION ETHEM ALPAYDIN The MIT Press, 2014 Lecure Slides for INTRODUCTION TO MACHINE LEARNING 3RD EDITION alpaydin@boun.edu.r hp://www.cmpe.boun.edu.r/~ehem/i2ml3e CHAPTER 2: SUPERVISED LEARNING Learning a Class

More information

Online Appendix to Solution Methods for Models with Rare Disasters

Online Appendix to Solution Methods for Models with Rare Disasters Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,

More information

Georey E. Hinton. University oftoronto. Technical Report CRG-TR February 22, Abstract

Georey E. Hinton. University oftoronto.   Technical Report CRG-TR February 22, Abstract Parameer Esimaion for Linear Dynamical Sysems Zoubin Ghahramani Georey E. Hinon Deparmen of Compuer Science Universiy oftorono 6 King's College Road Torono, Canada M5S A4 Email: zoubin@cs.orono.edu Technical

More information

ON THE BEAT PHENOMENON IN COUPLED SYSTEMS

ON THE BEAT PHENOMENON IN COUPLED SYSTEMS 8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,

More information

Random Walk with Anti-Correlated Steps

Random Walk with Anti-Correlated Steps Random Walk wih Ani-Correlaed Seps John Noga Dirk Wagner 2 Absrac We conjecure he expeced value of random walks wih ani-correlaed seps o be exacly. We suppor his conjecure wih 2 plausibiliy argumens and

More information

Final Spring 2007

Final Spring 2007 .615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o

More information

Notes on Kalman Filtering

Notes on Kalman Filtering Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren

More information

References are appeared in the last slide. Last update: (1393/08/19)

References are appeared in the last slide. Last update: (1393/08/19) SYSEM IDEIFICAIO Ali Karimpour Associae Professor Ferdowsi Universi of Mashhad References are appeared in he las slide. Las updae: 0..204 393/08/9 Lecure 5 lecure 5 Parameer Esimaion Mehods opics o be

More information

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data

Chapter 2. Models, Censoring, and Likelihood for Failure-Time Data Chaper 2 Models, Censoring, and Likelihood for Failure-Time Daa William Q. Meeker and Luis A. Escobar Iowa Sae Universiy and Louisiana Sae Universiy Copyrigh 1998-2008 W. Q. Meeker and L. A. Escobar. Based

More information

Demodulation of Digitally Modulated Signals

Demodulation of Digitally Modulated Signals Addiional maerial for TSKS1 Digial Communicaion and TSKS2 Telecommunicaion Demodulaion of Digially Modulaed Signals Mikael Olofsson Insiuionen för sysemeknik Linköpings universie, 581 83 Linköping November

More information

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED

0.1 MAXIMUM LIKELIHOOD ESTIMATION EXPLAINED 0.1 MAXIMUM LIKELIHOOD ESTIMATIO EXPLAIED Maximum likelihood esimaion is a bes-fi saisical mehod for he esimaion of he values of he parameers of a sysem, based on a se of observaions of a random variable

More information

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1. Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of

More information

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature

On Measuring Pro-Poor Growth. 1. On Various Ways of Measuring Pro-Poor Growth: A Short Review of the Literature On Measuring Pro-Poor Growh 1. On Various Ways of Measuring Pro-Poor Growh: A Shor eview of he Lieraure During he pas en years or so here have been various suggesions concerning he way one should check

More information

Article from. Predictive Analytics and Futurism. July 2016 Issue 13

Article from. Predictive Analytics and Futurism. July 2016 Issue 13 Aricle from Predicive Analyics and Fuurism July 6 Issue An Inroducion o Incremenal Learning By Qiang Wu and Dave Snell Machine learning provides useful ools for predicive analyics The ypical machine learning

More information

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits

Exponential Weighted Moving Average (EWMA) Chart Under The Assumption of Moderateness And Its 3 Control Limits DOI: 0.545/mjis.07.5009 Exponenial Weighed Moving Average (EWMA) Char Under The Assumpion of Moderaeness And Is 3 Conrol Limis KALPESH S TAILOR Assisan Professor, Deparmen of Saisics, M. K. Bhavnagar Universiy,

More information

The Potential Effectiveness of the Detection of Pulsed Signals in the Non-Uniform Sampling

The Potential Effectiveness of the Detection of Pulsed Signals in the Non-Uniform Sampling The Poenial Effeciveness of he Deecion of Pulsed Signals in he Non-Uniform Sampling Arhur Smirnov, Sanislav Vorobiev and Ajih Abraham 3, 4 Deparmen of Compuer Science, Universiy of Illinois a Chicago,

More information

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t

R t. C t P t. + u t. C t = αp t + βr t + v t. + β + w t Exercise 7 C P = α + β R P + u C = αp + βr + v (a) (b) C R = α P R + β + w (c) Assumpions abou he disurbances u, v, w : Classical assumions on he disurbance of one of he equaions, eg. on (b): E(v v s P,

More information

Impacts of both Tx and Rx IQ Imbalances on OFDM Systems - Analytical Approach

Impacts of both Tx and Rx IQ Imbalances on OFDM Systems - Analytical Approach mpacs of boh Tx and Rx Q mbalances on OFD Sysems - Analyical Approach Hassan Zareian, Vahid Tabaaba Vakili ran Universiy of Science and Technology UST, Tehran, ran Faculy of he slamic Republic of ran Broadcasing

More information

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi

Navneet Saini, Mayank Goyal, Vishal Bansal (2013); Term Project AML310; Indian Institute of Technology Delhi Creep in Viscoelasic Subsances Numerical mehods o calculae he coefficiens of he Prony equaion using creep es daa and Herediary Inegrals Mehod Navnee Saini, Mayank Goyal, Vishal Bansal (23); Term Projec

More information

Combined ML and QR Detection Algorithm for MIMO-OFDM Systems with Perfect ChanneI State Information

Combined ML and QR Detection Algorithm for MIMO-OFDM Systems with Perfect ChanneI State Information Combined ML and QR Deecion Algorihm for MIMO-OFDM Sysems wih Perfec ChanneI Sae Informaion Weizhi You, Lilin Yi, and Weisheng Hu An effecive signal deecion algorihm wih low complexiy is presened for iple-inpu

More information

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017 Two Popular Bayesian Esimaors: Paricle and Kalman Filers McGill COMP 765 Sep 14 h, 2017 1 1 1, dx x Bel x u x P x z P Recall: Bayes Filers,,,,,,, 1 1 1 1 u z u x P u z u x z P Bayes z = observaion u =

More information

GMM - Generalized Method of Moments

GMM - Generalized Method of Moments GMM - Generalized Mehod of Momens Conens GMM esimaion, shor inroducion 2 GMM inuiion: Maching momens 2 3 General overview of GMM esimaion. 3 3. Weighing marix...........................................

More information

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs

A Primal-Dual Type Algorithm with the O(1/t) Convergence Rate for Large Scale Constrained Convex Programs PROC. IEEE CONFERENCE ON DECISION AND CONTROL, 06 A Primal-Dual Type Algorihm wih he O(/) Convergence Rae for Large Scale Consrained Convex Programs Hao Yu and Michael J. Neely Absrac This paper considers

More information

Chapter Floating Point Representation

Chapter Floating Point Representation Chaper 01.05 Floaing Poin Represenaion Afer reading his chaper, you should be able o: 1. conver a base- number o a binary floaing poin represenaion,. conver a binary floaing poin number o is equivalen

More information

Notes for Lecture 17-18

Notes for Lecture 17-18 U.C. Berkeley CS278: Compuaional Complexiy Handou N7-8 Professor Luca Trevisan April 3-8, 2008 Noes for Lecure 7-8 In hese wo lecures we prove he firs half of he PCP Theorem, he Amplificaion Lemma, up

More information

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,

More information

Decentralized Stochastic Control with Partial History Sharing: A Common Information Approach

Decentralized Stochastic Control with Partial History Sharing: A Common Information Approach 1 Decenralized Sochasic Conrol wih Parial Hisory Sharing: A Common Informaion Approach Ashuosh Nayyar, Adiya Mahajan and Demoshenis Tenekezis arxiv:1209.1695v1 [cs.sy] 8 Sep 2012 Absrac A general model

More information

Longest Common Prefixes

Longest Common Prefixes Longes Common Prefixes The sandard ordering for srings is he lexicographical order. I is induced by an order over he alphabe. We will use he same symbols (,

More information

Quantum Computation Based Probability Density Function Estimation

Quantum Computation Based Probability Density Function Estimation June, 5:6 WSPC/NSRUCON LE ijqi nernaional Journal of uanum nformaion c World Scienific Publishing Company uanum Compuaion Based Probabiliy Densiy uncion Esimaion erenc Balázs, Sándor mre Mobile Communicaions

More information

10. State Space Methods

10. State Space Methods . Sae Space Mehods. Inroducion Sae space modelling was briefly inroduced in chaper. Here more coverage is provided of sae space mehods before some of heir uses in conrol sysem design are covered in he

More information

A Shooting Method for A Node Generation Algorithm

A Shooting Method for A Node Generation Algorithm A Shooing Mehod for A Node Generaion Algorihm Hiroaki Nishikawa W.M.Keck Foundaion Laboraory for Compuaional Fluid Dynamics Deparmen of Aerospace Engineering, Universiy of Michigan, Ann Arbor, Michigan

More information

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game

Sliding Mode Extremum Seeking Control for Linear Quadratic Dynamic Game Sliding Mode Exremum Seeking Conrol for Linear Quadraic Dynamic Game Yaodong Pan and Ümi Özgüner ITS Research Group, AIST Tsukuba Eas Namiki --, Tsukuba-shi,Ibaraki-ken 5-856, Japan e-mail: pan.yaodong@ais.go.jp

More information

Adaptive Noise Estimation Based on Non-negative Matrix Factorization

Adaptive Noise Estimation Based on Non-negative Matrix Factorization dvanced cience and Technology Leers Vol.3 (ICC 213), pp.159-163 hp://dx.doi.org/1.14257/asl.213 dapive Noise Esimaion ased on Non-negaive Marix Facorizaion Kwang Myung Jeon and Hong Kook Kim chool of Informaion

More information

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD HAN XIAO 1. Penalized Leas Squares Lasso solves he following opimizaion problem, ˆβ lasso = arg max β R p+1 1 N y i β 0 N x ij β j β j (1.1) for some 0.

More information

Two Coupled Oscillators / Normal Modes

Two Coupled Oscillators / Normal Modes Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own

More information

Comparing Means: t-tests for One Sample & Two Related Samples

Comparing Means: t-tests for One Sample & Two Related Samples Comparing Means: -Tess for One Sample & Two Relaed Samples Using he z-tes: Assumpions -Tess for One Sample & Two Relaed Samples The z-es (of a sample mean agains a populaion mean) is based on he assumpion

More information

Mean Square Projection Error Gradient-based Variable Forgetting Factor FAPI

Mean Square Projection Error Gradient-based Variable Forgetting Factor FAPI 3rd Inernaional Conference on Advances in Elecrical and Elecronics Engineering (ICAEE'4) Feb. -, 4 Singapore Mean Square Projecion Error Gradien-based Variable Forgeing Facor FAPI Young-Kwang Seo, Jong-Woo

More information

Particle Swarm Optimization Combining Diversification and Intensification for Nonlinear Integer Programming Problems

Particle Swarm Optimization Combining Diversification and Intensification for Nonlinear Integer Programming Problems Paricle Swarm Opimizaion Combining Diversificaion and Inensificaion for Nonlinear Ineger Programming Problems Takeshi Masui, Masaoshi Sakawa, Kosuke Kao and Koichi Masumoo Hiroshima Universiy 1-4-1, Kagamiyama,

More information

Introduction to Probability and Statistics Slides 4 Chapter 4

Introduction to Probability and Statistics Slides 4 Chapter 4 Inroducion o Probabiliy and Saisics Slides 4 Chaper 4 Ammar M. Sarhan, asarhan@mahsa.dal.ca Deparmen of Mahemaics and Saisics, Dalhousie Universiy Fall Semeser 8 Dr. Ammar Sarhan Chaper 4 Coninuous Random

More information

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006

2.160 System Identification, Estimation, and Learning. Lecture Notes No. 8. March 6, 2006 2.160 Sysem Idenificaion, Esimaion, and Learning Lecure Noes No. 8 March 6, 2006 4.9 Eended Kalman Filer In many pracical problems, he process dynamics are nonlinear. w Process Dynamics v y u Model (Linearized)

More information

Approximation Algorithms for Unique Games via Orthogonal Separators

Approximation Algorithms for Unique Games via Orthogonal Separators Approximaion Algorihms for Unique Games via Orhogonal Separaors Lecure noes by Konsanin Makarychev. Lecure noes are based on he papers [CMM06a, CMM06b, LM4]. Unique Games In hese lecure noes, we define

More information

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid

More information

Asymptotic Equipartition Property - Seminar 3, part 1

Asymptotic Equipartition Property - Seminar 3, part 1 Asympoic Equipariion Propery - Seminar 3, par 1 Ocober 22, 2013 Problem 1 (Calculaion of ypical se) To clarify he noion of a ypical se A (n) ε and he smalles se of high probabiliy B (n), we will calculae

More information

Some Ramsey results for the n-cube

Some Ramsey results for the n-cube Some Ramsey resuls for he n-cube Ron Graham Universiy of California, San Diego Jozsef Solymosi Universiy of Briish Columbia, Vancouver, Canada Absrac In his noe we esablish a Ramsey-ype resul for cerain

More information

Object tracking: Using HMMs to estimate the geographical location of fish

Object tracking: Using HMMs to estimate the geographical location of fish Objec racking: Using HMMs o esimae he geographical locaion of fish 02433 - Hidden Markov Models Marin Wæver Pedersen, Henrik Madsen Course week 13 MWP, compiled June 8, 2011 Objecive: Locae fish from agging

More information

Regularized Blind Detection for MIMO Communications

Regularized Blind Detection for MIMO Communications Regularized Blind Deecion for MIMO Communicaions Yuejie Chi, Yiyue Wu and Rober Calderbank Deparmen of Elecrical Engineering Princeon Universiy Princeon, NJ 08544, USA Email: ychi, yiyuewu, calderbk@princeon.edu

More information

Proposal of atomic clock in motion: Time in moving clock

Proposal of atomic clock in motion: Time in moving clock Proposal of aomic clock in moion: Time in moving clock Masanori Sao Honda Elecronics Co., d., 0 Oyamazuka, Oiwa-cho, Toyohashi, ichi 441-3193, Japan E-mail: msao@honda-el.co.jp bsrac: The ime in an aomic

More information

3.1 More on model selection

3.1 More on model selection 3. More on Model selecion 3. Comparing models AIC, BIC, Adjused R squared. 3. Over Fiing problem. 3.3 Sample spliing. 3. More on model selecion crieria Ofen afer model fiing you are lef wih a handful of

More information

Lecture 4. Goals: Be able to determine bandwidth of digital signals. Be able to convert a signal from baseband to passband and back IV-1

Lecture 4. Goals: Be able to determine bandwidth of digital signals. Be able to convert a signal from baseband to passband and back IV-1 Lecure 4 Goals: Be able o deermine bandwidh o digial signals Be able o conver a signal rom baseband o passband and back IV-1 Bandwidh o Digial Daa Signals A digial daa signal is modeled as a random process

More information

Chapter 2. First Order Scalar Equations

Chapter 2. First Order Scalar Equations Chaper. Firs Order Scalar Equaions We sar our sudy of differenial equaions in he same way he pioneers in his field did. We show paricular echniques o solve paricular ypes of firs order differenial equaions.

More information

Performance analysis of channel inversion over MIMO channels

Performance analysis of channel inversion over MIMO channels erformance analysis of channel inversion over MIMO channels Damih Senarane and Chinha Tellambura Deparmen of Elecrical and Compuer Engineering, Universiy of Albera, Edmonon, AB, Canada T6G V4 Email: {damih,

More information

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still. Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in

More information

BBP-type formulas, in general bases, for arctangents of real numbers

BBP-type formulas, in general bases, for arctangents of real numbers Noes on Number Theory and Discree Mahemaics Vol. 19, 13, No. 3, 33 54 BBP-ype formulas, in general bases, for arcangens of real numbers Kunle Adegoke 1 and Olawanle Layeni 2 1 Deparmen of Physics, Obafemi

More information

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1.

Timed Circuits. Asynchronous Circuit Design. Timing Relationships. A Simple Example. Timed States. Timing Sequences. ({r 6 },t6 = 1. Timed Circuis Asynchronous Circui Design Chris J. Myers Lecure 7: Timed Circuis Chaper 7 Previous mehods only use limied knowledge of delays. Very robus sysems, bu exremely conservaive. Large funcional

More information

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t... Mah 228- Fri Mar 24 5.6 Marix exponenials and linear sysems: The analogy beween firs order sysems of linear differenial equaions (Chaper 5) and scalar linear differenial equaions (Chaper ) is much sronger

More information

Low-complexity Algorithms for MIMO Multiplexing Systems

Low-complexity Algorithms for MIMO Multiplexing Systems Low-complexiy Algoihms fo MIMO Muliplexing Sysems Ouline Inoducion QRD-M M algoihm Algoihm I: : o educe he numbe of suviving pahs. Algoihm II: : o educe he numbe of candidaes fo each ansmied signal. :

More information

Joint Transmitter-Reciever Optimization for Multiple Input Multiple Output (MIMO) Systems

Joint Transmitter-Reciever Optimization for Multiple Input Multiple Output (MIMO) Systems Join Transmier-Reciever Opimizaion for Muliple Inpu Muliple Oupu (MIMO Sysems eun Chul WANG and wang Bo (Ed LEE School of Elecrical Engineering, Seoul Naional Universiy, OREA Absrac Muliple ransmi (Tx

More information

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation

Most Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,

More information

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol

Pade and Laguerre Approximations Applied. to the Active Queue Management Model. of Internet Protocol Applied Mahemaical Sciences, Vol. 7, 013, no. 16, 663-673 HIKARI Ld, www.m-hikari.com hp://dx.doi.org/10.1988/ams.013.39499 Pade and Laguerre Approximaions Applied o he Acive Queue Managemen Model of Inerne

More information

Evaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance

Evaluation of Mean Time to System Failure of a Repairable 3-out-of-4 System with Online Preventive Maintenance American Journal of Applied Mahemaics and Saisics, 0, Vol., No., 9- Available online a hp://pubs.sciepub.com/ajams/// Science and Educaion Publishing DOI:0.69/ajams--- Evaluaion of Mean Time o Sysem Failure

More information

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems

Single-Pass-Based Heuristic Algorithms for Group Flexible Flow-shop Scheduling Problems Single-Pass-Based Heurisic Algorihms for Group Flexible Flow-shop Scheduling Problems PEI-YING HUANG, TZUNG-PEI HONG 2 and CHENG-YAN KAO, 3 Deparmen of Compuer Science and Informaion Engineering Naional

More information

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H.

ACE 562 Fall Lecture 5: The Simple Linear Regression Model: Sampling Properties of the Least Squares Estimators. by Professor Scott H. ACE 56 Fall 005 Lecure 5: he Simple Linear Regression Model: Sampling Properies of he Leas Squares Esimaors by Professor Sco H. Irwin Required Reading: Griffihs, Hill and Judge. "Inference in he Simple

More information

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits

EECE251. Circuit Analysis I. Set 4: Capacitors, Inductors, and First-Order Linear Circuits EEE25 ircui Analysis I Se 4: apaciors, Inducors, and Firs-Order inear ircuis Shahriar Mirabbasi Deparmen of Elecrical and ompuer Engineering Universiy of Briish olumbia shahriar@ece.ubc.ca Overview Passive

More information

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary

More information

WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM

WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM WATER LEVEL TRACKING WITH CONDENSATION ALGORITHM Shinsuke KOBAYASHI, Shogo MURAMATSU, Hisakazu KIKUCHI, Masahiro IWAHASHI Dep. of Elecrical and Elecronic Eng., Niigaa Universiy, 8050 2-no-cho Igarashi,

More information

Hidden Markov Models. Adapted from. Dr Catherine Sweeney-Reed s slides

Hidden Markov Models. Adapted from. Dr Catherine Sweeney-Reed s slides Hidden Markov Models Adaped from Dr Caherine Sweeney-Reed s slides Summary Inroducion Descripion Cenral in HMM modelling Exensions Demonsraion Specificaion of an HMM Descripion N - number of saes Q = {q

More information

Chapter 4. Truncation Errors

Chapter 4. Truncation Errors Chaper 4. Truncaion Errors and he Taylor Series Truncaion Errors and he Taylor Series Non-elemenary funcions such as rigonomeric, eponenial, and ohers are epressed in an approimae fashion using Taylor

More information

Christos Papadimitriou & Luca Trevisan November 22, 2016

Christos Papadimitriou & Luca Trevisan November 22, 2016 U.C. Bereley CS170: Algorihms Handou LN-11-22 Chrisos Papadimiriou & Luca Trevisan November 22, 2016 Sreaming algorihms In his lecure and he nex one we sudy memory-efficien algorihms ha process a sream

More information

EKF SLAM vs. FastSLAM A Comparison

EKF SLAM vs. FastSLAM A Comparison vs. A Comparison Michael Calonder, Compuer Vision Lab Swiss Federal Insiue of Technology, Lausanne EPFL) michael.calonder@epfl.ch The wo algorihms are described wih a planar robo applicaion in mind. Generalizaion

More information

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course

Overview. COMP14112: Artificial Intelligence Fundamentals. Lecture 0 Very Brief Overview. Structure of this course OMP: Arificial Inelligence Fundamenals Lecure 0 Very Brief Overview Lecurer: Email: Xiao-Jun Zeng x.zeng@mancheser.ac.uk Overview This course will focus mainly on probabilisic mehods in AI We shall presen

More information

Stability and Bifurcation in a Neural Network Model with Two Delays

Stability and Bifurcation in a Neural Network Model with Two Delays Inernaional Mahemaical Forum, Vol. 6, 11, no. 35, 175-1731 Sabiliy and Bifurcaion in a Neural Nework Model wih Two Delays GuangPing Hu and XiaoLing Li School of Mahemaics and Physics, Nanjing Universiy

More information

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri.

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri. UNIVERSIY OF RENO DEPARMEN OF INFORMAION AND COMMUNICAION ECHNOLOGY 385 Povo reno Ialy Via Sommarive 4 hp://www.di.unin.i MEASUREMENS OF RANSIEN PHENOMENA WIH DIGIAL OSCILLOSCOPES Anonio Moschia Fabrizio

More information