Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới. Ứng dụng thông tin vị trí sử dụng các thành phần đa đường trong dự báo kênh cho hệ thống truyền thông vô tuyến thế hệ mới.
MINISTRY OF EDUCATION AND TRAINING HANOI UNIVERSITY OF SCIENCE AND TECHNOLOGY NGUYEN HONG ANH LOCATION-AWARE MULTIPATH-BASED CHANNEL PREDICTION FOR NEXT GENERATION WIRELESS COMMUNICATION SYSTEMS DOCTORAL DISSERTATION OF TELECOMMUNICATIONS ENGINEERING Hanoi−2022 MINISTRY OF EDUCATION AND TRAINING HANOI UNIVERSITY OF SCIENCE AND TECHNOLOGY NGUYEN HONG ANH LOCATION-AWARE MULTIPATH-BASED CHANNEL PREDICTION FOR NEXT GENERATION WIRELESS COMMUNICATION SYSTEMS Major: Telecommunication Engineering Code: 9520208 DOCTORAL DISSERTATION OF TELECOMMUNICATIONS ENGINEERING SUPERVISORS: 1.Assoc Prof Nguyen Van Khang 2.Assoc Prof Klaus Witrisal Hanoi−2022 DECLARATION OF AUTHORSHIP I declare that I have authored this thesis independently, that I have not used other than the declared sources/resources, and that I have explicitly indicated all material which have been quoted either literally or by content from the sources used Hanoi, / / 2022 PhD Student Nguyen Hong Anh SUPERVISORS Assoc.Prof Nguyen Van Khang i ACKNOWLEDGEMENT This dissertation was written during my doctoral course at School of Electronics and Telecommunications (SET), Hanoi University of Science and Technology (HUST) I was also received tremendous supports from the Signal Processing and Speech Com- munication Laboratory (SPSC), Graz University of Technology (TUGraz), Austria I am so grateful for all people who always support and encourage me for completing this study First, I would like to express my sincere gratitude to my advisors for their effective guidance, their patience, continuous support and encouragement, and their immense knowledge I would like to thank all members of SPSC, TUGraz They have been very kind and supportive during my visits to Graz They helped me a lot with their deep understand- ing of the group’s topics and researches I also would like to thank all my colleagues in SET, HUST They have always helped me with the research process and given helpful advice for me to overcome my own difficulties During my Ph.D course, I have received many supports from the Management Board of School of Electronics and Telecommunications Thanks to my employer, HUST for all necessary support and encouragement during my Ph.D journey I am also grateful to Vietnam’s Program 911, for their generous financial support Last but not least, I would like to thank OeAD and SPSC for giving funds for my research visits to Graz Special thanks to my family and relatives for their never-ending support and sacri- fice Hanoi, 2022 Ph.D Student iv CONTENTS DECLARATION OF AUTHORSHIP .i ACKNOWLEDGEMENT ii CONTENTS vi SYMBOLS vi SYMBOLS .ix LIST OF TABLES xiii LIST OF FIGURES xiv CHAPTER 1 INTRODUCTION AND MOTIVATION 1 1.1 .Lite rature review 1 1.1.1 Loc ation-awareness in mmWave beamforming 2 1.1.2 Loc ation-awareness in vehicular communications 3 1.1.3 Loc ation-awareness in adaptive mobile communications, scheduling and routing………… 4 1.1.4 Cha nnel quality metric (CQM) 5 1.2 Challenges and motivations 6 1.3 Purposes and objectives 6 1.4 Research hypotheses 7 1.4.1 Tow ards a site-specific radio propagation modeling .7 1.4.2 Tow ards a large-scale predicting of radio channel statistics 7 1.4.3 Towards a side information-aided single-anchor multipath-based localization 7 1.5 Contributions and outline .8 CHAPTER 2 SIGNAL AND SYSTEM MODELS .9 2.1 .Intr oduction 9 2.2 Syst em model .10 2.2.1 Repr esentation of reflectors using virtual anchors (VAs) 10 2.2.2 Floor plan/environment information for location-aware applications 13 2.3 Hybri d geometric/stochastic signal model 13 2.4 Chan nel quality indicators 15 2.4.1 SMC amplitude 15 2.4.2 Signa l-to-interference-plus-noise ratio (SINR) .16 2.4.3 Chan nel Capacity 17 2.4.4 Positi on error bound (PEB) .19 2.5 Disc ussion 20 2.5.1 Ener gy capture .20 2.5.2 Cont ribution of individual SMCs in the overall channel capacity 21 2.6 Chap ter conclusions .27 CHAPTER 3 GAUSSIAN PROCESS REGRESSION FOR SMC AMPLITUDES 28 3.1 Relat ed Work 28 3.2 SMC propagation model .29 3.3 GP Modeling (GPM) of the SMC Amplitudes .30 3.4 GPR 31 3.4.1 GP Model 31 3.4.2 Predi ction 32 3.4.3 Lear ning 32 3.4.4 Eval uate the quality of prediction 33 3.5 Expe riment and result 34 3.5.1 Expe riment .34 3.5.2 Meas urement pre-processing .34 3.5.3 GPR of SMC Amplitudes 35 3.5.4 GPR of SMC Phases 41 3.6 Chap ter conclusions .44 CHAPTER 4 RADIO ENVIRONMENT MAP FOR SITE-SPECIFIC PROPAGATION MODELING 45 4.1 Relat ed work 45 4.2 Radi o environment map (REM) using Gaussian Process regression (GPR) 47 4.3 SMC amplitudes 47 4.4 SINR 50 4.5 Positi on error bound 52 4.6 Chap ter conclusions .55 CHAPTER 5 APPLICATION OF GPR - ENABLED REMS TO ROBUST POSITIONING 57 5.1 Relat ed work 57 5.2 Probl em formulation .59 5.3 Prop osed algorithm .59 5.4 Resu lt 61 5.5 Chapter conclusions 64 PUBLICATIONS 67 BIBLIOGRAPHY 68 APPENDICES 81 A Description of channel measurement campaigns 81 A.1 Measurement campaign 1 .81 A.2 Measurement campaign 2 .82 B Variance of νk .88 C Predicted Variance 89 ABBREVIATIONS No Abbreviatio Meaning n 1 ACF AutoCorrelation Function 2 ADC Analog-to-Digital Converter 3 AOA Angle-Of-Arrival 4 AOD Angle-Of-Departure 5 AWGN Additive White Gaussian Noise 6 BER Bit Error Rate 7 BF Beam Forming 8 BS Base Station 9 CDF Cumulative Distribution Function 10 CIR Channel Impulse Response 11 CRLB Cramer Rao Lower Bound 12 CQM Channel Quality Metric 13 CSI Channel State Information 14 DM Diffuse Multipath 15 DMC Diffuse Multipath Component 16 EC Energy Capture 17 ECC European Communications Committee 18 EFIM Equivalent Fisher Information Matrix 19 EPB East Plaster Board 20 FCC Federal Communications Commission 21 FIM Fisher Information Matrix 22 GNSS Global Navigation Satellite System 23 GP Gaussian Process 24 GPM Gaussian Process Model 25 GPR Gaussian Process Regression 26 GPS Global Positioning System 27 GSCM Geometry-based Stochastic Channel Model 28 IoT Internet-of-Thing 29 LIDAR Light Detection And Ranging 30 LLHF Log LikeliHood Function 31 M2M Machine-to-Machine 32 MAC Media Access Control 33 MIMO Massive Input Massive Output 34 MINT Multipath-assisted Indoor Navigation and Tracking 35 ML Maximum Likelihood 36 MMSE Minimum Mean Square Error 37 MPC MultiPath Component 38 MRC Maximal Ratio Combining 39 MSLL Mean Square Log Loss 40 NLOS Non-Line-Of-Sight 41 NGW North Glass Window 42 LOS Line-Of-Sight 43 OFDM Orthogonal Frequency Division Multiplexing 44 PA Physical Anchor 45 PAM Pulse Amplitude Modulation 46 PDF Probability Distribution Function 47 PDP Power Delay Profile 48 PHY PHYsical Layer Protocol 49 PEB Position Error Bound 50 QAM Quadrature Amplitude Modulation 51 REM Radio Environment Map 52 RF Radio Frequency 53 RFID Radio Frequency IDentification 54 RRC Root Raised Cosine 55 RSS Received Signal Strength 56 RX Receiver 57 RV Random Variable 58 SALMA 59 SEP Single-Anchor Localization system using Multipath Assistance Symbol Error Probability 60 SIMO Single-Input-Multiple-Output 61 SINR Signal-to-Interference-plus-Noise Ratio 62 SLAM Simultaneous Localization And Mapping 63 SMC Specular Multipath Component 64 SMSE Standard Mean Square Error 65 SNR Signal-to-Noise Ratio Communications 2016.2569087 Letters, 5(4):pp 396–399 doi:10.1109/LWC [111] Han Y., Shen Y., Zhang X.P., Win M.Z., and Meng H (2016) Performance limits and geometric properties of array localization IEEE Transactions on Information Theory, 62(2):pp 1054–1075 doi:10.1109/TIT.2015.2511778 [112] Cadger F., Curran K., Santos J., and Moffett S (2013) A survey of geographical routing in wireless ad-hoc networks IEEE Communications Surveys Tutorials, 15(2):pp 621–653 doi:10.1109/SURV.2012.062612.00109 [113] Muppirisetty L., Svensson T., and Wymeersch H (Feb 2016) Spatial wireless channel prediction under location uncertainty 15(2):pp 1031–1044 ISSN 15361276 doi:10.1109/TWC.2015.2481879 [114] Meissner P., Leitinger E., Lafer M., and Witrisal K (June 2014) Real-time demonstration of multipath-assisted indoor navigation and tracking (mint) In 2014 IEEE International Conference on Communications Workshops (ICC), pp 144–149 ISSN 2164-7038 doi:10.1109/ICCW.2014.6881187 [115] Ulmschneider M (2021) Cooperative Multipath Assisted Positioning Ph.D thesis, Hamburg University of Technology Https://doi.org/10.15480/882.3299 [116] Rath M., Kulmer J., Leitinger E., and Witrisal K (2020) Singleanchor posi- tioning: Multipath processing with non-coherent directional measurements IEEE Access, 8:pp 88115–88132 doi:10.1109/ACCESS.2020.2993197 [117] Großwindhager B., Rath M., Kulmer J., Bakr M.S., Boano C.A., Witrisal K., and Römer K (2018) Salma: Uwb-based singleanchor localization system using multipath assistance Proceedings of the 16th ACM Conference on Embedded Networked Sensor Systems [118] P Stoica A.N (1995) On the concentrated stochastic likelihood function in array signal processing Circuits, Systems and Signal Processing [119] Venus A., Leitinger E., Tertinek S., and Witrisal K (May 2021) A message pass- ing based adaptive pda algorithm for robust radiobased localization and tracking In IEEE Radar Conference Atlanta, GA, USA doi:10.1109/RadarConf2147009 2021.9455311 [120] Wilding T., Leitinger E., and Witrisal K (2021) Multipath-based localization and tracking considering off-body channel effects ArXiv , abs/2110.09932 [121] Meissner P., Leitinger E., Fröhle M., and Witrisal K (2013) Accurate and robust indoor localization systems using ultra- wideband signals In European Navigation Conference European Navigation Conference, ENC ; Conference date: 23-04- 2013 Through 25-04-2013 [122] Santos T., Karedal J., Almers P., Tufvesson F., and Molisch A.F (2010) Model- ing the ultra-wideband outdoor channel: Measurements and parameter extraction method IEEE Transactions on Wireless Communications, 9(1):pp 282–290 doi: 10.1109/TWC.2010.01.090391 [123] Borish J (Mar 1984) Extension of the image model arbitrary polyhedra JASA, 75(6):pp 1827–1836 to [124] Sachs J., Herrmann R., Kmec M., Helbig M., and Schilling K (2007) Recent Advances and Applications of M-Sequence based Ultra-Wideband Sensors, pp 50–55 Singapore doi:10.1109/ICUWB.2007.4380914 [125] Molisch A.F (December 2010) Wireless communications Wiley-IEEE press, 2 edition [126] Cepeda R., Parker S.C.J., and Beach M (2007) The measurement of frequency dependent path loss in residential los environments using time domain uwb chan- nel sounding In 2007 IEEE International Conference on Ultra-Wideband , pp 328–333 doi:10.1109/ICUWB.2007.4380964 [127] Meissner P., Leitinger E., Lafer M., and Witrisal K (2013) MeasureMINT UWB database APPENDICES A Description of channel measurement campaigns A.1 Measurement campaign 1 Frequency Domain Measurements - Vector Network Analyzer Frequency-domain measurements have been obtained with a Rhode & Schwarz ZVA- 24 VNA The frequency range has been chosen as the full FCC bandwidth from 3.1 to 10.6 GHz (corresponding to a wavelength range of 9.67 cm to 2.83 cm), resulting in a delay resolution of 0.1333 ns and a spatial resolution of 4 cm At the l-th trajectory position, a sampled version Hl[k] of the CTF Hl(f ) with a frequency spacing of ∆f is measured The VNA has been calibrated up to (but not including) the antennas with a through-open-short-match (TOSM) calibration The FCC bandwidth has been measured for different discrete frequencies with a frequency resolution of 1.5 MHz The transmit power has been set to 15 dBm Measurement Post Processing For the VNA measurements, the major system influences on the measured CTF H(f ) have already been removed by the previously mentioned TOSM calibration This includes cables and connectors, but not the antennas, which are considered as part of the transmission channel The necessary post-processing tasks reduce to a filtering of the signal to select a desired frequency band out of the FCC range and to downconvert the signal transformed to time domain to obtain a baseband signal The filtering is done with a baseband pulse s(t) that covers the desired bandwidth The CTF is measured at Nf discrete frequencies fk = k∆f + fmin, k = 0, , Nf − 1, where fmin is the lowest measured frequency This sampled CTF H[k] corresponds to a Fourier series representation of the timedomain CIR h(τ ) [121], which is periodic with a period of τmax With f0 and fc denoting the lower band edge and the center frequency of the extracted band, respectively, and using an IFFT with size NFFT = ⌈(∆f ∆τ )−1⌉, where ∆τ is the desired delay resolution, the time domain equivalent baseband signal is obtained as r(t) = IFFTNFFT {H[k]S[k]}e−j2π(fc −f0 )t (A.1) Here, S[k] is the discrete frequency domain representation of the pulse s(t) in the desired frequency range This procedure is similar to [122] Measurement scenario We consider the simple scenario shown in Fig A.1, where one physical T anchor is present at position a1 = [4.2, 4] and the mobile agent is T placed at position p = [3.4, 1.4] UWB grid measurements are available for 484 (22x22) grid points pℓ with a spacing of 5 × 5 cm, surrounding p The measurements were performed us- ing a Rhode & Schwarz ZVA24 vector network analyzer with frequency range from 3.1 − 10.6 GHz, thereby covering the full FCC-regulated band for UWB Agent and anchor were equipped with dipole-like antennas made of Euro Cent coins mounted at a height of 1.5 m These antennas have an approximately uniform radiation pattern in the azimuth plane and zeros in the directions of floor and ceiling Within the total measured band, we selected the actual signal band using filtering with a raised cosine pulse s(t) with a roll-off-factor of 0.6 Varying values have been selected for the two- sided bandwidth, namely 100 MHz, 500 MHz, and 2 GHz, each at a carrier frequency of fc = 7 GHz A.2 Measurement campaign 2 Measurement scenario We consider indoor environments where fixed anchors communicate with a mobile agent by means of radio signals Figure A.3 illustrates the considered scenario, where two physical anchor positions at positions a(1) and a(2) are shown as blue crosses, the 1 1 agent positions p along a segmented trajectory, and some exemplary virtual anchors (VAs) are shown The VA positions are mirror images of the physical anchor positions that are induced by reflections at flat surfaces—typically walls—and thus depend on the surrounding environment (floor plan) [123] The position of the k-th VA of the jth physical anchor at position a(j) is denoted as a(j) Note that in this work we consider 1 k horizontal propagation only For brevity, we neglect the anchor index j from now on Also, we denote L as the set of measurement points Figure A.3 shows the laboratory room at Graz University of Technology that was used for the experimental validation The room consists of two plaster board walls and two reinforced concrete walls (shown as black outer lines), three glass windows at the north wall (shown as thick gray lines), one white board and one metal door at the south wall (indicated by A∗ and C∗, respectively) We introduce the following labels to refer to the involved reflection surfaces: EPB East plaster board SW South wall WW West wall NGW North glass wall 6 uwin 4 A2 2 p 0 y m −2 rwin lwall lwin −4 −6 A278 −8 −5 0 5 10 15 20 xm (a) Overview of floorplan 5 4 3 2 y m 1 0 −1 −2 2 4 6 xm 8 10 12 (b) Close-up of floorplan Figure A.1: Scenario floor-plan: a physical anchor is located at position a1 and an examplary VA is at position a2 The gray grid with positions pℓ indicates the × measurement grid with 5 5 cm spacing; the red dot indicates its center position p, the actual mobile agent position used in the illustration Blue lines depict specular reflections at wall segments Figure A.2: Photo of corridor scenario To conduct the channel measurements, an Ilmsens Ultra-Wide band Msequence de- vice [124] was used, c.f Figure A.4a The measurement principle is correlative channel sounding [125], i.e a binary code sequence with suitable autocorrelation properties is transmitted over the channel At the receiver, the channel impulse response is re- covered using a correlation with the known code sequence The channel sounder has one transmitter port and two receiver ports A 12-bit M-sequence has been employed, corresponding to a sequence length of 4095 samples This allows for an unambiguous delay window of 589.2 ns at a clock rate of 6.95 GHz The M-sequence is modulated onto a 6.95 GHz carrier, yielding a probing signal that covers a frequency band between approx 3.5 and 10.5 GHz Each of the ports was connected to a dipole coin antenna as shown in Figure A.4b According to [81], the coin antenna has a very wide bandwidth ranging from 3 to 9 GHz It also has with a nearly isotropic radiation pattern in the horizontal plane We used the two receiver ports as anchors and placed their antennas at fixed positions a(1) and a(2) The transmitter port is connected to another antenna that was 1 1 moved along a trajectory with 595 points p, as shown in the Figure A.3, to obtain the same number of channel measurements All antennas were mounted on tripods at the same height, therefore only the copolarized, azimuth radiation pattern of the antenna has an impact on the data The raw measurements at the receiver ports were filtered with an RRC pulse with center frequency 6.95 GHz, roll-off factor 0.5 and bandwidth A2 cpill p pl 11 segment 1 10 segment 2 segment 3 segment 4 9 segment 5 segment 6 8 segment 7 phys anchor VA a(1) 5 FE D CB A 7 ydir ect ion in me ter 6 D† a(1) 4 5 4 a(1) 2 a(1) p 1 C† 3 B† φk 2 1 A† a(2) 4 C∗ 0 a(2) 1 B∗ A∗ a(2) −1 −2 −13 a(2) 2 3 −12 −11 −10 −9 −8 −7 −6 −5 x-direction in meter −4 −3 −2 −1 0 1 Figure A.3: Floor plan of the evaluation scenario Bold black lines denote walls, thick gray lines represent glass windows, other lines illustrate other materials Two blue crosses represent the physical anchors; orange circles denote virtual anchors (VAs) which were considered in the experimental vali- dation An agent moves along a trajectory segmented into seven parts indicated with distinct colors Capital letters (with or without mark or ) refer to ∗ sub-segments of different materials along each wall (a) Ilmsens UWB M-sequence sounder and Ilmsens power supply (b) Europe coin antenna Figure A.4: A photo of the Ilmsens channel sounder, and a photo of the coin antenna used for transmit and receive 1/Tp = 2 GHz to obtain the received signals corresponding to the model in (2.3) The power spectral density of AWGN N0 is known and considered in the training and evaluation process Time domain measurement - M-Sequence Radar Time-domain measurements have been obtained with an Ilmsens UltraWide Band M-Sequence device [124] The measurement principle is correlative channel sounding [125] A binary code sequence with suitable autocorrelation properties (a large peak- to-off-peak-ratio) is transmitted over the channel At the receiver, the channel impulse response is recovered using a correlation with the known code sequence This Msequence radar has one transmitter and two receiver ports Hence, the mobile unit that has been moved along the measurement trajectories was the transmitter, and the two receiver ports have been used as anchors The transmit power of the M-sequence device in FCC mode is 18dBm The employed 12-bit M-sequence has a length of 4095 samples At the clock rate of 6.95GHz, this allows for a maximum delay of τmax = 589.2 ns Measurement post-processing Figure A.5 shows a block diagramm of the measurement setup using the M-Sequence radar As in the VNA measurements, the measurement system should be calibrated Figure A.5: Calibration setup for time domain measurements up to (but not including) the antennas Hence, the influence of the device internal transfer functions and the measurement cables and connectors, combined in the transfer function Hsys,i(f ) for the i-th RX channel, as well as the crosstalk between TX channel and i-th RX channel, Hcross,i(f ), have to be compensated For the further description, we will drop the channel index To achieve this, two types of measurements are necessary First, to determine the crosstalk, the TX antenna is unmounted and the TX port is terminated with a 50Ω match and the crosstalk signals are measured Second, also the RX antennas are unmounted and TX and RX cables are connected In this way, Hmeas(f ) = Hsys(f ) + Hcross(f ) are measured Using the measurement configuration with all the antennas as depicted in Figure A.5 yields Hmeas(f ) = H(f )Hsys(f )+Hcross(f ) Hence, a calibrated version of the radio channel transfer function is obtained as Hmeas(f ) − Hcross(f ) H(f ) = (A.2) Hsys(f ) − Hcross(f ) To avoid excessive noise gain, we use a thresholding on the time-domain representation of the denominator in (A.2) and set samples below the threshold to zero The time domain signal is obtained by an inverse Fourier transformation Finally, the time- domain signal within the ... reach their intended destination In its most basic form, given a destination d, a node i with neighbors Ni will choose to forward data to a neighbor closest to the destination 1.1.4 Channel quality... their limited coverage Last but not least, a well-known technique in network routing is geographic routing (georouting), which takes advantage of geographic information of nodes (actual geo- graphic... has to be searched exhaustively to align the beam pointing angles in current beam training solutions In [17], a database was built by collecting received power at every locations of the receiver