EECS 556 Winter 2007 Exam2 Solutions
|
|
- Diane McCormick
- 5 years ago
- Views:
Transcription
1 EECS 556 Winter 2007 Exam2 Solutions 1. [10pts] From Wiener filter derivation you need P g and P fg. g = x **b 2 + v = (f** b 1 + u)** b 2 + v = f** b 1 ** b 2 + u** b 2 +v, [n,m] dropped for convenience. P g = P f B 1 2 B P u B P v, (ω x,ω y ) dropped for convenience. R fg [n,m] = E{ f[n,m]g * [0,0] } = E{f[n,m](f**b 1 **b 2 ) * [0,0] + f[n,m](u** b 2 ) * [0,0] + f[n,m]v * [0,0]} = b 1 * ** b 2 * ** R f [n,m], thus P fg = B 1 * B 2 * P f. H = P fg / P g = { B 1 * B 2 * P f } / { P f B 1 2 B P u B P v }, (ω x,ω y ) dropped for convenience. Most of you did cascaded approach, which yields the same answer. Some of you had P u B 1 2 denominator (-2pts). in the 2. (a) [5pts] MSE = (NN), (linear), (cubic). Cubic interpolation works best for the cameraman image (img1). (b) [5pts] MSE = (NN), (linear), (cubic). Linear interpolation works best for this image (img2). (c) [10pts] The answer is NO. For part (b) linear is better than cubic. If you inspect the shape of the second image, you can observe linearly increasing intensity from the corners towards the center, a cone shape if you will. Upsampling via interpolation is filling in the missing intermediate values. As the intensities are varying linearly, the best interpolator to use is the linear interpolator for the second image. Cubic interpolator is also capable of interpolating linearly varying intensities. Under noisy conditions, cubic interpolator may deviate from a perfect line while linear interpolator always produces a straight line. You cannot generalize this result to other images. In general, higher order interpolator (popular choice is cubic) is likely to be better than lower order interpolator.
2 % exam2 prob2 clear; load e2p2.mat img = img1; img = img2; sub_img = sub_img1; sub_img = sub_img2; figure; subplot(331);imagesc(img);colormap gray; axis image;colorbar;title('original Image'); subplot(332);imagesc(sub_img);colormap gray; axis image;colorbar;title('downsampled Image'); interp_img1 = interp2(x, y, sub_img, x_int, y_int, 'nearest'); mse1 = mean2( (interp_img1 - img).^2 ); subplot(334);imagesc(interp_img1);colormap gray; axis image;colorbar;title('nn interp'); subplot(337);imagesc(interp_img1 - img);colormap gray; axis image;colorbar;title('difference'); title(sprintf('nn interp error, MSE=%g', mse1)) interp_img2 = interp2(x, y, sub_img, x_int, y_int, 'linear'); mse2 = mean2( (interp_img2 - img).^2 ); subplot(335);imagesc(interp_img2);colormap gray; axis image;colorbar;title('linear interp'); subplot(338);imagesc(interp_img2 - img);colormap gray; axis image;colorbar;title('difference'); title(sprintf('linear interp error, MSE=%g', mse2))
3 interp_img3 = interp2(x, y, sub_img, x_int, y_int, 'cubic'); mse3 = mean2( (interp_img3 - img).^2 ); subplot(336);imagesc(interp_img3);colormap gray; axis image;colorbar;title('cubic interp'); subplot(339);imagesc(interp_img3 - img);colormap gray; axis image;colorbar;title('difference'); title(sprintf('cubic interp error, MSE=%g', mse3)) 3. (a) [5pts] As per template, take yy(1:10,1:10) and compute sample variance. MSE = See the figures for the filtered output; <estimated noise, no blur> (b) [5pts] MSE = See the figures for the filtered output, <true noise, no blur> (c) [5pts] MSE = See the figures for the filtered output, <true noise, estimated blur> (d) [5pts] MSE = See the figures for the filtered output, < true noise, true blur> (e) [5pts] In part (c), blur is estimated through periodogram (P y ), which is known to be very noisy. Other components of the formula, P s and P v, are derived analytically so there are no errors in those components. If you plot B(ω x,ω y ) and estimated blur, B_est(ω x,ω y ), then you will observe that estimated blur is quite noisy. Using this noisy blur function led to worsening of the MSE in part (c) compared to part (b) where there is no blur. Basically, you are better off using no information (i.e., no blur) than using unreliable information (i.e., noisy estimated blur). The other reason is that the true blur is a 3x3 filter with high center value which resembles a delta function a lot. Thus, assuming that the blur is a delta function is quite reasonable. In general, estimating a blur function is a difficult problem to tackle. Some of you kindly suggested that using windowed periodogram helps, which is a good suggestion.
4 % exam2 prob3 % wiener filter problem with blur % reused code from hw7 prob5 rand('state', 0); randn('state', 0); nx = 62; ny = 50; h1 = zeros(27,27); h1(10:18, 11:17) = 1; p = 0.004; ss_tmp = (rand(nx,ny) < p) - p; ss = conv2(ss_tmp, h1, 'same') / sqrt(p*(1-p)); % random rect signal blur = 1/9*[0 1 0;1 5 1;0 1 0]; yy = conv2(ss, blur, 'same') + 3 * randn(size(ss)); % add gaussian noise - what variance? px = 2*nx; py = 2*ny; % zero padding to reduce space aliasing n1 = size(h1,1)-1; ix = [-n1:n1]; n2 = size(h1,2)-1; iy = [-n2:n2]; Rs = zeros(px,py); % true autocorr func h2 = xcorr2(h1, h1); Rs(ix + end/2+1,iy + end/2+1) = h2; wx = [-px/2:px/2-1]'*2*pi/px; wy = [-py/2:py/2-1]'*2*pi/py; Ps = fftshift ( (Diric(wx,9) * Diric(wy,7)').^2 ); % Method 1: analytical power spectrum Ps = abs(fft2(h1, px, py)).^2; % Method 2: H_1 (the simplest here) % estimate noise Pv from background sub_yy = yy(1:10,1:10); Pv_est = var(sub_yy(:)); % estimate blur Py = abs( fft2( xcorr2(yy,yy)/(nx*ny), px, py) ); Py2 = Py -9; Py2(Py2 < 0) = 0; B_est = sqrt(py2./ps); % B_est is real valued between 0 and 1 B_est(B_est > 1) = 1; B_est(B_est < 0) = 0; B = fftshift( (1/9)*(5 + 2*cos(wx)*ones(1,py) + 2*ones(px,1)*cos(wy)') ); % analytical method B = abs(fft2(blur, px, py)); % magnitude of fft of blur H1 = Ps./ (Ps + Pv_est); % Wiener filter wo blur w estimated noise H2 = Ps./ (Ps + 9); % Wiener filter wo blur w true noise H3 = (Ps.*B_est)./ (Ps.*(B_est.^2) + 9); % Wiener filter with estimated blur H4 = (Ps.*B)./ (Ps.*(B.^2) + 9); % Wiener filter with true blur
5 Y = fft2(yy, px, py); s_hat1 = real(ifft2(y.* H1)); s_hat1 = s_hat1(1:nx,1:ny); % un-zero-pad s_hat2 = real(ifft2(y.* H2)); s_hat2 = s_hat2(1:nx,1:ny); % un-zero-pad s_hat3 = real(ifft2(y.* H3)); s_hat3 = s_hat3(1:nx,1:ny); % un-zero-pad s_hat4 = real(ifft2(y.* H4)); s_hat4 = s_hat4(1:nx,1:ny); % un-zero-pad clf, pl = 230; colormap(1-gray(256)), subplot(pl+1), imagesc(1:nx,1:ny,ss'), axis xy, axis image xlabel n, ylabel m, title 'True: s(n,m)', colorbar horiz subplot(pl+2), imagesc(1:nx,1:ny,yy'), axis xy, axis image xlabel n, ylabel m, title 'Noisy: y(n,m)', colorbar horiz subplot(pl+3), imagesc(1:nx, 1:ny, s_hat1'), axis xy, axis image xlabel n, ylabel m, colorbar horiz wiener_mse = mean2((s_hat1-ss).^2); title(sprintf('no blur, est. noise, s.hat1(n,m), MSE=%g', wiener_mse)) subplot(pl+4), imagesc(1:nx, 1:ny, s_hat2'), axis xy, axis image xlabel n, ylabel m, colorbar horiz wiener_mse = mean2((s_hat2-ss).^2); title(sprintf('no blur, true noise, s.hat2(n,m), MSE=%g', wiener_mse)) subplot(pl+5), imagesc(1:nx, 1:ny, s_hat3'), axis xy, axis image xlabel n, ylabel m, colorbar horiz wiener_mse = mean2((s_hat3-ss).^2); title(sprintf('est blur, true noise, s.hat3(n,m), MSE=%g', wiener_mse)) subplot(pl+6), imagesc(1:nx, 1:ny, s_hat3'), axis xy, axis image xlabel n, ylabel m, colorbar horiz wiener_mse = mean2((s_hat4-ss).^2); title(sprintf('true blur, true noise, s.hat4(n,m), MSE=%g', wiener_mse))
6 4. (a) [5pts] NPLS (1 st order penalty), MSE= 2.399, beta = 2 3 = 8. See plots for filtered output and difference image. (b) [5pts] NPLS2 (2 nd order penalty), MSE= 3.879, beta = 2 2 = 4. See plots for filtered output and difference image. (C) [10pts] Combined block by block approach, the left potion of the image has piecewise constant (step edges) signal, which incurs no penalty under NPLS. Thus, you choose NPLS for the left potion of the image. The right potion has linearly increasing shape, which incurs no penalty under NPLS2. Thus, you choose NPLS2 for the right potion of the image. MSE= 2.145, beta1 (NPLS, left potion) = 2 2 = 4, beta2 (NPLS2, right potion) = 2 5 = 32. See plots for filtered output and difference image. You need to vary beta1 with NPLS to minimize MSE of just the left portion of the image and repeat the same process with NPLS2 to find out beta2 for the right potion of the image. Above was the original intended answer applying different NPSL depending on the profile of the signal. However, if you apply NPLS to both potions beta1 = 4 (NPLS, left potion) and beta2 = 16 (NPLS, right potion) your will actually get lower MSE = This happened because you are restricted to beta values between If you are allowed larger beta values minimum MSE = 2.12 occurs at beta1 = 4 (NPLS, left potion), beta2 = 2 16 or above (NPLS2, right potion). Answers applying NPLS to both potions are also correct. In this case, the rationale for choosing NPLS for both potions is to minimize MSE. Grading: both answers allowed. You need to provide justification for choosing NPLS/NPLS2 (signal profile dependent) or NPLS/NPLS (minimize MSE) combination. Not having the justification -3pts.
7 % exam2 prob4 % adaptive NPLS % code reused from hw7 and hw9 rand('state', 0); randn('state', 0); % 1st part of the signal nx = 62; ny = 50; h1 = zeros(27,27); h1(10:18, 11:17) = 1; p = 0.003; ss_tmp = (rand(nx,ny) < p) - p; ss = conv2(ss_tmp, h1, 'same') / sqrt(p*(1-p)); % random rect signal % 2nd part of the signal xtrue = zeros(nx,ny); ix = -(nx-1)/2:(nx-1)/2; iy = -(ny-1)/2:(ny-1)/2;
8 [ix, iy] = ndgrid(ix, iy); xtrue = (1 - min(abs(ix/(nx/3)),1)) * 30; xtrue = (1 - min(abs(iy/(ny/3)),1)).* xtrue; xtrue = xtrue - mean2(xtrue) + mean2(ss); nx2 = 2*nx; ny2 = ny; ss2 = [ss; xtrue]; yy = ss2 + 5 * randn(size(ss2)); subplot(421), imagesc(ss2',[-15 35]), axis xy, axis image, colormap gray xlabel n, ylabel m, title 'True: ss2(n,m)', colorbar subplot(422), imagesc(yy',[-15 35]), axis xy, axis image, xlabel n, ylabel m, title 'Noisy: yy(n,m)', colorbar % apply NPLS niter=200; delta=1; % fixed beta1 = 2^3; % play with this number xhat1 = npls_sps(yy, niter, beta1, delta); % apply NPLS2 beta2 = 2^2; % play with this number xhat2 = npls2_sps(yy, niter, beta2, delta); subplot(423), imagesc(xhat1',[-15 35]), axis xy, axis image xlabel n, ylabel m, colorbar, title( sprintf( 'NPLS output, beta=%g', beta1) ) subplot(424), imagesc((xhat1-ss2)',[-10 10]), axis xy, axis image xlabel n, ylabel m, colorbar title(sprintf('npls error, MSE=%g', mean2((xhat1-ss2).^2))) subplot(425), imagesc(xhat2',[-15 35]), axis xy, axis image xlabel n, ylabel m, colorbar, title( sprintf( 'NPLS2 output, beta=%g', beta2) ) subplot(426), imagesc((xhat2-ss2)',[-10 10]), axis xy, axis image xlabel n, ylabel m, colorbar title(sprintf('npls2 error, MSE=%g', mean2((xhat2-ss2).^2))) % apply NPLS and NPLS2 block by block yy_part1 = yy(1:nx,1:ny2); yy_part2 = yy(nx+1:2*nx,1:ny2); beta1 = 2^2; beta2 = 2^5; % play with this number xhat3_part1 = npls_sps(yy_part1, niter, beta1, delta); xhat3_part2 = npls2_sps(yy_part2, niter, beta2, delta);
9 xhat3 = [xhat3_part1; xhat3_part2]; subplot(427), imagesc(xhat3',[-15 35]), axis xy, axis image xlabel n, ylabel m, colorbar, title( sprintf( 'combined output, beta1=%g beta2=%g', beta1, beta2) ) subplot(428), imagesc((xhat3-ss2)',[-10 10]), axis xy, axis image xlabel n, ylabel m, colorbar title(sprintf('combined NPLS error, MSE=%g', mean2((xhat3-ss2).^2)))
Digital Image and Fourier Transform
Lab 5 Numerical Methods TNCG17 Digital Image and Fourier Transform Sasan Gooran (Autumn 2009) Before starting this lab you are supposed to do the preparation assignments of this lab. All functions and
More informationHandout 1 - Introduction to plots in Matlab 7
SPHSC 53 Speech Signal Processing UW Summer 6 Handout - Introduction to plots in Matlab 7 Signal analysis is an important part of signal processing. And signal analysis is not complete without signal visualization.
More information4.4 The FFT and MATLAB
4.4. THE FFT AND MATLAB 69 4.4 The FFT and MATLAB 4.4.1 The FFT and MATLAB MATLAB implements the Fourier transform with the following functions: fft, ifft, fftshift, ifftshift, fft2, ifft2. We describe
More informationThe Effect of Time-Domain Interpolation on Response Spectral Calculations. David M. Boore
The Effect of Time-Domain Interpolation on Response Spectral Calculations David M. Boore This note confirms Norm Abrahamson s finding that the straight line interpolation between sampled points used in
More informationUsing Multiple DMs for Increased Spatial Frequency Response
AN: Multiple DMs for Increased Spatial Frequency Response Using Multiple DMs for Increased Spatial Frequency Response AN Author: Justin Mansell Revision: // Abstract Some researchers have come to us suggesting
More informationAdaptive Resampling - Transforming From the Time to the Angle Domain
Adaptive Resampling - Transforming From the Time to the Angle Domain Jason R. Blough, Ph.D. Assistant Professor Mechanical Engineering-Engineering Mechanics Department Michigan Technological University
More informationDELTA MODULATION AND DPCM CODING OF COLOR SIGNALS
DELTA MODULATION AND DPCM CODING OF COLOR SIGNALS Item Type text; Proceedings Authors Habibi, A. Publisher International Foundation for Telemetering Journal International Telemetering Conference Proceedings
More informationUNIVERSITY OF BAHRAIN COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING
UNIVERSITY OF BAHRAIN COLLEGE OF ENGINEERING DEPARTMENT OF ELECTRICAL AND ELECTRONIC ENGINEERING EENG 373: DIGITAL COMMUNICATIONS EXPERIMENT NO. 3 BASEBAND DIGITAL TRANSMISSION Objective This experiment
More informationA NEW LOOK AT FREQUENCY RESOLUTION IN POWER SPECTRAL DENSITY ESTIMATION. Sudeshna Pal, Soosan Beheshti
A NEW LOOK AT FREQUENCY RESOLUTION IN POWER SPECTRAL DENSITY ESTIMATION Sudeshna Pal, Soosan Beheshti Electrical and Computer Engineering Department, Ryerson University, Toronto, Canada spal@ee.ryerson.ca
More informationExtraction Methods of Watermarks from Linearly-Distorted Images to Maximize Signal-to-Noise Ratio. Brandon Migdal. Advisors: Carl Salvaggio
Extraction Methods of Watermarks from Linearly-Distorted Images to Maximize Signal-to-Noise Ratio By Brandon Migdal Advisors: Carl Salvaggio Chris Honsinger A senior project submitted in partial fulfillment
More informationAudio Processing Exercise
Name: Date : Audio Processing Exercise In this exercise you will learn to load, playback, modify, and plot audio files. Commands for loading and characterizing an audio file To load an audio file (.wav)
More informationCalibrate, Characterize and Emulate Systems Using RFXpress in AWG Series
Calibrate, Characterize and Emulate Systems Using RFXpress in AWG Series Introduction System designers and device manufacturers so long have been using one set of instruments for creating digitally modulated
More informationAUTO-FOCUS USING PSD ESTIMATION EFFECTIVE BANDWIDTH By Laurence G. Hassebrook
AUTO-FOCUS USING PSD ESTIMATION EFFECTIVE BANDWIDTH By Laurence G. Hassebrook 2-22-2012 Please follow the tutorial and reproduce the figures with your own code. We demonstrate how to use the effective
More informationSampling Issues in Image and Video
Sampling Issues in Image and Video Spring 06 Instructor: K. J. Ray Liu ECE Department, Univ. of Maryland, College Park Overview and Logistics Last Time: Motion analysis Geometric relations and manipulations
More informationInverse Filtering by Signal Reconstruction from Phase. Megan M. Fuller
Inverse Filtering by Signal Reconstruction from Phase by Megan M. Fuller B.S. Electrical Engineering Brigham Young University, 2012 Submitted to the Department of Electrical Engineering and Computer Science
More informationCopyright Warning & Restrictions
Copyright Warning & Restrictions The copyright law of the United States (Title 17, United States Code) governs the making of photocopies or other reproductions of copyrighted material. Under certain conditions
More informationHands-on session on timing analysis
Amsterdam 2010 Hands-on session on timing analysis Introduction During this session, we ll approach some basic tasks in timing analysis of x-ray time series, with particular emphasis on the typical signals
More informationWE treat the problem of reconstructing a random signal
IEEE TRANSACTIONS ON SIGNAL PROCESSING, VOL. 57, NO. 3, MARCH 2009 977 High-Rate Interpolation of Random Signals From Nonideal Samples Tomer Michaeli and Yonina C. Eldar, Senior Member, IEEE Abstract We
More informationElasticity Imaging with Ultrasound JEE 4980 Final Report. George Michaels and Mary Watts
Elasticity Imaging with Ultrasound JEE 4980 Final Report George Michaels and Mary Watts University of Missouri, St. Louis Washington University Joint Engineering Undergraduate Program St. Louis, Missouri
More information1 Overview. 1.1 Digital Images GEORGIA INSTITUTE OF TECHNOLOGY. ECE 2026 Summer 2018 Lab #5: Sampling: A/D and D/A & Aliasing
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2026 Summer 2018 Lab #5: Sampling: A/D and D/A & Aliasing Date: 21 June 2018 Pre-Lab: You should read the Pre-Lab section
More information1 Overview. 1.1 Digital Images GEORGIA INSTITUTE OF TECHNOLOGY. ECE 2026 Summer 2016 Lab #6: Sampling: A/D and D/A & Aliasing
GEORGIA INSTITUTE OF TECHNOLOGY SCHOOL of ELECTRICAL and COMPUTER ENGINEERING ECE 2026 Summer 2016 Lab #6: Sampling: A/D and D/A & Aliasing Date: 30 June 2016 Pre-Lab: You should read the Pre-Lab section
More informationECE438 - Laboratory 1: Discrete and Continuous-Time Signals
Purdue University: ECE438 - Digital Signal Processing with Applications 1 ECE438 - Laboratory 1: Discrete and Continuous-Time Signals By Prof. Charles Bouman and Prof. Mireille Boutin Fall 2015 1 Introduction
More informationSingle Channel Speech Enhancement Using Spectral Subtraction Based on Minimum Statistics
Master Thesis Signal Processing Thesis no December 2011 Single Channel Speech Enhancement Using Spectral Subtraction Based on Minimum Statistics Md Zameari Islam GM Sabil Sajjad This thesis is presented
More informationAdaptive bilateral filtering of image signals using local phase characteristics
Signal Processing 88 (2008) 1615 1619 Fast communication Adaptive bilateral filtering of image signals using local phase characteristics Alexander Wong University of Waterloo, Canada Received 15 October
More informationECE438 - Laboratory 4: Sampling and Reconstruction of Continuous-Time Signals
Purdue University: ECE438 - Digital Signal Processing with Applications 1 ECE438 - Laboratory 4: Sampling and Reconstruction of Continuous-Time Signals October 6, 2010 1 Introduction It is often desired
More informationSupplemental Material: Color Compatibility From Large Datasets
Supplemental Material: Color Compatibility From Large Datasets Peter O Donovan, Aseem Agarwala, and Aaron Hertzmann Project URL: www.dgp.toronto.edu/ donovan/color/ 1 Unmixing color preferences In the
More informationGenerating Spectrally Rich Data Sets Using Adaptive Band Synthesis Interpolation
Generating Spectrally Rich Data Sets Using Adaptive Band Synthesis Interpolation James C. Rautio Sonnet Software, Inc. WFA: Microwave Component Design Using Optimization Techniques June 2003 Interpolation
More informationSystem Identification
System Identification Arun K. Tangirala Department of Chemical Engineering IIT Madras July 26, 2013 Module 9 Lecture 2 Arun K. Tangirala System Identification July 26, 2013 16 Contents of Lecture 2 In
More informationMultirate Signal Processing: Graphical Representation & Comparison of Decimation & Interpolation Identities using MATLAB
International Journal of Electronics and Communication Engineering. ISSN 0974-2166 Volume 4, Number 4 (2011), pp. 443-452 International Research Publication House http://www.irphouse.com Multirate Signal
More informationMultirate Digital Signal Processing
Multirate Digital Signal Processing Contents 1) What is multirate DSP? 2) Downsampling and Decimation 3) Upsampling and Interpolation 4) FIR filters 5) IIR filters a) Direct form filter b) Cascaded form
More informationAn Introduction to the Spectral Dynamics Rotating Machinery Analysis (RMA) package For PUMA and COUGAR
An Introduction to the Spectral Dynamics Rotating Machinery Analysis (RMA) package For PUMA and COUGAR Introduction: The RMA package is a PC-based system which operates with PUMA and COUGAR hardware to
More informationחלק מהשקפים מעובדים משקפים של פרדו דוראנד, טומס פנקהאוסר ודניאל כהן-אור קורס גרפיקה ממוחשבת 2009/2010 סמסטר א' Image Processing
חלק מהשקפים מעובדים משקפים של פרדו דוראנד, טומס פנקהאוסר ודניאל כהן-אור קורס גרפיקה ממוחשבת 2009/2010 סמסטר א' Image Processing 1 What is an image? An image is a discrete array of samples representing
More informationReconstruction of Ca 2+ dynamics from low frame rate Ca 2+ imaging data CS229 final project. Submitted by: Limor Bursztyn
Reconstruction of Ca 2+ dynamics from low frame rate Ca 2+ imaging data CS229 final project. Submitted by: Limor Bursztyn Introduction Active neurons communicate by action potential firing (spikes), accompanied
More informationDiscrete-time equivalent systems example from matlab: the c2d command
Discrete-time equivalent systems example from matlab: the cd command Create the -time system using the tf command and then generate the discrete-time equivalents using the cd command. We use this command
More informationDesign of Speech Signal Analysis and Processing System. Based on Matlab Gateway
1 Design of Speech Signal Analysis and Processing System Based on Matlab Gateway Weidong Li,Zhongwei Qin,Tongyu Xiao Electronic Information Institute, University of Science and Technology, Shaanxi, China
More informationSignal Stability Analyser
Signal Stability Analyser o Real Time Phase or Frequency Display o Real Time Data, Allan Variance and Phase Noise Plots o 1MHz to 65MHz medium resolution (12.5ps) o 5MHz and 10MHz high resolution (50fs)
More informationDigital Correction for Multibit D/A Converters
Digital Correction for Multibit D/A Converters José L. Ceballos 1, Jesper Steensgaard 2 and Gabor C. Temes 1 1 Dept. of Electrical Engineering and Computer Science, Oregon State University, Corvallis,
More informationFor the SIA. Applications of Propagation Delay & Skew tool. Introduction. Theory of Operation. Propagation Delay & Skew Tool
For the SIA Applications of Propagation Delay & Skew tool Determine signal propagation delay time Detect skewing between channels on rising or falling edges Create histograms of different edge relationships
More informationGPA for DigitalMicrograph
GPA for DigitalMicrograph Geometric Phase Analysis GPA Phase Manual 1.0 HREM Research Inc Conventions The typographic conventions used in this help are described below. Convention Bold Description Used
More informationDigital Lock-In Amplifiers SR850 DSP lock-in amplifier with graphical display
Digital Lock-In Amplifiers SR850 DSP lock-in amplifier with graphical display SR850 DSP Lock-In Amplifier 1 mhz to 102.4 khz frequency range >100 db dynamic reserve 0.001 degree phase resolution Time constants
More informationSpectrum Analyser Basics
Hands-On Learning Spectrum Analyser Basics Peter D. Hiscocks Syscomp Electronic Design Limited Email: phiscock@ee.ryerson.ca June 28, 2014 Introduction Figure 1: GUI Startup Screen In a previous exercise,
More informationTime series analysis
Time series analysis (July 12-13, 2011) Course Exercise Booklet MATLAB function reference 1 Introduction to time series analysis Exercise 1.1 Controlling frequency, amplitude and phase... 3 Exercise 1.2
More informationIntro to DSP: Sampling. with GNU Radio Jeff Long
Intro to DSP: Sampling with GNU Radio Jeff Long ADC SDR Hardware Reconfigurable Logic Front End Analog Bus USB2 USB3 GBE PCI Digital Data Control Analog Signals May include multiplesystem Typical SDR Radio
More informationRegion Adaptive Unsharp Masking based DCT Interpolation for Efficient Video Intra Frame Up-sampling
International Conference on Electronic Design and Signal Processing (ICEDSP) 0 Region Adaptive Unsharp Masking based DCT Interpolation for Efficient Video Intra Frame Up-sampling Aditya Acharya Dept. of
More informationME 565 HW 4 Solutions Winter Make image black and white. Compute the FFT of our image using fft2. clear all; close all; clc %%Exercise 4.
ME 565 HW 4 Solutions Winter 2017 clear all; close all; clc %%Exercise 4.1 %read in the image A= imread('recorder','jpg'); Make image black and white Abw2=rgb2gray(A); [nx,ny]=size(abw2); Compute the FFT
More informationEE369C: Assignment 1
EE369C Fall 17-18 Medical Image Reconstruction 1 EE369C: Assignment 1 Due Wednesday, Oct 4th Assignments This quarter the assignments will be partly matlab, and partly calculations you will need to work
More informationAgilent PN Time-Capture Capabilities of the Agilent Series Vector Signal Analyzers Product Note
Agilent PN 89400-10 Time-Capture Capabilities of the Agilent 89400 Series Vector Signal Analyzers Product Note Figure 1. Simplified block diagram showing basic signal flow in the Agilent 89400 Series VSAs
More informationLabView Exercises: Part II
Physics 3100 Electronics, Fall 2008, Digital Circuits 1 LabView Exercises: Part II The working VIs should be handed in to the TA at the end of the lab. Using LabView for Calculations and Simulations LabView
More informationAppendix D. UW DigiScope User s Manual. Willis J. Tompkins and Annie Foong
Appendix D UW DigiScope User s Manual Willis J. Tompkins and Annie Foong UW DigiScope is a program that gives the user a range of basic functions typical of a digital oscilloscope. Included are such features
More informationMusic Source Separation
Music Source Separation Hao-Wei Tseng Electrical and Engineering System University of Michigan Ann Arbor, Michigan Email: blakesen@umich.edu Abstract In popular music, a cover version or cover song, or
More informationQSched v0.96 Spring 2018) User Guide Pg 1 of 6
QSched v0.96 Spring 2018) User Guide Pg 1 of 6 QSched v0.96 D. Levi Craft; Virgina G. Rovnyak; D. Rovnyak Overview Cite Installation Disclaimer Disclaimer QSched generates 1D NUS or 2D NUS schedules using
More informationVideo Processing Applications Image and Video Processing Dr. Anil Kokaram
Video Processing Applications Image and Video Processing Dr. Anil Kokaram anil.kokaram@tcd.ie This section covers applications of video processing as follows Motion Adaptive video processing for noise
More informationUpgrading E-learning of basic measurement algorithms based on DSP and MATLAB Web Server. Milos Sedlacek 1, Ondrej Tomiska 2
Upgrading E-learning of basic measurement algorithms based on DSP and MATLAB Web Server Milos Sedlacek 1, Ondrej Tomiska 2 1 Czech Technical University in Prague, Faculty of Electrical Engineeiring, Technicka
More informationCourse Web site:
The University of Texas at Austin Spring 2018 EE 445S Real- Time Digital Signal Processing Laboratory Prof. Evans Solutions for Homework #1 on Sinusoids, Transforms and Transfer Functions 1. Transfer Functions.
More informationReproducibility Assessment of Independent Component Analysis of Expression Ratios from DNA microarrays.
Reproducibility Assessment of Independent Component Analysis of Expression Ratios from DNA microarrays. David Philip Kreil David J. C. MacKay Technical Report Revision 1., compiled 16th October 22 Department
More informationIntroduction: Overview. EECE 2510 Circuits and Signals: Biomedical Applications. ECG Circuit 2 Analog Filtering and A/D Conversion
EECE 2510 Circuits and Signals: Biomedical Applications ECG Circuit 2 Analog Filtering and A/D Conversion Introduction: Now that you have your basic instrumentation amplifier circuit running, in Lab ECG1,
More informationVocoder Reference Test TELECOMMUNICATIONS INDUSTRY ASSOCIATION
TIA/EIA STANDARD ANSI/TIA/EIA-102.BABC-1999 Approved: March 16, 1999 TIA/EIA-102.BABC Project 25 Vocoder Reference Test TIA/EIA-102.BABC (Upgrade and Revision of TIA/EIA/IS-102.BABC) APRIL 1999 TELECOMMUNICATIONS
More informationEVALUATION OF SIGNAL PROCESSING METHODS FOR SPEECH ENHANCEMENT MAHIKA DUBEY THESIS
c 2016 Mahika Dubey EVALUATION OF SIGNAL PROCESSING METHODS FOR SPEECH ENHANCEMENT BY MAHIKA DUBEY THESIS Submitted in partial fulfillment of the requirements for the degree of Bachelor of Science in Electrical
More informationIntroduction. Edge Enhancement (SEE( Advantages of Scalable SEE) Lijun Yin. Scalable Enhancement and Optimization. Case Study:
Case Study: Scalable Edge Enhancement Introduction Edge enhancement is a post processing for displaying radiologic images on the monitor to achieve as good visual quality as the film printing does. Edges
More informationInterface Practices Subcommittee SCTE STANDARD SCTE Measurement Procedure for Noise Power Ratio
Interface Practices Subcommittee SCTE STANDARD SCTE 119 2018 Measurement Procedure for Noise Power Ratio NOTICE The Society of Cable Telecommunications Engineers (SCTE) / International Society of Broadband
More informationSingle image super resolution with improved wavelet interpolation and iterative back-projection
IOSR Journal of VLSI and Signal Processing (IOSR-JVSP) Volume 5, Issue 6, Ver. II (Nov -Dec. 2015), PP 16-24 e-issn: 2319 4200, p-issn No. : 2319 4197 www.iosrjournals.org Single image super resolution
More informationAchieve Accurate Critical Display Performance With Professional and Consumer Level Displays
Achieve Accurate Critical Display Performance With Professional and Consumer Level Displays Display Accuracy to Industry Standards Reference quality monitors are able to very accurately reproduce video,
More informationAssessing and Measuring VCR Playback Image Quality, Part 1. Leo Backman/DigiOmmel & Co.
Assessing and Measuring VCR Playback Image Quality, Part 1. Leo Backman/DigiOmmel & Co. Assessing analog VCR image quality and stability requires dedicated measuring instruments. Still, standard metrics
More informationLecture 3, Opamps. Operational amplifiers, high-gain, high-speed
Lecture 3, Opamps Operational amplifiers, high-gain, high-speed What did we do last time? Multi-stage amplifiers Increases gain Increases number of poles Frequency domain Stability Phase margin 86 of 252
More informationDesign Approach of Colour Image Denoising Using Adaptive Wavelet
International Journal of Engineering Research and Development ISSN: 78-067X, Volume 1, Issue 7 (June 01), PP.01-05 www.ijerd.com Design Approach of Colour Image Denoising Using Adaptive Wavelet Pankaj
More informationSOA / PIN based OLT receiver update. David Piehler, Ruomei Mu 17 July 2007
SOA / based OLT receiver update David Piehler, Ruomei Mu 17 July 2007 dpiehler@alphion.com SOA/ OLT receiver New since last time (3av_0705_piehler_1.pdf): Calculations now use same assumptions as 3av_0705_takizawa_1.pdf
More informationProcessing data with Mestrelab Mnova
Processing data with Mestrelab Mnova This exercise has three parts: a 1D 1 H spectrum to baseline correct, integrate, peak-pick, and plot; a 2D spectrum to plot with a 1 H spectrum as a projection; and
More informationStudy of White Gaussian Noise with Varying Signal to Noise Ratio in Speech Signal using Wavelet
American International Journal of Research in Science, Technology, Engineering & Mathematics Available online at http://www.iasir.net ISSN (Print): 2328-3491, ISSN (Online): 2328-3580, ISSN (CD-ROM): 2328-3629
More informationArtifact rejection and running ICA
Artifact rejection and running ICA Task 1 Reject noisy data Task 2 Run ICA Task 3 Plot components Task 4 Remove components (i.e. back-projection) Exercise... Artifact rejection and running ICA Task 1 Reject
More informationE E Introduction to Wavelets & Filter Banks Spring Semester 2009
E E - 2 7 4 Introduction to Wavelets & Filter Banks Spring Semester 29 HOMEWORK 5 DENOISING SIGNALS USING GLOBAL THRESHOLDING One-Dimensional Analysis Using the Command Line This example involves a real-world
More informationPlease feel free to download the Demo application software from analogarts.com to help you follow this seminar.
Hello, welcome to Analog Arts spectrum analyzer tutorial. Please feel free to download the Demo application software from analogarts.com to help you follow this seminar. For this presentation, we use a
More informationAN030. DM Actuator Slaving. AN030: DM Actuator Slaving. Author: Justin Mansell Revision: 6/12/11. 1 Abstract
AN030 DM Actuator Slaving Author: Justin Mansell Revision: 6/12/11 1 Abstract In this application note, we present evidence that slaving of edge actuators creates a reduction in DM performance when using
More informationEE373B Project Report Can we predict general public s response by studying published sales data? A Statistical and adaptive approach
EE373B Project Report Can we predict general public s response by studying published sales data? A Statistical and adaptive approach Song Hui Chon Stanford University Everyone has different musical taste,
More informationAutomatic LP Digitalization Spring Group 6: Michael Sibley, Alexander Su, Daphne Tsatsoulis {msibley, ahs1,
Automatic LP Digitalization 18-551 Spring 2011 Group 6: Michael Sibley, Alexander Su, Daphne Tsatsoulis {msibley, ahs1, ptsatsou}@andrew.cmu.edu Introduction This project was originated from our interest
More informationMITOCW ocw f08-lec19_300k
MITOCW ocw-18-085-f08-lec19_300k The following content is provided under a Creative Commons license. Your support will help MIT OpenCourseWare continue to offer high quality educational resources for free.
More informationANALYSIS OF COMPUTED ORDER TRACKING
Mechanical Systems and Signal Processing (1997) 11(2), 187 205 ANALYSIS OF COMPUTED ORDER TRACKING K. R. FYFE AND E. D. S. MUNCK Department of Mechanical Engineering, University of Alberta, Edmonton, Alberta,
More informationAnalysis, Synthesis, and Perception of Musical Sounds
Analysis, Synthesis, and Perception of Musical Sounds The Sound of Music James W. Beauchamp Editor University of Illinois at Urbana, USA 4y Springer Contents Preface Acknowledgments vii xv 1. Analysis
More informationEnglish Project. Contents
English Project Contents Introduction 2. Many-Talker Prompt-File Distribution 3. Few-Talker Prompt-File Distribution 4. Very-Few-Talker Prompt Files Introduction This report documents the subjects, equipment,
More informationHello, welcome to the course on Digital Image Processing.
Digital Image Processing Prof. P. K. Biswas Department of Electronics and Electrical Communications Engineering Indian Institute of Technology, Kharagpur Module 01 Lecture Number 05 Signal Reconstruction
More informationFast Ethernet Consortium Clause 25 PMD-EEE Conformance Test Suite v1.1 Report
Fast Ethernet Consortium Clause 25 PMD-EEE Conformance Test Suite v1.1 Report UNH-IOL 121 Technology Drive, Suite 2 Durham, NH 03824 +1-603-862-0090 Consortium Manager: Peter Scruton pjs@iol.unh.edu +1-603-862-4534
More informationRounding Considerations SDTV-HDTV YCbCr Transforms 4:4:4 to 4:2:2 YCbCr Conversion
Digital it Video Processing 김태용 Contents Rounding Considerations SDTV-HDTV YCbCr Transforms 4:4:4 to 4:2:2 YCbCr Conversion Display Enhancement Video Mixing and Graphics Overlay Luma and Chroma Keying
More informationPHGN 480 Laser Physics Lab 4: HeNe resonator mode properties 1. Observation of higher-order modes:
PHGN 480 Laser Physics Lab 4: HeNe resonator mode properties Due Thursday, 2 Nov 2017 For this lab, you will explore the properties of the working HeNe laser. 1. Observation of higher-order modes: Realign
More informationQuiz #4 Thursday, April 25, 2002, 5:30-6:45 PM
Last (family) name: First (given) name: Student I.D. #: Circle section: Hu Saluja Department of Electrical and Computer Engineering University of Wisconsin - Madison ECE/CS 352 Digital System Fundamentals
More informationNOTICE: This document is for use only at UNSW. No copies can be made of this document without the permission of the authors.
Brüel & Kjær Pulse Primer University of New South Wales School of Mechanical and Manufacturing Engineering September 2005 Prepared by Michael Skeen and Geoff Lucas NOTICE: This document is for use only
More informationMeasurement of overtone frequencies of a toy piano and perception of its pitch
Measurement of overtone frequencies of a toy piano and perception of its pitch PACS: 43.75.Mn ABSTRACT Akira Nishimura Department of Media and Cultural Studies, Tokyo University of Information Sciences,
More informationDISTRIBUTED MOTION CONTROL
DISTRIBUTED OTION CONTROL Jacob Hefer Elmo otion Control Westford, A Abstract Distributed motion control is a reality; today's processing power, deterministic protocols and network technology make that
More informationTP1a mask, noise and jitter for SRn
TP1a mask, noise and jitter for SRn Piers Dawe Avago Technologies IEEE P802.3ba Quebec May 2009 TP1a mask, noise and jitter for SRn 1 Supporters Mike Dudek Jonathan King Brian Misek John Petrilla Independent*
More informationDepartment of Electrical & Electronic Engineering Imperial College of Science, Technology and Medicine. Project: Real-Time Speech Enhancement
Department of Electrical & Electronic Engineering Imperial College of Science, Technology and Medicine Project: Real-Time Speech Enhancement Introduction Telephones are increasingly being used in noisy
More informationPage I-ix / Lab Notebooks, Lab Reports, Graphs, Parts Per Thousand Information on Lab Notebooks, Lab Reports and Graphs
Page I-ix / Lab Notebooks, Lab Reports, Graphs, Parts Per Thousand Information on Lab Notebooks, Lab Reports and Graphs Lab Notebook: Each student is required to purchase a composition notebook (similar
More informationFFT Laboratory Experiments for the HP Series Oscilloscopes and HP 54657A/54658A Measurement Storage Modules
FFT Laboratory Experiments for the HP 54600 Series Oscilloscopes and HP 54657A/54658A Measurement Storage Modules By: Michael W. Thompson, PhD. EE Dept. of Electrical Engineering Colorado State University
More informationHomework 3 posted this week, due after Spring break Quiz #2 today Midterm project report due on Wednesday No office hour today
EE241 - Spring 2013 Advanced Digital Integrated Circuits Lecture 14: Statistical timing Latches Announcements Homework 3 posted this week, due after Spring break Quiz #2 today Midterm project report due
More informationDPD80 Visible Datasheet
Data Sheet v1.3 Datasheet Resolved Inc. www.resolvedinstruments.com info@resolvedinstruments.com 217 Resolved Inc. All rights reserved. General Description The DPD8 is a low noise digital photodetector
More informationMIE 402: WORKSHOP ON DATA ACQUISITION AND SIGNAL PROCESSING Spring 2003
MIE 402: WORKSHOP ON DATA ACQUISITION AND SIGNAL PROCESSING Spring 2003 OBJECTIVE To become familiar with state-of-the-art digital data acquisition hardware and software. To explore common data acquisition
More information10:15-11 am Digital signal processing
1 10:15-11 am Digital signal processing Data Conversion & Sampling Sampled Data Systems Data Converters Analog to Digital converters (A/D ) Digital to Analog converters (D/A) with Zero Order Hold Signal
More informationGeneral Index. bandpass 140 bandpass filter 141 bandstop 140 bandstop filters 141 bar plot 26 bars 26 bathymetry 154
General Index A accessible population 2 adaptive filtering 143 adaptive process 143 addition 18 Aitchisons log-ratio transformation 223 alternative hypothesis 51 amplitude 134 analog filters 119 analysis
More informationExperiment: Real Forces acting on a Falling Body
Phy 201: Fundamentals of Physics I Lab 1 Experiment: Real Forces acting on a Falling Body Objectives: o Observe and record the motion of a falling body o Use video analysis to analyze the motion of a falling
More informationMusic Database Retrieval Based on Spectral Similarity
Music Database Retrieval Based on Spectral Similarity Cheng Yang Department of Computer Science Stanford University yangc@cs.stanford.edu Abstract We present an efficient algorithm to retrieve similar
More informationMATLAB Programming. Visualization
Programming Copyright Software Carpentry 2011 This work is licensed under the Creative Commons Attribution License See http://software-carpentry.org/license.html for more information. Good science requires
More informationOpen loop tracking of radio occultation signals in the lower troposphere
Open loop tracking of radio occultation signals in the lower troposphere S. Sokolovskiy University Corporation for Atmospheric Research Boulder, CO Refractivity profiles used for simulations (1-3) high
More informationPredicting the immediate future with Recurrent Neural Networks: Pre-training and Applications
Predicting the immediate future with Recurrent Neural Networks: Pre-training and Applications Introduction Brandon Richardson December 16, 2011 Research preformed from the last 5 years has shown that the
More information