Signal Processing Summary
|
|
- Alyson Parks
- 5 years ago
- Views:
Transcription
1 Signal Processing Summary Jan Černocký, Valentina Hubeika DCGM FIT BUT Brno, FIT BUT Brno Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 1/41
2 Agenda Introduction Analog signals and sampling Frequency analysis of sampled signals. Random signals. Filters. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 2/41
3 Why digital signal processing? reproducibility (we don t have to think about component tolerance ). no changes due to material aging or temperature. easy setup (may be difficult tuning for analog signals). possible adaptive processing ( changing device functionality depending on the type of the signal ). simulation = application. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 3/41
4 Usual Approach x(t) A/D x (n) digital processing y (n) D/A y(t) (PC, DSP). storing transmision interpretation other processing... The input signal is continuous: function of a real-valued variable (in our case time) t, t is defined over the interval (, ), t. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 4/41
5 s(t) t To represent the signal in frequency domain (spectrum), Fourier transform is used: X(f) = x(t)e j2πft dt, (1) where the function X(f) is called the spectral function, or in short spectrum. The function is defined f from to and is complex. Its magnitude is X(f) and its phase is X(f). Here, we talk about magnitude and phase spectra. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 5/41
6 In case of real signals, we will need only the right side of the spectral function (f > ). The left side is complex adjoint to the right side: or X(f) = X( f) a arg X(f) = arg X( f). X(f) = X ( f), (2) S(f) -f max f f max Good signals are frequency constrained (energy within the band of (, f max ) ). Spectral function cant be calculated (infinite, integral,...). Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 6/41
7 Analog to Digital (AD) Converter x(t) antialias. filter sampling sampled signal xs(n) quantization quantized signal xq(n) Sampled signal is calculated by multiplying the original signal with a periodic function over time. s(t) D s(t) D x (t)=x(t)s(t) s 1/D T t t Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 7/41
8 Theoretically, sampling is multiplying the original analog signal with a series od Dirac deltas (unit impulse function) (infinite height, zero width, unit area). The result of multiplication is again a series of Dirac impulses, with the powers corresponding to the values of the original signal at times nt: s(t) s(t) 1 D 1 x (t)=x(t)s(t) s T t t T is the sampling period F s = 1 T is the sampling frequency Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 8/41
9 The spectrum of the sampled signal is periodized!!! X s (f) = 1 T + X ( f n T ) = 1 T + X (f nf s ) (3) n= n= From the point of view of the sampling theorem: 1) F s > 2f max : The copies of the spectrum of the original signal do not overlap so the original signal can be perfectly reconstructed from its sampled version by filtering it with low-pass filter with the frequency F s /2. 2) F s 2f max : the copies of the spectrum of the original signal do overlap so the resulting spectrum has a different shape than the original spectrum. The original signal cannot be reconstructed from its sampled version. We observe so called aliasing. Shanon Kotelnikov Nyquist sampling theorem F s > 2f max Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 9/41
10 1. Example: sampling and reconstruction OK F s = 8 Hz, f max = 3 Hz, so that Ω s = 16π rad/s, ω max = 6π rad/s. T = 1 8 s x x x x x 1 5 Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 1/41
11 2. Example: sampling and reconstruction BAD F s = 8 Hz, f max = 7 Hz, so that Ω s = 16π rad/s, ω max = 14π rad/s. T = 1 8 s x x x x x 1 5 Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 11/41
12 Antialiasing filter constrain to [ F s /2, F s /2] x x x x x x 1 5 Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 12/41
13 Denotation of the Sampled Signal x s (nt) or just x s [n] the sampled signal is a series of values. 1. Along with the sampled signal, we always need information on the sampling frequency (implicit or explicit) 2. When working with a sampled signal, we do not want to consider the sampling frequency in the calculations, therefore assume that the period is T = 1, thus F s = 1. The normalized time is then given as: The normalized frequency is: t = t T, thus n = nt T (4) f = f F s (5) Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 13/41
14 Example Define a function generating cosine with the frequency of 2 Hz with the sampling frequency F s = 8 Hz. In the continuous time: s(t) = cos(2πf t) = cos(2π2t). While sampling we will substitute the continuous time variable( t with the discrete time nt, where T is the sampling period: x(nt) = cos(2πf nt) = cos 2π f ) n. F s The frequency f ( is normalized. The desired signal is defined as: x(n) = cos 2π 1 ) F s 4 n. Generate 1 second of such a signal in Matlab: n = :7999; x = cos (2 * pi * 1 / 4 * n); wavwrite(x,8,16, sig.wav ); Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 14/41
15 Behaviour of the Sampled Signal in Frequency Domain Spectrum Discrete Fourier transform DTFT definition: X(k) = N 1 n= nk j2π x[n]e N pro k <, N 1 > (6) How to use it for discrete signals?: analyze a window of the length of N samples. result: when multiplying values X(k) with the sampling perioud T, the approximation of the spectral function at the points k f, where f = F s N f = 1 N (normalized frequency (usually not used) (true frequency) or ˆX(k f) = T N 1 n= nk j2π x[n]e N (7) in this definition, x[n] can be substitute for x(nt). Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 15/41
16 Comparing to the spectral function of the corresponding analog signal, the spectrum of the discrete signal is not the same thing! 1. we calculate spectrum of the sampled signal, so the spectrum is periodic with the period N units (corresponds to the sampling frequency F s ). If k (, + ), then ˆX(k f) is repeated after each N values. 2. The signal is windowed. The resulting spectrum carries information on the window: the signal is multiplied with the window function in time, the spectrum of the window is then convoluted with the spectrum of the original signal. This often results in blurring of the theoretically sharp lines of the spectrum (for instance during analysis of the harmonic signal). More in the lecture on preprocessing. 3. if the spectrum is discrete (we dispose of only N values from to F s ), then the signal must be periodic! time sampling periodization frequency periodization discretization Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 16/41
17 How can we use it in practice? We want to make frequency analysis of one voiced frame of speech s = wavread( test.wav ) ; sfr = frame (s,16,8); x = sfr(:,13); plot (x); Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 17/41
18 DTFT only careful with the frequency axes: Fs = 8; f = (:159) / 16 * Fs; X = fft(x); subplot (211); plot(f,abs(x)); subplot (212); plot(f,angle(x)); Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 18/41
19 DTFT only the left half is symmetric to the right half (we consider only on in analysis) Fs = 8; f = (:79) / 16 * Fs; X = fft(x); X = X(1:8); subplot (211); plot(f,abs(x)); subplot (212); plot(f,angle(x)); Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 19/41
20 If we want to smoother the spectrum (add more points shorten intervals between values), having only one frame (without a possibility of extension), we can use zero padding Fs = 8; f = (:511) / 124 * Fs; X = fft([x zeros(1,124-16)]); X = X(1:512); subplot (211); plot(f,abs(x)); subplot (212); plot(f,angle(x)); Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 2/41
21 FREQUENCY ANALYSIS OF RANDOM SIGNALS Speech signal is a random signal. Thus we use power spectral density (PSD) function for its frequency analysis. PSD is real and represent power density of the signal in frequency domain. One of the possible ways of computation uses DTFT: Which is just the power of the magnitude. Ĝ DTFT (k f) = 1 N X[k] 2. careful with X[k] 2. in Matlab power of the complex number is not the same as the power of its magnitude: X = fft(x); Gdft = X.^ 2; wrong! X = fft(x); Gdft = abs(x).^ 2; good! X = fft(x); Gdft = X.* conj(x); good and fast Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 21/41
22 Example: Fs = 8; f = (:511) / 124 * Fs X = fft([x zeros(1,124-16)]); X = X(1:512); Gdft= 1/16 *abs(x).^ 2; plot(f,gdft); Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 22/41
23 Dynamic performance of the power spectral density is higher than of DTFT (second power) which causes poor visualization of weak parts of the image. For this reason, often decibels are used (Matlab: function log1) Ĝ DTFT (k f) = 1 log 1 1 N X[k] Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 23/41
24 LINEAR FILTERS x(n) filter y(n) Usually, filters are linear hold linear combination : if x 1 (n) y 1 (n) a x 2 (n) y 2 (n), then ax 1 (n) + bx 2 (n) ay 1 (n) + by 2 (n), where a, b R. time invariant behaviour does not depend on time: if x(n) y(n), then as well x(n n ) y ( n n ), where n is any arbitrary shift. Sometimes however we desire filters with time variant behaviour adaptive systems, speech frames (change 1 ms). causal the filter does not see into future : y(n) y(m < n) a x(m n). Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 24/41
25 Impulse Response or impulse characteristics is the response of a filter to the Kroneck or unit impulse (do not confuse with the Dirac delta): pro n δ(n) = (8) 1 pro n = δ(n) h(n) filter Knowing the impulse response, we can calculate the response to any input signal. Each sample is perceived as a shifted and scaled unit impulse. We represent the output as the convolution: y(n) = x(n) h(n) = m= x(m)h(n m) = m= h(m)x(n m) (9) Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 25/41
26 Analysis of the impulse characteristics: if h(k) = for k <, then the filter is causal (all samples after n-th sample will be multiplied by zero). impulse response can be finite FIR (finite impulse response) or infinite IIR (infinite impulse response). the Fourier image in frequency then represents the complex characteristics of the filter: h(k) H(f) Convolution in time domain correspond to product in frequency domain, thus the spectrum of the resulting signal is: Y (f) = X(f)H(f) (1) We should keep in mind that we work with the discrete signals (so is the impulse response of a filter) the image in the frequency is therefore periodic with the period F s or 1 for the normalized frequency. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 26/41
27 Representation of a Filter y(n) x(n) z -1 z z -1 b Q Σ... z -1 z -1 z -1 b Q-1 b Q-2 -a -a b 1 -a P-1 b -a P The block z 1 denotes a delay of 1 sample. The behaviour of a filter can be expressed as a difference equation: y(n) = Q b k x(n k) k= P a k y(n k), (11) Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 27/41 k=1
28 where x(n k) are actual and delayed input values and y(n k) are delayed output values. Types of filters: FIR non-recursive: only b...b Q nonzero. Always stable. IIR solely recursive: only b, a 1...a P nonzero. IIR generally recursive: a i and b i nonzero. It is difficult to estimate the behaviour of the filter and analyze it from the point of view of stability from the difference equation. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 28/41
29 Z-Transform X(z) = n= x(n)z n (12) The z-transform is a generalization of the DTFT. The z-transform becomes the DTFT when z = exp jw (when DTFT is evaluated on the unit circle). In order to determine the frequency response of the system the z-transform must be evaluated on the circle (the region of convergence must contain a circle). Otherwise, the DTFT of the system does not exist. Lets define delay as: if x(n) X(z), then for y(n) = x(n n ) will hold: Y (z) = z n X(z) (13) for the delay of one time unit holds: x(n 1) z 1 X(z). Thus the delay of one time unit will be denoted as: z -1 Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 29/41
30 Transfer Function of a Filter We can rewrite the difference equation using z-transform Y (z) = Q b k X(z)z k k= P a k Y (z)z k, (14) k=1 The transfer function can be defined as the fraction: H(z) = Y (z) X(z) = 1 + Q b k z k k= = P a k z k k=1 B(z) A(z), (15) where B(z) and A(z) are two polynomes. The a coefficient of the denominator must be equal to 1. a does not physically present in the filter (denotes that the filter has the output). Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 3/41
31 Frequency Characteristics of a filter from to F s (or from to 1 in using normalized frequency) can be acquired from the transfer function by analyzing the unit circle. We are interested in the complex values of the function H(z): for the normalized frequency. Or: H(f) = H(z) z=e j2πf (16) H(f) = H(z) z=e j2πft (17) for the normal frequency. For each value f we need to calculate the position of the point on the unit circle: z = e j2πf (complex number), then for this number we need to compute the fraction of the polynomes B(z) and A(z) (complex number). In Matlab by using function freqz. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 31/41
32 +j Im "z" exp(j2 πf) 1 1/2~Fs/2 1 1~Fs Re j Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 32/41
33 Zeros and Poles of the Transmission Function... The other way how to represent the transfer function H(z) is using product: H(z) = B(z) A(z) = b + b 1 z b Q z Q 1 + a 1 z a P z P = z Q (b z Q + b 1 z Q b Q ) z P (z P + a 1 z P a P ) = = b z Q z P Q (z n k ) k=1 P (z p k ) k=1 = b z P Q Q (z n k ) k=1 P (z p k ), If a k, b k R, then poles p k and zeros n k can be either real or complex conjugate in pairs. Filter is stable when all poles lie within the unity circle: k=1 p k < 1 Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 33/41
34 Example of a Filter We want to implement a filter that will be simulating a telephone channel to filter a signal of the CD quality. Our filter will be a band-pass filter from 3 to 34 Hz. There are several functions available in Matlab to implement a filter. We will be using so called elliptic filters: Fs = 441; Fs2 = Fs/2; % has to be normalized by half Fs Wp = [3/Fs2 34/Fs2]; % pass-band Ws = [2/Fs2 35/Fs2]; % stop-band - approximately Rp = 3; % ripple in pass-bandu db Rs = 3; % muting stop-bandu db (obe hodnoty od % see the norms. [N, Wn] = ellipord(wp, Ws, Rp, Rs) % computing of the filter order [B,A] = ellip(n,rp,rs,wn) % computing of polynoms B a A...as the result we have 2 polynomes of the 12th order. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 34/41
35 Frequency characteristics: freqz (B,A,512,Fs); 5 Magnitude (db) Frequency (Hz) x 1 4 Phase (degrees) Frequency (Hz) x 1 4 Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 35/41
36 Poles and zeros: zplane (B,A); Imaginary Part Real Part Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 36/41
37 Implementation in C the basic implementation of the direct structure is very simple just rewrite the difference equation into C (see the file filter.c) in practice, usually optimal structures are used, which have only one delay bus and are not prone to the rounding errors. for more information on the filter theory see SXC/ISS lecture discrete systems : Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 37/41
38 Passing a Random Signal through a Filter Let the filter have a complex frequency characteristic H(f). For the input signal with the spectral power density G x (f) the output spectral power density is given as: G y (f) = H(f) 2 G x (f)...the input PSD is multiplied by the power of the magnitude of the complex frequency characteristic. Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 38/41
39 Example: Given a realization of a random signal. The signal is filtered by the filter defined as H(z) = 1.9z 1. The output ignal and its PSD are: Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 39/41
40 The magnitude of the frequency characteristic and its power are: H H Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 4/41
41 The output signal and its PSD: Signal Processing Summary Jan Černocký, Valentina Hubeika, CDGM FIT BUT Brno 41/41
Signal Processing for Speech Applications - Part 2-1. Signal Processing For Speech Applications - Part 2
Signal Processing for Speech Applications - Part 2-1 Signal Processing For Speech Applications - Part 2 May 14, 2013 Signal Processing for Speech Applications - Part 2-2 References Huang et al., Chapter
More informationSMS045 - DSP Systems in Practice. Lab 1 - Filter Design and Evaluation in MATLAB Due date: Thursday Nov 13, 2003
SMS045 - DSP Systems in Practice Lab 1 - Filter Design and Evaluation in MATLAB Due date: Thursday Nov 13, 2003 Lab Purpose This lab will introduce MATLAB as a tool for designing and evaluating digital
More informationECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015
Purdue University: ECE438 - Digital Signal Processing with Applications 1 ECE438 - Laboratory 7a: Digital Filter Design (Week 1) By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015 1 Introduction
More informationELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet
ELEC-C5230 Digitaalisen signaalinkäsittelyn perusteet Lecture 10: Summary Taneli Riihonen 16.05.2016 Lecture 10 in Course Book Sanjit K. Mitra, Digital Signal Processing: A Computer-Based Approach, 4th
More informationTopic 2. Signal Processing Review. (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music)
Topic 2 Signal Processing Review (Some slides are adapted from Bryan Pardo s course slides on Machine Perception of Music) Recording Sound Mechanical Vibration Pressure Waves Motion->Voltage Transducer
More informationy(n)= Aa n u(n)+bu(n) b m sin(2πmt)= b 1 sin(2πt)+b 2 sin(4πt)+b 3 sin(6πt)+ m=1 x(t)= x = 2 ( b b b b
Exam 1 February 3, 006 Each subquestion is worth 10 points. 1. Consider a periodic sawtooth waveform x(t) with period T 0 = 1 sec shown below: (c) x(n)= u(n). In this case, show that the output has the
More informationConcordia University. Discrete-Time Signal Processing. Lab Manual (ELEC442) Dr. Wei-Ping Zhu
Concordia University Discrete-Time Signal Processing Lab Manual (ELEC442) Course Instructor: Dr. Wei-Ping Zhu Fall 2012 Lab 1: Linear Constant Coefficient Difference Equations (LCCDE) Objective In this
More informationLecture 3, Multirate Signal Processing
Lecture 3, Multirate Signal Processing Frequency Response If we have coefficients of an Finite Impulse Response (FIR) filter h, or in general the impulse response, its frequency response becomes (using
More informationTeam proposals are due tomorrow at 6PM Homework 4 is due next thur. Proposal presentations are next mon in 1311EECS.
Lecture 8 Today: Announcements: References: FIR filter design IIR filter design Filter roundoff and overflow sensitivity Team proposals are due tomorrow at 6PM Homework 4 is due next thur. Proposal presentations
More informationEE 470 Signals and Systems
EE 470 Signals and Systems 9. Introduction to the Design of Discrete Filters Prof. Yasser Mostafa Kadah Textbook Luis Chapparo, Signals and Systems Using Matlab, 2 nd ed., Academic Press, 2015. Filters
More informationPROBLEM SET 6. Note: This version is preliminary in that it does not yet have instructions for uploading the MATLAB problems.
PROBLEM SET 6 Issued: 2/32/19 Due: 3/1/19 Reading: During the past week we discussed change of discrete-time sampling rate, introducing the techniques of decimation and interpolation, which is covered
More informationCS3291: Digital Signal Processing
CS39 Exam Jan 005 //08 /BMGC University of Manchester Department of Computer Science First Semester Year 3 Examination Paper CS39: Digital Signal Processing Date of Examination: January 005 Answer THREE
More informationLecture 3 Review of Signals and Systems: Part 2. EE4900/EE6720 Digital Communications
EE4900/EE6720: Digital Communications 1 Lecture 3 Review of Signals and Systems: Part 2 Block Diagrams of Communication System Digital Communication System 2 Informatio n (sound, video, text, data, ) Transducer
More informationOutline. Discrete time signals. Impulse sampling z-transform Frequency response Stability INF4420. Jørgen Andreas Michaelsen Spring / 37 2 / 37
INF4420 Discrete time signals Jørgen Andreas Michaelsen Spring 2013 1 / 37 Outline Impulse sampling z-transform Frequency response Stability Spring 2013 Discrete time signals 2 2 / 37 Introduction More
More informationSignals and Systems program and organization
Signals and Systems program and organization Valentina Hubeika, Jan Černocký DCGM FIT BUT {ihubeika cernocky}@fit.vutbr.cz organization goals motivation examples of signal processing program of the course
More informationLecture Schedule: Week Date Lecture Title
http://elec3004.org Sampling & More 2014 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date Lecture Title 1 2-Mar Introduction 3-Mar
More informationECE 429 / 529 Digital Signal Processing
ECE 429 / 529 Course Policy & Syllabus R. N. Strickland SYLLABUS ECE 429 / 529 Digital Signal Processing SPRING 2009 I. Introduction DSP is concerned with the digital representation of signals and the
More informationB.Tech III Year II Semester (R13) Regular & Supplementary Examinations May/June 2017 DIGITAL SIGNAL PROCESSING (Common to ECE and EIE)
Code: 13A04602 R13 B.Tech III Year II Semester (R13) Regular & Supplementary Examinations May/June 2017 (Common to ECE and EIE) PART A (Compulsory Question) 1 Answer the following: (10 X 02 = 20 Marks)
More information1. In the command window, type "help conv" and press [enter]. Read the information displayed.
ECE 317 Experiment 0 The purpose of this experiment is to understand how to represent signals in MATLAB, perform the convolution of signals, and study some simple LTI systems. Please answer all questions
More informationDigital Processing of Continuous-Time Signals
Chapter 4 Digital Processing of Continuous-Time Signals 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 Digital Processing of Continuous-Time Signals Digital
More informationSignals. Continuous valued or discrete valued Can the signal take any value or only discrete values?
Signals Continuous time or discrete time Is the signal continuous or sampled in time? Continuous valued or discrete valued Can the signal take any value or only discrete values? Deterministic versus random
More informationDigital Processing of
Chapter 4 Digital Processing of Continuous-Time Signals 清大電機系林嘉文 cwlin@ee.nthu.edu.tw 03-5731152 Original PowerPoint slides prepared by S. K. Mitra 4-1-1 Digital Processing of Continuous-Time Signals Digital
More informationIslamic University of Gaza. Faculty of Engineering Electrical Engineering Department Spring-2011
Islamic University of Gaza Faculty of Engineering Electrical Engineering Department Spring-2011 DSP Laboratory (EELE 4110) Lab#4 Sampling and Quantization OBJECTIVES: When you have completed this assignment,
More informationBiomedical Signals. Signals and Images in Medicine Dr Nabeel Anwar
Biomedical Signals Signals and Images in Medicine Dr Nabeel Anwar Noise Removal: Time Domain Techniques 1. Synchronized Averaging (covered in lecture 1) 2. Moving Average Filters (today s topic) 3. Derivative
More informationLecture 2 Review of Signals and Systems: Part 1. EE4900/EE6720 Digital Communications
EE4900/EE6420: Digital Communications 1 Lecture 2 Review of Signals and Systems: Part 1 Block Diagrams of Communication System Digital Communication System 2 Informatio n (sound, video, text, data, ) Transducer
More informationEECS 452 Practice Midterm Exam Solutions Fall 2014
EECS 452 Practice Midterm Exam Solutions Fall 2014 Name: unique name: Sign the honor code: I have neither given nor received aid on this exam nor observed anyone else doing so. Scores: # Points Section
More information1. Find the magnitude and phase response of an FIR filter represented by the difference equation y(n)= 0.5 x(n) x(n-1)
Lecture 5 1.8.1 FIR Filters FIR filters have impulse responses of finite lengths. In FIR filters the present output depends only on the past and present values of the input sequence but not on the previous
More informationEE 422G - Signals and Systems Laboratory
EE 422G - Signals and Systems Laboratory Lab 3 FIR Filters Written by Kevin D. Donohue Department of Electrical and Computer Engineering University of Kentucky Lexington, KY 40506 September 19, 2015 Objectives:
More informationSignals and Systems Lecture 6: Fourier Applications
Signals and Systems Lecture 6: Fourier Applications Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Winter 2012 arzaneh Abdollahi Signal and Systems Lecture 6
More informationZRE 01- Introductory Lab JanČernocký,ValentinaHubeika,FITBUTBrno
ZRE 01- Introductory Lab JanČernocký,ValentinaHubeika,FITBUTBrno Thegoalofthelabistointroducebasicoperationswithasoundsignal,mainlyusingMatlab.Thislab(aswell asallfollowinglabs)willberununderlinux.althoughmatlabisavailableunderwindowsosaswell,someofthe
More informationece 429/529 digital signal processing robin n. strickland ece dept, university of arizona ECE 429/529 RNS
ece 429/529 digital signal processing robin n. strickland ece dept, university of arizona 2007 SPRING 2007 SCHEDULE All dates are tentative. Lesson Day Date Learning outcomes to be Topics Textbook HW/PROJECT
More informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Basic Sampling Rate Alteration Devices Up-sampler - Used to increase the sampling rate by an integer factor Down-sampler - Used to increase the sampling rate by an integer
More informationUNIT-II MYcsvtu Notes agk
UNIT-II agk UNIT II Infinite Impulse Response Filter design (IIR): Analog & Digital Frequency transformation. Designing by impulse invariance & Bilinear method. Butterworth and Chebyshev Design Method.
More informationSignal processing preliminaries
Signal processing preliminaries ISMIR Graduate School, October 4th-9th, 2004 Contents: Digital audio signals Fourier transform Spectrum estimation Filters Signal Proc. 2 1 Digital signals Advantages of
More informationExperiment 8: Sampling
Prepared By: 1 Experiment 8: Sampling Objective The objective of this Lab is to understand concepts and observe the effects of periodically sampling a continuous signal at different sampling rates, changing
More informationWeek 1 Introduction of Digital Signal Processing with the review of SMJE 2053 Circuits & Signals for Filter Design
SMJE3163 DSP2016_Week1-04 Week 1 Introduction of Digital Signal Processing with the review of SMJE 2053 Circuits & Signals for Filter Design 1) Signals, Systems, and DSP 2) DSP system configuration 3)
More informationElectrical & Computer Engineering Technology
Electrical & Computer Engineering Technology EET 419C Digital Signal Processing Laboratory Experiments by Masood Ejaz Experiment # 1 Quantization of Analog Signals and Calculation of Quantized noise Objective:
More information2.1 BASIC CONCEPTS Basic Operations on Signals Time Shifting. Figure 2.2 Time shifting of a signal. Time Reversal.
1 2.1 BASIC CONCEPTS 2.1.1 Basic Operations on Signals Time Shifting. Figure 2.2 Time shifting of a signal. Time Reversal. 2 Time Scaling. Figure 2.4 Time scaling of a signal. 2.1.2 Classification of Signals
More informationSignals and Systems Lecture 6: Fourier Applications
Signals and Systems Lecture 6: Fourier Applications Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Winter 2012 arzaneh Abdollahi Signal and Systems Lecture 6
More informationBasic Signals and Systems
Chapter 2 Basic Signals and Systems A large part of this chapter is taken from: C.S. Burrus, J.H. McClellan, A.V. Oppenheim, T.W. Parks, R.W. Schafer, and H. W. Schüssler: Computer-based exercises for
More informationHW 1 is due on tuesday. PPI is due on Thurs ( to hero by 5PM) Lab starts next week.
EECS 452 Lecture 2 Today: Sampling and reconstruction review FIR and IIR filters C5515 ezdsp Direct digital synthesis Reminders: HW 1 is due on tuesday. PPI is due on Thurs (email to hero by 5PM) Lab starts
More informationPROBLEM SET 5. Reminder: Quiz 1will be on March 6, during the regular class hour. Details to follow. z = e jω h[n] H(e jω ) H(z) DTFT.
PROBLEM SET 5 Issued: 2/4/9 Due: 2/22/9 Reading: During the past week we continued our discussion of the impact of pole/zero locations on frequency response, focusing on allpass systems, minimum and maximum-phase
More informationThe University of Texas at Austin Dept. of Electrical and Computer Engineering Final Exam
The University of Texas at Austin Dept. of Electrical and Computer Engineering Final Exam Date: December 18, 2017 Course: EE 313 Evans Name: Last, First The exam is scheduled to last three hours. Open
More informationTABLE OF CONTENTS TOPIC NUMBER NAME OF THE TOPIC 1. OVERVIEW OF SIGNALS & SYSTEMS 2. ANALYSIS OF LTI SYSTEMS- Z TRANSFORM 3. ANALYSIS OF FT, DFT AND FFT SIGNALS 4. DIGITAL FILTERS CONCEPTS & DESIGN 5.
More informationDigital Filters IIR (& Their Corresponding Analog Filters) Week Date Lecture Title
http://elec3004.com Digital Filters IIR (& Their Corresponding Analog Filters) 2017 School of Information Technology and Electrical Engineering at The University of Queensland Lecture Schedule: Week Date
More informationCMPT 318: Lecture 4 Fundamentals of Digital Audio, Discrete-Time Signals
CMPT 318: Lecture 4 Fundamentals of Digital Audio, Discrete-Time Signals Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 16, 2006 1 Continuous vs. Discrete
More informationContinuous vs. Discrete signals. Sampling. Analog to Digital Conversion. CMPT 368: Lecture 4 Fundamentals of Digital Audio, Discrete-Time Signals
Continuous vs. Discrete signals CMPT 368: Lecture 4 Fundamentals of Digital Audio, Discrete-Time Signals Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 22,
More informationIt is the speed and discrete nature of the FFT that allows us to analyze a signal's spectrum with MATLAB.
MATLAB Addendum on Fourier Stuff 1. Getting to know the FFT What is the FFT? FFT = Fast Fourier Transform. The FFT is a faster version of the Discrete Fourier Transform(DFT). The FFT utilizes some clever
More informationDigital Video and Audio Processing. Winter term 2002/ 2003 Computer-based exercises
Digital Video and Audio Processing Winter term 2002/ 2003 Computer-based exercises Rudolf Mester Institut für Angewandte Physik Johann Wolfgang Goethe-Universität Frankfurt am Main 6th November 2002 Chapter
More informationFinal Exam. EE313 Signals and Systems. Fall 1999, Prof. Brian L. Evans, Unique No
Final Exam EE313 Signals and Systems Fall 1999, Prof. Brian L. Evans, Unique No. 14510 December 11, 1999 The exam is scheduled to last 50 minutes. Open books and open notes. You may refer to your homework
More informationObjectives. Presentation Outline. Digital Modulation Lecture 03
Digital Modulation Lecture 03 Inter-Symbol Interference Power Spectral Density Richard Harris Objectives To be able to discuss Inter-Symbol Interference (ISI), its causes and possible remedies. To be able
More informationExperiment 4- Finite Impulse Response Filters
Experiment 4- Finite Impulse Response Filters 18 February 2009 Abstract In this experiment we design different Finite Impulse Response filters and study their characteristics. 1 Introduction The transfer
More informationFinal Exam Solutions June 14, 2006
Name or 6-Digit Code: PSU Student ID Number: Final Exam Solutions June 14, 2006 ECE 223: Signals & Systems II Dr. McNames Keep your exam flat during the entire exam. If you have to leave the exam temporarily,
More informationBasics of Digital Filtering
4 Basics of Digital Filtering Willis J. Tompkins and Pradeep Tagare In this chapter we introduce the concept of digital filtering and look at the advantages, disadvantages, and differences between analog
More informationDSP Laboratory (EELE 4110) Lab#10 Finite Impulse Response (FIR) Filters
Islamic University of Gaza OBJECTIVES: Faculty of Engineering Electrical Engineering Department Spring-2011 DSP Laboratory (EELE 4110) Lab#10 Finite Impulse Response (FIR) Filters To demonstrate the concept
More informationDigital Filters. Linearity and Time Invariance. Implications of Linear Time Invariance (LTI) Music 270a: Introduction to Digital Filters
Digital Filters Music 7a: Introduction to Digital Filters Tamara Smyth, trsmyth@ucsd.edu Department of Music, University of California, San Diego (UCSD) November 7, 7 Any medium through which a signal
More informationEEM478-DSPHARDWARE. WEEK12:FIR & IIR Filter Design
EEM478-DSPHARDWARE WEEK12:FIR & IIR Filter Design PART-I : Filter Design/Realization Step-1 : define filter specs (pass-band, stop-band, optimization criterion, ) Step-2 : derive optimal transfer function
More informationEE 311 February 13 and 15, 2019 Lecture 10
EE 311 February 13 and 15, 219 Lecture 1 Figure 4.22 The top figure shows a quantized sinusoid as the darker stair stepped curve. The bottom figure shows the quantization error. The quantized signal to
More informationEE123 Digital Signal Processing
EE123 Digital Signal Processing Lecture 5A Time-Frequency Tiling Subtleties in filtering/processing with DFT x[n] H(e j! ) y[n] System is implemented by overlap-and-save Filtering using DFT H[k] π 2π Subtleties
More informationInfinite Impulse Response (IIR) Filter. Ikhwannul Kholis, ST., MT. Universitas 17 Agustus 1945 Jakarta
Infinite Impulse Response (IIR) Filter Ihwannul Kholis, ST., MT. Universitas 17 Agustus 1945 Jaarta The Outline 8.1 State-of-the-art 8.2 Coefficient Calculation Method for IIR Filter 8.2.1 Pole-Zero Placement
More informationDFT: Discrete Fourier Transform & Linear Signal Processing
DFT: Discrete Fourier Transform & Linear Signal Processing 2 nd Year Electronics Lab IMPERIAL COLLEGE LONDON Table of Contents Equipment... 2 Aims... 2 Objectives... 2 Recommended Textbooks... 3 Recommended
More informationNH 67, Karur Trichy Highways, Puliyur C.F, Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3
NH 67, Karur Trichy Highways, Puliyur C.F, 639 114 Karur District DEPARTMENT OF INFORMATION TECHNOLOGY DIGITAL SIGNAL PROCESSING UNIT 3 IIR FILTER DESIGN Structure of IIR System design of Discrete time
More informationMultirate Signal Processing Lecture 7, Sampling Gerald Schuller, TU Ilmenau
Multirate Signal Processing Lecture 7, Sampling Gerald Schuller, TU Ilmenau (Also see: Lecture ADSP, Slides 06) In discrete, digital signal we use the normalized frequency, T = / f s =: it is without a
More informationChapter-2 SAMPLING PROCESS
Chapter-2 SAMPLING PROCESS SAMPLING: A message signal may originate from a digital or analog source. If the message signal is analog in nature, then it has to be converted into digital form before it can
More information+ a(t) exp( 2πif t)dt (1.1) In order to go back to the independent variable t, we define the inverse transform as: + A(f) exp(2πif t)df (1.
Chapter Fourier analysis In this chapter we review some basic results from signal analysis and processing. We shall not go into detail and assume the reader has some basic background in signal analysis
More informationDISCRETE FOURIER TRANSFORM AND FILTER DESIGN
DISCRETE FOURIER TRANSFORM AND FILTER DESIGN N. C. State University CSC557 Multimedia Computing and Networking Fall 2001 Lecture # 03 Spectrum of a Square Wave 2 Results of Some Filters 3 Notation 4 x[n]
More informationSignals and Filtering
FILTERING OBJECTIVES The objectives of this lecture are to: Introduce signal filtering concepts Introduce filter performance criteria Introduce Finite Impulse Response (FIR) filters Introduce Infinite
More informationSampling and Reconstruction of Analog Signals
Sampling and Reconstruction of Analog Signals Chapter Intended Learning Outcomes: (i) Ability to convert an analog signal to a discrete-time sequence via sampling (ii) Ability to construct an analog signal
More informationTheory of Telecommunications Networks
Theory of Telecommunications Networks Anton Čižmár Ján Papaj Department of electronics and multimedia telecommunications CONTENTS Preface... 5 1 Introduction... 6 1.1 Mathematical models for communication
More informationDigital Filters FIR and IIR Systems
Digital Filters FIR and IIR Systems ELEC 3004: Systems: Signals & Controls Dr. Surya Singh (Some material adapted from courses by Russ Tedrake and Elena Punskaya) Lecture 16 elec3004@itee.uq.edu.au http://robotics.itee.uq.edu.au/~elec3004/
More informationMoving from continuous- to discrete-time
Moving from continuous- to discrete-time Sampling ideas Uniform, periodic sampling rate, e.g. CDs at 44.1KHz First we will need to consider periodic signals in order to appreciate how to interpret discrete-time
More informationFinal Exam Practice Questions for Music 421, with Solutions
Final Exam Practice Questions for Music 4, with Solutions Elementary Fourier Relationships. For the window w = [/,,/ ], what is (a) the dc magnitude of the window transform? + (b) the magnitude at half
More informationFrom Fourier Series to Analysis of Non-stationary Signals - VII
From Fourier Series to Analysis of Non-stationary Signals - VII prof. Miroslav Vlcek November 23, 2010 Contents Short Time Fourier Transform 1 Short Time Fourier Transform 2 Contents Short Time Fourier
More informationELECTRONOTES APPLICATION NOTE NO Hanshaw Road Ithaca, NY Nov 7, 2014 MORE CONCERNING NON-FLAT RANDOM FFT
ELECTRONOTES APPLICATION NOTE NO. 416 1016 Hanshaw Road Ithaca, NY 14850 Nov 7, 2014 MORE CONCERNING NON-FLAT RANDOM FFT INTRODUCTION A curiosity that has probably long been peripherally noted but which
More informationDigital Signal Processing
Digital Signal Processing System Analysis and Design Paulo S. R. Diniz Eduardo A. B. da Silva and Sergio L. Netto Federal University of Rio de Janeiro CAMBRIDGE UNIVERSITY PRESS Preface page xv Introduction
More informationECE503: Digital Filter Design Lecture 9
ECE503: Digital Filter Design Lecture 9 D. Richard Brown III WPI 26-March-2012 WPI D. Richard Brown III 26-March-2012 1 / 33 Lecture 9 Topics Within the broad topic of digital filter design, we are going
More informationECE 201: Introduction to Signal Analysis
ECE 201: Introduction to Signal Analysis Prof. Paris Last updated: October 9, 2007 Part I Spectrum Representation of Signals Lecture: Sums of Sinusoids (of different frequency) Introduction Sum of Sinusoidal
More informationAN2182 Application note
Application note Filters using the ST10 DSP library Introduction The ST10F2xx family provides a 16-bit multiply and accumulate unit (MAC) allowing control-oriented signal processing and filtering widely
More informationProject I: Phase Tracking and Baud Timing Correction Systems
Project I: Phase Tracking and Baud Timing Correction Systems ECES 631, Prof. John MacLaren Walsh, Ph. D. 1 Purpose In this lab you will encounter the utility of the fundamental Fourier and z-transform
More informationSubtractive Synthesis. Describing a Filter. Filters. CMPT 468: Subtractive Synthesis
Subtractive Synthesis CMPT 468: Subtractive Synthesis Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University November, 23 Additive synthesis involves building the sound by
More informationDIGITAL FILTERS. !! Finite Impulse Response (FIR) !! Infinite Impulse Response (IIR) !! Background. !! Matlab functions AGC DSP AGC DSP
DIGITAL FILTERS!! Finite Impulse Response (FIR)!! Infinite Impulse Response (IIR)!! Background!! Matlab functions 1!! Only the magnitude approximation problem!! Four basic types of ideal filters with magnitude
More information4. Design of Discrete-Time Filters
4. Design of Discrete-Time Filters 4.1. Introduction (7.0) 4.2. Frame of Design of IIR Filters (7.1) 4.3. Design of IIR Filters by Impulse Invariance (7.1) 4.4. Design of IIR Filters by Bilinear Transformation
More informationGEORGIA INSTITUTE OF TECHNOLOGY. SCHOOL of ELECTRICAL and COMPUTER ENGINEERING. ECE 2026 Summer 2018 Lab #8: Filter Design of FIR Filters
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2026 Summer 2018 Lab #8: Filter Design of FIR Filters Date: 19. Jul 2018 Pre-Lab: You should read the Pre-Lab section of
More informationDigital Filters IIR (& Their Corresponding Analog Filters) 4 April 2017 ELEC 3004: Systems 1. Week Date Lecture Title
http://elec3004.com Digital Filters IIR (& Their Corresponding Analog Filters) 4 April 017 ELEC 3004: Systems 1 017 School of Information Technology and Electrical Engineering at The University of Queensland
More informationFilters. Phani Chavali
Filters Phani Chavali Filters Filtering is the most common signal processing procedure. Used as echo cancellers, equalizers, front end processing in RF receivers Used for modifying input signals by passing
More informationEC6502 PRINCIPLES OF DIGITAL SIGNAL PROCESSING
1. State the properties of DFT? UNIT-I DISCRETE FOURIER TRANSFORM 1) Periodicity 2) Linearity and symmetry 3) Multiplication of two DFTs 4) Circular convolution 5) Time reversal 6) Circular time shift
More informationDigital Filtering: Realization
Digital Filtering: Realization Digital Filtering: Matlab Implementation: 3-tap (2 nd order) IIR filter 1 Transfer Function Differential Equation: z- Transform: Transfer Function: 2 Example: Transfer Function
More informationELECTRONICS AND COMMUNICATION ENGINEERING BOOLEAN ALGEBRA THE Z-TRANSFORM THE SUPERHETERODYNE RECEIVER
GATE - ELECTRONICS AND COMMUNICATION ENGINEERING SAMPLE THEORY BOOLEAN ALGEBRA THE Z-TRANSFORM THE SUPERHETERODYNE RECEIVER For IIT-JAM, JNU, GATE, NET, NIMCET and Other Entrance Exams -C-8, Sheela Chowdhary
More informationFourier Methods of Spectral Estimation
Department of Electrical Engineering IIT Madras Outline Definition of Power Spectrum Deterministic signal example Power Spectrum of a Random Process The Periodogram Estimator The Averaged Periodogram Blackman-Tukey
More informationBrief Introduction to Signals & Systems. Phani Chavali
Brief Introduction to Signals & Systems Phani Chavali Outline Signals & Systems Continuous and discrete time signals Properties of Systems Input- Output relation : Convolution Frequency domain representation
More informationLab 4 Digital Scope and Spectrum Analyzer
Lab 4 Digital Scope and Spectrum Analyzer Page 4.1 Lab 4 Digital Scope and Spectrum Analyzer Goals Review Starter files Interface a microphone and record sounds, Design and implement an analog HPF, LPF
More informationVolume 3 Signal Processing Reference Manual
Contents Volume 3 Signal Processing Reference Manual Contents 1 Sampling analogue signals 1.1 Introduction...1-1 1.2 Selecting a sampling speed...1-1 1.3 References...1-5 2 Digital filters 2.1 Introduction...2-1
More informationFundamental Frequency Detection
Fundamental Frequency Detection Jan Černocký, Valentina Hubeika {cernocky ihubeika}@fit.vutbr.cz DCGM FIT BUT Brno Fundamental Frequency Detection Jan Černocký, Valentina Hubeika, DCGM FIT BUT Brno 1/37
More informationSignal Processing. Introduction
Signal Processing 0 Introduction One of the premiere uses of MATLAB is in the analysis of signal processing and control systems. In this chapter we consider signal processing. The final chapter of the
More informationCorso di DATI e SEGNALI BIOMEDICI 1. Carmelina Ruggiero Laboratorio MedInfo
Corso di DATI e SEGNALI BIOMEDICI 1 Carmelina Ruggiero Laboratorio MedInfo Digital Filters Function of a Filter In signal processing, the functions of a filter are: to remove unwanted parts of the signal,
More informationModule 3 : Sampling and Reconstruction Problem Set 3
Module 3 : Sampling and Reconstruction Problem Set 3 Problem 1 Shown in figure below is a system in which the sampling signal is an impulse train with alternating sign. The sampling signal p(t), the Fourier
More informationDigital Signal Processing
Digital Signal Processing Lecture 9 Discrete-Time Processing of Continuous-Time Signals Alp Ertürk alp.erturk@kocaeli.edu.tr Analog to Digital Conversion Most real life signals are analog signals These
More informationESE531 Spring University of Pennsylvania Department of Electrical and System Engineering Digital Signal Processing
University of Pennsylvania Department of Electrical and System Engineering Digital Signal Processing ESE531, Spring 2017 Final Project: Audio Equalization Wednesday, Apr. 5 Due: Tuesday, April 25th, 11:59pm
More informationMusic 270a: Fundamentals of Digital Audio and Discrete-Time Signals
Music 270a: Fundamentals of Digital Audio and Discrete-Time Signals Tamara Smyth, trsmyth@ucsd.edu Department of Music, University of California, San Diego October 3, 2016 1 Continuous vs. Discrete signals
More informationPYKC 27 Feb 2017 EA2.3 Electronics 2 Lecture PYKC 27 Feb 2017 EA2.3 Electronics 2 Lecture 11-2
In this lecture, I will introduce the mathematical model for discrete time signals as sequence of samples. You will also take a first look at a useful alternative representation of discrete signals known
More information