Advanced Topics in Digital Communications Spezielle Methoden der digitalen Datenübertragung

Similar documents
Lecture 7: Support Vector Machine

Contents. I Basics 1. Copyright by SIAM. Unauthorized reproduction of this article is prohibited.

Wireless Networks Research Seminar April 22nd 2013

CPSC 340: Machine Learning and Data Mining. Principal Component Analysis Fall 2016

Multi-Hop Virtual MIMO Communication using STBC and Relay Selection

CS 6210 Fall 2016 Bei Wang. Review Lecture What have we learnt in Scientific Computing?

LECTURE 13: SOLUTION METHODS FOR CONSTRAINED OPTIMIZATION. 1. Primal approach 2. Penalty and barrier methods 3. Dual approach 4. Primal-dual approach

Two-view geometry Computer Vision Spring 2018, Lecture 10

CS545 Contents IX. Inverse Kinematics. Reading Assignment for Next Class. Analytical Methods Iterative (Differential) Methods

CPSC 340: Machine Learning and Data Mining. Kernel Trick Fall 2017

calibrated coordinates Linear transformation pixel coordinates

Principles of Wireless Sensor Networks. Fast-Lipschitz Optimization

Multiple View Geometry in Computer Vision

Advanced Techniques for Mobile Robotics Graph-based SLAM using Least Squares. Wolfram Burgard, Cyrill Stachniss, Kai Arras, Maren Bennewitz

Divide and Conquer Kernel Ridge Regression

CLASSIFICATION AND CHANGE DETECTION

ME 555: Distributed Optimization

Generalized trace ratio optimization and applications

Tracking system. Danica Kragic. Object Recognition & Model Based Tracking

Realization of Hardware Architectures for Householder Transformation based QR Decomposition using Xilinx System Generator Block Sets

IMAGE ANALYSIS, CLASSIFICATION, and CHANGE DETECTION in REMOTE SENSING

COMPUTATIONAL INTELLIGENCE (CS) (INTRODUCTION TO MACHINE LEARNING) SS16. Lecture 2: Linear Regression Gradient Descent Non-linear basis functions

A Two-phase Distributed Training Algorithm for Linear SVM in WSN

Multi-area Nonlinear State Estimation using Distributed Semidefinite Programming

Compressive Sensing for Multimedia. Communications in Wireless Sensor Networks

MATH 423 Linear Algebra II Lecture 17: Reduced row echelon form (continued). Determinant of a matrix.

Introduction to Deep Learning

Convex Optimization. Lijun Zhang Modification of

Image Analysis, Classification and Change Detection in Remote Sensing

Parallel Implementations of Gaussian Elimination

Lecture 8 Fitting and Matching

Algebraic Iterative Methods for Computed Tomography

Image Processing. Image Features

Deep Learning for Computer Vision

CS 565 Computer Vision. Nazar Khan PUCIT Lectures 15 and 16: Optic Flow

APP-PHY Interactions in Wireless Networks

Advanced phase retrieval: maximum likelihood technique with sparse regularization of phase and amplitude

Model Reduction for Variable-Fidelity Optimization Frameworks

ICRA 2016 Tutorial on SLAM. Graph-Based SLAM and Sparsity. Cyrill Stachniss

Robust Principal Component Analysis (RPCA)

Feature selection. Term 2011/2012 LSI - FIB. Javier Béjar cbea (LSI - FIB) Feature selection Term 2011/ / 22

A Course in Machine Learning

The Pre-Image Problem in Kernel Methods

CPSC 340: Machine Learning and Data Mining. Principal Component Analysis Fall 2017

Practical Least-Squares for Computer Graphics

Biometrics Technology: Image Processing & Pattern Recognition (by Dr. Dickson Tong)

Robust Face Recognition via Sparse Representation Authors: John Wright, Allen Y. Yang, Arvind Ganesh, S. Shankar Sastry, and Yi Ma

Probabilistic Graphical Models

Recent Developments in Model-based Derivative-free Optimization

COMPUTATIONAL INTELLIGENCE (INTRODUCTION TO MACHINE LEARNING) SS18. Lecture 2: Linear Regression Gradient Descent Non-linear basis functions

Multiple-View Object Recognition in Band-Limited Distributed Camera Networks

Exercise Set Decide whether each matrix below is an elementary matrix. (a) (b) (c) (d) Answer:

Lecture 9 Fitting and Matching

Bilevel Sparse Coding

New Approaches for EEG Source Localization and Dipole Moment Estimation. Shun Chi Wu, Yuchen Yao, A. Lee Swindlehurst University of California Irvine

Concept of Curve Fitting Difference with Interpolation

Sparse sampling in MRI: From basic theory to clinical application. R. Marc Lebel, PhD Department of Electrical Engineering Department of Radiology

Stability Analysis of the Muscl Method on General Unstructured Grids for Applications to Compressible Fluid Flow

Least-Squares Minimization Under Constraints EPFL Technical Report #

Robot Mapping. Least Squares Approach to SLAM. Cyrill Stachniss

Graphbased. Kalman filter. Particle filter. Three Main SLAM Paradigms. Robot Mapping. Least Squares Approach to SLAM. Least Squares in General

David G. Luenberger Yinyu Ye. Linear and Nonlinear. Programming. Fourth Edition. ö Springer

LEARNING-BASED ADAPTIVE TRANSMISSION FOR LIMITED FEEDBACK MULTIUSER MIMO-OFDM. Alberto Rico-Alvariño and Robert W. Heath Jr.

Lecture 3.3 Robust estimation with RANSAC. Thomas Opsahl

Two-View Geometry (Course 23, Lecture D)

Supplementary Material : Partial Sum Minimization of Singular Values in RPCA for Low-Level Vision

Learning Low-rank Transformations: Algorithms and Applications. Qiang Qiu Guillermo Sapiro

Implementation Of Quadratic Rotation Decomposition Based Recursive Least Squares Algorithm

Numerical schemes for Hamilton-Jacobi equations, control problems and games

Programming, numerics and optimization

Lecture 17 Sparse Convex Optimization

Machine Learning / Jan 27, 2010

THE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS. Alekseí Yu. Chekunov. 1. Introduction

King Fahd University of Petroleum & Minerals. Electrical Engineering Department. EE 575 Information Theory

x = 12 x = 12 1x = 16

Introduction to Topics in Machine Learning

CSE 5526: Introduction to Neural Networks Radial Basis Function (RBF) Networks

LAPACK. Linear Algebra PACKage. Janice Giudice David Knezevic 1

Fast-Lipschitz Optimization

Luca Schenato Workshop on cooperative multi agent systems Pisa, 6/12/2007

EECS 442 Computer vision. Fitting methods

Lecture 4 Duality and Decomposition Techniques

THE COMPUTER MODELLING OF GLUING FLAT IMAGES ALGORITHMS. Alekseí Yu. Chekunov. 1. Introduction

Image deblurring by multigrid methods. Department of Physics and Mathematics University of Insubria

Algebraic Iterative Methods for Computed Tomography

The end of affine cameras

Physical Layer Security from Inter-Session Interference in Large Wireless Networks

Basis Functions. Volker Tresp Summer 2017

Contents. Implementing the QR factorization The algebraic eigenvalue problem. Applied Linear Algebra in Geoscience Using MATLAB

(Creating Arrays & Matrices) Applied Linear Algebra in Geoscience Using MATLAB

IE598 Big Data Optimization Summary Nonconvex Optimization

An Improved Measurement Placement Algorithm for Network Observability

3D Computer Vision. Dense 3D Reconstruction II. Prof. Didier Stricker. Christiano Gava

Optimum Array Processing

Contents. I The Basic Framework for Stationary Problems 1

Outlier Pursuit: Robust PCA and Collaborative Filtering

FMA901F: Machine Learning Lecture 3: Linear Models for Regression. Cristian Sminchisescu

Detecting Burnscar from Hyperspectral Imagery via Sparse Representation with Low-Rank Interference

Capturing, Modeling, Rendering 3D Structures

Signal Reconstruction from Sparse Representations: An Introdu. Sensing

Transcription:

Advanced Topics in Digital Communications Spezielle Methoden der digitalen Datenübertragung Dr.-Ing. Carsten Bockelmann Institute for Telecommunications and High-Frequency Techniques Department of Communications Engineering Room: SPT C3160, Phone: 0421/218-62386 bockelmann@ant.uni-bremen.de Lecture Thursday, 10:00 12:00 in N3130 Exercise Wednesday, 14:00 16:00 in N1250 Dates for exercises will be announced during lectures. Tutor Tobias Monsees Room: SPT C3220 Phone 218-62407 tmonsees@ant.uni-bremen.de www.ant.uni-bremen.de/courses/atdc/

Outline Part 1: Linear Algebra Eigenvalues and eigenvectors, pseudo inverse Decompositions (QR, unitary matrices, singular value, Cholesky ) Part 2: Basics and Preliminaries Motivating systems with Multiple Inputs and Multiple Outputs (multiple access techniques) General classification and description of MIMO systems (SIMO, MISO, MIMO) Mobile Radio Channel Part 3: Information Theory for MIMO Systems Repetition of IT basics, channel capacity for SISO AWGN channel Extension to SISO fading channels Generalization for the MIMO case Part 4: Multiple Antenna Systems SIMO: diversity gain, beamforming at receiver MISO: space-time coding, beamforming at transmitter MIMO: BLAST with detection strategies Influence of channel (correlation) Part 5: Relaying Systems Basic relaying structures Relaying protocols and exemplary configurations Outline 2

Outline Part 6: In Network Processing Basic of distributed processing INP approach Part 7: Compressive Sensing Motivating Sampling below Nyquist Reconstruction principles and algorithms Applications Outline 3

Problem Statement Network of nodes perform different measurements of same quantity Measurements need to be combined to estimate unknown source Examples: SC1 SC2 (x j ; y j ; z j ) m j m 1 (x 1 ; y 1 ; z 1 ) F (x; y; z) m 2 (x 2 ; y 2 ; z 2 ) SC3 Wireless sensor network measuring temperature, humidity, salinity, etc of a medium estimate field distribution Network of base stations ( small cells ) receiving different replicas of radio signals recover data transmitted by users 4

Example: Linear Estimation Assume: network of J nodes connected by inter-node-links every node j makes different observation x j of the same quantity s observations are linearly distorted with additional noise term s 1 x 1 3 4 x j = H j s + n j 2 x j j Objective: Estimate s by fully utilizing all information available at J observations Note: H j can be any matrix, e.g. underdetermined system 5

Non-cooperative approach Every node performs an individual estimation For example: Least-Square estimation Pros: No inter-node cooperation required Cons: Per node the equation system might be underdetermined (H j is rank deficient), e.g Multi-User MIMO Estimations differ from node to node no consensus between nodes further processing required to reconstruct (the single) s 6

Centralized processing Every node passes on observations and channel H j to fusion center Fusion center performs centralized estimation using all observations and Pros: SotA algorithms yield optimum solution Cons: Communication effort routing protocols needed / single point of failure 7

In-Network-Processing Idea: Calculate a function within a network, e.g. average, MMSE/LS estimation Consensus based: Identical solution at all nodes Challenge: Desing of algorithms yielding optimum solution, but being signalling efficient and robust Our approach: Start from a mathematical optimization framework 8

In-Network-Processing Rewrite centralized optimization problem as Question: How to solve it by In-Network-Processing? Solution: It can be shown that Algorithmic approach: Augmented Lagrangian method / Method of Multipliers [Nocedal&Wright] Iterative algorithm with cooperation between neighbor nodes only For l 2 -norm: decomposition into local, but coupled problems 9

Distributed Consensus-Based Estimation (DiCE) Idea: Introduce node specific auxiliary variables z j and decouple variables s j Variables s j, z j : iteratively calculated and broadcast Lagrangian multipliers λλ jjjj : iteratively calculated and unicast Algorithm: Convergence: z j (k + 1) = f( j;i (k); s i (k)); Perfect links: DICE converges to optimum solution Noisy links: DICE converges in mean sense to optimum solution Link failures: DICE converges to optimum solution Main issue: Reduction of inter-node signaling i 2 N j [ j s j (k + 1) = f( i;j (k); z i (k + 1); x j ; H j ) i;j (k + 1) = f( i;j (k); s j (k + 1); z i (k + 1)) k: iteration index 10

DICE with reduced inter-node signaling Reduced Overhead DiCE (RO-DiCE) [Shin 13] Avoid exchange of edge specific Lagrange multipliers λλ jjii avoid unicast transmissions, broadcast transmissions only Fast-DiCE [Xu 14] Acceleration of DICE by using Nesterov s optimal gradient descend method Less #of iterations required for same estimation quality reduced inter-node communication [Shin 13] B.-S. Shin, H. Paul, D. Wübben, A. Dekorsy, Reduced Overhead Distributed Consensus-Based Estimation Algorithm, IWCPM 2013, Atlanta, GA, USA, December 2013. [Xu 14] G. Xu, H. Paul, D. Wübben, A. Dekorsy, Fast Distributed Consensus-based Estimation for Cooperative Wireless Sensor Networks, WSA 2014, Erlangen, Germany, March 2014. 11

4 xx 4, dd 4 Kernel-based distributed estimation xx 3, dd 3 3 1 xx 1, dd 1 ff(xx) 2 xx 2, dd 2 xx NN ss, dd NN ss NN ss Sensor network with NN ss nodes observes nonlinear function ff(xx) SN measures noisy value of function output dd jj at specific positions xx jj: dd jj = ff (xx jj ) + nn jj NN ss training samples {(xx jj, dd jj )} Objective: Estimate ff(xx) at arbitrary positions xx by INP within sensor network Distributed nonlinear regression on training samples { xx jj, dd jj } Use kernel methods for nonlinear regression: NN SS ff = arg min dd jj ff(xx jj ) ff H jj=1 12

Kernel DiCE Representer theorem states that optimal solution ff, i.e., the function estimate ff can be represented by linear superposition of positive semidefinite kernels κκ(, ): ff ( ) = NN ll=1 ww ll κκ xx ll, xx ll : dictionary (size NN) Kernel represents inner product between argument vectors in some higher dimensional (possibly infinite-dimensional) kernel Hilbert space Objective: Find weights ww ll for given kernel (e.g., Gaussian) and dictionary (e.g., sensor locations) estimation problem Approach: Calculate per-node weights using DiCE, enforcing consensus among nodes Reconstruction of field possible at every node using above equation 13

Example: Estimate diffusion by moving sensors ff(xx) describes a diffusion process with 2 sources sensors can move optimize sensor movement to improve estimation quality B.-S. Shin, H. Paul, A. Dekorsy: Spatial Field Reconstruction with Distributed Kernel Least Squares in Mobile Sensor Networks. Accepted for publication at SCC 2017 14