Solutions for examination in TSRT78 Digital Signal Processing,
|
|
- Madison Boone
- 6 years ago
- Views:
Transcription
1 Solutions for examination in TSRT78 Digital Signal Processing, s(t) is generated by s(t) = 1 w(t), q 1 Var(w(t)) = σ 2 w = 2. It is measured as y(t) = s(t) + n(t) where n(t) is white noise with Var(n(t)) = σ 2 n = 1. The noise processes are independent. a) 1 2 ( Φ ss (z) = σ 2 1 )( 1 ) z 1 w = z 1 = z (z + z 1 ) Since w(t) and e(t) are independent and white, it holds that Φ yy (z) = Φ ss (z) + Φ ee (z) = Φ ss (z) + σ 2 n = (z + z 1 ) (z + z 1 ). Furthermore, Φ sy (z) = Φ ss (z) follows from the independence of w(t) and e(t). The noncausal Wiener filter is given by H(z) = Φ sy(z) Φ yy (z) = (z + z 1 ) = 6.667z z z + 1 = 6.667z (z 10.20)(z 0.098) A separation into a causal and a anti-causal part is given below. b) Var(s(t) ŝ(t)) = R ss (0) i h(i)r sy(i) where R sy (i) = R ss (i). We begin with the covariance function R ss (k) = E [s(t)s(t k)]. For k = 0 and k = 1 we obtain R ss (0) = E [s(t)( 0.3s(t 1) + w(t))] = 0.3R ss (1) + 2, R ss (1) = E [( 0.3s(t 1) + w(t))s(t 1)] = 0.3R ss (0). Elimination of R ss (1) in the above equations yields R ss (0) = For general k > 0 we have R ss (k) = E [s(t + k)s(t)] = E [( 0.3s(t + k 1) + w(t + k))s(t)] = 0.3R ss (k 1) =... = 0.3 k R ss (0) and R ss ( k) = R ss (k) by symmetry. From this we can conclude that R ss (k) = ( 0.3) k. After decomposing H(z) into partial fractions we obtain H(z) = z z z z ( ) 0.098z = z z ( 1 ) = ( 0.098z) k + ( 0.098z 1 ) k = k=1 1 k= where we have used the geometric series results k=0 ( 0.098) k z k ( 0.098) k z k k=0 ( r) k = r k=0 and ( r) k = r 1 + r k=1 for r < 1. From the definition of z-transform, H(z) = k= h(k)z k, we get h(k) = ( 0.098) k. The error variance is then 1 Ver: 2014/05/06
2 Var(s(t) ŝ(t)) = R ss (0) k h(k)r sy (k) = k= = = ( ( ( 0.098) k )( ( 0.3) k ) k= ( ) k ( ) ( ) = = a) The zero-padded DTFT can be plotted by the following MATLAB code: % Generate the signal N = 128; x = 0:1:(N-1); A = 0.03; y = sin( 2 * pi * 0.29 * x ) + A * sin( 2 * pi * 0.25 * x); y1 = sin( 2 * pi * 0.29 * x ); y2 = A * sin( 2 * pi * 0.25 * x); %% Task 2a % Compute the zero-padded DFT and create the frequency grid yzeropadded = [y zeros(1,7*n)]; Yzeropadded = fft(yzeropadded); yzeropadded1 = [y1 zeros(1,7*n)]; Yzeropadded1 = fft(yzeropadded1); yzeropadded2 = [y2 zeros(1,7*n)]; Yzeropadded2 = fft(yzeropadded2); grid = 2 * pi / ( 8 * N * 1) * (0:1:(8*N-1) ); % Plot the DFT figure(1) subplot(3,2,[1 2]) semilogy(grid,abs(yzeropadded), k, LineWidth,2) set(gca, xlim,[2*pi*0.2 2*pi*0.35]) set(gca, ylim,[ ]) xlabel( Frequency (rad/s) ); ylabel( DTFT of signal ); hold on; plot(2*pi*0.29*[1 1],[ ], g, LineWidth,2) plot(2*pi*0.25*[1 1],[ ], g, LineWidth,2) subplot(3,2,3) semilogy(grid,abs(yzeropadded1), k, LineWidth,2) set(gca, xlim,[2*pi*0.2 2*pi*0.35]) set(gca, ylim,[ ]) xlabel( Frequency (rad/s) ); ylabel( DTFT of signal 1 ); hold on; plot(2*pi*0.29*[1 1],[ ], g, LineWidth,2) plot(2*pi*0.25*[1 1],[ ], g, LineWidth,2) subplot(3,2,4) semilogy(grid,abs(yzeropadded2), k, LineWidth,2) set(gca, xlim,[2*pi*0.2 2*pi*0.35]) set(gca, ylim,[ ]) xlabel( Frequency (rad/s) ); ylabel( DTFT of signal 2 ); ) ( ) k k=1 2 Ver: 2014/05/06
3 hold on; plot(2*pi*0.29*[1 1],[ ], g, LineWidth,2) plot(2*pi*0.25*[1 1],[ ], g, LineWidth,2) Running this code with A = 0.05 results in Figure 1, where we see that the peak of the weaker signal barely sticks out from the side lobes of the main peak. The limit of A to give a lobe that is detectable in the DTFT is around If the value is smaller than that it is difficult to distinguish from the side lobes of the main frequency. b) The difficulty in this case is mainly due to the leakage from the main lobe in the DTFT approximation. This can be solved to some extent by adding a window to the signal before computing the DTFT. c) We apply a Hanning window to the signal, zero-pad it and compute the DTFT using the following Matlab code: %% Task 2c % Windowing w = hanning(n); yw = w(:).*y(:) ; % Zero-pad and Approximate DTFT yzeropaddedw = [yw zeros(1,7*n)]; YzeropaddedW = fft(yzeropaddedw); subplot(3,2,[5 6]) semilogy(grid,abs(yzeropaddedw), k, LineWidth,2) set(gca, xlim,[2*pi*0.2 2*pi*0.35]) set(gca, ylim,[ ]) xlabel( Frequency (rad/s) ); ylabel( DTFT of signal ); hold on; plot(2*pi*0.29*[1 1],[ ], g, LineWidth,2) plot(2*pi*0.25*[1 1],[ ], g, LineWidth,2) In the lower part of Figure 1, we see that the windowing helps and both frequencies can now be seen in the DTFT. Here, we have also used A = 0.05 but it is possible to detect the small peak in the DTFT even if A = Therefore, we conclude that windowing helps and improves the frequency separation somewhat. 3. We start with the residual perspective and compute the loss function of the AR model at different model orders p on validation data. That is, we fit the model using the first 2/3 of the data and compute the prediction error on the last 1/3 of the data. This can be done in Matlab by: %% Generate the data N = 1000; w = randn(1,n); e = randn(1,n); y = filter([1 0.5],[1,-1,0.09],w) + filter([1 0.9],[ ],e); figure(1) subplot(3,2,[1 2]) plot(y); ylabel( signal ); xlabel( time ); %% Model order selection by residual analysis % Divide the data into estimation and validation sets yest = y(1:667); yval = y(668:1000); % Fit the AR(p) for different values of p. Make predictions of the % validation data and compute the prediction error variance resvar = zeros(20,1); for pp = 1:20 m = ar(yest,pp); yhat = predict(m,yval,1); resvar(pp) = mean( (yhat-yval ).^2 ); end subplot(3,2,[3 4]) 3 Ver: 2014/05/06
4 Figur 1: The DTFT of the signal with two sinusoids (upper), the two components (middle) and the windowed signal (lower). The green lines indicate the frequencies of the sinusoids. 4 Ver: 2014/05/06
5 % Plot the prediction error variance. Model order 5 seems nice! plot(1:20,resvar); axis([ ]); ylabel( Prediction error variance ); xlabel( model order (p) ); % Check whiteness of residuals m = ar(yest,5); yhat = predict(m,yval,1); acfresid = xcorr(yhat-yval,25); acfresid = acfresid(26:end); acfresid = acfresid./acfresid(1); subplot(3,2,5); plot(acfresid); % seems uncorrelated hold on; plot([0 30],2/sqrt(N)*[1 1], k: ); plot([0 30],-2/sqrt(N)*[1 1], k: ); xlabel( lag ); ylabel( autocorrelation of residuals ); subplot(3,2,6); qqplot(yhat); % look likes Gaussian (no deviation in the tail behaviour) In Figure 2, we present the data and the corresponding loss function for different values of p. We see that it has a knee point at around p = 5. We check the prediction residuals at this point and conclude that they are uncorrelated and Gaussian. Note that, in a QQ-plot the data (+) should not deviate to much from the theoretical quantiles (red dashed line) especially in the tails. Consequently, we have validated the model assumptions. The model order selection by frequency analysis can be done by the following Matlab code: figure(2) gam = 50; G = etfe(y,gam); [Per,wp]=spectrum(G); Per = squeeze(per); loglog(wp,per, k, LineWidth,2); [mag4,phase,w4] = bode(ar(y,4)); [mag6,phase,w6] = bode(ar(y,6)); [mag8,phase,w8] = bode(ar(y,8)); [mag10,phase,w10] = bode(ar(y,10)); hold on; plot(w4,squeeze(mag4(:)), r, LineWidth,2); plot(w6,squeeze(mag6(:)), b, LineWidth,2); plot(w8,squeeze(mag8(:)), c, LineWidth,2); plot(w10,squeeze(mag10(:)), g, LineWidth,2); legend({ B.-T., AR(4), AR(6), AR(8), AR(10) }, fontsize,16, location, best ) where we try some different model orders. The resulting output is presented in Figure 3, where it is quite difficult to see if any of the models are a better match than the other. As there are a number of peaks in the spectrum of the data, we require at least an AR(8) to cover the first four. Therefore, we conclude that an AR(8) or AR(10) is probably a good choice for this data. 4. a) We select a constant velocity model for this application (see p. 203 in the course book). In two dimensions the state vector is given by x = (p x, p y, v x, v y ) where p and v denote the position and the velocity in each dimension, respectively. The corresponding state space form is given by 1 0 T 0 T 2 / T 0 T 2 /2 x t+1 = x t ( ) y t = x t + e t, T 0 0 T v t, where T denotes the sampling time. Here, we assume that the noise sources are Gaussian, i.e. v t N (0, σ 2 v) and e t N (0, σ 2 e). For this application, we have T = 1 according to the data. b) We make use of a Kalman filter to estimate the position and velocity of the rescue leader given the constant velocity model. This can be done using the following Matlab code: 5 Ver: 2014/05/06
6 Figur 2: The generated data (upper), the loss function at different model orders (middle), the ACF (lower left) and the normal QQ-plot (lower right) of the prediction errors. 6 Ver: 2014/05/06
7 %% Exercise 4 load Apr14; % Noise variances and uncertainty in initial position Q0=0.1; R0=1; P0=[1 1 1e6 1e6]; T=1; % Setup the motion model and load the data N=length(Apr14Meas(:,2)); A=[1 0 T 0; T; ; ]; B=[T^2/2 0; 0 T^2/2; T 0; 0 T]; C=[ ; ]; Q=B*eye(2)*B *Q0; R=eye(2)*R0; P=diag(P0); % Set the initial conditions (zero velocity) xhat=zeros(n,4); yhat=zeros(n,2); xhat(1,:)=[apr14meas(1,1) Apr14Meas(1,2) 0 0]; % Kalman filter for k=2:n % Update Kalman Gain K=P*C /(C*P*C +R); % Eq % Estimate filtered state xhat(k,:)=(a-k*c*a)*xhat(k-1,:) +K*[Apr14Meas(k,1) Apr14Meas(k,2)] ; % Eq yhat(k,:)=c*xhat(k,:) ; % Update covariance estimate P=P-P*C /(C*P*C +R)*C*P; % Eq. 8.17e P=A*P*A +Q; % Eq. 8.17b end % Plot the results figure(4) subplot(2,2,[1 2]); plot(apr14meas(:,1),apr14meas(:,2), k ); hold on; plot(yhat(2:end,1),yhat(2:end,2), k-*, LineWidth,2); title( measured and filtered rescue leader position ); xlabel( x ); ylabel( y ); legend( measured, filtered ); subplot(2,2,3); plot(xhat(1:end,3), k, LineWidth,2); xlabel( time ); ylabel( est. velocity in x-direction ); subplot(2,2,4); plot(xhat(1:end,3), k, LineWidth,2); xlabel( time ); ylabel( est. velocity in y-direction ); The results obtained by running this code is presented in Figure 4. The Kalman filter manages to track the rescue leader s position quite well and gives reasonable estimates of the velocities in each direction. 7 Ver: 2014/05/06
8 Figur 3: The estimated spectrum from the data (black line) and for some different AR models (coloured lines). 8 Ver: 2014/05/06
9 Figur 4: The estimated and filtered position of the rescue leader (upper) and the estimated velocity in each direction (lower). 9 Ver: 2014/05/06
Solutions for examination in TSRT78 Digital Signal Processing,
Solutions for examination in TSRT78 Digital Signal Processing 212-12-2 1. (a) The forgetting factor balances the accuracy and the adaptation speed of the RLS algorithm. This as the forgetting factor controls
More informationω (rad/s)
1. (a) From the figure we see that the signal has energy content in frequencies up to about 15rad/s. According to the sampling theorem, we must therefore sample with at least twice that frequency: 3rad/s
More informationFast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding
Fast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding Fourier Series & Transform Summary x[n] = X[k] = 1 N k= n= X[k]e jkω
More informationFast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding
Fast Fourier Transform Discrete-time windowing Discrete Fourier Transform Relationship to DTFT Relationship to DTFS Zero padding J. McNames Portland State University ECE 223 FFT Ver. 1.03 1 Fourier Series
More informationEEO 401 Digital Signal Processing Prof. Mark Fowler
EEO 4 Digital Signal Processing Pro. Mark Fowler ote Set # Using the DFT or Spectral Analysis o Signals Reading Assignment: Sect. 7.4 o Proakis & Manolakis Ch. 6 o Porat s Book /9 Goal o Practical Spectral
More informationDigital Signal Processing, Lecture 2 Frequency description continued, DFT
Outline cture 2 2 Digital Signal Processing, cture 2 Frequency description continued, DFT Thomas Schön Division of Automatic Control Department of Electrical Engineering Linköping i University it E-mail:
More informationIntroduction to Signal Analysis Parts I and II
41614 Dynamics of Machinery 23/03/2005 IFS Introduction to Signal Analysis Parts I and II Contents 1 Topics of the Lecture 11/03/2005 (Part I) 2 2 Fourier Analysis Fourier Series, Integral and Complex
More informationStatistical methods. Mean value and standard deviations Standard statistical distributions Linear systems Matrix algebra
Statistical methods Mean value and standard deviations Standard statistical distributions Linear systems Matrix algebra Statistical methods Generating random numbers MATLAB has many built-in functions
More informationDiscrete-Time Gaussian Fourier Transform Pair, and Generating a Random Process with Gaussian PDF and Power Spectrum Mark A.
Discrete-Time Gaussian Fourier Transform Pair, and Generating a Random Process with Gaussian PDF and Power Spectrum Mark A. Richards October 3, 6 Updated April 5, Gaussian Transform Pair in Continuous
More informationLecture 17: variance in a band = log(s xx (f)) df (2) If we want to plot something that is more directly representative of variance, we can try this:
UCSD SIOC 221A: (Gille) 1 Lecture 17: Recap We ve now spent some time looking closely at coherence and how to assign uncertainties to coherence. Can we think about coherence in a different way? There are
More informationto have roots with negative real parts, the necessary and sufficient conditions are that:
THE UNIVERSITY OF TEXAS AT SAN ANTONIO EE 543 LINEAR SYSTEMS AND CONTROL H O M E W O R K # 7 Sebastian A. Nugroho November 6, 7 Due date of the homework is: Sunday, November 6th @ :59pm.. The following
More informationSolutions. Chapter 5. Problem 5.1. Solution. Consider the driven, two-well Duffing s oscillator. which can be written in state variable form as
Chapter 5 Solutions Problem 5.1 Consider the driven, two-well Duffing s oscillator which can be written in state variable form as ẍ + ɛγẋ x + x 3 = ɛf cos(ωt) ẋ = v v = x x 3 + ɛ( γv + F cos(ωt)). In the
More informationTemperature measurement
Luleå University of Technology Johan Carlson Last revision: July 22, 2009 Measurement Technology and Uncertainty Analysis - E7021E Lab 3 Temperature measurement Introduction In this lab you are given a
More informationLAB 2: DTFT, DFT, and DFT Spectral Analysis Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 2: DTFT, DFT, and DFT Spectral
More informationEach problem is worth 25 points, and you may solve the problems in any order.
EE 120: Signals & Systems Department of Electrical Engineering and Computer Sciences University of California, Berkeley Midterm Exam #2 April 11, 2016, 2:10-4:00pm Instructions: There are four questions
More informationSystem Identification & Parameter Estimation
System Identification & Parameter Estimation Wb3: SIPE lecture Correlation functions in time & frequency domain Alfred C. Schouten, Dept. of Biomechanical Engineering (BMechE), Fac. 3mE // Delft University
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for solving higher order ODEs and
More informationOverview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros)
Overview of Bode Plots Transfer function review Piece-wise linear approximations First-order terms Second-order terms (complex poles & zeros) J. McNames Portland State University ECE 222 Bode Plots Ver.
More informationBayesian Analysis - A First Example
Bayesian Analysis - A First Example This script works through the example in Hoff (29), section 1.2.1 We are interested in a single parameter: θ, the fraction of individuals in a city population with with
More informationCentre for Mathematical Sciences HT 2017 Mathematical Statistics
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2017 Mathematical Statistics Computer exercise 3 in Stationary stochastic processes, HT 17. The purpose of this exercise
More informationWindowing Overview Relevance Main lobe/side lobe tradeoff Popular windows Examples
ing Overview Relevance Main lobe/side lobe tradeoff Popular windows Examples ing In practice, we cannot observe a signal x(n) for n = to n = Thus we must truncate, or window, the signal: x(n) w(n) ing
More informationThe Kalman Filter. Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience. Sarah Dance
The Kalman Filter Data Assimilation & Inverse Problems from Weather Forecasting to Neuroscience Sarah Dance School of Mathematical and Physical Sciences, University of Reading s.l.dance@reading.ac.uk July
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MAT 75 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP. Use MATLAB solvers for solving higher order ODEs and systems
More information12. Prediction Error Methods (PEM)
12. Prediction Error Methods (PEM) EE531 (Semester II, 2010) description optimal prediction Kalman filter statistical results computational aspects 12-1 Description idea: determine the model parameter
More informationSolutions - Homework # 3
ECE-34: Signals and Systems Summer 23 PROBLEM One period of the DTFS coefficients is given by: X[] = (/3) 2, 8. Solutions - Homewor # 3 a) What is the fundamental period 'N' of the time-domain signal x[n]?
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:305:45 CBC C222 Lecture 8 Frequency Analysis 14/02/18 http://www.ee.unlv.edu/~b1morris/ee482/
More informationSolutions. Chapter 1. Problem 1.1. Solution. Simulate free response of damped harmonic oscillator
Chapter Solutions Problem. Simulate free response of damped harmonic oscillator ẍ + ζẋ + x = for different values of damping ration ζ and initial conditions. Plot the response x(t) and ẋ(t) versus t and
More informationMAE 143B - Homework 7
MAE 143B - Homework 7 6.7 Multiplying the first ODE by m u and subtracting the product of the second ODE with m s, we get m s m u (ẍ s ẍ i ) + m u b s (ẋ s ẋ u ) + m u k s (x s x u ) + m s b s (ẋ s ẋ u
More informationMAT 275 Laboratory 4 MATLAB solvers for First-Order IVP
MATLAB sessions: Laboratory 4 MAT 275 Laboratory 4 MATLAB solvers for First-Order IVP In this laboratory session we will learn how to. Use MATLAB solvers for solving scalar IVP 2. Use MATLAB solvers for
More informationAssignment 4 Solutions Continuous-Time Fourier Transform
Assignment 4 Solutions Continuous-Time Fourier Transform ECE 3 Signals and Systems II Version 1.01 Spring 006 1. Properties of complex numbers. Let c 1 α 1 + jβ 1 and c α + jβ be two complex numbers. a.
More informationProblem Set 1 Solution Sketches Time Series Analysis Spring 2010
Problem Set 1 Solution Sketches Time Series Analysis Spring 2010 1. Construct a martingale difference process that is not weakly stationary. Simplest e.g.: Let Y t be a sequence of independent, non-identically
More informationH 2 Optimal State Feedback Control Synthesis. Raktim Bhattacharya Aerospace Engineering, Texas A&M University
H 2 Optimal State Feedback Control Synthesis Raktim Bhattacharya Aerospace Engineering, Texas A&M University Motivation Motivation w(t) u(t) G K y(t) z(t) w(t) are exogenous signals reference, process
More information2. An Introduction to Moving Average Models and ARMA Models
. An Introduction to Moving Average Models and ARMA Models.1 White Noise. The MA(1) model.3 The MA(q) model..4 Estimation and forecasting of MA models..5 ARMA(p,q) models. The Moving Average (MA) models
More informationFrom Fourier Series to Analysis of Non-stationary Signals - X
From Fourier Series to Analysis of Non-stationary Signals - X prof. Miroslav Vlcek December 14, 21 Contents Stationary and non-stationary 1 Stationary and non-stationary 2 3 Contents Stationary and non-stationary
More informationLaboratory Project 2: Spectral Analysis and Optimal Filtering
Laboratory Project 2: Spectral Analysis and Optimal Filtering Random signals analysis (MVE136) Mats Viberg and Lennart Svensson Department of Signals and Systems Chalmers University of Technology 412 96
More informationDigital Signal Processing: Signal Transforms
Digital Signal Processing: Signal Transforms Aishy Amer, Mohammed Ghazal January 19, 1 Instructions: 1. This tutorial introduces frequency analysis in Matlab using the Fourier and z transforms.. More Matlab
More informationThe Z-Transform. For a phasor: X(k) = e jωk. We have previously derived: Y = H(z)X
The Z-Transform For a phasor: X(k) = e jωk We have previously derived: Y = H(z)X That is, the output of the filter (Y(k)) is derived by multiplying the input signal (X(k)) by the transfer function (H(z)).
More informationOverview of Discrete-Time Fourier Transform Topics Handy Equations Handy Limits Orthogonality Defined orthogonal
Overview of Discrete-Time Fourier Transform Topics Handy equations and its Definition Low- and high- discrete-time frequencies Convergence issues DTFT of complex and real sinusoids Relationship to LTI
More informationProblem Set 9 Solutions
Problem Set 9 Solutions EE23: Digital Signal Processing. From Figure below, we see that the DTFT of the windowed sequence approaches the actual DTFT as the window size increases. Gibb s phenomenon is absent
More informationEvaluation of Uncertainty in AC Power Calculation with Asynchronously Sampled Data
Journal of Physics: Conference Series OPEN ACCESS Evaluation of Uncertainty in AC Power Calculation with Asynchronously Sampled Data To cite this article: D Lindenthaler and H Zangl 01 J. Phys.: Conf.
More informationEE363 homework 5 solutions
EE363 Prof. S. Boyd EE363 homework 5 solutions 1. One-step ahead prediction of an autoregressive time series. We consider the following autoregressive (AR) system p t+1 = αp t + βp t 1 + γp t 2 + w t,
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 informationTTT4120 Digital Signal Processing Suggested Solutions for Problem Set 2
Norwegian University of Science and Technology Department of Electronics and Telecommunications TTT42 Digital Signal Processing Suggested Solutions for Problem Set 2 Problem (a) The spectrum X(ω) can be
More informationEE 565: Position, Navigation, and Timing
EE 565: Position, Navigation, and Timing Kalman Filtering Example Aly El-Osery Kevin Wedeward Electrical Engineering Department, New Mexico Tech Socorro, New Mexico, USA In Collaboration with Stephen Bruder
More information2.161 Signal Processing: Continuous and Discrete
MIT OpenCourseWare http://ocw.mit.edu 2.6 Signal Processing: Continuous and Discrete Fall 2008 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. MASSACHUSETTS
More informationSystem Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang
System Modeling and Identification CHBE 702 Korea University Prof. Dae Ryook Yang 1-1 Course Description Emphases Delivering concepts and Practice Programming Identification Methods using Matlab Class
More informationLABORATORY 1 DISCRETE-TIME SIGNALS
LABORATORY DISCRETE-TIME SIGNALS.. Introduction A discrete-time signal is represented as a sequence of numbers, called samples. A sample value of a typical discrete-time signal or sequence is denoted as:
More information3 rd -order Intercept Point
1 3 rd -order Intercept Point Anuranjan Jha, Columbia Integrated Systems Lab, Department of Electrical Engineering, Columbia University, New York, NY Last Revised: September 12, 2006 These derivations
More informationFourier Analysis of Signals Using the DFT
Fourier Analysis of Signals Using the DFT ECE 535 Lecture April 29, 23 Overview: Motivation Many applications require analyzing the frequency content of signals Speech processing study resonances of vocal
More informationOverview of Sampling Topics
Overview of Sampling Topics (Shannon) sampling theorem Impulse-train sampling Interpolation (continuous-time signal reconstruction) Aliasing Relationship of CTFT to DTFT DT processing of CT signals DT
More informationSummary of lecture 1. E x = E x =T. X T (e i!t ) which motivates us to define the energy spectrum Φ xx (!) = jx (i!)j 2 Z 1 Z =T. 2 d!
Summary of lecture I Continuous time: FS X FS [n] for periodic signals, FT X (i!) for non-periodic signals. II Discrete time: DTFT X T (e i!t ) III Poisson s summation formula: describes the relationship
More informationVelocity estimates from sampled position data
Velocity estimates from sampled position data It is often the case that position information has been sampled from continuous process and it is desired to estimate the velocity and acceleration. Example
More informationCorrelation, discrete Fourier transforms and the power spectral density
Correlation, discrete Fourier transforms and the power spectral density visuals to accompany lectures, notes and m-files by Tak Igusa tigusa@jhu.edu Department of Civil Engineering Johns Hopkins University
More informationDigital Signal Processing Module 1 Analysis of Discrete time Linear Time - Invariant Systems
Digital Signal Processing Module 1 Analysis of Discrete time Linear Time - Invariant Systems Objective: 1. To understand the representation of Discrete time signals 2. To analyze the causality and stability
More informationAspects of Continuous- and Discrete-Time Signals and Systems
Aspects of Continuous- and Discrete-Time Signals and Systems C.S. Ramalingam Department of Electrical Engineering IIT Madras C.S. Ramalingam (EE Dept., IIT Madras) Networks and Systems 1 / 45 Scaling the
More informationLecture 3: Statistical sampling uncertainty
Lecture 3: Statistical sampling uncertainty c Christopher S. Bretherton Winter 2015 3.1 Central limit theorem (CLT) Let X 1,..., X N be a sequence of N independent identically-distributed (IID) random
More informationSpectrograms Overview Introduction and motivation Windowed estimates as filter banks Uncertainty principle Examples Other estimates
Spectrograms Overview Introduction and motivation Windowed estimates as filter banks Uncertainty principle Examples Other estimates J. McNames Portland State University ECE 3 Spectrograms Ver. 1.1 1 Introduction
More informationModel-building and parameter estimation
Luleå University of Technology Johan Carlson Last revision: July 27, 2009 Measurement Technology and Uncertainty Analysis - E7021E MATLAB homework assignment Model-building and parameter estimation Introduction
More informationLab 4: Quantization, Oversampling, and Noise Shaping
Lab 4: Quantization, Oversampling, and Noise Shaping Due Friday 04/21/17 Overview: This assignment should be completed with your assigned lab partner(s). Each group must turn in a report composed using
More informationA time series is called strictly stationary if the joint distribution of every collection (Y t
5 Time series A time series is a set of observations recorded over time. You can think for example at the GDP of a country over the years (or quarters) or the hourly measurements of temperature over a
More informationINFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters 4. THE BUTTERWORTH ANALOG FILTER
INFINITE-IMPULSE RESPONSE DIGITAL FILTERS Classical analog filters and their conversion to digital filters. INTRODUCTION 2. IIR FILTER DESIGN 3. ANALOG FILTERS 4. THE BUTTERWORTH ANALOG FILTER 5. THE CHEBYSHEV-I
More informationCorrelator I. Basics. Chapter Introduction. 8.2 Digitization Sampling. D. Anish Roshi
Chapter 8 Correlator I. Basics D. Anish Roshi 8.1 Introduction A radio interferometer measures the mutual coherence function of the electric field due to a given source brightness distribution in the sky.
More informationLecture 16: State Space Model and Kalman Filter Bus 41910, Time Series Analysis, Mr. R. Tsay
Lecture 6: State Space Model and Kalman Filter Bus 490, Time Series Analysis, Mr R Tsay A state space model consists of two equations: S t+ F S t + Ge t+, () Z t HS t + ɛ t (2) where S t is a state vector
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 information3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE
3. ESTIMATION OF SIGNALS USING A LEAST SQUARES TECHNIQUE 3.0 INTRODUCTION The purpose of this chapter is to introduce estimators shortly. More elaborated courses on System Identification, which are given
More informationComplement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises
Complement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises Random Processes With Applications (MVE 135) Mats Viberg Department of Signals and Systems Chalmers University of Technology
More informationESSE 4020 / ESS 5020 Time Series and Spectral Analysis Notes 6 26 November 2018 State Space Models
ESSE 4020 / ESS 5020 Time Series and Spectral Analysis Notes 6 26 November 2018 State Space Models SS2011/2017 Chapter 6 State Space Models - Kalman Filtering or 2016 at https://www.stat.pitt.edu/stoffer/tsa4/chapter6.pdf
More informationSolutions - Homework #2
45 Scatterplot of Abundance vs. Relative Density Parasite Abundance 4 35 3 5 5 5 5 5 Relative Host Population Density Figure : 3 Scatterplot of Log Abundance vs. Log RD Log Parasite Abundance 3.5.5.5.5.5
More informationsinc function T=1 sec T=2 sec angle(f(w)) angle(f(w))
T=1 sec sinc function 3 angle(f(w)) T=2 sec angle(f(w)) 1 A quick script to plot mag & phase in MATLAB w=0:0.2:50; Real exponential func b=5; Fourier transform (filter) F=1.0./(b+j*w); subplot(211), plot(w,
More informationMinitab Project Report - Assignment 6
.. Sunspot data Minitab Project Report - Assignment Time Series Plot of y Time Series Plot of X y X 7 9 7 9 The data have a wavy pattern. However, they do not show any seasonality. There seem to be an
More informationExtended Kalman Filter Derivation Example application: frequency tracking
Extended Kalman Filter Derivation Example application: frequency tracking Nonlinear State Space Model Consider the nonlinear state-space model of a random process x n+ = f n (x n )+g n (x n )u n y n =
More informationCHAPTER 7. The Discrete Fourier Transform 436
CHAPTER 7. The Discrete Fourier Transform 36 hfa figconfg( P7a, long ); subplot() stem(k,abs(ck), filled );hold on stem(k,abs(ck_approx), filled, color, red ); xlabel( k, fontsize,lfs) title( Magnitude
More informationLearning Deep Broadband Hongjoo LEE
Learning Deep Broadband Network@HOME Hongjoo LEE Who am I? Machine Learning Engineer Software Engineer Fraud Detection System Software Defect Prediction Email Services (40+ mil. users) High traffic server
More informationElec4621 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis
Elec461 Advanced Digital Signal Processing Chapter 11: Time-Frequency Analysis Dr. D. S. Taubman May 3, 011 In this last chapter of your notes, we are interested in the problem of nding the instantaneous
More informationComplement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises
Complement on Digital Spectral Analysis and Optimal Filtering: Theory and Exercises Random Signals Analysis (MVE136) Mats Viberg Department of Signals and Systems Chalmers University of Technology 412
More informationA=randn(500,100); mu=mean(a); sigma_a=std(a); std_a=sigma_a/sqrt(500); [std(mu) mean(std_a)] % compare standard deviation of means % vs standard error
UCSD SIOC 221A: (Gille) 1 Reading: Bendat and Piersol, Ch. 5.2.1 Lecture 10: Recap Last time we looked at the sinc function, windowing, and detrending with an eye to reducing edge effects in our spectra.
More informationThe goal of the Wiener filter is to filter out noise that has corrupted a signal. It is based on a statistical approach.
Wiener filter From Wikipedia, the free encyclopedia In signal processing, the Wiener filter is a filter proposed by Norbert Wiener during the 1940s and published in 1949. [1] Its purpose is to reduce the
More informationCentre for Mathematical Sciences HT 2017 Mathematical Statistics. Study chapters 6.1, 6.2 and in the course book.
Lund University Stationary stochastic processes Centre for Mathematical Sciences HT 2017 Mathematical Statistics Computer exercise 2 in Stationary stochastic processes, HT 17. The purpose with this computer
More informationDSP Applications for Wireless Communications: Linear Equalisation of MIMO Channels
DSP Applications for Wireless Communications: Dr. Waleed Al-Hanafy waleed alhanafy@yahoo.com Faculty of Electronic Engineering, Menoufia Univ., Egypt Digital Signal Processing (ECE407) Lecture no. 5 August
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 1 Time-Dependent FT Announcements! Midterm: 2/22/216 Open everything... but cheat sheet recommended instead 1am-12pm How s the lab going? Frequency Analysis with
More informationComputer Exercise 1 Estimation and Model Validation
Lund University Time Series Analysis Mathematical Statistics Fall 2018 Centre for Mathematical Sciences Computer Exercise 1 Estimation and Model Validation This computer exercise treats identification,
More informationat least 50 and preferably 100 observations should be available to build a proper model
III Box-Jenkins Methods 1. Pros and Cons of ARIMA Forecasting a) need for data at least 50 and preferably 100 observations should be available to build a proper model used most frequently for hourly or
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 informationClassical Decomposition Model Revisited: I
Classical Decomposition Model Revisited: I recall classical decomposition model for time series Y t, namely, Y t = m t + s t + W t, where m t is trend; s t is periodic with known period s (i.e., s t s
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 informationCore Concepts Review. Orthogonality of Complex Sinusoids Consider two (possibly non-harmonic) complex sinusoids
Overview of Continuous-Time Fourier Transform Topics Definition Compare & contrast with Laplace transform Conditions for existence Relationship to LTI systems Examples Ideal lowpass filters Relationship
More informationEE 224 Lab: Graph Signal Processing
EE Lab: Graph Signal Processing In our previous labs, we dealt with discrete-time signals and images. We learned convolution, translation, sampling, and frequency domain representation of discrete-time
More informationNotes perso. - cb1 1. Version 3 - November u(t)e i2πft dt u(t) =
Notes perso. - cb1 1 Tutorial : fft, psd & coherence with Matlab Version 3 - November 01 1 Fourier transform by Fast Fourier Transform (FFT) Definition of the Fourier transform : û(f) = F[u(t)] = As an
More informationStochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents
Navtech Part #s Volume 1 #1277 Volume 2 #1278 Volume 3 #1279 3 Volume Set #1280 Stochastic Models, Estimation and Control Peter S. Maybeck Volumes 1, 2 & 3 Tables of Contents Volume 1 Preface Contents
More informationJ. McNames Portland State University ECE 223 Sampling Ver
Overview of Sampling Topics (Shannon) sampling theorem Impulse-train sampling Interpolation (continuous-time signal reconstruction) Aliasing Relationship of CTFT to DTFT DT processing of CT signals DT
More informationProf. Dr. Roland Füss Lecture Series in Applied Econometrics Summer Term Introduction to Time Series Analysis
Introduction to Time Series Analysis 1 Contents: I. Basics of Time Series Analysis... 4 I.1 Stationarity... 5 I.2 Autocorrelation Function... 9 I.3 Partial Autocorrelation Function (PACF)... 14 I.4 Transformation
More informationLesson 1. Optimal signalbehandling LTH. September Statistical Digital Signal Processing and Modeling, Hayes, M:
Lesson 1 Optimal Signal Processing Optimal signalbehandling LTH September 2013 Statistical Digital Signal Processing and Modeling, Hayes, M: John Wiley & Sons, 1996. ISBN 0471594318 Nedelko Grbic Mtrl
More informationLecture 4 - Spectral Estimation
Lecture 4 - Spectral Estimation The Discrete Fourier Transform The Discrete Fourier Transform (DFT) is the equivalent of the continuous Fourier Transform for signals known only at N instants separated
More informationReliability and Risk Analysis. Time Series, Types of Trend Functions and Estimates of Trends
Reliability and Risk Analysis Stochastic process The sequence of random variables {Y t, t = 0, ±1, ±2 } is called the stochastic process The mean function of a stochastic process {Y t} is the function
More informationEcon 424 Time Series Concepts
Econ 424 Time Series Concepts Eric Zivot January 20 2015 Time Series Processes Stochastic (Random) Process { 1 2 +1 } = { } = sequence of random variables indexed by time Observed time series of length
More informationUniversiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 3: DISCRETE TIME SYSTEM IN TIME DOMAIN
Universiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 3: DISCRETE TIME SYSTEM IN TIME DOMAIN Pusat Pengajian Kejuruteraan Komputer Dan Perhubungan Universiti Malaysia Perlis Discrete-Time
More informationLevinson Durbin Recursions: I
Levinson Durbin Recursions: I note: B&D and S&S say Durbin Levinson but Levinson Durbin is more commonly used (Levinson, 1947, and Durbin, 1960, are source articles sometimes just Levinson is used) recursions
More informationHW13 Solutions. Pr (a) The DTFT of c[n] is. C(e jω ) = 0.5 m e jωm e jω + 1
HW3 Solutions Pr..8 (a) The DTFT of c[n] is C(e jω ) = = =.5 n e jωn + n=.5e jω + n= m=.5 n e jωn.5 m e jωm.5e jω +.5e jω.75 =.5 cos(ω) }{{}.5e jω /(.5e jω ) C(e jω ) is the power spectral density. (b)
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More informationCCNY. BME I5100: Biomedical Signal Processing. Stochastic Processes. Lucas C. Parra Biomedical Engineering Department City College of New York
BME I5100: Biomedical Signal Processing Stochastic Processes Lucas C. Parra Biomedical Engineering Department CCNY 1 Schedule Week 1: Introduction Linear, stationary, normal - the stuff biology is not
More information