Stats

Locked Locked
Replies 12
Subscribers 125
Views 17519
Members are here 0

This discussion has been locked.

You can no longer post new replies to this discussion. If you have a question you can start a new discussion

Manipulate and export data from ViVa

BaaB over 9 years ago

I have two signals as in the picture below.

I would like to measure some information and export them to a table or a file something like this:

Pulse width fall time

614n 738n

696n 686n

581n 660n

660n 672n

... ...

I am wondering if there is a way to do that. It would be a great help. Otherwise, I have to do it manually with hundreds points.

Thank you.

Parents

Andrew Beckett over 9 years ago

This code should do it. See the example in the comments. Making this flexible is what took a bit of time. It probably could do with a bit more testing:

/* abGetCyclePeaks.il

Author     A.D.Beckett
Group      Custom IC (UK), Cadence Design Systems Ltd.
Language   SKILL
Date       Mar 23, 2016 
Modified   
By         

Function to calculate the peak values (or times of peak values)
within each period. Has some tolerances built in so that it should
detects when it first reaches the peak and filters out numerical
noise.

Example - measure fall time (peak to peak) in each period:

abGetCyclePeaks(v("ramp" ?result 'tran) "rising" "max" ?yValues "x")-
    abGetCyclePeaks(v("ramp" ?result 'tran) "rising" "min" ?yValues "x")

***************************************************

SCCS Info: @(#) abGetCyclePeaks.il 03/23/16.22:27:05 1.1

*/

/*******************************************************************
*                                                                  *
*      (abGetCyclePeaks wave crossType minMax [?mode "auto"]       *
*         [?threshold 0.0] [?yValues "y"] [?xName "cycle"]         *
*                  [?reltol 1e-3] [?abstol 1e-6])                  *
*                                                                  *
*     For each cycle - identified by the rising or falling (or     *
*     either, but that's not so likely to be used) crossing of     *
*     the threshold (which is controlled by the ?mode argument     *
*  "auto" means automatic, "user" means as specified) - find the   *
* "min" or "max" value of the signal within the cycle. In order to *
* filter out numerical noise on the signal influencing the precise *
*     position of the peak, only look at changes bigger than a     *
*   reltol+abstol type delta from the previously found peak. The   *
*   resulting waveform can be plotted versus "cycle" or "time",    *
*    and the yValues can be either the y or x values at the peak.  *
*                                                                  *
*******************************************************************/

(defun abGetCyclePeaks (wave crossType minMax @key (mode "auto") 
                             (threshold 0.0) (yValues "y") (xName "cycle") 
                             (reltol 1e-3) (abstol 1e-6) "gSSSfSSff")
  (cond
    ((drIsWaveform wave)
     (let (xVec yVec len outX outY newXVec newYVec newWave cycleNum lastY 
                thisY localPeakX localPeakY direction func absdelta tolerance)
       ;-----------------------------------------------------------------
       ; Validate the arguments to ensure they are valid strings
       ;-----------------------------------------------------------------
       (foreach validInfo
             `(
               (crossType ,(setq crossType (get_string crossType))
                          ("rising" "falling" "either"))
               (minMax ,(setq minMax (get_string minMax)) ("min" "max"))
               (?mode ,(setq mode (get_string mode)) ("auto" "user"))
               (?xName ,(setq xName (get_string xName)) ("cycle" "time"))
               (?yValues ,(setq yValues (get_string yValues)) ("y" "x"))
               )
             (unless (member (cadr validInfo) (caddr validInfo))
               (error "abGetCyclePeaks: %s argument must be one of \"%s\": %L\n"
                      (car validInfo)
                      (buildString (caddr validInfo) "\",\"")
                      (cadr validInfo))
               )
             )
       (when (equal mode "auto")
         (setq threshold (average wave)))
       (setq func (if (equal minMax "max") 'greaterp 'lessp))
       (setq xVec (drGetWaveformXVec wave))
       (setq yVec (drGetWaveformYVec wave))
       (setq len (drVectorLength xVec))
       (setq cycleNum 0)
       (setq lastY (drGetElem yVec 0))
       (setq localPeakY lastY)
       (setq localPeakX (drGetElem xVec 0))
       ;-----------------------------------------------------------------
       ; Iterate through all the points
       ;-----------------------------------------------------------------
       (for i 1 (sub1 len)
            (setq thisY (drGetElem yVec i))
            ;------------------------------------------------------------
            ; Determine if this is a transition, and if so, which direction
            ;------------------------------------------------------------
            (setq direction
                  (cond 
                    ((and (lessp lastY threshold)
                          (geqp thisY threshold)) "rising")
                    ((and (geqp lastY threshold)
                          (lessp thisY threshold)) "falling")
                    (t nil)))
            ;------------------------------------------------------------
            ; Assuming it's a crossing in the right direction, record
            ; the peak in that cycle (doesn't do it before the first
            ; transition though)
            ;------------------------------------------------------------
            (if (or (equal crossType direction)
                    (and (equal crossType "either") direction))
              (progn
                (when (greaterp cycleNum 0)
                  (setq outX 
                        (tconc outX (if (equal xName "cycle") cycleNum 
                                      localPeakX)))
                  (setq outY (tconc outY 
                                    (if (equal yValues "y") 
                                      localPeakY localPeakX)))
                  )
                (setq localPeakX (drGetElem xVec i))
                (setq localPeakY thisY)
                (postincrement cycleNum)
                )
              ;----------------------------------------------------------
              ; Otherwise see if the peak has been exceeded
              ;----------------------------------------------------------
              (when (funcall func thisY localPeakY)
                (setq absdelta (abs (difference thisY lastY)))
                (setq tolerance (plus (times (max (abs thisY) (abs lastY)) 
                                             reltol) abstol))
                (when (greaterp absdelta tolerance)
                  (setq localPeakX (drGetElem xVec i))
                  (setq localPeakY thisY)
                  )
                )
              )
            (setq lastY thisY)
            )
       ;-----------------------------------------------------------------
       ; Construct output vectors and waveforms with the right
       ; attributes
       ;-----------------------------------------------------------------
       (setq newXVec (drCreateVec 
                       (if (equal xName "cycle") 'intlong 'double) (car outX)))
       (setq newYVec (drCreateVec 'double (car outY)))
       (putpropq newXVec (if (equal xName "cycle") "" (getq xVec units)) units)
       (putpropq newXVec 
                 (if (equal xName "cycle") "cycle" (getq xVec name)) name)
       (putpropq newYVec 
                 (if (equal yValues "y") 
                   (getq yVec units) 
                   (getq xVec units)) 
                 units)
       (putpropq newYVec 
                 (if (equal yValues "y")
                   (getq yVec name)
                   (getq xVec name))
                 name)
       (setq newWave (drCreateWaveform newXVec newYVec))
       (famSetExpr newWave 
                   `(abGetCyclePeaks ,(famGetExpr wave) crossType minMax 
                                     ?mode ,mode ?threshold ,threshold 
                                     ?yValues yValues ?xName ,xName 
                                     ?reltol ,reltol ?abstol ,abstol))
       newWave
       )
     )
    ((famIsFamily wave)
     (famMap 'abGetCyclePeaks minMax ?mode mode ?threshold threshold 
             ?yValues yValues ?xName xName ?reltol reltol ?abstol abstol)
     )
    (t
      (error "abGetCyclePeaks: Cannot handle %L\n" wave)
      )
    )
  )

Reply

Andrew Beckett over 9 years ago

This code should do it. See the example in the comments. Making this flexible is what took a bit of time. It probably could do with a bit more testing:

/* abGetCyclePeaks.il

Author     A.D.Beckett
Group      Custom IC (UK), Cadence Design Systems Ltd.
Language   SKILL
Date       Mar 23, 2016 
Modified   
By         

Function to calculate the peak values (or times of peak values)
within each period. Has some tolerances built in so that it should
detects when it first reaches the peak and filters out numerical
noise.

Example - measure fall time (peak to peak) in each period:

abGetCyclePeaks(v("ramp" ?result 'tran) "rising" "max" ?yValues "x")-
    abGetCyclePeaks(v("ramp" ?result 'tran) "rising" "min" ?yValues "x")

***************************************************

SCCS Info: @(#) abGetCyclePeaks.il 03/23/16.22:27:05 1.1

*/

/*******************************************************************
*                                                                  *
*      (abGetCyclePeaks wave crossType minMax [?mode "auto"]       *
*         [?threshold 0.0] [?yValues "y"] [?xName "cycle"]         *
*                  [?reltol 1e-3] [?abstol 1e-6])                  *
*                                                                  *
*     For each cycle - identified by the rising or falling (or     *
*     either, but that's not so likely to be used) crossing of     *
*     the threshold (which is controlled by the ?mode argument     *
*  "auto" means automatic, "user" means as specified) - find the   *
* "min" or "max" value of the signal within the cycle. In order to *
* filter out numerical noise on the signal influencing the precise *
*     position of the peak, only look at changes bigger than a     *
*   reltol+abstol type delta from the previously found peak. The   *
*   resulting waveform can be plotted versus "cycle" or "time",    *
*    and the yValues can be either the y or x values at the peak.  *
*                                                                  *
*******************************************************************/

(defun abGetCyclePeaks (wave crossType minMax @key (mode "auto") 
                             (threshold 0.0) (yValues "y") (xName "cycle") 
                             (reltol 1e-3) (abstol 1e-6) "gSSSfSSff")
  (cond
    ((drIsWaveform wave)
     (let (xVec yVec len outX outY newXVec newYVec newWave cycleNum lastY 
                thisY localPeakX localPeakY direction func absdelta tolerance)
       ;-----------------------------------------------------------------
       ; Validate the arguments to ensure they are valid strings
       ;-----------------------------------------------------------------
       (foreach validInfo
             `(
               (crossType ,(setq crossType (get_string crossType))
                          ("rising" "falling" "either"))
               (minMax ,(setq minMax (get_string minMax)) ("min" "max"))
               (?mode ,(setq mode (get_string mode)) ("auto" "user"))
               (?xName ,(setq xName (get_string xName)) ("cycle" "time"))
               (?yValues ,(setq yValues (get_string yValues)) ("y" "x"))
               )
             (unless (member (cadr validInfo) (caddr validInfo))
               (error "abGetCyclePeaks: %s argument must be one of \"%s\": %L\n"
                      (car validInfo)
                      (buildString (caddr validInfo) "\",\"")
                      (cadr validInfo))
               )
             )
       (when (equal mode "auto")
         (setq threshold (average wave)))
       (setq func (if (equal minMax "max") 'greaterp 'lessp))
       (setq xVec (drGetWaveformXVec wave))
       (setq yVec (drGetWaveformYVec wave))
       (setq len (drVectorLength xVec))
       (setq cycleNum 0)
       (setq lastY (drGetElem yVec 0))
       (setq localPeakY lastY)
       (setq localPeakX (drGetElem xVec 0))
       ;-----------------------------------------------------------------
       ; Iterate through all the points
       ;-----------------------------------------------------------------
       (for i 1 (sub1 len)
            (setq thisY (drGetElem yVec i))
            ;------------------------------------------------------------
            ; Determine if this is a transition, and if so, which direction
            ;------------------------------------------------------------
            (setq direction
                  (cond 
                    ((and (lessp lastY threshold)
                          (geqp thisY threshold)) "rising")
                    ((and (geqp lastY threshold)
                          (lessp thisY threshold)) "falling")
                    (t nil)))
            ;------------------------------------------------------------
            ; Assuming it's a crossing in the right direction, record
            ; the peak in that cycle (doesn't do it before the first
            ; transition though)
            ;------------------------------------------------------------
            (if (or (equal crossType direction)
                    (and (equal crossType "either") direction))
              (progn
                (when (greaterp cycleNum 0)
                  (setq outX 
                        (tconc outX (if (equal xName "cycle") cycleNum 
                                      localPeakX)))
                  (setq outY (tconc outY 
                                    (if (equal yValues "y") 
                                      localPeakY localPeakX)))
                  )
                (setq localPeakX (drGetElem xVec i))
                (setq localPeakY thisY)
                (postincrement cycleNum)
                )
              ;----------------------------------------------------------
              ; Otherwise see if the peak has been exceeded
              ;----------------------------------------------------------
              (when (funcall func thisY localPeakY)
                (setq absdelta (abs (difference thisY lastY)))
                (setq tolerance (plus (times (max (abs thisY) (abs lastY)) 
                                             reltol) abstol))
                (when (greaterp absdelta tolerance)
                  (setq localPeakX (drGetElem xVec i))
                  (setq localPeakY thisY)
                  )
                )
              )
            (setq lastY thisY)
            )
       ;-----------------------------------------------------------------
       ; Construct output vectors and waveforms with the right
       ; attributes
       ;-----------------------------------------------------------------
       (setq newXVec (drCreateVec 
                       (if (equal xName "cycle") 'intlong 'double) (car outX)))
       (setq newYVec (drCreateVec 'double (car outY)))
       (putpropq newXVec (if (equal xName "cycle") "" (getq xVec units)) units)
       (putpropq newXVec 
                 (if (equal xName "cycle") "cycle" (getq xVec name)) name)
       (putpropq newYVec 
                 (if (equal yValues "y") 
                   (getq yVec units) 
                   (getq xVec units)) 
                 units)
       (putpropq newYVec 
                 (if (equal yValues "y")
                   (getq yVec name)
                   (getq xVec name))
                 name)
       (setq newWave (drCreateWaveform newXVec newYVec))
       (famSetExpr newWave 
                   `(abGetCyclePeaks ,(famGetExpr wave) crossType minMax 
                                     ?mode ,mode ?threshold ,threshold 
                                     ?yValues yValues ?xName ,xName 
                                     ?reltol ,reltol ?abstol ,abstol))
       newWave
       )
     )
    ((famIsFamily wave)
     (famMap 'abGetCyclePeaks minMax ?mode mode ?threshold threshold 
             ?yValues yValues ?xName xName ?reltol reltol ?abstol abstol)
     )
    (t
      (error "abGetCyclePeaks: Cannot handle %L\n" wave)
      )
    )
  )

Children

No Data

Community Guidelines

The Cadence Design Communities support Cadence users and technologists interacting to exchange ideas, news, technical information, and best practices to solve problems and get the most from Cadence technology. The community is open to everyone, and to provide the most value, we require participants to follow our Community Guidelines that facilitate a quality exchange of ideas and information. By accessing, contributing, using or downloading any materials from the site, you agree to be bound by the full Community Guidelines.