A T L A S E L E K T R O N I K G m b H

Bremen

VERFAHREN ZUM GENERIEREN VON ZIELMARKIERENDEN PEILTRACKS

Die Erfindung betrifft ein Verfahren zum Generieren von zielmarkierenden Peiltracks für mehrere schallabstrahlende Ziele, die in einem Seegebiet von mehreren, hydroakustische Empfangsantennen aufweisenden Sensoren gepeilt werden, nach dem Oberbegriff des Anspruchs 1.

In der Sonartechnik werden zur passiven Zielaufklärung eines Seegebiets von einem Trägerfahrzeug, z.B. U-Boot oder Oberflächenschiff, üblicherweise verschiedenartige Sensoren eingesetzt, die die Ziele über deren Schallabstrahlung ins Wasser auffassen und aufgefasste Ziele bis Zielverlust verfolgen und mehr oder weniger exakte Zieldaten liefern. (L. Kühnle „Classification and Identification - CAI - by Submarine Sonars", Naval Forces, Nov. 1987, S. 27 - 31) . Unter einem Sensor wird dabei eine hydroakustische Empfangsantenne, auch Unterwasserantenne genannt, wie z.B. eine Zylinderbasis (Cylindrical Array Sonar CAS) , eine Seitenantenne (Flank Array Sonar FAS), eine Schleppantenne (Towed Array Sonar TAS) , sowie deren nachgeschaltete Signalverarbeitungsalgorithmen, wie Breitband-Detektion (Broadband Detection BDT) Demon und Lofar, verstanden. Sensoren im hier gebrauchten Sinne sind somit CAS BDT, FAS BDT, die Peilungen zu Zielen liefern, CAS Demon, FAS Demon, TAS Demon, FAS Lofar und TAS Lofar, die Frequenzlinien der Ziele liefern. Von den eingesetzten Sensoren erfassen einige

Sensoren alle Ziele und einige Sensoren nur ein Teil der Ziele, z.B. wegen des durch die zugehörige Konfiguartion der Empfangsantenne vorgegebenen Frequenzbereichs und den damit einhergehenden unterschiedlichen Reichweiten oder wegen unterschiedlichen, tiefenabhängigen

Schallausbreitungsbedingungen. Die von den Sensoren CAS BDT, FAS BDT, TAS BDT gelieferten Peilungen werden üblicherweise in zeitlicher Reihenfolge für jedes Ziel und jeden Sensor dargestellt und dem Sonarbediener in einer sog. Wasserfalldarstellung angeboten, bei der die Zeitachse, denen die Peilungen zugeordnet sind, vertikal verläuft und am oberen Bildrand immer die zeitlich jüngsten Peilungen angezeigt sind. Bei einer solchen Aufbereitung und Darstellung der Sensordaten als von jedem Sensor für jedes von ihm gepeilte Ziel gelieferte Peilungsspuren, nachfolgend Sensortracks genannt, muss der Bediener fortlaufend eine Vielzahl von Sensortracks verfolgen. Insbesondere bei sich kreuzenden Sensortracks oder zeitweise unterbrochenen Sensortracks ist es für den Bediener schwierig bis nahezu unmöglich, zu entscheiden, ob der nach Unterbrechung oder Kreuzung sich fortsetzende Sensortrack das gleiche Ziel verfolgt oder der Sensor sich auf ein anderes Ziel gesetzt hat.

Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zur Generierung von zielmarkierenden Peiltracks oder zielmarkierenden Peilungsspuren anzugeben, das nur so viele zielmarkierende Peiltracks generiert, wie Ziele von der Gesamtheit der Sensoren aufgefasst werden, und das Trackstücke, die einunddasselbe Ziel markieren, zielrichtig einander zuordnet.

Die Aufgabe ist erfindungsgemäß durch die Merkmale des Anspruchs 1 gelöst.

Das erfindungsgemäße Verfahren hat den Vorteil, dass durch die Fusion der Sensortracks zu Kontakttracks eine Reduktion der Peilungsspuren, oder Peiltracks auf nur solche Peiltracks erreicht wird, die jeweils von einunddemselben detektierten Ziel herrühren - unabhängig davon von welchen Sensoren dieses Ziel verfolgt wird. Da damit die Zahl der Peilungsspuren deutlich abnimmt, wird der Sonarbediener bei der Zielϋberwachung, Zielverfolgung und Zielauswahl sehr stark entlastet. Das vom Bediener beobachtete Display, auf dem ausschließlich Kontakttracks und Zieltracks, vorzugsweise in unterschiedlicher Farbgebung, dargestellt werden, wird für den Bediener überschaubarer, und einzelne Ziele lassen sich auch nach im Minutenbereich liegenden, längeren Beobachtungspausen durch den Bediener leicht wiederfinden.

Darüber hinaus bietet das erfindungsgemäße Verfahren erhebliche Vorteile für die Lageerarbeitung in einem Seegebiet mittels der sog. TMA-Funktion (Target Motion Analysis), bei welcher aus passiv ermittelten Peilungen Positionen von Zielen sowie deren Bewegungsrichtung und Bewegungsgeschwindigkeit als Zieldaten geschätzt werden. Da das erfindungsgemäße Verfahren eine zuverlässige Zuordnung von endenden und zeitlich später wiederauftauchenden Kontakttracks, die Peilungen von einunddemselben Ziel repräsentieren, vornimmt, kann die TMA-Funktion sowohl auf Peilungen aus den endenden Kontakttracks als auch auf Peilungen aus den wiederauftauchenden Kontakttracks gestützt werden. Insbesondere kann zum Schätzen der Zieldaten auf Peilungen zum gleichen Ziel zurückgegriffen werden, die vor einem Eigenmanöver des Trägerfahrzeugs liegen, so dass

Peilungen von zwei vom Trägerfahrzeug durchfahrenen Wegstrecken zur Verfügung stehen. Durch den Zugriff auf Peilungen des gleichen Ziels aus der auch weit in der Vergangenheit zurückliegenden Messzeit können in sehr viel kürzerer Zeit zuverlässige Schätzungen der Zielposition, Zielgeschwindigkeit und Fahrtrichtung erreicht werden, da die TMA-Funktion dann zusammen mit den aktuellen Peilungen sehr viel schneller konvergiert. Zur TMA-Funktion wird auf die DE 101 29 726 Al und DE 103 52 738 Al verwiesen.

Durch das Zuordnen bzw. Zusammenführen zugehöriger Kontakttracks und das überbrücken fehlender Peilungen über längere Zeiträume sowie durch die richtige Zuordnung zu den nachfolgenden Peilungen wird es möglich, eine Zuordnung der Zielparameter für die Klassifizierung der Ziele auch aus der Vergangenheit heraus vorzunehmen und auf das Fahrverhalten der Ziele, z.B. Manöverfahrten oder Geschwindigkeits- und/oder Kursänderungen, zu schließen.

Zweckmäßige Ausführungsformen des erfindungsgemäßen Verfahrens mit vorteilhaften Weiterbildungen und Ausgestaltungen der Erfindung ergeben sich aus den weiteren Ansprüchen .

Die Erfindung ist anhand eines in der Zeichnung illustrierten Ausführungsbeispiels im folgenden näher beschrieben. Es zeigen :

Fig. 1 ein Blockdiagramm einer sequenziellen Track- to-Track-Fusion mit beispielhafter Darstellung von daraus resultierenden Kontakttracks als Funktionen der Peilwinkel B über der Zeit t,

Fig. 2 eine Kostenmatrix für die Track-to-Track-

Fusion von N _A Sensortracks eines Sensors A und N _B Sensortracks eines Sensors B,

Fig. 3 ein Zahlenbeispiel der Kostenmatrix in Fig. 2 für drei Sensortracks des Sensors A und drei Sensortracks des Sensors B,

Fig. 4 eine schematisierte, graphische Darstellung der jeweils drei Sensortracks der Sensoren A und B und der aus der Track-to-Track-Fusion dieser Sensortracks über die Kostenmatrix gemäß Fig. 3 zugeordneten und zu Kontakttracks (strichliniert dargestellt) fusionierten Sensortracks,

Fig. 5 eine Kostenmatrix für die Prüfung der

Zugehörigkeit von vier neuen Kontakttracks NT und drei alten Kontakttracks OT, wie sie im Display in Fig. 1 schematisiert dargestellt sind.

Mit dem in den Zeichnungen illustrierten Verfahren werden die fortlaufenden Peilungen (Peilwinkel B) von mehreren Zielen durch mehrere Sensoren, die als sog. Sensorspuren oder Sensortracks bezeichnet werden, so verarbeitet, dass nur noch Spuren oder Tracks vorhanden sind, die - unabhängig davon, welcher Sensor sie liefert - jeweils nur noch Peilungen einunddesselben Ziels enthalten. Dabei können die Sensortracks stückweise von verschiedenen Zielen stammen oder stückweise nur aus Rauschen bestehen und auch zeitlich unterbrochen sein. Hierzu werden in einem ersten Verfahrensschritt die Sensortracks der mehreren Sensoren 10,

die in Fig. 1 mit A, B, C ... Y gekennzeichnet sind, ausgehend von dem genauesten Sensor, z.B. Sensor A, in einer sequenziellen Spur-zu-Spur- oder Track-to-Track-Fusion zu Kontakttracks fusioniert. Diese Kontakttracks sind zielspezifisch, d. h. jeder Kontakttrack enthält nur Peilungen in Form von Winkelangaben gegen eine Bezugsrichtung zu einunddemselben Ziel. Da die Sensortracks stückweise unterbrochen sein können, z.B. durch Zielverlust und erneutes Auffassen des alten oder eines neuen Ziels, werden damit auch Kontakttracks entstehen, die zu einem Zeitpunkt enden, und Kontakttracks entstehen, die zu einem anderen Zeitpunkt beginnen. Solche zeitlich versetzten Kontakttracks werden in einem zweiten Verfahrensschritt paarweise daraufhin überprüft, ob sie zusammengehören, d.h. ob ein neu beginnender Kontakttrack die Fortsetzung eines zuvor geendeten Kontakttrack ist oder ein neu aufgefasstes Ziel repräsentiert. Ergibt die Prüfung ein hohes Maß an Wahrscheinlichkeit, dass die Kontakttracks zueinander gehören, so wird in einem dritten Verfahrensschritt die Zuordnung durchgeführt und die einander zugeordneten Kontakttracks als Zieltracks ausgewiesen. Bevorzugt werden dabei in der graphischen Darstellung Ende und Anfang zweier einander zugeordneter Kontakttracks durch eine Linie miteinander zu einem Zieltrack verbunden, so dass unterbrechungslose Zieltracks vorhanden sind, die dem Sonarbediener ein leichtes Verfolgen der einzelnen Zieltracks auch nach deren Kreuzung möglich machen. Solche Verbindungslinien sind im Display der Fig. 1 strichliniert symbolisiert. Zieltracks sind in Fig. 1 NT ₁ /OT3, NT ₂ /OT ₂ , NT ₃ ZOT ₁ und NT ₄ .

In Fig. 1 bis 4 ist beispielhaft eine sequenzielle Track-to- Track-Fusion illustriert, die als Ergebnis die Kontakttracks

liefert. Jeder der Sensoren 10 in Fig. 1 peilt ein oder mehrere Ziele im Seegebiet und liefert in zeitlicher Reihenfolge entsprechende Peilwinkel B, Peilungen genannt, die zeitlich aneinandergereiht Peilungsspuren der Sensoren 10, kurz Sensortracks genannt, ergeben. Ausgehend von dem genauesten Sensor, beispielsweise des Sensors A, der z.B. ein CAS BDT-Sensor ist, werden alle Sensortracks des Sensors A in eine zentrale Trackdatei 11 eingeschrieben. Die in der Trackdatei 11 enthaltenen Sensortracks werden im folgenden als zentrale Tracks bezeichnet. Zunächst werden die Sensortracks eines zweiten Sensors, z.B. des Sensors B, mit den in der Trackdatei 11 enthaltenen zentralen Tracks auf Zugehörigkeit oder Zuordnung zu einem gemeinsamen Ziel geprüft und bei einer ausreichend hohen Wahrscheinlichkeit an Zugehörigkeit dem jeweiligen zentralen Track zugeordnet und mit diesem fusioniert. Dies erfolgt im Block 12 „Zuordnung und Fusion". In der Trackdatei 11 werden die zentralen Tracks, aus denen die fusionierten Tracks hervorgegangen sind, mit den fusionierten Tracks überschrieben. Zentrale Tracks, denen keiner der Sensortracks des Sensors B zugeordnet worden sind, bleiben in der zentralen Trackdatei 11 erhalten. Sensortracks des Sensors B, die zu keiner Fusion mit einem der zentralen Tracks geführt haben, werden zusätzlich in die Trackdatei 11 eingeschrieben. Die Trackdatei 11 ist damit aktualisiert. Danach werden die Sensortracks des Sensors C mit den insgesamt in der aktualisierten Trackdatei 11 enthaltenen zentralen Tracks auf Zugehörigkeit oder auf Zuordnung zu einem gemeinsamen Ziel geprüft. Bei ausreichend hoher Wahrscheinlichkeit der Zugehörigkeit werden die Sensortracks des Sensors C den entsprechenden zentralen Tracks zugeordnet und mit diesen fusioniert. Mit den fusionierten Tracks werden in der Trackdatei 11 die zentralen Tracks, aus denen die

fusionierten Tracks hervorgegangen sind, überschrieben. Nicht fusionierte Sensortracks des Sensors C werden zusätzlich in die Trackdatei 11 eingeschrieben. Diese Vorgänge wiederholen sich mit allen Sensortracks der nächsten Sensoren bis hin zu den Sensortracks des Sensors Y solange, bis alle Sensortracks als zentrale Tracks in der aktualisierten Trackdatei 11 enthalten sind. Die in der zuletzt aktualisierten Trackdatei 11 enthaltenen zentralen Tracks werden nunmehr als Kontakttracks ausgegeben und dem Sonarbediener in einem Display 13 graphisch angezeigt, wobei je ein Kontakttrack ein Ziel markiert und als Funktion des Peilwinkels B über der Zeit t dargestellt wird. Die in Fig. 1 eingezeichneten Umschalter 14, 15, 16 symbolisieren lediglich die Reihenfolge der beschriebenen Verfahrensschritte, so dass für das Einschreiben der Sensortracks des Sensors A in die Trackdatei 11 der Umschalter 14 geschlossen und danach wieder geöffnet wird, der Umschalter 15 für das Auslesen der zentralen Tracks aus der Trackdatei 11 zwecks Prüfung mit den Sensortracks der übrigen Sensoren 10 mit dem Block 12 verbunden ist und für das Auslesen der zuletzt aktualisierten Trackdatei 11 auf das Display 13 umgeschaltet wird und der Umschalter 16 für das Zuführen der Sensortracks der Sensoren B, C ... Y zum Block 12 nacheinander den Block 12 mit den jeweiligen Sensoren B, C, ... Y verbindet. Block 17 ist für die entsprechende Steuerung der Umschalter 14, 15, 16 zuständig und symbolisiert die Einleitung der entsprechenden Verfahrensschritte.

Die Prüfung der Zuordnung der zentralen Tracks und der von einem Sensor 10 angelieferten Sensortracks erfolgt nach der Global-Nearest-Neighbor-Methode (GNN-Methode) . Das Ziel dieser GNN-Methode ist es, in jeder Zeitscheibe mit neuen Spur- oder Trackdaten die wahrscheinlichste Zuordnung der Spuren oder Tracks zweier Sensoren zu bestimmen. Hierzu wird

mit den i zentralen Tracks der aktualisierten, zentralen Trackdatei 11 und den j Sensortracks eines beliebigen Sensors 10 eine zweiteilige Kostenmatrix mit einem Zuordnungsteil und einem Nichtzuordnungsteil aufgestellt. In den Zuordnungsteil werden die Kosten der wechselseitigen Zuordnung der i zentralen Tracks und der j Sensortracks des entsprechenden Sensors 10 und in den Nichtzuordnungsteil die Kosten der NichtZuordnung der zentralen Tracks eingetragen. In Fig. 2 ist ein Beispiel für die Kostenmatrix dargestellt, in welchem angenommen ist, dass in der zentralen Trackdatei 11 ausschließlich die Sensortracks 1 bis N _A des Sensors A enthalten sind, die als zentrale Tracks im Block 12 „Zuordnung und Fusion" mit den Sensortracks 1 bis N _B des Sensors B geprüft und ggf. fusioniert werden.

Die Zuordnungskosten C ₁ ,-,, in Fig. 2 Cn, C ₁₂ ••• C _{N N} , werden zum Zeitpunkt t _k als Summe der Quadrate der normierten statistischen Abstande der Peilungen und der Summe der normierten statistischen Abstande der Peilraten berechnet, wobei wahlweise auf eine Hinzuziehung der Abstande der

Peilraten verzichtet werden kann. Ist θ die Peilung, θ die Peilrate, und σ die jeweilige Standardabweichung, so ergibt sich für die Zuordnungskosten C ₁ , _D

K ((S) _ (α ) ² K (ά) _ 0 Y

Wo - Z _{J ^} 2 , 2 + α ^• £_, 2 , _2 ^{[ 1} > ' n -k- bl + 1 ^σ θ _{A l} , _n + ^σ θ _{B j n} n -k-bl + 1 \ _{i n} + ^σ _0B , __

wobei die Summation über ein gleitendes Zeitfenster der Lange bl erfolgt. Wird α=0 gewählt, so wird auf die Einbeziehung der Peilraten in die Kostenrechnung verzichtet. Wird α=l gewählt, so werden die Peilraten in die Kostenrechnung mit einbezogen. Die Kosten für die NichtZuordnung aller i

zentralen Tracks (im Beispiel die Sensortracks des Sensors A) zu den Sensortracks des Sensors B werden für alle in der Trackdatei 11 vorhandenen i zentralen Tracks mit einem gleichen Wert C _N vorgegeben. Die in dem Nichtzuordnungsteil enthaltenen Kombinationsmöglichkeiten der Sensortracks, die verboten bzw. nicht möglich sind, sind mit einem X gekennzeichnet, das einen sehr großen Wert darstellt.

In Fig. 3 ist ein Zahlenbeispiel für die nach Gl.(l) berechneten Kosten bei der Prüfung von drei Sensortracks des Sensors A mit drei Sensortracks des Sensors B, wie sie graphisch in Fig. 4 dargestellt sind, angegeben. Die Kosten für die NichtZuordnung der drei Sensortracks des Sensors A zu jedem der Sensortracks des Sensors B, der sog. Gate-Wert, ist mit „9" vorgegeben. Diesem Gate-Wert liegt die überlegung zugrunde, dass die Wahrscheinlichkeit, dass zwei zusammengehörige Spuren oder Sensortracks einen Kostenwert verursachen, der größer „9" ist, nur 0,27% beträgt. Dieses Wahrscheinlichkeitsniveau wird in vielen Anwendungen verwendet. Selbstverständlich kann der Gate-Wert anders gewählt werden, z.B. „12". Allerdings werden bei einem solchen, höheren Gate-Wert auch mehr unwahrscheinlichere Zuordnungen zugelassen.

In der Kostenmatrix werden nun diejenigen Matrixelemente bestimmt, deren aufaddierte Kosten ein Minimum ergeben, oder anders gesagt, diejenigen Sensortrack-Paarungen bestimmt, deren Kosten Ci, j die kleinste Summe ergeben. Die Berechnung des Kostenminimums wird mittels des Auction-Algorithmus oder nach der JVC-(Jonker Vogelnant Castaflon) Methode durchgeführt. Diese Methoden sind in „D.P. Bertsekas, Auction algorithms for network flow problems : A tutorial introduction, Computational Optimization and Application,

Vol. 1, pp.7-66, 1992" und „O.E. Drummond, D. A. Castanon, M. S. Bellovin, Comparison of 2-D assignment algorithms for rectangular, floating point cost matrices, Proc . of SDI Panls on Tracking, No. 4, Dec. 1990 pp.81-97" beschrieben.

In dem Zahlenbeispiel der Fig. 3 ergibt sich das Kostenminimum für die in Fig. 3 eingekreisten Matrixelemente. Diese Matrixelemente mit den geringsten Kosten bestimmen die Zuordnungspaarungen der Sensortracks. Im Ausfuhrungsbeispiel der Fig. 3 werden somit die Sensortracks A2 und B2 mit den Kosten C ₂₂ =3 und die Sensortracks A3 und Bl mit den Kosten C ₃ i=8 einander paarweise zugeordnet.

Die Fusion der einander zugeordneten Sensortracks, im Beispiel die Sensortracks A2 und B2 einerseits und A3 und Bl anderseits, und damit die Generierung des fusionierten zentralen Tracks erfolgt durch konvexe Addition. Ist beispielsweise Xi der Zustandsvektor des einen Sensortracks und Pi die zugehörige Kovarianzmatrix und X ₂ der Zustandsvektor des anderen Sensortracks und P ₂ die zugehörige Kovarianzmatrix, so ergibt sich der Zustandsvektor X des fusionierten Tracks zu:

X = P ₂ (P ₁ + Pj- ¹ X ₁ + P ₁ (P ₁ + Pj- ¹ X ₂ (2) oder

X = P (P ₁ - ¹ X ₁ + P ₂ - ¹ X ₂ ) (3) und die Kovarianzmatrix P des fusionierten Tracks zu:

P = P ₁ (P ₁ + Pj- ¹ P ₂ (4), wobei für den Zustandsvektor gilt:

X = (5;

θ und für die Konvarianzmatrix gilt:

In Fig. 4 sind beispielhaft die in der Kostenmatrix der Fig. 3 enthaltenen Sensortracks Al, A2, A3 und Bl, B2 und B3 als Funktion des Peilwinkels B über der Zeit t schematisiert skizziert. Die durch Fusion der Sensoren A ₃ und B ₁ sowie A ₂ und B ₂ sich ergebenden Kontakttracks sind in Fig. 4 strichliniert angedeutet. In der Trackdatei 11 werden somit die als zentrale Tracks enthaltenen Sensortracks A2 und A3 mit den fusionierten Tracks überschrieben, der Sensortrack Al bleibt erhalten und der nicht zugeordnete Sensortrack B3 wird ebenfalls in die zentrale Trackdatei 11 als weiterer zentraler Track eingeschrieben. Wäre kein weiterer Sensor vorhanden, so würden nunmehr die in der zentralen Trackdatei 11 enthaltenen zentralen Tracks als Kontakttracks ausgegeben und im Display des Sonarbedieners in der in Fig. 4 dargestellten Weise angezeigt, so dass die Kontakttracks Al, B3 sowie die aus der Fusion hervorgegangenen Kontakttracks, die strichliniert dargestellt sind, sichtbar werden.

Bei Zielverlust, d.h. wenn einer oder mehrere Sensoren 10 das Ziel nicht mehr erfasst bzw. erfassen, enden auch die Kontakttracks. Bei erneuter Zielauffassung durch die Sensoren 10, sei es bei Wiederauffassung des alten Ziels oder Auffassung eines neuen Ziels, entstehen neue Kontakttracks, wie dies im Display 13 der Fig. 1 illustriert ist. Diese zeitlich versetzten Kontakttracks werden nunmehr auf die Wahrscheinlichkeit ihrer Zusammengehörigkeit geprüft, d.h. geprüft, ob ein neu beginnender Kontakttrack die Fortsetzung eines zuvor geendeten Kontakttracks ist oder der neu

beginnende Kontakttrack ein neu aufgefasstes Ziel markiert. Hierzu werden in vorgegebenen Zeitintervallen, z.B. alle 2 min, die zu einem aktuellen Zeitpunkt t _a (Fig. 1) vorhandenen Kontakttracks als neue Kontakttracks NT und alle Kontakttracks, deren Trackende vor dem aktuellen Zeitpunkt t _a liegt, als alte Kontakttracks OT aufgesucht und alle neuen und alten Kontakttracks NT, OT mittels Kostenrechnung auf paarweise Zugehörigkeit zueinander geprüft. Im Beispiel der Fig. 1 sind zum aktuellen Zeitpunkt t _a die neuen Tracks NTi, NT ₂ , NT ₃ und NT ₄ vorhanden und besitzen die alten Kontakttracks OTi, OT ₂ und OT ₃ ein Ende. Mit allen zum aktuellen Zeitpunkt t _a gefundenen j neuen Kontakttracks NT-, und i alten Kontakttracks OT ₁ wird eine zweiteilige Kostenmatrix mit einem Zuordnungsteil und einem Nichtzuordnungsteil aufgestellt, wie sie in Fig. 5 illustriert ist. In den Zuordnungsteil werden die Kosten C ₁ ,-, der wechselseitigen Zuordnung der j neuen und i alten Kontakttracks NT bzw. OT und in den Nichtzuordnungsteil die Kosten C _n der NichtZuordnung der i alten Kontakttrack OTi eingetragen. Zur Kostenberechnung werden für alle gefundenen neuen Kontakttracks NT ₃ für die Zeitlücken zwischen dem Zeitpunkt ihres Beginns und dem Zeitpunkt des Endes der alten

Kontakttracks OT ₁ Peilungen Bj. _etro retrodiziert und für alle gefundenen alten Kontakttracks OT ₁ für die Zeitlücken zwischen dem Zeitpunkt ihres Endes und dem Zeitpunkt des Beginns der neuen Kontakttracks NT ₃ Peilungen B _pra prädiziert.

Daraus werden die Zuordnungskosten C^ ₃ für einen Zeitpunkt tk gemäß

berechnet. σ _R und σ _R sind die Standardabweichungen der pradizierten bzw. retrodizierten Peilungen der Kontakttracks einer Kontakttrack-Paarung OT ₁ und NT-,. Die in die Kostenmatrix gemäß Fig. 2 einzutragenden Kosten C ₁ ,-, sind dann die gemittelten Kosten zu dem Zeitpunkten ti, der durch das zeitliche Ende des alten Kontakttrack OT ₁ bestimmt ist, für einen Zeitpunkt t _N , der durch den Beginn des neuen

Kontakttrack NT-, bestimmt ist und für einen Zeitpunkt t _{N /2} , der durch die halbe Zeitdifferenz zwischen dem Ende ti des alten und dem Beginn t _M ^w prä des neuen Kontakttracks bestimmt i s t , gemäß

/ 2

C W,3 = - — ₃ ( Vc- ¹ io + ⁺ C ^{U N} 1O ^pra c ^Np _" ( 8 )

Vor der Retrodiktion der Peilungen in die Zeitlucken werden die neuen Kontakttracks NT ₃ durch eine rückwärts gerichtete Glattung, vorzugsweise nach der Rauch-Tung-Striebel-Methode, geglättet. Der Kosteneintrag C _n für die NichtZuordnung der alten Kontakttracks OT ₁ zu den neuen Kontakttracks NT _D wird jeweils wieder als fester Wert vorgegeben, wobei beispielhaft C _n =9 gewählt wird.

Anstelle der Berechnung zu den drei angegebenen Zeitpunkten ti, t _N und t _{N /2} können auch die Zuordnungskosten für eine

Vielzahl von vorzugsweise aquidistanten Zeitpunkten t* zwischen dem Ende ti des alten Kontakttrack und dem Beginn t _H des neuen Kontakttrack berechnet werden und zum

Kosteneintrag C ₁ ,-, in den Zuordnungsteil der Kostenmatrix die gemittelte Summe der Zuordnungskosten C^. verwendet werden, gemäß

Sind die Kosten C ₁ ,-, für die Zuordnung ebenso wie die Kosten C _n für die NichtZuordnung in die Kostenmatrix gemäß Fig. 5 eingetragen, so werden wiederum diejenigen Matrixelemente bzw. Kontakttrack-Paarungen bestimmt, deren aufaddierte Kosten ein Minimum ergeben. Hierzu wird wiederum der Auction- Algorithmus oder der JVC-Algorithmus eingesetzt. Diejenigen Matrixelemente, deren aufaddierte Kosten ein Minimum ergeben, definieren die Kontakttrack-Paarungen OT ₁ , NT-,, die einander zugeordnet und als Zieltrack ausgewiesen werden. Wie bereits erwähnt, werden in der bildlichen Darstellung der Kontakttracks die Zeitlucken zwischen den einander zugeordneten alten und neuen Kontakttracks OT ₁ und NT ₃ durch Linien ausgefüllt, die das Ende des zugeordneten alten Kontakttrack und den Anfang des zugeordneten neuen Kontakttracks miteinander verbinden. Diese Linien sind vorzugsweise Geraden. Die miteinander verbundenen, zugeordneten Kontakttracks und die ununterbrochenen Kontakttracks sind Zieltracks, die jeweils eine einem Ziel zugehörige Spur von Peilungen ausweisen.

Um unwahrscheinliche Zuordnungen von neuen Kontakttracks NT und alten Kontakttracks OT auszuschließen und die Berechnungszeit für das Kostenminimum aus der Kostenmatrix zu reduzieren, werden die in die Kostenmatrix einzutragenden Kosten C ₁ , _j für die Zuordnung für eine Kontakttrack-Paarung eines alten und eines neuen Kontakttrack mit einem willkürlich vorgegebenen Maximalwert angesetzt, wenn die Zeitdifferenz zwischen dem Beginn des neuen Kontakttracks NT-,

und dem Ende des alten Kontakttrack OT ₁ großer ist als ein Vorgabewert, z.B. 10 min.

Zur Kostenberechnung können auch anstelle der pradizierten und retrodizierten Peilungen Peilraten herangezogen werden. In diesem Fall werden für alle gefundenen neuen Kontakttracks NT-, für die Zeitlucken zwischen dem Zeitpunkt ihres Beginns und dem Zeitpunkt des Endes der alten Kontakttracks Peilraten

BRj. _etro retrodiziert und für alle gefundenen alten

Kontakttracks OT ₁ für die Zeitlucken zwischen dem Zeitpunkt ihres Endes und dem Zeitpunkt des Beginns der neuen

Kontakttracks NT ₃ Peilraten BR _pra pradiziert und die Zuordnungskosten für einen Zeitpunkt t _k gemäß

berechnet. σ _BR und σ _BR sind die Standardabweichungen der pradizierten und retrodizxerten Peilraten. Zur Eintragung in die Kostenmatrix gelangen dann die Kosten C ₁ ,-, als Mittel aus den Zuordnungskosten C^ ₃ , die nach Gl. (10) für den durch das zeitliche Ende des alten Kontakttracks OT ₁ bestimmten Zeitpunkt k=ti und für den durch den Beginn des neuen Kontakttracks NT-, bestimmten Zeitpunkt k=t _N , berechnet werden, gemäß

Wie bei Verwendung der Peilungen wird auch hier vor der Retrodiktion der Peilraten in die Zeitlucken der jeweils

gefundene neue Kontakttrack; durch eine rückwärts gerichtete Glattung, vorzugsweise nach der Rauch-Tung-Striebel-Methode, geglättet .

Alternativ kann zur Kostenberechnung auch das Nutz-Stor- Verhaltms, im folgenden SNR (Signal-Noise-Ration) genannt, der z.B. von einem CAS BDT-Sensor abgenommenen Peilungen verwendet werden, dass als Attribut den Peilungen zugeordnet ist. In diesem Fall werden für alle gefundene neuen Kontakttracks NT-, für die Zeitlucken zwischen dem Zeitpunkt ihres Beginns und dem Zeitpunkt des Endes der alten Kontakttracks OT ₁ SNR-Werte SNR _retro retrodiziert und für alle gefundenen alten Kontakttracks OT ₁ für die Zeitlucken zwischen dem Zeitpunkt ihres Endes und dem Zeitpunkt des Beginns der neuen Kontakttracks NT-, SNR-Werte SNR _prä pradiziert. Die Zuordnungskosten C^. für einen Zeitpunkt t* werden aus der Differenz zwischen den für diesen Zeitpunkt t _k pradizierten und retrodizierten SNR-Werte SNR _pra und SNR _re tro und deren Standardabweichungen σ _SNR und σ _SNR der

Kontakttracks einer Kontakttrack-Paarung gemäß

C* = ^SNR p _ra (t _k ) ^~ SNR _retro (t _k ) σL _R _ . (t _k ) + σL _Mtrβ (t _k )

berechnet. Die in die Kostenmatrix gemäß Fig. 2 einzutragenden Kosten C ₁ , ₃ sind dann die gemittelte Summe gemäß Gl. (9) der Zuordnungskosten C^. , die für eine Vielzahl von vorzugsweise aquidistanten Zeitpunkten t _k zwischen dem Ende ti des alten Kontakttracks OT ₁ und dem Beginn t _{N a} des neuen Kontakttrack für eine Kontakttrack-Paarung gemäß Gl. (12) berechnet worden sind. Alternativ können zur

Eintragung in die Kostenmatrix auch die Kosten C ₁ ,-, herangezogen werden, die gemäß Gl. (11) das Mittel aus den

Zuordnungskosten C^ ₃ , bilden die wiederum aus Gl. (12) für den durch das zeitliche Ende des alten Kontakttracks OTi bestimmten Zeitpunkt k=ti und für den durch den Beginn des neuen Kontakttracks NT ₃ bestimmten Zeitpunkt k = t _N berechnet werden.

In einer alternativen Verfahrensvariante werden zur Kostenberechnung der Kosten C ₁ ,-, für die Zuordnung für alle gefundenen neuen und alten Kontakttracks NT _j bzw. OT ₁ deren Frequenzlinien FL _NT bzw. FL _0T aufgesucht, die üblicherweise den Peilungen, also den Kontakttracks, als Attribute zugeordnet sind. Die Zuordnungskosten werden nunmehr aus der Frequenzdifferenz der Frequenzlinien FL _0T und FL _NT und deren Standardabweichung σ _FLoτ und σ _FLm der Kontakttracks einer

Kontakttrack-Paarung berechnet, also gemäß

Bei einer gleichen Anzahl k von Frequenzlinien in jedem Kontakttrack eines Kontakttracks der Kontakttrack-Paarungen mit k=l,2 ... N _FL wird zum Kosteneintrag in den Zuordnungsteil der Kostenmatrix die Summe der für alle k Frequenzlinien berechneten Zuordnungskosten gemittelt, gemäß

<"

_

Besitzen dagegen die alten und neuen Kontakttracks eine unterschiedliche Anzahl von Frequenzlinien, so werden die

Zuordnungskosten C^ _: nur für diejenigen Kontakttrack- Paarungen von alten und neuen Kontakttracks OT ₁ , NT-, berechnet, bei denen Frequenzlinienpaare mit einer gegenüber einem Vorgabewert kleineren Frequenzdifferenz auftreten, k = 1, 2 ... N _Paare ist dabei die Numerierung der gefundenen Frequenzlinienpaare. Beispielhaft wird als Vorgabewert eine Frequenzdifferenz von 2Hz angenommen. Die in den Zuordnungsteil der Kostenmatrix einzutragenden Kosten C ₁ ,-, für die Zuordnung berechnen sich dann aus der gemittelten Summe der für alle N _Paare Frequenzlinienpaare berechneten Zuordnungskosten dividiert durch die quadrierte Anzahl von gefundenen Frequenzlinienpaaren gemäß

Außerdem werden in den Zuordnungsteil der Kostenmatrix die Kosten für alle Kontakttrack-Paarungen, die keine solche Frequenzlinienpaare aufweisen, auf ein willkürliches Maximum gesetzt .

Nach Aufstellen dieser Kostenmatrix wird in gleicher Weise das Minimum der aufaddierten Kosten der Matrixelemente nach dem Auction-Algorithmus oder dem JVC-Algorithmus bestimmt und wie beschrieben die Zuordnung der Kontakttracks vorgenommen.