CN107276937B - Signal collision separation method and device - Google Patents

Signal collision separation method and device Download PDF

Info

Publication number
CN107276937B
CN107276937B CN201710591476.0A CN201710591476A CN107276937B CN 107276937 B CN107276937 B CN 107276937B CN 201710591476 A CN201710591476 A CN 201710591476A CN 107276937 B CN107276937 B CN 107276937B
Authority
CN
China
Prior art keywords
signal
matrix
time
signals
beam forming
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.)
Active
Application number
CN201710591476.0A
Other languages
Chinese (zh)
Other versions
CN107276937A (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.)
Tianyu Aviation Data Technology Hefei Co ltd
Original Assignee
Beihang University
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 Beihang University filed Critical Beihang University
Priority to CN201710591476.0A priority Critical patent/CN107276937B/en
Publication of CN107276937A publication Critical patent/CN107276937A/en
Application granted granted Critical
Publication of CN107276937B publication Critical patent/CN107276937B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/03Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
    • H04L25/03891Spatial equalizers
    • H04L25/03898Spatial equalizers codebook-based design
    • H04L25/0391Spatial equalizers codebook-based design construction details of matrices
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Power Engineering (AREA)
  • Radio Relay Systems (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

The invention provides a signal collision separation method and a signal collision separation device. The invention provides a signal collision separation method, which comprises the following steps: receiving an aliasing signal formed by signal collision; 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 sequentially carrying out matrix transformation on the corresponding signal matrix in each time period to obtain a beam forming matrix, and obtaining separated signals according to the beam forming matrix. The signal collision separation method and the signal collision separation device provided by the invention reduce the complexity of an algorithm when the multi-signal collision condition is separated.

Description

Signal collision separation method and device
Technical Field
The invention relates to the aviation monitoring technology, in particular to a signal collision separation method and a signal collision separation device.
Background
The satellite-based Broadcast-based Automatic Dependent-Broadcast (ADS-B) technology is a novel aviation monitoring technology which is used for carrying an ADS-B receiver as a part of satellite load on a satellite platform and realizing the large-scale monitoring of aviation nodes by utilizing the large-scale coverage capability of the satellite. Because the coverage area of the satellite node is relatively wide, a plurality of airplanes often run simultaneously in the coverage area. When each aircraft broadcasts its own ADS-B signal, the ADS-B receiver will receive an aliased signal formed by many ADS-B signals colliding together. Proper separation of the aliased signals is required for the ADS-B receiver to obtain ADS-B information for each aircraft.
The existing projection algorithm is based on the fact that only two signals collide, the single-channel condition is converted into the multi-channel condition to separate collision signals, and therefore only the two signal collision conditions can be separated; the existing blind signal separation algorithm (such as constant coefficient constant separation algorithm and Manchester decoding algorithm) takes the zero return characteristic of a sampling point of a source signal as a constraint condition, separates the collided aliasing signals based on a solution matrix equation, and increases the operation complexity index of the algorithm along with the increase of the number of signal collisions.
The satellite-based ADS-B receiver has limited resources and capabilities for operation and processing, and needs an effective algorithm which can separate the collision situation of multiple signals and has low operation complexity.
Disclosure of Invention
The embodiment of the invention provides a signal collision separation method and a signal collision separation device, which are used for reducing the complexity of an algorithm when separation is carried out on the condition of multi-signal collision.
In a first aspect, the present invention provides a signal collision separation method, including:
receiving an aliasing signal formed by signal collision;
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 sequentially carrying out matrix transformation on the corresponding signal matrix in each time period to obtain a beam forming matrix, and obtaining separated signals according to the beam forming matrix.
In a second aspect, the present invention provides a signal collision separating apparatus comprising:
the receiving module is used for receiving an aliasing signal formed by signal collision;
the whitening module is used for 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 the separation module is used for carrying out matrix transformation on the corresponding signal matrix in each time period in sequence to obtain a beam forming matrix and obtaining the separated signals according to the beam forming matrix.
The invention estimates the number of signal collisions by whitening the received aliasing signals, divides a plurality of time periods according to the number of the signal collisions, sequentially carries out matrix transformation on corresponding signal matrixes in each time period to obtain a beamforming matrix, and obtains signals after separation according to the beamforming matrix, thereby realizing effective separation of a plurality of collision signals.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the description of the embodiments or the prior art will be briefly introduced below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to these drawings without creative efforts.
FIG. 1 is a flow chart of a first embodiment of a signal collision separation method according to the present invention;
FIG. 2 is a diagram illustrating time-domain segmentation of a signal according to a second embodiment of the signal collision separation method of the present invention;
FIG. 3 is a flow chart of a third embodiment of a signal collision separation method according to the present invention;
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;
fig. 7 is a schematic structural diagram of a signal collision separating apparatus according to a first embodiment of the present invention.
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
Figure BDA0001354850010000031
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:
Figure BDA0001354850010000051
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
Figure BDA0001354850010000052
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 ═ w1,…,wm-1]Wherein w isiAll satisfy
Figure BDA0001354850010000061
I.e. wiAnd m1Are orthogonal. The received signal X is filtered using the following equation,
Figure BDA0001354850010000062
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,t3) The time period is pure and contains only one signal. To X(1)At (t)2,t3) Corresponding signal matrix in time period
Figure BDA0001354850010000063
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=PHm2The second column m of the beam forming matrix can be obtained by inverse transformation2
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:
Figure BDA0001354850010000081
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
Figure BDA0001354850010000082
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,
Figure BDA0001354850010000092
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,t3) The time period is pure and contains only one signal. To X(1)At (t)2,t3) Corresponding signal matrix in time period
Figure BDA0001354850010000093
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=PHm2The second column m of the beam forming matrix can be obtained by inverse transformation2
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.

Claims (10)

1. A signal collision separation method, comprising:
receiving an aliasing signal formed by signal collision;
whitening the aliasing signal to obtain the number of signal collisions, and dividing a plurality of time periods according to the number of the signal collisions;
sequentially carrying out matrix transformation on the corresponding signal matrix in each time period to obtain a beam forming matrix, and obtaining separated signals according to the beam forming matrix;
the whitening processing the aliasing signal to obtain the number of signal collisions specifically includes:
sampling the aliasing signals at N points within time T to obtain received signals X, wherein X is an m multiplied by N dimensional matrix, and m is the number of antenna units of an array antenna for receiving the aliasing signals;
uniformly dividing the time T into a plurality of sections, performing singular value decomposition on a signal matrix corresponding to the received signal X in each time section, taking the number of singular values larger than a preset threshold value as the number of signals in each time section, and taking the maximum value of the number of the signals in the time sections as the number d of signal collisions in the time T;
in each time period, matrix transformation is performed on the corresponding signal matrix in sequence to obtain a beam forming matrix, and separated signals are obtained according to the beam forming matrix, which specifically comprises:
and in each time period in turn, respectively executing the following steps on the corresponding signal matrix:
step 1, matrix transformation is carried out to obtain a first column m of the beam forming matrix1
Step 2, constructing a matrix P ═ w1,…,wm-1]Wherein w isiAll satisfyBy PHFiltering the received signal to obtain X(1)
Step 3, for X(1)Repeating the step 1 and the step 2 on the matrix corresponding to the next time period to obtain a second column m of the beamforming matrix2And X(2)(ii) a Sequentially iterating until the d-th column of the beam forming matrix is obtained;
step 4, constructing the beam forming matrix M ═ M1,m2,…,md]Finding the generalized inverse M of M+By M+And filtering the received signal to obtain the separated signal.
2. The method according to claim 1, wherein the dividing a plurality of time segments according to the number of signal collisions specifically comprises:
dividing the adjacent time segments 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 in the time T is 1,2,3, … …, and d in sequence according to the sequence of the time segments.
3. The method according to claim 1 or 2, wherein the matrix transformation is performed to obtain the first column m of the beamforming matrix1The method specifically comprises the following steps:
for the corresponding signal X in the time periodt1Performing singular value decomposition, i.e.
Figure FDA0002289022590000021
Retention of X1Setting the rest values as 0 to remove the singular value, namely removing the influence of noise and interference to obtain X'1By passing
Figure FDA0002289022590000022
Construction X't1To X't1Performing a rank-one decomposition, i.e.
X′t1=m1st1
Wherein, Xt1Is mxL1Matrix of dimensions, U1Is a m x m dimensional square matrix, V1Is L1×L1Dimensional matrix, st1Is 1 XL1Matrix, m1I.e. the first column of the beamforming matrix.
4. Method according to claim 1 or 2, characterized in that said treatment with P is carried outHFiltering the received signal line to obtain X(1)The method specifically comprises the following steps:
by X(1)=PHX is calculated(1)
5. Method according to claim 1 or 2, characterized in that said treatment with M is carried out+Filtering the received signal to obtain a separated signal S, specifically including:
by S ═ M+X calculates the separation signal S.
6. A signal collision separating apparatus, comprising:
the receiving module is used for receiving an aliasing signal formed by signal collision;
the whitening module is used for 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;
the separation module is used for carrying out matrix transformation on the corresponding signal matrix in each time period in sequence to obtain a beam forming matrix and obtaining separated signals according to the beam forming matrix;
the whitening module is specifically configured to:
sampling the aliasing signals at N points within time T to obtain received signals X, wherein X is an m multiplied by N dimensional matrix, and m is the number of antenna units of an array antenna for receiving the aliasing signals;
uniformly dividing the time T into a plurality of sections, performing singular value decomposition on a signal matrix corresponding to the received signal X in each time section, taking the number of singular values larger than a preset threshold value as the number of signals in each time section, and taking the maximum value of the number of the signals in the time sections as the number d of signal collisions in the time T;
the separation module is specifically configured to:
sequentially and respectively executing the following steps on the corresponding signal matrix in each time period:
step 1, matrix transformation is carried out to obtain a first column m of the beam forming matrix1
Step 2, constructing a matrixP=[w1,…,wm-1]Wherein w isiAll satisfy
Figure FDA0002289022590000031
By PHFiltering the received signal to obtain X(1)
Step 3, for X(1)Repeating the step 1 and the step 2 on the matrix corresponding to the next time period to obtain a second column m of the beamforming matrix2And X(2)(ii) a Sequentially iterating until the d-th column of the beam forming matrix is obtained;
step 4, constructing the beam forming matrix M ═ M1,m2,…,md]Finding the generalized inverse M of M+By M+And filtering the received signal to obtain the separated signal.
7. The apparatus of claim 6, wherein the whitening module is further configured to:
dividing the adjacent time segments 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 in the time T is 1,2,3, … …, and d in sequence according to the sequence of the time segments.
8. The apparatus according to claim 6 or 7, wherein the separation module is further configured to:
for the corresponding signal X in the time periodt1Performing singular value decomposition, i.e.
Retention of X1Setting the rest values as 0 to remove the singular value, namely removing the influence of noise and interference to obtain X'1By passing
Figure FDA0002289022590000033
Construction X't1To X't1Performing a rank-one decomposition, i.e.
X′t1=m1st1
Wherein, Xt1Is mxL1Matrix of dimensions, U1Is a m x m dimensional square matrix, V1Is L1×L1Dimensional matrix, st1Is 1 XL1Matrix, m1I.e. the first column of the beamforming matrix.
9. The apparatus of claim 6 or 7, wherein the separation module is further configured to:
according to X(1)=PHX is calculated(1)
10. The apparatus of claim 6 or 7, wherein the separation module is further configured to:
according to S ═ M+X calculates the separation signal S.
CN201710591476.0A 2017-07-19 2017-07-19 Signal collision separation method and device Active CN107276937B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710591476.0A CN107276937B (en) 2017-07-19 2017-07-19 Signal collision separation method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710591476.0A CN107276937B (en) 2017-07-19 2017-07-19 Signal collision separation method and device

Publications (2)

Publication Number Publication Date
CN107276937A CN107276937A (en) 2017-10-20
CN107276937B true CN107276937B (en) 2020-02-18

Family

ID=60078291

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710591476.0A Active CN107276937B (en) 2017-07-19 2017-07-19 Signal collision separation method and device

Country Status (1)

Country Link
CN (1) CN107276937B (en)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004088878A1 (en) * 2003-03-31 2004-10-14 Matsushita Electric Industrial Co., Ltd. Radio communication method and radio communication device

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103957182B (en) * 2014-05-21 2017-03-08 苏州东奇信息科技股份有限公司 A kind of common channel full duplex system based on MPPSK modulation
US10264407B2 (en) * 2015-06-25 2019-04-16 The Board Of Trustees Of The University Of Alabama Intelligent multi-bean medium access control in ku-band for mission-oriented mobile mesh networks
CN106301687B (en) * 2016-08-19 2019-12-10 上海交通大学 satellite-borne AIS receiving method based on 4-bit CRC error correction and iterative interference elimination algorithm

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2004088878A1 (en) * 2003-03-31 2004-10-14 Matsushita Electric Industrial Co., Ltd. Radio communication method and radio communication device

Also Published As

Publication number Publication date
CN107276937A (en) 2017-10-20

Similar Documents

Publication Publication Date Title
US10652067B2 (en) Orthogonal differential vector signaling codes with embedded clock
Schniter Turbo reconstruction of structured sparse signals
CN110753937A (en) Data transmission network configuration
JP6506767B2 (en) Receiver, frame synchronization method, transmitter, transmission method, and program
CN106662645B (en) Adaptive selection STAP (static adaptive station) method for beam-Doppler channel
US10998885B2 (en) Method for equivalent high sampling rate FIR filtering based on FPGA
EP2242196A1 (en) Reception device, reception method, and communication system
US9729160B1 (en) Wideband analog to digital conversion by random or level crossing sampling
CN107037409B (en) MIMO radar waveform separation method based on compressed sensing
US20170064701A1 (en) Message coding for ncma-based multiple access networks
US8837572B2 (en) Receiver and method for equalizing signals
EP3001570B1 (en) Turbo equalization method and turbo equalization system
CN107276937B (en) Signal collision separation method and device
JP2013509753A (en) Channel estimation method and apparatus in multi-antenna system
GB2452319A (en) MIMO system with interpolation of precoder matrices from a subset of subcarriers
EP3490177A1 (en) Fast decoding method and device suitable for ovxdm system, and ovxdm system
JP2009055217A (en) Signal detector, signal detection method, program thereof, and recording medium
Mourad et al. Orthogonal matching pursuit with correction
CN102137050A (en) Processing method and equipment for data signal
KR101477664B1 (en) Method and apparatus for network coding using compressed sensing technique, method and apparatus for network decoding
KR101955102B1 (en) Method and Apparatus for Estimating Incidence Angle of Received Signal based on codebook in Multiple-Input Multiple-Output Communication Systems
Kaplan et al. A message transmission scheme for linear time-varying multipath channels
EP1578024A1 (en) Joint detection method using MMSE-LE for CDMA-based mobile communication systems
KR102246472B1 (en) Method for radar electrtonic surveillance based on mwc
KR102246461B1 (en) System for radar electrtonic surveillance based on mwc

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
TR01 Transfer of patent right

Effective date of registration: 20210120

Address after: 100088 room 1202, Weishi building, 39 Xueyuan Road, Haidian District, Beijing

Patentee after: Tibet Tianyu Aviation Data Technology Co.,Ltd.

Address before: 100191 f627, new main building of Beijing University of Aeronautics and Astronautics, 37 Xueyuan Road, Haidian District, Beijing

Patentee before: BEIHANG University

TR01 Transfer of patent right
CP01 Change in the name or title of a patent holder

Address after: 100088 room 1202, Weishi building, 39 Xueyuan Road, Haidian District, Beijing

Patentee after: Tianyu Aviation Data Technology (Hefei) Co.,Ltd.

Address before: 100088 room 1202, Weishi building, 39 Xueyuan Road, Haidian District, Beijing

Patentee before: Tibet Tianyu Aviation Data Technology Co.,Ltd.

CP01 Change in the name or title of a patent holder