A new method of time-frequency analysis, based on the Matching Pursuit (MP) algorithm, was used to extract and quantify EEG 'driving' or frequency-following responses produced in human primary somatosensory cortex (SI) by 33 Hz vibrotactile stimulation of the right index fingertip in a single subject. EEG signals were recorded from a 5 x 5 array of electrodes centered over the left hand area, time-locked to repeated presentations of four vibratory stimulus amplitudes. The MP algorithm was used to decompose the edited and and filtered EEG signals into waveforms selected from a large and redundant dictionary. Statistical discrimination of the vibratory stimulus amplitudes was then readily achieved in terms of trial-by-trial measures of response amplitude constructed in automated fashion from the calculated MP parameters. The results were orderly and physiologically coherent, and potentially open the way to correlation of psychophysical magnitude estimates with measures of neurophysiological response on a trial-by-trial basis. The approach developed here appears well suited to detection and characterisation of time dependent or transient target signals embedded in a noisy background.