Plans for stochastic hardware injection
Download
1 / 25

Plans for stochastic hardware injection - PowerPoint PPT Presentation


  • 96 Views
  • Uploaded on

Plans for stochastic hardware injection. Eric Thrane, Mandy Pihlaja and Vuk Mandic. Overview. Motivation Injection flow chart Relationship to reviewed code Testing and validation. Summary of injection. spherical harmonic coefficients and spectrum specified

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about ' Plans for stochastic hardware injection' - inoke


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript
Plans for stochastic hardware injection

Plans for stochastic hardware injection

Eric Thrane, Mandy Pihlaja and Vuk Mandic


Overview
Overview

  • Motivation

  • Injection flow chart

  • Relationship to reviewed code

  • Testing and validation


Summary of injection
Summary of injection

  • spherical harmonic coefficients and spectrum specified

  • User sets parameters of GPS, random seed, injection duration

  • Parent routine, SphMapToTimeSeries, calls simulateSkyMapTimeDomain, which returns time series [d1, d2] and overlapping time series [Ovrlp1, Ovrlp2] for {LHO,LLO}.

  • simulateSkyMapTimeDomain = reviewed stochastic code


Summary cont
Summary cont.

  • Then…

    • FFT both time series

    • Apply frequency pendulum transfer function

    • Window time series using sin function

    • Combine time series and overlapping time series

  • Write combined h(t) to text file

  • sbaniso.c (based on standard injection code, sb.c) reads in text file and converts to binary

  • Sent to psinject

  • Sent to awg


Injection flowchart
Injection flowchart

*new

user-specified params: random seeds, seg dur, etc.

SpH map file: l, m, Plm

spectrum file: f, H(f)

gets time series

gets overlapping time series

FFT time series

Apply pendulum T(f)

IFFT

Window

Combine overlapping series

Write combined h(t) to files.

SpHMapToTimeSeries*

d1,d2

simulateSkyMapTimeDomain

(reviewed)

pg. 6

binary out

10110010

awg

psinject

sbaniso.c

pg. 23


Simulateskymaptimeseries m
simulateSkyMapTimeSeries.m

  • Returns [d1(t),d2(t)] for LHO, LLO given user-specified P() and H(f).

  • Written by Stefan Balmer and Joe Romano

  • Reviewed during SpH review.


Validation and testing fft
Validation and testing- fft

  • The next four slides show the process of applying FFT, the pendulum transfer function, and the iFFT to the time series.

  • Then there is a consistency check to make sure this process could be reversed


Original Time Series

strain

time(s)


Original time series
Original time series

strain

time(s)


Fft of time series
FFT of time series

fft

frequency(Hz)



Pendulum transfer function
Pendulum Transfer Function function is given by sb.c

P(f)

frequency(Hz)


Fft multiplied by tfh
FFT multiplied by tfH function is given by sb.c

Time series

frequency(Hz)


Ifft back to time series
ifft back to time series function is given by sb.c

Time series

time(s)


  • In order to check that the pendulum transfer function was applied correctly, the transformed time series had the following inverse process applied

    • FFT, inverse pendulum transfer function, iFFT

  • The plot of the ratio of the ‘inverted’ time series to the original time series is on the next slide.

  • Within very small rounding, the original time series is recovered (~1.0)


Ratio of two time series
Ratio of two time series applied correctly, the transformed time series had the following inverse process applied

ratio

1.07

1

time(s)


Windowing
Windowing applied correctly, the transformed time series had the following inverse process applied

  • Each time series segment is multiplied by the sine function in order for beginning and end to scale to zero.


Windowed time series segment
Windowed time series segment applied correctly, the transformed time series had the following inverse process applied

counts

time(s)


Overlapping and combination
Overlapping and combination applied correctly, the transformed time series had the following inverse process applied

  • The overlapping time series lags the original time series by half of segment duration.

  • The second half of the current overlapping time series and the first half of the previous overlapping time series are concatenated to make the combined overlapping time series

  • The combined series is added to the original time series to form the final time series


Final time series
Final time series applied correctly, the transformed time series had the following inverse process applied

counts

time(s)


  • No modulation of the envelope applied correctly, the transformed time series had the following inverse process applied

  • Rms value of the pre-windowed time series = 2.2348

  • Rms value of the combined time series = 2.2334

  • The close rms values of the two time series shows apparent agreement


Checking code at lho and llo
Checking code at LHO and LLO applied correctly, the transformed time series had the following inverse process applied

  • The code was run at LLO and LHO given the same random seed to check that the time series were the same

  • Data in Xcel sheet (next page)


sb.c applied correctly, the transformed time series had the following inverse process applied

  • Originally generates time series, applies transfer function, overlaps segments and sends to AWG

  • Modified to read in text files created by SphMapToTimeSeries.m, send to AWG, and delete file after use (sbaniso.c)

  • Already used for S5 (w/o Vuk’s hacks) and reviewed


Handshake
Handshake applied correctly, the transformed time series had the following inverse process applied

  • Day long injection generated in one minute segments

  • In order to not create a 1TB injection file, SphMapToTimeSeries.m only generates three files at a time and waits for sbaniso.c to delete a file before continuing.


ad