Disclosure of Invention
Technical problem to be solved
The invention discloses a signal identification method and a signal identification device for a vehicle-mounted range radar, which aim to at least partially solve the technical problems.
(II) technical scheme
In order to achieve the above object, an embodiment of the present invention provides a signal identification method for a vehicle-mounted range radar, including: acquiring a plurality of sampling signals, wherein each sampling signal is a target echo signal or an interference signal; processing each sampling signal to obtain a plurality of processed signals; performing Fourier transform on each processed signal to obtain amplitude-frequency response information corresponding to each sampling signal; determining a feature vector corresponding to each of the sampled signals according to the amplitude-frequency response information; training a binary classification vector machine by using the feature vector; and identifying the received signal by using the trained binary vector machine to determine that the received signal belongs to the target echo signal or the interference signal.
According to an embodiment of the invention, the interfering signal comprises at least one of a noise amplitude modulated interfering signal and a sinusoidal amplitude modulated interfering signal.
According to an embodiment of the present invention, processing each of the sampled signals to obtain a plurality of processed signals includes: calculating the mean value of the time domain amplitude of each sampling signal; and subtracting the mean value of the time domain amplitude of each sampling signal from the time domain amplitude of each sampling signal to obtain the processed signal.
According to an embodiment of the present invention, the feature vector of each of the sampled signals includes at least one of a mean, a standard deviation, a kurtosis, a difference between a maximum and a minimum, and an average power spectral density of amplitude-frequency response information of each of the sampled signals.
According to an embodiment of the present invention, the kurtosis is calculated as follows:
wherein K is kurtosis, s (f)
iInformation representing the amplitude-frequency response of each of said sampled signals, σ represents the standard deviation of the amplitude-frequency response of each of said sampled signals,
the average value of the amplitude-frequency response information of each sampling signal is obtained, n is the number of sampling points of each sampling signal, and i is any one positive integer from 1 to n.
According to an embodiment of the present invention, the difference between the maximum value and the minimum value is calculated as follows: and pk is the difference between the maximum value and the minimum value, max is the maximum value of the amplitude-frequency response information of each sampling signal, and min is the minimum value of the amplitude-frequency response information of each sampling signal.
According to an embodiment of the invention, the average power spectral density is calculated by:
wherein E is an average power spectral density, f is a frequency of each of the sampling signals, s (f) is amplitude-frequency response information of each of the sampling signals, and N is the number of sampling points of each of the sampling signals.
According to the embodiment of the invention, a Gaussian radial basis function is preset in the two-classification vector machine, and the training of the two-classification vector machine by using the feature vector comprises the following steps: obtaining a feature vector sample; normalizing the feature vector sample to obtain a target feature vector sample, wherein the maximum value and the minimum value of each type of feature vector in the target feature vector sample are respectively kept consistent; and inputting the target feature vector sample into a two-classification vector machine, and outputting a classification decision function containing target parameters, wherein the target parameters comprise a penalty parameter C and a kernel function parameter gamma.
According to an embodiment of the present invention, identifying the received signal by using the trained binary vector machine to determine that the received signal belongs to the target echo signal or the interference signal includes: identifying the received signals by using the classification decision function to obtain interface information for distinguishing the target echo signals or the interference signals; and determining the target echo signal or the interference signal according to the interface information.
The invention further provides a signal identification device for the vehicle-mounted range radar, which is used for realizing the method.
(III) advantageous effects
According to the method and the device, the signal identification can be carried out according to different amplitude-frequency response characteristics of the target signal and the interference signal, so that the ranging capability of the vehicle-mounted ranging radar in a complex road environment is improved.
Detailed Description
In order that the objects, technical solutions and advantages of the present invention will become more apparent, the present invention will be further described in detail with reference to the accompanying drawings in conjunction with the following specific embodiments.
An embodiment of the present invention provides a signal identification method for a vehicle-mounted range radar, including: acquiring a plurality of sampling signals, wherein each sampling signal is a target echo signal or an interference signal; processing each sampling signal to obtain a plurality of processed signals; fourier transform is carried out on each processed signal to obtain amplitude-frequency response information corresponding to each sampling signal; determining a feature vector corresponding to each sampling signal according to the amplitude-frequency response information; training a two-classification vector machine by using the feature vectors; and identifying the received signal by using a trained two-classification vector machine to determine that the received signal belongs to a target echo signal or an interference signal.
Fig. 1 schematically shows a flowchart of a signal identification method for a vehicle-mounted ranging radar according to an embodiment of the present invention.
As shown in fig. 1, the flow includes operations S101 to S105.
In operation S101, a target echo signal and an interference signal are acquired.
According to an embodiment of the present invention, the target echo signal and the interference signal constitute, for example, the plurality of sampling signals, for example, sampling signals in the form of doppler frequencies. Taking a specific embodiment as an example, the sampling process of each target echo signal or interference signal may include: the method includes the steps of transmitting a sinusoidal signal by a modulation signal generator, forming a signal with a center frequency of 1GHz (for example only, which can be randomly adjusted according to actual application scenarios, the same applies hereinafter) after passing through an oscillator, wherein frequency changes are triangular waves, namely local oscillators, the offset of the local oscillators is 25MHz/V (for example only), generating the local oscillators by the modulation signal generator and the oscillator, generating a return signal by a circulator, mixing the local oscillators and the return signal to generate a difference frequency signal and a sum frequency signal, passing the difference frequency signal through a band-pass filter to obtain a required harmonic frequency, performing distance separation, performing secondary mixing on the harmonic signal, and filtering to obtain a Doppler frequency, namely the target echo signal.
According to an embodiment of the present invention, the interference signal may comprise at least one of a noise amplitude modulated interference signal and a sinusoidal amplitude modulated interference signal, for example. After the target echo signal is acquired in the foregoing manner, for example, the original circulator in the process for acquiring the target echo signal may be replaced by a sinusoidal signal transmitter and a noise transmitter, respectively, where the center frequency is aligned with the center frequency of the target echo signal, the frequency of the interference signal is, for example, 50Hz (by way of example only), and the step length (the time taken for each sweep of the sweep frequency, and after several sweeps, the coverage on all the frequencies to be swept is completed) is 0.001s, for 500 steps. By changing the interference frequency of the interference signal and increasing the distance, the two-way distance of the echo signal is changed into the one-way distance, and the acquisition of the noise amplitude modulation interference signal and the sine amplitude modulation interference signal can be realized.
According to the embodiment of the present invention, for example, 100 target echo signals under the undisturbed action, 100 noise amplitude modulation interference signals and 100 sinusoidal amplitude modulation interference signals may be collected as sampling signals, the sampling frequency is, for example, 100kHz, and the sampling time is, for example, 0.019 s.
In operation S102, a fast fourier transform is performed.
According to an embodiment of the present invention, the operation may further include, for example, processing each of the sampled signals to obtain a plurality of processed signals, and performing fourier transform on each of the processed signals to remove zero-frequency components, i.e., direct-current components, and obtain amplitude-frequency response information corresponding to each of the sampled signals.
Fig. 2 schematically shows a resultant graph of amplitude-frequency response information of a target echo signal according to an embodiment of the present invention.
Fig. 3 schematically shows a resulting graph of amplitude-frequency response information of a noise amplitude modulated jamming signal according to an embodiment of the invention.
Fig. 4 schematically shows a resulting graph of amplitude-frequency response information of a sinusoidal amplitude modulated jamming signal according to an embodiment of the invention.
As shown in fig. 2 to 4, it can be seen that the amplitude-frequency response peak of the target echo signal is relatively prominent, and the amplitude-frequency response peak of the noise signal is relatively small.
According to an embodiment of the present invention, the above-mentioned manner of processing each sampling signal may include, for example: calculating the mean value of the time domain amplitude of each sampling signal; and subtracting the mean value of the time domain amplitude of each sampling signal from the time domain amplitude of each sampling signal to obtain the processed signal.
In operation S103, a feature vector is extracted.
According to an embodiment of the present invention, the operation may include, for example, determining a feature vector corresponding to each of the sampled signals according to the amplitude-frequency response information, where the feature vector of each of the sampled signals may include, for example, at least one of a mean, a standard deviation, a kurtosis, a difference between maximum and minimum values, and an average power spectral density of the amplitude-frequency response information of each of the sampled signals.
Fig. 5 schematically shows a flow chart for obtaining feature vectors according to an embodiment of the invention.
As shown in fig. 5, the process includes calculating a mean value, a standard deviation, a kurtosis, a difference between a maximum value and a minimum value, an average power spectral density of a frequency spectrum of the magnitude-frequency response information obtained after fourier transform, and determining a feature vector.
According to the embodiment of the invention, the mean value, the standard deviation, the kurtosis, the difference between the maximum and the minimum values and the average power spectral density of the amplitude-frequency response of the processed signal can be obtained through a known function.
According to an embodiment of the present disclosure, the kurtosis may also be calculated, for example, by the following formula:
wherein s (f)
iInformation representing the amplitude-frequency response of each of said sampled signals, σ represents the standard deviation of the amplitude-frequency response of each of said sampled signals,
and taking the average value of the amplitude-frequency response information of each sampling signal, wherein n is the number of sampling points of each sampling signal, and i is any one positive integer from 1 to n.
According to an embodiment of the present disclosure, the difference between the maximum value and the minimum value may also be calculated, for example, as follows:
pk=max-min
and the method comprises the following steps of sampling signals, wherein pk is the difference between the maximum value and the minimum value, max is the maximum value of amplitude-frequency response information of each sampling signal, and min is the minimum value of amplitude-frequency response information of each sampling signal.
According to an embodiment of the present disclosure, the average power spectral density may also be calculated, for example, by:
wherein E is an average power spectral density, f is a frequency of each of the sampling signals, s (f) is amplitude-frequency response information of each of the sampling signals, and N is the number of sampling points of each of the sampling signals.
According to an embodiment of the present disclosure, referring to fig. 5, in the process of determining the feature vector, for example, the kurtosis, the difference between the maximum and minimum values, and the average power spectral density obtained through the above calculation may be used as a feature vector, and the feature vector may be expressed in a form of d ═ K, pk, E, for example, and used as an input for supporting a binary classification vector machine.
According to the embodiment of the invention, in order to better know whether the feature vector is effective, Kruskal-Wallis is used to test whether the distribution of a plurality of populations has significant difference analysis, and the original hypothesis H is0The p-value of the returned test result after the non-parametric analysis of the data is carried out on a plurality of independent samples from the same population or a plurality of populations generating the independent samples according to the same distribution Kruskal-Wallis, the p-value is an important parameter of the hypothesis test, the smaller the p-value, the more remarkable the result is, namely H is rejected0The more sufficient the reason for or the more confident we are to refuse to accept H0And accept alternative hypothesis H1The distribution of the various populations varies significantly and the results of the test are considered reliable for p-values less than 0.01.
In operation S104, a binary class vector machine is trained.
According to the embodiment of the invention, since the weight vector of the support vector machine must be mapped from a low dimension to a high dimension, the adoption of the kernel function becomes an indispensable condition, which makes the sample linearly separable, and based on this, the two-classification vector machine is preset with a gaussian radial basis kernel function, for example. Because some sample points may fall between the hyperplane and the boundary, a slack variable needs to be introduced, for each slack variable, a cost needs to be paid, a penalty parameter C is a penalty for the slack variable, the larger the C value is, the smaller the C value is, the misclassification point is as small as possible, so as to make the interval as large as possible, for example, C may be introduced into the two-classification vector machine to make a harmonic coefficient.
According to an embodiment of the present invention, the operation S401 may be represented by, for example, training a binary vector machine using the feature vector, and specifically, in the case of training samples with an interference signal, the operation may include: obtaining a feature vector sample; normalizing the feature vector sample to obtain a target feature vector sample, wherein the maximum value and the minimum value of each type of feature vector in the target feature vector sample are respectively kept consistent; and inputting the target feature vector sample into a two-classification vector machine, and outputting a classification decision function containing target parameters, wherein the target parameters comprise a penalty parameter C and a kernel function parameter gamma.
According to an embodiment of the present invention, before the obtaining the feature vector samples, for example, the method may further include: and establishing an interface between the interference signal and the target echo signal characteristic parameter value by adopting a two-classification SVM (vector machine), wherein the interference signal and the target echo signal are respectively positioned at two sides of the interface. By using a standard C-SVM and taking three parameters (such as the above-mentioned feature vector d) as input of the SVM, the kernel function selects a gaussian radial basis kernel function, and the planning problem can be solved under the condition that the parameters C i0 and γ 0.01, so that four support vectors and corresponding classification decision functions can be obtained, and the interface between the interference signal and the target signal can be obtained.
According to the embodiment of the invention, the target parameters can be determined by using a 10-fold cross-validation method, for example, by randomly dividing the original training set into 10 parts without repetition, selecting one part as a validation set, using the remaining 9 parts as the training set for model training, obtaining a model after training on the training set, testing on the validation set by using the model, and storing the evaluation indexes of the model for 10 times (ensuring that each subset has one opportunity to serve as the validation set).
According to an embodiment of the present invention, in order to further optimize the target parameter, for example, the feature vector may be divided into a training set and a test set by using a 10-fold cross-checking method, specifically, for example, the feature vector sample data after the normalization processing may be divided into 10 samples, where one sample is used as the test set, and the rest samples are used as the training set, and in order to further determine the parameter, for example, the obtained parameter value may be divided into a grid pattern, and each grid uses the 10-fold cross-checking method to calculate the verification classification accuracy so as to obtain the optimal target parameter value. After obtaining the optimal target parameter value, 3/4 samples (i.e., 75 sets of target echo signals and 150 sets of interference signals) are randomly taken out from the target signal samples and the interference signal samples respectively as training samples, and a learning training of the SVM classifier is performed, for example, the above-mentioned classification decision function containing the target parameter is obtained.
According to the embodiment of the present invention, after obtaining the classification decision function, for example, the remaining 1/4 sampling samples (i.e., 25 sets of target echo signals and 50 sets of interference signals) may be used as test samples, and input to train to obtain the SVM classifier, so as to obtain the target detection rate and the interference detection rate of the current experiment. By performing the random grouping test for 200 times, and finally averaging the target detection rate and the interference detection rate obtained each time, the final target detection rate and the final interference detection rate can be obtained, and the average value of the classification accuracy can be used as the performance index of the machine algorithm.
According to the embodiment of the invention, by continuously training the SVM model and performing simulation test, for example, a visual implementation manner of a classification decision function containing target parameters can be finally obtained.
In operation S105, a classification recognition capability is obtained.
According to an embodiment of the present invention, the operation may be represented as: and identifying the received signal by using a trained two-classification vector machine to determine that the received signal belongs to a target echo signal or an interference signal. Specifically, the implementation process may include, for example: identifying the received signals by using a classification decision function to obtain interface information for distinguishing target echo signals or the interference signals; and determining a target echo signal or an interference signal according to the interface information.
Another embodiment of the present invention provides a signal recognition apparatus for a vehicle-mounted ranging radar, which may implement the method as described above.
By the method and the device, signal identification can be performed according to different amplitude-frequency response characteristics of the target signal and the interference signal, so that the ranging capability of the vehicle-mounted ranging radar in a complex road environment is improved.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention and are not intended to limit the present invention, and any modifications, equivalents, improvements and the like made within the spirit and principle of the present invention should be included in the protection scope of the present invention.