Digital Signal Processing
|
|
- Alexandra Rogers
- 5 years ago
- Views:
Transcription
1 Digital Signal Processing Lab 1: FFT, Spectral Leakage, Zero Padding Moslem Amiri, Václav Přenosil Embedded Systems Laboratory Faculty of Informatics, Masaryk University Brno, Czech Republic Spring, 2015
2 Open MATLAB Create an m-file An m-file, or script file, is a simple text file where you can place MATLAB commands Choose New from the File menu and select Script Generate two separate 64 length buffers of the two sinusoids: Use a sampling frequency of 1 khz x 1 (t) = 100 sin(2π101.56t) x 2 (t) = 2 sin(2π156.25t) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
3 y = 100*sin(2*pi*101.56*t); figure(1);plot(fs*t(1:50),y(1:50)); title( Signal ); xlabel( time (milliseconds) ); Run the m-file Press F5 Or type filename at the MATLAB command window prompt Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
4 y = 2*sin(2*pi*156.25*t); figure(2);plot(fs*t(1:50),y(1:50)); title( Signal ); xlabel( time (milliseconds) ); Execute the commands Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
5 Evaluate the FFT of each of these signals using a rectangular window, and determine which frequency bins have the highest peaks Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
6 y = 100*sin(2*pi*101.56*t); Y = fft(y,l)/l; f = Fs/2*linspace(0,1,L/2+1); % Plot single-sided amplitude spectrum. figure(1);plot(f,2*abs(y(1:l/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
7 We should divide FFT magnitudes by N/2 for real inputs The division by L is done in Y = fft(y,l)/l; Multiplication by 2 is done in plot(f,2*abs(y(1:l/2+1))) Command y = linspace(a,b,n) Generates a row vector y of n points linearly spaced between and including a and b For n < 2, linspace returns b f = Fs/2*linspace(0,1,L/2+1); Frequency of index bins is mfs/n For L = 64, it must be Fs, i.e. one full interval in frequency domain, so for half of that (L/2+1), it is Fs/2 Execute the commands Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
8 y = 2*sin(2*pi*156.25*t); Y = fft(y,l)/l; f = Fs/2*linspace(0,1,L/2+1); figure(4);plot(f,2*abs(y(1:l/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
9 Repeat the FFT evaluations done previously, but using a 1024 point FFT Applying a 1024 point FFT to a 64 length data buffer has the effect of appending ( ) zeros to the end of the data, which increases the number of points at which the spectrum is evaluated Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
10 y = 100*sin(2*pi*101.56*t); Y = fft(y,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(1);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
11 Y = fft(x,n) returns the n-point DFT fft(x) is equivalent to fft(x,n) where n is the size of X in the first nonsingleton dimension If the length of X is less than n, X is padded with trailing zeros to length n If the length of X is greater than n, the sequence X is truncated When X is a matrix, the length of the columns are adjusted in the same manner Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
12 y = 2*sin(2*pi*156.25*t); Y = fft(y,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(4);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
13 The Hamming window is defined as: w(n) = cos(2πn/n) for n = 0, 1, 2,..., N 1 Generate and plot the Hamming window Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
14 n = (0:63); w = * cos(2 * pi * n / 64); plot(n,w); Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
15 Now repeat FFT evaluations (without and with zero padding), but applying first a Hamming window to the 64 length data buffers before evaluating the FFTs Observe the outputs of the FFTs Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
16 y = 100*sin(2*pi*101.56*t); w = * cos(2 * pi * (0 : 63) / 64); win = y.* w; Y = fft(win,l)/l; f = Fs/2*linspace(0,1,L/2+1); figure(2);plot(f,2*abs(y(1:l/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
17 y = 100*sin(2*pi*101.56*t); w = * cos(2 * pi * (0 : 63) / 64); win = y.* w; Y = fft(win,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(3);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
18 y = 2*sin(2*pi*156.25*t); w = * cos(2 * pi * (0 : 63) / 64); win = y.* w; Y = fft(win,l)/l; f = Fs/2*linspace(0,1,L/2+1); figure(1);plot(f,2*abs(y(1:l/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
19 y = 2*sin(2*pi*156.25*t); w = * cos(2 * pi * (0 : 63) / 64); win = y.* w; Y = fft(win,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(1);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
20 Now add these two sine waves together, and apply both 64-length and 1024-length FFTs, with both rectangular and Hamming windows Observe how spectral leakage can mask a small signal in a large one Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
21 y = 100*sin(2*pi*101.56*t) + 2*sin(2*pi*156.25*t); Y = fft(y,l)/l; f = Fs/2*linspace(0,1,L/2+1); figure(2);plot(f,2*abs(y(1:l/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
22 y = 100*sin(2*pi*101.56*t) + 2*sin(2*pi*156.25*t); w = * cos(2 * pi * (0 : 63) / 64); win = y.* w; Y = fft(win,l)/l; f = Fs/2*linspace(0,1,L/2+1); figure(2);plot(f,2*abs(y(1:l/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
23 y = 100*sin(2*pi*101.56*t) + 2*sin(2*pi*156.25*t); Y = fft(y,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(1);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
24 y = 100*sin(2*pi*101.56*t) + 2*sin(2*pi*156.25*t); w = * cos(2 * pi * (0 : 63) / 64); win = y.* w; Y = fft(win,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(3);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
25 Now add Gaussian noise to the combined signals, with a variance of 10 Compare the results of evaluating the 1024 point FFT on a 64-length data buffer (i.e. with zero padding), with a 1024 point FFT applied to a 1024-length data buffer (i.e. generate more data) In MATLAB, to generate the noise, try: noise = sqrt (variance) * randn(); Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
26 x = 100*sin(2*pi*101.56*t) + 2*sin(2*pi*156.25*t); y = x + sqrt(10)*randn(size(t)); % Sinusoids plus noise Y = fft(y,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(1);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
27 L = 1024; x = 100*sin(2*pi*101.56*t) + 2*sin(2*pi*156.25*t); y = x + sqrt(10)*randn(size(t)); % Sinusoids plus noise Y = fft(y,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(2);plot(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
28 Logarithmic axis scaling Log-log and semi-log plots are created with commands that act just like the plot command Command Name: loglog, Plot type: log(y) versus log(x) Command Name: semilogx, Plot type: y versus log(x) Command Name: semilogy, Plot type: log(y) versus x Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
29 L = 1024; x = 100*sin(2*pi*101.56*t) + 2*sin(2*pi*156.25*t); y = x + sqrt(10)*randn(size(t)); % Sinusoids plus noise Y = fft(y,1024)/l; f = Fs/2*linspace(0,1,1024/2+1); figure(2);semilogy(f,2*abs(y(1:1024/2+1))) Moslem Amiri, Václav Přenosil Digital Signal Processing Spring, / 29
Discrete Fourier Transform (DFT)
Amplitude Amplitude Discrete Fourier Transform (DFT) DFT transforms the time domain signal samples to the frequency domain components. DFT Signal Spectrum Time Frequency DFT is often used to do frequency
More informationL A B 3 : G E N E R A T I N G S I N U S O I D S
L A B 3 : G E N E R A T I N G S I N U S O I D S NAME: DATE OF EXPERIMENT: DATE REPORT SUBMITTED: 1/7 1 THEORY DIGITAL SIGNAL PROCESSING LABORATORY 1.1 GENERATION OF DISCRETE TIME SINUSOIDAL SIGNALS IN
More informationWireless Communication Systems Laboratory Lab#1: An introduction to basic digital baseband communication through MATLAB simulation Objective
Wireless Communication Systems Laboratory Lab#1: An introduction to basic digital baseband communication through MATLAB simulation Objective The objective is to teach students a basic digital communication
More informationLABORATORY - FREQUENCY ANALYSIS OF DISCRETE-TIME SIGNALS
LABORATORY - FREQUENCY ANALYSIS OF DISCRETE-TIME SIGNALS INTRODUCTION The objective of this lab is to explore many issues involved in sampling and reconstructing signals, including analysis of the frequency
More informationMemorial University of Newfoundland Faculty of Engineering and Applied Science. Lab Manual
Memorial University of Newfoundland Faculty of Engineering and Applied Science Engineering 6871 Communication Principles Lab Manual Fall 2014 Lab 1 AMPLITUDE MODULATION Purpose: 1. Learn how to use Matlab
More informationLog Booklet for EE2 Experiments
Log Booklet for EE2 Experiments Vasil Zlatanov DFT experiment Exercise 1 Code for sinegen.m function y = sinegen(fsamp, fsig, nsamp) tsamp = 1/fsamp; t = 0 : tsamp : (nsamp-1)*tsamp; y = sin(2*pi*fsig*t);
More informationELT COMMUNICATION THEORY
ELT 41307 COMMUNICATION THEORY Matlab Exercise #1 Sampling, Fourier transform, Spectral illustrations, and Linear filtering 1 SAMPLING The modeled signals and systems in this course are mostly analog (continuous
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 informationProblem Set 1 (Solutions are due Mon )
ECEN 242 Wireless Electronics for Communication Spring 212 1-23-12 P. Mathys Problem Set 1 (Solutions are due Mon. 1-3-12) 1 Introduction The goals of this problem set are to use Matlab to generate and
More informationMaking 2D Plots in Matlab
Making 2D Plots in Matlab Gerald W. Recktenwald Department of Mechanical Engineering Portland State University gerry@pdx.edu ME 350: Plotting with Matlab Overview Plotting in Matlab Plotting (x, y) data
More informationChapter 5 Window Functions. periodic with a period of N (number of samples). This is observed in table (3.1).
Chapter 5 Window Functions 5.1 Introduction As discussed in section (3.7.5), the DTFS assumes that the input waveform is periodic with a period of N (number of samples). This is observed in table (3.1).
More informationFrequency Domain Representation of Signals
Frequency Domain Representation of Signals The Discrete Fourier Transform (DFT) of a sampled time domain waveform x n x 0, x 1,..., x 1 is a set of Fourier Coefficients whose samples are 1 n0 X k X0, X
More informationSIGNALS AND SYSTEMS LABORATORY 3: Construction of Signals in MATLAB
SIGNALS AND SYSTEMS LABORATORY 3: Construction of Signals in MATLAB INTRODUCTION Signals are functions of time, denoted x(t). For simulation, with computers and digital signal processing hardware, one
More informationEE 215 Semester Project SPECTRAL ANALYSIS USING FOURIER TRANSFORM
EE 215 Semester Project SPECTRAL ANALYSIS USING FOURIER TRANSFORM Department of Electrical and Computer Engineering Missouri University of Science and Technology Page 1 Table of Contents Introduction...Page
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 informationTHE HONG KONG POLYTECHNIC UNIVERSITY Department of Electronic and Information Engineering. EIE2106 Signal and System Analysis Lab 2 Fourier series
THE HONG KONG POLYTECHNIC UNIVERSITY Department of Electronic and Information Engineering EIE2106 Signal and System Analysis Lab 2 Fourier series 1. Objective The goal of this laboratory exercise is to
More informationDistortion Analysis T S. 2 N for all k not defined above. THEOREM?: If N P is an integer and x(t) is band limited to f MAX, then
EE 505 Lecture 6 Spectral Analysis in Spectre - Standard transient analysis - Strobe period transient analysis Addressing Spectral Analysis Challenges Problem Awareness Windowing Post-processing . Review
More informationIntegrators, differentiators, and simple filters
BEE 233 Laboratory-4 Integrators, differentiators, and simple filters 1. Objectives Analyze and measure characteristics of circuits built with opamps. Design and test circuits with opamps. Plot gain vs.
More informationFourier Signal Analysis
Part 1B Experimental Engineering Integrated Coursework Location: Baker Building South Wing Mechanics Lab Experiment A4 Signal Processing Fourier Signal Analysis Please bring the lab sheet from 1A experiment
More informationINTRODUCTION TO MATLAB by. Introduction to Matlab
INTRODUCTION TO MATLAB by Mohamed Hussein Lecture 5 Introduction to Matlab More on XY Plotting Other Types of Plotting 3D Plot (XYZ Plotting) More on XY Plotting Other XY plotting commands are axis ([xmin
More informationELEC3104: Digital Signal Processing Session 1, 2013
ELEC3104: Digital Signal Processing Session 1, 2013 The University of New South Wales School of Electrical Engineering and Telecommunications LABORATORY 1: INTRODUCTION TO TIMS AND MATLAB INTRODUCTION
More informationNoise Measurements Using a Teledyne LeCroy Oscilloscope
Noise Measurements Using a Teledyne LeCroy Oscilloscope TECHNICAL BRIEF January 9, 2013 Summary Random noise arises from every electronic component comprising your circuits. The analysis of random electrical
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 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 informationAttia, John Okyere. Plotting Commands. Electronics and Circuit Analysis using MATLAB. Ed. John Okyere Attia Boca Raton: CRC Press LLC, 1999
Attia, John Okyere. Plotting Commands. Electronics and Circuit Analysis using MATLAB. Ed. John Okyere Attia Boca Raton: CRC Press LLC, 1999 1999 by CRC PRESS LLC CHAPTER TWO PLOTTING COMMANDS 2.1 GRAPH
More informationExperiment 1 Introduction to MATLAB and Simulink
Experiment 1 Introduction to MATLAB and Simulink INTRODUCTION MATLAB s Simulink is a powerful modeling tool capable of simulating complex digital communications systems under realistic conditions. It includes
More informationEEL 4350 Principles of Communication Project 2 Due Tuesday, February 10 at the Beginning of Class
EEL 4350 Principles of Communication Project 2 Due Tuesday, February 10 at the Beginning of Class Description In this project, MATLAB and Simulink are used to construct a system experiment. The experiment
More informationLaboratory Assignment 4. Fourier Sound Synthesis
Laboratory Assignment 4 Fourier Sound Synthesis PURPOSE This lab investigates how to use a computer to evaluate the Fourier series for periodic signals and to synthesize audio signals from Fourier series
More information9.1. Probability and Statistics
9. Probability and Statistics Measured signals exhibit deterministic (predictable) and random (unpredictable) behavior. The deterministic behavior is often governed by a differential equation, while the
More informationDESIGN OF A SIMPLE RELIABLE VOTER FOR MODULAR REDUNDANCY IMPLEMENTATIONS
DESIGN OF A SIMPLE RELIABLE VOTER FOR MODULAR REDUNDANCY IMPLEMENTATIONS Moslem Amiri, Václav Přenosil Faculty of Informatics, Masaryk University Brno, Czech Republic, amiri@mail.muni.cz, prenosil@fi.muni.cz
More informationAssignment 6: Solution to MATLAB code for BER generation of QPSK system over AWGN channel.
G. S. Sanyal School of Telecommunications Indian Institute of Technology Kharagpur MOOC: Spread Spectrum Communications & Jamming Assignment 6: Solution to MATLAB code for BER generation of QPSK system
More informationCSCD 409 Scientific Programming. Module 6: Plotting (Chpt 5)
CSCD 409 Scientific Programming Module 6: Plotting (Chpt 5) 2008-2012, Prentice Hall, Paul Schimpf All rights reserved. No portion of this presentation may be reproduced, in whole or in part, in any form
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 informationEE 464 Short-Time Fourier Transform Fall and Spectrogram. Many signals of importance have spectral content that
EE 464 Short-Time Fourier Transform Fall 2018 Read Text, Chapter 4.9. and Spectrogram Many signals of importance have spectral content that changes with time. Let xx(nn), nn = 0, 1,, NN 1 1 be a discrete-time
More informationDiscrete Fourier Transform, DFT Input: N time samples
EE445M/EE38L.6 Lecture. Lecture objectives are to: The Discrete Fourier Transform Windowing Use DFT to design a FIR digital filter Discrete Fourier Transform, DFT Input: time samples {a n = {a,a,a 2,,a
More informationOperational Amplifier Circuits
ECE VIII. Basic 5 Operational Amplifier Circuits Lab 8 In this lab we will verify the operation of inverting and noninverting amplifiers constructed using Operational Amplifiers. We will also observe the
More informationECE 2713 Homework 7 DUE: 05/1/2018, 11:59 PM
Spring 2018 What to Turn In: ECE 2713 Homework 7 DUE: 05/1/2018, 11:59 PM Dr. Havlicek Submit your solution for this assignment electronically on Canvas by uploading a file to ECE-2713-001 > Assignments
More informationENSC327 Communication Systems Fall 2011 Assignment #1 Due Wednesday, Sept. 28, 4:00 pm
ENSC327 Communication Systems Fall 2011 Assignment #1 Due Wednesday, Sept. 28, 4:00 pm All problem numbers below refer to those in Haykin & Moher s book. 1. (FT) Problem 2.20. 2. (Convolution) Problem
More informationEE 435. Lecture 34. Spectral Performance Windowing Quantization Noise
EE 435 Lecture 34 Spectral Performance Windowing Quantization Noise . Review from last lecture. Are there any strategies to address the problem of requiring precisely an integral number of periods to use
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 informationLab #2 First Order RC Circuits Week of 27 January 2015
ECE214: Electrical Circuits Laboratory Lab #2 First Order RC Circuits Week of 27 January 2015 1 Introduction In this lab you will investigate the magnitude and phase shift that occurs in an RC circuit
More informationUse Matlab Function pwelch to Find Power Spectral Density or Do It Yourself
Use Matlab Function pwelch to Find Power Spectral Density or Do It Yourself In my last post, we saw that finding the spectrum of a signal requires several steps beyond computing the discrete Fourier transform
More informationFall Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class
Fall 2018 2019 Music 320A Homework #2 Sinusoids, Complex Sinusoids 145 points Theory and Lab Problems Due Thursday 10/11/2018 before class Theory Problems 1. 15 pts) [Sinusoids] Define xt) as xt) = 2sin
More informationECE 5650/4650 MATLAB Project 1
This project is to be treated as a take-home exam, meaning each student is to due his/her own work. The project due date is 4:30 PM Tuesday, October 18, 2011. To work the project you will need access to
More informationSound synthesis with Pure Data
Sound synthesis with Pure Data 1. Start Pure Data from the programs menu in classroom TC307. You should get the following window: The DSP check box switches sound output on and off. Getting sound out First,
More informationEXPERIMENT 4 INTRODUCTION TO AMPLITUDE MODULATION SUBMITTED BY
EXPERIMENT 4 INTRODUCTION TO AMPLITUDE MODULATION SUBMITTED BY NAME:. STUDENT ID:.. ROOM: INTRODUCTION TO AMPLITUDE MODULATION Purpose: The objectives of this laboratory are:. To introduce the spectrum
More informationECEGR Lab #8: Introduction to Simulink
Page 1 ECEGR 317 - Lab #8: Introduction to Simulink Objective: By: Joe McMichael This lab is an introduction to Simulink. The student will become familiar with the Help menu, go through a short example,
More information(i) Understanding the basic concepts of signal modeling, correlation, maximum likelihood estimation, least squares and iterative numerical methods
Tools and Applications Chapter Intended Learning Outcomes: (i) Understanding the basic concepts of signal modeling, correlation, maximum likelihood estimation, least squares and iterative numerical methods
More informationReal Analog - Circuits 1 Chapter 11: Lab Projects
Real Analog - Circuits 1 Chapter 11: Lab Projects 11.2.1: Signals with Multiple Frequency Components Overview: In this lab project, we will calculate the magnitude response of an electrical circuit and
More informationLab 0: Introduction to TIMS AND MATLAB
TELE3013 TELECOMMUNICATION SYSTEMS 1 Lab 0: Introduction to TIMS AND MATLAB 1. INTRODUCTION The TIMS (Telecommunication Instructional Modelling System) system was first developed by Tim Hooper, then a
More informationEECS 452 Midterm Exam (solns) Fall 2012
EECS 452 Midterm Exam (solns) Fall 2012 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 I /40 Section
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 informationDSP First. Laboratory Exercise #11. Extracting Frequencies of Musical Tones
DSP First Laboratory Exercise #11 Extracting Frequencies of Musical Tones This lab is built around a single project that involves the implementation of a system for automatically writing a musical score
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 informationChapter Three. The Discrete Fourier Transform
Chapter Three. The Discrete Fourier Transform The discrete Fourier transform (DFT) is one of the two most common, and powerful, procedures encountered in the field of digital signal processing. (Digital
More informationComputer Programming ECIV 2303 Chapter 5 Two-Dimensional Plots Instructor: Dr. Talal Skaik Islamic University of Gaza Faculty of Engineering
Computer Programming ECIV 2303 Chapter 5 Two-Dimensional Plots Instructor: Dr. Talal Skaik Islamic University of Gaza Faculty of Engineering 1 Introduction Plots are a very useful tool for presenting information.
More informationUNIVERSITY OF WARWICK
UNIVERSITY OF WARWICK School of Engineering ES905 MSc Signal Processing Module (2010) AM SIGNALS AND FILTERING EXERCISE Deadline: This is NOT for credit. It is best done before the first assignment. You
More informationUNIVERSITY OF WARWICK
UNIVERSITY OF WARWICK School of Engineering ES905 MSc Signal Processing Module (2004) ASSIGNMENT 1 In this assignment, you will use the MATLAB package. In Part (A) you will design some FIR filters and
More informationGeorge Mason University Signals and Systems I Spring 2016
George Mason University Signals and Systems I Spring 2016 Laboratory Project #4 Assigned: Week of March 14, 2016 Due Date: Laboratory Section, Week of April 4, 2016 Report Format and Guidelines for Laboratory
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science Circuits & Electronics Spring 2005
Massachusetts Institute of Technology Department of Electrical Engineering and Computer Science 6.002 Circuits & Electronics Spring 2005 Lab #2: MOSFET Inverting Amplifiers & FirstOrder Circuits Introduction
More information######################################################################
Write a MATLAB program which asks the user to enter three numbers. - The program should figure out the median value and the average value and print these out. Do not use the predefined MATLAB functions
More information4 Experiment 4: DC Motor Voltage to Speed Transfer Function Estimation by Step Response and Frequency Response (Part 2)
4 Experiment 4: DC Motor Voltage to Speed Transfer Function Estimation by Step Response and Frequency Response (Part 2) 4.1 Introduction This lab introduces new methods for estimating the transfer function
More informationContents. An introduction to MATLAB for new and advanced users
An introduction to MATLAB for new and advanced users (Using Two-Dimensional Plots) Contents Getting Started Creating Arrays Mathematical Operations with Arrays Using Script Files and Managing Data Two-Dimensional
More informationEE 435/535: Error Correcting Codes Project 1, Fall 2009: Extended Hamming Code. 1 Introduction. 2 Extended Hamming Code: Encoding. 1.
EE 435/535: Error Correcting Codes Project 1, Fall 2009: Extended Hamming Code Project #1 is due on Tuesday, October 6, 2009, in class. You may turn the project report in early. Late projects are accepted
More informationTime-Frequency analysis of biophysical time series
Time-Frequency analysis of biophysical time series Sept 9 th 2010, NCTU, Taiwan Arnaud Delorme Frequency analysis synchronicity of cell excitation determines amplitude and rhythm of the EEG signal 30-60
More informationDigital Signal Processing Fourier Analysis of Continuous-Time Signals with the Discrete Fourier Transform
Digital Signal Processing Fourier Analysis of Continuous-Time Signals with the Discrete Fourier Transform D. Richard Brown III D. Richard Brown III 1 / 11 Fourier Analysis of CT Signals with the DFT Scenario:
More informationLaboratory Experiment #1 Introduction to Spectral Analysis
J.B.Francis College of Engineering Mechanical Engineering Department 22-403 Laboratory Experiment #1 Introduction to Spectral Analysis Introduction The quantification of electrical energy can be accomplished
More informationUniversiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 1: DISCRETE TIME SIGNALS IN THE TIME DOMAIN
Universiti Malaysia Perlis EKT430: DIGITAL SIGNAL PROCESSING LAB ASSIGNMENT 1: DISCRETE TIME SIGNALS IN THE TIME DOMAIN Pusat Pengajian Kejuruteraan Komputer Dan Perhubungan Universiti Malaysia Perlis
More informationFFT Analyzer. Gianfranco Miele, Ph.D
FFT Analyzer Gianfranco Miele, Ph.D www.eng.docente.unicas.it/gianfranco_miele g.miele@unicas.it Introduction It is a measurement instrument that evaluates the spectrum of a time domain signal applying
More informationENGR 210 Lab 12: Sampling and Aliasing
ENGR 21 Lab 12: Sampling and Aliasing In the previous lab you examined how A/D converters actually work. In this lab we will consider some of the consequences of how fast you sample and of the signal processing
More informationDiscrete Fourier Transform
Discrete Fourier Transform The DFT of a block of N time samples {a n } = {a,a,a 2,,a N- } is a set of N frequency bins {A m } = {A,A,A 2,,A N- } where: N- mn A m = S a n W N n= W N e j2p/n m =,,2,,N- EECS
More informationLab P-4: AM and FM Sinusoidal Signals. We have spent a lot of time learning about the properties of sinusoidal waveforms of the form: ) X
DSP First, 2e Signal Processing First Lab P-4: AM and FM Sinusoidal Signals Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises
More information6 Sampling. Sampling. The principles of sampling, especially the benefits of coherent sampling
Note: Printed Manuals 6 are not in Color Objectives This chapter explains the following: The principles of sampling, especially the benefits of coherent sampling How to apply sampling principles in a test
More informationREAL-TIME PROCESSING ALGORITHMS
CHAPTER 8 REAL-TIME PROCESSING ALGORITHMS In many applications including digital communications, spectral analysis, audio processing, and radar processing, data is received and must be processed in real-time.
More informationSignal Processing First Lab 20: Extracting Frequencies of Musical Tones
Signal Processing First Lab 20: Extracting Frequencies of Musical Tones Pre-Lab and Warm-Up: You should read at least the Pre-Lab and Warm-up sections of this lab assignment and go over all exercises in
More informationHCS / ACN 6389 Speech Perception Lab
HCS / ACN 6389 Speech Perception Lab Course Requirements Matlab problems & lab assignments (40%) Oral presentations (10%) Term project paper (50%) Dr. Peter Assmann Fall 2017 2 Term project: important
More information1.5 The voltage V is given as V=RI, where R and I are resistance matrix and I current vector. Evaluate V given that
Sheet (1) 1.1 The voltage across a discharging capacitor is v(t)=10(1 e 0.2t ) Generate a table of voltage, v(t), versus time, t, for t = 0 to 50 seconds with increment of 5 s. 1.2 Use MATLAB to evaluate
More informationIntroduction to Simulink
EE 460 Introduction to Communication Systems MATLAB Tutorial #3 Introduction to Simulink This tutorial provides an overview of Simulink. It also describes the use of the FFT Scope and the filter design
More informationReference Sources. Prelab. Proakis chapter 7.4.1, equations to as attached
Purpose The purpose of the lab is to demonstrate the signal analysis capabilities of Matlab. The oscilloscope will be used as an A/D converter to capture several signals we have examined in previous labs.
More informationWireless Communication Systems Laboratory #2. Understanding test equipments. The students will be familiar with the following items:
Wireless Communication Systems Laboratory #2 Understanding test equipments Objective The students will be familiar with the following items: Signal generation and analysis tools Description of the laboratory
More informationUnderstanding Digital Signal Processing
Understanding Digital Signal Processing Richard G. Lyons PRENTICE HALL PTR PRENTICE HALL Professional Technical Reference Upper Saddle River, New Jersey 07458 www.photr,com Contents Preface xi 1 DISCRETE
More informationUNIVERSITY OF UTAH ELECTRICAL AND COMPUTER ENGINEERING DEPARTMENT
UNIVERSITY OF UTAH ELECTRICAL AND COMPUTER ENGINEERING DEPARTMENT ECE1020 COMPUTING ASSIGNMENT 3 N. E. COTTER MATLAB ARRAYS: RECEIVED SIGNALS PLUS NOISE READING Matlab Student Version: learning Matlab
More informationPlotting in MATLAB. Trevor Spiteri
Functions and Special trevor.spiteri@um.edu.mt http://staff.um.edu.mt/trevor.spiteri Department of Communications and Computer Engineering Faculty of Information and Communication Technology University
More information6.S02 MRI Lab Acquire MR signals. 2.1 Free Induction decay (FID)
6.S02 MRI Lab 1 2. Acquire MR signals Connecting to the scanner Connect to VMware on the Lab Macs. Download and extract the following zip file in the MRI Lab dropbox folder: https://www.dropbox.com/s/ga8ga4a0sxwe62e/mit_download.zip
More informationTime Series/Data Processing and Analysis (MATH 587/GEOP 505)
Time Series/Data Processing and Analysis (MATH 587/GEOP 55) Rick Aster and Brian Borchers October 7, 28 Plotting Spectra Using the FFT Plotting the spectrum of a signal from its FFT is a very common activity.
More informationLAB 2 SPECTRUM ANALYSIS OF PERIODIC SIGNALS
Eastern Mediterranean University Faculty of Engineering Department of Electrical and Electronic Engineering EENG 360 Communication System I Laboratory LAB 2 SPECTRUM ANALYSIS OF PERIODIC SIGNALS General
More informationC.8 Comb filters 462 APPENDIX C. LABORATORY EXERCISES
462 APPENDIX C. LABORATORY EXERCISES C.8 Comb filters The purpose of this lab is to use a kind of filter called a comb filter to deeply explore concepts of impulse response and frequency response. The
More informationSwedish College of Engineering and Technology Rahim Yar Khan
PRACTICAL WORK BOOK Telecommunication Systems and Applications (TL-424) Name: Roll No.: Batch: Semester: Department: Swedish College of Engineering and Technology Rahim Yar Khan Introduction Telecommunication
More informationPART I: The questions in Part I refer to the aliasing portion of the procedure as outlined in the lab manual.
Lab. #1 Signal Processing & Spectral Analysis Name: Date: Section / Group: NOTE: To help you correctly answer many of the following questions, it may be useful to actually run the cases outlined in the
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 informationECEn 487 Digital Signal Processing Laboratory. Lab 3 FFT-based Spectrum Analyzer
ECEn 487 Digital Signal Processing Laboratory Lab 3 FFT-based Spectrum Analyzer Due Dates This is a three week lab. All TA check off must be completed by Friday, March 14, at 3 PM or the lab will be marked
More informationIntroduction. A Simple Example. 3. fo = 4; %frequency of the sine wave. 4. Fs = 100; %sampling rate. 5. Ts = 1/Fs; %sampling time interval
Introduction In this tutorial, we will discuss how to use the fft (Fast Fourier Transform) command within MATLAB. The fft command is in itself pretty simple, but takes a little bit of getting used to in
More informationME scope Application Note 02 Waveform Integration & Differentiation
ME scope Application Note 02 Waveform Integration & Differentiation The steps in this Application Note can be duplicated using any ME scope Package that includes the VES-3600 Advanced Signal Processing
More informationSignal Processing Toolbox
Signal Processing Toolbox Perform signal processing, analysis, and algorithm development Signal Processing Toolbox provides industry-standard algorithms for analog and digital signal processing (DSP).
More informationAPPENDIX A-5 THE CORRESPONDING FREQUENCY DOMAIN VALIDATION. APPENDIX A-5 shows plots for the corresponding time domain validation response
APPENDIX A-5 THE CORRESPONDING FREQUENCY DOMAIN VALIDATION APPENDIX A-5 shows plots for the corresponding time domain validation response records illustrated in chapter 7 for further validation. The following
More informationAdvanced Lab LAB 6: Signal Acquisition & Spectrum Analysis Using VirtualBench DSA Equipment: Objectives:
Advanced Lab LAB 6: Signal Acquisition & Spectrum Analysis Using VirtualBench DSA Equipment: Pentium PC with National Instruments PCI-MIO-16E-4 data-acquisition board (12-bit resolution; software-controlled
More informationTimbral Distortion in Inverse FFT Synthesis
Timbral Distortion in Inverse FFT Synthesis Mark Zadel Introduction Inverse FFT synthesis (FFT ) is a computationally efficient technique for performing additive synthesis []. Instead of summing partials
More informationThe Fundamentals of Mixed Signal Testing
The Fundamentals of Mixed Signal Testing Course Information The Fundamentals of Mixed Signal Testing course is designed to provide the foundation of knowledge that is required for testing modern mixed
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 informationArmstrong Atlantic State University Engineering Studies MATLAB Marina Sound Processing Primer
Armstrong Atlantic State University Engineering Studies MATLAB Marina Sound Processing Primer Prerequisites The Sound Processing Primer assumes knowledge of the MATLAB IDE, MATLAB help, arithmetic operations,
More information