US9391740B2 - Method for demodulating a signal - Google Patents
Method for demodulating a signal Download PDFInfo
- Publication number
- US9391740B2 US9391740B2 US14/769,538 US201414769538A US9391740B2 US 9391740 B2 US9391740 B2 US 9391740B2 US 201414769538 A US201414769538 A US 201414769538A US 9391740 B2 US9391740 B2 US 9391740B2
- Authority
- US
- United States
- Prior art keywords
- matrix
- signal
- eigenvectors
- eigenvalues
- weighting
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 230000002087 whitening effect Effects 0.000 claims abstract description 13
- 239000011159 matrix material Substances 0.000 claims description 336
- 239000013598 vector Substances 0.000 claims description 110
- 230000009467 reduction Effects 0.000 description 18
- 230000000875 corresponding effect Effects 0.000 description 14
- 230000005540 biological transmission Effects 0.000 description 8
- 238000004364 calculation method Methods 0.000 description 8
- 230000000717 retained effect Effects 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 230000004044 response Effects 0.000 description 3
- 238000001228 spectrum Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 229940050561 matrix product Drugs 0.000 description 1
- 238000003199 nucleic acid amplification method Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2649—Demodulators
- H04L27/26524—Fast Fourier transform [FFT] or discrete Fourier transform [DFT] demodulators in combination with other circuits for demodulation
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L1/00—Arrangements for detecting or preventing errors in the information received
- H04L1/004—Arrangements for detecting or preventing errors in the information received by using forward error control
- H04L1/0045—Arrangements at the receiver end
- H04L1/0055—MAP-decoding
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity 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/0842—Weighted combining
- H04B7/0848—Joint weighting
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity 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/0842—Weighted combining
- H04B7/0848—Joint weighting
- H04B7/0851—Joint weighting using training sequences or error signal
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B7/00—Radio transmission systems, i.e. using radiation field
- H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
- H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
- H04B7/08—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the receiving station
- H04B7/0837—Diversity 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/0842—Weighted combining
- H04B7/0848—Joint weighting
- H04B7/0854—Joint weighting using error minimizing algorithms, e.g. minimum mean squared error [MMSE], "cross-correlation" or matrix inversion
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/03—Shaping networks in transmitter or receiver, e.g. adaptive shaping networks
- H04L25/03006—Arrangements for removing intersymbol interference
- H04L25/03171—Arrangements involving maximum a posteriori probability [MAP] detection
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L27/00—Modulated-carrier systems
- H04L27/26—Systems using multi-frequency codes
- H04L27/2601—Multicarrier modulation systems
- H04L27/2647—Arrangements specific to the receiver only
- H04L27/2649—Demodulators
Definitions
- This invention relates to a method for demodulating a signal.
- It also relates to a multi-antenna receiving equipment and a corresponding computer programme.
- the field of the invention is that of the encoding of digital data, intended to be transmitted in particular in the presence of transmission noise, and of the decoding of said digital data after transmission.
- the invention more particularly but not in a limited manner relates to the field of optimising the transmission of digital data, for example via a broadband radio network.
- a multi-antenna receiver comprises at least two receiving antennas in order to receive replicas of the same emitted signal, introducing as such additional diversity having for effect to improve the quality of the reception.
- the additional diversity provided is either a space diversity if the antennas are sufficiently spaced apart, or a polarisation diversity if the antennas are collocated and polarised differently, or both a portion of one and of the other of these two types of diversity.
- the invention relates more particularly to multi-carrier transmission systems, for example of the OFDM type (Orthogonal Frequency Division Multiplexing).
- This multi-carrier transmission provides a frequency diversity that is separate from the space diversity and/or polarisation diversity provided by the use of several receiving antennas.
- LTE Long Term Evolution
- TEDS TETRA Enhanced Data Service
- DAB Digital Audio Broadcasting
- DVB-T Digital Video Broadcasting-Terrestrial
- the data is in general organised into frames.
- certain symbols distributed in the temps-frequency plane are inserted among the useful information intended for the receiver.
- These symbols referred to as “driver symbols” are known to the emitter and to the receiver. They are generally used for the purposes of synchronisation and estimation of the propagation channel.
- the noise and interferences that the signal is subjected to during its propagation between the emitter and the receiver degrade the reception of the signal.
- An insufficient processing of the noise component and of the interference then generates a high demodulation and decoding error rate.
- the noise component and the interference are processed by comparing a signal resulting from the combination of the signals received by the antennas of the receiver with an estimate of the signal emitted.
- This estimate of the signal emitted can be known a priori to the receiver by using for example the driver symbols.
- the estimate of the signal emitted can also be obtained through a prior treatment of the signals received on the various antennas, for example by means of a technique known as signal demodulation with a maximum combination of the signal ratio (MRC, for “Maximum Ratio Combining”).
- MRC maximum combination of the signal ratio
- Patent application FR 11 61062 describes a method for reducing interference. Although the method for reducing interference described in this document is of higher performance than the conventional methods for reducing interference in the case of the signals received with noise from strong interference, the method may not provide performance in certain cases wherein the noise component of the signals received contains little interference.
- the purpose of this invention is to overcome at lease this problem raised by prior art.
- Such an objective is achieved with a method for demodulating a signal in a receiver comprising at least two antennas each receiving a signal transmitted through an associated radio propagation channel, the signals received corresponding to the same emitted signal comprising time and frequency distributed symbol frames in which certain symbols, referred to as driver symbols, are known to the receiver, said method comprising the steps of:
- Driver symbol means any symbol known to the receiver, in particular both the driver symbols such as designated in the standards of known OFDM systems such as LTE and TEDS, and/or synchronisation symbols and/or symbols already demodulated and decided by the receiver in a prior processing step.
- this invention makes it possible to optimally reduce the level of the interferences in a multi-antenna receiver and as such improve the Signal to Interference plus Noise Ratio (SINR).
- SINR Signal to Interference plus Noise Ratio
- the method described hereinabove makes it possible to obtain two signals c 1 and c 2 of which the noise components are orthogonal between them.
- the signal c 1 has a level of interference that is optimally reduced and the signal c 2 of which the noise is orthogonal on each symbol therefore has a maximum level of interference. No information is lost when this combination is considered.
- the method performs in all of the cases, including when the noise component of the signals received contains little interference.
- the respective coefficients ( ⁇ 1 and ⁇ 2 ) of the channel seen by the useful signal associated with a respective antenna of the receiver and used in the step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, r′ and r′′, obtained at the end of the noise-whitening step can be determined using the maximum a posteriori criterion by maximising the probability of the occurring of the channel ( ⁇ 1 and ⁇ 2 ) conditionally with the error present on the one hand in the signal c 1 and on the other hand in the signal c 2 , by taking into account the propagation channel seen in the signal c 1 and in the signal c 2 .
- the step of determining the whitening step can comprise the sub-steps of:
- weighting vectors are homogeneous with the propagation channel.
- the covariance matrices G B and G D make it possible to represent the time and frequency constraints of this channel.
- the step of determining the whitening step can further comprise a sub-step of extracting a matrix C′ B from the matrix C B and a matrix ⁇ ′ B from the matrix ⁇ B , the matrix ⁇ ′ B comprising a determined number n B ′ of eigenvalues of ⁇ B and the matrix C′ B ′ comprising the eigenvectors associated with these n B ′ eigenvalues, and a sub-step of extracting a matrix C′ D from the matrix C D and a matrix ⁇ ′ D from the matrix ⁇ D , the matrix ⁇ ′ D comprising a determined number n D ′ of eigenvalues of ⁇ D and the matrix C′ D comprising the eigenvectors associated with these n D ′ eigenvalues, the weighting vectors then being determined using matrices of eigenvectors C′ B and C′ D and of eigenvalues ⁇ ′ B and ⁇ ′ D .
- the n B ′ eigenvalues retained are the largest eigenvalues of the matrix ⁇ B and the n′ D eigenvalues retained are the largest eigenvalues of the matrix ⁇ D .
- n′ B and n′ D are less than or equal to the number of driver symbols in each frame.
- the covariance matrix G B can be the Kronecker product of a covariance matrix G B,T in the time domain and of a covariance matrix G B,F in the frequency domain and the covariance matrix G D is the Kronecker product of a covariance matrix G D,T in the time domain and of a covariance matrix G D,F in the frequency domain and the step of determining ( 28 ) can comprise the sub-steps of:
- the weighting vectors then being determined using matrices of eigenvectors C′ B,T , C′ B,F , C′ D,F and of eigenvalues ⁇ ′ B,T , ⁇ B,F , ⁇ ′ D,T , ⁇ D,F .
- the Kronecker product of two matrices A and B, the matrix A having for components (a ij ), where i is an integer between 1 and m and j is an integer between 1 and n, is the matrix product noted as A ⁇ circle around ( ⁇ ) ⁇ B and defined by the following expression:
- a ⁇ B ( a 11 ⁇ B ... ... a 1 ⁇ n ⁇ B ⁇ ... ... ⁇ a m ⁇ ⁇ 1 ⁇ B ... ... a mn ⁇ B ) .
- step of determining the signal demodulating step can comprise the sub-steps of:
- the step of determining the signal demodulating step can further comprise a sub-step of extracting a matrix C′ B from the matrix C B and a matrix ⁇ ′ B from the matrix ⁇ B , the matrix ⁇ ′ B comprising a determined number n B ′ of eigenvalues of ⁇ B and the matrix C′ B comprising the eigenvectors associated with these n B ′ eigenvalues, the weighting vectors then being determined using the matrix of eigenvectors C′ B and of eigenvalues ⁇ ′ B .
- the covariance matrix G B can be the Kronecker product of a covariance matrix G B,T in the time domain and of a covariance matrix G B,F in the frequency domain and the step of determining the signal demodulating step comprises the sub-steps of:
- the weighting vectors then being determined using matrices of eigenvectors C′ B,T , C′ B,F and of eigenvalues ⁇ ′ B,T , ⁇ ′ B,F .
- the signal emitted is a multi-carrier signal, in particular an OFDM signal.
- the invention can also apply to single-carrier systems.
- a receiving equipment comprising at least two antennas each receiving a signal transmitted through an associated radio propagation channel, the signals received corresponding to the same emitted signal comprising time and frequency distributed symbol frames in which certain symbols, referred to as driver symbols, are known to the receiver, said method comprising means for:
- the receiving equipment according to the invention can furthermore comprise means for determining the respective coefficients ( ⁇ 1 and ⁇ 2 ) of the channel seen by the useful signal associated with a respective antenna of the receiver and used in the step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, r′ and r′′, obtained at the end of the noise-whitening step, can be determined using the maximum a posteriori criterion by maximising the probability of the occurring of the channel ( ⁇ 1 and ⁇ 2 ) conditionally with the error present on the one hand in the signal c 1 and on the other hand in the signal c 2 , by taking into account the propagation channel seen in the signal c 1 and in the signal c 2 .
- a computer programme comprising instructions for the implementation of the method according to the invention when the programme is executed by at least one processor.
- FIG. 1 is a diagram of a receiving equipment according to the invention
- FIG. 2 is a diagram showing the operating principle of the method of demodulation
- FIG. 3 is a flow chart showing a first part of the operation of the method of demodulation according to the invention.
- FIG. 4 shows a useful signal in a projection base
- FIG. 5 is a flow chart showing a second part of the operation of the method of demodulation according to the invention.
- alternatives of the invention can in particular be considered that comprise only a selection of the characteristics described in what follows separated from the other characteristics described (even if this selection is isolated within a sentence comprising these other characteristics), if this selection of characteristics is sufficient to confer a technical advantage or for differentiating the invention in relation to prior art.
- This selection comprises at least one more preferably functional characteristic without structural details, or with only one portion of the structural details if this portion only is sufficient to confer a technical advantage or to differentiate the invention in relation to prior art.
- FIG. 1 shows a receiving equipment 2 of an OFDM transmission system, such as for example a base station or a mobile terminal.
- the receiver 2 comprises two antennas 4 , 6 for the reception of two signals corresponding to the same OFDM signal transmitted from an emitter through two radio propagation channels associated with the antennas 4 , 6 .
- the channels associated with the antennas 4 , 6 are supposed to respond to identical physical constraints.
- the signal OFDM emitted by the emitter is organised into time and frequency distributed frames of symbols among which certain symbols, referred to as driver symbols, are known to the receiver 2 and are stored in a memory 8 of said receiver 2 .
- Each frame comprises as such n symbols with n f sub-carriers and n t temps-symbols, n being equal to the product of n f and of n t .
- the receiver 2 comprises an input module 10 comprising amplification, baseband formatting, sampling and guard interval suppression stages.
- the receiver 2 further comprises means for the time-frequency conversion 12 of the signals received by the antennas 4 , 6 and processed by the input module 10 in order to switch them from the time domain to the frequency domain.
- These means for time-frequency conversion 12 implement a Fast Fourier Transform (FFT).
- FFT Fast Fourier Transform
- the receiver 2 further comprises an interference reduction module 14 making it possible to reduce the level of interference in a useful signal resulting from the combination of the signals received by the two antennas 4 , 6 .
- interferences can be caused, by way of examples, by the presence of scramblers emitting on the same radio channel as the emitter.
- the interference reduction module 14 is able to provide data symbols wherein the contribution of the interference is minimised.
- a demodulator 16 of the receiver 2 makes it possible to demodulate these data symbols into demodulated bits according to the same modulation technique as that used in the emitter.
- the demodulator 16 is furthermore arranged in order to determine a similarity of each demodulated bit.
- the similarity of a bit has a negative or positive soft value, in comparison with a hard value such as the binary value “1” or “0”, in order to indicate that the demodulator 16 delivers actual floating values that each have a sign which imposes a later decision, by a decoder 18 of the receiver 2 , on the state of the corresponding bit, i.e. a decision on the “hard” value “0” or “1”,
- the decoder 18 makes it possible to decode the demodulated bits supplied by the demodulator 16 according to the previously determined similarities.
- the decoder 18 implements a recoding that corresponds to the encoding used when the signal is emitted, for example a convolutive decoding that corrects the errors using the Viterbi algorithm.
- FIG. 2 shows the principle of the method for reducing interference used by the interference reduction module 14 .
- the principle used by this invention consists in applying weightings w 1 and w 2 respectively to the signals s 1 and s 2 , then in combining the two weighted signals, for example by adding them together, in order to form a signal c 1 from which is subtracted an estimate of the emitted signal d weighted by a weighting w d .
- the resulting difference ⁇ shows a residual error.
- the method of the invention uses advantageously the maximum a posteriori approach in order to calculate the most probable weightings knowing this error.
- the interference reduction module 14 weights the signals s 1 and s 2 with first respective weighting vectors (w 1 ; w 2 ) associated with a respective antenna ( 4 , 6 ) of the receiver ( 2 ).
- the first weighting vectors (w 1 ; w 2 ) define a minimum direction, noted as DIRmin, for the interference as is shown in FIG. 4 .
- the content of the first weighting vectors (w 1 ; w 2 ) remains to be determined.
- the interference reduction module 14 combines, here by adding them together, the weighted signals s 1 and s 2 in order to form a first combined signal (c 1 ).
- the interference reduction module 14 weights a reference signal d, comprising the driver symbols, with another weighting vector w d of which the content is to be determined.
- the vectors w 1 , w 2 and w d are column vectors each containing as many lines as the signal received in a frame contains symbols, i.e. n lines.
- the interference reduction module 14 determines an error ⁇ that corresponds to the difference between the first combined signal c 1 and the weighted reference signal.
- the interference reduction module 14 calculates the vector w 1 , w 2 and w d by using the MAP approach.
- This approach consists in maximising the probability of the occurring of the weighting vector w. This probability if conditional to the observation of the error ⁇ .
- this probability is equal to the probability that the error ⁇ is observed conditionally to the probability that the weighting is equal to the vector w, this probability being multiplied by the probability that the weighting vector w is carried out.
- P ( W ) f ( w / ⁇ ) ⁇ f ( ⁇ / w ) ⁇ f ( w ),
- weighting vectors w 1 , w 2 and w d are linked to the propagation channel.
- This channel is conditioned by constraints relating to its maximum time spread, due to reflections on distant obstacles, and to its maximum frequency spread, due to the speed of the receiving equipment and to the carrier frequency, i.e. the Doppler spread.
- the frequency spread of the channel due to the reflections on near obstacles is limited. This spread is between ⁇ F D and +F D , wherein F D is the maximum Doppler frequency given by the relationship
- F D v c ⁇ F p , wherein v is the speed of the receiver 2 , c is the speed of light, and F p is the carrier frequency.
- the components of the frequency spectrum of the propagation channel according to the frequency axis are therefore between these limits ⁇ F D and +F D .
- the time spread of the channel due to the reflections on the distant obstacles, is limited. This time spread depends on the frequency band used and of the environment. By way of examples, at a carrier frequency of 400 MHz, in an urban environment the time spread is about 5 ⁇ s while in a mountainous environment, this spread is about 15 ⁇ s. The components of the time response of the channel are therefore between fixed limits for given environmental conditions.
- the limits of the frequency spectrum and of the time response of the channel are known to the receiver 2 and are stored in the memory 8 .
- G The global covariance matrix
- ⁇ is a constant and the notation X H indicates that it is a conjugate and transposed matrix X.
- ⁇ is a constant and ⁇ 2 represents the variance of the noise component in the signal corresponding to the signals received on the various weighted and combined antennas.
- ⁇ 2 ⁇ 1 2 ⁇ w 1 ⁇ 2 + ⁇ 2 2 ⁇ w 2 ⁇ 2 ,
- ⁇ 1 2 is the variance of the noise component on the first antenna 4 and ⁇ 2 2 is the variance of the noise component on the second antenna 6 .
- the interference reduction module 14 attempts to minimise this logarithm L(P(w)).
- the covariance matrix G is a diagonal matrix by blocks constituted by the concatenation of covariance matrices corresponding to each one of the weighting vectors w 1 , w 2 and w d .
- the covariance matrix of w 1 is the same as that of w 2 given that the two weighting vectors w 1 and w 2 are both homogeneous with a propagation channel, on the case of two receiving antennas.
- This covariance matrix is noted as G B and shows the time and frequency constraints relative to the propagation channel.
- the weighting vector w d is homogeneous with the product, symbol by symbol, of two propagation channels, on the case of two receiving antennas.
- the corresponding covariance matrix G D shows the time and frequency constraints relative to such a product.
- the covariance matrix G can therefore be written as:
- G ( G B 0 0 0 G B 0 0 0 G D ) .
- k 2 is set to 2. This makes it possible to consider that the weightings w 1 and w 2 are each of unitary power. Then the matrix G B is the normalised covariance matrix, i.e. obtained with a unitary average power of the channel, and the matrix G D is the covariance matrix corresponding to the product, symbol by symbol, of two channels of unitary power. Consequently, ⁇ 1 shows the inverse′ of the signal-to-noise ratio observed on any of the antennas.
- the minimisation problem of the step 28 is as such an optimisation problem with constraint, which is resolved with the Lagrange multipliers according to the following relationship:
- ⁇ w designates the gradient in relation to weighting vector w and the matrix N is the identity matrix for the two signals received at the antennas 4 , 6 and the zero matrix for the reference signal.
- the matrix N can be written as:
- N ( 1 0 ... ... ... ... ... ... 0 0 ⁇ ⁇ ⁇ ⁇ ⁇ 1 ⁇ ⁇ ⁇ ⁇ ⁇ 1 ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ ⁇ 1 ⁇ ⁇ ⁇ ⁇ 0 ⁇ ⁇ ⁇ ⁇ ⁇ 0 0 ... ... ... ... ... 0 ) .
- the weighting vector w solution of the minimisation problem is as such a generalised eigenvector both of the term
- the weighting vector w minimising L(P(w)) is the generalised eigenvector corresponding to the smallest generalised eigenvalue ⁇ .
- the covariance matrices G B and G D are decomposed in the step 28 into eigenvectors and eigenvalues in the following way:
- the matrices C B and C D are matrices of eigenvectors of the matrices G B and G D respectively and the matrices ⁇ B and ⁇ D are the diagonal matrices of corresponding eigenvalues respectively.
- the inverse matrix of the covariance matrix G is equal to:
- the complexity of the calculation is further reduced by retained only certainlichnvalues of the covariance matrix G as well as the corresponding eigenvalues.
- a matrix C′ is extracted from the matrix C and a matrix ⁇ ′ is extracted from the matrix ⁇ , the matrix ⁇ ′ comprising a determined number n′ of eigenvalues of ⁇ and the matrix C′ comprising the eigenvectors associated with these n′ eigenvalues.
- the weighting vector w is then determined using matrices of eigenvectors C′ and of eigenvalues ⁇ ′.
- the covariance matrix G is a diagonal matrix par bloc, constituted of the matrices G B , G B and G D , and for which each one of the matrices G B and G D is the Kronecker product of a covariance matrix G B,T , respectively G D,T , in the time domain and of a covariance matrix G B,F , respectively G D,F , in the frequency domain and the step of determining comprises the sub-steps of:
- the weighting vector w is then determined using matrices of eigenvectors C′ B,T , C′ B,F , C′ D,T , C′ D,F and of eigenvalues ⁇ ′ B,T , ⁇ ′ B,F ⁇ ′ D,T , ⁇ ′ D,F .
- the interference reduction module 14 weights the signals received s 1 and s 2 respectively by weighting vectors w* 2 and ⁇ w* 1 .
- the vector (w 1 , w 2 ) referred to as the first weighting vector, formed by the first weighting vectors (w 1 , w 2 ) is orthogonal (i.e. the hermitian product is zero) on each symbol to the vector (w* 2 , ⁇ w* 1 ), referred to as the second weighting vector, formed by the second weighting vectors (w* 2 , ⁇ w* 1 ).
- the second weighting vectors (w* 2 , ⁇ w* 1 ) define a maximum direction, noted as DIRmax, for the interference as is shown in FIG. 4 .
- the notation X* indicates that it is a conjugate matrix X.
- the interference reduction module 14 combines, here by adding them together, the weighted signals c 1 and c 2 by the second weighting vectors (w* 2 , ⁇ w* 1 ) in order to form a combined signal c 2 .
- the vectors ⁇ 1 , ⁇ 2 are calculated by using the MAP approach.
- the interference reduction module 14 weights the signals received s 1 and s 2 respectively by weighting vectors w 1 and w 2 determined previously.
- the interference reduction module 14 combines, here by adding them together, the weighted signals s 1 and s 2 in order to form a combined signal c 1 .
- the interference reduction module 14 weights a reference signal d, comprising the driver symbols, with another weighting vector w d .
- w d can be obtained by combining a weighting of the first weighting vector (w 1 , w 2 ) determined during the step of determining ( 28 ) of the step of whitening E 1 by the channel vector ( ⁇ 1 , ⁇ 2 ).
- w d w 1 ⁇ 1 +w 2 ⁇ 2 .
- the vectors ⁇ 1 , ⁇ 2 are column vectors each containing as many lines as the signal received in a frame contains symbols, i.e. n lines.
- the interference reduction module 14 calculates the vector ⁇ 1 , ⁇ 2 by using the MAP approach.
- This approach consists in maximising the probability of the occurring of the weighting vector ⁇ . This probability is conditional to the observation of the error ⁇ .
- this probability is equal to the probability that the error ⁇ is observed conditionally to the probability that the weighting is equal to the vector ⁇ , this probability being multiplied by the probability that the weighting vector ⁇ is carried out.
- P ( ⁇ ) f ( ⁇ / ⁇ ) ⁇ f ( ⁇ / ⁇ ) ⁇ f ( ⁇ )
- the limits of the frequency spectrum and of the time response of the channel are known to the receiver 2 and are stored in the memory 8 .
- G 2 ( G B 0 0 G B )
- G B C B , ⁇ B , C B H
- G B - 1 C B , ⁇ B - 1 , C B H
- ⁇ ′ is a constant
- ⁇ ′ is a constant and ⁇ ′ 2 represents the variance of the noise component in the signal corresponding to the signals received on the various weighted and combined antennas.
- the interference reduction module 14 attempts to minimise this logarithm L(P( ⁇ )), which reverts to minimising:
- the complexity of the calculation is further reduced by retained only certainlichnvalues of the covariance matrix G 2 as well as the corresponding eigenvalues.
- a matrix C′ is extracted from the matrix C 2 and a matrix ⁇ ′ is extracted from the matrix ⁇ 2 , the matrix ⁇ ′ comprising a determined number n′ of eigenvalues of ⁇ 2 and the matrix C′ comprising the eigenvectors associated with these n′ eigenvalues.
- the weighting vector ⁇ is then determined using matrices of eigenvectors C′ and of eigenvalues ⁇ ′.
- the covariance matrix G 2 is a diagonal matrix par bloc, constituted of the matrices G B , G B , and for which each one of the matrices G B is the Kronecker product of a covariance matrix G B,T , in the time domain and of a covariance matrix G B,F , in the frequency domain and the step of determining comprises the sub-steps of:
- the weighting vector ⁇ is then determined using matrices of eigenvectors C′ B,T , C′ B,F , and of eigenvalues ⁇ B,T and ⁇ B,F .
- the coefficients ⁇ 1 , ⁇ 2 are then used in the step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals (r′, r′′).
Landscapes
- Engineering & Computer Science (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Discrete Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Probability & Statistics with Applications (AREA)
- Power Engineering (AREA)
- Radio Transmission System (AREA)
- Noise Elimination (AREA)
Abstract
Description
-
- noise whitening in order to form two combined signals the noise components of which are separate;
- normalising the noise components of the combined signals in order to form two signals the noise components of which are separate and have equal average standards;
- performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, the noise components of which are separate and which have equal average standards
-
- weighting of each one of the signals received with respective first weighting vectors associated with a respective antenna of the receiver, the signal associated with the first antenna being weighted by a vector w1 and the signal associated with the second antenna being weighted by a vector w2,
- combining weighted received signals in order to form a first combined signal (c1),
- weighting a reference signal comprising said driver symbols with another weighting vector (Wd),
- comparing the first combined signal and the weighted reference signal in order to form an error (ε), and
- determining weighting vectors (w1, w2 and wd) using the maximum a posteriori criterion by maximising the probability of the occurring of said weighting vectors conditionally with the error obtained,
- weighting of each signal received with second weighting vectors, the signal received on the first antenna being weighted by the conjugate of the complex vector w2 and the signal received on the first antenna being weighted by the opposite of the conjugate of the complex vector w1,
- combining received signals weighted by the second weighting vectors in order to form a second combined signal (c2).
-
- calculating a covariance matrix GB of the propagation channel;
- calculating a covariance matrix GD of the product, symbol by symbol, of two propagation channels; and
- determining weighting vectors using covariance matrices GB and GD.
-
- decomposing the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,T H and GB,F=CB,FΛB,FCB,F H, wherein:
- CB,T is a matrix of eigenvectors of the matrix GB,T;
- ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
- CB,F is a matrix of eigenvectors of the matrix GB,F and
- ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
- extracting of a matrix from the matrix C′B,T and a matrix ΛB,T from the matrix ΛB,T the matrix ΛB,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues;
- extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
- decomposing of the covariance matrices GD,T and GD,F into eigenvectors according to the relationships GD,T=CD,TΛD,T H and GD,F=CD,FΛD,F H, wherein:
- CD,T is a matrix of eigenvectors of the matrix GD,T;
- ΛD,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,T;
- CD,F is a matrix of eigenvectors of the matrix GD,F; and
- ΛD,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,F;
- extracting of a matrix C′D,T from the matrix CD,T and a matrix Λ′D,T from the matrix ΛD,T, the matrix Λ′D,T comprising a determined number nDt′ of the largest eigenvalues of ΛD,T and the matrix C′D,T comprising the eigenvectors associated with these nDt′ largest eigenvalues; and
- extracting of a matrix C′D,F from the matrix CD,F and a matrix ΛD,F from the matrix ΛD,F, the matrix Λ′D,F comprising a determined number nDf′ of the largest eigenvalues of ΛD,F and the matrix C′D,F comprising the eigenvectors associated with these nDf′ largest eigenvalues,
-
- calculating a covariance matrix GB of the propagation channel;
- determining weighting vectors using the covariance matrix GB.
-
- decomposing of the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,T H and GB,F=CB,FΛB,FCB,F H, wherein:
- CB,T is a matrix of eigenvectors of the matrix GB,T;
- ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
- CB,F is a matrix of eigenvectors of the matrix GB,F; and
- ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
- extracting of a matrix C′B,T from the matrix CB,T and a matrix Λ′B,T from the matrix ΛB,T, the matrix Λ′B,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues;
- extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
-
- noise whitening in order to form two combined signals the noise components of which are separate;
- normalising the noise components of the combined signals in order to form two signals the noise components of which are separate and have equal average standards;
- step of performing signal demodulation with the maximum combination of the signal-to-noise ratio on the two signals, the noise components of which are separate and which have equal average standards;
-
- weighting of each one of the signals received with respective first weighting vectors associated with a respective antenna of the receiver, the signal associated with the first antenna being weighted by a vector w1 and the signal associated with the second antenna being weighted by a vector w2,
- combining weighted received signals in order to form a first combined signal (c1).
- weighting a reference signal comprising said driver symbols with another weighting vector (Wd),
- comparing the first combined signal and the weighted reference signal in order to form an error (ε), and
- determining weighting vectors (w1, W2 and Wd) using the maximum a posteriori criterion by maximising the probability of the occurring of said weighting vectors conditionally with the error obtained,
- weighting of each signal received with second weighting vectors, the signal received on the first antenna being weighted by the conjugate of the complex vector w2 and the signal received on the first antenna being weighted by the opposite of the conjugate of the complex vector w1, combining received signals weighted by the second weighting vectors in order to form a second combined signal (c2).
P(W)=f(w/ε)∝f(ε/w)·f(w),
wherein v is the speed of the
σ2=σ1 2 ∥w 1∥2+σ2 2 ∥w 2∥2,
σ2=σ1 2(∥w 1∥2 +∥w 2∥2)(.
k 2 =∥w 1∥2 +∥w 2∥2 =cst 2.
and of the term NHN and μ is the associated generalised eigenvalue.
-
- decomposing of the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,T H and GB,F=CB,FΛB,FCB,F H, wherein:
- CB,T is a matrix of eigenvectors of the matrix GB,T;
- ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
- CB,F is a matrix of eigenvectors of the matrix GB,F; and
- ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,F;
- extracting of a matrix C′B,T from the matrix CB,T and a matrix Λ′VB,T from the matrix ΛB,T, the matrix Λ′B,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues; and
- extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
- decomposing of covariance matrices GD,T and GD,F into eigenvectors according to the relationships GD,T=CD,TΛD,TCD,T H and GD,F=CD,FΛD,FCD,F H, wherein:
- CD,T is a matrix of eigenvectors of the matrix GD,T;
- ΛD,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,T;
- CD,F is a matrix of eigenvectors of the matrix GD,F; and
- ΛD,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GD,F;
- extracting of a matrix C′D,T from the matrix CD,T and a matrix Λ′D,T from the matrix ΛD,T, the matrix Λ′D,T comprising a determined number nDt′ of the largest eigenvalues of ΛD,T and the matrix C′D,T comprising the eigenvectors associated with these nDt′ largest eigenvalues; and
- extracting of a matrix C′D,F from the matrix CD,F and a matrix Λ′D,F from the matrix ΛD,F, the matrix Λ′D,F comprising a determined number nDf′ of the largest eigenvalues of ΛD,F and the matrix C′D,F comprising the eigenvectors associated with these nDf′ largest eigenvalues.
w=C·y.
P(α)=f(α/ε)∝f(ε/α)·f(α),
-
- decomposing of the covariance matrices GB,T and GB,F into eigenvectors according to the relationships GB,T=CB,TΛB,TCB,T H and GD,F=CB,FΛB,FCB,F H, wherein:
- CB,T, is a matrix of eigenvectors of the matrix GB,T;
- ΛB,T is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
- CB,F is a matrix of eigenvectors of the matrix GB,F; and
- ΛB,F is a diagonal matrix comprising eigenvalues associated with the eigenvectors of the matrix GB,T;
- extracting of a matrix C′B,T from the matrix CB,T and a matrix Λ′B,T from the matrix ΛB,T, the matrix Λ′B,T comprising a determined number nBt′ of the largest eigenvalues of ΛB,T and the matrix C′B,T comprising the eigenvectors associated with these nBt′ largest eigenvalues; and
- extracting of a matrix C′B,F from the matrix CB,F and a matrix Λ′B,F from the matrix ΛB,F, the matrix Λ′B,F comprising a determined number nBf′ of the largest eigenvalues of ΛB,F and the matrix C′B,F comprising the eigenvectors associated with these nBf′ largest eigenvalues,
Claims (13)
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
FR1351578A FR3002709B1 (en) | 2013-02-22 | 2013-02-22 | METHOD FOR DEMODULATING A SIGNAL |
FR1351578 | 2013-02-22 | ||
PCT/FR2014/000041 WO2014128369A1 (en) | 2013-02-22 | 2014-02-20 | Method for demodulating a signal |
Publications (2)
Publication Number | Publication Date |
---|---|
US20160013888A1 US20160013888A1 (en) | 2016-01-14 |
US9391740B2 true US9391740B2 (en) | 2016-07-12 |
Family
ID=49111295
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
US14/769,538 Active US9391740B2 (en) | 2013-02-22 | 2014-02-20 | Method for demodulating a signal |
Country Status (5)
Country | Link |
---|---|
US (1) | US9391740B2 (en) |
EP (1) | EP2959610B1 (en) |
ES (1) | ES2759857T3 (en) |
FR (1) | FR3002709B1 (en) |
WO (1) | WO2014128369A1 (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016127440A (en) * | 2015-01-05 | 2016-07-11 | 富士通株式会社 | Wireless communication system, wireless communication method, transmitter, and control method for transmitter |
CN109561041B (en) * | 2018-12-13 | 2021-08-13 | 西南交通大学 | Communication sequence construction method, system, equipment and computer storage medium |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6437737B1 (en) * | 1999-04-16 | 2002-08-20 | Science And Applied Technology, Inc. | Antenna data compression using multi-dipole antenna |
US20030142755A1 (en) | 2002-01-31 | 2003-07-31 | Richard Chi | Closed loop transmit diversity antenna verification using trellis decoding |
EP1530333A1 (en) | 2003-11-05 | 2005-05-11 | Siemens Mobile Communications S.p.A. | Method for channel estimation in a MIMO OFDM system |
US20080108314A1 (en) * | 2006-11-08 | 2008-05-08 | Sony Corporation | Wireless communication system, and apparatus and method for wireless communication |
FR2924884A1 (en) | 2007-12-11 | 2009-06-12 | Eads Secure Networks Soc Par A | REDUCTION OF INTERFERENCES IN AN ORTHOGONAL FREQUENCY-DISTRIBUTED SIGNAL |
US20090274235A1 (en) * | 2006-11-15 | 2009-11-05 | Wook Bong Lee | Data transmission method using dirty paper coding in mimo system |
FR2983667A1 (en) | 2011-12-01 | 2013-06-07 | Cassidian Sas | INTERFERENCE REDUCTION METHOD |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
FR1161062A (en) | 1956-06-09 | 1958-08-20 | Process for suppressing spontaneous combustion in agglomerated blocks of expanded cork |
-
2013
- 2013-02-22 FR FR1351578A patent/FR3002709B1/en not_active Expired - Fee Related
-
2014
- 2014-02-20 US US14/769,538 patent/US9391740B2/en active Active
- 2014-02-20 WO PCT/FR2014/000041 patent/WO2014128369A1/en active Application Filing
- 2014-02-20 ES ES14712699T patent/ES2759857T3/en active Active
- 2014-02-20 EP EP14712699.9A patent/EP2959610B1/en active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6437737B1 (en) * | 1999-04-16 | 2002-08-20 | Science And Applied Technology, Inc. | Antenna data compression using multi-dipole antenna |
US20030142755A1 (en) | 2002-01-31 | 2003-07-31 | Richard Chi | Closed loop transmit diversity antenna verification using trellis decoding |
EP1530333A1 (en) | 2003-11-05 | 2005-05-11 | Siemens Mobile Communications S.p.A. | Method for channel estimation in a MIMO OFDM system |
US20080108314A1 (en) * | 2006-11-08 | 2008-05-08 | Sony Corporation | Wireless communication system, and apparatus and method for wireless communication |
US20090274235A1 (en) * | 2006-11-15 | 2009-11-05 | Wook Bong Lee | Data transmission method using dirty paper coding in mimo system |
FR2924884A1 (en) | 2007-12-11 | 2009-06-12 | Eads Secure Networks Soc Par A | REDUCTION OF INTERFERENCES IN AN ORTHOGONAL FREQUENCY-DISTRIBUTED SIGNAL |
US20100322339A1 (en) * | 2007-12-11 | 2010-12-23 | Eads Secure Networks | Reduction of interference in an ofdm signal with weighting vectors evolving in vector subspaces |
FR2983667A1 (en) | 2011-12-01 | 2013-06-07 | Cassidian Sas | INTERFERENCE REDUCTION METHOD |
Non-Patent Citations (1)
Title |
---|
International Search Report as issued in International Patent Application No. PCT/FR2014/000041, dated May 9, 2014. |
Also Published As
Publication number | Publication date |
---|---|
WO2014128369A1 (en) | 2014-08-28 |
EP2959610A1 (en) | 2015-12-30 |
FR3002709B1 (en) | 2015-03-06 |
EP2959610B1 (en) | 2019-10-09 |
FR3002709A1 (en) | 2014-08-29 |
US20160013888A1 (en) | 2016-01-14 |
ES2759857T3 (en) | 2020-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8401131B2 (en) | Reduction of interference in an OFDM signal with weighting vectors evolving in vector subspaces | |
US8804884B2 (en) | Method and apparatus for suppressing co-channel interference | |
US8150326B2 (en) | Signal processing for multi-sectored wireless communications system and method thereof | |
US8428158B2 (en) | Method for channel estimation using recursive filtering and multicarrier receiver with interference-aware demodulation | |
US7539259B2 (en) | OFDM channel estimation and tracking for multiple transmit antennas | |
US8121210B2 (en) | Apparatus and method for signal detection in multiple input multiple output wireless communication system | |
EP1802000A1 (en) | Apparatus and method for cancelling interference from neighbor cells in broadband communication system | |
EP2533449A1 (en) | Transmitter apparatus, receiver apparatus, wireless communication system, transmission control method, reception control method, and processor | |
US8824605B2 (en) | Receiving device, receiving method, receiving program, and processor | |
EP1892908A1 (en) | Interference cancellation receiver and method | |
WO2014139550A1 (en) | Iterative interference cancellation method | |
US9391740B2 (en) | Method for demodulating a signal | |
CN109560850B (en) | MRC soft detection method, device, equipment and computer readable storage medium | |
KR100934170B1 (en) | Channel Estimation Apparatus and Method in Multi-antenna Wireless Communication System | |
US8982979B2 (en) | Methods and systems for reducing complexity of MIMO decoder | |
US9467218B2 (en) | Interference reduction method | |
US9166829B2 (en) | Method for estimating a radio channel | |
US9059751B2 (en) | Systems and methods for multiplicative superposition for MIMO broadcast channels | |
US9755769B2 (en) | Method for estimating a radioelectric propagation channel | |
KR101222130B1 (en) | Multiple input multiple output radio communication system with pre-equalizer and its mehtod | |
EP1932304A2 (en) | Ofdm reception in a multi-antenna mode | |
US9036461B2 (en) | LTE—advanced transmit diversity decoders | |
KR20210058654A (en) | Communication node performing repeated reception in moving network and operating method thereof | |
EP3139528B1 (en) | Communication device and method for signal determination in radio communication |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AS | Assignment |
Owner name: AIRBUS DS SAS, FRANCE Free format text: ASSIGNMENT OF ASSIGNORS INTEREST;ASSIGNORS:MEGE, PHILIPPE;MARTINOD, LAURENT;SIGNING DATES FROM 20150824 TO 20150826;REEL/FRAME:036544/0522 |
|
AS | Assignment |
Owner name: AIRBUS DS SAS, FRANCE Free format text: CHANGE OF NAME;ASSIGNOR:CASSIDIAN SAS;REEL/FRAME:037018/0596 Effective date: 20140717 |
|
STCF | Information on status: patent grant |
Free format text: PATENTED CASE |
|
MAFP | Maintenance fee payment |
Free format text: PAYMENT OF MAINTENANCE FEE, 4TH YEAR, LARGE ENTITY (ORIGINAL EVENT CODE: M1551); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY Year of fee payment: 4 |
|
FEPP | Fee payment procedure |
Free format text: MAINTENANCE FEE REMINDER MAILED (ORIGINAL EVENT CODE: REM.); ENTITY STATUS OF PATENT OWNER: LARGE ENTITY |