Sparse stochastic processes
|
|
- Brianna Richardson
- 5 years ago
- Views:
Transcription
1 BIOMEDICAL IMAGING GROUP (BIG) LABORATOIRE D IMAGERIE BIOMEDICALE EPFL LIB Bât. BM CH 115 Lausanne Switzerland Téléphone : Fax : Site web : Michael.Unser@epfl.ch Dr. Michael Liebling Sparse stochastic processes Biological Imaging Center California Inst. of Technology Mail Code Pasadena, CA 91125, USA Part I: Introduction Lausanne, August 19, 24 Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Dear Dr. Liebling, I am pleased to inform you that you were selected to the receive the 24 Research Award of the Swiss ofkorea Biomedical Engineering for your thesis work On Fresnelets, interference InvitedSociety seminar, Advanced Inst. of Science & Technology (KAIST), July 17-1, Seoul, Korea fringes, and digital holography. The award will be presented during the general assembly of the SSBE, September 3, Zurich, Switzerland. Please, lets us know if 1) you will be present to receive the award, 2) you would be willing to give a 1 minutes presentation of the work during the general assembly. 2th century statistical signal processing Hypothesis: Signal = stationary Gaussian process The award comes with a cash prize of 1.- CHF. Karhunen-Loève transform (KLT) is optimal for compression Would you please send your banking information to the treasurer of the SSBE, Uli Diermann ( uli.diermann@bfh.ch), so that he can transfer the cash prize to your account? I congratulate you on your achievement. DCT asymptotically equivalent to KLT With best regards, (Ahmed-Rao, 1975; U., 194) modulation the cross-modulationloss. improvement system can it conventional Michael Unser, Professor Chairman of the SSBE Award Committee cc: Ralph Mueller, president of the SSBE; Uli Diermann, treasurer Fig. 5. Rate versus distortion for various transforms for a second-order Markov process ( p =.95, N = 256). Gauss- (Pearl et al., IEEE Trans. Com 1972) 2
2 2th century statistical signal processing Hypothesis: Signal = Gaussian process y = Hs + n Noise: i.i.d. Gaussian with variance 2 Signal covariance: C s = E{s s T } Wiener filter is optimal for restoration/denoising s LMMSE = C s H T HC s H T + 2 I 1 y = FWiener y m L = Cs 1/2 : Whitening filter Wiener (LMMSE) solution = Gauss MMSE = Gauss MAP s MAP = arg min s 12 ky Hsk2 2 {z } Data Log likelihood + kc 1/2 s sk 2 2 {z } Gaussian prior likelihood, quadratic regularization (Tikhonov) klsk Then came wavelets... and sparsity Yves Meyer Alfred Haar 191 // 192 Stéphane Mallat 197- Ingrid Daubechies 1994 Sparsity 26 Applications Compressed sensing Martin Vetterli David Donoho Emmanuel Candès 4
3 Fact 1: Wavelets can outperform Wiener filter w = T (w) 2 w 5 Fact 2: Wavelet coding can outperform jpeg f(x) = X i,k i,k(x) w i,k Wavelet transform Inverse wavelet transform Discarding small coefficients (Shapiro, IEEE-IP 1993) 6
4 Fact 3: l1 schemes can outperform l2 optimization s? = argmin ky Hsk R(s) {z } {z } data consistency regularization Wavelet-domain regularization Wavelet expansion: s = Wv (typically, sparse) Wavelet-domain sparsity-constraint: R(s) =kvk`1 with v = W 1 s (Nowak et al., Daubechies et al. 24) `1 regularization (Total variation) (Rudin-Osher, 1992) R(s) =klsk`1 with L: gradient Compressed sensing/sampling (Candes-Romberg-Tao; Donoho, 26) 7 SPARSE STOCHASTIC MODELS: The step beyond Gaussianity Requirements for a comprehensive statistical framework Backward compatibility Continuous-domain formulation piecewise-smooth signals, translation and scale-invariance, sampling... Predictive power Can wavelets really outperform sinusoidal transforms (KLT)? Statistical justification and refinement of current algorithms Sparsity-promoting regularization, `1 norm minimization
5 1.3 FROM SPLINES TO STOCHASTIC PROCESSES Splines and Legos revisited Higher-order polynomial splines Random splines, innovations, Lévy processes Wavelet analysis of Lévy processes Lévy s synthesis of Brownian motion 9 Splines and Legos revisited Cardinal spline of degree : piecewise-contant f 1 (t) = X k2z f 1 [k] +(t k) +(t) = ( 1, for apple t<1, otherwise. Notion of D-spline: Df 1 (t) = X k2z a 1 [k] (t k) 1
6 B-spline and derivative operator Derivative Df(t) = df(t) dt D F! j! Finite difference operator D d f(t) =f(t) f(t 1) D d F! 1 e j! =( + Df)(t) B-spline of degree +(t) =D d D 1 (t) =D d + (t) l 1 +(t) = + (t) +(t 1) ˆ + (!) = 1 e j! j! Random splines and innovations cardinal Ds(t) = X n a n (t t n )=w(t) non-uniform Random weights {a n } i.i.d. and random knots {t n } (Poisson with rate ) Anti-derivative operators Shift-invariant solution: D 1 '(t) =( + ')(t) = Scale-invariant solution: D 1 '(t) = Z t '( )d Z t 1 '( )d 12
7 Innovation-based synthesis L= d dt =D ) L 1 : integrator (t) L 1 { } =L 1 (t t ) Translation invariance L 1 { } Impulse response (t t ) X a n (t t n ) n Linearity L 1 { } s(t) = X n a n (t t n ) 13 Compound Poisson process Stochastic differential equation Ds(t) =w(t) with boundary condition s() = Innovation: w(t) = X n a n (t t n ) Formal solution s(t) =D 1 w(t) =X n = X n a n D 1 { ( t n)}(t) a n + (t t n ) +( t n ) (impose boundary condition) 14
8 Generalization: Lévy processes Generalized innovations : white Lévy noise with E{w(t)w(t )} = 2 w (t t ) Ds = w (unstable SDE!) s =D 1 1 w, ' 2 S, h',si = hd ',wi White noise (innovation) Lévy process Gaussian Brownian motion (Wiener 1923) Integrator Impulsive w(t) Z t d s(t) Compound Poisson S S (Cauchy) Lévy flight (Paul Lévy circa 193) 15 Decoupling Lévy processes: increments Increment process: u(t) =D d s(t) =D d D 1 w(t) =( + w)(t). Increment process is stationary with autocorrelation function R u ( ) =E{u(t)u(t + )} = + ( +) _ R w (y) = 2 w 1 +( 1) * Discrete increments u[k] = s(k) s(k 1) = hw, +( k)i. u[k] are i.i.d. because { +( k)} are non-overlapping h, i w is independent at every point (white noise) 16
9 Wavelet analysis of Lévy processes Haar wavelets >< Haar(t) = >: 1, for apple t< 1 2 1, for 1 2 apple t<1, otherwise. >< i = >:,,2 i,k(t) =2 i/2 Haar t 2 i k 2 i >< i =1 >: 1, i =2 ( 2, 17 Wavelets as multi-scale derivatives >< i = >:,,2,,2 >< i = >: >< i =1 >: 1, >< i =1 >: 1, i =2 ( 2, i =2 ( 2, i,k =2 i/2 1 D i,k D 1 i,k =2 i/2 1 i,k. Wavelet coefficients of Lévy process Y i,k = hs, i,k i/hs, D i,k i /hd s, i,ki = hw, i,ki 1
10 M-term approximation: wavelets vs. KLT 1 Gaussian Identity KLT Haar Brownian motion Finite rate of innovation 1 1 Identity KLT Haar Compound Poisson Even sparser Identity KLT Haar Lévy flight (Cauchy) Wavelet-based synthesis of Brownian motion White Gaussian noise w = X X Z i,k i,k with Z i,k = hw, i,k i i2z k2z Brownian motion s(t) =D 1 w = X X Z i,k D 1 i,k(t) i2z k2z = X X 2 i/2 1 Z i,k i,k (t) i2z k2z 2
11 Unser and Tafti One puff to follow, one puff to follow, one puff to follow, one puff to follow, one puff to follow, one puff to follow, one puff to follow AUTHOR NAME, affiliation An Introduction to Sparse Stochastic Processes An Introduction to Sparse Stochastic Processes Unser and Tafti PPC. C M Y K Providing a novel approach to sparse stocastic processes, this comprehensive book presents the theory of stochastic process that are ruled by stochastic differential equations, and that admit a parsimonious representation in a matched wavelet-like basis. Two key themes are the statistical property of infinite divisibility, which leads to two distinct types of behavior Gaussian and sparse and the structural link between linear stochastic processes and spline functions, which is exploited to simplify the mathematical analysis. The core of the book is devoted to investigating sparse processes, including a complete description of their transform-domain statistics. The final part develops practical signal-processing algorithms that are based on these models, with special emphasis on biomedical image reconstruction. This is an ideal reference for graduate students and researchers with an interest in signal/image processing, compressed sensing, approximation theory, machine learning, or statistics. M I C H A E L U N S E R is Professor and Director of EPFL s Biomedical Imaging Group, Switzerland. He is a member of the Swiss Academy of Engineering Sciences, a fellow of EURASIP, and a fellow of the IEEE. Michael Unser and Pouya D. Tafti P O U YA D. TA F T I is a researcher at Qlaym BmbH, Düsseldort, and a former member of the Biomedical Imaging Group at EPFL, where he conducted research on the theory and applications of probablilistic models for data. Cover illustration: to follow ebook: web preprint EDEE Course 21 Extension to fractional orders: B-splines L = D (fractional derivative) Ld = + Fourier multiplier: (j!) (fractional differences) + (r) = + r+ Fourier multiplier: (1 1 F!... + (r) = +1 + r+ ( + 1) One-sided power function: e j! j! e j! j! e j! )... F! r+ = 1 ( +1 r, r, r< Degree = 1 22
12 Mandelbrot s fractional Brownian motion Solution of fractional SDE (Blu-U., 27) D s = w where w is Gaussian white noise ) s =I 2 w I 2 : scale-invariant, right-inverse of D Hurst exponent: H = Continuous-domain innovation model Shaping filter Generalized white noise w(x) L 1 { } (appropriate boundary conditions) Whitening operator L{ } Stochastic process s(x), x 2 R d Main outcome: non-gaussian solutions are necessarily sparse (infinitely divisible) Why?... as will explained in next chapters... (invoking powerful theorems in functional analysis: Bochner-Minlos, Gelfand, Schoenberg & Lévy-Khinchine) 24
13 2.1 DECOUPLING OF SPARSE PROCESSES s =L 1 w, w =Ls Discrete approximation of operator Operator-like wavelet analysis 25 Decoupling: Linear combination of samples Input: s(k), k 2 Z d (sampled values) s =L 1 w Discrete approximation of whitening operator: L d L d (x) = X k2z d d L [k] (x k) Discrete increment process: u[k] =L d s(x) x=k =( L w)(x) x=k = h L _ ( k),wi {z } ' Generalized B-spline: L(x) =L d L 1 (x) A-to-D translator 26
14 Decoupling: Wavelet analysis Ls = w Generalized operator-like wavelets: i(x) =L i(x) (Khalidov-U. 26, Ward-U. JFAA 213) Operator-like wavelet analysis of sparse process: h i ( x ),si =hl i( x ),si =h i ( x ), Ls)i =h i ( x ),wi =( _ i w)(x ) {z } ' 27
Sparse stochastic processes: A statistical framework for modern signal processing
Sparse stochastic processes: A statistical framework for modern signal processing Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Plenary talk, Int. Conf. Syst. Sig. Im. Proc. (IWSSIP),
More informationA Unified Formulation of Gaussian Versus Sparse Stochastic Processes
A Unified Formulation of Gaussian Versus Sparse Stochastic Processes Michael Unser, Pouya Tafti and Qiyu Sun EPFL, UCF Appears in IEEE Trans. on Information Theory, 2014 Presented by Liming Wang M. Unser,
More informationSparse stochastic processes with application to biomedical imaging
Sparse stochastic processes with application to biomedical imaging Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Plenary talk, ICIAR, October -, Vilamoura, Algarve, Portugal Variational
More informationSparse stochastic processes
Sparse stochastic processes Part II: Functional characterization of sparse stochastic processes Prof. Michael Unser, LIB Ecole CIMPA, Analyse et Probabilité, Abidjan 204 oad map for theory of sparse processes
More informationAn introduction to sparse stochastic
An introduction to sparse stochastic processes Michael Unser and Pouya Tafti Copyright 2012 Michael Unser and Pouya Tafti. These materials are protected by copyright under the Attribution-NonCommercial-NoDerivs
More informationSparse modeling and the resolution of inverse problems in biomedical imaging
Sparse modeling and the resolution of inverse problems in biomedical imaging Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Plenary talk, IEEE Int. Symp. Biomedical Imaging (ISBI 5),
More informationThe rise and fall of regularizations: l1 vs. l2 representer theorems
The rise and fall of regularizations: l1 vs. l2 representer theorems Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Joint work with Julien Fageot, John-Paul Ward, and Harshit Gupta
More informationStochastic Models for Sparse and Piecewise-Smooth Signals
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 59, NO. 3, MARCH 2011 989 Stochastic Models for Sparse and Piecewise-Smooth Signals Michael Unser, Fellow, IEEE, and Pouya Dehghani Tafti, Member, IEEE Abstract
More informationVersion beta 1.6 (December 6, 2013). This is a preproduction manuscript of a work protected by copyright. No reproduction of any part may take place
Version beta 1.6 (December 6, 2013). This is a preproduction manuscript of a work protected by copyright. No reproduction of any part may take place without the written permission of Cambridge University
More informationMMSEEstimationofSparseLévyProcesses
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 61, NO. 1, JANUARY 1, 2013 137 MMSEEstimationofSparseLévyProcesses Ulugbek S. Kamilov, Student Member, IEEE, Pedram Pad, Student Member, IEEE, Arash Amini,
More informationThe Analog Formulation of Sparsity Implies Infinite Divisibility and Rules Out Bernoulli-Gaussian Priors
The Analog Formulation of Sparsity Implies Infinite Divisibility and ules Out Bernoulli-Gaussian Priors Arash Amini, Ulugbek S. Kamilov, and Michael Unser Biomedical Imaging Group BIG) École polytechnique
More informationA unified formulation of Gaussian vs. sparse stochastic processes Part I: Continuous-domain theory
A unified formulation of Gaussian vs. 1 sparse stochastic processes Part I: Continuous-domain theory Michael Unser, Fellow, IEEE, Pouya Tafti, Member, IEEE, and Qiyu Sun arxiv:1108.6150v2 [cs.it] 5 Oct
More informationSparse linear models
Sparse linear models Optimization-Based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_spring16 Carlos Fernandez-Granda 2/22/2016 Introduction Linear transforms Frequency representation Short-time
More informationWavelets and differential operators: from fractals to Marr!s primal sketch
Wavelets and differential operators: from fractals to Marrs primal sketch Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Joint work with Pouya Tafti and Dimitri Van De Ville Plenary
More informationMMSE Denoising of 2-D Signals Using Consistent Cycle Spinning Algorithm
Denoising of 2-D Signals Using Consistent Cycle Spinning Algorithm Bodduluri Asha, B. Leela kumari Abstract: It is well known that in a real world signals do not exist without noise, which may be negligible
More informationWavelet Footprints: Theory, Algorithms, and Applications
1306 IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 51, NO. 5, MAY 2003 Wavelet Footprints: Theory, Algorithms, and Applications Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli, Fellow, IEEE Abstract
More informationNew representer theorems: From compressed sensing to deep learning
New representer theorems: From compressed sensing to deep learning Michael Unser Biomedical Imaging Group EPFL, Lausanne, Switzerland Joint work with Julien Fageot, John-Paul Ward and Kyong Jin Mathematisches
More informationStrengthened Sobolev inequalities for a random subspace of functions
Strengthened Sobolev inequalities for a random subspace of functions Rachel Ward University of Texas at Austin April 2013 2 Discrete Sobolev inequalities Proposition (Sobolev inequality for discrete images)
More informationSparse linear models and denoising
Lecture notes 4 February 22, 2016 Sparse linear models and denoising 1 Introduction 1.1 Definition and motivation Finding representations of signals that allow to process them more effectively is a central
More informationAN INTRODUCTION TO COMPRESSIVE SENSING
AN INTRODUCTION TO COMPRESSIVE SENSING Rodrigo B. Platte School of Mathematical and Statistical Sciences APM/EEE598 Reverse Engineering of Complex Dynamical Networks OUTLINE 1 INTRODUCTION 2 INCOHERENCE
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 12 Introduction to Wavelets Last Time Started with STFT Heisenberg Boxes Continue and move to wavelets Ham exam -- see Piazza post Please register at www.eastbayarc.org/form605.htm
More informationState Space Representation of Gaussian Processes
State Space Representation of Gaussian Processes Simo Särkkä Department of Biomedical Engineering and Computational Science (BECS) Aalto University, Espoo, Finland June 12th, 2013 Simo Särkkä (Aalto University)
More informationLessons in Estimation Theory for Signal Processing, Communications, and Control
Lessons in Estimation Theory for Signal Processing, Communications, and Control Jerry M. Mendel Department of Electrical Engineering University of Southern California Los Angeles, California PRENTICE HALL
More informationSparse & Redundant Signal Representation, and its Role in Image Processing
Sparse & Redundant Signal Representation, and its Role in Michael Elad The CS Department The Technion Israel Institute of technology Haifa 3000, Israel Wave 006 Wavelet and Applications Ecole Polytechnique
More informationSymmetric Wavelet Tight Frames with Two Generators
Symmetric Wavelet Tight Frames with Two Generators Ivan W. Selesnick Electrical and Computer Engineering Polytechnic University 6 Metrotech Center, Brooklyn, NY 11201, USA tel: 718 260-3416, fax: 718 260-3906
More informationOptimized Compact-support Interpolation Kernels
Optimized Compact-support Interpolation Kernels Ramtin Madani, Student Member, IEEE, Ali Ayremlou, Student Member, IEEE, Arash Amini, Farrokh Marvasti, Senior Member, IEEE, Abstract In this paper, we investigate
More informationMultichannel Sampling of Signals with Finite. Rate of Innovation
Multichannel Sampling of Signals with Finite 1 Rate of Innovation Hojjat Akhondi Asl, Pier Luigi Dragotti and Loic Baboulaz EDICS: DSP-TFSR, DSP-SAMP. Abstract In this paper we present a possible extension
More informationWhich wavelet bases are the best for image denoising?
Which wavelet bases are the best for image denoising? Florian Luisier a, Thierry Blu a, Brigitte Forster b and Michael Unser a a Biomedical Imaging Group (BIG), Ecole Polytechnique Fédérale de Lausanne
More informationDirectionlets. Anisotropic Multi-directional Representation of Images with Separable Filtering. Vladan Velisavljević Deutsche Telekom, Laboratories
Directionlets Anisotropic Multi-directional Representation of Images with Separable Filtering Vladan Velisavljević Deutsche Telekom, Laboratories Google Inc. Mountain View, CA October 2006 Collaborators
More informationGaussian and Sparse Processes Are Limits of Generalized Poisson Processes
1 Gaussian and Sparse Processes Are Limits of Generalized Poisson Processes Julien Fageot, Virginie Uhlmann, and Michael Unser arxiv:1702.05003v1 [math.pr] 16 Feb 2017 Abstract The theory of sparse stochastic
More informationAdaptive wavelet decompositions of stochastic processes and some applications
Adaptive wavelet decompositions of stochastic processes and some applications Vladas Pipiras University of North Carolina at Chapel Hill SCAM meeting, June 1, 2012 (joint work with G. Didier, P. Abry)
More information( nonlinear constraints)
Wavelet Design & Applications Basic requirements: Admissibility (single constraint) Orthogonality ( nonlinear constraints) Sparse Representation Smooth functions well approx. by Fourier High-frequency
More informationDirect Learning: Linear Classification. Donglin Zeng, Department of Biostatistics, University of North Carolina
Direct Learning: Linear Classification Logistic regression models for classification problem We consider two class problem: Y {0, 1}. The Bayes rule for the classification is I(P(Y = 1 X = x) > 1/2) so
More informationStatistical and Adaptive Signal Processing
r Statistical and Adaptive Signal Processing Spectral Estimation, Signal Modeling, Adaptive Filtering and Array Processing Dimitris G. Manolakis Massachusetts Institute of Technology Lincoln Laboratory
More informationSignal Recovery from Permuted Observations
EE381V Course Project Signal Recovery from Permuted Observations 1 Problem Shanshan Wu (sw33323) May 8th, 2015 We start with the following problem: let s R n be an unknown n-dimensional real-valued signal,
More informationCS 229r: Algorithms for Big Data Fall Lecture 19 Nov 5
CS 229r: Algorithms for Big Data Fall 215 Prof. Jelani Nelson Lecture 19 Nov 5 Scribe: Abdul Wasay 1 Overview In the last lecture, we started discussing the problem of compressed sensing where we are given
More informationContents. Acknowledgments
Table of Preface Acknowledgments Notation page xii xx xxi 1 Signals and systems 1 1.1 Continuous and discrete signals 1 1.2 Unit step and nascent delta functions 4 1.3 Relationship between complex exponentials
More informationWavelets, Filter Banks and Multiresolution Signal Processing
Wavelets, Filter Banks and Multiresolution Signal Processing It is with logic that one proves; it is with intuition that one invents. Henri Poincaré Introduction - 1 A bit of history: from Fourier to Haar
More informationNumerical Analysis for Statisticians
Kenneth Lange Numerical Analysis for Statisticians Springer Contents Preface v 1 Recurrence Relations 1 1.1 Introduction 1 1.2 Binomial CoefRcients 1 1.3 Number of Partitions of a Set 2 1.4 Horner's Method
More informationThe Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems
668 IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: ANALOG AND DIGITAL SIGNAL PROCESSING, VOL 49, NO 10, OCTOBER 2002 The Discrete Kalman Filtering of a Class of Dynamic Multiscale Systems Lei Zhang, Quan
More informationAn Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems
An Overview of Sparsity with Applications to Compression, Restoration, and Inverse Problems Justin Romberg Georgia Tech, School of ECE ENS Winter School January 9, 2012 Lyon, France Applied and Computational
More informationAdaptive Filtering. Squares. Alexander D. Poularikas. Fundamentals of. Least Mean. with MATLABR. University of Alabama, Huntsville, AL.
Adaptive Filtering Fundamentals of Least Mean Squares with MATLABR Alexander D. Poularikas University of Alabama, Huntsville, AL CRC Press Taylor & Francis Croup Boca Raton London New York CRC Press is
More informationInvariant Scattering Convolution Networks
Invariant Scattering Convolution Networks Joan Bruna and Stephane Mallat Submitted to PAMI, Feb. 2012 Presented by Bo Chen Other important related papers: [1] S. Mallat, A Theory for Multiresolution Signal
More informationIntroduction to Wavelet. Based on A. Mukherjee s lecture notes
Introduction to Wavelet Based on A. Mukherjee s lecture notes Contents History of Wavelet Problems of Fourier Transform Uncertainty Principle The Short-time Fourier Transform Continuous Wavelet Transform
More informationCambridge University Press The Mathematics of Signal Processing Steven B. Damelin and Willard Miller Excerpt More information
Introduction Consider a linear system y = Φx where Φ can be taken as an m n matrix acting on Euclidean space or more generally, a linear operator on a Hilbert space. We call the vector x a signal or input,
More informationData Compression and Harmonic Analysis
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 44, NO. 6, OCTOBER 1998 2435 Data Compression and Harmonic Analysis David L. Donoho, Martin Vetterli, Fellow, IEEE, R. A. DeVore, and Ingrid Daubechies, Senior
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 218. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationVariations. ECE 6540, Lecture 10 Maximum Likelihood Estimation
Variations ECE 6540, Lecture 10 Last Time BLUE (Best Linear Unbiased Estimator) Formulation Advantages Disadvantages 2 The BLUE A simplification Assume the estimator is a linear system For a single parameter
More informationAn Introduction to Sparse Approximation
An Introduction to Sparse Approximation Anna C. Gilbert Department of Mathematics University of Michigan Basic image/signal/data compression: transform coding Approximate signals sparsely Compress images,
More informationProbability and Statistics for Final Year Engineering Students
Probability and Statistics for Final Year Engineering Students By Yoni Nazarathy, Last Updated: May 24, 2011. Lecture 6p: Spectral Density, Passing Random Processes through LTI Systems, Filtering Terms
More informationLarge-Scale L1-Related Minimization in Compressive Sensing and Beyond
Large-Scale L1-Related Minimization in Compressive Sensing and Beyond Yin Zhang Department of Computational and Applied Mathematics Rice University, Houston, Texas, U.S.A. Arizona State University March
More informationSIMON FRASER UNIVERSITY School of Engineering Science
SIMON FRASER UNIVERSITY School of Engineering Science Course Outline ENSC 810-3 Digital Signal Processing Calendar Description This course covers advanced digital signal processing techniques. The main
More informationKarhunen-Loève Expansions of Lévy Processes
Karhunen-Loève Expansions of Lévy Processes Daniel Hackmann June 2016, Barcelona supported by the Austrian Science Fund (FWF), Project F5509-N26 Daniel Hackmann (JKU Linz) KLE s of Lévy Processes 0 / 40
More informationEE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6)
EE 367 / CS 448I Computational Imaging and Display Notes: Image Deconvolution (lecture 6) Gordon Wetzstein gordon.wetzstein@stanford.edu This document serves as a supplement to the material discussed in
More informationWE begin by briefly reviewing the fundamentals of the dual-tree transform. The transform involves a pair of
Gabor wavelet analysis and the fractional Hilbert transform Kunal Narayan Chaudhury and Michael Unser Biomedical Imaging Group, Ecole Polytechnique Fédérale de Lausanne (EPFL), Switzerland ABSTRACT We
More informationSparsity Measure and the Detection of Significant Data
Sparsity Measure and the Detection of Significant Data Abdourrahmane Atto, Dominique Pastor, Grégoire Mercier To cite this version: Abdourrahmane Atto, Dominique Pastor, Grégoire Mercier. Sparsity Measure
More informationProperties of the Autocorrelation Function
Properties of the Autocorrelation Function I The autocorrelation function of a (real-valued) random process satisfies the following properties: 1. R X (t, t) 0 2. R X (t, u) =R X (u, t) (symmetry) 3. R
More informationStatistics of Stochastic Processes
Prof. Dr. J. Franke All of Statistics 4.1 Statistics of Stochastic Processes discrete time: sequence of r.v...., X 1, X 0, X 1, X 2,... X t R d in general. Here: d = 1. continuous time: random function
More informationAn Introduction to Sparse Representations and Compressive Sensing. Part I
An Introduction to Sparse Representations and Compressive Sensing Part I Paulo Gonçalves CPE Lyon - 4ETI - Cours Semi-Optionnel Méthodes Avancées pour le Traitement des Signaux 2014 Objectifs Part I The
More informationDigital Object Identifier /MSP
DIGITAL VISION Sampling Signals from a Union of Subspaces [A new perspective for the extension of this theory] [ Yue M. Lu and Minh N. Do ] Our entire digital revolution depends on the sampling process,
More informationCourse content (will be adapted to the background knowledge of the class):
Biomedical Signal Processing and Signal Modeling Lucas C Parra, parra@ccny.cuny.edu Departamento the Fisica, UBA Synopsis This course introduces two fundamental concepts of signal processing: linear systems
More informationNear Optimal Signal Recovery from Random Projections
1 Near Optimal Signal Recovery from Random Projections Emmanuel Candès, California Institute of Technology Multiscale Geometric Analysis in High Dimensions: Workshop # 2 IPAM, UCLA, October 2004 Collaborators:
More informationWavelet Based Image Restoration Using Cross-Band Operators
1 Wavelet Based Image Restoration Using Cross-Band Operators Erez Cohen Electrical Engineering Department Technion - Israel Institute of Technology Supervised by Prof. Israel Cohen 2 Layout Introduction
More informationCompressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery
Compressibility of Infinite Sequences and its Interplay with Compressed Sensing Recovery Jorge F. Silva and Eduardo Pavez Department of Electrical Engineering Information and Decision Systems Group Universidad
More informationIEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 52, NO. 12, DECEMBER
IEEE TRANSACTIONS ON INFORMATION THEORY, VOL. 52, NO. 12, DECEMBER 2006 5177 The Distributed Karhunen Loève Transform Michael Gastpar, Member, IEEE, Pier Luigi Dragotti, Member, IEEE, and Martin Vetterli,
More informationCompressed Sensing and Sparse Recovery
ELE 538B: Sparsity, Structure and Inference Compressed Sensing and Sparse Recovery Yuxin Chen Princeton University, Spring 217 Outline Restricted isometry property (RIP) A RIPless theory Compressed sensing
More informationSuper-resolution via Convex Programming
Super-resolution via Convex Programming Carlos Fernandez-Granda (Joint work with Emmanuel Candès) Structure and Randomness in System Identication and Learning, IPAM 1/17/2013 1/17/2013 1 / 44 Index 1 Motivation
More informationFRACTIONAL BROWNIAN MODELS FOR VECTOR FIELD DATA. Pouya Dehghani Tafti, Member, IEEE, and Michael Unser, Fellow, IEEE
FRACTIONAL BROWNIAN MODELS FOR VECTOR FIELD DATA Pouya Dehghani Tafti, Member, IEEE, and Michael Unser, Fellow, IEEE Laboratoire d imagerie biomédicale, Ecole Polytechnique Fédérale de Lausanne (EPFL),
More informationBayesian Methods for Sparse Signal Recovery
Bayesian Methods for Sparse Signal Recovery Bhaskar D Rao 1 University of California, San Diego 1 Thanks to David Wipf, Jason Palmer, Zhilin Zhang and Ritwik Giri Motivation Motivation Sparse Signal Recovery
More informationPROBABILITY: LIMIT THEOREMS II, SPRING HOMEWORK PROBLEMS
PROBABILITY: LIMIT THEOREMS II, SPRING 15. HOMEWORK PROBLEMS PROF. YURI BAKHTIN Instructions. You are allowed to work on solutions in groups, but you are required to write up solutions on your own. Please
More informationDo probabilists flip coins all the time?
Include Only If Paper Has a Subtitle Department of Mathematics and Statistics Do probabilists flip coins all the time? Math Graduate Seminar February 17, 2009 Outline 1 Bernoulli random variables and Rademacher
More informationLinear Models for Regression CS534
Linear Models for Regression CS534 Prediction Problems Predict housing price based on House size, lot size, Location, # of rooms Predict stock price based on Price history of the past month Predict the
More informationCOMPRESSED SENSING IN PYTHON
COMPRESSED SENSING IN PYTHON Sercan Yıldız syildiz@samsi.info February 27, 2017 OUTLINE A BRIEF INTRODUCTION TO COMPRESSED SENSING A BRIEF INTRODUCTION TO CVXOPT EXAMPLES A Brief Introduction to Compressed
More informationThe Sparsity Gap. Joel A. Tropp. Computing & Mathematical Sciences California Institute of Technology
The Sparsity Gap Joel A. Tropp Computing & Mathematical Sciences California Institute of Technology jtropp@acm.caltech.edu Research supported in part by ONR 1 Introduction The Sparsity Gap (Casazza Birthday
More informationCompressed Sensing in Astronomy
Compressed Sensing in Astronomy J.-L. Starck CEA, IRFU, Service d'astrophysique, France jstarck@cea.fr http://jstarck.free.fr Collaborators: J. Bobin, CEA. Introduction: Compressed Sensing (CS) Sparse
More informationLearning MMSE Optimal Thresholds for FISTA
MITSUBISHI ELECTRIC RESEARCH LABORATORIES http://www.merl.com Learning MMSE Optimal Thresholds for FISTA Kamilov, U.; Mansour, H. TR2016-111 August 2016 Abstract Fast iterative shrinkage/thresholding algorithm
More informationBrownian Motion. 1 Definition Brownian Motion Wiener measure... 3
Brownian Motion Contents 1 Definition 2 1.1 Brownian Motion................................. 2 1.2 Wiener measure.................................. 3 2 Construction 4 2.1 Gaussian process.................................
More informationGabor wavelet analysis and the fractional Hilbert transform
Gabor wavelet analysis and the fractional Hilbert transform Kunal Narayan Chaudhury and Michael Unser (presented by Dimitri Van De Ville) Biomedical Imaging Group, Ecole Polytechnique Fédérale de Lausanne
More informationLecture 9. d N(0, 1). Now we fix n and think of a SRW on [0,1]. We take the k th step at time k n. and our increments are ± 1
Random Walks and Brownian Motion Tel Aviv University Spring 011 Lecture date: May 0, 011 Lecture 9 Instructor: Ron Peled Scribe: Jonathan Hermon In today s lecture we present the Brownian motion (BM).
More informationColor Scheme. swright/pcmi/ M. Figueiredo and S. Wright () Inference and Optimization PCMI, July / 14
Color Scheme www.cs.wisc.edu/ swright/pcmi/ M. Figueiredo and S. Wright () Inference and Optimization PCMI, July 2016 1 / 14 Statistical Inference via Optimization Many problems in statistical inference
More informationMGA Tutorial, September 08, 2004 Construction of Wavelets. Jan-Olov Strömberg
MGA Tutorial, September 08, 2004 Construction of Wavelets Jan-Olov Strömberg Department of Mathematics Royal Institute of Technology (KTH) Stockholm, Sweden Department of Numerical Analysis and Computer
More informationSummary of the ISBI 2013 Grand Challenge on 3D Deconvolution Microscopy
Summary of the ISBI 23 Grand Challenge on 3D Deconvolution Microscopy Cédric Vonesch & Stamatios Lefkimmiatis Biomedical Imaging Group EPFL, Lausanne, Switzerland April 7, 23 Thanks to our sponsors! ISBI
More informationModule 4. Multi-Resolution Analysis. Version 2 ECE IIT, Kharagpur
Module 4 Multi-Resolution Analysis Lesson Multi-resolution Analysis: Discrete avelet Transforms Instructional Objectives At the end of this lesson, the students should be able to:. Define Discrete avelet
More informationIntroduction to Biomedical Engineering
Introduction to Biomedical Engineering Biosignal processing Kung-Bin Sung 6/11/2007 1 Outline Chapter 10: Biosignal processing Characteristics of biosignals Frequency domain representation and analysis
More informationLinear Models for Regression CS534
Linear Models for Regression CS534 Example Regression Problems Predict housing price based on House size, lot size, Location, # of rooms Predict stock price based on Price history of the past month Predict
More informationFinding Increment Statistics on various types of Wavelets under 1-D Fractional Brownian Motion Synthesis
Finding Increment Statistics on various types of Wavelets under 1-D Fractional Brownian Motion Synthesis Abhishek Sharma 1, Bhubneshwar Sharma 2 1 M.Tech student, Department of ECE Eternal University,
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationWavelets and Multiresolution Processing
Wavelets and Multiresolution Processing Wavelets Fourier transform has it basis functions in sinusoids Wavelets based on small waves of varying frequency and limited duration In addition to frequency,
More informationVarious signal sampling and reconstruction methods
Various signal sampling and reconstruction methods Rolands Shavelis, Modris Greitans 14 Dzerbenes str., Riga LV-1006, Latvia Contents Classical uniform sampling and reconstruction Advanced sampling and
More informationSparse Legendre expansions via l 1 minimization
Sparse Legendre expansions via l 1 minimization Rachel Ward, Courant Institute, NYU Joint work with Holger Rauhut, Hausdorff Center for Mathematics, Bonn, Germany. June 8, 2010 Outline Sparse recovery
More informationLong-Range Dependence and Self-Similarity. c Vladas Pipiras and Murad S. Taqqu
Long-Range Dependence and Self-Similarity c Vladas Pipiras and Murad S. Taqqu January 24, 2016 Contents Contents 2 Preface 8 List of abbreviations 10 Notation 11 1 A brief overview of times series and
More informationMultiplicative and Additive Perturbation Effects on the Recovery of Sparse Signals on the Sphere using Compressed Sensing
Multiplicative and Additive Perturbation Effects on the Recovery of Sparse Signals on the Sphere using Compressed Sensing ibeltal F. Alem, Daniel H. Chae, and Rodney A. Kennedy The Australian National
More informationDigital Image Processing
Digital Image Processing, 2nd ed. Digital Image Processing Chapter 7 Wavelets and Multiresolution Processing Dr. Kai Shuang Department of Electronic Engineering China University of Petroleum shuangkai@cup.edu.cn
More informationProbability Models in Electrical and Computer Engineering Mathematical models as tools in analysis and design Deterministic models Probability models
Probability Models in Electrical and Computer Engineering Mathematical models as tools in analysis and design Deterministic models Probability models Statistical regularity Properties of relative frequency
More informationCompressed Sensing and Neural Networks
and Jan Vybíral (Charles University & Czech Technical University Prague, Czech Republic) NOMAD Summer Berlin, September 25-29, 2017 1 / 31 Outline Lasso & Introduction Notation Training the network Applications
More informationOptimal series representations of continuous Gaussian random fields
Optimal series representations of continuous Gaussian random fields Antoine AYACHE Université Lille 1 - Laboratoire Paul Painlevé A. Ayache (Lille 1) Optimality of continuous Gaussian series 04/25/2012
More informationGoing off the grid. Benjamin Recht Department of Computer Sciences University of Wisconsin-Madison
Going off the grid Benjamin Recht Department of Computer Sciences University of Wisconsin-Madison Joint work with Badri Bhaskar Parikshit Shah Gonnguo Tang We live in a continuous world... But we work
More informationWavelets and multiresolution representations. Time meets frequency
Wavelets and multiresolution representations Time meets frequency Time-Frequency resolution Depends on the time-frequency spread of the wavelet atoms Assuming that ψ is centred in t=0 Signal domain + t
More information2D Wavelets. Hints on advanced Concepts
2D Wavelets Hints on advanced Concepts 1 Advanced concepts Wavelet packets Laplacian pyramid Overcomplete bases Discrete wavelet frames (DWF) Algorithme à trous Discrete dyadic wavelet frames (DDWF) Overview
More informationINTRODUCTION TO. Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR)
INTRODUCTION TO WAVELETS Adapted from CS474/674 Prof. George Bebis Department of Computer Science & Engineering University of Nevada (UNR) CRITICISM OF FOURIER SPECTRUM It gives us the spectrum of the
More information