Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Fig. 1 is a flowchart of a first embodiment of a signal collision separation method according to the present invention, and as shown in fig. 1, the method of this embodiment may include:
step 101, receiving an aliasing signal formed by signal collision;
the present invention can be based on various forms of signals to be separated, and the present embodiment takes the separation of the ADS-B signal as an example, but does not limit the present invention.
The duration of an ADS-B baseband signal is 120 mus, wherein the first 8 mus is a leading pulse which is composed of signals with fixed formats and is 8bit, the leading pulse is composed of four pulses which are respectively present at 0 mus, 1 mus, 3.5 mus and 4.5 mus, and the leading pulse is used for a receiver to find signals and a synchronous sampling clock; next 112 μ s is the data portion and the ADS-B signal data portion is 112 bits. The ADS-B signal is Manchester encoded, i.e., symbol 0 is encoded as [0,1], symbol 1 is encoded as [1,0], and thus the total code length is 240 bits.
An ADS-B baseband signal can be expressed as
Wherein b [ n ]]B is the nth code element of B, B is a sequence formed by Manchester encoding ADS-B signals and containing 240 code elements, p (T) is a square wave pulse function, and the pulse width Ts=0.5μs。
In this embodiment, the received alias signal is formed by a plurality of ADS-B baseband signals s (t) after being subjected to up-conversion and then colliding at a receiving end.
In order to utilize the difference brought by the signals of different incident angles to the output, the array antenna is adopted at the receiving front end of the receiver to receive the aliasing signals, the specific form of the array antenna is not specially required, in the embodiment, the uniform linear array antenna is taken as an example for explanation, the multipath interference of the signals can be ignored, and the antenna units of the uniform linear array antenna are separated by half wavelength.
102, whitening the aliasing signals to obtain the number of signal collisions, and dividing a plurality of time periods according to the number of the signal collisions;
and whitening the received aliasing signals, estimating the number of signal collisions in the period of time, and dividing the period of time into a plurality of relatively short time periods according to the number of the signal collisions.
103, sequentially performing matrix transformation on corresponding signal matrixes in each time period to obtain a beam forming matrix, and obtaining separated signals according to the beam forming matrix;
the received signal can be expressed as the matrix product of the beam forming matrix of the array antenna and the transmitted signal, so that the separated transmitted signal can be obtained only by estimating the beam forming matrix of the array antenna and filtering the received signal by adopting the beam forming matrix.
And if the time periods obtained by the division in the step 102 are relatively short, matrix transformation is performed on the corresponding signal matrix to obtain a beam forming matrix, so that the separated signals can be obtained.
In this embodiment, the number of signal collisions is estimated by performing whitening processing on received aliasing signals, a plurality of time periods are divided according to the number of signal collisions, matrix transformation is performed on corresponding signal matrices in each time period in sequence to obtain a beamforming matrix, and separated signals are obtained according to the beamforming matrix, so that effective separation of a plurality of collision signals is achieved. The problem of the prior art when separating the condition of many signal collisions algorithm complexity is too high is solved.
On the basis of the foregoing embodiment, the following embodiment further provides a process of performing whitening processing on the aliasing signal in step 102 to obtain the number of signal collisions, and dividing a plurality of time segments according to the number of signal collisions, which includes the following specific steps:
and N-point sampling is carried out on the received aliasing signals within the time T, the sampling frequency meets the Nyquist law, and the array antenna for receiving the aliasing signals is provided with m antenna units, so that a matrix with the dimension of m multiplied by N is obtained.
The received signal X is a matrix of dimensions m × N, as follows:
the relationship between the received signal X and the original signal S can be expressed as:
X=MS+Z
wherein, M is a beamforming matrix with dimension of M multiplied by d, S is an original signal with dimension of d multiplied by N, and Z is a noise matrix with dimension of M multiplied by N.
The time T is divided evenly into n segments, i.e. T ═ n × Δ T, and a signal matrix X corresponding to the received signal X in each Δ T segmentΔT,iAnd (i-1, 2, … …, n) performing singular value decomposition, wherein the number of singular values larger than a preset threshold is the number of signals in the time period delta T, the maximum value of the number of n time period delta T signals is the number d of signal collisions in the time period T, and the number d of signal collisions is less than or equal to the number m of antenna elements of the array antenna.
Dividing adjacent delta T with the same number of signals into the same time segment, and dividing the time segment (T) in the time T1,t2),(t2,t3),(t3,t4),……,(td,td+1) (ii) a The number of signals included in each time segment is 1,2,3, … …, d in order of time segment. I.e. the d signals of a collision arrive one after the other. For time domain segmentation of the signal, reference may be made to the schematic diagram of time domain segmentation of the signal shown in fig. 3.
The embodiment divides the received signal to obtain a plurality of relatively small signal matrixes XΔT,iAnd (i is 1,2, … …, n) singular value decomposition is carried out, so that the complexity of the signal whitening processing is reduced. And the number of collision signals in each segment can be estimated according to the relationship between the singular value obtained by decomposition and a preset threshold value, so that the number d of signal collisions in the time T is estimated and the time T is segmented. The time domain segmentation is carried out on the signals, and a larger signal matrix is divided into a plurality of smaller signal matrices for matrix operation respectively, so that the time domain segmentation can be greatly carried outThe algorithm complexity of the matrix operation is reduced.
Fig. 2 is a flowchart of a third embodiment of the signal collision separation method of the present invention, and on the basis of the above embodiment, this embodiment further provides specific steps of performing matrix transformation on the corresponding signal matrix in each time period in sequence in step 103 to obtain a beamforming matrix, and obtaining separated signals according to the beamforming matrix, where the specific steps are as follows:
step 1031:
from the above time-domain segmentation of the received signal, at (t)1,t2) Contains only 1 signal in a time period during which the signal is clean, pair (t)1,t2) Signal matrix X corresponding to time periodt1Singular value decomposition is carried out, namely:
wherein, Xt1Is mxL1Matrix of dimensions, L1=fs×(t2-t1),fsIs the sampling frequency; u shape1Is a m x m dimensional square matrix, V1Is L1×L1Dimensional matrix, X1Is mxL1A matrix of dimensions.
Retaining only X1Setting the rest values as 0 to remove the singular value, namely removing the influence of noise and interference to obtain X'1Is thus
Through to X't1Is subjected to rank-one decomposition to obtain
X′t1=m1st1
m1Is a matrix of m x 1 dimensions, st1Is 1 XL1And (4) matrix. m is1I.e. the first column of the desired beamforming matrix.
Step 1032:
the construction matrix P ═ w
1,…,w
m-1]Wherein w is
iAll satisfy
I.e. w
iAnd m
1Are orthogonal. The received signal X is filtered using the following equation,
from the above formula, the signal X(1)Is the result of separating the first signal from X and no longer contains any information of the first signal. And after the conversion, the received signal X is converted into a matrix X of (m-1) xN dimensions(1)The matrix dimension is reduced from m to m-1, and m'2=PHm2Is a (m-1) × 1 dimensional matrix.
Step 1033:
due to X
(1)No longer contains any information of the first signal, so X
(1)At (t)
2,t
3) The time period is pure and contains only one signal. To X
(1)At (t)
2,t
3) Corresponding signal matrix in time period
Repeating the
steps 1031 and 1032 to obtain m'
2And a further dimension-reduced (m-2) X N-dimensional matrix X no longer containing any information of the first signal and the second signal
(2)Is through to m'
2=P
Hm
2The second column m of the beam forming matrix can be obtained by inverse transformation
2。
And sequentially iterating until the d-th column of the beam forming matrix is obtained.
Step 1034:
forming beam forming matrix M ═ M1,m2,…,md]And finding the generalized inverse M of M+By S ═ M+And filtering the received signal X by the X to obtain a separation signal S.
In this embodiment, each column of the beamforming matrix is sequentially obtained in an iterative manner, and in the iterative process, the dimension of the matrix is reduced by one dimension once per iteration, thereby further reducing the complexity of the algorithm. The signal collision separation method provided by the embodiment can meet the requirement on the operation complexity when the satellite-based ADS-B receiver separates the multi-signal collision condition.
FIG. 4 is a schematic diagram of a bit error rate simulation result of the signal collision separation method of the present invention; FIG. 5 is a schematic diagram of simulation results of signal-to-noise ratio output by the signal collision separation method according to the present invention; FIG. 6 is a schematic diagram of a simulation result of algorithm complexity of the signal collision separation method of the present invention. Wherein, AZCMA represents a constant coefficient constant separation algorithm, MDA represents a manchester decoding algorithm, and IPA represents an algorithm corresponding to an embodiment of the present invention. Fig. 4 shows the error rates output by three different algorithms with different signal-to-noise ratios, and it can be seen from the figure that the error rate performance of the IPA algorithm is similar to that of the AZCMA algorithm and better than that of the MDA algorithm. Fig. 5 shows the output snr of the three algorithms under different snr environments, and it can be seen from the graph that the IPA algorithm performs similar to the AZCMA algorithm in performance better than the MDA algorithm under the condition of low snr. In an environment with low signal-to-noise ratio of satellite-based ADS-B, the IPA algorithm is more suitable for ADS-B signal separation in the environment than the AZCMA algorithm and the IPA algorithm. Fig. 6 shows the computation time required by the three algorithms to separate different numbers of collision signals, and it can be seen from the graph that the computation amounts of the three algorithms are not obviously different when the number of collision signals is small, but when the number of collision signals is large, the computation amount is obviously reduced by the IPA algorithm compared with the other two algorithms, and the IPA algorithm is more suitable for a common high-order collision scenario in the satellite-based ADS-B reception. In conclusion, the signal collision separation method can greatly reduce the operation complexity of the algorithm and shorten the time required for separating collision signals under the condition of not losing the performance, and is more suitable for the satellite-based ADS-B receiver with limited resources and capacity of operation and processing.
Fig. 7 is a schematic structural diagram of a first embodiment of a signal collision separation apparatus according to the present invention, and as shown in fig. 7, an apparatus 20 of the present embodiment may include: a receiving module 201, a whitening module 202 and a separating module 203.
A receiving module 201, configured to receive an alias signal formed by signal collision;
in order to utilize the difference brought by the signals of different incident angles to the output, the receiving module adopts the array antenna to receive the aliasing signals, the specific form of the array antenna is not required to be specially required, in the embodiment, the uniform linear array antenna is taken as an example for explanation, and the antenna units of the uniform linear array antenna are separated by half wavelength.
The whitening module 202 is configured to perform whitening processing on the aliasing signals to obtain the number of signal collisions, and divide a plurality of time periods according to the number of signal collisions;
the separation module 203 is configured to perform matrix transformation on the corresponding signal matrix in each time period in sequence to obtain a beamforming matrix, and obtain a separated signal according to the beamforming matrix.
The apparatus of this embodiment may be used to implement the technical solution of the method embodiment shown in fig. 1, and the implementation principle and the technical effect are similar, which are not described herein again.
The embodiment is implemented on the basis of the embodiment of fig. 7, and specifically includes the following steps:
optionally, the whitening module is specifically configured to:
and N-point sampling is carried out on the received aliasing signals within the time T, the sampling frequency meets the Nyquist law, and the array antenna for receiving the aliasing signals is provided with m antenna units, so that a matrix with the dimension of m multiplied by N is obtained.
The time T is divided evenly into n segments, i.e. T ═ n × Δ T, and a signal matrix X corresponding to the received signal X in each Δ T segmentΔT,iAnd (i-1, 2, … …, n) performing singular value decomposition, wherein the number of singular values larger than a preset threshold is the number of signals in the time period delta T, the maximum value of the number of n time period delta T signals is the number d of signal collisions in the time period T, and the number d of signal collisions is less than or equal to the number m of antenna elements of the array antenna.
Dividing adjacent delta T with the same number of signals into the same time segment, and dividing the time segment (T) in the time T1,t2),(t2,t3),(t3,t4),……,(td,td+1) (ii) a The number of signals included in each time segment is 1,2,3, … …, d in order of time segment. I.e. the d signals of a collision arrive one after the other.
Optionally, the separation module is specifically configured to:
separating the collision signal to obtain a separated signal S according to the following steps:
step 1:
from the above time-domain segmentation of the received signal, at (t)1,t2) Contains only 1 signal in a time period during which the signal is clean, pair (t)1,t2) Signal matrix X corresponding to time periodt1Singular value decomposition is carried out, namely:
wherein, Xt1Is mxL1Matrix of dimensions, L1=fs×(t2-t1),fsIs the sampling frequency; u shape1Is a m x m dimensional square matrix, V1Is L1×L1Dimensional matrix, X1Is mxL1A matrix of dimensions.
Retaining only X1Setting the rest values as 0 to remove the singular value, namely removing the influence of noise and interference to obtain X'1Is thus
Through to X't1Is subjected to rank-one decomposition to obtain
X′t1=m1st1
m1Is a matrix of m x 1 dimensions, st1Is 1 XL1And (4) matrix. m is1I.e. the first column of the desired beamforming matrix.
Step 2:
the construction matrix P ═ w1,…,wm-1]Wherein w isiAll satisfyI.e. wiAnd m1Are orthogonal. The received signal X is filtered using the following equation,
from the above formula, the signal X(1)Is the result of separating the first signal from X and no longer contains any information of the first signal. And after the conversion, the received signal X is converted into a matrix X of (m-1) xN dimensions(1)The matrix dimension is reduced from m to m-1, and m'2=PHm2Is a (m-1) × 1 dimensional matrix.
And step 3:
due to X
(1)No longer contains any information of the first signal, so X
(1)At (t)
2,t
3) The time period is pure and contains only one signal. To X
(1)At (t)
2,t
3) Corresponding signal matrix in time period
Repeating the
step 1 and the step 2 to obtain m'
2And a further dimension-reduced (m-2) X N-dimensional matrix X no longer containing any information of the first signal and the second signal
(2)Is through to m'
2=P
Hm
2The second column m of the beam forming matrix can be obtained by inverse transformation
2。
And sequentially iterating until the d-th column of the beam forming matrix is obtained.
And 4, step 4:
forming beam forming matrix M ═ M1,m2,…,md]And finding the generalized inverse M of M+By S ═ M+And filtering the received signal X by the X to obtain a separation signal S.
In the several embodiments provided in the present application, it should be understood that the disclosed method and apparatus may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, a division of modules is merely a division of logical functions, and an actual implementation may have another division, for example, a plurality of modules or components may be combined or integrated into another system, or some features may be omitted, or not executed.
Those of ordinary skill in the art will understand that: all or a portion of the steps of implementing the above-described method embodiments may be performed by hardware associated with program instructions. The program may be stored in a computer-readable storage medium. When executed, the program performs steps comprising the method embodiments described above; and the aforementioned storage medium includes: various media that can store program codes, such as ROM, RAM, magnetic or optical disks.
Finally, it should be noted that: the above embodiments are only used to illustrate the technical solution of the present invention, and not to limit the same; while the invention has been described in detail and with reference to the foregoing embodiments, it will be understood by those skilled in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some or all of the technical features may be equivalently replaced; and the modifications or the substitutions do not make the essence of the corresponding technical solutions depart from the scope of the technical solutions of the embodiments of the present invention.