DE10325147A1 - Signal analysis method for time frequency analysis of signal sequences uses atoms in a dictionary in trials for a signal sequence - Google Patents
Signal analysis method for time frequency analysis of signal sequences uses atoms in a dictionary in trials for a signal sequence Download PDFInfo
- Publication number
- DE10325147A1 DE10325147A1 DE2003125147 DE10325147A DE10325147A1 DE 10325147 A1 DE10325147 A1 DE 10325147A1 DE 2003125147 DE2003125147 DE 2003125147 DE 10325147 A DE10325147 A DE 10325147A DE 10325147 A1 DE10325147 A1 DE 10325147A1
- Authority
- DE
- Germany
- Prior art keywords
- signal
- atoms
- signal sequence
- dictionary
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Withdrawn
Links
Classifications
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R23/00—Arrangements for measuring frequencies; Arrangements for analysing frequency spectra
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/14—Fourier, Walsh or analogous domain transformations, e.g. Laplace, Hilbert, Karhunen-Loeve, transforms
- G06F17/148—Wavelet transforms
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/24—Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
- A61B5/316—Modalities, i.e. specific diagnostic methods
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7232—Signal processing specially adapted for physiological signals or for diagnostic purposes involving compression of the physiological signal, e.g. to extend the signal recording period
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
Description
Das Anwendungsgebiet der Erfindung ist die Zeit-Frequenz-Analyse von Signalfolgen. Dabei können diese Signalfolgen aus beliebigen Bereichen von Wissenschaft und Technik entstammen. Diese Signalfolgen entstehen bei jeder Messwertaufnahme und sind häufig durch lange Folgen von Messwerten gekennzeichnet. Biomedizinische Beispiele für diese Signalfolgen sind das menschliche Elektroenzephalogramm (EEG), das Magnetoenzephalogramm (MEG), das Elektrokardiogramm (EKG) und das Magnetokardiogramm (MKG). Andere Beispiele sind Audiodaten, wie Sprache, meteorologische Daten, geologische Daten oder Messdaten aus industriellen Prozessen.The Field of application of the invention is the time-frequency analysis of Signal sequences. You can do this Signal sequences from any area of science and technology come. These signal sequences arise with every measurement value recording and are common characterized by long sequences of measured values. Biomedical examples for these signal sequences are the human electroencephalogram (EEG), the magnetoencephalogram (MEG), the electrocardiogram (EKG) and the magnetocardiogram (MKG). Other examples are audio data, such as language, meteorological data, geological data or measurement data from industrial processes.
Bei diesen Signalfolgen werden zum Zwecke der Analyse, Klassifizierung, Filterung und/oder Komprimierung Zeit-Frequenz-Analysen durchgeführt. Derartige Zeit-Frequenz-Analysen ermöglichen über die klassische Signalamplitudendarstellung hinaus die Charakterisierung des Frequenzgehaltes zu bestimmten Zeiten. Das Frequenzspektrum eines Signals bezeichnet dabei einfach die Frequenzkomponenten des Signals. Dies kann durch die Fourier-Transformation (FT) des Signals beschrieben werden. Für das Signal x(t) ist die Fourier-Transformation definiert durch: Time-frequency analyzes are carried out on these signal sequences for the purpose of analysis, classification, filtering and / or compression. Such time-frequency analyzes enable the characterization of the frequency content at certain times in addition to the classic signal amplitude representation. The frequency spectrum of a signal simply denotes the frequency components of the signal. This can be described by the Fourier transform (FT) of the signal. For the signal x (t) the Fourier transformation is defined by:
Die Fourier-Transformation ist die Darstellung des Signals als gewichtete Summe der verschiedenen sinusförmigen Frequenzen. Die Information des Integrals beinhaltet alle Zeitmomente, da die Integration von minus unendlich bis plus unendlich ist. Daraus folgt, dass, ganz gleich wann, in der Zeit die Komponente mit der Frequenz "f" erscheint, das Ergebnis der Integration identisch beeinflusst wird. Die Transformation liefert die Information über den spektralen Inhalt des Signals, aber nicht darüber, wann spezielle Komponenten auftreten. Zusammengefasst bedeutet das, dass die Fourier-Transformation eine ausgezeichneten Frequenzauflösung, aber keine Zeitauflösung hat.The Fourier transform is the representation of the signal as weighted Sum of the different sinusoidal Frequencies. The information of the integral includes all moments of time, because the integration is from minus infinity to plus infinity. from that follows that, no matter when, in time the component with the Frequency "f" appears, the result the integration is influenced identically. The transformation delivers the information about the spectral content of the signal, but not when special components occur. In summary, that means that the Fourier transform has excellent frequency resolution, however no time resolution Has.
Um nun eine Auflösung in Zeit und Frequenz zu erhalten, werden sogenannte Multi-Resolution-Verfahren eingesetzt. Neben der Wavelet-Analyse, die aber eine schlechte Auflösung bei höheren Frequenzen besitzt, wird dazu das Matching Pursuit eingesetzt (z. B. S. Mallat and Z. Zhang, Matching pursuit with time-frequency dictionaries, IEEE Transactions on Signal Processing, 41(12), 3397–3414, 1993; G. Davis, S. G. Mallat & Z. Zhang: Adaptive time-frequency approximations with matching pursuit, in Wavelets: Theory, Algorithms and Applications, 1994; C. K. Chui, L. Montefusco and L. Puccio (Eds.), San Diego: Academic Press, 271–293; P. J. Durka: Time-frequency analyses of EEG, PhD dissertation, Warsaw University, Warsaw, Poland, 1996; P. J. Durka and K. J. Blinowska: Analysis of EEG transients by means of matching pursuit, Ann. Biomed. Eng., vol. 23, 1995, 608–611; Walden, A. T. & Contreras Cristan, A: Matching pursuit by undecimated discrete wavelet transform for non-stationary time series of arbitrary length. Statistics and Computing 8, 1998, 205–219; P. Vandergheynst and P. Frossard: Efficient image representation by anisotropic refinement in matching pursuit, in Proceedings of IEEE, volume 3, Salt Lake City UT, May 2001, ICASSP. 4; P. J. Durka, K. J. Blinowska: A unified time-frequency parametrisation of EEGs, IEEE Eng. Med. Biol. Mag. 20, 2001, 47–53; B. Shen, R. Zapulla and H. W. Harper: Matching pursuit analysis of epileptic EEG triansients, in Soc. Neurosci. Abst., 1996, vol. 22, 185; J. J. Allen, B. Shen, J. G. Williams and J. B. Angevine: Matching pursuit on ERP analysis, in Soc. Neurosci, Abst., 1995, vol. 21, 1927).Around now a resolution So-called multi-resolution procedures are obtained in time and frequency used. In addition to the wavelet analysis, which, however, poor resolution higher Frequencies, the matching pursuit is used (e.g. B. S. Mallat and Z. Zhang, Matching pursuit with time-frequency dictionaries, IEEE Transactions on Signal Processing, 41 (12), 3397-3414, 1993; G. Davis, S.G. Mallat & Z. Zhang: Adaptive time-frequency approximations with matching pursuit, in Wavelets: Theory, Algorithms and Applications, 1994; C.K. Chui, L. Montefusco and L. Puccio (Eds.), San Diego: Academic Press, 271-293; P. J. Durka: Time-frequency analyzes of EEG, PhD dissertation, Warsaw University, Warsaw, Poland, 1996; P. J. Durka and K. J. Blinowska: Analysis of EEG transients by means of matching pursuit, Ann. Biomed. Eng., Vol. 23, 1995, 608-611; Walden, A.T. & Contreras Cristan, A: Matching pursuit by undecimated discrete wavelet transform for non-stationary time series of arbitrary length. Statistics and Computing 8, 1998, 205-219; P. Vandergheynst and P. Frossard: Efficient image representation by anisotropic refinement in matching pursuit, in Proceedings of IEEE, volume 3, Salt Lake City UT, May 2001, ICASSP. 4; P. J. Durka, K. J. Blinowska: A unified time-frequency parametrization of EEGs, IEEE Eng. Med. Biol. Mag. 20, 2001, 47-53; B. Shen, R. Zapulla and H. W. Harper: Matching pursuit analysis of epileptic EEG triansients, in Soc. Neurosci. Abst., 1996, vol. 22, 185; J. J. Allen, B. Shen, J.G. Williams and J.B. Angevine: Matching pursuit on ERP analysis, in Soc. Neurosci, Abst., 1995, vol. 21, 1927).
Das
Matching Pursuit ist ein adaptiver Algorithmus. Es basiert auf einer
periodischen Entwicklung des Signals mit sogenannten (Zeit-Frequenz-)Atomen
gγ als
Ansatzfunktionen. Diese Atome werden durch die Skalierung, Verschiebung
und Modulation der Gauss-Funktion g(t) kreiert:
u – Verschiebung,
ξ – Modulation.The matching pursuit is an adaptive algorithm. It is based on a periodic development of the signal with so-called (time-frequency) atoms g γ as approach functions. These atoms are created by scaling, shifting and modulating the Gauss function g (t):
u - displacement,
ξ - modulation.
Das Wörterbuch (dictionary) ist redundant (enthält mehr Atome als zur Beschreibung eines spezifischen Signals notwendig). Daher kann das Signal wie folgt beschrieben werden: The dictionary is redundant (contains more atoms than to describe a specific one) signal necessary). Therefore, the signal can be described as follows:
Der Algorithmus sucht im Wörterbuch das Atom (Index γ0 = (s0, u0, ξ0)), welches am besten mit dem Signal korreliert. Dann kann das Signal als Summe des gewichteten Atoms und des Residuums Rf dargestellt werden: The algorithm searches for the atom in the dictionary (Index γ 0 = (s 0 , u 0 , ξ 0 )), which correlates best with the signal. Then the signal can be the sum of the weighted atom and the residual Rf are shown:
Anschließend sucht der Algorithmus das am besten korrelierte Atom für das Residuum: The algorithm then looks for the best correlated atom for the residue:
Dieser
letzte Schritt wird nun so lange wiederholt bis die Approximation
das erforderliche Niveau (geforderte Energie des zu beschreibenden
Signals oder Zahl der Iterationen (vgl.
Dabei ist die Energie des Signals wie folgt definiert: The energy of the signal is defined as follows:
Für die Zeit-Frequenz-Darstellung
des Signals (vgl.
Zusammenfassend lässt sich sagen, dass dieses Verfahren des Matching Pursuit bei sehr kurzen Signalfolgen (bis einige hundert Abtastwerte) sehr gute Ergebnisse liefert. Wird das Matching Pursuit aber auf lange Signalfolgen angewendet, wie sie in der Technik häufig vorkommen, dann ist die notwendige Rechenzeit so enorm groß, dass eine Analyse praktisch undurchführbar ist.In summary let yourself say that this method of matching pursuit at very short Signal sequences (up to a few hundred samples) have very good results supplies. If the matching pursuit is applied to long signal sequences, as they are common in technology then the necessary computing time is so enormous that an analysis practically impracticable is.
Der Erfindung liegt deshalb die Aufgabe zu Grunde, die Rechenzeit für die Auswahl der im Wörterbuch (dictionary) enthaltenen Atome zur Anwendung auf die Signalfolge wesentlich zu verkürzen, so dass auch für die Zeit-Frequenz-Analyse großer Impulsfolgen ein vertretbarer Zeitaufwand resultiert.The The invention is therefore based on the object of computing time for the selection the one in the dictionary contained atoms for application to the signal sequence significantly shorten, so that also for the Time-frequency analysis large Pulse sequences result in a reasonable amount of time.
Diese Aufgabe wird erfindungsgemäß dadurch gelöst, dass für eine Auswahl der im Wörterbuch enthaltenen Atome zur Anwendung auf die Signalfolge im Gegensatz zum bekannten Stand der Technik nicht jeweils die gesamte Impulsfolge jeweils anhand der einzelnen vorhandenen Atome auf deren Eignung und Auswahl untersucht wird, sondern dass jedes Atom im Wörterbuch jeweils nacheinander mit allen Sequenzen (trials) der gesamten Impulsfolge auf Korrelationsfähigkeit geprüft wird.This The object of the invention solved, that for a selection of those in the dictionary contained atoms for application to the signal sequence in contrast not the entire pulse sequence to the known prior art in each case based on the individual atoms present for their suitability and selection is examined, but that every atom in the dictionary one after the other with all sequences (trials) of the entire pulse sequence on correlation ability checked becomes.
Überraschend
hat sich gezeigt, dass sich im Vergleich zur bekannten und eingangs
genannten Auswahl der Atome für
die Korrelationsfähigkeit
mit dem vorgeschlagenen Verfahrensschritt, insbesondere bei größeren Impulsfolgen,
drastische Verkürzungen
der Auswahlzeit ergeben (vgl.
Ein weiterer Vorteil des erfindungsgemäßen Verfahrens liegt in seiner geringen Anforderung an den Hauptspeicher des verwendeten Rechnersystems.On Another advantage of the method according to the invention is its low demands on the main memory of the computer system used.
Die Erfindung soll nachstehend anhand eines in der Zeichnung dargestellten Ausführungsbeispiels näher erläutert werden.The Invention is intended below with reference to one shown in the drawing embodiment are explained in more detail.
Es zeigen:It demonstrate:
a)
Signal, b) Approximation mit 9 Atomen,
c) Approximation mit
14 Atomen, d) Approximation mit 20 Atomen
a) signal, b) approximation with 9 atoms,
c) Approximation with 14 atoms, d) Approximation with 20 atoms
Im
ersten Schritt werde die im dictionary enthaltenen Atome gefunden,
die am besten mit den Sequenzen der darzustellenden Signalfolge
korrelieren: fi – Sequenzen
(i = 1...N),
N – Zahl
der Sequenzen,
Rfi – Residuum für die Sequenz
fi,
– Atom,
welches am besten mit der Sequenz fi korreliert.In the first step, the atoms contained in the dictionary are found that best correlate with the sequences of the signal sequence to be displayed: f i - sequences (i = 1 ... N),
N - number of sequences,
Rf i - residue for the sequence f i ,
- atom which correlates best with the sequence f i .
Für diesen Schritt wird ein erstes Atom des Wörterbuchs (dictionary) verwendet und mit allen Sequenzen der darzustellenden Signalfolge verglichen. Der Vergleich geschieht unter Zuhilfenahme des inneren Produktes Die inneren Produkte für jede Sequenz und die Indexe γmax,i = γa werden abgespeichert. Dann wird das nächste Atom mit allen Sequenzen der Signalfolge verglichen. For this step, a first atom of the dictionary (dictionary) used and compared with all sequences of the signal sequence to be displayed. The comparison is done with the help of the inner product The inner products for each sequence and the indexes γ max, i = γ a are stored. Then the next atom compared with all sequences of the signal sequence.
Wenn
giltdann wird γmax,i = γb gesetzt.
Diese Prozedur wird für
alle Atome des Wörterbuchs
wiederholt. Daraus ergibt sich ein Set (vgl.
Anschließend werden die Residuen Rfi jeweils durch das Subtrahieren der Atome von den entsprechenden Sequenzen fi, berechnet: Then the residuals Rf i each by subtracting the atoms from the corresponding sequences f i , calculated:
In der zweiten Iteration werden, analog zum Prozess der ersten Iteration, diejenigen Atome gesucht, die am besten mit den Residuen Rfi korrelieren: In the second iteration, analogously to the process of the first iteration, those atoms are searched which best correlate with the residuals Rf i :
Allgemein können wir daher das Residuum wie folgt definieren: In general, we can therefore define the residual as follows:
Die Iterationen werden beendet, nachdem entweder eine vorgegebene Genauigkeit für das Residuum erreicht wurde oder eine maximal vorgegebene Anzahl von Iterationen erreicht wurde.The Iterations are ended after either a given accuracy for the Residual was reached or a maximum predetermined number of Iterations was achieved.
Die Erfindung ist zur Untersuchung der Korrelationsfähigkeit und Auswahl der im Wörterbuch enthaltenen Atome für die Darstellung der Impulsfolge nicht auf die Berechnung des Korrelationskoeffizienten beschränkt. Eine andere Methode ist beispielsweise die ebenfalls an sich bekannte Berechnung von RMS-Werten. Auch können mehrere, insbesondere eigenständige, Wörterbücher gleichzeitig oder separat zur Zeit-Frequenz-Analyse herangezogen werden.The invention is not limited to the calculation of the correlation coefficient for examining the correlation ability and selection of the atoms contained in the dictionary for the representation of the pulse train. Another method is, for example, the RMS calculation, which is also known per se Several, in particular independent, dictionaries can also be used simultaneously or separately for time-frequency analysis.
Claims (10)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE2003125147 DE10325147A1 (en) | 2003-05-28 | 2003-05-28 | Signal analysis method for time frequency analysis of signal sequences uses atoms in a dictionary in trials for a signal sequence |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
DE2003125147 DE10325147A1 (en) | 2003-05-28 | 2003-05-28 | Signal analysis method for time frequency analysis of signal sequences uses atoms in a dictionary in trials for a signal sequence |
Publications (1)
Publication Number | Publication Date |
---|---|
DE10325147A1 true DE10325147A1 (en) | 2004-12-16 |
Family
ID=33441567
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
DE2003125147 Withdrawn DE10325147A1 (en) | 2003-05-28 | 2003-05-28 | Signal analysis method for time frequency analysis of signal sequences uses atoms in a dictionary in trials for a signal sequence |
Country Status (1)
Country | Link |
---|---|
DE (1) | DE10325147A1 (en) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101895297A (en) * | 2010-07-30 | 2010-11-24 | 哈尔滨工业大学 | Compressed sensing-oriented block-sparse signal reconfiguring method |
CN101533045B (en) * | 2008-12-31 | 2011-04-13 | 重庆大学 | Spectral analysis method for neutron pulse sequences |
CN102279390A (en) * | 2011-05-06 | 2011-12-14 | 西南交通大学 | Intra-pulse modulation and recognition method of low signal-to-noise radar radiation source signal |
CN103675444A (en) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | High-precision time-frequency analysis method |
WO2014209439A1 (en) * | 2013-06-26 | 2014-12-31 | Thomson Licensing | System and method for predicting audience responses to content from electro-dermal activity signals |
-
2003
- 2003-05-28 DE DE2003125147 patent/DE10325147A1/en not_active Withdrawn
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101533045B (en) * | 2008-12-31 | 2011-04-13 | 重庆大学 | Spectral analysis method for neutron pulse sequences |
CN101895297A (en) * | 2010-07-30 | 2010-11-24 | 哈尔滨工业大学 | Compressed sensing-oriented block-sparse signal reconfiguring method |
CN101895297B (en) * | 2010-07-30 | 2012-12-26 | 哈尔滨工业大学 | Compressed sensing-oriented block-sparse signal reconfiguring method |
CN102279390A (en) * | 2011-05-06 | 2011-12-14 | 西南交通大学 | Intra-pulse modulation and recognition method of low signal-to-noise radar radiation source signal |
CN103675444A (en) * | 2012-08-30 | 2014-03-26 | 中国石油化工股份有限公司 | High-precision time-frequency analysis method |
CN103675444B (en) * | 2012-08-30 | 2016-09-14 | 中国石油化工股份有限公司 | A kind of high-precision Time-Frequency Analysis Method |
WO2014209439A1 (en) * | 2013-06-26 | 2014-12-31 | Thomson Licensing | System and method for predicting audience responses to content from electro-dermal activity signals |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CH716863B1 (en) | Method for building a depression detection classification model based on channel selection of multichannel electroencephalography. | |
CN107260166A (en) | A kind of electric artefact elimination method of practical online brain | |
CN106236080A (en) | Based on the removing method of myoelectricity noise in multichannel EEG signals | |
Effern et al. | Time adaptive denoising of single trial event-related potentials in the wavelet domain | |
CN114081503A (en) | Method for removing ocular artifacts in electroencephalogram signals | |
Uthayakumar et al. | Multifractal-wavelet based denoising in the classification of healthy and epileptic EEG signals | |
CN112998704A (en) | Wearable device blood oxygen saturation calculation method | |
Lamoš et al. | Spatial-temporal-spectral EEG patterns of BOLD functional network connectivity dynamics | |
Rosso et al. | Evidence of self-organization in brain electrical activity using wavelet-based informational tools | |
DE4327429C2 (en) | Method and device for brain wave analysis | |
Jervis et al. | Spectral analysis of EEG responses | |
DE10325147A1 (en) | Signal analysis method for time frequency analysis of signal sequences uses atoms in a dictionary in trials for a signal sequence | |
Principe et al. | Representing and decomposing neural potential signals | |
Pereyra et al. | Wavelet Jensen–Shannon divergence as a tool for studying the dynamics of frequency band components in EEG epileptic seizures | |
Shen et al. | Method for extracting time-varying rhythms of electroencephalography via wavelet packet analysis | |
DE10008792A1 (en) | Device for processing body signals | |
DE19608733C1 (en) | Classification of time series with sampling values for evaluating EEG or ECG | |
DE102015118226A1 (en) | Apparatus for the examination of materials by acoustic spectroscopy | |
EP1273007A1 (en) | Method for determining a characteristic data record for a data signal | |
Zhang et al. | Method of removing noise from EEG signals based on HHT method | |
Fraiwan et al. | Multiclass Epilepsy Classification Using Wavelet Decomposition, Direct Quadrature, and Shannon Entropy | |
Efe et al. | Comparison of Time-Frequency Analyzes for a Sleep Staging Application with CNN | |
Dhiman et al. | Wavelet analysis of electrical signals from brain: the electroencephalogram | |
Yadav et al. | Comparative study of various techniques for elimination of noise in EMG signal | |
Benkherrat et al. | Individual evoked potential extraction by multiresolution wavelets decomposition |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
8139 | Disposal/non-payment of the annual fee |