WO2004029825A1 - Rekonstruktion von originalsignalen aus relativmessungen - Google Patents

Rekonstruktion von originalsignalen aus relativmessungen Download PDF

Info

Publication number
WO2004029825A1
WO2004029825A1 PCT/EP2003/009161 EP0309161W WO2004029825A1 WO 2004029825 A1 WO2004029825 A1 WO 2004029825A1 EP 0309161 W EP0309161 W EP 0309161W WO 2004029825 A1 WO2004029825 A1 WO 2004029825A1
Authority
WO
WIPO (PCT)
Prior art keywords
reconstruction
unit circle
relative measurements
signals
original signals
Prior art date
Application number
PCT/EP2003/009161
Other languages
English (en)
French (fr)
Inventor
Klaus-Ulrich Wolter
Original Assignee
Db Netz Ag
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Db Netz Ag filed Critical Db Netz Ag
Priority to AU2003258632A priority Critical patent/AU2003258632A1/en
Priority to EP03798108A priority patent/EP1543439A1/de
Publication of WO2004029825A1 publication Critical patent/WO2004029825A1/de

Links

Classifications

    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B61—RAILWAYS
    • B61K—AUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
    • B61K9/00—Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
    • B61K9/08—Measuring installations for surveying permanent way
    • E—FIXED CONSTRUCTIONS
    • E01—CONSTRUCTION OF ROADS, RAILWAYS, OR BRIDGES
    • E01B—PERMANENT WAY; PERMANENT-WAY TOOLS; MACHINES FOR MAKING RAILWAYS OF ALL KINDS
    • E01B35/00—Applications of measuring apparatus or devices for track-building purposes
    • E01B35/12—Applications of measuring apparatus or devices for track-building purposes for measuring movement of the track or of the components thereof under rolling loads, e.g. depression of sleepers, increase of gauge
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/0211—Frequency selective networks using specific transformation algorithms, e.g. WALSH functions, Fermat transforms, Mersenne transforms, polynomial transforms, Hilbert transforms
    • H03H17/0213—Frequency domain filters using Fourier transforms
    • H—ELECTRICITY
    • H03—ELECTRONIC CIRCUITRY
    • H03H—IMPEDANCE NETWORKS, e.g. RESONANT CIRCUITS; RESONATORS
    • H03H17/00—Networks using digital techniques
    • H03H17/02—Frequency selective networks
    • H03H17/04—Recursive filters

Definitions

  • Patent specification DE 44 29 517 describes a method for calculating an inverse transfer function or for correcting measurement signals.
  • a measurement signal recorded by a measurement system is corrected by adding a high-frequency component of the measurement signal, which is determined by an approximate deconvolution, to a solution of a differential equation of the low-frequency component of the measurement signal, the integration constants of which are determined by the boundary conditions given by the measurement technology or by the true signal curve are.
  • this method is mathematically complex and does not represent a target-oriented solution for all tasks. Furthermore, it is an approximation method in which the original measurement signal is only approximately / approximately reconstructed.
  • Claims 2 to 6 indicate advantageous application examples of the method according to the invention.
  • the method of claim 1 is used for the reconstruction of signals according to a chord measurement method.
  • a symmetrical chord division is used in the chord measuring method. This results in a symmetrical impulse response.
  • the symmetrical impulse response has a linear phase response and therefore a constant group delay. From the knowledge of the phase response or the group delay, the phase distortion of the measuring system is advantageously also exactly compensated.
  • FIG. 3 shows an impulse response h d (z) of the measuring system of the exemplary embodiment in the locally discrete system
  • FIG. 4 shows a transfer function H d (z) of the measuring system of the exemplary embodiment in the locally discrete system
  • the measuring system of the GMTZ does not reproduce the exact physical track positions or the track position deviations. This distortion is due to a systematic error caused by the measurement method. It is imperative to exactly compensate for this systematic error in order to maintain the track position and, in particular, to assess the deviations from the track position using true-to-shape track position signals.
  • the task to be solved is to design a compensation filter h c [n] - additional LTI system - which compensates for the transmission behavior of the GMTZ.
  • 5 shows the overall system including the compensation filter with the designations: g * [n] g [n] (2.12)
  • Equation (2.17) can also be used in the form
  • H d (z) H dn ⁇ n (z) - H ap (z) (2.20)
  • H dmin (z) and H d (z) have the same amplitude distortion, but the same amplitude response but a different phase response
  • the numerator polynomial Z '(z) for H dmin (z) results from the zeros within the unit circle and the mirrored zeros. All poles and zeros of H dmin (z) are inside the unit circle.
  • H dm j n (z) and H ap (z) are shown in FIGS. 7 and 8.
  • H dm is minimal in (z)
  • H c (z) can be calculated according to equation (2.19)
  • the mechanical structure of a measuring system is considered using the wandering tendon method.
  • This measuring system has a free chord length and a free chord division.
  • the impulse response h (x) or transfer function H (j ⁇ ) or system function H (z) of the measuring system becomes so by a suitable choice of the chord length and the chord division optimizes that there is a minimal compensation filter H c (z). With this compensation filter, all zeros of H (z) are at a large distance from the unit circle or are already close to the center of the unit circle without reflection within the unit circle.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Civil Engineering (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Architecture (AREA)
  • Mathematical Analysis (AREA)
  • Structural Engineering (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mechanical Engineering (AREA)
  • Position Fixing By Use Of Radio Waves (AREA)
  • Transmission And Conversion Of Sensor Element Output (AREA)
  • Materials For Photolithography (AREA)

Abstract

Die Erfindung betrifft ein Verfahren zur Rekonstruktion von Originalsignalen aus Relativmessungen. Die Erfindung beschreibt ein verfahren, bei dem die Amplitudenverzerrung eines Messsystems exakt durch ein Kompensationsfilter kompensiert wird, wobei das Kompensationsfilter eindentig durch das Messsystem bestimmt ist. Die beschreibende Systemfunktion des Messsystems wird minimalphasig und damit invertierbar, indem alle ihre Nullstellen, die ausserhalb des Einheitskreises liegen am Einheistkreis gespiegelt werden und somit in das Innere des Einheitskreises gelangen. Vorteilshafte Anwendungsbeispiele sind die Rekonstruktion von Messungen aus einem Wandersehnenverfahren sowie zur Rekonstruktion von Gleislagemessungen des Schienenverkehrs.

Description

Rekonstruktion von Oriqinalsiqnalen aus Relativmessunqen
Beschreibung
Die Erfindung betrifft ein Verfahren zur Rekonstruktion von Originalsignalen aus Relativmessungen.
Bedingt durch das dynamische Verhalten und durch den konstruktiven Aufbau eines Messsystems ist im Allgemeinen das ausgegebene Messsignal nicht i- dentisch mit der zu messenden Größe. Bei vielen Anwendungen ist jedoch die- ser durch das Messsystem systematisch bedinge Fehler nicht tolerierbar.
Dieser Fehler kann bei der Kenntnis des dynamischen Verhaltens und/oder des konstruktiven Aufbaus des Messsystems kompensiert werden, so dass das vom Messsystem ausgebebene Messsignal exakt der zu messenden Größe entspricht. Allerdings ist dies nur bei einer bestimmten Klasse von Messsystemen möglich.
In der Patentschrift DE 44 29 517 wird ein Verfahren zur Berechnung einer in- versen Übertragungsfunktion bzw. zur Korrektur von Messsignalen beschrieben. Hierbei wird ein von einem Messsystem aufgenommenes Messsignal korrigiert durch Addition eines hochfrequenten Anteils des Messsignals, der durch eine approximative Entfaltung ermittelt wird, zu einer Lösung einer Differentialgleichung des tieffrequenten Anteils des Messsignals, wobei deren Integrationskonstanten durch die messtechnischen oder durch den wahren Signalverlauf gegebenen Randbedingungen bestimmt sind. Dieses Verfahren ist jedoch mathematisch aufwendig und stellt nicht für alle Aufgabenstellungen einen zielfüh- renden Lösungsansatz dar. Desweiteren handelt es sich um ein Näherungsverfahren, bei dem das originale Messignal nur näherungsweise/approximativ rekonstruiert wird.
Es ist daher Aufgabe der Erfindung, ein Verfahren zur Rekonstruktion von Originalsignalen aus Relativmessungen anzugeben, bei dem die geschilderten Nachteile des Standes der Technik gelöst werden.
Diese Aufgabe wird in Verbindung mit dem Oberbegriff des Hauptanspruches erfindungsgemäß durch die in Anspruch 1 angegebenen Merkmale gelöst.
Vorteil der Erfindung gegenüber dem Stand der Technik ist, dass - die Verzerrung der Amplitude des Messsignals exakt (mathematisch nachweisbar) kompensiert wird, es handelt sich somit nicht um ein Näherungsverfahren,
- das Kompensationsfilter einfach zu realisieren ist. - nur einfache Numerik (wenige Multiplikationen und Additionen) verwendet wird und daher das Verfahren auch für Echtzeitanwendungen geeignet ist.
- das Kompensationsfilter eindeutig zu bestimmen ist, es ist somit kein Optimierungsverfahren notwendig.
Ansprüche 2 bis 6 geben vorteilhafte Anwendungsbeispiele des erfindungsge- mäßen Verfahrens an.
Nach Anspruch 2 wird das Verfahren aus Anspruch 1 zur Rekonstruktion von Signalen nach einem Wandersehnenmessverfahren verwendet.
Nach Anspruch 3 wird bei dem Wandersehnenmessverfahren eine symmetrische Sehnenteilung verwendet. Hierdurch ergibt sich eine symmetrische Im- pulsantwort. Die symmetrische Impulsantwort hat einen linearen Phasengang und damit eine konstante Gruppenlaufzeit zur Folge. Aus der Kenntnis des Phasenganges bzw. der Gruppenlaufzeit wird vorteilhaft zusätzlich die Phasenverzerrung des Messssystems exakt kompensiert.
Nach Anspruch 4 wird das Verfahren aus Anspruch 1 und/oder 2 zur Rekon- struktion von Gleislagesignalen verwendet.
Die Erfindung wird nachstehend anhand zweier Ausführungsbeispiele und den Figuren 1 bis 13 näher erläutert. Die Figuren zeigen in
- Fig. 1 ein prizipielles Schaltbild mit einem originalen Signal g(x), das von einem Messsystem mit einer Impulsantwort h(x) erfasst und als Messsignal m(x) ausgegeben wird,
- Fig. 2 wie Fig. 1 , jedoch in einem ortsdiskreten System,
- Fig. 3 eine Impulsantwort hd(z) des Messsystems des Ausführungsbeispiels im ortsdiskreten System, - Fig. 4 eine Übertragungsfunktion Hd(z) des Messsystems des Ausführungsbeispiels im ortsdiskreten System,
- Fig. 5 wie Fig. 2, jedoch ergänzt um ein Kompensationsfilter mit der inver- sen Impulsantwort hc[n], so dass sich ein rekonstruiertes originales Signal g*[n] ergibt, - Fig. 6 Pole und Nullstellen der Übertragungsfunktion Hd(z) des Messsystems des Ausführungsbeispiels in der komplexen Ebene,
- Fig. 7 Pole und Nullstellen der minimalphasigen Übertragungsfunktion Hdmin(z) des Ausführungsbeispiels in der komplexen Ebene, - Fig. 8 Pole und Nullstellen der Allpass-Übertragungsfunktion Hap(z) des Ausführungsbeispiels in der komplexen Ebene,
- Fig. 9 eine minimalphasige Übertragungsfunktion Hdmin(z) des Ausführungsbeispiels nach Gleichung (2.20) im ortsdiskreten System,
- Fig. 10 eine Allpass-Übertragungsfunktion Hap(z) des Ausführungsbeispiels nach Gleichung (2.20) im ortsdiskreten System,
- Fig. 11 Pole und Nullstellen der inversen Übertragungsfunktionen Hc(z) des Kompensationsfilters des Ausführungsbeispiels in der komplexen Ebene,
- Fig. 12 eine inverse Übertragungsfunktionen H0(z) des Kompensationsfilters des Ausführungsbeispiels im ortsdiskreten System, - Fig. 13 eine kanonische Realisierung einer Übertragungsfunktion mittels digitaler Filter nach Gleichung (2.29).
In dem ersten Ausführungsbeispiel wird ein Gleismesstriebzug (GMTZ) der Deutschen Bahn AG betrachtet, der das Messverfahren des Wandersehnen- messverfahrens verwendet. Es wird beschrieben, wie eine Verzerrung der Amplitude des Messsignals des GMTZ exakt kompensiert wird.
Das Messsystem des GMTZ gibt hierbei nicht die exakten physikalischen Gleislagen bzw. die Gleislageabweichungen wieder. Diese Verzerrung beruht auf einem durch das Messverfahren bedingten systematischen Fehler. Für die Gleislageinstandhaltung und insbesondere zur Beurteilung der Gleislageabweichungen anhand formtreuer Gleislagesignale ist es zwingend erforderlich diesen systematischen Fehler exakt zu kompensieren.
Das Messsystem des GMTZ kann als lineares zeitinvariantes System (LTI- System) beschrieben werden. Die Gleislage g(x) ist hierbei die Eingangsgröße für das LTI-System. Das LTI-System wird vollständig durch die Impulsantwort h(x) bzw. Übertragungsfunktion ( ω) beschrieben. Hierbei ist H(jw) die Fou- riertransformierte der Impulsantwort h(x). m(x) ist das vom Messsystem gelieferte Signal, siehe Fig. 1 : m(x) = g(x) * h(x) - h(x-τ)dτ (2.1)
Figure imgf000006_0001
M(jω) = G(jω) • H(jω) (2.2)
Beim Messsystem GMTZ sind die Impulsantworten h(x) und die Übertragungsfunktionen H(jω) eindeutig aufgrund der Geometrie des Wandersehenmessver- fahrens gegeben. Es ist zu berücksichtigen, dass für die Messung der Längshöhen- und Richtungsabweichungen jeweils unterschiedliche Impulsantworten und Übertragungsfunktionen existieren. Des Weiteren ist auch die Messrichtung - Fahrtrichtung des GMTZ zu berücksichtigen. Somit ergeben sich für die Messung der Längshöhen- und der Richtungsabweichungen insgesamt vier unter- schiedliche Impulsantworten bzw. Übertragungsfunktionen.
Die Impulsantworten h(x) und die Übertragungsfunktionen H(jα/) des GMTZ sind durch folgende Gleichungen vollständig beschrieben:
h(x) = δ{x) — —δ(x-b) — —δ(x+a) (2.3) a + b a + b
H(jώ) = \ — Q J -~- e (2.4) a +b ^ a +b
Die Parameter a und b bezeichnen die Sehnenteilungen.
GMTZ Impulsantwort im ortsdiskreten Fall
Für die Berechnung der inversen Übertragungsfunktion wird von der z- Transformation Gebrauch gemacht, so dass bereits an dieser Stelle die kontinuierlichen Übertragungsfunktionen als ortsdiskrete Funktionen dargestellt werden.
Für den Übergang von der kontinuierlichen zur diskreten Funktion müssen folgende Parameter ersetzt werden: x -> n - Ax n : Abtastwert (2.5)
Δx : Abtastrate
Zu beachten ist hierbei, dass auch die Sehnenteilungen diskretisiert werden müssen:
— - d -— = d Ax " Δx ~ 2 dl;d2 e Z (2.6)
ie diskreten Sehnenteilungen ad und bd können direkt berechnet werden: a d = Δx (2.7) bΛ
Δx ad;bd e Z (2.8)
Die ortsdiskrete Impulsantwort des Messsystems lautet: hd[n] = δ[n] — δ[n -b] — δ[n + a] (2.9) a+b a +b
(die Diskretisierung wird zusätzlich durch eckige Klammern kenntlich gemacht). In Gleichung (2.9) bezeichnen die Parameter a und b die diskreten Sehnenteilungen.
Das Messsignal m[n] ergibt sich aus der ortsdiskreten Faltung der Gleislage g[n] mit der Impulsantwort hd[n], siehe Fig. 2:
m[n] = ∑g[k]- hd[n -k] (2.10) k=-∞
Die diskrete Impulsantwort hd[n] kann durch die diskrete Fouriertransformation (DFT) auch im Frequenzbereich dargestellt werden. In Fig. 3 ist die Impulsant- wort und in Fig. 4 die Übertragungsfunktion dargestellt.
Die Sehnenteilungen betragen in diesem Beispiel 2,6 m und 6 m, die Abtastrate wurde auf 0,2 m festgelegt:
Δx = 0,2m α = ^L = 13 (2.11)
Δx b = ^i = 30
Δx
Kompensation des Übertragungsverhaltens Die zu lösende Aufgabe besteht darin, ein Kompensationsfilter hc[n] - zusätzliches LTI-System- zu entwerfen, welches das Übertragungsverhalten des GMTZ kompensiert. In Fig. 5 ist das Gesamtsystem inklusive Kompensationsfilter dargestellt mit den Bezeichnungen: g * [n] g[n] (2.12)
m[n] = 2 Sik] - hd[n - k] (2.13)
g *[n] = m[k] - hc[n - k] (2.14)
4=-co
Berechnung des Kompensationsfilters h£fnl
Für die Berechnung des Kompensationsfilters hc[n] wird die Impulsantwort hd[n] mittels der komplexen z-Transformation transformiert. Der Vorteil dieser Heran- gehensweise liegt darin, dass wesentliche Eigenschaften von LTI-Systemen in der komplexen z-Ebene in einfacher Weise darstellbar sind. hd[ \-^→Hd (z) <-2-15)
Unter Berücksichtigung der Eigenschaften der z-Transformation und der Verwendung bekannter z-Transformationspaare kann Hd(z) direkt aus hd[n] be- stimmt werden:
^[»] = ^[«] δ[n -b] δ[n + a] (2.16) a + b a + b
Hä(z) = l ?—Z-b ~z" (2.17) a+b a+b
Durch algebraische Umformung kann Gleichung (2.17) auch in der Form
BM.i_______^_±__.m a + b N(z) dargestellt werden.
Allgemein gilt:
- Die Νullstellen von Hd(z) sind die Νullstellen von Z(z)
- Die Polstellen von Hd(z) sind die Νullstellen von Ν(z) - Damit das LTI-System Hd(z) kausal und stabil ist müssen alle Pole von Hd(z) im Inneren des Einheitskreises in der komplexen z-Ebene liegen
- Soll das System auch invertierbar sein, so müssen auch die Nullstellen von Hd(z) im Inneren des Einheitskreises in der komplexen z-Ebene liegen Liegen alle Pol- und Nullstellen von Hd(z) im Inneren des Einheitskreises in der komplexen z-Ebene nennt man Hd(z) minimalphasig.
ie zu Hd(z) inverse Systemfunktionen Hc(z) minimalphasiger Systeme kann ie folgt angegeben werden:
N(z)
Hc(z) = Hd 1 (z) = - (2.19)
Hd (z) Z(z)
In Fig. 6 sind die Pole und Nullstellen von Hd(z) entsprechend der Gleichung (2.18) in der komplexen z-Ebene eingezeichnet. Da Nullstellen von Hd(z) außerhalb des Einheitskreises liegen ist dieses System nicht minimalphasig. Da Hd(z) nicht minimalphasig ist kann Hc(z) auch nicht entsprechend der Gleichung (2.19) berechnet werden.
Zur Berechnung von Hc(z) muss Hd(z) in zwei Systemfunktionen aufgeteilt werden.
Hd(z) = Hdnάn (z) - Hap(z) (2.20)
Eigenschaften von Hdmin(z) und Hap(z):
Hdmin(z)
- Alle Pole und Nullstellen von Hdmιn(z) liegen im Inneren des Einheitskreises
- Hdmin(z) und Hd(z) haben die gleiche Amplitudenverzerrung, den glei- chen Amplitudengang jedoch unterschiedlichen Phasengang
- Hdn,in(z) ist minimalphasig
- Hap(z) ist ein Allpasssystem
- Amplitudengang konstant eins, jedoch Phasengang, Phasenverzer- rung Berechnung von Hriminfz) und Hap(z):
(a + b) -a -z-b -b- za -a -b-z(cΛb) +(a + b) -zb _ Z(z)
Hd{z) = - a + b ; = ( τa— + b~)-z-b b — - — = - Nττ(rzτ) <2-21>
- Berechnen der Nullstellen von Hd(z):
Z(z) = 0 =.> -a - b ■ z(fl+4) + (a + b) - z" = 0 (2.22)
- Berechnen der Polstellen von Hd(z):
N(Z) = 0 => (α +b) - Z4 = 0 (2.23)
- Z(z) von Hd(z) umformen:
Z(Z) = -b -(z(β+i) -~Z" + ) = -b[(z- N1)(Z-N2)(z - N3)...(z -N(a+4))] (2.24)
0 b
- Alle Νullstellen Νx von Z(z) die außerhalb des Einheitskreises liegen werden am Einheitskreis gespiegelt. Spiegelung der entsprechenden Nullstellen an ihre konjugiert reziproke Position innerhalb des Einheitskreises.
NJ > 1 => Z -» — (2.25)
Z
- Aus den Νullstellen innerhalb des Einheitskreises und den gespiegelten Νullstellen ergibt sich das Zählerpolynom Z'(z) für Hdmin(z). Alle Pole und Νullstellen von Hdmin(z) liegen im Inneren des Einheitskreises.
Hdmin(z) ist minimalphasig.
Figure imgf000010_0001
- Hap(z) besteht aus allen Νullstellen von Hd(z) die außerhalb des Einheitskreises liegen zusammen mit den Polen die die gespiegelten konjugiert re- ziproken Νullstellen von Hdmin(z) wieder aufheben.
In Fig. 7 und Fig. 8 sind die Pole und Νullstellen von Hdmjn(z) und Hap(z) eingezeichnet.
Fig. 9 und Fig. 10 stellen die Übertragungsfunktionen von Hdmin(jω) und Hap(jω) dar. Berechnung von Hs(z):
Da Hdmin(z) minimalphasig ist kann Hc(z) entsprechend der Gleichung (2.19) berechnet werden
HC(Z) = H ! ( ) = ---— = !^
Hdvb_ {z Z z) (2 27)
In Fig. 11 und Fig. 12 sind die Pole und die Nullstellen sowie die Übertragungsfunktionen von Ηo(jω) dargestellt.
Übertragungsverhalten H^z) des gesamten Systems:
Das Übertragungsverhalten des gesamten Systems Hges(z) lautet:
Hges(z) = Hd(z) - Hc(z) = Hd(z) - l __ Hd(z) -^ß- ^ Hap z) (2.28)
Der Amplitudengang von Hd(z) wird exakt kompensiert. Der Phasengang des gesamten Systems ist gleich dem Phasengang von Hap(z).
Bestimmung der Filterkoeffizienten a, b aus Hc(z)
Zur Realisierung der Übertragungsfunktion mittels digitaler Filter bietet sich die Realisierung mit einer kanonischen Grundstruktur an. Der Vorteil besteht darin, dass die Filterkoeffizienten a und b direkt aus der Systemfunktion abgelesen werden können. Gegebenenfalls ist ein linearer Faktor zu berücksichtigen.
Allgemeine Form, siehe Fig. 13:
rr, . h +b2 -z~l +b3 -z~2 +b4 -z~3 + --- + +bn,, -z'n
H(z) = — ! 2— 3 ; * : *** — (2.29) al +a2 -z +az -z +a - z +--- + + m - z
Das Übertragungsverhalten des Kompensationsfilter Hc(z) kann entsprechend der Gleichung (2.29) dargestellt werden. Somit können auch für Hc(z) die Filterkoeffizienten a und b direkt abgelesen werden.
Im zweiten Ausführungsbeispiel wird der mechanische Aufbau eines Messsystems nach dem Wandersehnenverfahren betrachtet. Dieses Messystem weist eine freie Sehnenlänge und eine freie Sehnenteilung auf. Durch eine geeignete Wahl der Sehnenlänge und der Sehnenteilung wird die Impulsantwort h(x) bzw. Übertragungsfunktion H(jω) bzw. Systemfunktion H(z) des Messsystems so optimiert, dass sich ein minimales Kompensationsfilter Hc(z) ergibt. Bei diesem Kompensationsfilter weisen alle Nullstellen von H(z) einen großen Abstand vom Einheitskreis auf oder befinden sich bereits ohne Spiegelung innerhalb des Einheitskreises nahe am Mittelpunkt des Einheitskreises.
Hierdurch wird erreicht, dass sich das Kompensationsverfahren wesentlich vereinfacht bzw. eine optimale Rekonstruktion erhalten wird.
Verzeichnis der verwendeten Formelzeichen
m(x) Messsignal h(x) Impulsantwort des Messsystems g(x) originales Signal g*(x) rekonstruiertes Signal m[n] Messsignal im diskreten System h[n] Impulsantwort des Messsystems im diskreten System g[n] originales Signal im diskreten System g*[n] rekonstruiertes Signal im diskreten System M(jω) Fouriertransformierte von m(x)
H(jω) Fouriertransformierte von h(x), Übertragungsfunktion
G(jω) Fouriertransformierte von g(x)
M(z) z-Transformierte von m(x)
H(z) z-Transformierte von h(x), Systemfunktion G(z) z-Transformierte von g(x) x Ort j imaginäre Einheit ω Kreisfrequenz δ(x) Dirac-Funktion a, b Parameter der Sehnenteilung bzw. Filterkoeffizienten n Abtastwert
Δx Abtastrate z Ort der komplexen z-Ebene, komplexe Zahl z* konjugiert komplexe von z d Index des ortsdiskreten Systems
Z(z) Zählerpolynom
N(z) Nennerpolynom
Hc(z) Systemfunktion des Kompensationsfilters Hmin(z) minimalphasige Systemfunktion zu Hd(z)
Hap(z) Allpass-Systemfunktion
Nx Nullstelle
Hges(z) Systemfunktion des gesamten Systems

Claims

Patentansprüche
1. Verfahren zur Rekonstruktion von Originalsignalen aus Relativmessungen, bei dem ein Messsignal m(x) mit seiner Fourriertransformierten M(jω) durch ein Messsystem mit einer bekannten und eindeutigen Impulsantwort h(x) bzw. Übertragungsfunktion H(jω) bzw. Systemfunktion H(z) aus einem originalen Signal g(x) mit seiner Fourriertransformierten G(jω) ermittelt wird mit dem Zusammenhang M(jω) = G(jω) -H(jω), dadurch gekennzeichnet, dass
- die Impulsantwort mit der komplexen z-Transformation transformiert wird h(x) → H(z), - H(z) in zwei Systemfunktionen aufgeteilt wird H(z) = Hmin(z) -Hap(z), wobei Hmin(z) ein minimalphasiges System und Hap(z) ein Allpasssystem beschreibt,
- ein Kompensationsfilter Hc(z) = 1/ Hmin(z) verwendet wird, durch das die Amplitudenverzerrung bzw. der Amplitudengang des Messsystems exakt kompensiert wird, wobei Hmin(z) berechnet wird indem
- Nullstellen von H(z), die außerhalb des Einheitskreises liegen am Ein- heistkreis gespiegelt werden und somit in das Innere des Einheitskreises gelangen, indem die entsprechende Nullstelle mit z → 1/z* an ihre kon- jugiert reziproke Position innerhalb des Einheitskreises gespiegelt wird, wobei z* konjugiert komplex zu z ist,
- trotz dieser Spiegelung der Amplitudengang von H(z) gleich dem Amplitudengang von Hmin(z) gelassen wird,
- die neue Systemfunktion Hmin(z) dadurch minimalphasig wird und somit Hc(z) berechnet wird.
2. Verfahren zur Rekonstruktion von Originalsignalen aus Relativmessungen nach Anspruch 1 , dadurch gekennzeichnet, dass das Verfahren zur Rekonstruktion von Signalen nach einem Wandersehnenmessverfahren ver- wendet wird.
3. Verfahren zur Rekonstruktion von Originalsignalen aus Relativmessungen nach Anspruch 2, dadurch gekennzeichnet, dass eine symmetrische Sehnenteilung verwendet wird. Verfahren zur Rekonstruktion von Originalsignalen aus Relativmessungen nach einem oder mehreren der Ansprüche 1 bis 3, dadurch gekennzeichnet, dass das Verfahren zur Rekonstruktion von Gleislagesignalen verwendet wird.
PCT/EP2003/009161 2002-09-24 2003-08-19 Rekonstruktion von originalsignalen aus relativmessungen WO2004029825A1 (de)

Priority Applications (2)

Application Number Priority Date Filing Date Title
AU2003258632A AU2003258632A1 (en) 2002-09-24 2003-08-19 Reconstruction of original signals from relative measurements
EP03798108A EP1543439A1 (de) 2002-09-24 2003-08-19 Rekonstruktion von originalsignalen aus relativmessungen

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE10244357.2 2002-09-24
DE10244357 2002-09-24

Publications (1)

Publication Number Publication Date
WO2004029825A1 true WO2004029825A1 (de) 2004-04-08

Family

ID=31969489

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2003/009161 WO2004029825A1 (de) 2002-09-24 2003-08-19 Rekonstruktion von originalsignalen aus relativmessungen

Country Status (4)

Country Link
EP (1) EP1543439A1 (de)
AU (1) AU2003258632A1 (de)
DE (1) DE10337976B4 (de)
WO (1) WO2004029825A1 (de)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3005321A1 (fr) * 2013-05-05 2014-11-07 Leyfa Measurement Dispositif de mesure de la geometrie d'une voie ferree et procede d'estimation des profils de nivellement et de dressage de ladite voie ferree

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE102009031819B4 (de) * 2009-07-03 2016-05-04 Deutsche Bahn Ag Verfahren zur Ermittlung von kurzwelliger Gleislagegeometrie und der Schieneneinsenkungen unter Last
DE102011101226A1 (de) 2011-05-11 2012-11-15 Deutsche Bahn Ag Verfahren zur Beschreibung von Gleislageabweichungen
EP2957674B1 (de) 2014-06-18 2017-10-11 HP3 Real GmbH Verfahren zum Betreiben einer auf einer Gleisanlage verfahrbaren Oberbaumaschine
EP2960371B1 (de) 2014-06-27 2017-08-09 HP3 Real GmbH Vorrichtung zum Vermessen von Gleisen
CN104501755A (zh) * 2014-12-30 2015-04-08 苏州路云机电设备有限公司 一种便携式尖轨心轨测量仪
CN106553086B (zh) * 2016-10-27 2018-06-08 大连理工大学 快速高精度的复杂曲面制孔点法矢量测量方法
AT520795B1 (de) 2017-12-21 2020-03-15 Plasser & Theurer Export Von Bahnbaumaschinen Gmbh Gleisbaumaschine und Verfahren zum Nivellieren eines Gleises
AT520894B1 (de) 2018-01-22 2021-01-15 Hp3 Real Gmbh Verfahren zur Gleislageverbesserung durch eine gleisfahrbare Gleisstopfmaschine
AT523717B1 (de) 2020-06-18 2021-11-15 Hp3 Real Gmbh Verfahren zum Vermessen einer Gleislage
AT523627B1 (de) 2020-09-16 2021-10-15 Plasser & Theurer Export Von Bahnbaumaschinen Gmbh Verfahren und System zur Ermittlung eines Soll-Gleisverlaufs für eine Lagekorrektur
AT524435B1 (de) 2020-11-25 2022-06-15 Plasser & Theurer Export Von Bahnbaumaschinen Gmbh Verfahren und System zur Ermittlung von Korrekturwerten für eine Lagekorrektur eines Gleises

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4429517A1 (de) * 1994-08-19 1996-02-22 Man Technologie Gmbh Verfahren zur Korrektur einer Meßkurve oder eines Signalverlaufs, eine zugehörige Vorrichtung und dessen bzw. deren Anwendung zur Rekonstruktion von Lagefehlern bei Bahngleisen aus geometrischen Relativmessungen

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE4429517A1 (de) * 1994-08-19 1996-02-22 Man Technologie Gmbh Verfahren zur Korrektur einer Meßkurve oder eines Signalverlaufs, eine zugehörige Vorrichtung und dessen bzw. deren Anwendung zur Rekonstruktion von Lagefehlern bei Bahngleisen aus geometrischen Relativmessungen

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
ESVELD C: "DAS MESSEN UND KORRIGIEREN DER GLEISGEOMETRIE", ETR EISENBAHNTECHNISCHE RUNDSCHAU, HESTRA-VERLAG. DARMSTADT, DE, vol. 29, no. 5, 1980, pages 351 - 354,356,35, XP002061334 *
OPPENHEIM A V; SCHAFER R W: "Zeitdiskrete Signalverarbeitung", 1999, R. OLDENBOURG VERLAG, MÜNCHEN, XP002262235 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR3005321A1 (fr) * 2013-05-05 2014-11-07 Leyfa Measurement Dispositif de mesure de la geometrie d'une voie ferree et procede d'estimation des profils de nivellement et de dressage de ladite voie ferree
EP2806065A1 (de) 2013-05-05 2014-11-26 Leyfa Measurement Schätzverfahren der Profile zum Ausgleichen und Abrichten besagter Eisenbahnstrecke, und entsprechende Messvorrichtung
EP3162958A1 (de) * 2013-05-05 2017-05-03 Leyfa Measurement Mess- und schätzverfahren der geometrie einer eisenbahnschiene

Also Published As

Publication number Publication date
AU2003258632A1 (en) 2004-04-19
DE10337976A1 (de) 2004-04-01
EP1543439A1 (de) 2005-06-22
DE10337976B4 (de) 2007-12-27

Similar Documents

Publication Publication Date Title
EP1543439A1 (de) Rekonstruktion von originalsignalen aus relativmessungen
EP2696924B1 (de) Vorrichtung und verfahren zur datenverarbeitung physiologischer signale
DE102012214351B4 (de) Angiographisches Untersuchungsverfahren eines Gefäßsystems
DE69627755T2 (de) Verfahren und Vorrichtung zur temporären Rauschfilterung einer Bildfolge
WO2001086819A9 (de) Verfahren und vorrichtung zur kompression und/oder dekompression sowie zur analyse und darstellung von daten
DE2436966C3 (de) Bandfilter
EP0574465A1 (de) Anordnung zum herausfiltern von grundlinienschwankungen aus physiologischen messsignalen.
DE2059728B2 (de) Spulenloser daempfungs- und laufzeitentzerrer
Jäger et al. Psychiatrische Versorgungsepidemiologie.
DE1762157B2 (de) Vorrichtung zum umsetzen von analogen signalen, die einer monoton zunehmenden funktion entsprechen, in digitale werte
DE1562215B2 (de) Variabler entzerrer
DE3709532A1 (de) Verfahren zur pruefung von anordnungen
DE928967C (de) Leiternetzwerk mit einer innerhalb des Durchlassbereichs geebneten Betriebsdaempfungskurve
DE4429517C2 (de) Vorrichtung und Verfahren zur Korrektur einer Meßkurve oder eines Signalverlaufs und deren bzw. dessen Anwendung zur Rekonstruktion von Lagefehlern bei Bahngleisen aus geometrischen Relativmessungen
DE69401727T2 (de) Vorrichtung zur Beseitigung von Nachschwingungen in gefilterten EKG-Signalen
DE2421033A1 (de) Bandpassfilterschaltung
DE617116C (de) Vierpolige Siebschaltung
Altenhöner et al. Occupational health management interventions: evaluation of a seminar series for shift workers
DE4327275C2 (de) Verfahren zur digitalen Verarbeitung von schmalbandigen Signalen mit variabler Mittenfrequenz
DE102009016013B4 (de) Verfahren zur Aufhärtungskorrektur
DE452495C (de) Kettenleiter
DE2806921C2 (de) Filterschaltung mit einer Übertragungsfunktion dritten Grades
DE2052842C3 (de) Schmalbandige schnelleinschwingende Filteranordnung
DE562317C (de) Drehkondensator, insbesondere fuer Messzwecke
DE3150238C2 (de) Mit kontinuierlicher Drehung abtastender Fächerstrahl-Tomograph

Legal Events

Date Code Title Description
AK Designated states

Kind code of ref document: A1

Designated state(s): AE AG AL AM AT AU AZ BA BB BG BR BY BZ CA CH CN CO CR CU CZ DE DK DM DZ EC EE ES FI GB GD GE GH GM HR HU ID IL IN IS JP KE KG KP KR KZ LC LK LR LS LT LU LV MA MD MG MK MN MW MX MZ NI NO NZ OM PG PH PL PT RO RU SC SD SE SG SK SL SY TJ TM TN TR TT TZ UA UG US UZ VC VN YU ZA ZM ZW

AL Designated countries for regional patents

Kind code of ref document: A1

Designated state(s): GH GM KE LS MW MZ SD SL SZ TZ UG ZM ZW AM AZ BY KG KZ MD RU TJ TM AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LU MC NL PT RO SE SI SK TR BF BJ CF CG CI CM GA GN GQ GW ML MR NE SN TD TG

121 Ep: the epo has been informed by wipo that ep was designated in this application
WWE Wipo information: entry into national phase

Ref document number: 2003798108

Country of ref document: EP

WWP Wipo information: published in national office

Ref document number: 2003798108

Country of ref document: EP

NENP Non-entry into the national phase

Ref country code: JP

WWW Wipo information: withdrawn in national office

Country of ref document: JP

WWR Wipo information: refused in national office

Ref document number: 2003798108

Country of ref document: EP

WWW Wipo information: withdrawn in national office

Ref document number: 2003798108

Country of ref document: EP