a. Use (at least) window lengths of 256, 1024, and 4096 samples to compute the average spectrum using a window overlap of 0.5.
|
|
- Vincent Wheeler
- 6 years ago
- Views:
Transcription
1 1. Download the file signal.mat from the website. This is continuous 10 second recording of a signal sampled at 1 khz. Assume the noise is ergodic in time and that it is white. I used the MATLAB Signal Processing Toolbox function sptool, a GUI for analyzing signals to create the plots below. This is a very convenient tool for creating spectra and comparing window types, lengths, spectral estimation methods, etc. a. Use (at least) window lengths of 256, 1024, and 4096 samples to compute the average spectrum using a window overlap of 0.5. Figure 1 compares the spectrum of our signal calculated with a rectangular window of 256, 1024, and 4096 points (blue, green, red, lines, respectively). Lengthening the window increases the spectral resolution and increases the energy in correct frequency channel as the resolution increases. Spectral lines in the signal that were not seen with 256-point windows become apparent with 4096 point windows. Figure 1 b. Use Matlab s sptool (Figure 1) to perform spectral analysis using a sub-optimal window based on the performance characteristics discussed and the table and graph handed out in class. Use (at least) window lengths of 256, 1024, and 4096 samples to compute the average spectrum using a window overlap of 0.5. c. Repeat the analysis with a more optimal window with the same 0.5 overlap.
2 Figure 2 demonstrates that, given a 256-point window, the window type affects the spectral resolution and sidelobe level. Shown are rectangle, Hann, Chebyshev(blue, green, red lines, respectively). Figure 2 Figure 3 demonstrates that when long windows are used, side-lobe suppression becomes the main issue point windows (rectangle, Hann, Chebyshev; blue, green, red lines, respectively) are shown.
3 Figure 3 d. Tabulate the spectral content of the signal. The correct answer, i.e. the spectral frequencies that were added to the original sequence were at: f1 = ; % all frequencies in Hz. f2 = ; f3 = ; f4 = ; f5 = ; e. Discuss how your results differ for the window types and durations used. See discussion above. Extra Credit: Design a bandpass filter using the fdatool to pass the non-zero frequencies in the spectrum of the signal. This is worth 10 extra points. Design a bandpass filter with bandwidth Hz to suppress noise from outside this range. An appropriate selection of filter type would favor low ripple in the pass band over narrow transition bands in order to minimize ringing artifact.
4 2. In blood oxygen dependent contrast (BOLD) imaging the flow of blood to the brain increases for voxels that are activated during nervous system stimulation, e.g. visual or motor cortex. Our goal is to isolate the pixels that are responding to the stimulation using spectral analysis techniques. These voxels are said to be activated. Download the series of 80 brain images (acquired at a temporal resolution of 2 sec/image) from a functional MRI data set (titled fmri_images.mat ) available on the website. Assume these are the data acquired using a simple finger tapping paradigm recorded in paradigm.mat. a. The measured signal is the hemodynamic response to the paradigm (Figure 2) plus noise. Note that one feature of BOLD data is that the signal changes due to increased blood flow are relatively small. First characterize the noise in the fmri images by calculating the 2D noise power assuming each image represents an independent realization of the underlying noise process. What type of noise appears to be represented in this problem? The approach I took to evaluate the noise process in this experiment was to assume that 2 consecutive images were independent noise realizations. This allows one to remove the signal from the process to evaluate the noise using the autocovariance noise power in directly in 2-dimensions. >> load fmri_images >> I_sub = I{2}-I{1}; >> I_sub = I_sub./sqrt(2); % Assume no covariance in time thus noise variance is root 2 % larger for subtracted image. >> I_ac = xcorr2(i_sub); Figure 4: Raw noise process appears to be structured. Recall the effects of the magnitude operation on a complex noise process. % Matrix is complex but imaginary terms are zero; need to scale the xcorrelation by the matrix size. >> figure;imagesc(abs(i_ac)./(128.^2));axis('image');colormap('gray');colorbar % >> figure;imagesc(i_sub);axis('image');colormap('gray');colorbar
5 % Calculate the noise power using the DFT of the autocovariance >> i_ac = fftshift(fft2(fftshift(i_ac))); >> figure;imagesc(abs(i_ac)./(128.^2));axis('image');colormap('gray');colorbar >> brighten(0.5) % Display 1D profile through the 2D noise power image >> figure;plot(abs(i{1}(128,:))) Figure 5: The noise process in space appears to be white with a variance of ±2.5 gray levels. (a) The 2D autocovariance is a 2D delta-function with amplitude ~2.5 (b) White noise power b. Now estimate the noise process in time by calculating the noise power along the time dimension. Remember that the MRI data have undergone a magnitude operation. The approach I took here assumed that the time-dependent noise process was well represented by a single sample within tissue that was unlikely to show activation (i.e. fat or skull signal around the brain in the Shepp-Logan phantom). I then subtracted the mean value from this time-dependent signal: % Obtain the time-dependent signal and store in array t1(). >> for i = 1:80, t1(i) = I{i}(8,64); end % Subtract mean from t1 >> tnoise = t1-mean(t1); >> figure; plot(tnoise) % Calculate autocovariance >> tac = xcorr(tnoise)./80; >> figure; plot(abs(tac)) % Calculate noise power
6 >> Tac = fftshift(fft(fftshift(tac))); >> figure; plot(abs(tac)) Figure 6: The noise process in time appears to be colored and has stronger components and specific frequencies. (a) Raw noise in time (b) Noise power using Rect/Boxcar Window To further analyze, try windowing the original signal and repeating the noise power calculation. % Apply a Kaiser window >> W = kaiser(80,8); >> tac = xcorr(tnoise.*w')./80; >> figure; plot(abs(tac)) % Calculate the noise power of the windowed function >> tac = xcorr(tnoise)./80; >> Tac = fftshift(fft(fftshift(tac))); >> figure; plot(-(79/160):(1/160):(79/160),abs(tac)') Figure 7: The noise process in the windowed function identifies a prominent peak at 0.25 Hz suggestive of a small periodicity of the noise in time. This is something we need to be aware of when cross-correlating with the paradigm function. The paradigm is at much lower frequency than this small periodicity, so we will ignore it in this problem.
7 c. Calculate the cross-correlation function for each voxel in the data set using the xcorr() function in Matlab. Judicious use of image thresholding may be used to reduce computation time. Your result should be an image that represents a 2D map of correlation coefficients in the brain. Figure 8: Resulting 2D correlation map showing region of activation in the posterior (occipital) lobe. Figure 9: After applying a Kaiser window Figure 10 d. Repeat part c. with the use of a satisfactory window of your choice (see window() function in Matlab). Does windowing improve the quality of your results? I have to admit that I intended for the temporal windowing to improve the results. However, it is clear from the figures above that the application of the Kaiser window, which may lead to more temporal noise correlations, degrades the intensity of the activated region relative to the background. e. Repeat part c. after applying a 2D low pass filter to the image data. Does filtering improve the quality of your results? This clearly improves the results although there is clearly degradation in the spatial resolution. f. No activation is expected in the frontal lobe of the brain. Therefore, we can use the region as a representative of the background correlation. For you best correlation result, use hypothesis testing (i.e. T-test) to select the threshold for activated voxels assuming the null hypothesis, H o, and a twotailed test with α = 0.05.
8 Now we use the concepts of hypothesis testing to determine which voxels are significant. Where BW is a mask that identifies regions in the frontal lobe versus the center of the region of activation. >> mean(corrmap(find(double(bw).*corrmap))) ans = >> std(corrmap(find(double(bw).*corrmap))) ans = for i = 1:128*128, if(((corrmap(i) )/0.3014)<1.99) % 80-1 degrees of freedom corrmap(i) = 0; end end >> figure; imagesc(1000.*corrmap); axis('image');colorbar Figure 11: Region of significant activation.
Signal 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 informationImage Quality/Artifacts Frequency (MHz)
The Larmor Relation 84 Image Quality/Artifacts (MHz) 42 ω = γ X B = 2πf 84 0.0 1.0 2.0 Magnetic Field (Tesla) 1 A 1D Image Magnetic Field Gradients Magnet Field Strength Field Strength / Gradient Coil
More information2014 M.S. Cohen all rights reserved
2014 M.S. Cohen all rights reserved mscohen@g.ucla.edu IMAGE QUALITY / ARTIFACTS SYRINGOMYELIA Source http://gait.aidi.udel.edu/res695/homepage/pd_ortho/educate/clincase/syrsco.htm Surgery is usually recommended
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 informationFirst-level fmri modeling. UCLA Advanced NeuroImaging Summer School, 2010
First-level fmri modeling UCLA Advanced NeuroImaging Summer School, 2010 Task on Goal in fmri analysis Find voxels with BOLD time series that look like this Delay of BOLD response Voxel with signal Voxel
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 informationEE469B: Assignment 4 Solutions
EE469B Fall 26-7 RF Pulse Design for MRI EE469B: Assignment 4 Solutions Due Thursday Oct 27. True Null/Flyback Spectral-Spatial Pulses True null and flyback designs are very closely related. In this problem
More informationHarmonic Analysis. Purpose of Time Series Analysis. What Does Each Harmonic Mean? Part 3: Time Series I
Part 3: Time Series I Harmonic Analysis Spectrum Analysis Autocorrelation Function Degree of Freedom Data Window (Figure from Panofsky and Brier 1968) Significance Tests Harmonic Analysis Harmonic analysis
More informationMRI Summer Course Lab 2: Gradient Echo T1 & T2* Curves
MRI Summer Course Lab 2: Gradient Echo T1 & T2* Curves Experiment 1 Goal: Examine the effect caused by changing flip angle on image contrast in a simple gradient echo sequence and derive T1-curves. Image
More information1 Introduction. 2 The basic principles of NMR
1 Introduction Since 1977 when the first clinical MRI scanner was patented nuclear magnetic resonance imaging is increasingly being used for medical diagnosis and in scientific research and application
More informationCoE4TN4 Image Processing. Chapter 3: Intensity Transformation and Spatial Filtering
CoE4TN4 Image Processing Chapter 3: Intensity Transformation and Spatial Filtering Image Enhancement Enhancement techniques: to process an image so that the result is more suitable than the original image
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 informationPHYS225 Lecture 15. Electronic Circuits
PHYS225 Lecture 15 Electronic Circuits Last lecture Difference amplifier Differential input; single output Good CMRR, accurate gain, moderate input impedance Instrumentation amplifier Differential input;
More informationEE469B: Assignment 1 Solutions
EE469B Fall 26-7 RF Pulse Design for MRI EE469B: Assignment Solutions Due Thursday Oct 6 Introduction This assignment concerns typical Fourier transform designs of excitation pulses. This includes designing
More informationRemoval of Line Noise Component from EEG Signal
1 Removal of Line Noise Component from EEG Signal Removal of Line Noise Component from EEG Signal When carrying out time-frequency analysis, if one is interested in analysing frequencies above 30Hz (i.e.
More informationIntroduction to Computational Neuroscience
Introduction to Computational Neuroscience Lecture 4: Data analysis I Lesson Title 1 Introduction 2 Structure and Function of the NS 3 Windows to the Brain 4 Data analysis 5 Data analysis II 6 Single neuron
More informationPHYSIOLOGICAL DE-NOISING FMRI DATA. Katie Dickerson & Jeff MacInnes February 11th, 2013
PHYSIOLOGICAL DE-NOISING FMRI DATA Katie Dickerson & Jeff MacInnes February 11th, 2013 OUTLINE OUTLINE Theoretical overview OUTLINE Theoretical overview OUTLINE Theoretical overview Tutorial in FSL OVERVIEW
More informationEE482: Digital Signal Processing Applications
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu EE482: Digital Signal Processing Applications Spring 2014 TTh 14:30-15:45 CBC C222 Lecture 12 Speech Signal Processing 14/03/25 http://www.ee.unlv.edu/~b1morris/ee482/
More informationPulse Sequence Design and Image Procedures
Pulse Sequence Design and Image Procedures 1 Gregory L. Wheeler, BSRT(R)(MR) MRI Consultant 2 A pulse sequence is a timing diagram designed with a series of RF pulses, gradients switching, and signal readout
More informationTunable Multi Notch Digital Filters A MATLAB demonstration using real data
Tunable Multi Notch Digital Filters A MATLAB demonstration using real data Jon Bell CSIRO ATNF 27 Sep 2 1 Introduction Many people are investigating a wide range of interference suppression techniques.
More informationSupplementary Material
Supplementary Material Orthogonal representation of sound dimensions in the primate midbrain Simon Baumann, Timothy D. Griffiths, Li Sun, Christopher I. Petkov, Alex Thiele & Adrian Rees Methods: Animals
More information(i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters
FIR Filter Design Chapter Intended Learning Outcomes: (i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters (ii) Ability to design linear-phase FIR filters according
More informationDesign Digital Non-Recursive FIR Filter by Using Exponential Window
International Journal of Emerging Engineering Research and Technology Volume 3, Issue 3, March 2015, PP 51-61 ISSN 2349-4395 (Print) & ISSN 2349-4409 (Online) Design Digital Non-Recursive FIR Filter by
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 informationMethods. Experimental Stimuli: We selected 24 animals, 24 tools, and 24
Methods Experimental Stimuli: We selected 24 animals, 24 tools, and 24 nonmanipulable object concepts following the criteria described in a previous study. For each item, a black and white grayscale photo
More informationSUPPLEMENTARY INFORMATION
Supplementary Information S1. Theory of TPQI in a lossy directional coupler Following Barnett, et al. [24], we start with the probability of detecting one photon in each output of a lossy, symmetric beam
More information(i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters
FIR Filter Design Chapter Intended Learning Outcomes: (i) Understanding of the characteristics of linear-phase finite impulse response (FIR) filters (ii) Ability to design linear-phase FIR filters according
More informationExamples of image processing
Examples of image processing Example 1: We would like to automatically detect and count rings in the image 3 Detection by correlation Correlation = degree of similarity Correlation between f(x, y) and
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 informationCHAPTER 2 FIR ARCHITECTURE FOR THE FILTER BANK OF SPEECH PROCESSOR
22 CHAPTER 2 FIR ARCHITECTURE FOR THE FILTER BANK OF SPEECH PROCESSOR 2.1 INTRODUCTION A CI is a device that can provide a sense of sound to people who are deaf or profoundly hearing-impaired. Filters
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 informationSupporting Text Signal Conditioning.
Supporting Text Signal Conditioning. Electrode impedances in physiological saline were typically 1 M! at 10 Hz for both reactive and resistive components. All electrical signals, i.e., those for the mystacial
More informationEE 791 EEG-5 Measures of EEG Dynamic Properties
EE 791 EEG-5 Measures of EEG Dynamic Properties Computer analysis of EEG EEG scientists must be especially wary of mathematics in search of applications after all the number of ways to transform data is
More informationSupplementary Figure 1
Supplementary Figure 1 Left aspl Right aspl Detailed description of the fmri activation during allocentric action observation in the aspl. Averaged activation (N=13) during observation of the allocentric
More informationDesign of FIR Filter for Efficient Utilization of Speech Signal Akanksha. Raj 1 Arshiyanaz. Khateeb 2 Fakrunnisa.Balaganur 3
IJSRD - International Journal for Scientific Research & Development Vol. 3, Issue 03, 2015 ISSN (online): 2321-0613 Design of FIR Filter for Efficient Utilization of Speech Signal Akanksha. Raj 1 Arshiyanaz.
More informationInstruction Manual for Concept Simulators. Signals and Systems. M. J. Roberts
Instruction Manual for Concept Simulators that accompany the book Signals and Systems by M. J. Roberts March 2004 - All Rights Reserved Table of Contents I. Loading and Running the Simulators II. Continuous-Time
More informationThis content has been downloaded from IOPscience. Please scroll down to see the full text.
This content has been downloaded from IOPscience. Please scroll down to see the full text. Download details: IP Address: 148.251.232.83 This content was downloaded on 10/07/2018 at 03:39 Please note that
More informationSignal segmentation and waveform characterization. Biosignal processing, S Autumn 2012
Signal segmentation and waveform characterization Biosignal processing, 5173S Autumn 01 Short-time analysis of signals Signal statistics may vary in time: nonstationary how to compute signal characterizations?
More informationSUPPLEMENTARY INFORMATION
a b STS IOS IOS STS c "#$"% "%' STS posterior IOS dorsal anterior ventral d "( "& )* e f "( "#$"% "%' "& )* Supplementary Figure 1. Retinotopic mapping of the non-lesioned hemisphere. a. Inflated 3D representation
More informationPractice 2. Baseband Communication
PRACTICE : Practice. Baseband Communication.. Objectives To learn to use the software Simulink of MATLAB so as to analyze baseband communication systems... Practical development... Unipolar NRZ signal
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 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 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 informationThe Scientist and Engineer's Guide to Digital Signal Processing By Steven W. Smith, Ph.D.
The Scientist and Engineer's Guide to Digital Signal Processing By Steven W. Smith, Ph.D. Home The Book by Chapters About the Book Steven W. Smith Blog Contact Book Search Download this chapter in PDF
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 informationList and Description of MATLAB Script Files. add_2(n1,n2,b), n1 and n2 are data samples to be added with b bits of precision.
List and Description of MATLAB Script Files 1. add_2(n1,n2,b) add_2(n1,n2,b), n1 and n2 are data samples to be added with b bits of precision. Script file forms sum using 2-compl arithmetic with b bits
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 informationEENG473 Mobile Communications Module 3 : Week # (12) Mobile Radio Propagation: Small-Scale Path Loss
EENG473 Mobile Communications Module 3 : Week # (12) Mobile Radio Propagation: Small-Scale Path Loss Introduction Small-scale fading is used to describe the rapid fluctuation of the amplitude of a radio
More informationijdsp Workshop: Exercise 2012 DSP Exercise Objectives
Objectives DSP Exercise The objective of this exercise is to provide hands-on experiences on ijdsp. It consists of three parts covering frequency response of LTI systems, pole/zero locations with the frequency
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 informationTime and Frequency Domain Windowing of LFM Pulses Mark A. Richards
Time and Frequency Domain Mark A. Richards September 29, 26 1 Frequency Domain Windowing of LFM Waveforms in Fundamentals of Radar Signal Processing Section 4.7.1 of [1] discusses the reduction of time
More informationFFT 1 /n octave analysis wavelet
06/16 For most acoustic examinations, a simple sound level analysis is insufficient, as not only the overall sound pressure level, but also the frequency-dependent distribution of the level has a significant
More informationWindows Connections. Preliminaries
Windows Connections Dale B. Dalrymple Next Annual comp.dsp Conference 21425 Corrections Preliminaries The approach in this presentation Take aways Window types Window relationships Windows tables of information
More informationEstimation of cross coupling of receiver noise between the EoR fat-dipole antennas
Estimation of cross coupling of receiver noise between the EoR fat-dipole antennas Due to the proximity of the fat dipoles in the EoR receiver configuration, the receiver noise of individual antennas may
More informationPREDICTION OF FINGER FLEXION FROM ELECTROCORTICOGRAPHY DATA
University of Tartu Institute of Computer Science Course Introduction to Computational Neuroscience Roberts Mencis PREDICTION OF FINGER FLEXION FROM ELECTROCORTICOGRAPHY DATA Abstract This project aims
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 informationDesign of FIR Filters
Design of FIR Filters Elena Punskaya www-sigproc.eng.cam.ac.uk/~op205 Some material adapted from courses by Prof. Simon Godsill, Dr. Arnaud Doucet, Dr. Malcolm Macleod and Prof. Peter Rayner 1 FIR as a
More informationExperiment 2 Effects of Filtering
Experiment 2 Effects of Filtering INTRODUCTION This experiment demonstrates the relationship between the time and frequency domains. A basic rule of thumb is that the wider the bandwidth allowed for the
More informationEMBEDDED DOPPLER ULTRASOUND SIGNAL PROCESSING USING FIELD PROGRAMMABLE GATE ARRAYS
EMBEDDED DOPPLER ULTRASOUND SIGNAL PROCESSING USING FIELD PROGRAMMABLE GATE ARRAYS Diaa ElRahman Mahmoud, Abou-Bakr M. Youssef and Yasser M. Kadah Biomedical Engineering Department, Cairo University, Giza,
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 informationSIEMENS MAGNETOM Skyra syngo MR D13
Page 1 of 12 SIEMENS MAGNETOM Skyra syngo MR D13 \\USER\CIND\StudyProtocols\PTSA\*ep2d_M0Map_p2_TE15 TA:7.9 s PAT:2 Voxel size:2.5 2.5 3.0 mm Rel. SNR:1.00 :epfid Properties Routine Contrast Prio Recon
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 informationChapter 4 SPEECH ENHANCEMENT
44 Chapter 4 SPEECH ENHANCEMENT 4.1 INTRODUCTION: Enhancement is defined as improvement in the value or Quality of something. Speech enhancement is defined as the improvement in intelligibility and/or
More informationspeech signal S(n). This involves a transformation of S(n) into another signal or a set of signals
16 3. SPEECH ANALYSIS 3.1 INTRODUCTION TO SPEECH ANALYSIS Many speech processing [22] applications exploits speech production and perception to accomplish speech analysis. By speech analysis we extract
More informationThe Polyphase Filter Bank Technique
CASPER Memo 41 The Polyphase Filter Bank Technique Jayanth Chennamangalam Original: 2011.08.06 Modified: 2014.04.24 Introduction to the PFB In digital signal processing, an instrument or software that
More information21/01/2014. Fundamentals of the analysis of neuronal oscillations. Separating sources
21/1/214 Separating sources Fundamentals of the analysis of neuronal oscillations Robert Oostenveld Donders Institute for Brain, Cognition and Behaviour Radboud University Nijmegen, The Netherlands Use
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 informationNoise estimation and power spectrum analysis using different window techniques
IOSR Journal of Electrical and Electronics Engineering (IOSR-JEEE) e-issn: 78-1676,p-ISSN: 30-3331, Volume 11, Issue 3 Ver. II (May. Jun. 016), PP 33-39 www.iosrjournals.org Noise estimation and power
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 informationInitial ARGUS Measurement Results
Initial ARGUS Measurement Results Grant Hampson October 8, Introduction This report illustrates some initial measurement results from the new ARGUS system []. Its main focus is on simple measurements of
More informationA Closer Look at 2-Stage Digital Filtering in the. Proposed WIDAR Correlator for the EVLA
NRC-EVLA Memo# 1 A Closer Look at 2-Stage Digital Filtering in the Proposed WIDAR Correlator for the EVLA NRC-EVLA Memo# Brent Carlson, June 2, 2 ABSTRACT The proposed WIDAR correlator for the EVLA that
More informationPrewhitening. 1. Make the ACF of the time series appear more like a delta function. 2. Make the spectrum appear flat.
Prewhitening What is Prewhitening? Prewhitening is an operation that processes a time series (or some other data sequence) to make it behave statistically like white noise. The pre means that whitening
More informationQuestion 1 Draw a block diagram to illustrate how the data was acquired. Be sure to include important parameter values
Data acquisition Question 1 Draw a block diagram to illustrate how the data was acquired. Be sure to include important parameter values The block diagram illustrating how the signal was acquired is shown
More informationDigital Image Processing 3/e
Laboratory Projects for Digital Image Processing 3/e by Gonzalez and Woods 2008 Prentice Hall Upper Saddle River, NJ 07458 USA www.imageprocessingplace.com The following sample laboratory projects are
More informationEECE 301 Signals & Systems Prof. Mark Fowler
EECE 301 Signals & Systems Prof. Mark Fowler Note Set #16 C-T Signals: Using FT Properties 1/12 Recall that FT Properties can be used for: 1. Expanding use of the FT table 2. Understanding real-world concepts
More informationImage Enhancement in the Spatial Domain (Part 1)
Image Enhancement in the Spatial Domain (Part 1) Lecturer: Dr. Hossam Hassan Email : hossameldin.hassan@eng.asu.edu.eg Computers and Systems Engineering Principle Objective of Enhancement Process an image
More informationApplication of pulse compression technique to generate IEEE a-compliant UWB IR pulse with increased energy per bit
Application of pulse compression technique to generate IEEE 82.15.4a-compliant UWB IR pulse with increased energy per bit Tamás István Krébesz Dept. of Measurement and Inf. Systems Budapest Univ. of Tech.
More informationSimulation Based Design Analysis of an Adjustable Window Function
Journal of Signal and Information Processing, 216, 7, 214-226 http://www.scirp.org/journal/jsip ISSN Online: 2159-4481 ISSN Print: 2159-4465 Simulation Based Design Analysis of an Adjustable Window Function
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 informationInternational Journal of Modern Trends in Engineering and Research e-issn No.: , Date: 2-4 July, 2015
International Journal of Modern Trends in Engineering and Research www.ijmter.com e-issn No.:2349-9745, Date: 2-4 July, 2015 Analysis of Speech Signal Using Graphic User Interface Solly Joy 1, Savitha
More informationWhat is image enhancement? Point operation
IMAGE ENHANCEMENT 1 What is image enhancement? Image enhancement techniques Point operation 2 What is Image Enhancement? Image enhancement is to process an image so that the result is more suitable than
More informationFIR window method: A comparative Analysis
IOSR Journal of Electronics and Communication Engineering (IOSR-JECE) e-issn: 2278-2834,p- ISSN: 2278-8735.Volume 1, Issue 4, Ver. III (Jul - Aug.215), PP 15-2 www.iosrjournals.org FIR window method: A
More informationA Closer Look at 2-Stage Digital Filtering in the. Proposed WIDAR Correlator for the EVLA. NRC-EVLA Memo# 003. Brent Carlson, June 29, 2000 ABSTRACT
MC GMIC NRC-EVLA Memo# 003 1 A Closer Look at 2-Stage Digital Filtering in the Proposed WIDAR Correlator for the EVLA NRC-EVLA Memo# 003 Brent Carlson, June 29, 2000 ABSTRACT The proposed WIDAR correlator
More informationEncoding a Hidden Digital Signature onto an Audio Signal Using Psychoacoustic Masking
The 7th International Conference on Signal Processing Applications & Technology, Boston MA, pp. 476-480, 7-10 October 1996. Encoding a Hidden Digital Signature onto an Audio Signal Using Psychoacoustic
More informationExam in 1TT850, 1E275. Modulation, Demodulation and Coding course
Exam in 1TT850, 1E275 Modulation, Demodulation and Coding course EI, TF, IT programs 16th of August 2004, 14:00-19:00 Signals and systems, Uppsala university Examiner Sorour Falahati office: 018-471 3071
More informationCorrelation, Interference. Kalle Ruttik Department of Communications and Networking School of Electrical Engineering Aalto University
Correlation, Interference Kalle Ruttik Department of Communications and Networking School of Electrical Engineering Aalto University Correlation Correlation Digital communication uses extensively signals
More informationLecture 2 Exercise 1a. Lecture 2 Exercise 1b
Lecture 2 Exercise 1a 1 Design a converter that converts a speed of 60 miles per hour to kilometers per hour. Make the following format changes to your blocks: All text should be displayed in bold. Constant
More informationPulse Compression. Since each part of the pulse has unique frequency, the returns can be completely separated.
Pulse Compression Pulse compression is a generic term that is used to describe a waveshaping process that is produced as a propagating waveform is modified by the electrical network properties of the transmission
More informationPulse Sequence Design Made Easier
Pulse Sequence Design Made Easier Gregory L. Wheeler, BSRT(R)(MR) MRI Consultant gurumri@gmail.com 1 2 Pulse Sequences generally have the following characteristics: An RF line characterizing RF Pulse applications
More informationLab 8. Signal Analysis Using Matlab Simulink
E E 2 7 5 Lab June 30, 2006 Lab 8. Signal Analysis Using Matlab Simulink Introduction The Matlab Simulink software allows you to model digital signals, examine power spectra of digital signals, represent
More informationVISUAL NEURAL SIMULATOR
VISUAL NEURAL SIMULATOR Tutorial for the Receptive Fields Module Copyright: Dr. Dario Ringach, 2015-02-24 Editors: Natalie Schottler & Dr. William Grisham 2 page 2 of 38 3 Introduction. The goal of this
More informationADAPTIVE ANTENNAS. TYPES OF BEAMFORMING
ADAPTIVE ANTENNAS TYPES OF BEAMFORMING 1 1- Outlines This chapter will introduce : Essential terminologies for beamforming; BF Demonstrating the function of the complex weights and how the phase and amplitude
More informationImage Enhancement in Spatial Domain
Image Enhancement in Spatial Domain 2 Image enhancement is a process, rather a preprocessing step, through which an original image is made suitable for a specific application. The application scenarios
More information25 CP Generalize Concepts in Abstract Multi-dimensional Image Model Component Semantics Page 1
25 CP-1390 - Generalize Concepts in Abstract Multi-dimensional Image Model Component Semantics Page 1 1 STATUS Letter Ballot 2 Date of Last Update 2014/09/08 3 Person Assigned David Clunie 4 mailto:dclunie@dclunie.com
More information1 PeZ: Introduction. 1.1 Controls for PeZ using pezdemo. Lab 15b: FIR Filter Design and PeZ: The z, n, and O! Domains
DSP First, 2e Signal Processing First Lab 5b: FIR Filter Design and PeZ: The z, n, and O! Domains The lab report/verification will be done by filling in the last page of this handout which addresses a
More informationLakehead University. Department of Electrical Engineering
Lakehead University Department of Electrical Engineering Lab Manual Engr. 053 (Digital Signal Processing) Instructor: Dr. M. Nasir Uddin Last updated on January 16, 003 1 Contents: Item Page # Guidelines
More informationUsing the DFT as a Filter: Correcting a Misconception by Richard G. Lyons
Using the DFT as a Filter: Correcting a Misconception by Richard G. Lyons I have read, in some of the literature of DSP, that when the discrete Fourier transform (DFT) is used as a filter the process of
More informationCapacitive MEMS accelerometer for condition monitoring
Capacitive MEMS accelerometer for condition monitoring Alessandra Di Pietro, Giuseppe Rotondo, Alessandro Faulisi. STMicroelectronics 1. Introduction Predictive maintenance (PdM) is a key component of
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 informationCommunications over Sparse Channels:
Communications over Sparse Channels: Fundamental limits and practical design Phil Schniter (With support from NSF grant CCF-1018368, NSF grant CCF-1218754, and DARPA/ONR grant N66001-10-1-4090) Intl. Zürich
More information