Disclosure of Invention
The invention aims to provide an axle weight estimation method of a vehicle dynamic weighing system, which can reduce calculation and requirements on a hardware system and improve weighing precision.
The technical scheme for achieving the aim of the invention is as follows: an axle weight estimation method of a vehicle dynamic weighing system is characterized in that:
(1) when the vehicle passes through the dynamic axle weight scale, the weighing sensor collects original axle weight signals of each axle, the original axle weight signals comprise static axle weight signals, low-frequency interference signals and high-frequency interference signals, and the number of sampling points is counted;
(2) filtering out high-frequency interference signals with the frequency being more than or equal to 50 Hz;
(3) calculating the speed of the axle according to the number of sampling points, the set sampling frequency and the stroke of the vehicle passing through the dynamic axle weight scale;
(4) intercepting the effective sampling points when the vehicle passes through the dynamic axle weight scale according to the axle speed to obtain an effective axle weight signal;
(5) extracting 10-150 data points on the basis of the effective axle weight signal to form an axle weight signal for simplifying calculation;
(6) obtaining parameters required by nonlinear fitting by using the axle weight signal, and obtaining an approximate low-frequency interference signal by using a nonlinear fitting algorithm;
(7) and subtracting the approximate low-frequency interference signal from the effective axle weight signal or the axle weight signal to obtain an estimated static axle weight signal.
The invention adopts the technical scheme and has the advantages that:
1. according to the invention, the speed of each axle when the vehicle passes through the dynamic axle weight scale is dynamically estimated according to the sampling points of the collected signals, and the effective sampling points are intercepted by a speed-based method to obtain effective axle weight signals, so that the intercepted effective axle weight signals can be ensured to be accurate and reliable.
2. Aiming at the conditions that the nonlinear fitting algorithm, such as the Levenberg-Marquardt fitting algorithm, has high requirements on a hardware platform and cannot meet the practical application, the method for extracting and processing the effective axle weight signal is adopted, so that the calculation precision is ensured, the calculation complexity is simplified, and the method can be applied to the conventional dynamic weighing system.
3. Aiming at the condition that the frequency of the low-frequency interference signal is uncertain, the method adopts a frequency-variable nonlinear fitting algorithm to fit the interference signal, can obtain an approximate low-frequency interference signal, and improves the weighing precision of the dynamic axle load scale due to the improvement of the accuracy of removing the low-frequency interference signal.
Detailed Description
When a vehicle passes through a weighing platform 1 of the dynamic axle weight scale along the direction shown in figure 1, whether the vehicle passes through is identified by the connection and disconnection of light curtains 2 on the two sides of the weighing platform 1, then the wheel axle is identified through a tire identifier 3, a ground induction coil 4 is used as a standby of the light curtain 2, four weighing sensors or more weighing sensors arranged below the weighing platform 1 collect axle weight data, a signal processing platform based on an RAM is adopted, a 51 chip is used as a signal control platform, and the two are communicated through a parallel port.
The axle weight estimation method of the vehicle dynamic weighing system of the invention is shown in figure 2, a weighing sensor arranged below a weighing platform 1 collects an original axle weight signal of each axle, the original axle weight signal comprises a static axle weight signal, a low-frequency interference signal and a high-frequency interference signal, the number of sampling points is counted, and a discrete mathematical model of the original axle weight signal can be as follows:
in the above formula, y (n) is a discrete mathematical model of the original axle weight signal, w is the static axle weight of the vehicle, A is the amplitude of the low-frequency interference signal, f is the frequency of the low-frequency interference signal, and FS is set for the weighing systemThe sampling frequency of (a) is determined,
is the phase of the low-frequency interference signal, A
iIs the amplitude of the high-frequency interference signal, f
iIs the frequency of the high-frequency interference signal,
iis the phase of the high frequency interference signal.
The initial axle weight signal is preprocessed by an FIR (finite Impulse response) or IIR (infinite impulse response) type low-pass filter, high-frequency interference signals contained in the initial axle weight signal are filtered, the high-frequency interference signals with the frequency being more than or equal to 50Hz are filtered by the filter, and at the moment, the axle weight signal basically only contains static axle weight and low-frequency interference signals.
Calculating the speed of the axle according to the number of sampling points, the set sampling frequency and the stroke of the vehicle passing through the dynamic axle weight scale; the mathematical expression for the axle speed is: v is FS (L + delta)/length, wherein v is the speed of each axle when the vehicle passes through the dynamic axle weight scale, FS is the sampling frequency, L + delta is the stroke of the dynamic axle weight scale up and down, L is the platform width of the dynamic axle weight scale, delta is the compensation coefficient, delta is between 0.2 and 0.8, and length is the number of sampling points.
Intercepting effective sampling points when a vehicle passes through the dynamic axle weight scale according to the calculated axle speed to obtain an effective axle weight signal, as shown in figure 3, which is a waveform diagram of the effective axle weight signal, wherein the collected axle weight signal comprises three sections including an upper dynamic axle weight scale section of the vehicle, an effective weighing section of the vehicle on the dynamic axle weight scale and a lower dynamic axle weight scale section of the vehicle, so that the effective sampling points can be obtained by removing the sampling points of the upper dynamic axle weight scale section and the lower dynamic axle weight scale section of the vehicle, the effective sampling points can be obtained by calculation, the mathematical expression of the effective sampling points is L2 length-2L1, wherein L2 is the effective sampling points, L1 is the sampling points of the upper dynamic axle weight scale section and the lower dynamic axle weight scale section of the vehicle, and the sampling points L1 of the upper dynamic axle weight scale section and the lower dynamic axle weight scale section of the vehicle are inversely proportional to the axle speed v, so that the effective axle weight signal can be obtained by calculation, the mathematical expression is L1 ═ FS × S/v, wherein S is the passing displacement of the upper dynamic axle weight scale section and the lower dynamic axle weight scale section of the vehicle, and the fluctuation part of the upper dynamic axle weight scale and the lower dynamic axle weight scale of the vehicle is removed, so that the accuracy of vehicle weight estimation can be improved.
Extracting 10-150 data points on the basis of the effective axle weight signal, in the sampling process, adopting the modes of moving average, simple extraction and the like, wherein the sampled effective axle weight signal forms an axle weight signal for simplifying calculation, and the mathematical expression of the axle weight signal is as follows:
x(n)=w+A′*sin(2*pi*f*n*/FS′+
'), where x (n) is a discrete mathematical model, A ' is the amplitude of the decimated jammer signal, f is the frequency of the jammer signal, FS ' is the decimated sampling frequency,
the extracted phase can obtain effective axle weight signals with fewer points, so that the requirement that parameters required by nonlinear fitting are obtained by processing the effective axle weight signals by using the axle weight signals by using a nonlinear algorithm in embedded systems such as RAM (random access memory), DSP (digital signal processor) and the like can be met, the condition that four parameters, namely an axle weight initial value w, are required by fitting the model by using the nonlinear fitting algorithm can be known from a mathematical expression of the axle weight signals
0Amplitude initial value A of low-frequency interference signal
0Frequency initial value f of low-frequency interference signal
0And phase initialization value of low frequency interference signal
0. Wherein the initial value w of the axle weight
0Can be obtained by calculating the average value of the axle weight signal, and N is the number of data points after sampling, or the maximum value or the minimum value of the axle weight signal is adopted,
the mathematical expression is as follows: <math><mrow><msub><mi>w</mi><mn>0</mn></msub><mo>=</mo><mfenced open='{' close='' separators=' '><mtable><mtr><mtd><mfrac><mn>1</mn><mi>N</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>max</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr><mtr><mtd><mi>min</mi><mrow><mo>(</mo><mi>x</mi><mrow><mo>(</mo><mi>i</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd></mtr></mtable><mo>.</mo></mfenced></mrow></math>
and the amplitude initial value A of the low-frequency interference signal
0Obtaining the weight by using half of the maximum value Max and the minimum value Min in the axle weight signals, wherein the numerical expression is as follows:
and the frequency initial value f of the low-frequency interference signal
0The frequency of the sampling frequency FS 'can be calculated by selecting 1 Hz-5 Hz or by using the maximum value coordinate Maxp, the minimum value coordinate Minp and the extracted sampling frequency FS' of the corresponding axle weight signal, and the mathematical expression is f
0FS'/(Maxp-Minp |), the maximum value and the minimum value can be ensured within a certain range, and thus a range satisfying the frequency of the low-frequency interference signal can be ensured. Phase initialization value of low frequency interference signal
00.5-2 can be selected, or the first extreme point coordinate pp corresponding to the axle weight signal and the initial value f of the low-frequency interference signal frequency can be used
0Calculating the extracted sampling frequency FS ', and converting the sampling frequency FS' into a range of 0-2 pi, wherein the mathematical expression is as follows:
obtaining the initial phase value of the low-frequency interference signal
0Obtaining approximate low-frequency interference signal by using nonlinear fitting algorithmThe mathematical expression for approximating a low frequency interference signal is:
wherein,
to fit the resulting amplitude of the glitch,
for fitting the obtained frequency of the low-frequency interference
The phase obtained by fitting.
Finally, the approximate low-frequency interference signal is subtracted from the effective axle weight signal or the axle weight signal to obtain an estimated static axle weight signal, and as can be seen from the waveform diagram in fig. 4, the waveform of the estimated static axle weight signal is smoother than that of the axle weight signal.
The axle weight signal acquired by the axle weight estimation method of the vehicle dynamic weighing system is processed, the vehicle speed is less than or equal to 20km/h, the axle total weight estimation error is controlled within +/-2.5 percent, and the weighing precision of the dynamic axle weight scale is improved.