Create a Code Snippet by Logging In or Creating an Account

Wave Folding

Thu, 01/17/2019 - 10:29 am

This small snippet allows one to efficiently fold and average a 1D wave containing periodic data, for instances, a wave of 100 points containing ten cycles of 10-point data.

The end result is a 10-point 1D wave containing the averaged results of the ten cycles originally contained by the 100-point 1D wave.

function foldWave(foldedWave, waveToFold)

    wave foldedWave, waveToFold

    variable numRowfW = numpnts(foldedWave)

    variable numRowwTF = numpnts(waveToFold)

    variable numFold = floor(numRowwTF/numRowfw)

    duplicate/o/R=[0,numFold*numRowfW] waveToFold, fWMatrix

    redimension/N=(numRowfW, numFold,0,0) fWMatrix

    multithread foldedWave = meanRow(fwMatrix,p)

end

threadsafe static Function MeanRow (wave2D, whichRow)

    wave wave2D

    variable whichRow       

    MatrixOP/O/free w=row(wave2D,whichRow)^t

   return mean(w)

End