EP1056063A1 - Method and device for the determination of the traffic conditions of a roadsegment - Google Patents

Method and device for the determination of the traffic conditions of a roadsegment Download PDF

Info

Publication number
EP1056063A1
EP1056063A1 EP00108359A EP00108359A EP1056063A1 EP 1056063 A1 EP1056063 A1 EP 1056063A1 EP 00108359 A EP00108359 A EP 00108359A EP 00108359 A EP00108359 A EP 00108359A EP 1056063 A1 EP1056063 A1 EP 1056063A1
Authority
EP
European Patent Office
Prior art keywords
traffic
model
section
segment
vector
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.)
Granted
Application number
EP00108359A
Other languages
German (de)
French (fr)
Other versions
EP1056063B1 (en
Inventor
Andreas Lagger
Kai Müller
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Siemens Schweiz AG
Original Assignee
Siemens Schweiz 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 Siemens Schweiz AG filed Critical Siemens Schweiz AG
Publication of EP1056063A1 publication Critical patent/EP1056063A1/en
Application granted granted Critical
Publication of EP1056063B1 publication Critical patent/EP1056063B1/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/01Detecting movement of traffic to be counted or controlled
    • G08G1/0104Measuring and analyzing of parameters relative to traffic conditions

Definitions

  • the present invention relates to a method and a device according to the preamble of the claim 1 or 6.
  • Measuring devices are often used for automatic monitoring and control of road traffic. which measure the traffic volume or the number of vehicles that are within a certain Time pass a measuring point. Based on the measurement results, however, only a useful statement about the local traffic volume can be made. About the condition It is not possible to specify the traffic in the traffic section before the measurement point mentioned.
  • traffic measurements are therefore made at its entrance and exit serving sensors are provided. By comparing the measurement results of the two sensors can make a statement about the traffic condition within the monitored traffic segment be made. Because the measured traffic volume shows strong fluctuations, the stochastic Faults around a "true" value occur, the measured values are filtered.
  • Kalman filters or so-called observers use a model as well as the input and output information of the object whose condition is to be determined or tracked. Basically, the The more information, the faster and more accurately the actual object states can be determined of the object can be used for the model, which shows the relationships between the inputs and Output variables should reflect correctly. From [1] the use of a model is known that is based on the use of eight parameters, some of which are difficult to determine. A part the parameter also depends on the weather. The results achieved with this system, Despite the increased effort, therefore, deviate significantly from the results, which ideally let achieve.
  • the present invention is therefore based on the object of a method and an apparatus for to quickly determine the state of a traffic segment.
  • a Kalman filter with a model of the monitored traffic section is used.
  • the traffic condition can be determined quickly and precisely using a Kalman KF filter.
  • Faults are recognized almost without delay, which means that the necessary measures are taken (Traffic jam warning, etc.) can be issued quickly. Further advantages of the invention or of the selected model are explained below.
  • Models differ in their input and output sizes and their interdependency as well as the choice of state variables that reflect the "inner state” of the model. Out For reasons of measurement technology there is a location and time discretization of the continuous traffic instead of. Measured values are only available at certain locations and at certain times (sampling times) Available.
  • a transport-oriented model is used for the state assessment.
  • the monitored traffic segment is accordingly regarded as a transport medium, which, analogous to a conveyor belt, elements or vehicles are transported at a certain speed or let it happen. To simplify matters, it is provided that all vehicles are in the traffic section moving at a constant speed.
  • a Kalman filter a quick and reliable determination of traffic disruptions and other parameters from strongly fluctuating measurements of traffic volume q and speed v.
  • the model used describes the conditions of an "ideal" street, the one has infinite capacity. This property in particular enables faults to be very quickly and safely recognizable, since measured values in the event of a fault differ from the "predicted" values of the model Differentiate values strongly.
  • the model is not suitable for simulating road traffic. The use of this model is therefore only useful in connection with a Kalman filter. Because of however, the linearity of the transport model is optimally suited for this purpose.
  • a major advantage of the method according to the invention is that all parameters of the Model results in a simple manner from the geometric arrangement of the measuring points. An expensive one Identification of model parameters is therefore no longer necessary.
  • the parameters depend on the used Model does not depend on weather influences such as wetness or smoothness, which means that it can be adapted to different weather conditions Ambient conditions are superfluous.
  • the model shows a traffic model corresponding to a traffic segment, to which the traffic volume q1 measured at the entrance to the traffic segment and the speed v1 measured at the entrance to the traffic segment as well as the naturally constant length I s of the traffic segment are supplied as an input variable.
  • the model provides an output quantity q2 m , which, taking into account the calculated running time T L of the vehicles Fz, indicates the traffic volume at the exit of the traffic section. For the sake of simplicity, it is assumed that all vehicles in the section travel at the same constant speed v1.
  • the change in traffic volume q is thus analogous to the transport of material (e.g. sand) an assembly line moving at a constant speed v.
  • material e.g. sand
  • the material appears after a time delay that depends on the transport speed and the length of the conveyor belt dependent, unchanged at the point of delivery.
  • the vehicles F1, F2 and F3 appear according to the idealized one shown in FIG. 1 Model without differences in transit time, at the same mutual distance at the exit of the traffic section.
  • a transport process of this kind corresponds to the partial differential equation (1):
  • Kalman filter also requires a description of the process in Form of an ordinary and not a partial differential equation, which partial derivatives an unknown function with multiple variables (in relation to the difference from ordinary and partial differential equations see [3], page 435).
  • an ordinary differential equation is obtained by only considering discrete locations in the traffic section.
  • the traffic section is divided into n segments s1, s2, ..., s n , the inner limits of which do not have sensors or sensors.
  • the speed is assumed to be constant in each of these segments.
  • the first value of the trowel q 1, 0 is the value q1 at the input q1; the last value of the q 1, n is identical to the q2 value at the exit of the model.
  • the individual matrices are as follows:
  • the selected model is now used to describe the transport process and thus the undisturbed traffic suitable in a traffic segment.
  • This shape is also very suitable for a Kalman filter because it is linear in the state variables.
  • the model serves as a building block for a traffic section with two Measuring points. With the help of this module, models for arbitrarily complex topologies can be easily set up, which have entrances, exits, branches and junctions.
  • the measured values are available as discrete-time values for a specific measuring interval (sampling time Ts). Therefore, the discrete version of the model is chosen, which results from the continuous model as follows (the matrix C does not change):
  • the scanning step k was written as an index.
  • the index k denotes the interval for the time t kT S ⁇ t ⁇ ( k + 1) T S
  • the deviations between the predicted values of the model q2 m and the actually measured values q2 are corrected by the Kalman filter, which calculates the estimated values q2 e , in the undisturbed case. However, if malfunctions occur that were not taken into account in the modeling, the Kalman filter also supplies incorrect values. The accident is now taken into account by extending the model, as shown in FIG. 4, by virtual entrances and exits to a “parking space with infinite capacity”.
  • the traffic flows on the virtual exits which are designated q2 d , are a measure of the effects of a disruption that has occurred.
  • the Kalman filter estimates the values q2 d (value q2 de ).
  • the state vector x provided for four segments s1, s2, s3, s4 is expanded as follows:
  • the Kalman filter can use the virtual traffic flows to compare measured values and estimated values of the model.
  • the vector x (k) is corrected by the Kalman filter KF (see FIG. 5), which results in the vector x e (k) with corrected values x 1e , ..., x 4e and x de .
  • the corrected (internal) state values x 4e plus x de or x de thus correspond to the estimated values q2e and q2 de of the Kalman filter KF.
  • the values q2 de are an immediate measure of the degree of disturbance that has occurred in the monitored section (see FIG. 7).
  • the corrected model value q2 * m (k) thus corresponds to the value q2 (k) actually measured in the interval k.
  • the assumed disturbances result from the traffic process and correspond to the form of a frequency-independent noise with a Gaussian-distributed amplitude spectrum. If one looks at measured courses of traffic intensity q, one can see that the strong fluctuations can be represented quite well as a noise process around an average.
  • the Kalman KF filter is an optimal filter that measures the variance of the estimation error x err ( k ) ⁇ x ( k ) - x e ( k ) minimized.
  • the "true" state vector is not known. Therefore x e (k) is used as an estimate for the unknown state vector x (k).
  • the matrix G (k) which is exactly like the matrices A (k), B (k) by discretization arises, indicates how the (unknown) disturbance variables are distributed among the internal state variables of the process. In the simplest case, it is assumed that these disturbance variables have an even effect on all internal state variables.
  • the matrix G (k) therefore indicates where, but not with what intensity, the disturbance variables affect the model.
  • the intensity of the disturbances is described by the covariance matrices Q and R described below. If, for example, there are three state variables and two disturbance variables and the first disturbance variable only affects the first state variable and the second disturbance variable only affects the following two state variables, this would be taken into account with a matrix G (k) of the form given below:
  • the matrix G (k) is selected in such a way that process disturbances p all sizes of the state vector influence x (k) in the same way.
  • the control parameters of the Kalman filter KF are the covariance matrices Q and R of the (assumed) Noise processes that relate to the process itself or to the measured values output by the sensors act.
  • the matrix Q describes the intensity of the process disturbances p.
  • the matrix R describes the Intensity of sensor disturbances see Although the matrices Q and R may theoretically be time-variant, assumed that the fluctuations in the measured values regardless of the traffic situation and thus are independent of time.
  • the matrices Q and R thus take on constant values.
  • the elements of the matrices Q and R are the design parameters of the Kalman KF filter, with which the settling time and sensitivity to interference of the Kalman filter KF can be set.
  • the signal flow in the Kalman filter KF according to the invention is shown in FIG. 5.
  • the products of the input value q1 (k) times matrix B (k) and the estimated state vector x e (k) times matrix A (k) are added in the addition stage ADD1.
  • the difference level DIFF is then used to form the difference between the traffic intensity q2 m (k) newly determined by the model and the actually measured traffic intensity q2 (k) (at the exit of the traffic segment).
  • the resulting difference (q2 m (k) - q2 (k)) is multiplied by the Kalman matrix L (k), whereby values are formed with which the state vector x (k) is corrected using the addition stage ADD3 or in the estimated values State vector x e (k) is converted.
  • C 1 the sum of the estimated values x4 e and x de is formed, which gives the estimated value q2 e (k).
  • the matrices C 1 and C 2 are constant.
  • the last element of the matrix C 2 is negative, so that a positive value for q2 de (k) is obtained in the event of a jam.
  • a computer system known to the person skilled in the art (for example: processor, signal processor).
  • FIG. 7 shows the traffic intensity q2 e estimated by the Kalman filter KF and the estimated traffic intensity q2 de on the virtual entrance and exit, influenced by a disturbance that occurs at a time 11 and during the time T dist up to one Time t2 continues.
  • the vehicles can pass through the traffic section unhindered, so that the traffic volumes determined at the entrance and exit correspond to the value q1, neglecting the travel time differences.
  • a disturbance occurs, which creates a bottleneck within the traffic section, which only allows traffic to pass with the reduced traffic volume q2 dist .
  • the fault in the traffic control center must be recognized immediately.
  • the device according to the invention and the method by means of which input and output variables q1, q2 are processed using the Kalman filter KF now allow the occurrence of the fault to be detected with practically no delay.
  • the value q2 de starts up immediately after the fault occurs. This corresponds to an increase in the estimated traffic volume q2 de on the entrance to the virtual parking lot.
  • the traffic volume q2 de increases again from the virtual parking lot back into the traffic section and then falls back to zero when the traffic volumes q1 and q2 are adjusted at the entrance and exit of the traffic section.
  • the threshold values are chosen according to the size of the fault to be detected. Preferably several threshold values are provided which correspond to the states “light traffic obstruction” ,, “slow traffic” or “traffic jam”. When the corresponding conditions occur therefore the appropriate measures (e.g. traffic jam warning) are initialized.
  • the expected ones can also be expected Calculate travel times quickly and precisely. This can also be optimal for road users Itineraries to be set.

Abstract

The method involves using a Kalman filter based on a model that can take into account the length of the traffic section and measured levels of traffic at the entry to and exit from the section. A transport-oriented model is selected in which a vector i formed contg. a correction value and a traffic level value for each of at least two segments using the length of the traffic sec and the sampled levels of traffic at the entry to and exit from the section. The vector is corrected using the Kalman filter to provide a precise estimate of the exit traffic level and of the traffic state. An Independent claim is also included for an arrangement for determining the traffic state of a section of traffic.

Description

Die vorliegende Erfindung betrifft ein Verfahren und eine Vorrichtung nach dem Oberbegriff des Patentanspruchs 1 bzw. 6.The present invention relates to a method and a device according to the preamble of the claim 1 or 6.

Zur automatischen Überwachung und Steuerung des Strassenverkehrs werden oft Messgeräte eingesetzt, welche die Verkehrsstärke bzw. die Zahl der Fahrzeuge messen, die innerhalb einer bestimmten Zeitdauer einen Messpunkt passieren. Anhand der dabei anfallenden Messresultate kann jedoch nur eine brauchbare Aussage über die lokal auftretende Verkehrsstärke gemacht werden. Über den Zustand des Verkehrs in dem vor dem genannten Messpunkt liegenden Verkehrsabschnitt ist keine Angabe möglich.Measuring devices are often used for automatic monitoring and control of road traffic. which measure the traffic volume or the number of vehicles that are within a certain Time pass a measuring point. Based on the measurement results, however, only a useful statement about the local traffic volume can be made. About the condition It is not possible to specify the traffic in the traffic section before the measurement point mentioned.

Auch anhand der Messung der Fahrzeuggeschwindigkeit bzw. Geschwindigkeitstendenz ist eine Bestimmung des Verkehrszustands möglich. Störungen im Verkehrsablauf sind aufgrund der Geschwindigkeitsmesswerte jedoch erst so spät feststellbar, dass entsprechende Massnahmen zur Verkehrsregulierung der Störung nicht mehr entgegenwirken können.A determination is also based on the measurement of the vehicle speed or speed tendency of the traffic condition possible. Traffic disturbances are due to the speed measurements However, it can only be determined so late that appropriate measures to regulate traffic can no longer counteract the disturbance.

Zur Überwachung eines Verkehrsabschnittes werden daher an dessen Ein- und Ausgang der Verkehrsmessung dienende Messfühler vorgesehen. Durch einen Vergleich der Messresultate der beiden Messfühler kann eine Aussage über den Verkehrszustand innerhalb dem überwachten Verkehrsabschnitt gemacht werden. Da die gemessene Verkehrsstärke starke Schwankungen aufweist, die als stochastische Störungen um einen "wahren" Wert auftreten, werden die gemessenen Werte gefiltert.To monitor a traffic segment, traffic measurements are therefore made at its entrance and exit serving sensors are provided. By comparing the measurement results of the two sensors can make a statement about the traffic condition within the monitored traffic segment be made. Because the measured traffic volume shows strong fluctuations, the stochastic Faults around a "true" value occur, the measured values are filtered.

Bei der Filterung anhand von klassischen Tiefpass-, Bandpass- oder Hochpassfiltern entstehen Verzögerungen, die verhindem, dass eine rasche Aussage über den Verkehrszustand gemacht werden kann. Bekannt ist, dass Messwerte anhand eines Kalman-Filters praktisch verzögerungsfrei gefiltert werden können.Delays occur when filtering using classic low-pass, band-pass or high-pass filters, which prevent a quick statement about the traffic condition from being made. It is known that measured values are filtered practically without delay using a Kalman filter can.

Kalman-Filter oder sogenannte Beobachter verwenden ein Modell sowie die Ein- und Ausgangsinformationen des Objektes, dessen Zustand festgestellt bzw. verfolgt werden soll. Grundsätzlich gilt, dass die tatsächlichen Objektzustände umso schneller und genauer bestimmt werden können, je mehr Informationen des Objektes für das Modell verwendet werden, welches die Zusammenhänge zwischen den Ein- und Ausgangsgrössen korrekt widergeben soll. Aus [1] ist die Verwendung eines Modells bekannt, das auf der Verwendung von acht Parametern basiert, die sich zum Teil nur schwer ermitteln lassen. Ein Teil der Parameter hängt zudem von der Witterung ab. Die erreichten Resultate können bei diesem System, trotz des erhöhten Aufwandes, daher erheblich von den Resultaten abweichen, die sich idealerweise erzielen lassen.Kalman filters or so-called observers use a model as well as the input and output information of the object whose condition is to be determined or tracked. Basically, the The more information, the faster and more accurately the actual object states can be determined of the object can be used for the model, which shows the relationships between the inputs and Output variables should reflect correctly. From [1] the use of a model is known that is based on the use of eight parameters, some of which are difficult to determine. A part the parameter also depends on the weather. The results achieved with this system, Despite the increased effort, therefore, deviate significantly from the results, which ideally let achieve.

Der vorliegenden Erfindung liegt daher die Aufgabe zugrunde, ein Verfahren und eine Vorrichtung zur schnellen Ermittlung des Zustandes eines Verkehrsabschnittes anzugeben. The present invention is therefore based on the object of a method and an apparatus for to quickly determine the state of a traffic segment.

Diese Aufgabe wird durch die im kennzeichnenden Teil des Patentanspruchs 1 bzw. 6 angegebenen Massnahmen gelöst. Vorteilhafte Ausgestaltungen der Erfindung sind in weiteren Ansprüchen angegeben.This object is achieved by those specified in the characterizing part of patent claim 1 and 6, respectively Measures solved. Advantageous embodiments of the invention are specified in further claims.

Erfindungsgemäss wird ein Kalman-Filter mit einem Modell des überwachten Verkehrsabschnittes verwendet. Durch die vorteilhafte Wahl eines einfachen Modells für den störungsfreien Fall in Verbindung mit einem Kalman-Filter KF kann der Verkehrszustand jeweils schnell und präzise ermittelt werden. Insbesondere werden Störungen nahezu verzögerungsfrei erkannt, wodurch jeweils die notwendigen Massnahmen (Stauwarnung, etc.) schnell abgegeben werden können. Weitere Vorteile der Erfindung bzw. des gewählten Modells werden nachfolgend erläutert.According to the invention, a Kalman filter with a model of the monitored traffic section is used. By the advantageous choice of a simple model for the trouble-free case in connection The traffic condition can be determined quickly and precisely using a Kalman KF filter. In particular Faults are recognized almost without delay, which means that the necessary measures are taken (Traffic jam warning, etc.) can be issued quickly. Further advantages of the invention or of the selected model are explained below.

Modelle unterscheiden sich durch ihre Ein- und Ausgangsgrössen und deren gegenseitige Abhängigkeit sowie durch die Wahl der Zustandsgrössen, die den "inneren Zustand" des Modells widerspiegeln. Aus messtechnischen Gründen findet eine Orts- und eine Zeitdiskretisierung des kontinuierlichen Verkehrs statt. Messwerte stehen nur an bestimmten Orten sowie zu bestimmten Zeiten (Abtastzeitpunkten) zur Verfügung. Erfindungsgemäss wird für die Zustandsabschätzung ein transportorientiertes Modell verwendet. Der überwachte Verkehrsabschnitt wird demgemäss als Transportmedium betrachtet, welches, analog zu einem Förderband, Elemente bzw. Fahrzeuge mit einer bestimmten Geschwindigkeit transportiert bzw. passieren lässt. Vereinfachend wird vorgesehen, dass sich alle Fahrzeuge in dem Verkehrsabschnitt mit einer konstanten Geschwindigkeit fortbewegen. Trotz des einfachen Aufbaus des verwendeten Transportmodells ermöglicht das erfindungsgemässe Verfahren mit Hilfe eines Kalman-Filters eine schnelle und zuverlässige Bestimmung von Verkehrsstörungen und anderer Kenngrössen aus stark schwankenden Messwerten der Verkehrsstärke q und der Geschwindigkeit v.Models differ in their input and output sizes and their interdependency as well as the choice of state variables that reflect the "inner state" of the model. Out For reasons of measurement technology there is a location and time discretization of the continuous traffic instead of. Measured values are only available at certain locations and at certain times (sampling times) Available. According to the invention, a transport-oriented model is used for the state assessment. The monitored traffic segment is accordingly regarded as a transport medium, which, analogous to a conveyor belt, elements or vehicles are transported at a certain speed or let it happen. To simplify matters, it is provided that all vehicles are in the traffic section moving at a constant speed. Despite the simple structure of the used transport model enables the inventive method with the help of a Kalman filter a quick and reliable determination of traffic disruptions and other parameters from strongly fluctuating measurements of traffic volume q and speed v.

Wesentlich ist, dass das verwendete Modell die Verhältnisse einer "idealen" Strasse beschreibt, die eine unendliche Kapazität aufweist. Gerade diese Eigenschaft ermöglicht es, Störungen sehr schnell und sicher zu erkennen, da sich gemessene Werte im Störungsfall von den durch das Modell "vorhergesagten" Werten stark unterscheiden. Das Modell ist hingegen für die Simulation des Strassenverkehrs nicht geeignet. Der Einsatz dieses Modells ist daher nur in Verbindung mit einem Kalman-Filter sinnvoll. Aufgrund dessen Linearität ist das Transportmodell zu diesem Zweck hingegen optimal geeignet.It is essential that the model used describes the conditions of an "ideal" street, the one has infinite capacity. This property in particular enables faults to be very quickly and safely recognizable, since measured values in the event of a fault differ from the "predicted" values of the model Differentiate values strongly. However, the model is not suitable for simulating road traffic. The use of this model is therefore only useful in connection with a Kalman filter. Because of however, the linearity of the transport model is optimally suited for this purpose.

Ein wesentlicher Vorteil des erfindungsgemässen Verfahrens besteht darin, dass sich alle Parameter des Modells in einfacher Weise aus der geometrischen Anordnung der Messstellen ergibt. Eine aufwendige Identifikation von Modellparametem entfällt damit. Insbesondere hängen die Parameter des verwendeten Modells nicht von Witterungseinflüssen wie Nässe oder Glätte ab, wodurch sich eine Anpassung an unterschiedliche Umgebungsbedingungen erübrigt.A major advantage of the method according to the invention is that all parameters of the Model results in a simple manner from the geometric arrangement of the measuring points. An expensive one Identification of model parameters is therefore no longer necessary. In particular, the parameters depend on the used Model does not depend on weather influences such as wetness or smoothness, which means that it can be adapted to different weather conditions Ambient conditions are superfluous.

Die Erfindung wird nachfolgend anhand einer Zeichnung beispielsweise näher erläutert. Dabei zeigt:

Fig. 1
ein Verkehrsmodell mit Ein- und Ausgangsgrössen q1, q2,
Fig. 2
das Verkehrsmodell gemäss Fig. 1 für einen Verkehrsabschnitt, den Fahrzeuge mit einer bestimmten Gruppenlaufzeit durchfahren,
Fig. 3
den Verlauf der Verkehrsstärke q2 am Ausgang des Verkehrsabschnittes als Reaktion eines synthetischen Verlaufs der Verkehrsstärke q1 am Eingang des Verkehrsabschnittes für das kontinuierliche Modell,
Fig. 4
das Verkehrsmodell mit virtueller Ein- und Ausfahrt,
Fig. 5
das verwendete Kalman-Filter KF mit einem erweiterten diskreten Modell MOD,
Fig. 6
den Verlauf der Verkehrsstärke q2 am Ausgang des Verkehrsabschnittes als Reaktion eines synthetischen Verlaufs der Verkehrsstärke q1 am Eingang des Verkehrsabschnittes für das diskrete Modell und
Fig. 7
den vom Kalman-Filter KF geschätzten Verlauf der Verkehrsstärke q2e (verzögerungsfreie Glättung von q2) sowie die Verkehrsstärke q2d auf den virtuellen Ein- und Ausfahrten beim Auftreten einer Störung (Wert q2d durch das Kalman-Filter KF geschätzt, wird mit q2de bezeichnet).
The invention is explained in more detail below using a drawing, for example. It shows:
Fig. 1
a traffic model with input and output variables q1, q2,
Fig. 2
1 for a traffic segment that vehicles with a certain group transit time are traveling through,
Fig. 3
the course of the traffic volume q2 at the exit of the traffic section in response to a synthetic course of the traffic strength q1 at the entrance of the traffic section for the continuous model,
Fig. 4
the traffic model with virtual entry and exit,
Fig. 5
the Kalman filter KF used with an extended discrete model MOD,
Fig. 6
the course of the traffic volume q2 at the exit of the traffic section in response to a synthetic course of the traffic intensity q1 at the entrance of the traffic section for the discrete model and
Fig. 7
the course of the traffic intensity q2 e estimated by the Kalman filter KF (instantaneous smoothing of q2) and the traffic intensity q2 d on the virtual entrances and exits when a fault occurs (value q2 d estimated by the Kalman filter KF is estimated with q2 de designated).

Einleitend werden die verwendeten Bezeichnungen tabellarisch aufgelistet:

Figure 00030001
To begin with, the terms used are listed in a table:
Figure 00030001

Fig. 1 zeigt ein zu einem Verkehrsabschnitt korrespondierendes Verkehrsmodell, dem als Eingangsgrösse die am Eingang des Verkehrsabschnittes gemessene Verkehrsstärke q1 sowie als Parameter die am Eingang des Verkehrsabschnittes gemessene Geschwindigkeit v1 sowie die selbstverständlich konstante Länge Is des Verkehrsabschnittes zugeführt werden. Das Modell liefert eine Ausgangsgrösse q2m, welche, unter Berücksichtigung der errechneten Laufzeit TL der Fahrzeuge Fz, die Verkehrsstärke am Ausgang des Verkehrsabschnittes angibt. Für das Modell wird vereinfachend angenommen, dass sich alle Fahrzeuge in dem Abschnitt mit der gleichen konstanten Geschwindigkeit v1 fortbewegen.1 shows a traffic model corresponding to a traffic segment, to which the traffic volume q1 measured at the entrance to the traffic segment and the speed v1 measured at the entrance to the traffic segment as well as the naturally constant length I s of the traffic segment are supplied as an input variable. The model provides an output quantity q2 m , which, taking into account the calculated running time T L of the vehicles Fz, indicates the traffic volume at the exit of the traffic section. For the sake of simplicity, it is assumed that all vehicles in the section travel at the same constant speed v1.

Die Änderung der Verkehrsstärke q wird somit in Analogie zum Transport von Material (z.B. Sand) auf einem mit konstanter Geschwindigkeit v bewegten Fliessband beschrieben. Das Material erscheint nach einer zeitlichen Verzögerung, die von der Transportgeschwindigkeit und der Länge des Fliessbandes abhängig ist, unverändert am Abgabeort. Es bestehen keine Laufzeitunterschiede der transportierten Materialpartikel. Deshalb erscheinen die Fahrzeuge F1, F2 und F3 gemäss dem in Fig. 1 gezeigten idealisierten Modell ohne Laufzeitunterschiede, in gleichem gegenseitigen Abstand am Ausgang des Verkehrsabschnittes.The change in traffic volume q is thus analogous to the transport of material (e.g. sand) an assembly line moving at a constant speed v. The material appears after a time delay that depends on the transport speed and the length of the conveyor belt dependent, unchanged at the point of delivery. There are no differences in the transit times of the transported Material particles. Therefore, the vehicles F1, F2 and F3 appear according to the idealized one shown in FIG. 1 Model without differences in transit time, at the same mutual distance at the exit of the traffic section.

Ein Transportvorgang dieser Art entspricht der partiellen Differentialgleichung (1):

Figure 00040001
A transport process of this kind corresponds to the partial differential equation (1):
Figure 00040001

Bei konstanter Geschwindigkeit v wiederholt sich also am Ende des Förderbandes nach der Laufzeit:

Figure 00040002
exakt der Veglauf von q1, d.h. es gilt: q2(t) = q1(t-TL).At constant speed v, the following is repeated at the end of the conveyor belt after the running time:
Figure 00040002
the exact run of q1, ie: q2 (t) = q1 (tT L ) .

Dieser Sachverhalt wird dem Verhaften der Fahrzeuge jedoch nicht gerecht, da die Fahrzeuge, die in einem Zeitintervall gemessen werden, nicht dieselbe Geschwindigkeit besitzen und folglich, wie in Fig. 2 gezeigt, nach unterschiedlichen Laufzeiten den zweiten Messquerschnitt passieren.This fact does not do justice to the arrest of the vehicles, since the vehicles that are in measured at a time interval do not have the same speed and consequently as in FIG. 2 shown to pass the second measurement cross-section after different running times.

Der vorgesehene Einsatz eines Kalman-Filters erfordert zudem eine Beschreibung des Prozesses in Form einer gewöhnlichen und nicht einer partiellen Differentialgleichung, welche partielle Ableitungen einer unbekannten Funktion mit mehreren Variablen enthält (bezüglich dem Unterschied von gewöhnlichen und partiellen Differentialgleichungen siehe [3], Seite 435).The intended use of a Kalman filter also requires a description of the process in Form of an ordinary and not a partial differential equation, which partial derivatives an unknown function with multiple variables (in relation to the difference from ordinary and partial differential equations see [3], page 435).

Man gelangt erfindungsgemäss zu einer gewöhnlichen Differentialgleichung, indem man nur noch diskrete Orte in dem Verkehrsabschnitt betrachtet. Der Verkehrsabschnitt wird hierzu in n Segmente s1, s2, ..., sn unterteilt, deren innere Grenzen nicht über Sensoren bzw. Messfühler verfügen. In jedem dieser Segmente wird die Geschwindigkeit als konstant angenommen. According to the invention, an ordinary differential equation is obtained by only considering discrete locations in the traffic section. For this purpose, the traffic section is divided into n segments s1, s2, ..., s n , the inner limits of which do not have sensors or sensors. The speed is assumed to be constant in each of these segments.

Beschreibt man das Verhalten eines der n Untersegmente durch eine diskrete Verzögerung 1. Ordnung, so verhält sich das gesamte Segment wie ein Laufzeitelement le, dessen Zeitkonstante wie folgt lautet:

Figure 00050001
If the behavior of one of the n subsegments is described by a discrete 1st order delay, the entire segment behaves like a runtime element le, the time constant of which is as follows:
Figure 00050001

Je grösser die Zahl n der Segmente ist, desto genauer lässt sich das Verhalten beschreiben. Der Transportvorgang wird somit abgebildet, indem n Verzögerungen 1. Ordnung hintereinander geschaltet werden:

Figure 00050002
The larger the number n of segments, the more precisely the behavior can be described. The transport process is thus mapped by connecting n first order delays:
Figure 00050002

Durch diese Massnahme erfolgt für gewählte Segmente s1, s2, ... (siehe Fig. 1) eine Ortsdiskretisierung, wodurch eine gewöhnliche Differentialgleichung erreicht wird. Der erste Wert der Kelle q1, 0 ist der Wert q1 am Eingang q1; der letzte Wert der Keife q1, n ist identisch mit dem Wert q2 am Ausgang des Modells.This measure results in a location discretization for selected segments s1, s2, ... (see FIG. 1), whereby an ordinary differential equation is achieved. The first value of the trowel q 1, 0 is the value q1 at the input q1; the last value of the q 1, n is identical to the q2 value at the exit of the model.

Fasst man die Werte q1, i zu einem den inneren Zustand des Systems beschreibenden Zustandsvektor X zusammen, so lässt sich das Differentialgleichungssystem in Zustandsform angeben. Ein System, das aus vier Segmenten besteht, besitzt mit Tie = 4 / TL folgende Darstellung:

Figure 00050003
If the values q 1, i are combined to form a state vector X describing the internal state of the system, the system of differential equations can be stated in the form of a state. A system consisting of four segments also has T ie = 4 / T L following representation:
Figure 00050003

Die einzelnen Matrizen lauten wie folgt:

Figure 00050004
The individual matrices are as follows:
Figure 00050004

In Fig. 3 ist die Verkehrsstärke q2 als Reaktion auf einen synthetischen Verlauf von q1 für n=4 Segmente (bzw. Verzögerungen) dargestellt. Man erkennt, dass nun ein Teil der Fahrzeuge bereits vor Ablauf der Zeit TL eintrifft, während ein anderer Teil aufgrund einer etwas tieferen Geschwindigkeit mehr Zeit benötigt (siehe Fig. 2). 3 shows the traffic volume q2 in response to a synthetic course of q1 for n = 4 segments (or delays). It can be seen that some of the vehicles now arrive before the time T L has elapsed, while another part requires more time due to a somewhat lower speed (see FIG. 2).

Das gewählte Modell ist nun zur Beschreibung des Transportvorganges und damit des ungestörten Verkehrs in einem Verkehrsabschnitt geeignet. Diese Form ist zudem für ein Kalman-Filter sehr geeignet, da es linear in den Zustandsgrössen ist. Das Modell dient als Baustein für einen Verkehrsabschnitt mit zwei Messstellen. Mit Hilfe dieses Bausteins lassen sich leicht Modelle für beliebig komplexe Topologien aufstellen, die Einfahrten, Ausfahrten, Verzweigungen und Zusammenführungen aufweisen.The selected model is now used to describe the transport process and thus the undisturbed traffic suitable in a traffic segment. This shape is also very suitable for a Kalman filter because it is linear in the state variables. The model serves as a building block for a traffic section with two Measuring points. With the help of this module, models for arbitrarily complex topologies can be easily set up, which have entrances, exits, branches and junctions.

Die Messwerte liegen als zeitdiskrete Werte für ein bestimmtes Messintervall (Abtastzeit Ts) vor. Daher wird die diskrete Version des Modells gewählt, die sich aus dem kontinuierlichen Modell wie folgt ergibt (die Matrix C ändert sich nicht):

Figure 00060001
The measured values are available as discrete-time values for a specific measuring interval (sampling time Ts). Therefore, the discrete version of the model is chosen, which results from the continuous model as follows (the matrix C does not change):
Figure 00060001

Die unendliche Reihe konvergiert aufgrund der Fakultät im Nenner sehr schnell, so dass keine sehr hohe Rechenleistung notwendig ist. Allerdings muss die Diskretisierung in jedem Messintervall k neu durchgeführt werden, da die Matrizen

Figure 00060002
von der gemessenen Geschwindigkeit abhängen und somit nicht konstant sind. Die Diskretisierung kontinuierlicher Systeme ist u.a. auch aus [4], Seiten 74-80, Kapitel 3.1 (Diskretisierung der Regelstrecke) bekannt. Für das diskrete Modell gehen somit die Gleichungen : x(k + 1) = A(k) xe (k) + B(k) q1(k), q2 m (k) = C 1 x(k) + D(k) q1(k). The infinite series converges very quickly due to the factorial in the denominator, so that very high computing power is not necessary. However, the discretization must be carried out anew in every measuring interval k, since the matrices
Figure 00060002
depend on the measured speed and are therefore not constant. The discretization of continuous systems is also known from [4], pages 74-80, chapter 3.1 (discretization of the control system). The equations for the discrete model are therefore: x ( k + 1) = A ( k ) x e ( k ) + B ( k ) q 1( k ), q 2nd m ( k ) = C. 1 x ( k ) + D ( k ) q 1( k ).

Vereinfachend wurde der Abtastschritt k als Index geschrieben. Der Index k kennzeichnet für die Zeit t das Intervall kTS ≤ t < (k + 1) TS To simplify matters, the scanning step k was written as an index. The index k denotes the interval for the time t kT S ≤ t <( k + 1) T S

Erfindungsgemäss werden in einem Intervall alle Werte als konstant angenommen, wodurch eine einfache Berechnung möglich wird. Es ist zu beachten, dass die Parameter des Prozessmodells jedoch nicht konstant sind. Nach jedem Abtastintervall TS ergibt sich ein neuer Wert für die am Eingang des Verkehrsabschnittes gemessene Geschwindigkeit v1, mit der sich, gemäss Annahme, alle Fahrzeuge in dem überwachten Abschnitt fortbewegen.According to the invention, all values are assumed to be constant in one interval, which makes simple calculation possible. It should be noted that the parameters of the process model are not constant. After each sampling interval T S , there is a new value for the speed v1 measured at the entrance of the traffic section, with which, according to the assumption, all vehicles move in the monitored section.

Die Abweichungen zwischen den Vorhersagewerten des Modells q2m und den tatsächlich gemessenen Werten q2 werden vom Kalman-Filter, das die Schätzwerte q2e berechnet, im ungestörten Fall korrigiert. Falls jedoch Störungen auftreten, die bei der Modellierung nicht berücksichtigt wurden, liefert auch das Kalman-Filter fehlerhafte Werte. Der Störfall wird nun berücksichtigt, indem das Modell, wie in Fig. 4 gezeigt, um virtuelle Ein- und Ausfahrten zu einem "Parkplatz mit unendlicher Kapazität" erweitert wird. The deviations between the predicted values of the model q2 m and the actually measured values q2 are corrected by the Kalman filter, which calculates the estimated values q2 e , in the undisturbed case. However, if malfunctions occur that were not taken into account in the modeling, the Kalman filter also supplies incorrect values. The accident is now taken into account by extending the model, as shown in FIG. 4, by virtual entrances and exits to a “parking space with infinite capacity”.

Die Verkehrsströme auf den virtuellen Ausfahrten, die mit q2d bezeichnet sind, sind ein Mass für die Auswirkungen einer aufgetretenen Störung. Durch das Kalman-Filter werden die Werte q2d geschätzt (Wert q2de).The traffic flows on the virtual exits, which are designated q2 d , are a measure of the effects of a disruption that has occurred. The Kalman filter estimates the values q2 d (value q2 de ).

Da das Modell, wie eingangs erwähnt, eine unbeschränkte Kapazität aufweist, kann aus der Eingangsinformation q1 und dem Messwert q2 nicht auf den Wert q2d geschlossen werden. Da das Kalman-Filter optimale Schätzungen nur basierend auf einem Modell durchführen kann, muss auch der Wert q2d durch das Modell berücksichtigt werden. Da das Modell für den störungsfreien Betrieb ausgelegt ist, wird konsequenterweise angenommen, dass der Anfangswert von q2d gleich Null ist und sich anschliessend von Abtastperiode zu Abtastperiode wie folgt verhält: q2d (k+1) = q2d(k). Das Modell sieht daher, wie modelliert, keine Störung vor, so dass der Wert q2d konstant ist bzw. mit jedem Intervall k nur durch das Kalman-Filter KF verändert werden kann (q2d (k+1) = q2d(k) + corrKF(k)). Dazu wird der für vier Segmente s1, s2, s3, s4 vorgesehene Zustandsvektor x wie folgt erweitert:

Figure 00070001
Since the model has an unlimited capacity, as mentioned at the beginning, the value q2 d cannot be concluded from the input information q1 and the measured value q2. Since the Kalman filter can only make optimal estimates based on a model, the value q2 d must also be taken into account by the model. Since the model is designed for trouble-free operation, it is consequently assumed that the initial value of q2 d is zero and then behaves as follows from sampling period to sampling period: q2 d (k + 1) = q2 d (k) . As modeled, the model therefore does not provide for any disturbance, so that the value q2d is constant or can only be changed with the Kalman filter KF at every interval k ( q2 d (k + 1) = q2 d (k) + corr KF (k) ). For this purpose, the state vector x provided for four segments s1, s2, s3, s4 is expanded as follows:
Figure 00070001

Im ungestörten Fall treten keine virtuellen Verkehrsflüsse (siehe Fig. 4) auf. Kommt es zu Störungen, so kann das Kalman-Filter über die virtuellen Verkehrsflüsse Messwerte und Schätzwerte des Modells abgleichen. Der Vektor x(k) wird durch das Kalman-Filter KF korrigiert (siehe Fig. 5), wodurch der Vektor xe(k) mit korrigierten Werten x1e, ..., x4e sowie xde entsteht. Die korrigierten (internen) Zustandswerte x4e plus xde bzw. xde entsprechen somit den Schätzwerten q2e und q2de des Kalman-Filters KF. Die Werte q2de sind ein unmittelbares Mass für den Grad der im überwachten Abschnitt aufgetretenen Störung (siehe Fig. 7).In the undisturbed case, there are no virtual traffic flows (see FIG. 4). If malfunctions occur, the Kalman filter can use the virtual traffic flows to compare measured values and estimated values of the model. The vector x (k) is corrected by the Kalman filter KF (see FIG. 5), which results in the vector x e (k) with corrected values x 1e , ..., x 4e and x de . The corrected (internal) state values x 4e plus x de or x de thus correspond to the estimated values q2e and q2 de of the Kalman filter KF. The values q2 de are an immediate measure of the degree of disturbance that has occurred in the monitored section (see FIG. 7).

Grundlage des in Fig. 5 gezeigten Kalman-Filters KF ist das oben beschriebene, erweiterte diskrete Modell MOD, das den Verhältnissen im Idealzustand genügt. Zur Anpassung an die realen Verhältnisse wird zusätzlich angenommen, dass Prozess- und Sensorstörungen p bzw. s an den in Fig. 5 gezeigten Stellen wie folgt wirksam sind (Annahme D(k) =[0]): x(k + 1) = A(k)x(k) + B(k)q1(k) + G(k) p(k), q2* m (k) = C(k)x(k) + D(k)q1(k) + s(k) ≅ q2(k) The basis of the Kalman filter KF shown in FIG. 5 is the expanded discrete model MOD described above, which satisfies the conditions in the ideal state. To adapt to the real conditions, it is additionally assumed that process and sensor faults p and s at the points shown in FIG. 5 are effective as follows (assumption D (k) = [0]): x ( k + 1) = A ( k ) x ( k ) + B ( k ) q 1( k ) + G ( k ) p ( k ), q 2 * m ( k ) = C. ( k ) x ( k ) + D ( k ) q 1( k ) + s ( k ) ≅ q 2 ( k )

Der korrigierte Modellwert q2*m(k) entspricht somit dem im Intervall k tatsächlich gemessenen Wert q2(k). Die angenommenen Störungen resultieren aus dem Verkehrsprozess und entsprechen der Form eines frequenzunabhängigen Rauschens mit gaussverteiltem Amplitudenspektrum. Betrachtet man gemessene Verläufe der Verkehrsstärke q, so erkennt man, dass sich die starken Schwankungen recht gut als Rauschprozess um einen Mittelwert darstellen lassen. Das Kalman-Filter KF ist ein Optimalfilter, das die Varianz des Schätzfehlers xerr (k) ≅ x(k) - xe (k) minimiert. Der "wahre" Zustandsvektor ist nicht bekannt. Man verwendet deshalb xe(k) als Schätzwert für den unbekannten Zustandsvektor x(k).The corrected model value q2 * m (k) thus corresponds to the value q2 (k) actually measured in the interval k. The assumed disturbances result from the traffic process and correspond to the form of a frequency-independent noise with a Gaussian-distributed amplitude spectrum. If one looks at measured courses of traffic intensity q, one can see that the strong fluctuations can be represented quite well as a noise process around an average. The Kalman KF filter is an optimal filter that measures the variance of the estimation error x err ( k ) ≅ x ( k ) - x e ( k ) minimized. The "true" state vector is not known. Therefore x e (k) is used as an estimate for the unknown state vector x (k).

Die Matrix G(k), die genau wie die Matrizen A(k), B(k) durch Diskretisierung

Figure 00080001
entsteht, gibt an, wie sich die (unbekannten) Störgrössen auf die internen Zustandsgrössen des Prozesses verteilen. Im einfachsten Fall nimmt man an, dass diese Störgrössen gleichmässig auf alle internen Zustandsgrössen einwirken. Durch die Matrix G(k) wird daher angegeben wo nicht aber mit welcher Intensität die Störgrössen auf das Modell einwirken. Die Intensität der Störungen wird durch die nachstehend beschriebenen Kovarianzmatrizen Q und R beschrieben. Sofern zum Beispiel drei Zustandsgrössen und zwei Störgrössen vorhanden sind und die erste Störgrösse nur auf die erste Zustandsgrösse und die zweite Störgrösse nur auf die beiden folgenden Zustandsgrössen einwirkt, würde dieser Sachverhalt mit einer Matrix G(k) der nachstehend angegebenen Form berücksichtigt:
Figure 00080002
The matrix G (k), which is exactly like the matrices A (k), B (k) by discretization
Figure 00080001
arises, indicates how the (unknown) disturbance variables are distributed among the internal state variables of the process. In the simplest case, it is assumed that these disturbance variables have an even effect on all internal state variables. The matrix G (k) therefore indicates where, but not with what intensity, the disturbance variables affect the model. The intensity of the disturbances is described by the covariance matrices Q and R described below. If, for example, there are three state variables and two disturbance variables and the first disturbance variable only affects the first state variable and the second disturbance variable only affects the following two state variables, this would be taken into account with a matrix G (k) of the form given below:
Figure 00080002

Erfindungsgemäss wird die Matrix G(k) derart gewählt, dass Prozessstörungen p alle Grössen des Zustandsvektors x(k) in gleicher Weise beeinflussen.According to the invention, the matrix G (k) is selected in such a way that process disturbances p all sizes of the state vector influence x (k) in the same way.

Die Steuerparameter des Kalman-Filters KF sind die Kovarianzmatrizen Q und R der (angenommenen) Rauschprozesse, die auf den Prozess selbst bzw. auf die von den Sensoren abgegebenen Messwerte einwirken. Die Matrix Q beschreibt die Intensität der Prozessstörungen p. Die Matrix R beschreibt die Intensität der Sensorstörungen s. Obwohl die Matrizen Q und R theoretisch zeitvariant sein dürfen, wird angenommen, dass die Schwankungen der Messwerte unabhängig von der Verkehrssituation und somit zeitunabhängig sind. Die Matrizen Q und R nehmen damit konstante Werte an. Die Elemente der Matrizen Q und R sind die Entwurfsparameter des Kalman-Filters KF, mit denen die Einschwingzeit und Störempfindlichkeit des Kalman-Filters KF eingestellt werden.The control parameters of the Kalman filter KF are the covariance matrices Q and R of the (assumed) Noise processes that relate to the process itself or to the measured values output by the sensors act. The matrix Q describes the intensity of the process disturbances p. The matrix R describes the Intensity of sensor disturbances see Although the matrices Q and R may theoretically be time-variant, assumed that the fluctuations in the measured values regardless of the traffic situation and thus are independent of time. The matrices Q and R thus take on constant values. The elements of the matrices Q and R are the design parameters of the Kalman KF filter, with which the settling time and sensitivity to interference of the Kalman filter KF can be set.

Mit den Initialisierungswerten P(0) = G(0)QGT (0) und x(0) = 0 sind folgende Gleichungen zu berechnen. Man unterscheidet zwischen den Gleichungen zur Bestimmung der optimalen Schätzwerte und den Gleichungen zur Berechnung des Folgezustandes.With the initialization values P (0) = G (0) QG T (0) and x (0) = 0 the following equations must be calculated. A distinction is made between the equations for determining the optimal estimated values and the equations for calculating the subsequent state.

Die Gleichungen zur Bestimmung der optimalen Schätzwerte lauten wie folgt (I ist die Einheitsmatrix) :

Figure 00090001
The equations for determining the optimal estimates are as follows (I is the unit matrix):
Figure 00090001

Die Gleichungen zur Berechnung des Folgezustandes lauten wie folgt: x(k + 1) = A(k)xe (k) + B(k)q1(k), P(k + 1) = A(k)Pe (k)AT (k) + G(k)QGT (k) The equations for calculating the subsequent state are as follows: x ( k + 1) = A ( k ) x e ( k ) + B ( k ) q 1( k ), P ( k + 1) = A ( k ) P e ( k ) A T ( k ) + G ( k ) QG T ( k )

Die obengenannten Gleichungen werden für jeden Abtastschritt berechnet. Die Reihenfolge, in der die Gleichungen berechnet werden, darf nicht vertauscht werden, da die Ergebnisse teilweise voneinander abhängen.The above equations are calculated for each scan step. The order in which the Equations should not be mixed up because the results are partly different depend.

Der Signalfluss im erfindungsgemässen Kalman-Filter KF ist in Fig. 5 gezeigt. Für jedes Intervall erfolgt die Addition der Produkte aus Eingangswert q1(k) mal Matrix B(k) sowie geschätzter Zustandsvektor xe(k) mal Matrix A(k) in der Additionsstufe ADD1. Die Summe ergibt nach der zeitlichen Verschiebung den neuen unkorrigierten internen Zustandsvektor x(k) der mit der Matrix C1 = [00011] multipliziert wird, wodurch die Summe der Werte x4 plus xd gebildet wird (q2m = x4 + xd ). Anschliessend wird anhand der Differenzstufe DIFF die Differenz zwischen der vom Modell neu bestimmten Verkehrsstärke q2m(k) und der tatsächlich gemessenen Verkehrsstärke q2(k) (am Ausgang des Verkehrsabschnittes) gebildet. Die resultierende Differenz (q2m(k) - q2(k)) wird mit der Kalman-Matrix L(k) multipliziert, wodurch Werte gebildet werden, mit denen der Zustandsvektor x(k) anhand der Additionsstufe ADD3 korrigiert bzw. in den geschätzten Zustandsvektor xe(k) umgewandelt wird. Die geschätzten Werte x4e und xde werden anhand der Matrizen C1=[00011] und C2=[0000-1] ausgelesen. Anhand von C1 wird die Summe der geschätzten Werte x4e und xde gebildet, welche den Schätzwert q2e(k) ergibt. Wie bereits erwähnt, sind die Matrizen C1 und C2 konstant. Das letzte Element der Matrix C2 ist negativ, damit man bei Stau einen positiven Wert für q2de(k) erhält. Selbstverständlich können alle Rechenoperationen durch ein dem Fachmann bekanntes Rechnersystem (z.B.: Prozessor, Signalprozessor) durchgeführt werden.The signal flow in the Kalman filter KF according to the invention is shown in FIG. 5. For each interval, the products of the input value q1 (k) times matrix B (k) and the estimated state vector x e (k) times matrix A (k) are added in the addition stage ADD1. After the time shift, the sum results in the new uncorrected internal state vector x (k) which is multiplied by the matrix C 1 = [00011], whereby the sum of the values x 4 plus x d is formed ( q2 m = x 4th + x d ). The difference level DIFF is then used to form the difference between the traffic intensity q2 m (k) newly determined by the model and the actually measured traffic intensity q2 (k) (at the exit of the traffic segment). The resulting difference (q2 m (k) - q2 (k)) is multiplied by the Kalman matrix L (k), whereby values are formed with which the state vector x (k) is corrected using the addition stage ADD3 or in the estimated values State vector x e (k) is converted. The estimated values x4 e and x de are read out using the matrices C 1 = [00011] and C 2 = [0000-1]. Using C 1 , the sum of the estimated values x4 e and x de is formed, which gives the estimated value q2 e (k). As already mentioned, the matrices C 1 and C 2 are constant. The last element of the matrix C 2 is negative, so that a positive value for q2 de (k) is obtained in the event of a jam. Of course, all arithmetic operations can be carried out by a computer system known to the person skilled in the art (for example: processor, signal processor).

Fig. 7 zeigt den vom Kalman-Filter KF geschätzten Vertauf der Verkehrsstärke q2e sowie die geschätzte Verkehrsstärke q2de auf der virtuellen Ein- und Ausfahrt, beeinflusst durch eine Störung, die zu einem Zeitpunkt 11 auftritt und während der Zeit Tdist bis zu einem Zeitpunkt t2 anhält. Bis zum Zeitpunkt t1 können die Fahrzeuge den Verkehrsabschnitt ungehindert passieren, so dass die am Ein- und Ausgang festgestellten Verkehrsstärken unter Vernachlässigung der Laufzeitunterschiede dem Wert q1 entsprechen. Zum Zeitpunkt t1 tritt eine Störung auf, durch die innerhalb dem Verkehrsabschnitt ein Engpass entsteht, der den Verkehr nur noch mit der reduzierten Verkehrsstärke q2dist passieren lässt. Damit entsprechende Korrektur- und Sicherheitsmassnahmen rechtzeitig ergriffen werden können, muss die Störung in der Verkehrsleitzentrale unverzüglich erkannt werden. In Fig. 7 ist zusätzlich ein für die Messgrössen q1 und q2 typischer Verlauf angegeben. Dabei ist ersichtlich, dass die Messgrössen q1, q2 auch im ungestörten Fall grosse statistische Schwankungen aufweisen und dass die Auswirkung einer Störung innerhalb des Verkehrsabschnittes erst nach einer erheblichen Verzögerung feststellbar ist. Die erfindungsgemässe Vorrichtung und das Verfahren, durch die Ein- und Ausgangsgrössen q1, q2 anhand des Kalman-Filters KF verarbeitet werden, erlauben nun, das Auftreten der Störung praktisch verzögerungsfrei zu erfassen. Unmittelbar nach Auftreten der Störung läuft der Wert q2de hoch. Dies entspricht einem Anstieg der geschätzten Verkehrsstärke q2de auf der Zufahrt zum virtuellen Parkplatz. Nach Aufhebung der Störung zum Zeitpunkt t2 steigt die Verkehrsstärke q2de vom virtuellen Parkplatz zurück in den Verkehrsabschnitt wieder an und fällt anschliessend mit der Angleichung der Verkehrsstärken q1 und q2 am Ein- und Ausgang des Verkehrsabschnittes wieder auf Null zurück. Durch Vergleich der geschätzten Verkehrsstärken q2de zwischen dem virtuellen Parkplatz und dem überwachten Verkehrsabschnitt mit positiven und negativen Schwellwerten thp; thm kann daher jeweils schnell festgestellt werden, ob eine Störung aufgetreten ist oder aufgehoben wurde.FIG. 7 shows the traffic intensity q2 e estimated by the Kalman filter KF and the estimated traffic intensity q2 de on the virtual entrance and exit, influenced by a disturbance that occurs at a time 11 and during the time T dist up to one Time t2 continues. Up to time t1, the vehicles can pass through the traffic section unhindered, so that the traffic volumes determined at the entrance and exit correspond to the value q1, neglecting the travel time differences. At time t1, a disturbance occurs, which creates a bottleneck within the traffic section, which only allows traffic to pass with the reduced traffic volume q2 dist . In order for corrective and safety measures to be taken in good time, the fault in the traffic control center must be recognized immediately. FIG. 7 also shows a typical course for the measured variables q1 and q2. It can be seen that the measured variables q1, q2 have large statistical fluctuations even in the undisturbed case and that the effect of a disturbance within the traffic section can only be determined after a considerable delay. The device according to the invention and the method by means of which input and output variables q1, q2 are processed using the Kalman filter KF now allow the occurrence of the fault to be detected with practically no delay. The value q2 de starts up immediately after the fault occurs. This corresponds to an increase in the estimated traffic volume q2 de on the entrance to the virtual parking lot. After the disturbance has been eliminated at time t2, the traffic volume q2 de increases again from the virtual parking lot back into the traffic section and then falls back to zero when the traffic volumes q1 and q2 are adjusted at the entrance and exit of the traffic section. By comparing the estimated traffic volumes q2 de between the virtual parking lot and the monitored traffic section with positive and negative threshold values thp; Therefore, thm can quickly determine whether a fault has occurred or has been remedied.

Die Schwellwerte werden entsprechend der Grösse der zu erfassenden Störung gewählt. Vorzugsweise werden mehrere Schwellwerte vorgesehen, welche den Zuständen "leichte Verkehrsbehinderung",, "zähflüssiger Verkehr" oder "Stau" entsprechen. Beim Eintreten der entsprechenden Zustände können daher die geeigneten Massnahmen (z.B. Stauwarnung) initialisiert werden.The threshold values are chosen according to the size of the fault to be detected. Preferably several threshold values are provided which correspond to the states "light traffic obstruction" ,, "slow traffic" or "traffic jam". When the corresponding conditions occur therefore the appropriate measures (e.g. traffic jam warning) are initialized.

Zudem können verschiedene Strassenabschnitte miteinander verbunden werden. Die Modelle unterschiedlicher Strassenabschnitte lassen sich einfach aneinander reihen, wobei für jedes Modell wiederum Länge, Ausgangs- und Eingangsgrössen des betreffenden Verkehrsabschnittes zu berücksichtigen sind. Messwerte am Ausgang eines Verkehrsabschnittes können daher gleichzeitig als Eingangsgrössen für den anschliessenden Verkehrsabschnitt verwendet werden. Durch Verbindung mehrerer Verkehrsabschnitte kommt eine verzögerungsfreie Glättungswirkung zum tragen.In addition, different road sections can be connected to each other. The models different Street sections can easily be strung together, again for each model Length, output and input variables of the relevant traffic section must be taken into account. Measured values at the exit of a traffic segment can therefore also serve as input variables for the subsequent traffic section can be used. By connecting several traffic sections there is an instantaneous smoothing effect.

Anhand der gewonnen Kenntnisse über die Verkehrsstärken lassen sich ferner auch die zu erwartenden Reisezeiten schnell und präzise berechnen. Dadurch können für die Verkehrsteilnehmer auch optimale Reiserouten festgelegt werden.Based on the knowledge gained about the traffic levels, the expected ones can also be expected Calculate travel times quickly and precisely. This can also be optimal for road users Itineraries to be set.

Literaturverzeichnisbibliography

  • [1] Michael Cremer, Der Verkehrsfluss auf Schnellstrassen, Springer Verlag, Berlin 1979[1] Michael Cremer, The flow of traffic on expressways, Springer Verlag, Berlin 1979
  • [2] Kai Müller, Entwurf robuster Regelungen, Teubner Verlag, Stuttgart 1996[2] Kai Müller, draft robust regulations, Teubner Verlag, Stuttgart 1996
  • [3] Lothar Papula, Mathematik für Ingenieure und Naturwissenschaftler, Vieweg Verlag, Wiesbaden 1997, 8. Auflage[3] Lothar Papula, mathematics for engineers and scientists, Vieweg Verlag, Wiesbaden 1997, 8th edition
  • [4] Jürgen Ackermann, Abtastregelung, Springer Verlag, 3. Auflage, Berlin 1988[4] Jürgen Ackermann, scanning control, Springer Verlag, 3rd edition, Berlin 1988
  • Claims (10)

    Verfahren zur Ermittlung des Verkehrszustandes innerhalb eines Verkehrsabschnittes mit einem auf einem Modell (MOD) basierenden Kalman-Fifter (KF), durch das die Länge ls des Verkehrsabschnittes sowie die am Ein- und Ausgang des Verkehrsabschnittes gemessenen Verkehrsstärken q1 bzw. q2 berücksichtigt werden, dadurch gekennzeichnet, dass ein transportorientiertes Modell (MOD) gewählt wird, in dem für wenigstens zwei Segmente (s1; ...; sn) des Verkehrsabschnittes anhand der Länge ls sowie der am Eingang des Verkehrsabschnittes gemessenen Verkehrsstärke q1 und der Fahrzeuggeschwindigkeit v1(k) in Abtastintervallen k jeweils ein einen Korrekturwert (xd) sowie für jedes Segment (s1; ...; sn) einen Verkehrsstärkewert (x1; ...; xn) enthaltenden Vektor x(k) gebildet wird, der durch Vergleich mit der am Ausgang des Verkehrsabschnittes gemessenen Verkehrsstärke q2 anhand des Kalman-Filters (KF) korrigiert wird, so dass dem korrigierten Vektor xe(k) präzisierte Schätzwerte q2e für die Verkehrsstärke am Ausgang des Verkehrsabschnittes sowie ein Schätzwert q2de für eine den inneren Zustand des Verkehrsabschnittes repräsentierende Differenzverkehrsstärke q2d entnehmbar sind.Method for determining the traffic condition within a traffic segment using a Kalman Fifter (KF) based on a model (MOD), by means of which the length l s of the traffic segment and the traffic volumes q1 and q2 measured at the entrance and exit of the traffic segment are taken into account, characterized in that a transport-oriented model (MOD) is selected in which for at least two segments (s1; ...; s n ) of the traffic section based on the length l s and the traffic intensity q1 measured at the entrance of the traffic section and the vehicle speed v1 ( k) at sampling intervals k, a vector x (k) containing a correction value (x d ) and for each segment (s1; ...; s n ) a traffic intensity value (x1; ...; x n ) is formed, which is formed by Comparison with the traffic volume q2 measured at the exit of the traffic section is corrected using the Kalman filter (KF), so that the corrected vector x e (k) makes the estimated value more precise e q2 e for the traffic intensity at the exit of the traffic segment and an estimate q2 de for a differential traffic intensity q2 d representing the inner state of the traffic segment can be gathered. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass ein neuer Vektor x(k+1) anhand von in Abhängigkeit der Geschwindigkeit v1(k) für jedes Abtastintervall k neu festgelegten Matrizen A(k) und B(k) des erweiterten Verkehrsmodells (MOD) wie folgt berechnet wird: x(k + 1) = A(k) xe (k) + B(k) q1(k), wobei der alte Vektor x(k) durch das Kalman-Filter (KF) in den korrigierten Vektor xe(k) überführt wird:
    Figure 00110001
    indem die Differenz zwischen dem Messwert q2(k) und der Summe der Werte xn(k) und xd(k) mit einer für jedes Abtastintervall neu festgelegten Kalman-Matrix L(k) multipliziert wird.
    Method according to Claim 1, characterized in that a new vector x (k + 1) on the basis of matrices A (k) and B (k) of the extended traffic model (MOD) newly defined for each sampling interval k as a function of the speed v1 (k) is calculated as follows: x ( k + 1) = A ( k ) x e ( k ) + B ( k ) q 1( k ), whereby the old vector x (k) is converted into the corrected vector x e (k) by the Kalman filter (KF):
    Figure 00110001
    by multiplying the difference between the measured value q2 (k) and the sum of the values x n (k) and x d (k) by a Kalman matrix L (k) newly defined for each sampling interval.
    Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass die Schätzwerte q2de für die den inneren Zustand des Verkehrsabschnittes repräsentierenden Differenzverkehrsstärken q2d mit wenigstens einem Schwellwert thp; thm verglichen werden, nach dessen Überschreiten eine Zustandsänderung innerhalb des Verkehrsabschnittes festgestellt wird.Method according to Claim 1 or 2, characterized in that the estimated values q2 de for the differential traffic volumes q2 d representing the inner state of the traffic section have at least one threshold value thp; thm are compared, after which a change in state within the traffic section is determined. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass in Abhängigkeit der Polarität und Amplitude der Schätzwerte q2de schrittweise Reduktionen und/oder Steigerungen der Kapazität des Verkehrsabschnittes festgestellt und in Abhängigkeit davon Massnahmen zur Information der Verkehrsbeteiligten und/oder zur Steuerung des Verkehrs getroffen werden. Method according to Claim 3, characterized in that, depending on the polarity and amplitude of the estimated values q2, de gradual reductions and / or increases in the capacity of the traffic section are determined and, depending on this, measures are taken to inform the traffic participants and / or to control the traffic. Verfahren nach einem der Ansprüche 1-4, dadurch gekennzeichnet, dass mehrere Verkehrsabschnitte in einem Modell (MOD) verknüpft und von einem Kalman-Filter (KF) überwacht werden und/oder dass in Abhängigkeit der festgestellten Zustände Reisezeiten und/oder optimale Verkehrsrouten berechnet werden.Method according to one of claims 1-4, characterized in that several traffic sections are linked in a model (MOD) and monitored by a Kalman filter (KF) and / or that travel times and / or optimal traffic routes are calculated depending on the determined conditions . Vorrichtung zur Ermittlung des Verkehrszustandes innerhalb eines Verkehrsabschnittes mit einem auf einem Modell (MOD) basierenden Kalman-Filter (KF), durch das die Länge ls sowie die mittels Messfühlem am Ein- und Ausgang des Verkehrsabschnittes gemessenen Verkehrsstärken q1 bzw. q2 berücksichtigt werden, dadurch gekennzeichnet, dass ein transportorientiertes Modell (MOD) vorgesehen ist, in dem für wenigstens zwei Segmente (s1; ...; sn) des Verkehrsabschnittes anhand der Länge ls sowie der am Eingang des Verkehrsabschnittes gemessenen Verkehrsstärke q1 und der Fahrzeuggeschwindigkeit v1(k) in Abtastintervallen k jeweils ein einen Korrekturwert (xd) sowie für jedes Segment (s1; ...; sn) einen Verkehrsstärkewert (x1; ...; xn) enthaltenden Vektor x(k) erzeugbar ist, der durch Vergleich mit der am Ausgang des Verkehrsabschnittes gemessenen Verkehrsstärke q2 anhand des Kalman-Filters (KF) korrigierbar ist, so dass dem korrigierten Vektor xe(k) präzisierte Schätzwerte q2e für die Verkehrsstärke am Ausgang des Verkehrsabschnittes sowie ein Schätzwert q2de für eine den inneren Zustand des Verkehrsabschnittes repräsentierende Differenzverkehrsstärke q2d entnehmbar sind.Device for determining the traffic condition within a traffic segment with a Kalman filter (KF) based on a model (MOD), by which the length l s and the traffic volumes q1 and q2 measured by means of sensors at the entrance and exit of the traffic segment are taken into account, characterized in that a transport-oriented model (MOD) is provided, in which for at least two segments (s1; ...; s n ) of the traffic section based on the length l s and the traffic intensity q1 measured at the entrance of the traffic section and the vehicle speed v1 ( k) at sampling intervals k a vector x (k) containing a correction value (x d ) and for each segment (s1; ...; s n ) a traffic intensity value (x1; ...; x n ) can be generated Comparison with the traffic volume q2 measured at the exit of the traffic section can be corrected using the Kalman filter (KF), so that the corrected vector x e (k) makes the estimated value more precise e q2 e for the traffic intensity at the exit of the traffic segment and an estimate q2 de for a differential traffic intensity q2 d representing the inner state of the traffic segment can be gathered. Vorrichtung nach Anspruch 6, dadurch gekennzeichnet, dass ein neuer Vektor x(k+1) anhand von in Abhängigkeit der Geschwindigkeit v1(k) für jedes Abtastintervall k neu festgelegten Matrizen A(k) und B(k) des erweiterten Verkehrsmodells (MOD) wie folgt berechenbar ist: x(k + 1) = A(k) xe (k) + B(k) q1(k), wobei der alte Vektor x(k) durch das Kalman-Filter (KF) in den korrigierten Vektor xe(k) überführbar ist
    Figure 00120001
    indem die Differenz zwischen dem Messwert q2(k) und der Summe der Werte xn(k) und xd(k) mit einer für jedes Abtastintervall neu festgelegten Kalman-Matrix L(k) multiplizierbar ist.
    Apparatus according to claim 6, characterized in that a new vector x (k + 1) on the basis of matrices A (k) and B (k) of the extended traffic model (MOD) newly defined for each sampling interval k as a function of the speed v1 (k) can be calculated as follows: x ( k + 1) = A ( k ) x e ( k ) + B ( k ) q 1( k ), the old vector x (k) can be converted into the corrected vector x e (k) by the Kalman filter (KF)
    Figure 00120001
    in that the difference between the measured value q2 (k) and the sum of the values x n (k) and x d (k) can be multiplied by a Kalman matrix L (k) newly defined for each sampling interval.
    Vorrichtung nach Anspruch 6 oder 7, dadurch gekennzeichnet, dass die Schätzwerte q2de für die den inneren Zustand des Verkehrsabschnittes repräsentierenden Differenzverkehrsstärken q2d mit wenigstens einem Schwellwert thp; thm vergleichbar sind, nach dessen Überschreiten eine Zustandsänderung innerhalb des Verkehrsabschnittes feststellbar ist. Apparatus according to claim 6 or 7, characterized in that the estimated values q2 de for the differential traffic volumes q2 d representing the inner state of the traffic section with at least one threshold value thp; thm are comparable, after which a change in state within the traffic section can be determined. Vorrichtung nach Anspruch 8, dadurch gekennzeichnet, dass in Abhängigkeit der Polarität und Amplitude der Schätzwerte q2de schrittweise Reduktionen und/oder Steigerungen der Kapazität des Verkehrsabschnittes feststellbar und in Abhängigkeit davon Massnahmen zur Information der Verkehrsbeteiligten und/oder zur Steuerung des Verkehrs realisierbar sind.Apparatus according to claim 8, characterized in that depending on the polarity and amplitude of the estimated values q2 de gradual reductions and / or increases in the capacity of the traffic section can be determined and, depending on this, measures can be implemented to inform the traffic participants and / or to control the traffic. Vorrichtung nach einem der Ansprüche 6-9, dadurch gekennzeichnet, dass mehrere Verkehrsabschnitte in einem Modell (MOD) verknüpft und durch ein Kalman-Filter (KF) überwacht sind und/oder dass in Abhängigkeit der festgestellten Zustände Reisezeiten und/oder optimale Verkehrsrouten berechenbar sind.Device according to one of claims 6-9, characterized in that several traffic sections are linked in a model (MOD) and monitored by a Kalman filter (KF) and / or that travel times and / or optimal traffic routes can be calculated depending on the determined conditions .
    EP00108359A 1999-05-17 2000-04-17 Method and device for the determination of the traffic conditions of a roadsegment Expired - Lifetime EP1056063B1 (en)

    Applications Claiming Priority (2)

    Application Number Priority Date Filing Date Title
    CH93199 1999-05-17
    CH93199 1999-05-17

    Publications (2)

    Publication Number Publication Date
    EP1056063A1 true EP1056063A1 (en) 2000-11-29
    EP1056063B1 EP1056063B1 (en) 2007-06-06

    Family

    ID=4198504

    Family Applications (1)

    Application Number Title Priority Date Filing Date
    EP00108359A Expired - Lifetime EP1056063B1 (en) 1999-05-17 2000-04-17 Method and device for the determination of the traffic conditions of a roadsegment

    Country Status (3)

    Country Link
    EP (1) EP1056063B1 (en)
    AT (1) ATE364211T1 (en)
    DE (1) DE50014383D1 (en)

    Cited By (6)

    * Cited by examiner, † Cited by third party
    Publication number Priority date Publication date Assignee Title
    FR2822576A1 (en) * 2001-03-22 2002-09-27 Renault Method of monitoring of road traffic involves detecting obstacles adjacent to vehicle and producing map with binary trace of road condition data
    EP1280118A1 (en) * 2001-07-25 2003-01-29 Siemens Schweiz AG Method for determining road traffic states
    DE102004009898A1 (en) * 2004-02-26 2005-09-22 Siemens Ag Determining traffic situation on section of road network involves determining traffic situation from measurement values using analysis technique with which errors in measurement evaluation process are also detected
    DE102005043471A1 (en) * 2005-09-13 2007-03-15 Daimlerchrysler Ag Vehicle-sided traffic-adaptive assistance system controlling method for use in control device, involves evaluating two spatially and/or temporally sections of road from environment information and selecting parameters for controlling system
    DE102006032162A1 (en) * 2006-07-12 2007-11-29 Daimlerchrysler Ag Control of traffic-adaptive assistance system comprises measuring e.g. speed gradient for sections of road ahead of vehicle and associating with time constants for controlling system which decrease with proximity of stretch to vehicle
    CN113160555A (en) * 2021-03-01 2021-07-23 武汉理工大学 Road-state ripple processing method and system based on road side sensing equipment and storage medium

    Citations (1)

    * Cited by examiner, † Cited by third party
    Publication number Priority date Publication date Assignee Title
    US5801943A (en) * 1993-07-23 1998-09-01 Condition Monitoring Systems Traffic surveillance and simulation apparatus

    Patent Citations (1)

    * Cited by examiner, † Cited by third party
    Publication number Priority date Publication date Assignee Title
    US5801943A (en) * 1993-07-23 1998-09-01 Condition Monitoring Systems Traffic surveillance and simulation apparatus

    Cited By (8)

    * Cited by examiner, † Cited by third party
    Publication number Priority date Publication date Assignee Title
    FR2822576A1 (en) * 2001-03-22 2002-09-27 Renault Method of monitoring of road traffic involves detecting obstacles adjacent to vehicle and producing map with binary trace of road condition data
    EP1280118A1 (en) * 2001-07-25 2003-01-29 Siemens Schweiz AG Method for determining road traffic states
    DE102004009898A1 (en) * 2004-02-26 2005-09-22 Siemens Ag Determining traffic situation on section of road network involves determining traffic situation from measurement values using analysis technique with which errors in measurement evaluation process are also detected
    AT502459B1 (en) * 2004-02-26 2008-03-15 Siemens Ag METHOD FOR DETERMINING TRANSPORT CONDITION
    DE102004009898B4 (en) * 2004-02-26 2009-05-20 Siemens Ag A method for determining the traffic condition on a section of a road network and traffic management center for performing the method
    DE102005043471A1 (en) * 2005-09-13 2007-03-15 Daimlerchrysler Ag Vehicle-sided traffic-adaptive assistance system controlling method for use in control device, involves evaluating two spatially and/or temporally sections of road from environment information and selecting parameters for controlling system
    DE102006032162A1 (en) * 2006-07-12 2007-11-29 Daimlerchrysler Ag Control of traffic-adaptive assistance system comprises measuring e.g. speed gradient for sections of road ahead of vehicle and associating with time constants for controlling system which decrease with proximity of stretch to vehicle
    CN113160555A (en) * 2021-03-01 2021-07-23 武汉理工大学 Road-state ripple processing method and system based on road side sensing equipment and storage medium

    Also Published As

    Publication number Publication date
    ATE364211T1 (en) 2007-06-15
    DE50014383D1 (en) 2007-07-19
    EP1056063B1 (en) 2007-06-06

    Similar Documents

    Publication Publication Date Title
    DE69734474T2 (en) Vehicle detection system and signal processing system therefor
    EP2280383B1 (en) Method for determining traffic information for a section of a road network and traffic calculator to implement the method
    DE10133945A1 (en) Method and device for exchanging and processing data
    DE1798182A1 (en) Method and device for measuring the amount of grain conveyed per unit of time by means of a flow device
    EP2044406B1 (en) Method for the diagnosis of a blockage of an impulse line in a pressure measurement transducer, and pressure measurement transducer
    DE19647127A1 (en) Road traffic jam monitoring and analysis method
    DE102008003039A1 (en) Vehicle traffic condition recognizing method for traffic assistance system, involves differentiating traffic jam conditions on basis of measured speed, and recognizing traffic jam condition if another traffic jam condition is not recognized
    EP1136826A2 (en) Device for speed measurement
    EP3068643A1 (en) Method for controlling an actuator
    EP1056063B1 (en) Method and device for the determination of the traffic conditions of a roadsegment
    DE102008021380B4 (en) Method and device for predicting a course of a roadway and driver assistance system
    DE102012204123A1 (en) Method for determining queue length of traffic objects and/or vehicle groups, involves determining numerically identified traffic density by interpolation and core approximations at specific order, with extreme and/or turning points
    EP0501193B1 (en) Method for the automatic coordination of an independent intersection traffic light control device with one or more neighbouring intersections
    EP0821334A1 (en) Traffic control method and device
    EP1276085A1 (en) Method for determining a traffic jam index and for determining tailback lengths
    AT502459B1 (en) METHOD FOR DETERMINING TRANSPORT CONDITION
    EP0035960A1 (en) Method and device for the distinction of vehicles of different constructions in road traffic
    DE102019209255A1 (en) Determining a road condition
    DE102019218111A1 (en) Method for generating location information on relevant road sections for the detection of wrong-way trips and method for detecting wrong-way trips
    DE4228934C2 (en) Device for determining the confidence interval of percentile measured values of continuous stochastic sound signals
    DE102022201199B3 (en) Method and device for monitoring a lane detection system, lane detection system
    EP2550540A1 (en) Method for detecting magnetically marked objects and corresponding device
    EP3782868A2 (en) Method for calibrating a speed sensor of a railway vehicle
    DE102022111179A1 (en) METHOD AND DEVICE FOR GENERATING AN ARTIFICIAL INTELLIGENCE-BASED PREDICTOR AND USE THEREOF, AND COMPUTER PROGRAM
    EP3712786A1 (en) Method for evaluating at least one signal

    Legal Events

    Date Code Title Description
    PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

    Free format text: ORIGINAL CODE: 0009012

    AK Designated contracting states

    Kind code of ref document: A1

    Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LI LU MC NL PT SE

    AX Request for extension of the european patent

    Free format text: AL;LT;LV;MK;RO;SI

    17P Request for examination filed

    Effective date: 20001218

    AKX Designation fees paid

    Free format text: AT BE CH CY DE DK ES FI FR GB GR IE IT LI LU MC NL PT SE

    GRAP Despatch of communication of intention to grant a patent

    Free format text: ORIGINAL CODE: EPIDOSNIGR1

    GRAS Grant fee paid

    Free format text: ORIGINAL CODE: EPIDOSNIGR3

    GRAA (expected) grant

    Free format text: ORIGINAL CODE: 0009210

    AK Designated contracting states

    Kind code of ref document: B1

    Designated state(s): AT BE CH CY DE DK ES FI FR GB GR IE IT LI LU MC NL PT SE

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: FI

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070606

    REG Reference to a national code

    Ref country code: GB

    Ref legal event code: FG4D

    Free format text: NOT ENGLISH

    REG Reference to a national code

    Ref country code: CH

    Ref legal event code: EP

    REG Reference to a national code

    Ref country code: IE

    Ref legal event code: FG4D

    Free format text: LANGUAGE OF EP DOCUMENT: GERMAN

    REF Corresponds to:

    Ref document number: 50014383

    Country of ref document: DE

    Date of ref document: 20070719

    Kind code of ref document: P

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: SE

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070906

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: ES

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070917

    ET Fr: translation filed
    GBV Gb: ep patent (uk) treated as always having been void in accordance with gb section 77(7)/1977 [no translation filed]

    Effective date: 20070606

    REG Reference to a national code

    Ref country code: IE

    Ref legal event code: FD4D

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: PT

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20071106

    Ref country code: IE

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070606

    PLBE No opposition filed within time limit

    Free format text: ORIGINAL CODE: 0009261

    STAA Information on the status of an ep patent application or granted ep patent

    Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: GR

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070907

    Ref country code: DK

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070606

    Ref country code: GB

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070606

    26N No opposition filed

    Effective date: 20080307

    PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

    Ref country code: AT

    Payment date: 20080312

    Year of fee payment: 9

    PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

    Ref country code: IT

    Payment date: 20080424

    Year of fee payment: 9

    Ref country code: BE

    Payment date: 20080422

    Year of fee payment: 9

    PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

    Ref country code: CH

    Payment date: 20080708

    Year of fee payment: 9

    Ref country code: DE

    Payment date: 20080620

    Year of fee payment: 9

    Ref country code: NL

    Payment date: 20080417

    Year of fee payment: 9

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: MC

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20080430

    REG Reference to a national code

    Ref country code: CH

    Ref legal event code: PCOW

    Free format text: SIEMENS SCHWEIZ AG;FREILAGERSTRASSE 40;8047 ZUERICH (CH)

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: CY

    Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

    Effective date: 20070606

    BERE Be: lapsed

    Owner name: SIEMENS SCHWEIZ A.G.

    Effective date: 20090430

    REG Reference to a national code

    Ref country code: CH

    Ref legal event code: PL

    NLV4 Nl: lapsed or anulled due to non-payment of the annual fee

    Effective date: 20091101

    REG Reference to a national code

    Ref country code: FR

    Ref legal event code: ST

    Effective date: 20091231

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: AT

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20090417

    Ref country code: DE

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20091103

    Ref country code: CH

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20090430

    Ref country code: LI

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20090430

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: NL

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20091101

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: FR

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20091222

    PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

    Ref country code: FR

    Payment date: 20080422

    Year of fee payment: 9

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: BE

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20090430

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: LU

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20080417

    PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

    Ref country code: IT

    Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

    Effective date: 20090417