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 | |
---|