Adaptive Linear Filtering Using Interior Point. Optimization Techniques. Lecturer: Tom Luo
|
|
- Ann Underwood
- 5 years ago
- Views:
Transcription
1 Adaptive Linear Filtering Using Interior Point Optimization Techniques Lecturer: Tom Luo
2 Overview A. Interior Point Least Squares (IPLS) Filtering Introduction to IPLS Recursive update of IPLS Convergence/transient analysis of IPLS B. Applications System identification Beamforming Channel equalization in a CDMA forward link 1
3 Interior Point Optimization for Optimal Linear Filtering A discrete-time linear system can be described by yi = x T i w + vi, i = 1, 2,... Using input output pairs {xi, yi} the linear least-squares problem is then to estimate a filter w that minimizes the mean-squared error Fn(w) = 1 n n ( ) 2 yi x T i w = 1 n yt n y n 2w T n p xy(n) + w T Rxx(n)w, (1) i=1 where yn = [y1, y2,..., yn] T, pxy(n) = 1 n n i=1 x iyi, Rxx(n) = 1 n n i=1 x ix T i. Note: pxy(n) and Rxx(n) are both recursively updatable with per-sample complexity of O(M 2 ). The optimum linear filter then satisfies Fn(w) = 0, or Rxx(n)w pxy(n) = 0. 2
4 One Motivation: Transient Convergence The RLS algorithm estimates (see e.g., Sayed and Kailath 96) w rls n := [ δ n I + R xx(n)] 1 pxy(n), where δ n I is a regularization term to improve conditioning. Problem: The regularization depends entirely on the constant δ. If, for example, SNR is underestimated = slower asymptotic convergence SNR is overestimated = bad transient behaviour Remedy: wn := [αni + Rxx(n)] 1 pxy(n), αn adjusted adaptively 3
5 The Analytic Center Approach Formulate a convex feasibility problem at each iteration. w is a feasible filter only if it is contained in Ωn = {w R M Fn(w) τn, w 2 R 2 }, (2) 1 st constraint: minimize the mean-squared error Fn(w). 2 nd constraint: make Ωn a bounded region. The analytic center w a n of Ω n is the minimizer of φn(w) = log(τn Fn(w)) log(r 2 w 2 ) which can be found by solving φn(w) = 0 ( Fn sn(w n a ) + 2wa n s tn(w n a ) = 0, therefore wa n = n(w a n ) ) 1 tn(w n a )I + R xx(n) pxy(n). where sn(w) := τn Fn(w) and tn(w) := R 2 w a n 2. 4
6 Definition of τn We are concerned with the behaviour of αn = sn/tn. 1. The goal is to make αn Fn(w). Fn(w) = 2(Rxx(n)w pxy(n)) Thus, if Fn(w) is large = need αn large for regularization. Fn(w) is small = need only a small αn. 2. Define sn := β R 2 Fn(w a n 1 ) R/ 2: a normalization constant required to show asymptotic convergence β: also required to show convergence 3. The definition of τn follows simply from τn = Fn(w a n 1 ) + s n 5
7 6 w(1) w(1) w(1) w(2) 0 w(2) 0 w(2) w(1) w(1) w(1) w(2) 0 w(2) 0 w(2) Stanford University EE392o Z.Q. Luo An Example
8 Asymptotic Convergence Analysis Condition 1. (Bounded Autocorrelation matrix) There exist n0 > 0, σ1 > 0, σ2 > 0 such that σ1i 1 n xix T i σ2i, n n0. n i=1 Condition 2. (Bounded Outputs) There exists a fixed ρy such that for all n > n0 there holds 1 n y 2 i n ρ y. i=1 The left inequality in Condition 1 is known as weak persistent excitation condition. Theorem 1. Let the sequence of estimates {wn, n = 1, 2, 3,...} be generated by the IPLS algorithm. Then Fn(wn) = 2 Rxx(n)wn pxy(n) = O (1/n). 7
9 Transient Convergence Analysis Condition 1. (Bounded Autocorrelation matrix) There exist n0 > 0, σ1 > 0, σ2 > 0 such that σ1i 1 n xix T i σ2i, n n0. n i=1 Condition 2. (No Noise) The system is free from measurement noise, i.e., yi = x T i w, i = 1, 2,.... We assume that the data has no statistical fluctuations. Convergence then implies the phasing out of effects of initialization and thus is dictated entirely by the transient behaviour of the algorithm. Theorem 2. Let the sequence of estimates {wn, n M} be generated by the IPLS algorithm. If the observations are free of noise, then wn w = O(R 1 ) wn 1 w. 8
10 Transient Analysis a simple example Example. yi = w xi + vi, with xi = w = 1, i = 1, 2,... RLS Assuming no statistical averaging (e.g., no noise), Rxx(n) = E(xx T ) = 1, pxy(n) = E(xy) = 1, n The RLS estimator then reduces to w rls n = (1 + δ/n) 1. IPLS Now, Fn(w) = (w 1) 2, and Fn(w) = 2(w 1). Evaluating τn, the condition φn(w a n ) = 0 for the analytic center becomes, 2(w a n 1) a βr 2 w n 1 a 1 + (wa n 1 1)2 (w n a + 2w n 1)2 R 2 (w = 0 n a )2 which implies w a n 1 = O(R 1 ) w a n 1 1, i.e., exponential decay of the transient error. 9
11 Direct Comparison of RLS, IPLS Property RLS IPLS Asymptotic Convergence O(1/n) O(1/n) Computational Complexity O(M 2 ) O(M 2.2 ) Transient Convergence O(1/n) O(R n ) Robustness to Initialization no yes Additional constraints req. new algorithm easily accommodated Numerical Stability (in constrained case) problems occur when λ is small (λ : forgetting factor) stable even at small values of λ, and in cases of limited precision calculations Sliding window implementation can be accommodated easily accommodated 10
12 The Interior Point Least Squares (IPLS) algorithm 1. We don t need the exact analytic center of Ωn, an approximate center is sufficient. 2. Such an approximate center is found by taking just a single Newton iteration in the minimization of φn(w). wn := wn 1 ( 2 φn(wn 1)) 1 φn(wn 1), (3) To compute (3) we need φn(w) = F n sn(w) + 2w tn(w), (4) 2 φn(w) = ( F n) T Fn s 2 n (w) + 2 Fn sn(w) + 4wwT 2I + t 2 n (w) tn(w) (5) where Fn = 2pxy(n) + 2Rxx(n)w and 2 Fn = 2Rxx(n). 3. To compute the Newton direction, an O(M 2.2 ) recursive update procedure has been devised (using the work of Powell, 1997). 11
13 Interior Point Least Squares (IPLS) Algorithm Step 1: Initialization. Let β, R be given. Set w0 = 0, pxy(0) = 0, Rxx(0) = 0, F0(0) = 0. Step 2: Updating. For n 1, acquire new data xn, yn. Then recursively update pxy(n) = n 1 n p xy(n 1) + 1 n x nyn, Rxx(n) = n 1 n R xx(n 1) + 1 n x nx T n. Update 2 Fn(wn 1) and Fn(wn 1) sn(wn 1) and tn(wn 1) ( 2 φn(wn 1)) 1 φn(wn 1) using the update procedure (or using (4), (5)) Step 3: Recentering. The new center of Ωn is obtained by taking just one Newton iteration starting at wn 1: wn := wn 1 ( 2 φn(wn 1)) 1 φn(wn 1) Set n := n + 1, and return to Step 2. 12
14 Summary Contributions provided a new look at (recursive) adaptive filtering first application of interior point optimization to a dynamic problem Features of IPLS converges asymptotically at the rate O(1/n) exhibits fast transient convergence, and is robust to initialization easily accommodates additional linear or convex quadratic constraints, and is numerically stable O(M 2.2 ) complexity 13
15 Application: System Identification Performance Measure εip(n) = wn w 2 and εrls(n) = w rls n w 2 Sources (i) White Gaussian noise, (ii) White Gaussian noise filtered through H(z) = 1 + 2z 1 + 3z 2 ( z z 2 )( z 1 ). SNRs (i) SNR1 = 40dB, (ii) SNR2 = 10dB Nominal Parameter Settings RLS λ = 1, δ = 10 4 IPLS β = 2, R = 1000 Experiment 1 w R 20, w(i) [ 1, +1], 500 independent Monte Carlo trials Experiment 2 Comparing sliding window versions of RLS (Liu & He 95) and IPLS: w R 10, Tl = 15 14
16 System Identification: Experiment 1 IPAF RLS IPAF RLS Iteration k Iteration k (a) SNR = 40dB, Source: White Gaussian Noise (b) SNR = 40dB, Source: Correlated IPAF RLS IPAF RLS Iteration k Iteration k (c) SNR = 10dB, Source: White Gaussian Noise (d) SNR = 10dB, Source: Correlated 15
17 System Identification: Experiment SW RLS SW IPAF Iteration k Figure 1: Comparison of sliding-window versions of IPLS and RLS when channel characteristics change abruptly (at iteration 100). 16
18 Application: Minimum Variance Beamforming Sensor Array Adjustable Weights Target Signal Interferer Σ Output y(n) Adaptive Filtering Algorithm Steering Vectors 17
19 Minimum Variance Beamforming By adaptively adjusting the tap weights hi(n) the beamformer must 1. Steering Capability: protect the target signal c H (θ)h(n) = 1, n, θ = θ1, θ2,... c H (θ) = [1, e jθ,..., e j(m 1)θ ], where M : number of tap weights hi, θi : Electrical Angle determined by the direction of the target i with respect to the first sensor 2. Minimize the effects of the interferers i.e., minimize the Output Power E( y 2 ) of the beamformer This beamforming problem can be cast in the framework of a Constrained Adaptive Estimation Problem. 18
20 Beamforming: Constrained Adaptive Estimation minimize Fn := 1 n n λ n i d(i) x T i h 2, i=1 subject to C T h = f, h R M, (6) λ: forgetting factor d( ): desired response x T i : vector input sequence h: vector of tap weights C, f: define the linear constraints on h In the Minimum Variance Beamforming problem the reference signal d( ) is zero. During the adaption process we assume that no target is present. The rows of C correspond to steering vector constraints. 19
21 Beamforming: Numerical Simulation Input: (interference at 0.3, and 0.7) x(n) = sin(0.3nπ) + sin(0.325nπ) + sin(0.7nπ) + b(n). b(n): white Gaussian noise at 40dB. Constraints: (desired response at freq. 0.2 and 0.5) C T = 1 cos(0.2π)... cos((m 1)0.2π) 1 cos(0.5π)... cos((m 1)0.5π) 0 sin(0.2π)... sin((m 1)0.2π) 0 sin(0.5π)... sin((m 1)0.5π) 1, f =
22 21 LCMV µ = 0.1 LCFLS λ = 0.99, Eo = 0.1 IPLS λ = 0.99, ε = 0.01, R = 100, β = 2 Precision 4 digits for LCFLS and IPLS Figure 2: (a) Freq. response at Iteration 4000, (b) Mean-squared error in h(n) Normalized frequency Time (Iteration) LCMV LC FLS IPM1 Optimum Filter 10 1 abs H(f) [db] Norm[ h opt h(n) ] LCMV LC FLS IPM Stanford University EE392o Z.Q. Luo Beamforming: Numerical Simulation
23 Application: Channel Equalization in a CDMA Downlink u (k) 1 Code Filters L C (z) 1 u (k) 2 L C (z) 2 Σ Channel noise H (z) Σ chip rate sampler Received Signal y(k) u (k) K L C (z) K Figure 3: Discrete-time model of CDMA downlink 22
24 Equalizer/Decoder Structure code matched filter Received signal y(k) FFF -1 C 1 (z ) r(k) x(k) L Σ ^ u 1 (k) x(k) FBF Figure 4: Code Matched Filter Chip rate DFE 23
25 24 Static Channel Fading Channel sampled at a random instant. Symbol Index LOS Path First multi path Second Multi Path Individual path amplitudes Fading Channel LOS component is 5 db higher than 2 multipath components, fading rate fd = 0.005, delay spread: 6Tc Sources QPSK with uniform probabilities for each symbol (i.i.d.) Stanford University EE392o Z.Q. Luo CDMA Downlink: System Description
26 25 Message Signal N = 200, NT = 10, CL = 16, Users = 1 Algorithms λ = 1.0, δ = R = 10 4, β = 2 Equalizer Mff = 14, Mfb = 2, delay = 1, pfr = SNR SNR SQRT RLS IPLS BER Prob(packet arrival) SQRT RLS IPLS Stanford University EE392o Z.Q. Luo Experiment 1: Static Channel, Single user
27 26 Message Signal N = 200, NT = 10, CL = 16, Users = 4 Algorithms λ = 1.0, δ = R = 10 4, β = 2 Equalizer Mff = 14, Mfb = 2, delay = 1, pfr = SNR SNR SQRT RLS IPLS BER 10 2 Prob(packet arrival) SQRT RLS IPLS Stanford University EE392o Z.Q. Luo Experiment 1: Static Channel, 4 Users
28 27 Message Signal N = 200, SNR = 12dB, CL = 16, Users = 1 Algorithms λ = 1.0, δ = R = 10 4, β = 2 Equalizer Mff = 14, Mfb = 2, delay = 1, pfr = N T N T 0.1 SQRT RLS IPLS BER Prob(packet arrival) SQRT RLS IPLS Stanford University EE392o Z.Q. Luo Experiment 1: Dependence on Training Length
29 28 Message Signal N = 200, NT = 2/10, CL = 16, Users = 1 Algorithms λ = 0.85, δ = R = 10 4, β = 2 Equalizer Mff = 14, Mfb = 2, delay = 1, pfr = 10 2 SNR SNR SQRT RLS, N =2 T SQRT RLS, N T =10 IPLS, N =2 T IPLS, N =10 T BER 10 2 Prob(packet arrival) 0.5 SQRT RLS, N T =2 SQRT RLS, N T =10 IPLS, N T =2 IPLS, N T = Stanford University EE392o Z.Q. Luo Experiment 2: Time-Varying Channel, Single user
30 29 Message Signal N = 200, NT = 2/10, CL = 16, Users = 4 Algorithms λ = 0.85, δ = R = 10 4, β = 2 Equalizer Mff = 14, Mfb = 2, delay = 1, pfr = 10 2 SNR SNR SQRT RLS, N T =10 IPLS, N =10 T BER Prob(packet arrival) SQRT RLS, N =10 T IPLS, N =10 T Stanford University EE392o Z.Q. Luo Experiment 2: Time-Varying Channel, 4 users
31 30 Message Signal N = 200, SNR = 16dB, CL = 16, Users = 1 Algorithms λ = 0.85, δ = R = 10 4, β = 2 Equalizer Mff = 14, Mfb = 2, delay = 1, pfr = N T N T 0.2 SQRT RLS IPLS BER Prob(packet arrival) SQRT RLS IPLS Stanford University EE392o Z.Q. Luo Experiment 2: Dependence on Training Length
32 Conclusions Transient convergence of IPLS is O(1/R n ), n M Superior transient convergence to RLS even when n < M Gain of using IPLS over the RLS algorithm can range from 5-6 db to well over 10 db 31
Lecture 7: Wireless Channels and Diversity Advanced Digital Communications (EQ2410) 1
Wireless : Wireless Advanced Digital Communications (EQ2410) 1 Thursday, Feb. 11, 2016 10:00-12:00, B24 1 Textbook: U. Madhow, Fundamentals of Digital Communications, 2008 1 / 15 Wireless Lecture 1-6 Equalization
More informationAdaptive Filter Theory
0 Adaptive Filter heory Sung Ho Cho Hanyang University Seoul, Korea (Office) +8--0-0390 (Mobile) +8-10-541-5178 dragon@hanyang.ac.kr able of Contents 1 Wiener Filters Gradient Search by Steepest Descent
More informationBeamforming. A brief introduction. Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University
Beamforming A brief introduction Brian D. Jeffs Associate Professor Dept. of Electrical and Computer Engineering Brigham Young University March 2008 References Barry D. Van Veen and Kevin Buckley, Beamforming:
More informationDecision Weighted Adaptive Algorithms with Applications to Wireless Channel Estimation
Decision Weighted Adaptive Algorithms with Applications to Wireless Channel Estimation Shane Martin Haas April 12, 1999 Thesis Defense for the Degree of Master of Science in Electrical Engineering Department
More informationJoint Channel Estimation and Co-Channel Interference Mitigation in Wireless Networks Using Belief Propagation
Joint Channel Estimation and Co-Channel Interference Mitigation in Wireless Networks Using Belief Propagation Yan Zhu, Dongning Guo and Michael L. Honig Northwestern University May. 21, 2008 Y. Zhu, D.
More informationLecture 2. Fading Channel
1 Lecture 2. Fading Channel Characteristics of Fading Channels Modeling of Fading Channels Discrete-time Input/Output Model 2 Radio Propagation in Free Space Speed: c = 299,792,458 m/s Isotropic Received
More informationTensor approach for blind FIR channel identification using 4th-order cumulants
Tensor approach for blind FIR channel identification using 4th-order cumulants Carlos Estêvão R Fernandes Gérard Favier and João Cesar M Mota contact: cfernand@i3s.unice.fr June 8, 2006 Outline 1. HOS
More informationRobust Capon Beamforming
Robust Capon Beamforming Yi Jiang Petre Stoica Zhisong Wang Jian Li University of Florida Uppsala University University of Florida University of Florida March 11, 2003 ASAP Workshop 2003 1 Outline Standard
More informationDiversity Combining Techniques
Diversity Combining Techniques When the required signal is a combination of several plane waves (multipath), the total signal amplitude may experience deep fades (Rayleigh fading), over time or space.
More informationLinear Optimum Filtering: Statement
Ch2: Wiener Filters Optimal filters for stationary stochastic models are reviewed and derived in this presentation. Contents: Linear optimal filtering Principle of orthogonality Minimum mean squared error
More informationAdaptive Filtering Part II
Adaptive Filtering Part II In previous Lecture we saw that: Setting the gradient of cost function equal to zero, we obtain the optimum values of filter coefficients: (Wiener-Hopf equation) Adaptive Filtering,
More informationTransmitting. Base-station. Interfering. Base-station Mobile Receiver of Desired User
ReducedRank Linear Equalization for the Forward Link in HighSpeed CDMA PhD Preliminary Examination, Electrical and Computer Engineering, Purdue University. Major Professor : Michael D. Zoltowski November
More informationOverview of Beamforming
Overview of Beamforming Arye Nehorai Preston M. Green Department of Electrical and Systems Engineering Washington University in St. Louis March 14, 2012 CSSIP Lab 1 Outline Introduction Spatial and temporal
More informationBlind Source Separation with a Time-Varying Mixing Matrix
Blind Source Separation with a Time-Varying Mixing Matrix Marcus R DeYoung and Brian L Evans Department of Electrical and Computer Engineering The University of Texas at Austin 1 University Station, Austin,
More informationNOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION. M. Schwab, P. Noll, and T. Sikora. Technical University Berlin, Germany Communication System Group
NOISE ROBUST RELATIVE TRANSFER FUNCTION ESTIMATION M. Schwab, P. Noll, and T. Sikora Technical University Berlin, Germany Communication System Group Einsteinufer 17, 1557 Berlin (Germany) {schwab noll
More informationBlind Channel Equalization in Impulse Noise
Blind Channel Equalization in Impulse Noise Rubaiyat Yasmin and Tetsuya Shimamura Graduate School of Science and Engineering, Saitama University 255 Shimo-okubo, Sakura-ku, Saitama 338-8570, Japan yasmin@sie.ics.saitama-u.ac.jp
More informationEE 5407 Part II: Spatial Based Wireless Communications
EE 5407 Part II: Spatial Based Wireless Communications Instructor: Prof. Rui Zhang E-mail: rzhang@i2r.a-star.edu.sg Website: http://www.ece.nus.edu.sg/stfpage/elezhang/ Lecture II: Receive Beamforming
More informationParametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d)
Parametric Signal Modeling and Linear Prediction Theory 1. Discrete-time Stochastic Processes (cont d) Electrical & Computer Engineering North Carolina State University Acknowledgment: ECE792-41 slides
More informationSquare Root Raised Cosine Filter
Wireless Information Transmission System Lab. Square Root Raised Cosine Filter Institute of Communications Engineering National Sun Yat-sen University Introduction We consider the problem of signal design
More informationOptimization in Wireless Communication
Zhi-Quan (Tom) Luo Department of Electrical and Computer Engineering University of Minnesota 200 Union Street SE Minneapolis, MN 55455 2007 NSF Workshop Challenges Optimization problems from wireless applications
More informationPerformance Analysis of Spread Spectrum CDMA systems
1 Performance Analysis of Spread Spectrum CDMA systems 16:33:546 Wireless Communication Technologies Spring 5 Instructor: Dr. Narayan Mandayam Summary by Liang Xiao lxiao@winlab.rutgers.edu WINLAB, Department
More informationSliding Window Recursive Quadratic Optimization with Variable Regularization
11 American Control Conference on O'Farrell Street, San Francisco, CA, USA June 29 - July 1, 11 Sliding Window Recursive Quadratic Optimization with Variable Regularization Jesse B. Hoagg, Asad A. Ali,
More informationExpectation propagation for signal detection in flat-fading channels
Expectation propagation for signal detection in flat-fading channels Yuan Qi MIT Media Lab Cambridge, MA, 02139 USA yuanqi@media.mit.edu Thomas Minka CMU Statistics Department Pittsburgh, PA 15213 USA
More informationLeast Squares Estimation Namrata Vaswani,
Least Squares Estimation Namrata Vaswani, namrata@iastate.edu Least Squares Estimation 1 Recall: Geometric Intuition for Least Squares Minimize J(x) = y Hx 2 Solution satisfies: H T H ˆx = H T y, i.e.
More information!P x. !E x. Bad Things Happen to Good Signals Spring 2011 Lecture #6. Signal-to-Noise Ratio (SNR) Definition of Mean, Power, Energy ( ) 2.
Bad Things Happen to Good Signals oise, broadly construed, is any change to the signal from its expected value, x[n] h[n], when it arrives at the receiver. We ll look at additive noise and assume the noise
More informationAdvanced 3 G and 4 G Wireless Communication Prof. Aditya K Jagannathan Department of Electrical Engineering Indian Institute of Technology, Kanpur
Advanced 3 G and 4 G Wireless Communication Prof. Aditya K Jagannathan Department of Electrical Engineering Indian Institute of Technology, Kanpur Lecture - 19 Multi-User CDMA Uplink and Asynchronous CDMA
More informationECE531 Lecture 10b: Dynamic Parameter Estimation: System Model
ECE531 Lecture 10b: Dynamic Parameter Estimation: System Model D. Richard Brown III Worcester Polytechnic Institute 02-Apr-2009 Worcester Polytechnic Institute D. Richard Brown III 02-Apr-2009 1 / 14 Introduction
More informationAn Adaptive Sensor Array Using an Affine Combination of Two Filters
An Adaptive Sensor Array Using an Affine Combination of Two Filters Tõnu Trump Tallinn University of Technology Department of Radio and Telecommunication Engineering Ehitajate tee 5, 19086 Tallinn Estonia
More informationDirect-Sequence Spread-Spectrum
Chapter 3 Direct-Sequence Spread-Spectrum In this chapter we consider direct-sequence spread-spectrum systems. Unlike frequency-hopping, a direct-sequence signal occupies the entire bandwidth continuously.
More information5 Kalman filters. 5.1 Scalar Kalman filter. Unit delay Signal model. System model
5 Kalman filters 5.1 Scalar Kalman filter 5.1.1 Signal model System model {Y (n)} is an unobservable sequence which is described by the following state or system equation: Y (n) = h(n)y (n 1) + Z(n), n
More informationLECTURE 18. Lecture outline Gaussian channels: parallel colored noise inter-symbol interference general case: multiple inputs and outputs
LECTURE 18 Last time: White Gaussian noise Bandlimited WGN Additive White Gaussian Noise (AWGN) channel Capacity of AWGN channel Application: DS-CDMA systems Spreading Coding theorem Lecture outline Gaussian
More informationFundamentals of Statistical Signal Processing Volume II Detection Theory
Fundamentals of Statistical Signal Processing Volume II Detection Theory Steven M. Kay University of Rhode Island PH PTR Prentice Hall PTR Upper Saddle River, New Jersey 07458 http://www.phptr.com Contents
More informationEstimation of the Capacity of Multipath Infrared Channels
Estimation of the Capacity of Multipath Infrared Channels Jeffrey B. Carruthers Department of Electrical and Computer Engineering Boston University jbc@bu.edu Sachin Padma Department of Electrical and
More informationAdaptive Filters. un [ ] yn [ ] w. yn n wun k. - Adaptive filter (FIR): yn n n w nun k. (1) Identification. Unknown System + (2) Inverse modeling
Adaptive Filters - Statistical digital signal processing: in many problems of interest, the signals exhibit some inherent variability plus additive noise we use probabilistic laws to model the statistical
More informationComputer Intensive Methods in Mathematical Statistics
Computer Intensive Methods in Mathematical Statistics Department of mathematics johawes@kth.se Lecture 16 Advanced topics in computational statistics 18 May 2017 Computer Intensive Methods (1) Plan of
More informationChapter 2 Fundamentals of Adaptive Filter Theory
Chapter 2 Fundamentals of Adaptive Filter Theory In this chapter we will treat some fundamentals of the adaptive filtering theory highlighting the system identification problem We will introduce a signal
More informationJ. Liang School of Automation & Information Engineering Xi an University of Technology, China
Progress In Electromagnetics Research C, Vol. 18, 245 255, 211 A NOVEL DIAGONAL LOADING METHOD FOR ROBUST ADAPTIVE BEAMFORMING W. Wang and R. Wu Tianjin Key Lab for Advanced Signal Processing Civil Aviation
More informationImpulsive Noise Filtering In Biomedical Signals With Application of New Myriad Filter
BIOSIGAL 21 Impulsive oise Filtering In Biomedical Signals With Application of ew Myriad Filter Tomasz Pander 1 1 Division of Biomedical Electronics, Institute of Electronics, Silesian University of Technology,
More informationImproved channel estimation for massive MIMO systems using hybrid pilots with pilot anchoring
Improved channel estimation for massive MIMO systems using hybrid pilots with pilot anchoring Karthik Upadhya, Sergiy A. Vorobyov, Mikko Vehkapera Department of Signal Processing and Acoustics Aalto University,
More informationPassive Sonar Detection Performance Prediction of a Moving Source in an Uncertain Environment
Acoustical Society of America Meeting Fall 2005 Passive Sonar Detection Performance Prediction of a Moving Source in an Uncertain Environment Vivek Varadarajan and Jeffrey Krolik Duke University Department
More informationELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications
ELE539A: Optimization of Communication Systems Lecture 15: Semidefinite Programming, Detection and Estimation Applications Professor M. Chiang Electrical Engineering Department, Princeton University March
More informationMultiple Carrier Frequency Offset and Channel State Estimation in the Fading Channel
Multiple Carrier Frequency Offset and Channel State Estimation in the Fading Channel Brad W Zarikoff, James K Cavers School of Engineering Science, Faculty of Applied Sciences Simon Fraser University Burnaby,
More informationECS455: Chapter 5 OFDM. ECS455: Chapter 5 OFDM. OFDM: Overview. OFDM Applications. Dr.Prapun Suksompong prapun.com/ecs455
ECS455: Chapter 5 OFDM OFDM: Overview Let S = (S 1, S 2,, S ) contains the information symbols. S IFFT FFT Inverse fast Fourier transform Fast Fourier transform 1 Dr.Prapun Suksompong prapun.com/ecs455
More informationDesign of MMSE Multiuser Detectors using Random Matrix Techniques
Design of MMSE Multiuser Detectors using Random Matrix Techniques Linbo Li and Antonia M Tulino and Sergio Verdú Department of Electrical Engineering Princeton University Princeton, New Jersey 08544 Email:
More informationEE538 Final Exam Fall :20 pm -5:20 pm PHYS 223 Dec. 17, Cover Sheet
EE538 Final Exam Fall 005 3:0 pm -5:0 pm PHYS 3 Dec. 17, 005 Cover Sheet Test Duration: 10 minutes. Open Book but Closed Notes. Calculators ARE allowed!! This test contains five problems. Each of the five
More informationZERO-FORCING MULTIUSER DETECTION IN CDMA SYSTEMS USING LONG CODES
ZERO-FORCING MUTIUSER DETECTION IN CDMA SYSTEMS USING ONG CODES Cássio B Ribeiro, Marcello R de Campos, and Paulo S R Diniz Electrical Engineering Program and Department of Electronics and Computer Engineering
More informationSemi-Definite Programming (SDP) Relaxation Based Semi-Blind Channel Estimation for Frequency-Selective MIMO MC-CDMA Systems
Semi-Definite Programming (SDP) Relaxation Based Semi-Blind Channel Estimation for Frequency-Selective MIMO MC-CDMA Systems Naveen K. D. Venkategowda Department of Electrical Engineering Indian Institute
More informationAdaptive Channel Modeling for MIMO Wireless Communications
Adaptive Channel Modeling for MIMO Wireless Communications Chengjin Zhang Department of Electrical and Computer Engineering University of California, San Diego San Diego, CA 99- Email: zhangc@ucsdedu Robert
More informationLecture Notes 4 Vector Detection and Estimation. Vector Detection Reconstruction Problem Detection for Vector AGN Channel
Lecture Notes 4 Vector Detection and Estimation Vector Detection Reconstruction Problem Detection for Vector AGN Channel Vector Linear Estimation Linear Innovation Sequence Kalman Filter EE 278B: Random
More informationA ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT
Progress In Electromagnetics Research Letters, Vol. 16, 53 60, 2010 A ROBUST BEAMFORMER BASED ON WEIGHTED SPARSE CONSTRAINT Y. P. Liu and Q. Wan School of Electronic Engineering University of Electronic
More informationADAPTIVE DETECTION FOR A PERMUTATION-BASED MULTIPLE-ACCESS SYSTEM ON TIME-VARYING MULTIPATH CHANNELS WITH UNKNOWN DELAYS AND COEFFICIENTS
ADAPTIVE DETECTION FOR A PERMUTATION-BASED MULTIPLE-ACCESS SYSTEM ON TIME-VARYING MULTIPATH CHANNELS WITH UNKNOWN DELAYS AND COEFFICIENTS Martial COULON and Daniel ROVIRAS University of Toulouse INP-ENSEEIHT
More informationPenalty and Barrier Methods General classical constrained minimization problem minimize f(x) subject to g(x) 0 h(x) =0 Penalty methods are motivated by the desire to use unconstrained optimization techniques
More informationLecture 19 IIR Filters
Lecture 19 IIR Filters Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/5/10 1 General IIR Difference Equation IIR system: infinite-impulse response system The most general class
More informationADAPTIVE ANTENNAS. SPATIAL BF
ADAPTIVE ANTENNAS SPATIAL BF 1 1-Spatial reference BF -Spatial reference beamforming may not use of embedded training sequences. Instead, the directions of arrival (DoA) of the impinging waves are used
More informationChapter 10. Timing Recovery. March 12, 2008
Chapter 10 Timing Recovery March 12, 2008 b[n] coder bit/ symbol transmit filter, pt(t) Modulator Channel, c(t) noise interference form other users LNA/ AGC Demodulator receive/matched filter, p R(t) sampler
More informationRobust Range-rate Estimation of Passive Narrowband Sources in Shallow Water
Robust Range-rate Estimation of Passive Narrowband Sources in Shallow Water p. 1/23 Robust Range-rate Estimation of Passive Narrowband Sources in Shallow Water Hailiang Tao and Jeffrey Krolik Department
More informationEfficient Semi-Blind Channel Estimation and Equalization Based on a Parametric Channel Representation
Efficient Semi-Blind Channel Estimation and Equalization Based on a Parametric Channel Representation Presenter: Kostas Berberidis University of Patras Computer Engineering & Informatics Department Signal
More informationCh. 7: Z-transform Reading
c J. Fessler, June 9, 3, 6:3 (student version) 7. Ch. 7: Z-transform Definition Properties linearity / superposition time shift convolution: y[n] =h[n] x[n] Y (z) =H(z) X(z) Inverse z-transform by coefficient
More informationSparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels
Sparse Least Mean Square Algorithm for Estimation of Truncated Volterra Kernels Bijit Kumar Das 1, Mrityunjoy Chakraborty 2 Department of Electronics and Electrical Communication Engineering Indian Institute
More informationSignal Design for Band-Limited Channels
Wireless Information Transmission System Lab. Signal Design for Band-Limited Channels Institute of Communications Engineering National Sun Yat-sen University Introduction We consider the problem of signal
More informationLecture 6: Modeling of MIMO Channels Theoretical Foundations of Wireless Communications 1
Fading : Theoretical Foundations of Wireless Communications 1 Thursday, May 3, 2018 9:30-12:00, Conference Room SIP 1 Textbook: D. Tse and P. Viswanath, Fundamentals of Wireless Communication 1 / 23 Overview
More informationState-Space Methods for Inferring Spike Trains from Calcium Imaging
State-Space Methods for Inferring Spike Trains from Calcium Imaging Joshua Vogelstein Johns Hopkins April 23, 2009 Joshua Vogelstein (Johns Hopkins) State-Space Calcium Imaging April 23, 2009 1 / 78 Outline
More informationLecture 6: Modeling of MIMO Channels Theoretical Foundations of Wireless Communications 1. Overview. CommTh/EES/KTH
: Theoretical Foundations of Wireless Communications 1 Wednesday, May 11, 2016 9:00-12:00, Conference Room SIP 1 Textbook: D. Tse and P. Viswanath, Fundamentals of Wireless Communication 1 / 1 Overview
More informationOptimal and Adaptive Filtering
Optimal and Adaptive Filtering Murat Üney M.Uney@ed.ac.uk Institute for Digital Communications (IDCOM) 26/06/2017 Murat Üney (IDCOM) Optimal and Adaptive Filtering 26/06/2017 1 / 69 Table of Contents 1
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 11 Adaptive Filtering 14/03/04 http://www.ee.unlv.edu/~b1morris/ee482/
More informationELEC546 Review of Information Theory
ELEC546 Review of Information Theory Vincent Lau 1/1/004 1 Review of Information Theory Entropy: Measure of uncertainty of a random variable X. The entropy of X, H(X), is given by: If X is a discrete random
More informationMedian Filter Based Realizations of the Robust Time-Frequency Distributions
TIME-FREQUENCY SIGNAL ANALYSIS 547 Median Filter Based Realizations of the Robust Time-Frequency Distributions Igor Djurović, Vladimir Katkovnik, LJubiša Stanković Abstract Recently, somenewefficient tools
More informationEXTENSION OF NESTED ARRAYS WITH THE FOURTH-ORDER DIFFERENCE CO-ARRAY ENHANCEMENT
EXTENSION OF NESTED ARRAYS WITH THE FOURTH-ORDER DIFFERENCE CO-ARRAY ENHANCEMENT Qing Shen,2, Wei Liu 2, Wei Cui, Siliang Wu School of Information and Electronics, Beijing Institute of Technology Beijing,
More informationMinimum BER Linear Transceivers for Block. Communication Systems. Lecturer: Tom Luo
Minimum BER Linear Transceivers for Block Communication Systems Lecturer: Tom Luo Outline Block-by-block communication Abstract model Applications Current design techniques Minimum BER precoders for zero-forcing
More informationMulti User Detection I
January 12, 2005 Outline Overview Multiple Access Communication Motivation: What is MU Detection? Overview of DS/CDMA systems Concept and Codes used in CDMA CDMA Channels Models Synchronous and Asynchronous
More informationDynamic Multipath Estimation by Sequential Monte Carlo Methods
Dynamic Multipath Estimation by Sequential Monte Carlo Methods M. Lentmaier, B. Krach, P. Robertson, and T. Thiasiriphet German Aerospace Center (DLR) Slide 1 Outline Multipath problem and signal model
More informationLecture 4: Linear and quadratic problems
Lecture 4: Linear and quadratic problems linear programming examples and applications linear fractional programming quadratic optimization problems (quadratically constrained) quadratic programming second-order
More informationEffect of Multipath Propagation on the Noise Performance of Zero Crossing Digital Phase Locked Loop
International Journal of Electronics and Communication Engineering. ISSN 0974-2166 Volume 9, Number 2 (2016), pp. 97-104 International Research Publication House http://www.irphouse.com Effect of Multipath
More informationADAPTIVE FILTER ALGORITHMS. Prepared by Deepa.T, Asst.Prof. /TCE
ADAPTIVE FILTER ALGORITHMS Prepared by Deepa.T, Asst.Prof. /TCE Equalization Techniques Fig.3 Classification of equalizers Equalizer Techniques Linear transversal equalizer (LTE, made up of tapped delay
More informationDigital Band-pass Modulation PROF. MICHAEL TSAI 2011/11/10
Digital Band-pass Modulation PROF. MICHAEL TSAI 211/11/1 Band-pass Signal Representation a t g t General form: 2πf c t + φ t g t = a t cos 2πf c t + φ t Envelope Phase Envelope is always non-negative,
More informationTheoretical Development for Blind Identification of Non Linear Communication Channels
, pp.1-8 http://dx.doi.org/10.14257/ijeic.2015.6.6.01 Theoretical Development for Blind Identification of Non Linear Communication Channels Mohammed Zidane 1, Said Safi 2, Mohamed Sabri 1 and Ahmed Boumezzough
More informationDiscrete-Time Signals & Systems
Chapter 2 Discrete-Time Signals & Systems 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 2-1-1 Discrete-Time Signals: Time-Domain Representation (1/10) Signals
More informationExamples. 2-input, 1-output discrete-time systems: 1-input, 1-output discrete-time systems:
Discrete-Time s - I Time-Domain Representation CHAPTER 4 These lecture slides are based on "Digital Signal Processing: A Computer-Based Approach, 4th ed." textbook by S.K. Mitra and its instructor materials.
More informationLecture 4 Capacity of Wireless Channels
Lecture 4 Capacity of Wireless Channels I-Hsiang Wang ihwang@ntu.edu.tw 3/0, 014 What we have learned So far: looked at specific schemes and techniques Lecture : point-to-point wireless channel - Diversity:
More informationDigital Communications: A Discrete-Time Approach M. Rice. Errata. Page xiii, first paragraph, bare witness should be bear witness
Digital Communications: A Discrete-Time Approach M. Rice Errata Foreword Page xiii, first paragraph, bare witness should be bear witness Page xxi, last paragraph, You know who you. should be You know who
More informationROBUST ADAPTIVE BEAMFORMING BASED ON CO- VARIANCE MATRIX RECONSTRUCTION FOR LOOK DIRECTION MISMATCH
Progress In Electromagnetics Research Letters, Vol. 25, 37 46, 2011 ROBUST ADAPTIVE BEAMFORMING BASED ON CO- VARIANCE MATRIX RECONSTRUCTION FOR LOOK DIRECTION MISMATCH R. Mallipeddi 1, J. P. Lie 2, S.
More informationRobust Adaptive Beam-forming Algorithm Using Rotation and Iteration
Robust Adaptive Beam-forming Algorithm Using Rotation and Iteration Guanling Wang, Zhengchuan Fang, Fengsui Wang School of Electrical Engineering, Anhui Polytechnic University, Wuhu, Anhui, 241000, China
More informationKNOWN approaches for improving the performance of
IEEE TRANSACTIONS ON CIRCUITS AND SYSTEMS II: EXPRESS BRIEFS, VOL. 58, NO. 8, AUGUST 2011 537 Robust Quasi-Newton Adaptive Filtering Algorithms Md. Zulfiquar Ali Bhotto, Student Member, IEEE, and Andreas
More informationBlind Identification of FIR Systems and Deconvolution of White Input Sequences
Blind Identification of FIR Systems and Deconvolution of White Input Sequences U. SOVERINI, P. CASTALDI, R. DIVERSI and R. GUIDORZI Dipartimento di Elettronica, Informatica e Sistemistica Università di
More informationAutomatic Speech Recognition (CS753)
Automatic Speech Recognition (CS753) Lecture 12: Acoustic Feature Extraction for ASR Instructor: Preethi Jyothi Feb 13, 2017 Speech Signal Analysis Generate discrete samples A frame Need to focus on short
More informationLecture 8: Policy Gradient
Lecture 8: Policy Gradient Hado van Hasselt Outline 1 Introduction 2 Finite Difference Policy Gradient 3 Monte-Carlo Policy Gradient 4 Actor-Critic Policy Gradient Introduction Vapnik s rule Never solve
More informationUNIT 1. SIGNALS AND SYSTEM
Page no: 1 UNIT 1. SIGNALS AND SYSTEM INTRODUCTION A SIGNAL is defined as any physical quantity that changes with time, distance, speed, position, pressure, temperature or some other quantity. A SIGNAL
More informationCHAPTER 14. Based on the info about the scattering function we know that the multipath spread is T m =1ms, and the Doppler spread is B d =0.2 Hz.
CHAPTER 4 Problem 4. : Based on the info about the scattering function we know that the multipath spread is T m =ms, and the Doppler spread is B d =. Hz. (a) (i) T m = 3 sec (ii) B d =. Hz (iii) ( t) c
More informationCHAPTER 3 ROBUST ADAPTIVE BEAMFORMING
50 CHAPTER 3 ROBUST ADAPTIVE BEAMFORMING 3.1 INTRODUCTION Adaptive beamforming is used for enhancing a desired signal while suppressing noise and interference at the output of an array of sensors. It is
More informationCoherentDetectionof OFDM
Telematics Lab IITK p. 1/50 CoherentDetectionof OFDM Indo-UK Advanced Technology Centre Supported by DST-EPSRC K Vasudevan Associate Professor vasu@iitk.ac.in Telematics Lab Department of EE Indian Institute
More informationAdaptive beamforming for uniform linear arrays with unknown mutual coupling. IEEE Antennas and Wireless Propagation Letters.
Title Adaptive beamforming for uniform linear arrays with unknown mutual coupling Author(s) Liao, B; Chan, SC Citation IEEE Antennas And Wireless Propagation Letters, 2012, v. 11, p. 464-467 Issued Date
More informationEECE Adaptive Control
EECE 574 - Adaptive Control Basics of System Identification Guy Dumont Department of Electrical and Computer Engineering University of British Columbia January 2010 Guy Dumont (UBC) EECE574 - Basics of
More informationRoot-MUSIC Time Delay Estimation Based on Propagator Method Bin Ba, Yun Long Wang, Na E Zheng & Han Ying Hu
International Conference on Automation, Mechanical Control and Computational Engineering (AMCCE 15) Root-MUSIC ime Delay Estimation Based on ropagator Method Bin Ba, Yun Long Wang, Na E Zheng & an Ying
More informationCh5: Least Mean-Square Adaptive Filtering
Ch5: Least Mean-Square Adaptive Filtering Introduction - approximating steepest-descent algorithm Least-mean-square algorithm Stability and performance of the LMS algorithm Robustness of the LMS algorithm
More informationMultiple Bits Distributed Moving Horizon State Estimation for Wireless Sensor Networks. Ji an Luo
Multiple Bits Distributed Moving Horizon State Estimation for Wireless Sensor Networks Ji an Luo 2008.6.6 Outline Background Problem Statement Main Results Simulation Study Conclusion Background Wireless
More informationComparative Performance Analysis of Three Algorithms for Principal Component Analysis
84 R. LANDQVIST, A. MOHAMMED, COMPARATIVE PERFORMANCE ANALYSIS OF THR ALGORITHMS Comparative Performance Analysis of Three Algorithms for Principal Component Analysis Ronnie LANDQVIST, Abbas MOHAMMED Dept.
More informationDetermining the Optimal Decision Delay Parameter for a Linear Equalizer
International Journal of Automation and Computing 1 (2005) 20-24 Determining the Optimal Decision Delay Parameter for a Linear Equalizer Eng Siong Chng School of Computer Engineering, Nanyang Technological
More informationSlide Set Data Converters. Digital Enhancement Techniques
0 Slide Set Data Converters Digital Enhancement Techniques Introduction Summary Error Measurement Trimming of Elements Foreground Calibration Background Calibration Dynamic Matching Decimation and Interpolation
More informationRADIO SYSTEMS ETIN15. Lecture no: Equalization. Ove Edfors, Department of Electrical and Information Technology
RADIO SYSTEMS ETIN15 Lecture no: 8 Equalization Ove Edfors, Department of Electrical and Information Technology Ove.Edfors@eit.lth.se Contents Inter-symbol interference Linear equalizers Decision-feedback
More informationECE534, Spring 2018: Solutions for Problem Set #5
ECE534, Spring 08: s for Problem Set #5 Mean Value and Autocorrelation Functions Consider a random process X(t) such that (i) X(t) ± (ii) The number of zero crossings, N(t), in the interval (0, t) is described
More information