Time series analysis
|
|
- Felicia Stephens
- 6 years ago
- Views:
Transcription
1 Time series analysis (July 12-13, 2011) Course Exercise Booklet MATLAB function reference 1
2 Introduction to time series analysis Exercise 1.1 Controlling frequency, amplitude and phase... 3 Exercise 1.2 Signal representation in the Frequency Domain... 4 Analysis in the Time Domain Exercise 2.1 Event Spacing of the Cave Creek runoff data... 5 Exercise 2.2 Autocorrelation of the Cave Creek runoff data... 6 Exercise 2.3 Autocorrelation of the SPECMAP Stack... 7 Signal Filtering Exercise 3.1 Removing long-term trends from the Mauna data set... 8 Exercise 3.2 Calculating the first difference curve of the Mauna data... 9 Exercise 3.3 Using a moving average Exercise 3.4 Filtering in the frequency domain The Frequency Domain and Spectral Analysis Exercise 4.1 Calculation of a periodogram Exercise 4.2 Time in the Frequency domain Exercise 4.3a Trends in the Frequency domain Exercise 4.3b Trends in the Frequency domain(2) Exercise 4.4 Processing unequally spaced data Exercise 4.5 The CLEAN algorithm Exercise 4.6 White noise in the Frequency Domain Exercise 4.7 Welch-Overlapped-Segment-Averaging Exercise 4.8 Red noise in the Time and Frequency Domains The Time-Frequency Plane Exercise 5.1 A Stationary Signal Exercise 5.2 A nonstationary signal Exercise 5.3 Evolutionary spectral analysis of a nonstationary signal Exercise 5.4: Evolutionary spectral analysis of the ODP O record PCA / EOF Analysis Exercise 7.1 PCA of a collection of time series 25 Exercise 7.2 PCA and the Hockey Stick".26 Correlation of two time series Exercise 8.1 Correlation and autocorrelation 27 Exercise 8.2 A semi-emperical ice volume model 28 2
3 Exercise 1.1 Controlling frequency, amplitude and phase. The first task of this exercise is to calculate and plot a sinusoid with a given frequency. To do this we need to make time points, define the frequency of the sinusoid and then calculate the cycle itself. >> time=[0:1:800] ; Time points between 0 and 800 with a spacing of 1 >> f=0.01 Define the frequency of the sinusoid as 0.01 >> signal=sin(2.*pi.*time.*f); Calculate the sinusoid >> plot(time,signal) Plot the sinusoid against time Now we can repeat the calculation but also define the amplitude of the sinusoid. >> time=[0:1:800] ; Time points between 0 and 800 with a spacing of 1 >> f=0.01 Define the frequency of the sinusoid as 0.01 >> A=2.5 Define the amplitude of the sinusoid as 2.5 >> signal=sin(2.*pi.*time.*f).*a; Calculate the sinusoid >> plot(time,signal) Plot the sinusoid against time Finally, calculate a sinusoid with a defined phase. >> time=[0:1:800] ; Time points between 0 and 800 with a spacing of 1 >> f=0.01 Define the frequency of the sinusoid as 0.01 >> phi=pi./2 Define the phase of the sinusoid as /2 >> signal=sin(2.*pi.*time.*f+phi); Calculate the sinusoid >> plot(time,signal) Plot the sinusoid against time It is possible to calculate three separate sinusoids with the same frequencies as the main Milankovitch cycles. These sinusoids can then be added together and plotted. >> time=[0:1:800] ; Time points between 0 and 800 with a spacing of 1 >> fe=1./100; Define the eccentricity frequency >> eccen=sin(2.*pi.*time.*fe); Calculate a sinusoid with frequency fe >> fo=1./41; Define the obliquity frequency >> obliq=sin(2.*pi.*time.*fo); Calculate a sinusoid with frequency fo >> fp=1./21; Define the precession frequency >> prec=sin(2.*pi.*time.*fp); Calculate a sinusoid with frequency fp >> final=eccen+obliq+prec; Sum the 3 sinusoids together to make a final signal >> plot(time,final); Plot the signal against time Does your signal which contains the Milankovitch periods look like known patterns of orbital scale climate change? 3
4 Exercise 1.2 Signal representation in the Frequency Domain The aim of this exercise is for you to investigate how a signal made from a mixture of sine waves is represented in the Frequency Domain (in other words what happens when you perform the Fourier transform. Using MATLAB calculate 3 different sine waves (i.e. with different amplitudes, frequencies and phases) and add them together to give a composite signal. 1. Produce a plot which shows your composite signal 2. Use fft_plot to produce a plot that shows the frequency spectrum of your composite signal. So, the next part is the important bit; 3. Provide an interpretation of the frequency spectrum that you obtained from your composite signal. Think about the positions of the peaks in the diagram and their relative heights. 4
5 Exercise 2.1 Event Spacing of the Cave Creek runoff data The time series we will study shows the monthly amount of runoff water (measured in inches) from Cave Creek in Kentucky. Use the event spacing method to estimate the period of each runoff cycle. In MATLAB we first have to load the data into the memory. You can load the Cave Creek data using: >> load cave_creek >> whos This will tell you which variables are in the memory At the bottom of the screen you will see some text which looks like: Name Size Bytes Class month 216x double array runoff 216x double array Grand total is 432 elements using 3456 bytes This shows us that the variables month and runoff have been loaded into the memory. We can now make a plot of the data: >> plot(month,runoff) >> xlabel('month') Add a label to the x-axis >> ylabel('runoff (inches)') Add a label to the y-axis We will use the function ginput.m to mark the points on the chart which we think are events: >> [x,y]=ginput Obtain data from the figure You can now use the mouse to record the positions that you think correspond to events in the runoff data. When you have clicked on all the points of interest, press Enter to stop the function and return to the MATLAB prompt. The variables x and y now correspond to the coordinates of the positions you clicked on. You can mark these points on your chart in the following way: >> hold on This allows you to add more data to the existing plot. >> plot(x,y, o ) Add the click positions to the plot as circles. We just need to find the difference between the values in x, in MATLAB this is simple: >> x_diff=diff(x) >> x_mean=mean(x_diff); mean is the command for calculating a mean >> x_std=std(x_diff); std is the command for calculating a standard deviation The event spacing of the Cave Creek runoff data = ± months What is your interpretation of the mean event spacing that you calculated? 5
6 Exercise 2.2 Autocorrelation of the Cave Creek runoff data The time series we will study shows the monthly amount of runoff water (measured in inches) from Cave Creek in Kentucky. Use the event spacing method to estimate the period of each runoff cycle. In MATLAB we first have to remove any data which might still be in the memory using the clear command and then load the Cave Creek data; >> clear all, close all Remove all data from the memory >> load cave_creek Load the data set Our first task is to detrend the data. To do this we can use the function detrend_signal.m >> runoff_d=detrend_signal(month,runoff,1); This will produce the detrended data runoff_d, which has had a straight-line trend removed. A figure is produced by the function that shows the original data, the fitted line and the detrended data. The data is now prepared for the autocorrelation. To do the analysis we use the function autocorr.m >> [rt,lags]=autocorr(runoff_d); This function has two outputs rt and lags >> figure Make a new figure window >> plot(lags,rt) Plot the final autocorrelogram >> xlabel('lag number'); Add a label to the x-axis >> ylabel('r'); Add a label to the y-axis Again, here is the important part, what is your interpretation of the autocorrelogram and how does it compare to the Event Spacing method? 6
7 Exercise 2.3 Autocorrelation of the SPECMAP Stack Load the SPECMAP file into MATLAB, you will find it contains two variables; age and data (the units of age are ka and the data units are normalised oxygen isotope ). Use MATLAB to produce a plot of the SPECMAP record. Add appropriate labels to both the x and the y axes. Detrend the SPECMAP signal. Perform an autocorrelation including significance levels on the detrended series. Make an interpretation of the autocorrelation in terms of known climate variation. 7
8 Exercise 3.1 Removing long-term trends from the Mauna data set We can use the function detrend_signal to remove long-term change from the Mauna data and keep only the higher frequency variation. >> load mauna load the data into MATLAB There are two variables time (units of days) and mauna (CO 2 data, units of ppm) First, try removing a simple linear trend (a straight line) of the form: y a1x a0 >> xc=detrend_signal(time,mauna,1) In this case the 1 tells MATLAB to remove a 1 st order polynomial from the data, i.e. a straight line. 2 Next, try removing a quadratic function of the form: y a2 x a1x a0 >> xc=detrend_signal(time,mauna,2) In this case the 2 tells MATLAB to remove a 2 nd order polynomial from the data, i.e. a quadratic function. 3 2 Finally, try removing a cubic function of the form: y a3x a2 x a1x a0 >> xc=detrend_signal(time,mauna,3) In this case the 3 tells MATLAB to remove a 3 rd order polynomial from the data, i.e. a cubic function. 8
9 Exercise 3.2 Calculating the first difference curve of the Mauna data Calculate the first difference curve for the Mauna CO 2 data. Then plot both the original data and the first-difference curve. >> load mauna load the data into MATLAB There are two variables time (units of days) and mauna (CO 2 data, units of ppm) >> xd=diff(mauna); Calculate the first-difference of the CO 2 data >> xt=(diff(time)./2)+time(1:229); Calculate the mid-points of the time data Now we can start the plotting, we want to plot two sets of axes in one window, so we can use the subplot command. >> subplot(2,1,1) This tells MATLAB to set the figure window into a 2 x 1 matrix of plotting areas, and chooses the 1 st area to be active. >> plot(time,mauna) Plot the original data >> xlabel('time (yrs)') Puts a title onto the x-axis >> ylabel('co2 content (ppm)') Puts a title onto the y-axis Now we want to make the second plot so we give the subplot command again, but now choose the 2 nd area to be active. >> subplot(2,1,2) >> plot(xt,xd) Plot the first-difference data >> xlabel('time (yrs)') Puts a title onto the x-axis >> ylabel('first Difference (ppm)') Puts a title onto the y-axis 9
10 Exercise 3.3 Using a moving average The MATLAB file porosity contains normalised porosity data from core GeoB taken from the equatorial Atlantic. Use the maverage function to produce a smooth record of the data. What is a good smoothing window to use? >> load porosity load the porosity data First look at the data by making a simple plot >> plot(age,porosity,'g') This plots the data as a green line >> xlabel('age (kyr)') Add a x-axis title >> ylabel('normalised porosity') Add a y-axis title >> hold on Allow lines to be added to the plot You can read the instructions on the use of maverage. One required input is a matrix that describes the moving average window to be used. For example; >> f=[ ] 5pt moving average with equal weights >> f=[ ] 5pt weighted moving average >> [Xout,Yout]=maverage(age,porosity,f); Calculate the smoothed data >> plot(xout,yout,'r') Plot the smoothed data as a red line If you want to try different smoothing windows in a new plot give the command; >> figure this will produce a new plot window in which you can plot the new data. 10
11 Exercise 3.4 Filtering in the frequency domain Use the function filter_signal to filter the normalised porosity from core GeoB4311. Use the subplot command to make plots which compare the filtered porosity record to the SPECMAP stack. The normalised oxygen isotope data for the SPECMAP stack can be found in the MATLAB file SPECMAP. 11
12 Exercise 4.1 Calculation of a periodogram In this exercise we will construct two signals with different frequencies and amplitudes and determine their variance. The signals are then combined together and the periodogram (unsmoothed variance spectrum) is calculated using the function pdg. We can then check that the variance returned in the spectrum is the same as in the individual input components. >> time=[0:1:1000]'; Form a time array between 0 and 1000 >> f1=1./100; Frequency of the first signal >> f2=1./40 Frequency of the second signal >> data1=sin(2.*pi.*time.*f1); Calculate the first sinusoid >> data2=sin(2.*pi.*time.*f2); Calculate the second sinusoid Now we have the sinusoids but they have equal amplitudes and equal variances. If we change the amplitude then that will also change the variance. >> data1=data1.*4.5; Set the amplitude of the first signal to 4.5 >> data2=data2.*7.0; Set the amplitude of the second signal to 7.0 >> v1=var(data1) Calculate the variance of the first signal >> v2=var(data2) Calculate the variance of the second signal >> input_data=data1+data2; Combine the signals to give the final time series >> subplot(2,1,1) Create a plot in the upper half of the figure >> plot(time,data1,time,data2) Plot the individual sinusoids >> xlabel('time') Add a x-axis title >> ylabel('signal') Add a y-axis title >> subplot(2,1,2) Create a plot in the lower half of the figure >> plot(time,input_data) Plot the final input time series >> xlabel('time') Add a x-axis title >> ylabel('signal') Add a y-axis title The final part of this exercise is to calculate the periodogram and plot the results. >> [f,pyy]=pdg(time,input_data); Calculate the periodogram >> figure Make a new figure window >> plot(f,pyy) Plot frequency against the variance >> xlim([0 0.05]) Set the maximum plotted frequency to 0.05 >> xlabel('frequency') Add a x-axis title >> ylabel('variance') Add a y-axis title You should now be able to compare the peak heights in the periodogram (i.e. the variances of the different spectral components) to the variances you calculated for the individual input sinusoids, are they the same? 12
13 Exercise 4.2 Time in the Frequency domain We can compare the frequency spectra of the original SPECMAP stack and a modified version. These signals are basically the opposite of each other, but how are their frequency spectra different? >> clear all Remove all existing data from the memory >> close all Close all the existing figure windows >> load SPECMAP2 Load the data into the memory There should now be 3 variables in MATLAB s memory. The first is age, then data1 which is the normal SPECMAP record and data2 which is the flipped and reversed stack. First you can plot the data to check they are okay. >> subplot(2,1,1) Create a plot in the upper half of the figure >> plot(age,data1); Plot the original SPECMAP data >> subplot(2,1,2) Create a plot in the lower half of the figure >> plot(age,data2); Plot the reversed SPECMAP data Now a periodogram can be calculated for each signal and they can be compared in the frequency domain. >> [f1,pyy1]=pdg(time,data1) Periodogram of the original SPECMAP data >> figure Create a new plot window >> subplot(2,1,1) Create a plot in the upper half of the figure >> plot(f1,pyy1) Plot the periodogram >> xlabel('frequency') Add a x-axis title >> ylabel('variance (original)') Add a y-axis title >> [f2,pyy2]=pdg(time,data2) Periodogram of the reversed SPECMAP data >> subplot(2,1,2) Create a plot in the lower half of the figure >> plot(f2,pyy2) Plot the periodogram >> xlabel('frequency') Add a x-axis title >> ylabel('variance (flipped)') Add a y-axis title You should now be able to compare the periodograms of the two signals. What differences can you see? 13
14 Exercise 4.3a Trends in the Frequency domain Using the raw Mauna Loa CO 2 data we can investigate the effect long-term trends have in the frequency domain >> clear all Remove all existing data from the memory >> close all Close all the existing figure windows >> load mauna Load the Mauna Loa data into the memory There are two variables in the data set, time (measured in days) and mauna which is the CO 2 data (measured in ppm). >> subplot(2,1,1) Create a plot in the upper half of the figure >> plot(time,mauna); Plot the original CO 2 data >> xlabel('time (days)') >> ylabel('co2 cotent (ppm)') >> subplot(2,1,2) Create a plot in the lower half of the figure >> [f,pyy]=pdg(time,mauna) Periodogram of the CO 2 data >> plot(f,pyy); Plot the periodogram >> xlabel('frequency (1/days)') >> ylabel('variance') What structure do you see in the periodogram? 14
15 Exercise 4.3b Trends in the Frequency domain(2) You can now try detrending the Mauna Loa CO 2 data using the function detrend_signal. What effect does the detrending have on the periodogram. >> xc=detrend_signal(time,mauna) Detrend the data with a 1 st order polynomial We can no repeat the previous part of the exercise, but using the detrended rather than the raw data. >> figure Make a new figure window >> subplot(2,1,1) Create a plot in the upper half of the figure >> plot(time,xc); Plot the detrended CO 2 data >> xlabel('time (days)') >> ylabel('detrended signal') >> subplot(2,1,2) Create a plot in the lower half of the figure >> [f,pyy]=pdg(time,xc) Periodogram of the detrended data >> plot(f,pyy); Plot the periodogram >> xlabel('frequency (1/days)') >> ylabel('variance') How has the structure of the periodogram changed now the data has been detrended? 15
16 Exercise 4.4 Processing unequally spaced data If we have data that are unequally spaced in the time domain then we must interpolate them onto an equally spaced time axis before calculating a traditional periodogram. In the file interpolation there are two versions of the same signal, one equally space and one unequally spaced in time. >> load interpolation The variables time_reg and data_reg make up the equally spaced data set, whilst time_irreg and data_irreg are the unequally spaced data set. First we can plot the signal. >> subplot(3,1,1) Activate the top third of the plot window >> plot(time_reg,data_reg); Plot the regularly spaced signal >> xlabel('time') >> ylabel('signal') You will see this is a complicated signal and to understand it we should study the periodogram of the equally space data. >> [f,pyy]=pdg(time_reg,data_reg); Periodogram of the regularly spaced signal >> subplot(3,1,2) Activate the middle third of the plot window >> plot(f,pyy) Plot the Periodogram >> xlabel('frequency') >> ylabel('variance') Examination of the periodogram shows that the signal is made up of 8 sinusoids with equal amplitudes and frequencies at spacings of To calculate a periodogram for the irregularly spaced signal we must interpolate the data onto a regularly spaced time array. The simplest way to do this is to interpolate it onto the original time_reg points. >> data_interp=interp1(time_irreg,data_irreg,time_reg,'linear') The above command performs a one-dimensional linear interpolation of the irregular signal (time_irreg and data_irreg) onto the regular time array (time_reg). We can now use the original time_reg points with the data_interp array to calculate a periodogram for the interpolated signal. >> [f,pyy]=pdg(time_reg,data_interp); Periodogram of the interpolated signal >> subplot(3,1,3) Activate the bottom third of the plot window >> plot(f,pyy) Plot the Periodogram >> xlabel('frequency') >> ylabel('variance') What effect does the interpolation appear to have on the signal? 16
17 Exercise 4.5 The CLEAN algorithm The CLEAN algorithm is just one method where frequency spectra can be calculated directly from unequally spaced data without a need for interpolation. To apply CLEAN to the data simply give the command: >> [f,pyy]=clean(time_irreg,data_irreg) You can see from the input variables that we are putting the irregularly spaced data directly into the algorithm without any interpolation. >> figure >> plot(f,pyy) Plot the CLEAN spectrum >> xlabel('frequency') >> ylabel('variance') The high frequency parts of the signal are not attenuated and CLEAN has successfully produced an accurate spectrum for the signal. There are a number of alternative methods available for spectral analysis of unevenly spaced data. A few references are: Roberts D.H., Lehar J., Dreher J.W., (1987) Time Series Analysis with Clean - I - Derivation of a Spectrum. Astronomical J. 93, (4) 968 J.D. Scargle, (1982) Studies in astronomical time series analysis II. Statistical aspects of spectral analysis of unevenly sampled data. Astrophysical J. 263, And with specific reference to palaeoclimatic data: Schulz, M. and Stattegger, K. (1997): SPECTRUM: Spectral analysis of unevenly spaced paleoclimatic time series. Comput. Geosci., 23, Heslop D. and Dekkers M.J. (2002). Spectral analysis of unevenly spaced climatic time series using CLEAN: signal recovery and derivation of significance levels using a Monte Carlo simulation, Phys. Earth Planet. Inter., 130,
18 Exercise 4.6 White noise in the Frequency Domain We will construct 3 white noise signals with different lengths (64, 256 and 1024 data points) and study the form of their periodograms. We can generate normally distributed random numbers (with a mean = 0 and variance = 1.0) in MATLAB using the randn function. >> time1=[0:1:63]; Generate a time array ranging from 0 to 63 >> signal1=randn(1,64); 64 normally distributed random numbers for the signal >> time2=[0:1:255]; Generate a time array ranging from 0 to 255 >> signal2=randn(1,256); 256 normally distributed random numbers for the signal >> time3=[0:1:1023]; Generate a time array ranging from 0 to 1023 >> signal3=randn(1,1024);1024 normally distributed random numbers for the signal It is a good idea to first plot the signals so you have some idea what they look like. >> subplot(3,1,1) Activate the top one-third of the plot window >> plot(time1,signal1) Plot signal1 >> subplot(3,1,2) Activate the middle one-third of the plot window >> plot(time2,signal2) Plot signal2 >> subplot(3,1,3) Activate the bottom one-third of the plot window >> plot(time3,signal3) Plot signal3 We can now calculate the periodogram for each signal and plot them in a new figure window. >> figure Make a new plot window >> [f,pyy]=pdg(time1,signal1) Calculate the periodogram of signal1 >> subplot(3,1,1) Activate the top one-third of the plot window >> plot(f,pyy) Plot the periodogram of signal1 >> xlabel('frequency') >> ylabel('variance') >> [f,pyy]=pdg(time2,signal2) Calculate the periodogram of signal2 >> subplot(3,1,2) Activate the middle one-third of the plot window >> plot(f,pyy) Plot the periodogram of signal2 >> xlabel('frequency') >> ylabel('variance') >> [f,pyy]=pdg(time3,signal3) Calculate the periodogram of signal3 >> subplot(3,1,3) Activate the bottom one-third of the plot window >> plot(f,pyy) Plot the periodogram of signal3 >> xlabel('frequency') >> ylabel('variance') 18
19 Exercise 4.7 Welch-Overlapped-Segment-Averaging Using the pwelch function it is possible to investigate how the power spectra of a white noise signal changes as the WOSA method is applied >> clear all Remove existing variables from the memory >> close all Close all the existing figure windows >> signal=randn(1,1024); 1024 normally distributed random numbers for the signal We can now calculate the spectra for the signal with different segment lengths (1024, 256 and 64) and then plot the results. >> [Pyy,f]=pwelch(signal,1024,[],[],1) Calculate with a 1024 point segment >> subplot(3,1,1) Activate the top one-third of the plot window >> plot(f,pyy) Plot the spectrum >> xlabel('frequency') >> ylabel('power') >> [Pyy,f]=pwelch(signal,256,[],[],1) Calculate with 256 point segments >> subplot(3,1,2) Activate the middle one-third of the plot window >> plot(f,pyy) Plot the spectrum >> xlabel('frequency') >> ylabel('power') >> [Pyy,f]=pwelch(signal,64,[],[],1) Calculate with 64 point segments >> subplot(3,1,3) Activate the bottom one-third of the plot window >> plot(f,pyy) Plot the spectrum >> xlabel('frequency') >> ylabel('power') Compare the frequency spectra. It is important to consider how close they are to the theoretical spectra for white noise and their relative frequency resolutions. 19
20 Exercise 4.8 Red noise in the Time and Frequency Domains The function AR1n can be used to generate red noise series with a given value of. We will generate 3 different series and see what they look like in both the time and frequency domains. >> clear all Remove existing variables from the memory >> close all Close all the existing figure windows >> rho=0.0 Define the lag-one autocorrelation coefficient >> [Rn1,Wn1]=AR1n(rho,2048); Create a series of length 2048 The AR1 function has two outputs; Rn1 is the red noise series and Wn1 is the white noise series that was used in the construction of Rn1. >> figure Generate a new plot window >> subplot(2,1,1) Activate the upper half of the plot window >> plot([0:1:2047],wn1, b ) Plot the white noise as a blue line >> hold on Allow more data to be added to the chart >> plot([0:1:2047],rn1, r ) Plot the red noise as a red line We saw in the previous exercise that to obtain a consistent spectrum we should use the WOSA method included in the function pwelch. >> [Pyy,f]=pwelch(Rn1,64,[],[],1) Calculate with a 64 point segment >> subplot(2,1,2) Activate the lower half of the plot window >> hold off Switch off the hold command >> plot(f,pyy) Plot the spectrum >> xlabel('frequency') >> ylabel('power') Because we set = 0 the calculated series is simply white noise and should have a flat spectrum. What happens if you have a non-zero value for the lag-one autocorrelation coefficient? Repeat the exercise with = 0.8 and =
21 Exercise 5.1 A Stationary Signal First construct a vector containing the points in time, in this case starting at 0, and finishing at 3199 in intervals of 1 unit. >> time=[0:1:3199]'; We will calculate 4 sine waves with different frequencies, so first we define the frequency values: >> f1=1./400; Long period Eccentricity >> f2=1./100; Eccentricity >> f3=1./41; Obliquity >> f4=1./23; Precession Now calculate a sine wave for each of the 4 frequencies: >> signal1=sin(2.*pi.*f1.*time); >> signal2=sin(2.*pi.*f2.*time); >> signal3=sin(2.*pi.*f3.*time); >> signal4=sin(2.*pi.*f4.*time); Finally we add all 4 signals together to produce a composite signal: >> signal=signal1+signal2+signal3+signal4; >> plot(time,signal) Now determine the periodogram using pdg and plot the resulting frequency spectrum and set the x-axis limit between 0 and 0.1 >> [f,pyy]=pdg(time,signal); The frequency data is given in the variable f and the power of the signal at each frequency is stored in the variable Pyy. >> figure >> plot(f,pyy) >> xlim([0 0.1]) >> xlabel('frequency') >> ylabel('variance') 21
22 Exercise 5.2 A nonstationary signal Now we ll construct a signal that changes through time. We will work with 4 different time segments and each one will contain a cycle with a different frequency. Finally, we will combine the segments and calculate a periodogram. Segment 1: spans 0 to 799 kyr with a frequency of 1/400 kyr -1 Segment 2: spans 800 to 1599 kyr with a frequency of 1/100 kyr -1 Segment 3: spans 1600 to 2399 kyr with a frequency of 1/41 kyr -1 Segment 4: spans 2400 to 3199 kyr with a frequency of 1/23 kyr -1 As a first step we construct our 4 different time segments: >> clear all, close all >> time1=[0:1:799]'; Time segment spanning 0 to 799 >> time2=[800:1:1599]'; Time segment spanning 800 to 1599 >> time3=[1600:1:2399]'; Time segment spanning 1600 to 2399 >> time4=[2400:1:3199]'; Time segment spanning 2400 to 3199 Then define the 4 different frequencies: >> f1=1./400; Long period Eccentricity >> f2=1./100; Eccentricity >> f3=1./41; Obliquity >> f4=1./23; Precession Now calculate the sine waves for each time segment: >> signal1=sin(2.*pi.*f1.*time1); Long eccentricity spanning time 0 to 799 >> signal2=sin(2.*pi.*f2.*time2); Eccentricity spanning time 800 to 1599 >> signal3=sin(2.*pi.*f3.*time3); Obliquity spanning time 1600 to 2399 >> signal4=sin(2.*pi.*f4.*time4); Precession spanning time 2400 to 3199 Now we must combine the time segments and signal segments together. We do this using the square brackets operator, which allows us to glue existing variables together: >> time=[time1;time2;time3;time4]; This tells MATLAB to make a new variable called time that starts with the vector time1. Underneath time1 it places vector time2, then underneath that places time3 and finally places time4 at the bottom. Now we do the same procedure for the signal, plot the result and perform the spectral analysis: >> signal=[signal1;signal2;signal3;signal4]; >> plot(time,signal) >> [f,pyy]=pdg(time,signal); >> plot(f,pyy) >> xlim([0 0.1]) 22
23 Exercise 5.3 Evolutionary spectral analysis of a nonstationary signal Perform evolutionary spectral analysis on the nonstationary signal you produced in Exercise 5.2. Perform the analysis with a spacing of 10 kyr and try some different window lengths (for example 200 kyr, 500 kyr, 800 kyr, 1200 kyr). e.g. >> evolpsd(time,signal,200,10); >> evolpsd(time,signal,500,10); >> evolpsd(time,signal,800,10); >> evolpsd(time,signal,1200,10); What influence does varying the window length have? What size window lengths do you need to obtain a good resolution in time and frequency? 23
24 Exercise 5.4: Evolutionary spectral analysis of the ODP O record Load the data into MATLAB from the file ODP677 >> load odp677 There are two variables: age and data The ODP O record is nonstationary, make an interpretation of its spectral content and climatic information using evolutive spectral analysis. Do the relative contributions of the different Milankovitch cycles change through time? >> evolpsd(age,data,window,spacing) adjust the window size to get the best balance between resolution time and frequency. Extra Hint: Before you start the analysis plot the data, do you think it is necessary to detrend the signal first? 24
25 Exercise 7.1 PCA of a collection of time series The data are stored in the file pca_climate.mat. There are 4 variables, age (a common time scale for the data), ms (magnetic susceptibility), gs (grain size) and d18o (oxygen isotopes). We ll now plot the 3 time series. >> clear all, close all Clear the memory and close all figures >> load pca_climate Load the data file >> figure Create a new figure >> subplot(3,1,1) axes in 3 rows & 1 column, activate first set of axes >> plot(age,zscore(ms)) plot the standardized magnetic susceptibility data >> ylabel('normalized Susceptibility') label the y-axis >> subplot(3,1,2) axes in 3 rows & 1 column, activate second set of axes >> plot(age,zscore(gs)) plot the standardized grain size data >> ylabel('normalized Grain size') label the y-axis >> subplot(3,1,3) axes in 3 rows & 1 column, activate second set of axes >> plot(age,zscore(d18o)) plot the standardized oxygen isotope data >> ylabel('normalized d18o') label the y-axis >> xlabel('age [ka]') label the x-axis We ll combine 3 time series which are on a common timescale into a single matrix and then standardize the columns. If the time series were not on a common timescale then they would have to be interpolated first to obtain versions that are on a common timescale. PCA is performed using the function princomp. With the results we can find out what proportion of the variance the 1st PC explains and plot the scores as a function of age. >> X=[ms,gs,d18O]; form a matrix composed of the time series >> X=zscore(X); standardize the columns of X >> [loadings,scores,latent]=princomp(x); perform the PCA >> latent=latent./sum(latent) normalize the PC contributions >> latent(1) the variance explained by the 1st PC >> figure Generate a new figure >> plot(age,scores(:,1)) plot the scores of the 1st PC >> ylabel('pc Score') label the y-axis >> xlabel('age [ka]') label the x-axis Using the example above plot the scores of the 2nd PC as a function of age. What proportion of the total variance does the 2nd PC account for? 25
26 Exercise 7.2 PCA and the Hockey Stick We ll run a test with a collection of red noise time series to analyze the PCA technique employed by Mann et al. (1999) and see how it compares to the classical PCA approach. To start we ll create a collection of age points between 1400 and 2000 AD and then generate 112 red noise time series with =0.8. >> clear all, close all >> age=[1400:1:2000]'; Create the age points >> Rn1=AR1n(0.8,numel(age),112); Generate 112 red time series >> figure New figure >> plot(age,rn1) Plot the time series >> xlabel('age [Yr]') Label the x axis >> ylabel('red Noise input') Label the y axis Now calculate the 1 st principal component of the data using the classical approach where the standard deviation of each record is set to 1 and the mean is set to 0. >> A=bsxfun(@minus,Rn1,mean(Rn1)); Subtract the mean of each record >> A=bsxfun(@rdivide,A,std(A)); Divide by the std of each record >> [coeffa,scorea] = princomp_hs(a); Perform the PCA on the data in A >> figure New figure >> plot(age,zscore(scorea(:,1)),'b') Plot the scores of the 1 st component >> xlabel('age [Yr]') Label the x axis >> ylabel('pca Score') Label the y axis Now we ll calculate the 1 st principal component of the data using the Mann approach where the standard deviation of each record is set to 1 during the period AD and the mean is set to 0 during the same period. >> idx=find(age>=1900 & age<=2000); Find points in the calibration period >> B=bsxfun(@minus,Rn1,mean(Rn1(idx,:))); Normalized the mean >> B=bsxfun(@rdivide,B,std(B(idx,:))); Normalized the standard deviation >> [coeffb,scoreb] = princomp_hs(b); Perform the PCA on the data in B >> hold on Add to the current plot >> plot(age,zscore(scoreb(:,1)),'r') Plot the scores of the 1 st component >> legend('classical PCA','Mann PCA',0) Add a legend Look at the final plot and compare the time series generated by the classical and Mann approaches. Remember, we are working with random data, that will show no preference to increasing or decreasing in the calibration period. 26
27 Exercise 8.1 Correlation and autocorrelation We ll generate two random red noise time series consisting of 200 points and test the significance of their correlation with and without employing the effective sample size. >> clear all, close all Clear the memory, close all existing figures >> N=200; Define the length of the time series >> X=AR1n(0.95,N); 200 point red noise series with =0.95 >> Y=AR1n(0.99,N); 200 point red noise series with =0.99 >> figure Generate a new figure >> subplot(2,1,1) 2 x 1 subplot, activate the first plot >> plot(1:200,x,'k') Plot the X time series >> ylabel('x Series') label the y-axis >> subplot(2,1,2) 2 x 1 subplot, activate the second plot >> plot(1:200,y,'k') Plot the X time series >> ylabel('y Series') label the y-axis >> xlabel('time') label the x-axis We can now form a bivariate plot and calculate the correlation >> figure, plot(x,y,'.') Form an XY plot of the data >> xlabel('x'), ylabel('y') Label the axis >> r = corrcoef(x,y) Calculate the correlation matrix >> r = abs(r(2,1)); Find the absolute value of r First we ll assess if the correlation is significant whilst ignoring the existence of any autocorrelation. >> t0 = r.*sqrt((n-2)./(1-r.^2)) Calculate the t statistic >> t_crit0=tinv( /2,n-2) Calculate the critical value of t Is t0 > t_crit0, if so we can state that the two time series are significantly correlated at the = 0.05 level. If not then the correlation is not significant. Now we ll assess if the correlation is significant whilst taking the existence of autocorrelation into account. >> Neff=N.*(1-0.95*0.99)./(1+0.95*0.99) Effective number of samples >> t1 = r.*sqrt((neff-2)./(1-r.^2)) Calculate the t statistic >> t_crit1=tinv( /2,neff-2) Calculate the critical value of t Is t0 > t_crit0, if so we can state that the two time series are significantly correlated at the = 0.05 level. If not then the correlation is not significant. 27
28 Exercise 8.2 A semi-empirical ice volume model We ll test the significance of the correlation upon which ice volume predictions are made. First we ll load the data and plot the time series. >> clear all, close all Clear the memory >> load sealevel_data load the data >> temp=detrend(temp) detrend the temperature rise >> rate=detrend(rate) detrend the rate >> subplot(2,1,1) 2 x 1 plots, use the 1st plot >> plot(temp,'.k-') plot the temperature rise >> ylabel('rate of Change (cm/yr)') label the y-axis >> subplot(2,1,2) 2 x 1, use the 2nd plot >> plot(rate,'.k-') plot the rate >> ylabel('warming above mean') label the y-axis >> xlabel('time') label the x-axis Now test the significance of the relationship whilst not taking the autocorrelation into account. >> r = corrcoef(temp,rate); Calculate the correlation matrix >> r = abs(r(2,1)); The absolute correlation coefficient >> N=length(rate) Number of values in the time series >> t0 = r.*sqrt((n-2)./(1-r.^2)) calculate the t statistic >> t_crit0=tinv( /2,n-2) calculate the critical value of t Now we ll repeat the test, but first estimating the effective number of samples based on the autocorrelation coefficients. >> r = corrcoef(temp,rate); Calculate the correlation matrix >> r = abs(r(2,1)); The absolute correlation coefficient >> N=length(rate) Number of values in the time series >> rho_temp=ar1(temp); AR1 coefficient of temp >> rho_rate=ar1(rate); AR1 coefficient of rate >> Neff=N.*(1-rho_temp*rho_rate)./(1+rho_temp*rho_rate) effective N >> t1 = r.*sqrt((neff-2)./(1-r.^2)) calculate the t statistic >> t_crit1=tinv( /2,neff-2) calculate the critical value of t Draw conclusions concerning the validity of the model, given the values in t0, t_crit0, t1 and t_crit1. 28
29 Using the maverage.m function [Xout,Yout]=maverage(x,y,f) This function applies a moving average to your data Inputs x: x-axis data, this will normally be depth or time y: y-axis data, this is the signal data f: the moving average window. Xout: the smoothed version of x Yout: the smoothed version of y. Output Examples For a traditional 3-point moving average use: [Xout,Yout]=maverage(x,y,[1 1 1]) For a traditional 5-point moving average use: [Xout,Yout]=maverage(x,y,[ ]) To apply different weights to each position just enter them into f For a 5-point weighted moving average use: [Xout,Yout]=maverage(x,y,[ ]) (this applies a weight of 1 to the first point, 3 to the second point, 5 to the third point, 3 to the forth point and 1 to the fifth point). You can use any window length and weightings you want. 29
30 Using the filter_signal.m function [Xout,Yout]=filter_signal(x,y,f,type) This function applies a frequency domain filter to your data. Inputs x: x-axis data, this will normally be depth or time y: y-axis data, this is the signal data f = cut-off frequency, for a low-pass or high-pass filter this is a single cut-off value. For a band-pass filter you must provide two frequencies which are the minimum and maximum frequencies that will be allowed through the filter, i.e. [fmin fmax]. type: what kind of filter should the function perform you can enter: 'low' for a low-pass filter 'high' for a high-pass filter 'band' for a band-pass filter Xout: the output of x. Yout: the filtered version of y. Output Examples [Xout,Yout]=filter_signal(x,y,f,'low') for a low-pass filter. [Xout,Yout]=filter_signal(x,y,f,'high') for a high-pass filter. [Xout,Yout]=filter_signal(x,y,[fmin fmax],'band') for a band-pass filter. 30
LabView 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 informationPS User Guide Series Seismic-Data Display
PS User Guide Series 2015 Seismic-Data Display Prepared By Choon B. Park, Ph.D. January 2015 Table of Contents Page 1. File 2 2. Data 2 2.1 Resample 3 3. Edit 4 3.1 Export Data 4 3.2 Cut/Append Records
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 informationNENS 230 Assignment #2 Data Import, Manipulation, and Basic Plotting
NENS 230 Assignment #2 Data Import, Manipulation, and Basic Plotting Compound Action Potential Due: Tuesday, October 6th, 2015 Goals Become comfortable reading data into Matlab from several common formats
More informationDigital 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 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 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 informationInvestigation of Digital Signal Processing of High-speed DACs Signals for Settling Time Testing
Universal Journal of Electrical and Electronic Engineering 4(2): 67-72, 2016 DOI: 10.13189/ujeee.2016.040204 http://www.hrpub.org Investigation of Digital Signal Processing of High-speed DACs Signals for
More informationLaboratory Assignment 3. Digital Music Synthesis: Beethoven s Fifth Symphony Using MATLAB
Laboratory Assignment 3 Digital Music Synthesis: Beethoven s Fifth Symphony Using MATLAB PURPOSE In this laboratory assignment, you will use MATLAB to synthesize the audio tones that make up a well-known
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 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 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 informationLab 5 Linear Predictive Coding
Lab 5 Linear Predictive Coding 1 of 1 Idea When plain speech audio is recorded and needs to be transmitted over a channel with limited bandwidth it is often necessary to either compress or encode the audio
More informationMDF Exporter This new function allows output in the ASAM MDF 4.0 file format, which is now widely used in the auto industry. Data import from the DS-3000 series Data Station Recorded data is transferred
More informationRestoration of Hyperspectral Push-Broom Scanner Data
Restoration of Hyperspectral Push-Broom Scanner Data Rasmus Larsen, Allan Aasbjerg Nielsen & Knut Conradsen Department of Mathematical Modelling, Technical University of Denmark ABSTRACT: Several effects
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 informationSwept-tuned spectrum analyzer. Gianfranco Miele, Ph.D
Swept-tuned spectrum analyzer Gianfranco Miele, Ph.D www.eng.docente.unicas.it/gianfranco_miele g.miele@unicas.it Video section Up until the mid-1970s, spectrum analyzers were purely analog. The displayed
More informationSpeech and Speaker Recognition for the Command of an Industrial Robot
Speech and Speaker Recognition for the Command of an Industrial Robot CLAUDIA MOISA*, HELGA SILAGHI*, ANDREI SILAGHI** *Dept. of Electric Drives and Automation University of Oradea University Street, nr.
More informationON THE INTERPOLATION OF ULTRASONIC GUIDED WAVE SIGNALS
ON THE INTERPOLATION OF ULTRASONIC GUIDED WAVE SIGNALS Jennifer E. Michaels 1, Ren-Jean Liou 2, Jason P. Zutty 1, and Thomas E. Michaels 1 1 School of Electrical & Computer Engineering, Georgia Institute
More informationAlgebra I Module 2 Lessons 1 19
Eureka Math 2015 2016 Algebra I Module 2 Lessons 1 19 Eureka Math, Published by the non-profit Great Minds. Copyright 2015 Great Minds. No part of this work may be reproduced, distributed, modified, sold,
More informationRemoving the Pattern Noise from all STIS Side-2 CCD data
The 2010 STScI Calibration Workshop Space Telescope Science Institute, 2010 Susana Deustua and Cristina Oliveira, eds. Removing the Pattern Noise from all STIS Side-2 CCD data Rolf A. Jansen, Rogier Windhorst,
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 informationGetting started with Spike Recorder on PC/Mac/Linux
Getting started with Spike Recorder on PC/Mac/Linux You can connect your SpikerBox to your computer using either the blue laptop cable, or the green smartphone cable. How do I connect SpikerBox to computer
More information2. AN INTROSPECTION OF THE MORPHING PROCESS
1. INTRODUCTION Voice morphing means the transition of one speech signal into another. Like image morphing, speech morphing aims to preserve the shared characteristics of the starting and final signals,
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 informationBitWise (V2.1 and later) includes features for determining AP240 settings and measuring the Single Ion Area.
BitWise. Instructions for New Features in ToF-AMS DAQ V2.1 Prepared by Joel Kimmel University of Colorado at Boulder & Aerodyne Research Inc. Last Revised 15-Jun-07 BitWise (V2.1 and later) includes features
More informationMore About Regression
Regression Line for the Sample Chapter 14 More About Regression is spoken as y-hat, and it is also referred to either as predicted y or estimated y. b 0 is the intercept of the straight line. The intercept
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 information1 Ver.mob Brief guide
1 Ver.mob 14.02.2017 Brief guide 2 Contents Introduction... 3 Main features... 3 Hardware and software requirements... 3 The installation of the program... 3 Description of the main Windows of the program...
More informationMachinery Diagnostic Plots Part 2 ORBIT Back-to-Basics: What does the data really tell us?
Machinery Diagnostic Plots Part 2 ORBIT Back-to-Basics: What does the data really tell us? Gaston Desimone Latin America Technical Leader Bently Nevada* Machinery Diagnostic Services (MDS) Buenos Aires
More informationMATLAB Basics 6 plotting
1 MATLAB Basics 6 plotting Anthony Rossiter University of Sheffield For a neat organisation of all videos and resources http://controleducation.group.shef.ac.uk/indexwebbook.html Introduction 2 1. The
More informationThe Measurement Tools and What They Do
2 The Measurement Tools The Measurement Tools and What They Do JITTERWIZARD The JitterWizard is a unique capability of the JitterPro package that performs the requisite scope setup chores while simplifying
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 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 informationDATA COMPRESSION USING THE FFT
EEE 407/591 PROJECT DUE: NOVEMBER 21, 2001 DATA COMPRESSION USING THE FFT INSTRUCTOR: DR. ANDREAS SPANIAS TEAM MEMBERS: IMTIAZ NIZAMI - 993 21 6600 HASSAN MANSOOR - 993 69 3137 Contents TECHNICAL BACKGROUND...
More informationDISPLAY WEEK 2015 REVIEW AND METROLOGY ISSUE
DISPLAY WEEK 2015 REVIEW AND METROLOGY ISSUE Official Publication of the Society for Information Display www.informationdisplay.org Sept./Oct. 2015 Vol. 31, No. 5 frontline technology Advanced Imaging
More informationA Parametric Autoregressive Model for the Extraction of Electric Network Frequency Fluctuations in Audio Forensic Authentication
Proceedings of the 3 rd International Conference on Control, Dynamic Systems, and Robotics (CDSR 16) Ottawa, Canada May 9 10, 2016 Paper No. 110 DOI: 10.11159/cdsr16.110 A Parametric Autoregressive Model
More informationFrequencies. Chapter 2. Descriptive statistics and charts
An analyst usually does not concentrate on each individual data values but would like to have a whole picture of how the variables distributed. In this chapter, we will introduce some tools to tabulate
More informationCS229 Project Report Polyphonic Piano Transcription
CS229 Project Report Polyphonic Piano Transcription Mohammad Sadegh Ebrahimi Stanford University Jean-Baptiste Boin Stanford University sadegh@stanford.edu jbboin@stanford.edu 1. Introduction In this project
More informationNanoGiant Oscilloscope/Function-Generator Program. Getting Started
Getting Started Page 1 of 17 NanoGiant Oscilloscope/Function-Generator Program Getting Started This NanoGiant Oscilloscope program gives you a small impression of the capabilities of the NanoGiant multi-purpose
More informationPulseCounter Neutron & Gamma Spectrometry Software Manual
PulseCounter Neutron & Gamma Spectrometry Software Manual MAXIMUS ENERGY CORPORATION Written by Dr. Max I. Fomitchev-Zamilov Web: maximus.energy TABLE OF CONTENTS 0. GENERAL INFORMATION 1. DEFAULT SCREEN
More informationEE 261 The Fourier Transform and its Applications Fall 2007 Problem Set Two Due Wednesday, October 10
EE 6 The Fourier Transform and its Applications Fall 007 Problem Set Two Due Wednesday, October 0. (5 points) A periodic, quadratic function and some surprising applications Let f(t) be a function of period
More informationMPEGTool: An X Window Based MPEG Encoder and Statistics Tool 1
MPEGTool: An X Window Based MPEG Encoder and Statistics Tool 1 Toshiyuki Urabe Hassan Afzal Grace Ho Pramod Pancha Magda El Zarki Department of Electrical Engineering University of Pennsylvania Philadelphia,
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 informationEASY-MCS. Multichannel Scaler. Profiling Counting Rates up to 150 MHz with 15 ppm Time Resolution.
Multichannel Scaler Profiling Counting Rates up to 150 MHz with 15 ppm Time Resolution. The ideal solution for: Time-resolved single-photon counting Phosphorescence lifetime spectrometry Atmospheric and
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 informationE X P E R I M E N T 1
E X P E R I M E N T 1 Getting to Know Data Studio Produced by the Physics Staff at Collin College Copyright Collin College Physics Department. All Rights Reserved. University Physics, Exp 1: Getting to
More informationThe Research of Controlling Loudness in the Timbre Subjective Perception Experiment of Sheng
The Research of Controlling Loudness in the Timbre Subjective Perception Experiment of Sheng S. Zhu, P. Ji, W. Kuang and J. Yang Institute of Acoustics, CAS, O.21, Bei-Si-huan-Xi Road, 100190 Beijing,
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 informationm RSC Chromatographie Integration Methods Second Edition CHROMATOGRAPHY MONOGRAPHS Norman Dyson Dyson Instruments Ltd., UK
m RSC CHROMATOGRAPHY MONOGRAPHS Chromatographie Integration Methods Second Edition Norman Dyson Dyson Instruments Ltd., UK THE ROYAL SOCIETY OF CHEMISTRY Chapter 1 Measurements and Models The Basic Measurements
More informationAn Effective Filtering Algorithm to Mitigate Transient Decaying DC Offset
An Effective Filtering Algorithm to Mitigate Transient Decaying DC Offset By: Abouzar Rahmati Authors: Abouzar Rahmati IS-International Services LLC Reza Adhami University of Alabama in Huntsville April
More informationStaMPS Persistent Scatterer Practical
StaMPS Persistent Scatterer Practical ESA Land Training Course, Leicester, 10-14 th September, 2018 Andy Hooper, University of Leeds a.hooper@leeds.ac.uk This practical exercise consists of working through
More informationSupplemental Material for Gamma-band Synchronization in the Macaque Hippocampus and Memory Formation
Supplemental Material for Gamma-band Synchronization in the Macaque Hippocampus and Memory Formation Michael J. Jutras, Pascal Fries, Elizabeth A. Buffalo * *To whom correspondence should be addressed.
More informationNoise. CHEM 411L Instrumental Analysis Laboratory Revision 2.0
CHEM 411L Instrumental Analysis Laboratory Revision 2.0 Noise In this laboratory exercise we will determine the Signal-to-Noise (S/N) ratio for an IR spectrum of Air using a Thermo Nicolet Avatar 360 Fourier
More informationMultiple-Window Spectrogram of Peaks due to Transients in the Electroencephalogram
284 IEEE TRANSACTIONS ON BIOMEDICAL ENGINEERING, VOL. 48, NO. 3, MARCH 2001 Multiple-Window Spectrogram of Peaks due to Transients in the Electroencephalogram Maria Hansson*, Member, IEEE, and Magnus Lindgren
More informationGetting Started. Connect green audio output of SpikerBox/SpikerShield using green cable to your headphones input on iphone/ipad.
Getting Started First thing you should do is to connect your iphone or ipad to SpikerBox with a green smartphone cable. Green cable comes with designators on each end of the cable ( Smartphone and SpikerBox
More informationSetting Up the Warp System File: Warp Theater Set-up.doc 25 MAY 04
Setting Up the Warp System File: Warp Theater Set-up.doc 25 MAY 04 Initial Assumptions: Theater geometry has been calculated and the screens have been marked with fiducial points that represent the limits
More informationTempo and Beat Analysis
Advanced Course Computer Science Music Processing Summer Term 2010 Meinard Müller, Peter Grosche Saarland University and MPI Informatik meinard@mpi-inf.mpg.de Tempo and Beat Analysis Musical Properties:
More informationDetection and demodulation of non-cooperative burst signal Feng Yue 1, Wu Guangzhi 1, Tao Min 1
International Conference on Applied Science and Engineering Innovation (ASEI 2015) Detection and demodulation of non-cooperative burst signal Feng Yue 1, Wu Guangzhi 1, Tao Min 1 1 China Satellite Maritime
More informationQuick-Start for READ30
Quick-Start for READ30 The program READ30 was written for the purpose of reading and configuring the digital pressure-transmitter of the series 30. The two features are divided into the following parts:
More informationAnalysis of WFS Measurements from first half of 2004
Analysis of WFS Measurements from first half of 24 (Report4) Graham Cox August 19, 24 1 Abstract Described in this report is the results of wavefront sensor measurements taken during the first seven months
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 informationPre-Processing of ERP Data. Peter J. Molfese, Ph.D. Yale University
Pre-Processing of ERP Data Peter J. Molfese, Ph.D. Yale University Before Statistical Analyses, Pre-Process the ERP data Planning Analyses Waveform Tools Types of Tools Filter Segmentation Visual Review
More informationNormalization Methods for Two-Color Microarray Data
Normalization Methods for Two-Color Microarray Data 1/13/2009 Copyright 2009 Dan Nettleton What is Normalization? Normalization describes the process of removing (or minimizing) non-biological variation
More informationWhat is Statistics? 13.1 What is Statistics? Statistics
13.1 What is Statistics? What is Statistics? The collection of all outcomes, responses, measurements, or counts that are of interest. A portion or subset of the population. Statistics Is the science of
More informationA Parametric Autoregressive Model for the Extraction of Electric Network Frequency Fluctuations in Audio Forensic Authentication
Journal of Energy and Power Engineering 10 (2016) 504-512 doi: 10.17265/1934-8975/2016.08.007 D DAVID PUBLISHING A Parametric Autoregressive Model for the Extraction of Electric Network Frequency Fluctuations
More informationA Matlab toolbox for. Characterisation Of Recorded Underwater Sound (CHORUS) USER S GUIDE
Centre for Marine Science and Technology A Matlab toolbox for Characterisation Of Recorded Underwater Sound (CHORUS) USER S GUIDE Version 5.0b Prepared for: Centre for Marine Science and Technology Prepared
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 informationPole Zero Correction using OBSPY and PSN Data
Pole Zero Correction using OBSPY and PSN Data Obspy provides the possibility of instrument response correction. WinSDR and WinQuake already have capability to embed the required information into the event
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 informationECE 4220 Real Time Embedded Systems Final Project Spectrum Analyzer
ECE 4220 Real Time Embedded Systems Final Project Spectrum Analyzer by: Matt Mazzola 12222670 Abstract The design of a spectrum analyzer on an embedded device is presented. The device achieves minimum
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 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 informationBlueline, Linefree, Accuracy Ratio, & Moving Absolute Mean Ratio Charts
INTRODUCTION This instruction manual describes for users of the Excel Standard Celeration Template(s) the features of each page or worksheet in the template, allowing the user to set up and generate charts
More informationExperiment 2: Sampling and Quantization
ECE431, Experiment 2, 2016 Communications Lab, University of Toronto Experiment 2: Sampling and Quantization Bruno Korst - bkf@comm.utoronto.ca Abstract In this experiment, you will see the effects caused
More informationEE391 Special Report (Spring 2005) Automatic Chord Recognition Using A Summary Autocorrelation Function
EE391 Special Report (Spring 25) Automatic Chord Recognition Using A Summary Autocorrelation Function Advisor: Professor Julius Smith Kyogu Lee Center for Computer Research in Music and Acoustics (CCRMA)
More informationStaMPS Persistent Scatterer Exercise
StaMPS Persistent Scatterer Exercise ESA Land Training Course, Bucharest, 14-18 th September, 2015 Andy Hooper, University of Leeds a.hooper@leeds.ac.uk This exercise consists of working through an example
More informationUNIVERSAL SPATIAL UP-SCALER WITH NONLINEAR EDGE ENHANCEMENT
UNIVERSAL SPATIAL UP-SCALER WITH NONLINEAR EDGE ENHANCEMENT Stefan Schiemenz, Christian Hentschel Brandenburg University of Technology, Cottbus, Germany ABSTRACT Spatial image resizing is an important
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 informationpotentiostat/galvanostat
potentiostat/galvanostat Rev. 12-2012 potentiostat/galvanostat A battery-powered, handheld instrument which allows the application of most of the relevant voltammetric and amperometric techniques. The
More informationMusic Segmentation Using Markov Chain Methods
Music Segmentation Using Markov Chain Methods Paul Finkelstein March 8, 2011 Abstract This paper will present just how far the use of Markov Chains has spread in the 21 st century. We will explain some
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 informationTERRESTRIAL broadcasting of digital television (DTV)
IEEE TRANSACTIONS ON BROADCASTING, VOL 51, NO 1, MARCH 2005 133 Fast Initialization of Equalizers for VSB-Based DTV Transceivers in Multipath Channel Jong-Moon Kim and Yong-Hwan Lee Abstract This paper
More informationGG450 4/12/2010. Today s material comes from p in the text book. Please read and understand all of this material!
GG450 April 13, 2010 Seismic Reflection III Data Processing Today s material comes from p. 163-198 in the text book. Please read and understand all of this material! Reflection Processing We've been talking
More informationElectrospray-MS Charge Deconvolutions without Compromise an Enhanced Data Reconstruction Algorithm utilising Variable Peak Modelling
Electrospray-MS Charge Deconvolutions without Compromise an Enhanced Data Reconstruction Algorithm utilising Variable Peak Modelling Overview A.Ferrige1, S.Ray1, R.Alecio1, S.Ye2 and K.Waddell2 1 PPL,
More informationCSC475 Music Information Retrieval
CSC475 Music Information Retrieval Monophonic pitch extraction George Tzanetakis University of Victoria 2014 G. Tzanetakis 1 / 32 Table of Contents I 1 Motivation and Terminology 2 Psychacoustics 3 F0
More informationRobert Alexandru Dobre, Cristian Negrescu
ECAI 2016 - International Conference 8th Edition Electronics, Computers and Artificial Intelligence 30 June -02 July, 2016, Ploiesti, ROMÂNIA Automatic Music Transcription Software Based on Constant Q
More informationPracticum 3, Fall 2010
A. F. Miller 2010 T1 Measurement 1 Practicum 3, Fall 2010 Measuring the longitudinal relaxation time: T1. Strychnine, dissolved CDCl3 The T1 is the characteristic time of relaxation of Z magnetization
More informationUSB Mini Spectrum Analyzer User s Guide TSA5G35
USB Mini Spectrum Analyzer User s Guide TSA5G35 Triarchy Technologies, Corp. Page 1 of 21 USB Mini Spectrum Analyzer User s Guide Copyright Notice Copyright 2011 Triarchy Technologies, Corp. All rights
More informationVisual Encoding Design
CSE 442 - Data Visualization Visual Encoding Design Jeffrey Heer University of Washington A Design Space of Visual Encodings Mapping Data to Visual Variables Assign data fields (e.g., with N, O, Q types)
More informationThe XYZ Colour Space. 26 January 2011 WHITE PAPER. IMAGE PROCESSING TECHNIQUES
www.omnitek.tv IMAE POESSIN TEHNIQUES The olour Space The colour space has the unique property of being able to express every colour that the human eye can see which in turn means that it can express every
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 informationMestReNova A quick Guide. Adjust signal intensity Use scroll wheel. Zoomen Z
MestReNova A quick Guide page 1 MNova is a program to analyze 1D- and 2D NMR data. Start of MNova Start All Programs Chemie NMR MNova The MNova Menu 1. 2. Create expanded regions Adjust signal intensity
More informationResearch Article. ISSN (Print) *Corresponding author Shireen Fathima
Scholars Journal of Engineering and Technology (SJET) Sch. J. Eng. Tech., 2014; 2(4C):613-620 Scholars Academic and Scientific Publisher (An International Publisher for Academic and Scientific Resources)
More informationResults of the June 2000 NICMOS+NCS EMI Test
Results of the June 2 NICMOS+NCS EMI Test S. T. Holfeltz & Torsten Böker September 28, 2 ABSTRACT We summarize the findings of the NICMOS+NCS EMI Tests conducted at Goddard Space Flight Center in June
More informationAnalysis. mapans MAP ANalysis Single; map viewer, opens and modifies a map file saved by iman.
Analysis Analysis routines (run on LINUX): iman IMage ANalysis; makes maps out of raw data files saved be the acquisition program (ContImage), can make movies, pictures of green, compresses and decompresses
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 informationMoving on from MSTAT. March The University of Reading Statistical Services Centre Biometrics Advisory and Support Service to DFID
Moving on from MSTAT March 2000 The University of Reading Statistical Services Centre Biometrics Advisory and Support Service to DFID Contents 1. Introduction 3 2. Moving from MSTAT to Genstat 4 2.1 Analysis
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 informationLecture 2 Video Formation and Representation
2013 Spring Term 1 Lecture 2 Video Formation and Representation Wen-Hsiao Peng ( 彭文孝 ) Multimedia Architecture and Processing Lab (MAPL) Department of Computer Science National Chiao Tung University 1
More information