METHOD FOR SYNCHRONIZING SIGNALS

Abstract

A method for synchronizing signals of a plurality of participants. A relationship between the signals is given by a physical relation. The signals are each filtered with a first filter in order to determine a shift between the signals. The determined shift is a measure of the phase shift between the signals. The shift is subsequently eliminated by filtering the signals respectively with a second filter.

Claims

1-13. (canceled)

14. A method for synchronizing signals of at least one participant, a relationship between the signals being given by a mathematical relation, the method comprising the following steps: filtering each of the signals with a first filter to determine a shift between the signals, the determined shift being a measure of the phase shift between the signals; and eliminating the determined shift by filtering the signals respectively with a second filter.

15. The method as recited in claim 14, wherein a symmetry of the first filter and the second filter is used to determine and to eliminate both a positive shift and a negative shift.

16. The method as recited in claim 14, wherein the first filter is a Kalman filter, a measurement equation being realized as a first-order fractional delay filter.

17. The method as recited in claim 14, wherein the first-order fractional delay filter is a Lagrange filter.

18. The method as recited claim 14, wherein the second filter is a first-order or higher-order fractional delay filter.

19. The method as recited in claim 14, wherein a Kalman filter for estimating the shift and a Lagrange filter for synchronizing the signals are combined in one filter.

20. The method as recited in claim 19, wherein after the synchronization of the signals, an estimation of at least one parameter is carried based on the synchronized signals.

21. The method as recited in claim 20, wherein the at least one parameter is estimated with a Kalman filter.

22. The method as recited in claim 14, wherein the method is carried out in an onboard electrical network of a vehicle.

23. The method as recited in claim 14, wherein the signals have different sampling rates, and the signals are treated with a resampling filter.

24. A system for synchronizing signals of a plurality of participants, a relationship between the signals being given by a mathematical relation, the system configured to: filter each of the signals with a first filter to determine a shift between the signals, the determined shift being a measure of the phase shift between the signals; and eliminate the determined shift by filtering the signals respectively with a second filter.

25. The system as recited in claim 24, wherein the system is a control device of a vehicle.

26. A non-transitory machine-readable storage medium on which is stored a computer program for synchronizing signals of at least one participant, a relationship between the signals being given by a mathematical relation, the computer program, when executed by a computer, causing the computer to perform the following steps: filtering each of the signals with a first filter to determine a shift between the signals, the determined shift being a measure of the phase shift between the signals; and eliminating the determined shift by filtering the signals respectively with a second filter.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0024] FIG. 1 shows, in a diagram, a method for estimating parameters in accordance with an example embodiment of the present invention.

[0025] FIG. 2 shows an example embodiment of the method of the present invention.

[0026] FIG. 3 shows a graph illustrating the filter characteristic.

[0027] FIG. 4 shows a modeling of the paths before the filter.

[0028] FIG. 5 shows, in a graph, the principle of functioning of a synchronization.

[0029] FIG. 6 shows, in a schematic representation, a system for carrying out the method in accordance with an example embodiment of the present invention.

DETAILED DESCRIPTION OF EXAMPLE EMBODIMENTS

[0030] In the following, the problem with the asynchronicity of two signals is explained for the case of application of an onboard electrical network of a motor vehicle. It is to be kept in mind that the presented method is not limited to this case of application, but rather can be used whenever two asynchronous signals are to be synchronized.

[0031] FIG. 1 shows, in a diagram, a method for estimating parameters. The illustration shows a first block 50 for time updating or prediction, having a state prediction 52 and a prediction 54 of the error covariance, as well as a second block 60 for updating the measurement, or correction, having a calculation 62 of the Kalman gain, an estimation updating 64 with measurement and an updating 66 of the error covariance. An initial estimation is provided at input 70.

[0032] FIG. 2 shows, in a diagram, a possible embodiment of the presented method. The representation shows a physical system 10, a first Kalman filter 12 for estimating the delay, a fractional delay filter 14 for synchronizing the input and measurement variables, and a second Kalman filter 16 for estimating the parameters of physical system 10. Measurement variables z(k) 20 and input variables u(k) 22 are coupled to one another via physical system 10. Physical system 10 is given for example by the equation U=R*I, where the voltage U is the measured variable z(k) and the current I is the input variable u(k). The physical coupling of these two variables is provided via the resistor R.

[0033] Input variables u(k) go into physical system 10, first Kalman filter 12, and fractional delay filter 14. First Kalman filter 12 outputs a shift, or delay, D 24. Fractional delay filter 14 outputs u(k+D/2*t.sub.s), z(k−D/2*t.sub.s) 26, which represent the synchronized input and measured variables. Second Kalman filter 16 outputs estimated parameters 28 of physical system 12.

[0034] As already stated, in an onboard electrical or energy network of a vehicle diagnostic results are strongly influenced by the asynchronicity of the measured variables. If the asynchronicity can be determined and eliminated, this has a positive effect on the result of the diagnosis. The diagnostic design is based on time-discrete measured variables that have the same sampling rate but have a temporal shift relative to one another, the temporal shift being D*t.sub.s. Here, D is the delay factor between the signals as a linear factor relative to the sampling time, and t.sub.s is the sampling rate.

[0035] Therefore, under the assumption that the time shift is identical for both measured variables of a particular participant, of the named measured variables the following measured variables are available at sampling time k:

U.sub.v(k), I.sub.v(k), U.sub.Batt(k+D*t.sub.s), I.sub.Batt(k+D*t.sub.s)

[0036] For the determination of delay factor D, two designs can be distinguished. In the first design, for all signals of one of the participants a filtering is carried out that carries out a time shift by an initial factor D. The disadvantage of this design, compared to the following second design, is that only positive delays can be estimated and synchronized, and the signals of only one participant undergo the filter damping. Thus, it has to be known ahead of time which of the signals is advanced in order to make it possible to use the design.

[0037] According to the second design, for the determination of delay factor D the input and measured variables u(k), z(k) of all participants are filtered, as shown in FIG. 4 and FIG. 5. Here, for the signals of participant 1 a time shift of N/2+D/2 is carried out, and for the signals of participant 2 a time shift of N/2−D/2 is carried out. The difference of the time shifts of the two participants thus yields the overall shift factor D. The time shift of the signals can be carried out using a so-called fractional delay filter.

[0038] In the following, the realization of the fractional delay filter is described by a first-order Lagrange filter (N=1). Higher-order Lagrange filters (N>1) are also possible, but then the filter coefficients are different.

u(k−D*t.sub.s)=(u(k)*(1−D)+u(k−1)*D)

[0039] In this way, the estimated measurement variable z(k−(1/2−D/2) t.sub.s) and the estimated measured variable h (k−(1/2+D/2) t.sub.s) are calculated according to the procedure described above.

[0040] Here:

[0041] k: sampling time

[0042] z: measured variable

[0043] h: measured variable estimated from u via the model equation

[0044] u: input variable (measured variable, not synchronous with measured variable z)

[0045] For the parameter estimation of the supply line and contact resistances of the cable harness, the already-calculated values for z(k−(1/2−D/2) t.sub.s) and u(k−(1/2+D/2) t.sub.s) can then be used. Alternatively, as shown in FIG. 2, using the calculated factor D 24 an additional filtering 14 having higher quality can be realized, for example a higher-order Lagrange filter, whose filtered variables are then provided to the parameter estimation. This has the advantage that, with a low computing expense for the estimation of factor D, a high signal quality of the delayed signals can be achieved. The method can also be set up using only Kalman filter 12, in which case the advantages named above are not present.

[0046] The synchronization is thus made up of two components: Kalman filter 12 for the estimation of the shift, and fractional delay filter 14 for the synchronization of the signals. The shift D is a linear factor that describes the time shift between the signals as D*t.sub.s, where t.sub.s is the sample rate. The synchronized signals are then used by Kalman filter 16 to estimate the parameters, for example of a resistor.

[0047] For the synchronization of the signals of two subscribers, the relative shift D*t.sub.s to one another is the critical factor. In order to determine this shift and to enable the named advantages to be realized, a fractional delay filter, for example an n-order Lagrange filter, is newly modeled. On this, see FIG. 3.

[0048] FIG. 3 shows, in a graph 100 on whose abscissa 102 delay D is plotted and on whose ordinate 104 the damping |A| is plotted, the characteristic of a first-order filter, i.e. N=1, on the basis of a curve 106. A double arrow 110 indicates a ΔD of 0.4. A first arrow 112 indicates z(k−(0.5−D/2). A second arrow 114 indicates h(k−(0.5+D/2). A third arrow 116 indicates z(k−(0.5−D/2). A fourth arrow 118 indicates h(k−(0.5+D/2).

[0049] A D of 0.4 means that the signals of the first participant are advanced ahead of the signals of the second participant by 0.4*t.sub.s, where t.sub.s is equal to the sample time of the signals. A D of −0.4 means that the signals of the first participant run behind the signals of the second participant by 0.4*t.sub.s.

[0050] In the new modeling, the symmetry of the filter is used, which means that the damping |A| is equal for a delay of D and N−D. N corresponds here to the order of the filter. From this new modeling and the simultaneous filtering of the signals of the two participants, there result the two advantages, namely that an equal damping is achieved, and that both positive and negative delays are possible.

[0051] In the following, the filtering of the measured variables z(k) and of the estimated measured variables, calculated via the filtered input variables u(k) and the model equation h(k), are shown for the example of a first-order Lagrange filter (N=1):

z(k−(1/2−D/2) t.sub.s)=z(k)(1−(1/2−D/2))+z(k−1)(1/2−D/2)

h(k−(1/2+D/2) t.sub.s)=h(k)(1−(1/2+D/2))+h(k−1)(1/2+D/2)

[0052] If this measurement equation is implemented in a parameter estimator, then, using the noisy signals u(k) and z(k) and the model equation h(k), which is stated below, the delay D can be estimated.

D(k)=D(k−1)+[z(k−(1/2−D/2))−h(k−(1/2+D/2)t.sub.s)]*[h(k−1−(1/2+D/2) t.sub.s)−h(k−(1/2+D/2) t.sub.s)].sup.−1

[0053] FIG. 4 shows, in a diagram, the modeling of the paths before the fractional delay filter (reference character 14 in FIG. 2). The representation shows an addition unit 150, a subtraction unit 152, a first fractional delay filter 154, and a second fractional delay filter 156. At the input of addition unit 150 there are provided the value N/2 160 and D/2 162, which results from a multiplication of the delay D 164 by a factor of 0.5 166. D/2 162 and a value N/2 168 are also provided at the input of subtraction unit 152.

[0054] Thus, the addition unit outputs N/2+D/2 170. Subtraction unit 152 outputs N/2−D/2 172. h(k) 180 and N/2+D/2 170 are inputted to first fractional delay filter 154. z(k) 182 and N/2−D/2 172 are inputted to second fractional delay filter 156.

[0055] First fractional delay filter 154 outputs:

h(k−(N/2+D/2)*t.sub.s)

[0056] Second fractional delay filter 156 outputs:

z(k−(N/2−D/2)*t.sub.s)

[0057] The modeling of the shifts before fractional delay filters 154, 156 makes it possible to map positive and negative shifts between the components. Through the symmetry of fractional delay filters 154, 156, all signals undergo the same filter damping.

[0058] FIG. 5 shows, in a graph 300 on whose abscissa 302 k is plotted and on whose ordinate 304 the voltage U [V] is plotted, the principle of functioning of a synchronization as interpolation between two measured variables, represented by a first signal 310 and a second signal 312. These signals 310, 312 are asynchronous to one another. The relations of dependence of the measurements are given by a mathematical relationship, in this case

U.sub.1=U.sub.2+I.sub.2*R,

z(k) h(k)

[0059] where the left side of the equation is represented by first signal 310 and the right part of the equation is represented by second signal 312. The filtering now means an interpolation between measurement points that are identified by points in the Figure. The interpolations are illustrated by the straight lines in the Figure. Deviations in the mathematical relationship are then corrected by shifting the measurements. The interpolation of the measured variables of the two components, namely U.sub.1, U.sub.2, I.sub.2, enables the signal shift. Here there take place a shift of the interpolated measurement variables (U.sub.1) of the first component by the factor D/2+N/2 and a shift of the interpolated measurement variables (U.sub.2, I.sub.2) of the second component by the factor N/2−D/2.

[0060] FIG. 6 shows, in a schematic, highly simplified representation, a system for carrying out the method, designated as a whole by reference character 200. In this case, this system 200 is fashioned as a control device of a vehicle.

[0061] System 200 is connected to a first participant, or control device 202, and to a second participant, or control device 204, first control device 202 sending a first signal 206 to system 200, and second control device 204 sending a second signal 208 to system 200. The two signals 206, 208, which each carry measured values as information, are to be combined for evaluation in system 200; here it is to be kept in mind that the two signals 206, 208 are asynchronous to one another. In system 200, using the method presented herein a synchronization of the two signals 206, 208 can now be carried out, so that subsequently an evaluation of signals 206, 208 can be carried out, in this embodiment also in system 200.

[0062] Of course, the method can also be carried out with more than two participants or control devices. Here, the control devices can be synchronized with one another. However, a synchronization can also be carried out between one or more control devices and system 200.

[0063] The method can be applied in many ways if the following requirements are met:

[0064] The signals of the participants must have a mathematical relation.

[0065] The signals of the participants must have the same sampling rate, which can be resolved by a resampling filter.