Biophysical Techniques (BPHS 4090/PHYS 5800)
|
|
- Roxanne Gibson
- 5 years ago
- Views:
Transcription
1 Biophysical Techniques (BPHS 49/PHYS 58) Instructors: Prof. Christopher Bergevin Schedule: MWF :3-2:3 (CB 22) Website: York University Winter 27 Lec.7
2 Goal Develop knowledge and intuiaon dealing w/ 2-D Fourier transforms 2-D case (e.g., image) Why Fourier (i.e., spectral) analysis? And why 2-D?
3 Reminder Franklin & Gosling (953) Sodium thymonucleate fibres give two disanct types of X-ray diagram. The first corresponds to a crystalline form... Watson & Crick (953)
4 Crystal What is a crystal? Something w/ a periodic structure (we ll return to this later in the semester) Halliday & Resnick
5 Crystals made up of complex materials In a nutshell: Thus, the diffracaon payern of a protein crystal is the Fourier transform of the unit cell Ames the Fourier transform of the crystal lazce. The layer is called reciprocal la*ce à Should make sense that we might want to use a theoreacal foundaaon with a set of periodic basis funcaons... NolAng
6 Goal Develop knowledge and intuiaon dealing w/ 2D Fourier transforms 2-D case (e.g., image).8 -D case time waveform Ame waveform recorded from mic Pressure [arb] Time [s] à Will first explore (from a computaaonal viewpoint) -D Fourier transforms... Note: Following background notes taken from PHYS 23 W6
7 % ### EXbuildImpulse.m ###.3.4! % Code to visually build up a signal by successively adding higher and! % higher frequency terms from corresponding FFT! clear; clf;! % ! SR= 44; Npoints= 892; % sample rate [Hz]! % length of fft window (# of points) [should ideally be 2^N]! % [time window will be the same length]! INDXon= ; % index at which click turns 'on' (i.e., go from to )! INDXoff= ; % index at which click turns 'off' (i.e., go from to )! % !! dt= /SR; % spacing of time steps! freq= [:Npoints/2]; % create a freq. array (for FFT bin labeling)! freq= SR*freq./Npoints;! t=[:/sr:(npoints-)/sr]; % create an appropriate array of time points! % build signal! clktemp= zeros(,npoints); clktemp2= ones(,indxoff-indxon);! signal= [clktemp(:indxon-) clktemp2 clktemp(indxoff:end)];! % ! % *******! % plot time waveform of signal! if ==! figure(); clf; plot(t*,signal,'ko-','markersize',5)! grid on; hold on; xlabel('time [ms]'); ylabel('signal'); title('time Waveform')! end! % *******! % now compute/plot FFT of the signal! sigspec= rfft(signal);! % MAGNITUDE! figure(2); clf;! subplot(2); plot(freq/,db(sigspec),'ko-','markersize',3)! hold on; grid on; ylabel('magnitude [db]'); title('spectrum')! % PHASE! subplot(22); plot(freq/,cycs(sigspec),'ko-','markersize',3)! xlabel('frequency [khz]'); ylabel('phase [cycles]'); grid on;! % *******! % now make animation of click getting built up, using the info from the FFT! sum= zeros(,numel(t)); % (initial) array for reconstructed waveform! figure(3); clf;! for nn=:numel(freq)! sum= sum+ abs(sigspec(nn))*cos(2*pi*freq(nn)*t + angle(sigspec(nn)));! plot(t,sum); xlabel('time [s]');! legend(['highest freq= ',num2str(freq(nn)/),' khz'])! pause(2/(nn))! end! EXbuildImpulse.m
8 6 x 4 Highest freq= khz 4 x 3 Highest freq= khz Time Waveform Temporal 7 Spectrum EXbuildImpulse.m Spectral Signal is an impulse (i.e., a delta function) Magnitude [db] Signal Time [ms] Phase [cycles] Spectral representation has flat amplitude and a group delay Frequency [khz] 5 4 Reconstruct waveform by adding sinusoids (only lowest frequency here) Now the first 5 terms are included Time [s] Time [s] à Eventually all the sinusoids add up such that things cancel out everywhere except at the point of the impulse!
9 % ### EXspecREP3.m ###.29.4! % Example code to just fiddle with basics of discrete FFTs and connections! % back to common real-valued time waveforms! % --> Demonstrates several useful concepts such as 'quantizing' the frequency! % Requires: rfft.m, irfft.m, cycs.m, db.m, cyc.m! % ! % Stimulus Type Legend! % stimt= - non-quantized sinusoid! % stimt= - quantized sinusoid! % stimt= 2 - one quantized sinusoid, one un-quantized sinusoid! % stimt= 3 - two quantized sinusoids! % stimt= 4 - click I.e., an impulse)! % stimt= 5 - noise (uniform in time)! % stimt= 6 - chirp (flat mag.)! % stimt= 7 - noise (Gaussian; flat spectrum, random phase)! % stimt= 8 - exponentially decaying sinusoid (i.e., HO impulse response)!! clear; clf;! % ! SR= 44; % sample rate [Hz]! Npoints= 892; % length of fft window (# of points) [should ideally be 2^N]! % [time window will be the same length]! stimt= 8; % Stimulus Type (see legend above)! f= 258.; % Frequency (for waveforms w/ tones) [Hz]! ratio=.22; % specify f2/f2 ratio (for waveforms w/ two tones)! % Note: Other stimulus parameters can be changed below! % ! dt= /SR; % spacing of time steps! freq= [:Npoints/2]; % create a freq. array (for FFT bin labeling)! freq= SR*freq./Npoints;! % quantize the freq. (so to have an integral # of cycles in time window)! df = SR/Npoints;! fq= ceil(f/df)*df; % quantized natural freq.! t=[:/sr:(npoints-)/sr]; % create an array of time points, Npoints long! % ----! % compute stimulus! if stimt== % non-quantized sinusoid! signal= cos(2*pi*f*t);! disp(sprintf(' \n *Stimulus* - (non-quantized) sinusoid, f = %g Hz \n', f));! disp(sprintf('specified freq. = %g Hz', f));! elseif stimt== % quantized sinusoid! signal= cos(2*pi*fq*t);! disp(sprintf(' \n *Stimulus* - quantized sinusoid, f = %g Hz \n', fq));! disp(sprintf('specified freq. = %g Hz', f));! disp(sprintf('quantized freq. = %g Hz', fq));! elseif stimt==2 % one quantized sinusoid, one un-quantized sinusoid! signal= cos(2*pi*fq*t) + cos(2*pi*ratio*fq*t);! disp(sprintf(' \n *Stimulus* - two sinusoids (one quantized, one not) \n'));! elseif stimt==3 % two quantized sinusoids! fq2= ceil(ratio*f/df)*df;! signal= cos(2*pi*fq*t) + cos(2*pi*fq2*t);! disp(sprintf(' \n *Stimulus* - two sinusoids (both quantized) \n'));! elseif stimt==4 % click! CLKon= ; % index at which click turns 'on' (starts at )! CLKoff= ; % index at which click turns 'off'! clktemp= zeros(,npoints);! clktemp2= ones(,clkoff-clkon);! signal= [clktemp(:clkon-) clktemp2 clktemp(clkoff:end)];! disp(sprintf(' \n *Stimulus* - Click \n'));! elseif stimt==5 % noise (flat)! signal= rand(,npoints);! disp(sprintf(' \n *Stimulus* - Noise \n'));! elseif stimt==6 % chirp (flat)! fs= 2.; % if a chirp (stimt=2) starting freq. [Hz] [freq. swept linearly w/ time]! fe= 4.; % ending freq. (energy usually extends twice this far out)! fsq= ceil(fs/df)*df; %quantize the start/end freqs. (necessary?)! feq= ceil(fe/df)*df;! % LINEAR sweep rate! fswp= fsq + (feq-fsq)*(sr/npoints)*t;! signal = sin(2*pi*fswp.*t)';! disp(sprintf(' \n *Stimulus* - Chirp \n'));! EXspecREP3.m elseif stimt==7 % noise (Gaussian)! Asize=Npoints/2 +;! % create array of complex numbers w/ random phase and unit magnitude! for n=:asize! theta= rand*2*pi;! N2(n)= exp(i*theta);! end! N2=N2';! % now take the inverse FFT of that using Chris' irfft.m code! tnoise=irfft(n2);! % scale it down so #s are between - and (i.e. normalize)! if (abs(min(tnoise)) > max(tnoise))! tnoise= tnoise/abs(min(tnoise));! else! tnoise= tnoise/max(tnoise);! end! signal= tnoise;! disp(sprintf(' \n *Noise* - Gaussian, flat-spectrum \n'));! elseif stimt==8 % exponentially decaying cos! alpha= 5;! signal= exp(-alpha*t).*sin(2*pi*fq*t);! disp(sprintf(' \n *Exponentially decaying (quantized) sinusoid* \n'));! end!! % ! % *******! figure(); clf % plot time waveform of signal! plot(t*,signal,'k.-','markersize',5); grid on; hold on;! xlabel('time [ms]'); ylabel('signal'); title('time Waveform')! % *******! % now plot rfft of the signal! % NOTE: rfft just takes /2 of fft.m output and nomalizes! sigspec= rfft(signal);! figure(2); clf; % MAGNITUDE! subplot(2)! plot(freq/,db(sigspec),'ko-','markersize',3)! hold on; grid on;! ylabel('magnitude [db]')! title('spectrum')! subplot(22) % PHASE! plot(freq/,cycs(sigspec),'ko-','markersize',3)! xlabel('frequency [khz]'); ylabel('phase [cycles]'); grid on;! % ! % play the stimuli as an output sound?! if (==), sound(signal,sr); end! % ! % compute inverse Fourier transform and plot?! if ==! figure();! signalinv= irfft(sigspec);! plot(t*,signalinv,'rx','markersize',4)! legend('original waveform','inverse transformed')! end
10 Fourier transforms of basic (-D) waveforms stimt= - non-quantized sinusoid! EXspecREP3.m SR= 44; % sample rate [Hz]! Npoints= 892; % length of fft window! Time domain Spectral domain Time Waveform Spectrum Magnitude [db] Signal Time [ms] Phase [cycles] Frequency [khz] Ø Magnitude shows a peak at the sinusoid s frequency Note: The phase is unwrapped in all the spectral plots
11 Fourier transforms of basic (-D) waveforms EXspecREP3.m stimt= 3 - two quantized sinusoids! Time domain Spectral domain 2 Time Waveform Spectrum.5.5 Magnitude [db] 2 3 Signal Phase [cycles] Time [ms] Frequency [khz] Ø Magnitude shows two peaks (note the beating in the time domain)
12 Fourier transforms of basic (-D) waveforms EXspecREP3.m stimt= 4 - click (i.e., an impulse)! Time domain Spectral domain Time Waveform 7 Spectrum Magnitude [db] Signal Phase [cycles] Time [ms] Frequency [khz] Ø Click has a flat magnitude (This is also a good place to mention the concept of a group delay )
13 Fourier transforms of basic (-D) waveforms EXspecREP3.m stimt= 5 - noise (uniform distribution)! Time domain Spectral domain Time Waveform 2 Spectrum.9 Signal Magnitude [db] Time [ms] Phase [cycles] Frequency [khz] Ø Magnitude is flat-ish (on log scale), but actually noisy. Phase is noisy too.
14 Fourier transforms of basic (-D) waveforms EXspecREP3.m stimt= 7 - noise (Gaussian distribution)! Time domain Spectral domain Time Waveform 44 Spectrum Magnitude [db] Signal Phase [cycles] Time [ms] Frequency [khz] Ø Magnitude is flat just like an impulse (i.e., flat), but the phase is random
15 Fourier transforms of basic (-D) waveforms Impulse Noise Time Waveform Time Waveform Time domain Signal.5.4 Signal Time [ms] Time [ms] 7 Spectrum 44 Spectrum Magnitude [db] Magnitude [db] Spectral domain Phase [cycles] Frequency [khz] Phase [cycles] Frequency [khz] à Remarkable that the magnitudes are idenacal (more or less) between two signals with such different properaes. The key difference here is the phase: Timing is a cri/cal piece of the puzzle!
16 Fourier transforms of basic (-D) waveforms EXspecREP3.m stimt= 6 - chirp (flat mag.)! Time domain Spectral domain Time Waveform 2 Spectrum Magnitude [db] Signal Time [ms] Hard to see on this timescale, but frequency is changing (increasing) with time Phase [cycles] Frequency [khz]
17 Fourier transforms of basic (-D) waveforms EXspecREP3.m stimt= 8 - exponentially decaying sinusoid! Time domain Spectral domain Time Waveform 2 Spectrum Magnitude [db] Signal Phase [cycles] Time [ms] Frequency [khz] Ø This seems to look familiar...
18 ConnecAon back to the harmonic oscillator! Time Waveform Signal Time [ms] 2 Spectrum Magnitude [db] Phase [cycles] Frequency [khz] Ø The steady-state response of the sinusoidally-driven harmonic harmonic oscillator acts like a band-pass filter Ø DisAncAon between steady-state response & impulse response [we ll come back to this]
19 2-D Fourier analysis.8 -D case time waveform Ame waveform recorded from mic.6.4 Pressure [arb] D case (e.g., image) Time [s] Ø Ø Same basic idea between -D and 2-D (though math becomes a bit more complicated) Just as a waveform can have temporal frequencies, an image can have spa/al frequencies
20 2-D Fourier analysis Note: Independent variables (x and y here) can represent any physical quanaty, but for 2-D we typically use posiaon in Cartesian plane (rather than Ame)
21 NotaAon re Fourier analysis -D 2-D
22 2-D Fourier analysis A. Zisserman (Oxford)
23 2-D Fourier analysis Note: Only ½ of the informaaon is being shown for the spectral domain (i.e., just the magnitude) A. Zisserman (Oxford)
24 2-D Fourier analysis Spectral domain SpaAal domain Note: Only ½ of the informaaon is being shown for the spectral domain (i.e., just the magnitude) Hobbie & Roth
25 2-D Fourier analysis A. Zisserman (Oxford)
26 2-D Fourier analysis Swap phases, then inverse transform to get the image back A. Zisserman (Oxford)
Complex Sounds. Reading: Yost Ch. 4
Complex Sounds Reading: Yost Ch. 4 Natural Sounds Most sounds in our everyday lives are not simple sinusoidal sounds, but are complex sounds, consisting of a sum of many sinusoids. The amplitude and frequency
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 informationSAMPLING THEORY. Representing continuous signals with discrete numbers
SAMPLING THEORY Representing continuous signals with discrete numbers Roger B. Dannenberg Professor of Computer Science, Art, and Music Carnegie Mellon University ICM Week 3 Copyright 2002-2013 by Roger
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 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 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 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 informationPHYC 500: Introduction to LabView. Exercise 9 (v 1.1) Spectral content of waveforms. M.P. Hasselbeck, University of New Mexico
PHYC 500: Introduction to LabView M.P. Hasselbeck, University of New Mexico Exercise 9 (v 1.1) Spectral content of waveforms This exercise provides additional experience with the Waveform palette, along
More informationDiscrete 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 informationStructure of Speech. Physical acoustics Time-domain representation Frequency domain representation Sound shaping
Structure of Speech Physical acoustics Time-domain representation Frequency domain representation Sound shaping Speech acoustics Source-Filter Theory Speech Source characteristics Speech Filter characteristics
More informationThe Fundamentals of FFT-Based Signal Analysis and Measurement Michael Cerna and Audrey F. Harvey
Application ote 041 The Fundamentals of FFT-Based Signal Analysis and Measurement Michael Cerna and Audrey F. Harvey Introduction The Fast Fourier Transform (FFT) and the power spectrum are powerful tools
More informationLecture 3 Complex Exponential Signals
Lecture 3 Complex Exponential Signals Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/3/1 1 Review of Complex Numbers Using Euler s famous formula for the complex exponential The
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 informationSignal Processing. Naureen Ghani. December 9, 2017
Signal Processing Naureen Ghani December 9, 27 Introduction Signal processing is used to enhance signal components in noisy measurements. It is especially important in analyzing time-series data in neuroscience.
More informationELT Receiver Architectures and Signal Processing Fall Mandatory homework exercises
ELT-44006 Receiver Architectures and Signal Processing Fall 2014 1 Mandatory homework exercises - Individual solutions to be returned to Markku Renfors by email or in paper format. - Solutions are expected
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 informationHideo Okawara s Mixed Signal Lecture Series. DSP-Based Testing Fundamentals 6 Spectrum Analysis -- FFT
Hideo Okawara s Mixed Signal Lecture Series DSP-Based Testing Fundamentals 6 Spectrum Analysis -- FFT Verigy Japan October 008 Preface to the Series ADC and DAC are the most typical mixed signal devices.
More informationFigure 1: Block diagram of Digital signal processing
Experiment 3. Digital Process of Continuous Time Signal. Introduction Discrete time signal processing algorithms are being used to process naturally occurring analog signals (like speech, music and images).
More informationFourier Series and Gibbs Phenomenon
Fourier Series and Gibbs Phenomenon University Of Washington, Department of Electrical Engineering This work is produced by The Connexions Project and licensed under the Creative Commons Attribution License
More informationSignal Processing for Digitizers
Signal Processing for Digitizers Modular digitizers allow accurate, high resolution data acquisition that can be quickly transferred to a host computer. Signal processing functions, applied in the digitizer
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 informationFourier Theory & Practice, Part I: Theory (HP Product Note )
Fourier Theory & Practice, Part I: Theory (HP Product Note 54600-4) By: Robert Witte Hewlett-Packard Co. Introduction: This product note provides a brief review of Fourier theory, especially the unique
More informationLab S-8: Spectrograms: Harmonic Lines & Chirp Aliasing
DSP First, 2e Signal Processing First Lab S-8: Spectrograms: Harmonic Lines & Chirp Aliasing Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification:
More informationFourier Analysis. Chapter Introduction Distortion Harmonic Distortion
Chapter 5 Fourier Analysis 5.1 Introduction The theory, practice, and application of Fourier analysis are presented in the three major sections of this chapter. The theory includes a discussion of Fourier
More informationWindow Method. designates the window function. Commonly used window functions in FIR filters. are: 1. Rectangular Window:
Window Method We have seen that in the design of FIR filters, Gibbs oscillations are produced in the passband and stopband, which are not desirable features of the FIR filter. To solve this problem, window
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 informationSignals & Systems for Speech & Hearing. Week 6. Practical spectral analysis. Bandpass filters & filterbanks. Try this out on an old friend
Signals & Systems for Speech & Hearing Week 6 Bandpass filters & filterbanks Practical spectral analysis Most analogue signals of interest are not easily mathematically specified so applying a Fourier
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 information14 fasttest. Multitone Audio Analyzer. Multitone and Synchronous FFT Concepts
Multitone Audio Analyzer The Multitone Audio Analyzer (FASTTEST.AZ2) is an FFT-based analysis program furnished with System Two for use with both analog and digital audio signals. Multitone and Synchronous
More informationSpectrum Analysis: The FFT Display
Spectrum Analysis: The FFT Display Equipment: Capstone, voltage sensor 1 Introduction It is often useful to represent a function by a series expansion, such as a Taylor series. There are other series representations
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 informationFilter Banks I. Prof. Dr. Gerald Schuller. Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany. Fraunhofer IDMT
Filter Banks I Prof. Dr. Gerald Schuller Fraunhofer IDMT & Ilmenau University of Technology Ilmenau, Germany 1 Structure of perceptual Audio Coders Encoder Decoder 2 Filter Banks essential element of most
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 informationThe Fast Fourier Transform
The Fast Fourier Transform Basic FFT Stuff That s s Good to Know Dave Typinski, Radio Jove Meeting, July 2, 2014, NRAO Green Bank Ever wonder how an SDR-14 or Dongle produces the spectra that it does?
More informationTransfer Function (TRF)
(TRF) Module of the KLIPPEL R&D SYSTEM S7 FEATURES Combines linear and nonlinear measurements Provides impulse response and energy-time curve (ETC) Measures linear transfer function and harmonic distortions
More informationNotes on Fourier transforms
Fourier Transforms 1 Notes on Fourier transforms The Fourier transform is something we all toss around like we understand it, but it is often discussed in an offhand way that leads to confusion for those
More informationExperiment One: Generating Frequency Modulation (FM) Using Voltage Controlled Oscillator (VCO)
Experiment One: Generating Frequency Modulation (FM) Using Voltage Controlled Oscillator (VCO) Modified from original TIMS Manual experiment by Mr. Faisel Tubbal. Objectives 1) Learn about VCO and how
More informationPractical Applications of the Wavelet Analysis
Practical Applications of the Wavelet Analysis M. Bigi, M. Jacchia, D. Ponteggia ALMA International Europe (6- - Frankfurt) Summary Impulse and Frequency Response Classical Time and Frequency Analysis
More informationDefinition of Sound. Sound. Vibration. Period - Frequency. Waveform. Parameters. SPA Lundeen
Definition of Sound Sound Psychologist's = that which is heard Physicist's = a propagated disturbance in the density of an elastic medium Vibrator serves as the sound source Medium = air 2 Vibration Periodic
More informationAdvanced Audiovisual Processing Expected Background
Advanced Audiovisual Processing Expected Background As an advanced module, we will not cover introductory topics in lecture. You are expected to already be proficient with all of the following topics,
More informationSECTION 7: FREQUENCY DOMAIN ANALYSIS. MAE 3401 Modeling and Simulation
SECTION 7: FREQUENCY DOMAIN ANALYSIS MAE 3401 Modeling and Simulation 2 Response to Sinusoidal Inputs Frequency Domain Analysis Introduction 3 We ve looked at system impulse and step responses Also interested
More informationDepartment of Electronic Engineering NED University of Engineering & Technology. LABORATORY WORKBOOK For the Course SIGNALS & SYSTEMS (TC-202)
Department of Electronic Engineering NED University of Engineering & Technology LABORATORY WORKBOOK For the Course SIGNALS & SYSTEMS (TC-202) Instructor Name: Student Name: Roll Number: Semester: Batch:
More informationDigital Signal Processing Lecture 1 - Introduction
Digital Signal Processing - Electrical Engineering and Computer Science University of Tennessee, Knoxville August 20, 2015 Overview 1 2 3 4 Basic building blocks in DSP Frequency analysis Sampling Filtering
More informationDCSP-10: DFT and PSD. Jianfeng Feng. Department of Computer Science Warwick Univ., UK
DCSP-10: DFT and PSD Jianfeng Feng Department of Computer Science Warwick Univ., UK Jianfeng.feng@warwick.ac.uk http://www.dcs.warwick.ac.uk/~feng/dcsp.html DFT Definition: The discrete Fourier transform
More informationThe quality of the transmission signal The characteristics of the transmission medium. Some type of transmission medium is required for transmission:
Data Transmission The successful transmission of data depends upon two factors: The quality of the transmission signal The characteristics of the transmission medium Some type of transmission medium is
More informationSet-up. Equipment required: Your issued Laptop MATLAB ( if you don t already have it on your laptop)
All signals found in nature are analog they re smooth and continuously varying, from the sound of an orchestra to the acceleration of your car to the clouds moving through the sky. An excerpt from http://www.netguru.net/ntc/ntcc5.htm
More informationAcoustics, signals & systems for audiology. Week 4. Signals through Systems
Acoustics, signals & systems for audiology Week 4 Signals through Systems Crucial ideas Any signal can be constructed as a sum of sine waves In a linear time-invariant (LTI) system, the response to a sinusoid
More informationEE 3302 LAB 1 EQIUPMENT ORIENTATION
EE 3302 LAB 1 EQIUPMENT ORIENTATION Pre Lab: Calculate the theoretical gain of the 4 th order Butterworth filter (using the formula provided. Record your answers in Table 1 before you come to class. Introduction:
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 informationDIGITAL SIGNAL PROCESSING CCC-INAOE AUTUMN 2015
DIGITAL SIGNAL PROCESSING CCC-INAOE AUTUMN 2015 Fourier Transform Properties Claudia Feregrino-Uribe & Alicia Morales Reyes Original material: Rene Cumplido "The Scientist and Engineer's Guide to Digital
More informationAppendix. Harmonic Balance Simulator. Page 1
Appendix Harmonic Balance Simulator Page 1 Harmonic Balance for Large Signal AC and S-parameter Simulation Harmonic Balance is a frequency domain analysis technique for simulating distortion in nonlinear
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 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 informationME scope Application Note 01 The FFT, Leakage, and Windowing
INTRODUCTION ME scope Application Note 01 The FFT, Leakage, and Windowing NOTE: The steps in this Application Note can be duplicated using any Package that includes the VES-3600 Advanced Signal Processing
More informationLLS - Introduction to Equipment
Published on Advanced Lab (http://experimentationlab.berkeley.edu) Home > LLS - Introduction to Equipment LLS - Introduction to Equipment All pages in this lab 1. Low Light Signal Measurements [1] 2. Introduction
More informationLab S-7: Spectrograms of AM and FM Signals. 2. Study the frequency resolution of the spectrogram for two closely spaced sinusoids.
DSP First, 2e Signal Processing First Lab S-7: Spectrograms of AM and FM Signals Pre-Lab: Read the Pre-Lab and do all the exercises in the Pre-Lab section prior to attending lab. Verification: The Exercise
More informationWavelets and wavelet convolution and brain music. Dr. Frederike Petzschner Translational Neuromodeling Unit
Wavelets and wavelet convolution and brain music Dr. Frederike Petzschner Translational Neuromodeling Unit 06.03.2015 Recap Why are we doing this? We know that EEG data contain oscillations. Or goal is
More informationMulti-Path Fading Channel
Instructor: Prof. Dr. Noor M. Khan Department of Electronic Engineering, Muhammad Ali Jinnah University, Islamabad Campus, Islamabad, PAKISTAN Ph: +9 (51) 111-878787, Ext. 19 (Office), 186 (Lab) Fax: +9
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 informationCross-Frequency Coupling. Naureen Ghani. April 28, 2018
Cross-Frequency Coupling Naureen Ghani April 28, 2018 Introduction The interplay of excitation and inhibition is a fundamental feature of cortical information processing. The opposing actions of pyramidal
More informationADC Clock Jitter Model, Part 2 Random Jitter
db ADC Clock Jitter Model, Part 2 Random Jitter In Part 1, I presented a Matlab function to model an ADC with jitter on the sample clock, and applied it to examples with deterministic jitter. Now we ll
More informationSynthesis: From Frequency to Time-Domain
Synthesis: From Frequency to Time-Domain I Synthesis is a straightforward process; it is a lot like following a recipe. I Ingredients are given by the spectrum X (f )={(X 0, 0), (X 1, f 1 ), (X 1, f 1),...,
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 informationECE 440L. Experiment 1: Signals and Noise (1 week)
ECE 440L Experiment 1: Signals and Noise (1 week) I. OBJECTIVES Upon completion of this experiment, you should be able to: 1. Use the signal generators and filters in the lab to generate and filter noise
More informationDigital Signal Processing +
Digital Signal Processing + Nikil Dutt UC Irvine ICS 212 Winter 2005 + Material adapted from Tony Givargis & Rajesh Gupta Templates from Prabhat Mishra ICS212 WQ05 (Dutt) DSP 1 Introduction Any interesting
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 informationFourier transforms, SIM
Fourier transforms, SIM Last class More STED Minflux Fourier transforms This class More FTs 2D FTs SIM 1 Intensity.5 -.5 FT -1.5 1 1.5 2 2.5 3 3.5 4 4.5 5 6 Time (s) IFT 4 2 5 1 15 Frequency (Hz) ff tt
More informationSignal Characteristics
Data Transmission The successful transmission of data depends upon two factors:» The quality of the transmission signal» The characteristics of the transmission medium Some type of transmission medium
More informationPulsed S-Parameter Measurements using the ZVA network Analyzer
Pulsed S-Parameter Measurements using the ZVA network Analyzer 1 Pulse Profile measurements ZVA Advanced Network Analyser 3 Motivation for Pulsed Measurements Typical Applications Avoid destruction of
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 informationSpectrum Analyzer TEN MINUTE TUTORIAL
Spectrum Analyzer TEN MINUTE TUTORIAL November 4, 2011 Summary The Spectrum Analyzer option allows users who are familiar with RF spectrum analyzers to start using the FFT with little or no concern about
More informationVU Signal and Image Processing. Torsten Möller + Hrvoje Bogunović + Raphael Sahann
052600 VU Signal and Image Processing Torsten Möller + Hrvoje Bogunović + Raphael Sahann torsten.moeller@univie.ac.at hrvoje.bogunovic@meduniwien.ac.at raphael.sahann@univie.ac.at vda.cs.univie.ac.at/teaching/sip/17s/
More informationMassachusetts Institute of Technology Dept. of Electrical Engineering and Computer Science Spring Semester, Introduction to EECS 2
Massachusetts Institute of Technology Dept. of Electrical Engineering and Computer Science Spring Semester, 2007 6.082 Introduction to EECS 2 Lab #3: Modulation and Filtering Goal:... 2 Instructions:...
More informationMATLAB Assignment. The Fourier Series
MATLAB Assignment The Fourier Series Read this carefully! Submit paper copy only. This project could be long if you are not very familiar with Matlab! Start as early as possible. This is an individual
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 informationMATRIX TECHNICAL NOTES MTN-109
200 WOOD AVENUE, MIDDLESEX, NJ 08846 PHONE (732) 469-9510 E-mail sales@matrixtest.com MATRIX TECHNICAL NOTES MTN-109 THE RELATIONSHIP OF INTERCEPT POINTS COMPOSITE DISTORTIONS AND NOISE POWER RATIOS Amplifiers,
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 informationData Acquisition Systems. Signal DAQ System The Answer?
Outline Analysis of Waveforms and Transforms How many Samples to Take Aliasing Negative Spectrum Frequency Resolution Synchronizing Sampling Non-repetitive Waveforms Picket Fencing A Sampled Data System
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 informationDigital Signal Processing ETI
2012 Digital Signal Processing ETI265 2012 Introduction In the course we have 2 laboratory works for 2012. Each laboratory work is a 3 hours lesson. We will use MATLAB for illustrate some features in digital
More informationDigital Signal Processing ETI
2011 Digital Signal Processing ETI265 2011 Introduction In the course we have 2 laboratory works for 2011. Each laboratory work is a 3 hours lesson. We will use MATLAB for illustrate some features in digital
More informationLecture 7 Frequency Modulation
Lecture 7 Frequency Modulation Fundamentals of Digital Signal Processing Spring, 2012 Wei-Ta Chu 2012/3/15 1 Time-Frequency Spectrum We have seen that a wide range of interesting waveforms can be synthesized
More informationReading: Johnson Ch , Ch.5.5 (today); Liljencrants & Lindblom; Stevens (Tues) reminder: no class on Thursday.
L105/205 Phonetics Scarborough Handout 7 10/18/05 Reading: Johnson Ch.2.3.3-2.3.6, Ch.5.5 (today); Liljencrants & Lindblom; Stevens (Tues) reminder: no class on Thursday Spectral Analysis 1. There are
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 informationADC Clock Jitter Model, Part 1 Deterministic Jitter
ADC Clock Jitter Model, Part 1 Deterministic Jitter Analog to digital converters (ADC s) have several imperfections that effect communications signals, including thermal noise, differential nonlinearity,
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 informationIMAC 27 - Orlando, FL Shaker Excitation
IMAC 27 - Orlando, FL - 2009 Peter Avitabile UMASS Lowell Marco Peres The Modal Shop 1 Dr. Peter Avitabile Objectives of this lecture: Overview some shaker excitation techniques commonly employed in modal
More informationEE 442 Homework #3 Solutions (Spring 2016 Due February 13, 2017 ) Print out homework and do work on the printed pages.
NAME Solutions EE 44 Homework #3 Solutions (Spring 06 Due February 3, 07 ) Print out homework and do work on the printed pages. Textbook: B. P. Lathi & Zhi Ding, Modern Digital and Analog Communication
More informationSignals, sampling & filtering
Signals, sampling & filtering Scientific Computing Fall, 2018 Paul Gribble 1 Time domain representation of signals 1 2 Frequency domain representation of signals 2 3 Fast Fourier transform (FFT) 2 4 Sampling
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 informationChapter 1: Introduction. EET-223: RF Communication Circuits Walter Lara
Chapter 1: Introduction EET-223: RF Communication Circuits Walter Lara Introduction Electronic communication involves transmission over medium from source to destination Information can contain voice,
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 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 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 informationThe 29 th Annual ARRL and TAPR Digital Communications Conference. DSP Short Course Session 1: DSP Intro and Basics. Rick Muething, KN6KB/AAA9WK
The 29 th Annual ARRL and TAPR Digital Communications Conference DSP Short Course Session 1: DSP Intro and Basics Rick Muething, KN6KB/AAA9WK Session 1 Overview What is DSP? Why is DSP better/different
More informationSpur Detection, Analysis and Removal Stable32 W.J. Riley Hamilton Technical Services
Introduction Spur Detection, Analysis and Removal Stable32 W.J. Riley Hamilton Technical Services Stable32 Version 1.54 and higher has the capability to detect, analyze and remove discrete spectral components
More informationLAB 2 Machine Perception of Music Computer Science 395, Winter Quarter 2005
1.0 Lab overview and objectives This lab will introduce you to displaying and analyzing sounds with spectrograms, with an emphasis on getting a feel for the relationship between harmonicity, pitch, and
More informationDSP First. Laboratory Exercise #7. Everyday Sinusoidal Signals
DSP First Laboratory Exercise #7 Everyday Sinusoidal Signals This lab introduces two practical applications where sinusoidal signals are used to transmit information: a touch-tone dialer and amplitude
More informationCommunication Systems Theory for Undergraduate Students using Matlab
Paper ID #16213 Communication Systems Theory for Undergraduate Students using Matlab Dr. Chandana K.K. Jayasooriya, Engineering Technology Division, University of Pittsburgh at Johnstown Chandana K.K.
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 information