CN111245490A - Broadband signal extraction method and device and electronic equipment - Google Patents

Broadband signal extraction method and device and electronic equipment Download PDF

Info

Publication number
CN111245490A
CN111245490A CN201911090146.9A CN201911090146A CN111245490A CN 111245490 A CN111245490 A CN 111245490A CN 201911090146 A CN201911090146 A CN 201911090146A CN 111245490 A CN111245490 A CN 111245490A
Authority
CN
China
Prior art keywords
signal
initial
time
interference
broadband
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
CN201911090146.9A
Other languages
Chinese (zh)
Other versions
CN111245490B (en
Inventor
潘峰
刘亚奇
师一帅
范鹏程
吕万里
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.)
Individual
Original Assignee
Individual
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 Individual filed Critical Individual
Priority to CN201911090146.9A priority Critical patent/CN111245490B/en
Publication of CN111245490A publication Critical patent/CN111245490A/en
Application granted granted Critical
Publication of CN111245490B publication Critical patent/CN111245490B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/0848Joint weighting
    • H04B7/0854Joint weighting using error minimizing algorithms, e.g. minimum mean squared error [MMSE], "cross-correlation" or matrix inversion
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/08Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
    • H04B7/0837Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station using pre-detection combining
    • H04B7/0842Weighted combining
    • H04B7/086Weighted combining using weights depending on external parameters, e.g. direction of arrival [DOA], predetermined weights or beamforming
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02DCLIMATE CHANGE MITIGATION TECHNOLOGIES IN INFORMATION AND COMMUNICATION TECHNOLOGIES [ICT], I.E. INFORMATION AND COMMUNICATION TECHNOLOGIES AIMING AT THE REDUCTION OF THEIR OWN ENERGY USE
    • Y02D30/00Reducing energy consumption in communication networks
    • Y02D30/70Reducing energy consumption in communication networks in wireless communication networks

Abstract

The embodiment of the disclosure provides a broadband signal extraction method, a broadband signal extraction device and electronic equipment, belonging to the technical field of signal processing, wherein the method comprises the following steps: receiving an initial signal to be extracted, wherein the initial signal comprises an initial expected signal and an interference signal; transforming the initial signal to a time-frequency domain by utilizing bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix on each time-frequency point; calculating to obtain a single-source self-item set of the expected signal and a single-source self-item set of the interference signal according to the space-time frequency distribution matrix; estimating a steering vector of the expected signal by utilizing two types of single-source self-entry sets, and reconstructing a broadband interference and noise covariance matrix; and inputting the expected signal steering vector and the broadband interference and noise covariance matrix into a broadband beam former to obtain a broadband expected signal corresponding to the initial signal. By the broadband signal extraction scheme, the stability and the output performance of extracting broadband signals in a complex environment are improved.

Description

Broadband signal extraction method and device and electronic equipment
Technical Field
The present disclosure relates to the field of signal processing technologies, and in particular, to a method and an apparatus for extracting a broadband signal, and an electronic device.
Background
At present, the method has entered into the comprehensive information era, and with the development of technologies such as wireless communication and radar detection, the space electromagnetic environment is becoming more and more complex, and the interference signal is usually overlapped with the desired signal in the time domain and is located in the same frequency band in the frequency domain, so that the conventional time domain and frequency domain processing method is difficult to realize the extraction of the desired signal, and therefore, the spatial filtering technology based on the array antenna is rapidly growing up and is widely applied. At the present stage, in practical application, a fixed beam forming mode is mostly adopted in the medium-space filtering technology, and the mode has strong stability and can suppress interference and noise to a certain extent. However, the performance advantages of the array antenna cannot be fully exerted, the spatial sensing capability is not provided, and the capability of inhibiting accidental strong interference is poor. Compared with the prior art, the self-adaptive beam forming technology can automatically adjust the weighting coefficient according to the electromagnetic environment, thereby effectively inhibiting various interferences and noises and realizing the optimal output of the expected signal. However, the adaptive beamforming technology has strict requirements on information such as array element position and target arrival angle, and when there is an error in the a priori information, the adaptive beamformer suppresses the desired signal as interference, which causes a sharp drop in output performance, which is also a main reason for the small practical application. Therefore, how to implement robustness of adaptive beamforming becomes a hot spot of current research.
Therefore, the existing broadband signal extraction method has the technical problems of poor stability and poor output performance of extracting broadband signals in a complex environment.
Disclosure of Invention
In view of the above, the disclosed embodiments provide a wideband signal extraction method, which at least partially solves the problems in the prior art.
In a first aspect, an embodiment of the present disclosure provides a wideband signal extraction method, where the method includes:
receiving an initial signal to be extracted, wherein the initial signal comprises an initial expected signal and an interference signal;
transforming the initial signal to a time-frequency domain by utilizing bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix on each time-frequency point;
calculating according to the space-time frequency distribution matrix to obtain a single-source self-item set of an initial expected signal and a single-source self-item set of an interference signal;
estimating a steering vector of the initial expected signal by using the set of single-source self-terms of the initial expected signal and the set of single-source self-terms of the interference signal, and reconstructing a broadband interference and noise covariance matrix;
and inputting the steering vector of the initial expected signal and the broadband interference and noise covariance matrix into a broadband beam former to obtain a broadband expected signal corresponding to the initial expected signal of the initial signal.
According to a specific implementation manner of the embodiment of the present disclosure, the space-time frequency distribution matrix includes a delay-free array received signal vector, a conjugate transpose and a kernel function.
According to a specific implementation manner of the embodiment of the present disclosure, the step of obtaining a single source self-entry set of an initial desired signal and a single source self-entry set of an interference signal by calculating according to the space-time frequency distribution matrix includes:
calculating a time delay vector according to the space-time-frequency distribution matrix;
and respectively calculating to obtain a single source origin item set of the initial expected signal and a single source origin item set of the interference signal by using the time delay vector and the initial signal.
According to a specific implementation manner of the embodiment of the present disclosure, the step of respectively calculating a single source self-entry set of the initial desired signal and a single source self-entry set of the interference signal by using the delay vector and the initial signal includes:
taking the time delay vector as a feature, extracting a time frequency point corresponding to the initial expected signal from the initial signal to form a time frequency point set corresponding to the initial expected signal, and taking the time frequency points left after the time frequency point set of the initial expected signal is extracted from the initial signal as the time frequency point set of an interference signal and a noise signal;
and removing the mutual item point, the mutual item point and the noise point in the time-frequency point set of the initial expected signal and the time-frequency point set of the interference signal and the noise signal by utilizing the single-source item time-frequency point to obtain the initial expected signal single-source item point set and the interference signal single-source item point set.
According to a specific implementation manner of the embodiment of the present disclosure, the step of extracting, with the delay vector as a feature, a time-frequency point corresponding to the initial desired signal from the initial signal, as a time-frequency point set of the initial desired signal, and taking remaining time-frequency points as time-frequency point sets of an interference signal and a noise signal includes:
obtaining a calculation result by using the time delay vector and a preset time delay vector of the initial expected signal according to a preset algorithm;
and comparing the calculation result with a first threshold, taking the time frequency points less than or equal to the first threshold as a time frequency point set of the initial expected signal, and taking the time frequency points more than the first threshold as a time frequency point set of the interference and noise signal.
According to a specific implementation manner of the embodiment of the present disclosure, the step of obtaining the single source self-item set of the initial desired signal and the single source self-item set of the interference signal by removing the mutual item, the mutual self-item, and the noise point in the time-frequency point set of the initial desired signal and the time-frequency point set of the interference signal and the noise signal by using the single source self-item time-frequency point includes:
and according to a second threshold value and the maximum characteristic value, eliminating mutual item points, interactive self item points and noise points from the time frequency point set of the initial expected signal to obtain a single-source self item point set of the expected signal, and eliminating mutual item points, interactive self item points and noise points from the time frequency point set of the interference signal and the noise signal to obtain a single-source self item point set of the interference signal.
According to a specific implementation of the embodiment of the present disclosure, the estimating a steering vector of the initial desired signal by using the set of mono-origin terms of the initial desired signal and the set of mono-origin terms of the interference signal, and reconstructing a wideband interference-plus-noise covariance matrix includes:
taking the time delay vector as a characteristic, and dividing the single-source item point set of the interference signal into a single-source item point set of each interference signal by using a mean value clustering method;
calculating a steering vector of the initial signal on different frequencies by using the single-origin-term set of the initial expected signal and the single-origin-term set of each interference signal, wherein the steering vector comprises the steering vector of the initial expected signal and the steering vector of the interference signal;
and reconstructing the broadband interference and noise covariance matrix according to the guide vector of the interference signal.
According to a specific implementation manner of the embodiment of the present disclosure, the step of inputting the steering vector of the initial desired signal and the wideband interference-plus-noise covariance matrix into a wideband beamformer to obtain a wideband desired signal corresponding to the initial desired signal of the initial signal includes:
calculating a constraint matrix by using the guide vector of the initial expected signal;
calculating a weighting coefficient according to the constraint matrix and the broadband interference and noise covariance matrix;
and performing spatial filtering on the initial signal by using the weighting coefficient to extract a broadband expected signal.
According to a specific implementation of the embodiments of the present disclosure, the constraint matrix includes steering vectors of the initial desired signal at different frequencies.
In a second aspect, an embodiment of the present disclosure provides a wideband signal extracting apparatus, including:
the device comprises a receiving module, a processing module and a processing module, wherein the receiving module is used for receiving an initial signal to be extracted, and the initial signal comprises an initial desired signal and an interference signal;
the transformation module is used for transforming the initial signal to a time-frequency domain by utilizing bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix on each time-frequency point;
the calculation module is used for calculating to obtain a single-source self-item set of an initial expected signal and a single-source self-item set of an interference signal according to the space-time frequency distribution matrix;
an estimation module, configured to estimate a steering vector of the initial desired signal using the set of single-source self-terms of the initial desired signal and the set of single-source self-terms of the interference signal, and reconstruct a wideband interference-plus-noise covariance matrix;
and the acquisition module is used for inputting the steering vector of the initial expected signal and the broadband interference and noise covariance matrix into a broadband beam former to obtain a broadband expected signal corresponding to the initial expected signal of the initial signal.
In a third aspect, an embodiment of the present disclosure further provides an electronic device, where the electronic device includes:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of extracting a broadband signal of any implementation of the first aspect or the first aspect.
The broadband signal extraction scheme in the embodiment of the disclosure comprises receiving an initial signal to be extracted, wherein the initial signal comprises an initial expected signal and an interference signal; transforming the initial signal to a time-frequency domain by utilizing bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix on each time-frequency point; calculating to obtain a single-source self-item set of the expected signal and a single-source self-item set of the interference signal according to the space-time frequency distribution matrix; estimating a steering vector of the desired signal by using the set of single-source self-terms of the desired signal and the set of single-source self-terms of the interference signal, and reconstructing a broadband interference-plus-noise covariance matrix; and inputting the steering vector of the expected signal and the broadband interference and noise covariance matrix into a broadband beam former to obtain a broadband expected signal corresponding to the initial expected signal of the initial signals. Through the scheme disclosed by the invention, the stability and the output performance of extracting the broadband expected signal in a complex environment are improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present disclosure, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present disclosure, and it is obvious for those skilled in the art to obtain other drawings without creative efforts.
Fig. 1 is a schematic flowchart of a wideband signal extraction method provided in an embodiment of the present disclosure;
fig. 2 is a block diagram of a wideband beamformer provided by an embodiment of the present disclosure;
fig. 3 is a diagram of a delay vector estimation result provided by the embodiment of the present disclosure;
fig. 4 is a diagram of a spatial filtering result provided by an embodiment of the present disclosure;
FIG. 5 is a graph of output performance versus number of sample points for embodiments of the present disclosure and other methods;
FIG. 6 is a graph of output performance versus input signal-to-noise ratio for embodiments of the present disclosure and other algorithms;
fig. 7 is a schematic structural diagram of a wideband signal extraction apparatus provided in an embodiment of the present disclosure;
fig. 8 is a schematic view of an electronic device provided in an embodiment of the present disclosure.
Detailed Description
The embodiments of the present disclosure are described in detail below with reference to the accompanying drawings.
The embodiments of the present disclosure are described below with specific examples, and other advantages and effects of the present disclosure will be readily apparent to those skilled in the art from the disclosure in the specification. It is to be understood that the described embodiments are merely illustrative of some, and not restrictive, of the embodiments of the disclosure. The disclosure may be embodied or carried out in various other specific embodiments, and various modifications and changes may be made in the details within the description without departing from the spirit of the disclosure. It is to be noted that the features in the following embodiments and examples may be combined with each other without conflict. All other embodiments, which can be derived by a person skilled in the art from the embodiments disclosed herein without making any creative effort, shall fall within the protection scope of the present disclosure.
It is noted that various aspects of the embodiments are described below within the scope of the appended claims. It should be apparent that the aspects described herein may be embodied in a wide variety of forms and that any specific structure and/or function described herein is merely illustrative. Based on the disclosure, one skilled in the art should appreciate that one aspect described herein may be implemented independently of any other aspects and that two or more of these aspects may be combined in various ways. For example, an apparatus may be implemented and/or a method practiced using any number of the aspects set forth herein. Additionally, such an apparatus may be implemented and/or such a method may be practiced using other structure and/or functionality in addition to one or more of the aspects set forth herein.
It should be noted that the drawings provided in the following embodiments are only for illustrating the basic idea of the present disclosure, and the drawings only show the components related to the present disclosure rather than the number, shape and size of the components in actual implementation, and the type, number and ratio of the components in actual implementation may be changed arbitrarily, and the layout of the components may be more complicated.
In addition, in the following description, specific details are provided to facilitate a thorough understanding of the examples. However, it will be understood by those skilled in the art that the aspects may be practiced without these specific details.
The embodiment of the disclosure provides a broadband signal extraction method. The broadband signal extraction method provided by the present embodiment may be executed by a computing device, which may be implemented as software, or implemented as a combination of software and hardware, and may be integrally provided in a server, a terminal device, or the like.
Referring to fig. 1, a wideband signal extraction method provided by the embodiment of the present disclosure includes:
s101, receiving an initial signal to be extracted, wherein the initial signal comprises an initial expected signal and an interference signal.
The broadband signal extraction method provided by the embodiment of the disclosure can be applied to the extraction process of broadband signals in the environments such as radio stations, houses and the like.
The electronic device may be internally provided with a signal receiving module or externally connected with a signal receiving device for receiving an initial signal to be extracted. Where the received initial signal typically comprises an initial desired signal corresponding to the final output requirement and other interfering and noise signals than the initial desired signal, the purpose of the wideband signal extraction is to extract a final optimized signal corresponding to the initial desired signal from the mixed initial signal. After receiving the initial signal, the signal receiving module or the signal receiving device may directly send the initial signal to the processor for subsequent analysis processing operation, or may store the received initial signal in a predetermined storage space, and call the initial signal from the predetermined storage space when the analysis processing is required.
And S102, transforming the initial signal to a time-frequency domain by utilizing bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix on each time-frequency point.
The electronic device may transform the received initial signal to a time-frequency domain through bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix at each time-frequency point through calculation. For example, the following calculation formula can be adopted for the space-time-frequency domain conversion:
Figure RE-GDA0002455628270000071
in the above-mentioned formula,
Figure RE-GDA0002455628270000072
representing a non-delayed array received signal vector, (-)HThe expression is to take conjugate transpose, phi (i, l) is a kernel function, and different kernel functions can be selected according to different requirements.
And S103, calculating to obtain a single source self-item set of the initial expected signal and a single source self-item set of the interference signal according to the space-time frequency distribution matrix.
After the space-time frequency distribution matrix is obtained through calculation, the time-frequency points corresponding to the initial signals can be classified by using the space-time frequency distribution matrix, and a single-source self-item set of the initial expected signals and a single-source self-item set of the interference signals are obtained.
And S104, estimating a guide vector of the initial expected signal by using the single-source term set of the initial expected signal and the single-source term set of the interference signal, and reconstructing a broadband interference and noise covariance matrix.
And calculating a guide vector of the expected signal and a guide vector of the interference signal according to a single-source self-item set of the initial expected signal and a single-source self-item set of the interference signal obtained by classifying the time-frequency points corresponding to the initial signal, and reconstructing a broadband interference and noise covariance matrix according to the obtained guide vector of the interference signal.
And S105, inputting the steering vector of the initial expected signal and the broadband interference and noise covariance matrix into a broadband beam former to obtain a broadband expected signal corresponding to the initial expected signal of the initial signal.
The wideband beamformer is a device that spatially filters the received original signal, as shown in fig. 2.
And utilizing a broadband beam former, and calculating the steering vector of the initial expected signal obtained by the steps and the broadband interference and noise covariance matrix to obtain a weighting coefficient, and being characterized by being capable of carrying out spatial filtering without delay compensation. And performing space domain filtering on the initial signal by using the broadband beam former to obtain a broadband expected signal.
The wideband signal extraction method provided by the embodiment of the disclosure is directed at classifying received signals in a complex electromagnetic signal environment, respectively obtaining a single-source self-point set of an initial desired signal and a single-source self-point set of an interference signal, then accurately estimating an initial desired signal steering vector and reconstructing an interference-plus-noise covariance matrix by using the two sets, and stably and efficiently extracting a wideband desired signal through a wideband beam former without delay compensation. The stability and the output performance of extracting the broadband signal are improved.
On the basis of the foregoing embodiment of the present disclosure, the step of calculating, according to the space-time frequency distribution matrix, to obtain a single source self-entry set of the initial desired signal and a single source self-entry set of the interference signal in step S103 may further include:
calculating a time delay vector according to the space-time-frequency distribution matrix;
and respectively calculating to obtain a single source origin item set of the initial expected signal and a single source origin item set of the interference signal by using the time delay vector and the initial signal.
Further, the step of respectively calculating a single origin and self-item set of the initial desired signal and a single origin and self-item set of the interference signal by using the delay vector and the initial signal may further include:
taking the time delay vector as a feature, extracting a time frequency point corresponding to an initial expected signal from the initial signal to form a time frequency point set corresponding to the initial expected signal, and taking the time frequency points left after the time frequency point set corresponding to the initial expected signal is extracted from the initial signal as a time frequency point set of an interference signal and a noise signal;
and removing the mutual item point, the mutual item point and the noise point in the time-frequency point set of the initial expected signal and the time-frequency point set of the interference signal and the noise signal by utilizing the single-source item time-frequency point to obtain the initial expected signal single-source item point set and the interference signal single-source item point set.
Further, the step of extracting, with the delay vector as a feature, a time-frequency point corresponding to the initial desired signal from the initial signal, as a time-frequency point set of the initial desired signal, and taking the remaining time-frequency points as a time-frequency point set of an interference signal and a noise signal may further include:
obtaining a calculation result by utilizing the time delay vector and a preset time delay vector of the initial expected signal according to a preset algorithm;
and comparing the calculation result with a first threshold, taking the time frequency points less than or equal to the first threshold as a time frequency point set of the initial expected signal, and taking the time frequency points more than the first threshold as a time frequency point set of the interference and noise signal.
In a specific calculation process, the time-frequency points related to the initial expectation satisfy:
Figure BDA0002266606050000091
wherein
Figure RE-GDA0002455628270000092
A predetermined delay vector for the initial desired signal (determined using the known DOA with error and the position of the elements). Epsilon1Is a constant threshold, usually ε1Is less than 0.5. By using the formula, a time frequency point set corresponding to the initial expected signal can be obtained, and the rest time frequency points correspond to an interference signal and noise time frequency point set.
In another specific implementation manner of the embodiment of the present disclosure, the step of obtaining the single source self-item set of the initial desired signal and the single source self-item set of the interference signal by removing the mutual item, the mutual self-item, and the noise point in the time-frequency point set of the initial desired signal and the time-frequency point set of the interference signal and the noise signal by using the single source self-item time-frequency point includes:
and according to a second threshold value and the maximum characteristic value, eliminating mutual item points, interactive self item points and noise points from the time frequency point set of the initial expected signal to obtain a single-source self item point set of the expected signal, and eliminating mutual item points, interactive self item points and noise points from the time frequency point set of the interference signal and the noise signal to obtain a single-source self item point set of the interference signal.
For example, the delay vector on each time frequency point is calculated according to the space-time frequency distribution matrix, then the delay vector on each time frequency point is used as a feature to extract the time frequency point corresponding to the initial expected signal, the time frequency point is defined as the time frequency point set of the initial expected signal, and the time frequency points left after the time frequency point corresponding to the initial expected signal is extracted are defined as the time frequency point set of the interference signal and the noise signal. And then according to a preset algorithm, obtaining a calculation result by using the delay vector and a preset delay vector of the expected signal, wherein the algorithm formula can be expressed as:
Figure BDA0002266606050000101
wherein
Figure BDA0002266606050000102
Is the predetermined delay vector of the desired signal (found using the known DOA with error and the array element position). Epsilon1Is a constant threshold, usually ε1<0.5. Comparing a calculation result obtained by using the delay vector and a preset delay vector of the expected signal with the first threshold, taking the time frequency points smaller than or equal to the first threshold as a time frequency point set of the initial expected signal, and taking the time frequency points larger than the first threshold as a time frequency point set of the interference and noise signal. Then extracting formula according to time-frequency point of single derived term
Figure BDA0002266606050000103
Wherein λ is0(t,f),λ1(t,f),L,λM-1(t, f) is Dx%x%M characteristic values of (t, f), λmax(t, f) is the maximum eigenvalue. Epsilon2Is a constant threshold, usually 0.8 ≦ ε2<1. Initializing the initial value according to a second threshold value and a maximum eigenvalueAnd rejecting mutual item points, interactive self item points and noise points in the time frequency point set of the expected signal and the time frequency point set of the interference signal and the noise signal to obtain a time frequency point set of the initial expected signal to serve as the single-source self item set of the initial expected signal, and obtain a time frequency point set of the interference signal and the noise signal to serve as the single-source self item set of the interference signal.
On the basis of the foregoing embodiment of the present disclosure, the step of estimating a steering vector of the initial desired signal by using the set of single-source terms of the initial desired signal and the set of single-source terms of the interference signal, and reconstructing a wideband interference-plus-noise covariance matrix in S104 may further include:
taking the time delay vector as a characteristic, and dividing the single-source item point set of the interference signal into a single-source item point set of each interference signal by using a mean value clustering method;
calculating a steering vector of the initial signal on different frequencies by using the single-origin-term set of the initial expected signal and the single-origin-term set of each interference signal, wherein the steering vector comprises the steering vector of the initial expected signal and the steering vector of the interference signal;
and reconstructing the broadband interference and noise covariance matrix according to the guide vector of the interference signal.
For example, the time-frequency points in the interference signal single-source item point set are divided into K-1 omega types by using a K-means clustering method12,L,ΩK-1Defined as a single source per interfering signal from a set of terms. Calculating a steering vector of the initial signal on different frequencies according to the single origin self-term set of each interference signal, wherein the steering vector comprises a steering vector of the initial expected signal and a steering vector of the interference signal, and the formula can be expressed as:
Figure RE-GDA0002455628270000111
wherein
Figure BDA0002266606050000112
Which represents the product of the Kronecker reaction,
Figure BDA0002266606050000113
delay T only with tapssThis can be directly obtained.
And obtaining an expression of interference plus noise covariance matrix reconstruction according to the calculated steering vector of the interference signal and the minimum characteristic value of the covariance matrix of the initial signal, wherein the expression is as follows:
Figure BDA0002266606050000114
on the basis of the foregoing embodiment of the present disclosure, the step S105 of inputting the steering vector of the desired signal and the wideband interference-plus-noise covariance matrix into a wideband beamformer to obtain a wideband desired signal corresponding to an initial desired signal of the initial signal includes:
calculating a constraint matrix by using the guide vector of the initial expected signal;
calculating a weighting coefficient according to the constraint matrix and the broadband interference and noise covariance matrix;
and performing spatial filtering on the initial signal by using the weighting coefficient to extract a broadband expected signal.
Optionally, the constraint matrix includes steering vectors of the initial desired signal at different frequencies.
For example, a constraint matrix is calculated from the steering vectors of the desired signals
Figure BDA0002266606050000115
Then, the reconstructed interference and noise covariance matrix is used to calculate the weight coefficient of beam forming, and the formula can be expressed as
Figure BDA0002266606050000116
And performing spatial filtering on the initial signal according to the weighting coefficient, extracting and enhancing the initial period signal, suppressing interference signals and noise, and finally outputting a broadband expected signal.
In the above calculation process, as shown in fig. 3, the time delay vector of the signal can be accurately estimated under the condition of the position error of the angular error array element, the estimated time delay vector is used to construct the steering vector of the desired signal and the interference-plus-noise covariance matrix, and then the final weighting vector is obtained to perform spatial filtering on the initial signal, and the beam forming diagram of the spatial filtering result is shown in fig. 4, which can form a main lobe in the direction of the desired signal and a null notch in the direction of the interference, so that the beam former can effectively enhance the desired signal and suppress the interference signal, thereby realizing the extraction of the broadband desired signal. Further, fig. 5 and fig. 6 show a comparison between the performance of the embodiment of the present disclosure and other algorithms, and it can be seen that the output signal to interference plus noise ratio of the method of the present disclosure is always close to an optimal value, and a higher output signal to interference plus noise ratio can be obtained with a small number of fast beats, and the performance is superior to that of other algorithms.
Corresponding to the above method embodiment, referring to fig. 7, the present disclosure also provides a wideband signal extracting apparatus 70, including:
a receiving module 701, configured to receive an initial signal to be extracted, where the initial signal includes an initial desired signal and an interference signal;
a transformation module 702, configured to transform the initial signal to a time-frequency domain by using bilinear time-frequency transformation, so as to obtain a space-time-frequency distribution matrix at each time-frequency point;
a calculating module 703, configured to calculate, according to the space-time frequency distribution matrix, a single source self-entry set of the initial expected signal and a single source self-entry set of the interference signal;
an estimating module 704, configured to estimate a steering vector of the initial desired signal using the set of single-source self-terms of the initial desired signal and the set of single-source self-terms of the interference signal, and reconstruct a wideband interference-plus-noise covariance matrix;
an obtaining module 705, configured to input the steering vector of the initial desired signal and the wideband interference-plus-noise covariance matrix into a wideband beamformer, and obtain a wideband desired signal corresponding to the initial desired signal of the initial signal.
The apparatus shown in fig. 7 may correspondingly execute the content in the above method embodiment, and details of the part not described in detail in this embodiment refer to the content described in the above method embodiment, which is not described again here.
Referring to fig. 8, an embodiment of the present disclosure also provides an electronic device 80, which includes:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the method of extracting a broadband signal of the method embodiments described above.
The disclosed embodiments also provide a non-transitory computer-readable storage medium storing computer instructions for causing a computer to execute the broadband signal extraction method in the foregoing method embodiments.
The disclosed embodiments also provide a computer program product comprising a computer program stored on a non-transitory computer readable storage medium, the computer program comprising program instructions which, when executed by a computer, cause the computer to perform the wideband signal extraction method in the aforementioned method embodiments.
Referring now to FIG. 8, a block diagram of an electronic device 80 suitable for use in implementing embodiments of the present disclosure is shown. The electronic devices in the embodiments of the present disclosure may include, but are not limited to, mobile terminals such as mobile phones, notebook computers, digital broadcast receivers, PDAs (personal digital assistants), PADs (tablet computers), PMPs (portable multimedia players), in-vehicle terminals (e.g., in-vehicle navigation terminals), and the like, and fixed terminals such as digital TVs, desktop computers, and the like. The electronic device shown in fig. 8 is only an example, and should not bring any limitation to the functions and the scope of use of the embodiments of the present disclosure.
As shown in fig. 8, the electronic device 80 may include a processing means (e.g., a central processing unit, a graphics processor, etc.) 801 that may perform various appropriate actions and processes in accordance with a program stored in a Read Only Memory (ROM)802 or a program loaded from a storage means 808 into a Random Access Memory (RAM) 803. In the RAM 803, various programs and data necessary for the operation of the electronic apparatus 80 are also stored. The processing device 801, the ROM 802, and the RAM 803 are connected to each other by a bus 804. An input/output (I/O) interface 805 is also connected to bus 804.
Generally, the following devices may be connected to the I/O interface 805: input devices 806 including, for example, a touch screen, touch pad, keyboard, mouse, image sensor, microphone, accelerometer, gyroscope, or the like; output devices 807 including, for example, a Liquid Crystal Display (LCD), speakers, vibrators, and the like; storage 808 including, for example, magnetic tape, hard disk, etc.; and a communication device 809. The communication means 809 may allow the electronic device 80 to communicate wirelessly or by wire with other devices to exchange data. While the figures illustrate an electronic device 80 having various means, it is to be understood that not all illustrated means are required to be implemented or provided. More or fewer devices may alternatively be implemented or provided.
In particular, according to an embodiment of the present disclosure, the processes described above with reference to the flowcharts may be implemented as computer software programs. For example, embodiments of the present disclosure include a computer program product comprising a computer program embodied on a computer readable medium, the computer program comprising program code for performing the method illustrated in the flow chart. In such an embodiment, the computer program may be downloaded and installed from a network through the communication means 809, or installed from the storage means 808, or installed from the ROM 802. The computer program, when executed by the processing apparatus 801, performs the above-described functions defined in the methods of the embodiments of the present disclosure.
It should be noted that the computer readable medium in the present disclosure can be a computer readable signal medium or a computer readable storage medium or any combination of the two. A computer readable storage medium may be, for example, but not limited to, an electronic, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the foregoing. More specific examples of the computer readable storage medium may include, but are not limited to: an electrical connection having one or more wires, a portable computer diskette, a hard disk, a Random Access Memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the present disclosure, a computer readable storage medium may be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device. In contrast, in the present disclosure, a computer readable signal medium may comprise a propagated data signal with computer readable program code embodied therein, either in baseband or as part of a carrier wave. Such a propagated data signal may take many forms, including, but not limited to, electro-magnetic, optical, or any suitable combination thereof. A computer readable signal medium may also be any computer readable medium that is not a computer readable storage medium and that can communicate, propagate, or transport a program for use by or in connection with an instruction execution system, apparatus, or device. Program code embodied on a computer readable medium may be transmitted using any appropriate medium, including but not limited to: electrical wires, optical cables, RF (radio frequency), etc., or any suitable combination of the foregoing.
The computer readable medium may be embodied in the electronic device; or may exist separately and not be incorporated into the electronic device.
The computer readable medium carries one or more programs which, when executed by the electronic device, cause the electronic device to perform the steps associated with the method embodiments.
Alternatively, the computer readable medium carries one or more programs which, when executed by the electronic device, enable the electronic device to perform the steps associated with the method embodiments.
Computer program code for carrying out operations for aspects of the present disclosure may be written in any combination of one or more programming languages, including an object oriented programming language such as Java, Smalltalk, C + +, and conventional procedural programming languages, such as the "C" programming language or similar programming languages. The program code may execute entirely on the user's computer, partly on the user's computer, as a stand-alone software package, partly on the user's computer and partly on a remote computer or entirely on the remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a Local Area Network (LAN) or a Wide Area Network (WAN), or the connection may be made to an external computer (for example, through the Internet using an Internet service provider).
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present disclosure. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
The units described in the embodiments of the present disclosure may be implemented by software, and may also be implemented by hardware. Where the name of a unit does not in some cases constitute a limitation of the unit itself, for example, the first retrieving unit may also be described as a "unit for retrieving at least two internet protocol addresses".
It should be understood that portions of the present disclosure may be implemented in hardware, software, firmware, or a combination thereof.
The above description is only for the specific embodiments of the present disclosure, but the scope of the present disclosure is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present disclosure should be covered within the scope of the present disclosure. Therefore, the protection scope of the present disclosure shall be subject to the protection scope of the claims.

Claims (11)

1. A method of wideband signal extraction, the method comprising:
receiving an initial signal to be extracted, wherein the initial signal comprises an initial expected signal and an interference signal;
transforming the initial signal to a time-frequency domain by utilizing bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix on each time-frequency point;
calculating according to the space-time frequency distribution matrix to obtain a single-source self-item set of an initial expected signal and a single-source self-item set of an interference signal;
estimating a steering vector of the initial expected signal by using the single-source self-term set of the initial expected signal and the single-source self-term set of the interference signal, and reconstructing a broadband interference and noise covariance matrix;
and inputting the steering vector of the initial expected signal and the broadband interference and noise covariance matrix into a broadband beam former to obtain a broadband expected signal corresponding to the initial expected signal of the initial signal.
2. The method of claim 1, wherein the space-time distribution matrix comprises a non-delayed array received signal vector, a conjugate transpose, and a kernel function.
3. The method according to claim 1, wherein the step of calculating a set of single-source terms of the initial desired signal and a set of single-source terms of the interference signal according to the space-time frequency distribution matrix comprises:
calculating a time delay vector according to the space-time-frequency distribution matrix;
and respectively calculating to obtain a single-source self-item set of the initial expected signal and a single-source self-item set of the interference signal by using the time delay vector and the initial signal.
4. The method of claim 3, wherein the step of calculating the set of terms from which the initial desired signal originates and the set of terms from which the interfering signal originates from using the delay vector and the initial signal respectively comprises:
taking the time delay vector as a feature, extracting a time frequency point corresponding to the initial expected signal from the initial signal to form a time frequency point set corresponding to the initial expected signal, and taking the time frequency points left after the time frequency point set of the initial expected signal is extracted from the initial signal as the time frequency point set of an interference signal and a noise signal;
and removing the mutual item point, the mutual item point and the noise point in the time-frequency point set of the initial expected signal and the time-frequency point set of the interference signal and the noise signal by utilizing the single-source item-derived time-frequency point to obtain the single-source item point set of the initial expected signal and the single-source item point set of the interference signal.
5. The method according to claim 4, wherein the step of extracting time-frequency points corresponding to the initial desired signal from the initial signal as a time-frequency point set of the initial desired signal and using the remaining time-frequency points as a time-frequency point set of interference signals and noise signals by using the delay vector as a feature comprises:
obtaining a calculation result by utilizing the time delay vector and a preset time delay vector of the initial expected signal according to a preset algorithm;
and comparing the calculation result with a first threshold, taking the time frequency points less than or equal to the first threshold as a time frequency point set of the initial expected signal, and taking the time frequency points more than the first threshold as a time frequency point set of interference and noise signals.
6. The method according to claim 4, wherein the step of obtaining the initial desired signal mono-source term set and the interference signal mono-source term set by removing the inter-term point, the inter-term point and the noise point in the time-frequency point set of the initial desired signal and the time-frequency point set of the interference signal and the noise signal by using the mono-source term time-frequency point comprises:
and according to a second threshold value and the maximum characteristic value, eliminating mutual item points, interactive self item points and noise points from the time frequency point set of the initial expected signal to obtain a single-source self item point set of the expected signal, and eliminating mutual item points, interactive self item points and noise points from the time frequency point set of the interference signal and the noise signal to obtain a single-source self item point set of the interference signal.
7. The method of claim 1, wherein the step of estimating a steering vector of the initial desired signal using the set of single-sourced terms of the initial desired signal and the set of single-sourced terms of the interfering signal, and reconstructing a wideband interference-plus-noise covariance matrix comprises:
taking the time delay vector as a characteristic, and dividing the single-source self-item set of the interference signal into a set of single-source self-items of each interference signal by using a mean value clustering method;
calculating a steering vector of the initial signal on different frequencies by using the single-origin-term set of the initial expected signal and the single-origin-term set of each interference signal, wherein the steering vector comprises the steering vector of the initial expected signal and the steering vector of the interference signal;
and reconstructing the broadband interference and noise covariance matrix according to the guide vector of the interference signal.
8. The method of claim 7, wherein the step of inputting the steering vector of the initial desired signal and the wideband interference-plus-noise covariance matrix to a wideband beamformer to obtain a wideband desired signal corresponding to the initial desired signal of the initial signals comprises:
calculating a constraint matrix by using the guide vector of the initial expected signal;
calculating a weighting coefficient according to the constraint matrix and the broadband interference and noise covariance matrix;
and performing spatial filtering on the initial signal by using the weighting coefficient to extract a broadband expected signal.
9. The method of claim 8, wherein the constraint matrix comprises steering vectors of the initial desired signal at different frequencies.
10. A wideband signal extraction apparatus, comprising:
the device comprises a receiving module, a processing module and a processing module, wherein the receiving module is used for receiving an initial signal to be extracted, and the initial signal comprises an initial expected signal and an interference signal;
the transformation module is used for transforming the initial signal to a time-frequency domain by utilizing bilinear time-frequency transformation to obtain a space-time-frequency distribution matrix on each time-frequency point;
the calculation module is used for calculating to obtain a single-source self-item set of an initial expected signal and a single-source self-item set of an interference signal according to the space-time frequency distribution matrix;
an estimation module, configured to estimate a steering vector of the initial desired signal using the set of single-source self-terms of the initial desired signal and the set of single-source self-terms of the interference signal, and reconstruct a wideband interference-plus-noise covariance matrix;
and the acquisition module is used for inputting the steering vector of the initial expected signal and the broadband interference and noise covariance matrix into a broadband beam former to obtain a broadband expected signal corresponding to the initial expected signal of the initial signal.
11. An electronic device, characterized in that the electronic device comprises:
at least one processor; and the number of the first and second groups,
a memory communicatively coupled to the at least one processor; wherein the content of the first and second substances,
the memory stores instructions executable by the at least one processor to enable the at least one processor to perform the wideband signal extraction method of any of the preceding claims 1-9.
CN201911090146.9A 2019-11-08 2019-11-08 Broadband signal extraction method and device and electronic equipment Active CN111245490B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911090146.9A CN111245490B (en) 2019-11-08 2019-11-08 Broadband signal extraction method and device and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911090146.9A CN111245490B (en) 2019-11-08 2019-11-08 Broadband signal extraction method and device and electronic equipment

Publications (2)

Publication Number Publication Date
CN111245490A true CN111245490A (en) 2020-06-05
CN111245490B CN111245490B (en) 2023-09-01

Family

ID=70864013

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911090146.9A Active CN111245490B (en) 2019-11-08 2019-11-08 Broadband signal extraction method and device and electronic equipment

Country Status (1)

Country Link
CN (1) CN111245490B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112073076A (en) * 2020-07-31 2020-12-11 南京航空航天大学 Space-frequency adaptive processing anti-interference method based on short-time Fourier transform

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070164902A1 (en) * 2005-12-02 2007-07-19 Samsung Electronics Co., Ltd. Smart antenna beamforming device in communication system and method thereof
CN103091661A (en) * 2013-02-01 2013-05-08 西安科技大学 Broadband signal arriving direction estimation method based on iteration spectral reconfiguration
CN107167809A (en) * 2017-06-14 2017-09-15 哈尔滨工程大学 It is a kind of that array beamforming method is blocked based on the broadband that signal subspace is focused on
CN108880586A (en) * 2018-06-28 2018-11-23 中国人民解放军战略支援部队信息工程大学 A kind of broadband weak signal enhancement method and apparatus
CN108872946A (en) * 2018-04-20 2018-11-23 西安电子科技大学 The robust ada- ptive beamformer method of steering vector and covariance matrix Joint iteration
CN109743078A (en) * 2019-01-11 2019-05-10 河南宏泰控飞信息技术有限公司 Unmanned plane frequency hopping remote signal intercept receiving method and equipment based on array antenna

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20070164902A1 (en) * 2005-12-02 2007-07-19 Samsung Electronics Co., Ltd. Smart antenna beamforming device in communication system and method thereof
CN103091661A (en) * 2013-02-01 2013-05-08 西安科技大学 Broadband signal arriving direction estimation method based on iteration spectral reconfiguration
CN107167809A (en) * 2017-06-14 2017-09-15 哈尔滨工程大学 It is a kind of that array beamforming method is blocked based on the broadband that signal subspace is focused on
CN108872946A (en) * 2018-04-20 2018-11-23 西安电子科技大学 The robust ada- ptive beamformer method of steering vector and covariance matrix Joint iteration
CN108880586A (en) * 2018-06-28 2018-11-23 中国人民解放军战略支援部队信息工程大学 A kind of broadband weak signal enhancement method and apparatus
CN109743078A (en) * 2019-01-11 2019-05-10 河南宏泰控飞信息技术有限公司 Unmanned plane frequency hopping remote signal intercept receiving method and equipment based on array antenna

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112073076A (en) * 2020-07-31 2020-12-11 南京航空航天大学 Space-frequency adaptive processing anti-interference method based on short-time Fourier transform
CN112073076B (en) * 2020-07-31 2021-10-22 南京航空航天大学 Space-frequency adaptive processing anti-interference method based on short-time Fourier transform

Also Published As

Publication number Publication date
CN111245490B (en) 2023-09-01

Similar Documents

Publication Publication Date Title
WO2020207174A1 (en) Method and apparatus for generating quantized neural network
US11222646B2 (en) Apparatus and method for generating audio signal with noise attenuated based on phase change rate
CN107240396B (en) Speaker self-adaptation method, device, equipment and storage medium
CN109756818B (en) Dual-microphone noise reduction method and device, storage medium and electronic equipment
EP4266308A1 (en) Voice extraction method and apparatus, and electronic device
CN111308414A (en) Method and system for estimating direction of arrival, intelligent terminal and storage medium
CN110188782B (en) Image similarity determining method and device, electronic equipment and readable storage medium
CN110009101B (en) Method and apparatus for generating a quantized neural network
CN115208720A (en) Electronic device
CN111245490B (en) Broadband signal extraction method and device and electronic equipment
CN111211821B (en) System and method for analog beamforming for a single connection antenna array
CN112087235A (en) Sparsity self-adaptive DOA estimation method and system based on pseudo-inverse perception dictionary
CN110555861A (en) optical flow calculation method and device and electronic equipment
WO2022222922A1 (en) Voice signal processing method and apparatus
CN111626044B (en) Text generation method, text generation device, electronic equipment and computer readable storage medium
CN111755021B (en) Voice enhancement method and device based on binary microphone array
CN113808606B (en) Voice signal processing method and device
US20240046948A1 (en) Device and method for removing wind noise and electronic device comprising wind noise removing device
US10911096B1 (en) System and method for providing sub-band whitening in the presence of partial-band interference
KR20190086339A (en) Apparatus and method of non-iterative singular-value decomposition
CN111090758B (en) Media data processing method, device and storage medium
CN110263797B (en) Method, device and equipment for estimating key points of skeleton and readable storage medium
US11069373B2 (en) Speech processing method, speech processing apparatus, and non-transitory computer-readable storage medium for storing speech processing computer program
CN116342785A (en) Image processing method, device, equipment and medium
CN116167933A (en) Image noise reduction method, device, electronic equipment and computer readable medium

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant