| 1 | command SPIKING <trc> <lo-bound> <hi-bound> <reg> |
|---|
| 2 | =============== |
|---|
| 3 | |
|---|
| 4 | key: compute spiking filter |
|---|
| 5 | |
|---|
| 6 | Creates spiking filter (inverse wavelet) of a trace window given by |
|---|
| 7 | <trc>, <lo-bound> and <hi-bound>. The time response function of the |
|---|
| 8 | filter is added to the display. You can use the "fold" command to |
|---|
| 9 | apply this filter to other traces. |
|---|
| 10 | |
|---|
| 11 | |
|---|
| 12 | parameters |
|---|
| 13 | ---------- |
|---|
| 14 | |
|---|
| 15 | <trc> --- parameter type: trace |
|---|
| 16 | Trace to compute spiking filter |
|---|
| 17 | |
|---|
| 18 | <lo-bound>, <hi-bound> --- parameter type: real |
|---|
| 19 | Time window on trace <trc> for computation of spiking filter |
|---|
| 20 | |
|---|
| 21 | <reg> --- parameter type: real |
|---|
| 22 | Control parameter of spiking. The computation is done in time |
|---|
| 23 | domain using the toeplitz matrix made from the autocorrelation |
|---|
| 24 | of the specified window. Before evaluation the diagonal |
|---|
| 25 | elements of the toeplitz matrix are multiplied by <reg>+1. |
|---|
| 26 | <reg>=0 means ordinary inversion. <reg> > 0 means low-frequent |
|---|
| 27 | inversion (prewhitening). |
|---|
| 28 | |
|---|
| 29 | |
|---|
| 30 | qualifiers |
|---|
| 31 | ---------- |
|---|
| 32 | |
|---|
| 33 | /pos=<spike-pos> |
|---|
| 34 | Specifies the spike position of the desired trace in seconds |
|---|
| 35 | relative to the start of the time window <lo-bound>. As a |
|---|
| 36 | default, the spike is positioned at the center of the signal |
|---|
| 37 | <spike-pos> = sum_i{i*dt*a[i]} / sum_i{a[i]*a[i]}, |
|---|
| 38 | where a[i] are the amplitude values in the window. |
|---|
| 39 | You also can specify "/pos=max", then the spike is positioned |
|---|
| 40 | at the maximum amplitude inside the window. |
|---|
| 41 | |
|---|
| 42 | /lth=<lth-fac> |
|---|
| 43 | Determines the length of the filter trace in units of the |
|---|
| 44 | length of the specified window. Default is 0.5, i.e. the |
|---|
| 45 | filter has half the length of the signal window. |
|---|
| 46 | |
|---|
| 47 | |
|---|
| 48 | examples |
|---|
| 49 | -------- |
|---|
| 50 | |
|---|
| 51 | spiking 1 5.0 50.0 1.0 ! computes spiking filter of trace 1 |
|---|
| 52 | ! in window from 5.0 sec to 50.0 sec. |
|---|
| 53 | ! A control parameter of 1.0 is used |
|---|
| 54 | |
|---|
| 55 | spiking/lth=1.0 2 *;;0.0 ! computes spiking filter of trace 2; |
|---|
| 56 | ! the time window must be selected by |
|---|
| 57 | ! graphic cursor. The output trace |
|---|
| 58 | ! has the same length as the selected |
|---|
| 59 | ! window. Control parameter is 0.0. |
|---|
| 60 | |
|---|
