WO2012079348A1 - Methods for estimating broadband co-channel interference and noise and suppressing interference, and systems thereof - Google Patents

Methods for estimating broadband co-channel interference and noise and suppressing interference, and systems thereof Download PDF

Info

Publication number
WO2012079348A1
WO2012079348A1 PCT/CN2011/075517 CN2011075517W WO2012079348A1 WO 2012079348 A1 WO2012079348 A1 WO 2012079348A1 CN 2011075517 W CN2011075517 W CN 2011075517W WO 2012079348 A1 WO2012079348 A1 WO 2012079348A1
Authority
WO
WIPO (PCT)
Prior art keywords
data
data stream
subcarrier
interference
pilot
Prior art date
Application number
PCT/CN2011/075517
Other languages
French (fr)
Chinese (zh)
Inventor
朱登魁
宁迪浩
肖华华
鲁照华
Original Assignee
中兴通讯股份有限公司
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 中兴通讯股份有限公司 filed Critical 中兴通讯股份有限公司
Publication of WO2012079348A1 publication Critical patent/WO2012079348A1/en

Links

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/03006Arrangements for removing intersymbol interference
    • H04L25/03159Arrangements for removing intersymbol interference operating in the frequency domain

Definitions

  • the present invention relates to the field of communications, and in particular, to a method and a corresponding system for wideband co-channel interference noise estimation and interference suppression.
  • Multi-antenna technology is a major breakthrough in smart antenna technology in the field of wireless mobile communications. This technology can multiply the capacity and spectrum utilization of communication systems without increasing bandwidth. It can also use multipath to mitigate multipath fading. It can effectively eliminate channel interference, improve channel reliability, and reduce bit error rate. It is a key technology for a new generation of mobile communication systems. Multi-antenna technology has been widely used in many wireless broadband systems such as Long Term Evolution (LTE) and World Interoperability for Microwave Access (WiMAX). For wireless communications in which the network is arranged in a cellular structure, co-channel interference between adjacent cells is one of the most important factors leading to a degradation in communication quality, as shown in FIG.
  • LTE Long Term Evolution
  • WiMAX World Interoperability for Microwave Access
  • the receiving end must accurately estimate the channel coefficient, the interference channel coefficient or the interference feature of the desired data to perform more accurate data detection.
  • the data pilot subcarriers between adjacent cells coincide at the time-frequency position, it brings great difficulty to the interference estimation. Since the coincidence of the interference pilots leads to a degradation of the channel estimation quality, that is, the channel estimation itself carries the interference information, making the estimation of the interference noise characteristics very difficult or very inaccurate.
  • the performance of common interference suppression receiving algorithms such as Minimum Mean-Square Error (MMSE) or Interference Rejection Combination (IRC) algorithm, will be greatly reduced.
  • MMSE Minimum Mean-Square Error
  • IRC Interference Rejection Combination
  • the present invention provides a method for wideband co-channel interference noise estimation, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference
  • OFDM Orthogonal Frequency Division Multiplexing
  • OFDMA Orthogonal Frequency Division Multiple Access
  • the method includes: receiving, on each pilot subcarrier corresponding to the data stream, the reception on each pilot subcarrier corresponding to the data stream a weighted average of the product of the signal and the conjugate transpose of the received signal as an interference noise covariance matrix estimate for the data subcarrier position; and for each data subcarrier corresponding to the data stream, the data subcarrier will be The interference noise covariance matrix estimated value
  • the weighted average of the product of the received signal on each pilot subcarrier corresponding to the data stream and the conjugate transposed on the received signal for each data subcarrier corresponding to the data stream In the step of estimating the interference noise covariance matrix of the data subcarrier position, the equation (a) is used to calculate the interference noise covariance matrix estimate of the data subcarrier position:
  • the number of carriers; 1, ⁇ ⁇ ⁇ , /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region;
  • is to calculate the _th data corresponding to the data stream in the interference suppression region
  • the interference noise covariance matrix estimate for the subcarrier position is given the weight of ( ( i ) ) ( y p (i )f ,
  • the method further includes: The interference suppression area is divided into one or more interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least one pilot subcarrier and one data subcarrier therein;
  • the same weight is given to the product of the received signal on each pilot subcarrier in the same interference noise estimation unit and its conjugate transposed.
  • the method further includes:
  • the interference suppression region is further divided into an interference noise estimation unit, and each interference noise estimation unit is a time domain two.
  • a set of indices ⁇ of the pilot subcarriers corresponding to the data stream included in the first interference noise estimation unit 1, ⁇ ⁇ ⁇ , /, / is the pilot corresponding to the data stream in the interference suppression region
  • is the received signal of the receiving end on the first pilot subcarrier corresponding to the data stream; (a wide () conjugate transposed matrix;
  • ⁇ ⁇ ⁇ is calculated as ic M — D , it is assigned to each pilot subcarrier corresponding to the data stream in the third interference noise estimation unit (the weight of the product of (,
  • the data interference result obtained by diagonally loading the interference noise covariance matrix estimation value of the data subcarrier position is used as the data subcarrier position for each data subcarrier corresponding to the data stream.
  • the number of data subcarriers; ⁇ 0; is the estimated value of the interference noise covariance matrix of the _; data subcarrier positions corresponding to the data stream in the interference suppression region; ⁇ >0 ⁇ represents the diagonal matrix of N xN 3 ⁇ 4 And 3 ⁇ 4 denote the number of receiving antennas at the receiving end.
  • the interference noise covariance matrix of the first data subcarrier position is:
  • each pilot subcarrier position corresponding to the data stream is calculated in the following manner. Estimated value of the channel coefficient: Multiply the received signal on the pilot subcarrier with the conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier to obtain an estimated channel coefficient of the pilot subcarrier position.
  • the present invention also provides a system for wideband co-channel interference noise estimation, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression region.
  • OFDM Orthogonal Frequency Division Multiplexing
  • OFDMA Orthogonal Frequency Division Multiple Access
  • the interference suppression area being a time-frequency two-dimensional resource block in the received data bearer area
  • the system comprising: a first device, configured to: each corresponding to the data stream a data subcarrier, a weighted average of a product of a received signal on each pilot subcarrier corresponding to the data stream and a conjugate transposed of the received signal, as an interference noise covariance matrix estimation value of the data subcarrier position; And a second device, configured to: diagonally load the estimated value of the interference noise covariance matrix of the data subcarrier position obtained by the first device for each data subcarrier corresponding to the data stream As a result, the interference noise covariance matrix is the position of the data subcarrier.
  • the first device is configured to calculate an interference noise covariance matrix estimate of the data subcarrier position using equation (a):
  • the number of carriers; 1, ⁇ ⁇ ⁇ , /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region;
  • is to calculate the first data subcarrier corresponding to the data stream in the interference suppression region imparting ((I)) weight (y p (i) f the noise covariance matrix estimate interference position,
  • a third device is included, where: the third device is configured to: divide the interference suppression region into one or more interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least one a pilot subcarrier and a data subcarrier; the first device calculates the interference noise covariance matrix estimation value of the data subcarrier position according to the equation (a), and is the reception on each pilot subcarrier in the same interference noise estimation unit Signal and its conjugate Set the product, giving the same weight.
  • the system further includes a fourth device, where: the fourth device is configured to: divide the interference suppression region into interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least a pilot subcarrier and a data subcarrier are positive integers; for each data subcarrier corresponding to the data stream, the first device is set to be on each pilot subcarrier corresponding to the data stream according to the following formula A weighted average of the product of the received signal and the conjugate transpose of the received signal as the estimated value of the interference noise covariance matrix for the data subcarrier position: among them,
  • a set of index I of pilot subcarriers corresponding to the data stream included in the first interference noise estimation unit 1, . . . , /, / is the pilot subcarrier corresponding to the data stream in the interference suppression region Number of
  • ⁇ ;) is the received signal of the receiving end on the first pilot subcarrier corresponding to the data stream; y p (if is the conjugate transposed matrix of (); fi ml is calculated - D , the first Each channel corresponding to the data stream in the interference noise estimation unit
  • the second device for each data subcarrier corresponding to the data stream, is configured to diagonally load the estimated value of the interference noise covariance matrix of the data subcarrier position according to the following formula: As a result, the interference noise covariance matrix as the data subcarrier position: NI-D
  • R W _ D / is the interference noise covariance matrix of the _th data subcarrier position corresponding to the data stream in the interference suppression region
  • 1, ⁇
  • J is the interference suppression region
  • the number of data subcarriers corresponding to the data stream; ⁇ 0; is the estimated value of the interference noise covariance matrix of the _; data subcarrier positions corresponding to the data stream in the interference suppression region; ⁇ >0 ⁇ represents N xN 3 ⁇ 4
  • the diagonal matrix, 3 ⁇ 4, represents the number of receiving antennas at the receiving end.
  • the interference noise covariance matrix of the first data subcarrier position calculated by the second device is:
  • the unit matrix further includes: a fifth device, wherein: the fifth device is configured to obtain, at the following manner, a channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream at the transmitting end and output the The first device: multiplies the received signal on the pilot subcarrier by a conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier to obtain an estimated channel coefficient of the pilot subcarrier position.
  • the above-mentioned wideband co-channel interference estimation method can obtain more accurate interference noise characteristics, which is beneficial to improve the performance of interference suppression and the accuracy of data detection.
  • Another object of the present invention is to provide a method for suppressing broadband co-channel interference and a corresponding system to solve the problem of poor interference suppression performance when adjacent cells have co-channel interference.
  • the present invention provides a method for wideband co-channel interference suppression, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression manner.
  • OFDM Orthogonal Frequency Division Multiplexing
  • OFDMA Orthogonal Frequency Division Multiple Access
  • the method includes: obtaining, according to the interference noise estimation method, the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream, and the interference noise covariance matrix of each data subcarrier position; corresponding to the data stream a weighted average of the channel coefficient estimates of the respective pilot subcarrier positions corresponding to the data stream, as the channel coefficient estimation value of the data subcarrier position; and each data substring corresponding to the data stream Transmitting, according to the received signal on the data subcarrier, and the channel coefficient estimation value and the interference noise covariance matrix of the data subcarrier position, calculating a data signal estimate on the data subcarrier; wherein the interference suppression region is receiving A time-frequency two-dimensional resource block in the data bearer area.
  • F is the number of bits included in the data signal on the Jth data carrier corresponding to the data stream in the interference suppression area, and indicates the channel gain corresponding to the _/th data carrier corresponding to the data stream in the interference suppression area. 4 ( ) (R W — (7) or signal to interference and noise ratio
  • each channel estimation unit is a time domain two-dimensional a resource block and including at least one pilot subcarrier and one data subcarrier, where f is a positive integer; and for each data subcarrier corresponding to the data stream, a channel of each pilot subcarrier location corresponding to the data stream.
  • O kl is calculated, giving the first / channel estimation unit for each pilot sub-carrier position ⁇ ⁇ () of
  • 1 , 0 ⁇ ⁇ 1 ,
  • indicates the number of pilot subcarriers included, and is in the right
  • the present invention also provides a system for wideband co-channel interference suppression, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression region.
  • OFDM Orthogonal Frequency Division Multiplexing
  • OFDMA Orthogonal Frequency Division Multiple Access
  • Interference suppression is performed on a data stream carried in the received data bearer region
  • the system includes: a first subsystem, configured to: In the same manner as the system for interference noise estimation, the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream and the interference noise covariance matrix of each data subcarrier position are obtained; the second subsystem is set to: Each data subcarrier corresponding to the data stream, a weighted average of channel coefficient estimates of respective pilot subcarrier positions corresponding to the data stream, as a channel coefficient estimate of the data subcarrier position; and a third subsystem, Set to: each data subcarrier corresponding to the data stream, according to the received signal on the data subcarrier, and the location of the data subcarrier Channel coefficient estimates and the covariance matrix of interference noise, the data signal is calculated on the data subcarrier estimated.
  • the system further includes a fourth subsystem, wherein: the fourth subsystem is configured to: calculate, for each data subcarrier corresponding to the data stream, a data signal estimate on the data carrier calculated by the third subsystem After that, the estimated data signal is demodulated, and the soft information corresponding to each bit of the signal obtained after demodulation is adjusted, and the soft information corresponding to each bit after adjustment is respectively
  • ' l, ⁇ ⁇ ⁇ , J , J are the number of data subcarriers corresponding to the data stream in the interference suppression region;
  • F is the number of bits included in the data signal on the Jth data carrier corresponding to the data stream in the interference suppression area, and indicates the _/ ⁇ th data carrier pair corresponding to the data stream in the interference suppression area.
  • M _ P W is the interference noise covariance matrix of the first pilot subcarrier position, where 1 is the inverse matrix of ⁇ ; () is the number corresponding to the data stream in the interference suppression region of the receiving end, The received signal on the pilot
  • the system further includes a fifth subsystem and a sixth subsystem, wherein: the fifth subsystem is configured to: divide the interference suppression region into K channel estimation units, each channel estimation unit is a time domain two-dimensional a resource block and including at least one pilot subcarrier and one data subcarrier, where f is a positive integer; the sixth subsystem is set to: for each data subcarrier corresponding to the data stream, using the following formula to the data The weighted average of the channel coefficient estimates for the respective pilot subcarrier locations corresponding to the stream, as the channel coefficient estimates for the data subcarrier locations:
  • the number of pilot subcarriers; the channel of the first pilot subcarrier position corresponding to the data stream in the interference suppression region Coefficient estimate; a kl is the ⁇ ⁇ () assigned to the position of each pilot subcarrier in the channel estimation unit when calculating
  • 1 , 0 ⁇ ⁇ 1 ,
  • indicates the number of pilot subcarriers included, and is in the right
  • the above-mentioned method for suppressing broadband co-channel interference is based on relatively accurate interference noise characteristics, which can improve the performance of interference suppression and the accuracy of data detection.
  • FIG. 1 is a schematic diagram of a neighboring multi-cell in the prior art
  • FIG. 2 is a flowchart of a method for estimating and suppressing interference of wideband co-channel interference according to an embodiment of the present invention
  • FIG. 3 to FIG. A schematic diagram of three ways of dividing, the bold coil represents the pilot subcarrier, the thin coil represents the data subcarrier, and the following FIG. 6 and FIG. 7 are the same
  • FIG. 6 is a division of the interference suppression region pattern 2
  • FIG. 7 is a schematic diagram of a manner of dividing the interference suppression region pattern three
  • FIG. 8 is a schematic diagram of a manner of dividing the interference suppression region pattern four, where the thick coil indicates the data stream 1 corresponding The pilot subcarrier, the dotted coil represents the pilot subcarrier corresponding to the data stream 2, and the thin coil represents the data subcarrier, which is the same as FIG. 9;
  • FIG. 9 is a manner of dividing the interference suppression region pattern five. Schematic diagram.
  • the transmitting end in the text may be a control device such as a base station and a relay station, or may be a terminal device such as a mobile phone, a notebook computer, or a handheld computer.
  • the receiving end is configured to receive the data signal of the transmitting end, and the receiving end may be a terminal device such as a mobile phone, a notebook computer, and a handheld computer, or may be a control device such as a base station and a relay station.
  • the receiving end divides the received data bearer area into one or more interference suppression areas, and each interference suppression area is a time-frequency two-dimensional resource block in a frame or a field structure, that is, each interference suppression area includes multiple times in time.
  • a continuous OFDM/OFDMA symbol comprising a plurality of consecutive subcarriers in the frequency domain.
  • the receiving data bearer area may include a time-frequency two-dimensional resource block, and may also include a plurality of separate time-frequency two-dimensional resource blocks.
  • each of the independent time-frequency two-dimensional resource blocks is used as an interference. Suppress area.
  • the relatively independent time-frequency two maintenance resources in the received data bearer area may be further divided into multiple interference suppression areas.
  • the interference suppression area may carry one or more data streams, and each data stream corresponds to one or more data subcarriers and pilot subcarriers, and different pilot streams corresponding to different data streams are different.
  • Step 20 for each data subcarrier corresponding to the data stream, the location of the data subcarrier will be The result obtained by diagonally loading the estimated value of the interference noise covariance matrix as the interference noise covariance matrix of the data subcarrier position;
  • ⁇ 0, ⁇ 0, ⁇ represents N xN 3 ⁇ 4 diagonal matrix, i.e., in addition to the non-zero diagonal elements, the elements of the other positions are 0, N 3 ⁇ 4 denotes the number of receive antennas.
  • N 3 ⁇ 4 denotes the number of receive antennas.
  • R w — D O) R M _ D O) + 7 Spotify ( 3 )
  • the estimation of the broadband co-channel interference noise in the interference suppression region is performed. After the interference suppression regions in the data bearer region are calculated according to the above method, the broadband co-channel interference noise estimation for the data bearer region is completed.
  • the operation of this step is a conventional operation.
  • the data signal corresponding to the data subcarrier DsC() corresponding to the data stream in the interference suppression region is calculated by the following method: When ⁇ (represented as column direction) Time,
  • 4( ⁇ ) is the estimated channel coefficient corresponding to the data subcarrier DsC(/), which is a conjugated transposition of 4( ⁇ )
  • ′′ indicates that the element of 4( ⁇ ) is ⁇
  • (R w —.(_ /)) - 1 is the inverse matrix of R w — )
  • (_/ ⁇ ) is the received signal on DsC().
  • (_/) is represented as a column vector, and if (7 ⁇ ) is expressed as a row vector, the above formula needs to be adaptively changed, and will not be described again.
  • the corresponding data calculation can be obtained according to the above steps.
  • the specific weight selection can be different.
  • the pilot subcarriers and the data subcarriers corresponding to the data stream in the foregoing steps all refer to the pilot subcarriers and the data subcarriers in the current interference suppression region.
  • the data signal estimation on each of the data subcarriers obtained above can be sent to a demodulation decoding device to complete the detection of the data.
  • the soft information metric L[,...'! corresponding to each bit obtained after demodulation can also be performed.
  • the soft information metric corresponding to each bit after adjustment is G(j)Li, --,, G(j)U FO where F is the number of bits contained in ( ⁇ ), indicating the _/ ⁇ Data
  • Corresponding channel coefficient estimated value of the first pilot subcarrier position, h p H () is a conjugate transposed matrix of p (), and () is the number corresponding to the data stream of the receiving end in the interference suppression region
  • the received signal on the pilot subcarriers is the transmit signal on the first pilot subcarrier corresponding to the data stream in the interference suppression region.
  • the estimated channel coefficient estimates of the pilot subcarriers and the data subcarrier positions used in the steps of the wideband co-channel interference noise estimation and interference suppression method are calculated by the following methods: get:
  • Step 1 For each pilot subcarrier corresponding to the data stream in the interference suppression area, the receiving end shares the received signal on the pilot subcarrier with the pilot signal sent by the transmitting end on the pilot subcarrier. Multiplying the yoke to obtain a channel coefficient estimation value of the pilot subcarrier position; the channel coefficient estimation value h P () corresponding to the first pilot subcarrier PsC(z) corresponding to the data stream in the interference suppression region is Get:
  • ⁇ ( ⁇ ( ( 6 )
  • ⁇ , ⁇ , /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region
  • ( ⁇ ) is the interference suppression at the receiving end
  • the received signal on the first pilot subcarrier corresponding to the data stream in the area is a pilot signal sent by the transmitting end on the first pilot subcarrier corresponding to the data stream in the interference suppression area. Conventional), indicating that the pair is conjugated. Because the correlation of the pilot signals on the same pilot subcarrier is relatively low, the interference signal of the adjacent cell pilot band on the pilot subcarrier can be filtered out by the above operation. Obtaining a relatively accurate channel coefficient estimation value.
  • Step 2 Corresponding to the data stream For each data subcarrier, the receiving end weights the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream in the interference suppression region as the channel coefficient estimation value of the data subcarrier position.
  • the interference suppression area the data stream corresponding to a first-data subcarriers denoted DsC (/), DsC () the position of the channel coefficient estimates /) according to equation (7) to give:
  • DsC the data stream corresponding to a first-data subcarriers
  • DsC the position of the channel coefficient estimates /) according to equation (7) to give:
  • the weight of ⁇ () is given, and the part (the weight of 0 " y can be 0, and other parameters have the meanings as described above.
  • the channel estimation unit partitioning when the channel coefficient estimation value of a certain data subcarrier position is calculated according to the formula (7), the pilot subcarrier positions corresponding to the data stream in the same channel estimation unit are The channel coefficient estimates give the same weight.
  • the manner of the above two embodiments may be combined. as follows:
  • the number of pilot subcarriers included in the set At the time-frequency, the closer the pilot subcarriers are to a certain data subcarrier, the stronger the channel correlation. Therefore, preferably, in calculating the weight used, " w ,
  • the weighted average of step 20 may be performed based on the interference noise estimation unit.
  • Each time-frequency two-dimensional resource block is used as an interference noise estimation unit, and each interference noise estimation unit includes at least one pilot sub-carrier.
  • the division of the channel estimation unit and the interference noise estimation unit in the same interference suppression region may be the same.
  • interference noise estimation unit division when the interference noise covariance matrix estimation value of a certain data subcarrier position is calculated according to formula (1), each pilot in the same interference noise estimation unit is used. The interference noise covariance matrix of the subcarrier position is given the same weight.
  • a further embodiment of performing interference noise estimation unit partitioning may be combined with the manner of the above two embodiments.
  • the estimated noise interference covariance matrix of each data subcarrier position corresponding to the data stream in the w interference noise estimation unit is equal, denoted as M m — D , and the receiving end calculates according to the following formula:
  • 1 , 0 ⁇ ⁇ 1 , where
  • the number of pilot subcarriers included the interference noise covariance matrix estimation value of a data subcarrier position according to formula (1), the interference noise covariance matrix of each pilot subcarrier position in the same interference noise estimation unit, Taking the same weight; and calculating the interference noise covariance matrix estimation value of each data subcarrier position in the same interference noise estimation unit, by taking the same set of weights, making the interference noise covariance of each data subcarrier position
  • the matrix estimates are the same. In the time-frequency region, the closer the pilot subcarriers are to a certain data subcarrier, the stronger the channel correlation.
  • the embodiment further provides a system for wideband co-channel interference noise estimation, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression manner.
  • OFDM Orthogonal Frequency Division Multiplexing
  • OFDMA Orthogonal Frequency Division Multiple Access
  • Interference noise estimation is performed on a data stream carried in the area, where the interference suppression area is a time-frequency two-dimensional resource block in the received data bearer area
  • the system includes: a first device, configured to: correspond to the data stream a weighted average of the product of the received signal on each pilot subcarrier corresponding to the data stream and the conjugate transposed of the received signal, as the estimated value of the interference noise covariance matrix of the data subcarrier position
  • a second device configured to: perform a diagonal loading on the estimated value of the interference noise covariance matrix of the data subcarrier position obtained by the first device for each data subcarrier corresponding to the data stream, As the interference noise covariance matrix of the data subcarrier position.
  • the first device is configured to calculate an interference noise covariance matrix estimate of the data subcarrier position according to equation (a):
  • the system for frequency interference noise estimation may further include a third device, wherein: the third device is configured to: divide the interference suppression region into one or more interference noise estimation units, each interference noise estimation unit is a time domain two a dimension resource block and including at least one pilot subcarrier and one data subcarrier; correspondingly, the first device calculates the interference noise covariance matrix estimation value of the data subcarrier position according to formula (a), and is the same interference noise
  • the product of the received signal on each pilot subcarrier in the estimation unit and its conjugate transpose is given the same weight.
  • the system for wideband co-channel interference noise estimation may further comprise a fourth device, wherein: the fourth device is set The interference suppression area is divided into interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least one pilot. And a data subcarrier, which is a positive integer; correspondingly, for each data subcarrier corresponding to the data stream, the first device is set to be on each pilot subcarrier corresponding to the data stream according to the following formula A weighted average of the product of the received signal and the conjugate transpose of the received signal as an estimate of the interference noise covariance matrix for the data subcarrier position: among them,
  • ⁇ ;) is the received signal of the receiving end on the first pilot subcarrier corresponding to the data stream; y p (if is the conjugate transposed matrix of (); fi ml is calculated - D , the first / interference noise estimation unit corresponding to each pilot subcarrier corresponding to the data stream (the weight of the product of (,
  • the second device is configured to diagonally load the interference noise covariance matrix estimation value of the data subcarrier position as the data, as the data Interference noise covariance matrix for subcarrier position:
  • R W _ D / is the interference noise covariance matrix of the _th data subcarrier position corresponding to the data stream in the interference suppression region
  • 1, ⁇
  • J is the interference suppression region
  • the number of data subcarriers corresponding to the data stream; ⁇ 0; is the estimated value of the interference noise covariance matrix of the _; data subcarrier positions corresponding to the data stream in the interference suppression region; ⁇ >0 ⁇ represents N xN 3 ⁇ 4
  • the diagonal matrix, 3 ⁇ 4, represents the number of receiving antennas at the receiving end.
  • the embodiment further provides a system for wideband co-channel interference suppression, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM)/Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression region.
  • OFDM Orthogonal Frequency Division Multiplexing
  • OFDMA Orthogonal Frequency Division Multiple Access
  • the system comprising: a first subsystem, configured to: according to the above-mentioned broadband co-channel interference noise The estimated system obtains the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream and the interference noise covariance matrix of each data subcarrier position; the second subsystem is set to: the data stream Corresponding to each data subcarrier, a weighted average of channel coefficient estimates of respective pilot subcarrier positions corresponding to the data stream is used as an estimated channel coefficient of the data subcarrier position; and a third subsystem is set to : for each data subcarrier corresponding to the data stream, according to the received signal on the data subcarrier, and the channel of the data subcarrier position Estimated value covariance matrix and the interference noise, the data signal estimate is calculated on the data subcarriers.
  • the system for the broadband inter-frequency interference suppression may further include a fourth subsystem: the fourth subsystem is configured to: calculate, for each data subcarrier corresponding to the data stream, the third subsystem After estimating the data signal on the data carrier, the estimated data signal is demodulated, and the soft information corresponding to each bit of the signal obtained after demodulation is adjusted, and the soft information corresponding to each bit after adjustment is respectively
  • ⁇ ⁇ , ⁇ , J J is the number of data subcarriers corresponding to the data stream in the interference suppression region
  • F is the jth data carrier corresponding to the data stream in the interference suppression region
  • the number of bits included in the data signal indicates the _/ ⁇ th data carrier pair corresponding to the data stream in the interference suppression region
  • M _ P W is the interference noise covariance matrix of the first pilot subcarrier position, where 1 is the inverse matrix of ⁇ ; () is the number corresponding to the data stream in the interference suppression
  • the system for wideband co-channel interference suppression may further include a fifth subsystem and a sixth subsystem: the fifth subsystem is configured to: divide the interference suppression region into K channel estimation units, where each channel estimation unit is a time domain two-dimensional resource block and including at least one pilot subcarrier and one data subcarrier, where f is a positive integer; the sixth subsystem is set to: for each data subcarrier corresponding to the data stream, The weighting average of the channel coefficient estimation values of the respective pilot subcarrier positions corresponding to the data stream is used as a step of estimating the channel coefficient of the data subcarrier position as follows:
  • the number of pilot subcarriers; the channel of the first pilot subcarrier position corresponding to the data stream in the interference suppression region Coefficient estimated value; and a kl is the ( ) given to the position of each pilot subcarrier in the channel estimation unit when calculating
  • ⁇ ” 1 , 0 ⁇ ⁇ 1 ,
  • indicates the number of pilot subcarriers included, and is in the right
  • each parameter is the same as that of the above embodiment, and it is assumed that the receiving end has obtained the received signal on each data subcarrier.
  • the example mainly shows how to further calculate the interference noise covariance matrix of the data subcarrier position in the case of different interference suppression region patterns and interference noise estimation unit partitioning. For the data signal estimation, see above, it will not be repeated.
  • Application example one is the same as that of the above embodiment, and it is assumed that the receiving end has obtained the received signal on each data subcarrier.
  • the example mainly shows how to further calculate the interference noise covariance matrix of the data subcarrier position in the case of different interference suppression region patterns and interference noise estimation unit partitioning. For the data signal estimation, see above, it will not be repeated.
  • the interference suppression region is divided into one interference noise estimation unit, and the index of the 20 pilot subcarriers included in the interference suppression region belongs to one pilot index set, and 1 to 20 belong to .
  • the estimated value of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit is ⁇ , with:
  • the suppression region is divided into two interference noise estimation units, and the 20 pilot subcarriers included in the interference suppression region belong to two pilot index sets, wherein 1 to 8 belong to ⁇ , and 9 to 20 belong to ⁇ 2 .
  • the estimated value of the interference noise covariance matrix of each data subcarrier position in the first interference noise estimation unit is ⁇ , with:
  • the suppression region is divided into two interference noise estimation units, and the 20 pilot subcarriers included in the interference suppression region belong to two pilot index sets, where: 1, 2, 5, 6, 9, 10, 13, 14 17, 18, 18 belong to ⁇ , and the rest of the pilot index belongs to ⁇ 2 .
  • the interference noise covariance matrix of each pilot subcarrier position in the first interference noise estimation unit is ⁇
  • the interference noise covariance matrix of each pilot subcarrier position in the second interference noise estimation unit is:
  • the interference noise covariance matrix estimation of the jth data subcarrier position in the interference suppression region The value is:
  • the interference suppression area in this example is the interference suppression area pattern 2. It contains 12 consecutive OFDM/OFDMA symbols in the time domain and 4 consecutive subcarriers in the frequency domain, which carries one data stream.
  • the interference suppression region is divided into one interference noise estimation unit, and the index of the 16 pilot subcarriers included in the interference suppression region belongs to one pilot index set, and all of 1 to 16 belong.
  • the interference suppression area in this example is the interference suppression area pattern 3. It contains 9 consecutive OFDM/OFDMA symbols in the time domain and 4 consecutive subcarriers in the frequency domain, which carries one data stream.
  • the interference suppression region is divided into one interference noise estimation unit, and the index of the 12 pilot subcarriers included in the interference suppression region belongs to one pilot index set, and 1 to 12 belong to each other.
  • the following uses the application examples 6 and 7 to illustrate the interference noise estimation method for carrying two data streams in an interference suppression region.
  • the channel coefficient estimate on the first pilot subcarrier corresponding to the first data stream is fi p ( ); the channel coefficient estimate on the first pilot subcarrier corresponding to the second data stream is fi p2 ⁇ >;
  • the pilot signal sent by the transmitting end on the first pilot subcarrier corresponding to the first data stream is the pilot signal sent by the transmitting end on the first pilot subcarrier corresponding to the second data stream is ⁇ receiving
  • the received signal received by the terminal on the first pilot subcarrier corresponding to the first data stream and the received signal received by the receiving end on the first pilot subcarrier corresponding to the second data stream are 2 ( ).
  • Application Example 6 As shown in FIG.
  • the interference suppression region in this application example is an interference suppression region pattern 4, which includes 15 consecutive OFDM/OFDMA symbols in the time domain and 4 consecutive subcarriers in the frequency domain, where Carry two data streams.
  • the interference suppression region is divided into an interference noise estimation unit, and the indexes of the 10 pilot subcarriers corresponding to each data stream included in the interference suppression region belong to two pilot index sets, where: 1 ⁇ 10 belongs to ⁇ , and 11 ⁇ 20 belongs to ⁇ 3 ⁇ 4.
  • the estimated value of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit is ⁇ 11 , having:
  • the above method for performing interference noise estimation for each data stream is the same.
  • the division of the interference noise estimation unit for different data streams may also be different.
  • Application Example 7 As shown in FIG. 9, the interference suppression region in this application example is an interference suppression region pattern five, which includes 6 consecutive OFDM/OFDMA symbols in the time domain and 6 consecutive subcarriers in the frequency domain, where Carry two data streams.
  • the interference suppression region is divided into one interference noise estimation unit, and the four pilot subcarrier indexes corresponding to each data stream included in the interference suppression region belong to one pilot index set ⁇ , that is, 1 ⁇ 4 belongs to ⁇ ,.
  • the estimated value of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit is ⁇ 11
  • the estimated values of the interference noise covariance matrix for each data subcarrier position within the interference noise estimation unit are ⁇ 12 , with:
  • the first data subcarrier position's interference noise covariance matrix is estimated to be UJ, for any size of the stem
  • the interference suppression region the interference noise covariance matrix estimated value of the _/ ⁇ data subcarrier positions is obtained after the diagonal loading, and the interference noise covariance matrix R M — D /;) of the data subcarrier position is obtained,
  • / represents the number of pilot subcarriers that the current data stream contains in the interference suppression region.
  • is a constant and satisfies the condition 0 ⁇ ⁇ 1.
  • the interference noise covariance matrix estimated value of the _/ ⁇ data subcarrier positions is subjected to diagonal loading to obtain the interference noise covariance matrix R M _ D C/;> of the data subcarrier position, and there are:
  • the following application examples mainly describe the soft information adjustment steps performed by the receiving end after demodulation.
  • Application Example 10 After the receiving end demodulates the transmitted signal on each data subcarrier, it is assumed that the transmitting end transmits F soft information in the signal transmitted on the first data subcarrier, respectively, A, ...
  • the soft information can be adjusted by using the channel gain, and the soft information metric corresponding to each bit after adjustment is Gm " , wherein G (_7) is calculated as C /) (R M - (_ /)) - Of course, the soft information is adjusted, namely:
  • the following application examples primarily illustrate broadband inter-frequency interference noise estimation and interference suppression using the method of the present invention. It should be noted that the following describes the case where only one data stream is carried in one interference suppression area. For the case of carrying multiple data streams simultaneously, for each data stream, interference is performed by the same method as the following application example. Suppress the reception of the signal.
  • Application Example 11 In this application example, the interference suppression region pattern 1 is taken as an example.
  • the channel estimation unit division method and the interference noise estimation unit are divided in the same manner, that is, included in each channel estimation unit.
  • the subcarriers are identical to the subcarriers included in each interference noise estimation unit.
  • the channel coefficient estimates for all data subcarrier locations within the channel estimation unit are:
  • the estimated value of the interference noise covariance matrix of all data subcarrier positions in the interference noise estimation unit is ⁇ , with:
  • R L ⁇ (0( ⁇ (0) Completion of channel estimation and interference noise covariance matrix estimation corresponding to data subcarrier DsC(/)
  • the receiver performs data detection, which has:
  • the soft information corresponding to each bit included in each data symbol s is obtained: , optionally, the soft information is adjusted accordingly, and the adjusted
  • the soft information metrics corresponding to the bits are (? ⁇ /; ⁇ ,... ⁇ /;) ⁇ . among them, Of course, the above soft information adjustment process can also be neglected, that is, the information after demodulation is directly sent. Go to the decoding module to complete the entire data receiving process.
  • the method and system for estimating the wideband co-channel interference of the present invention can obtain more accurate interference noise characteristics, which is advantageous for improving the performance of interference suppression and the accuracy of data detection.
  • the wideband co-channel interference suppression method and system of the present invention is based on relatively accurate interference noise characteristics, and can improve the performance of interference suppression and the accuracy of data detection.

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Radio Transmission System (AREA)
  • Noise Elimination (AREA)

Abstract

Methods for estimating broadband co-channel interference and noise and suppressing interference, and systems thereof, are disclosed. The method for estimating broadband co-channel interference and noise includes the following steps: for each data sub-carrier corresponding to a data flow, a weighted average for a product of a received signal on each pilot sub-carrier corresponding to the data flow and the conjugate transpose of the received signal is determined as an estimated value of an interference and noise covariance matrix in the position of the data sub-carrier; and for each data sub-carrier corresponding to the data flow, the result obtained by performing diagonal loading to the estimated value of the interference and noise covariance matrix in the position of the data sub-carrier is determined as the interference and noise covariance matrix in the position of the data sub-carrier. The present invention can acquire a more accurate interference and noise characteristics, thus being advantageous to improving the performance of interference rejection and the accuracy of data detection.

Description

一种宽带同频干扰噪声估计和干扰抑制的方法及相应系统 技术领域 本发明涉及通信领域, 尤其涉及到一种宽带同频干扰噪声估计和干扰抑 制的方法及相应系统。  TECHNICAL FIELD The present invention relates to the field of communications, and in particular, to a method and a corresponding system for wideband co-channel interference noise estimation and interference suppression.
背景技术 Background technique
多天线技术是无线移动通信领域中智能天线技术的一个重大突破, 该技 术可以在不增加带宽的情况下成倍地提高通信系统的容量和频谱利用率, 还 可以利用多径来减轻多径衰落, 并能有效地消除信道干扰, 提高信道的可靠 性, 降低误码率, 是新一代移动通信系统的关键技术。 多天线技术已经被广 泛地应用于长期演进 ( Long Term Evolution, LTE )和全球微波接入互操作 性 ( World Interoperability for Microwave Access, WiMAX )等多种无线宽带 系统中。 对于以蜂窝结构布置网络的无线通信来说, 相邻小区间的同频干扰是导 致通信质量下降的最重要因素之一, 如图 1所示。 因为干扰源为邻小区用户 在同时同频资源上发送的数据信号, 使得接收端必须要较为准确的估计出期 望数据的信道系数、 干扰信道系数或者干扰的特征才能够进行较为准确的数 据检测, 但是当相邻小区间的数据导频子载波在时频位置上重合时, 就给干 扰估计带来了巨大的困难。 因为干扰导频的重合会导致信道估计质量下降, 即信道估计本身携带干扰信息, 从而使得干扰噪声特征估计变得十分困难或 者说是非常不准确。 在此种情况下,常见的干扰抑制接收算法,比如最小均方误差(Minimum Mean- Squared Error , MMSE ) 或者干扰 p制合并 ( Interference Rejection Combination, IRC ) 算法的性能都会大打折扣。  Multi-antenna technology is a major breakthrough in smart antenna technology in the field of wireless mobile communications. This technology can multiply the capacity and spectrum utilization of communication systems without increasing bandwidth. It can also use multipath to mitigate multipath fading. It can effectively eliminate channel interference, improve channel reliability, and reduce bit error rate. It is a key technology for a new generation of mobile communication systems. Multi-antenna technology has been widely used in many wireless broadband systems such as Long Term Evolution (LTE) and World Interoperability for Microwave Access (WiMAX). For wireless communications in which the network is arranged in a cellular structure, co-channel interference between adjacent cells is one of the most important factors leading to a degradation in communication quality, as shown in FIG. Because the interference source is the data signal sent by the neighboring cell user on the same frequency resource, the receiving end must accurately estimate the channel coefficient, the interference channel coefficient or the interference feature of the desired data to perform more accurate data detection. However, when the data pilot subcarriers between adjacent cells coincide at the time-frequency position, it brings great difficulty to the interference estimation. Since the coincidence of the interference pilots leads to a degradation of the channel estimation quality, that is, the channel estimation itself carries the interference information, making the estimation of the interference noise characteristics very difficult or very inaccurate. In this case, the performance of common interference suppression receiving algorithms, such as Minimum Mean-Square Error (MMSE) or Interference Rejection Combination (IRC) algorithm, will be greatly reduced.
发明内容 本发明的目的是提供一种宽带同频干扰噪声估计的方法及相应系统, 以 解决相邻小区存在同频干扰时干扰噪声特征估计不准确的问题。 为了解决上述技术问题, 本发明提供了一种宽带同频干扰噪声估计的方 法, 应用于正交频分复用(OFDM)或正交频分多址 (OFDMA)系统的接收端, 在一干扰抑制区域内, 用该方法对其中承载的一个数据流进行干扰噪声估计 时, 该方法包括: 对该数据流对应的每一数据子载波, 将该数据流对应的各导频子载波上 的接收信号与该接收信号的共轭转置的乘积的加权平均, 作为该数据子载波 位置的干扰噪声协方差矩阵估计值; 以及 对该数据流对应的每一数据子载波, 将对该数据子载波位置的干扰噪声 协方差矩阵估计值进行对角加载后得到的结果, 作为该数据子载波位置的干 扰噪声协方差矩阵; 其中, 该干扰抑制区域为接收数据承载区域中的一时频二维资源块。 本发明的方法中, 所述对该数据流对应的每一数据子载波, 将该数据流 对应的各导频子载波上的接收信号与该接收信号的共轭转置的乘积的加权平 均, 作为该数据子载波位置的干扰噪声协方差矩阵估计值的步骤中, 釆用式 ( a )计算数据子载波位置的干扰噪声协方差矩阵估计值:
Figure imgf000004_0001
SUMMARY OF THE INVENTION It is an object of the present invention to provide a method and corresponding system for wideband co-channel interference noise estimation, The problem of inaccurate estimation of interference noise characteristics when there is co-channel interference in adjacent cells is solved. In order to solve the above technical problem, the present invention provides a method for wideband co-channel interference noise estimation, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference In the suppression region, when the method performs interference noise estimation on a data stream carried in the method, the method includes: receiving, on each pilot subcarrier corresponding to the data stream, the reception on each pilot subcarrier corresponding to the data stream a weighted average of the product of the signal and the conjugate transpose of the received signal as an interference noise covariance matrix estimate for the data subcarrier position; and for each data subcarrier corresponding to the data stream, the data subcarrier will be The interference noise covariance matrix estimated value of the position is diagonally loaded as the interference noise covariance matrix of the data subcarrier position; wherein the interference suppression region is a time-frequency two-dimensional resource block in the received data bearer region . In the method of the present invention, the weighted average of the product of the received signal on each pilot subcarrier corresponding to the data stream and the conjugate transposed on the received signal for each data subcarrier corresponding to the data stream, In the step of estimating the interference noise covariance matrix of the data subcarrier position, the equation (a) is used to calculate the interference noise covariance matrix estimate of the data subcarrier position:
Figure imgf000004_0001
其中, 为该干扰抑制区域中该数据流对应的第 个数据子载波位 置的干扰噪声协方差矩阵估计值, · = 1,· · · , J为该干扰抑制区域中该数据流 对应的数据子载波个数; = 1, · · ·,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; ^为计算该干扰抑制区域中该数据流对应的第 _;个数据子载 波位置的干扰噪声协方差矩阵估计值时赋予 ( ( i)) (yp (i )f 的权值,Wherein, the estimated value of the interference noise covariance matrix of the first data subcarrier position corresponding to the data stream in the interference suppression region, · = 1, · · · , J is the data corresponding to the data stream in the interference suppression region The number of carriers; = 1, · · ·, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region; ^ is to calculate the _th data corresponding to the data stream in the interference suppression region The interference noise covariance matrix estimate for the subcarrier position is given the weight of ( ( i ) ) ( y p (i )f ,
∑βϋ = \ ; W为所述接收端在该干扰抑制区域中该数据流对应的第 ,个导 频子载波上的接收信号, 以及 ( 为 ( 的共轭转置矩阵。 按式 (a)计算数据子载波位置的干扰噪声协方差矩阵估计值的步骤之前, 所述方法还包括: 将该干扰抑制区域划分为一个或多个干扰噪声估计单元, 每一干扰噪声 估计单元为一个时域二维资源块且其中包含至少一个导频子载波和一个数据 子载波; Σβ ϋ = \; W to the receiving end in the area of the interference suppression data stream corresponding to a first received signal on pilot subcarriers, and (for the (conjugate transposed matrix according to formula (a. Before the step of calculating an interference noise covariance matrix estimate of the data subcarrier position, the method further includes: The interference suppression area is divided into one or more interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least one pilot subcarrier and one data subcarrier therein;
按式 (a)计算数据子载波位置的干扰噪声协方差矩阵估计值时, 为同一干 扰噪声估计单元中各个导频子载波上的接收信号与其共轭转置的乘积, 赋予 相同的权值。 所述方法还包括:  When the interference noise covariance matrix estimation value of the data subcarrier position is calculated according to equation (a), the same weight is given to the product of the received signal on each pilot subcarrier in the same interference noise estimation unit and its conjugate transposed. The method further includes:
所述在一干扰抑制区域内, 用该方法对其中承载的一个数据流进行干扰 噪声估计时, 还将该干扰抑制区域划分为 个干扰噪声估计单元,每一干扰 噪声估计单元为一个时域二维资源块且其中包含至少一个导频子载波和一数 据子载波, 为正整数; 所述对该数据流对应的每一数据子载波, 将该数据流对应的各导频子载 波上的接收信号与该接收信号的共轭转置的乘积的加权平均, 作为该数据子 载波位置的干扰噪声协方差矩阵估计值的步骤中, 釆用式( b )计算数据子载 波位置的干扰噪声协方差矩阵估计值:  When the interference noise estimation is performed on a data stream carried by the method in an interference suppression region, the interference suppression region is further divided into an interference noise estimation unit, and each interference noise estimation unit is a time domain two. a dimension resource block and including at least one pilot subcarrier and a data subcarrier, which are positive integers; and each data subcarrier corresponding to the data stream, receiving on each pilot subcarrier corresponding to the data stream A weighted average of the product of the signal and the conjugate transpose of the received signal as a step of estimating the interference noise covariance matrix of the data subcarrier position, using equation (b) to calculate the interference noise covariance of the data subcarrier position Matrix estimate:
RR
Figure imgf000005_0001
其中,
Figure imgf000005_0001
among them,
O 第 m 个干扰噪声估计单元中该数据流对应的每一数据子载波位 置的干扰噪声协方差矩阵估计值, m = 1, 2, · · · ,Μ; /为一循环变量, / = 1, 2, · · · , ;  O The estimated value of the interference noise covariance matrix for each data subcarrier position corresponding to the data stream in the mth interference noise estimation unit, m = 1, 2, · · · , Μ; / is a cyclic variable, / = 1 , 2, · · · , ;
为第 /个干扰噪声估计单元包含的该数据流对应的导频子载波的索引 ζ· 的集合, = 1, · · ·, /, /为该干扰抑制区域中该数据流对应的导频子载波的个 数; a set of indices ζ· of the pilot subcarriers corresponding to the data stream included in the first interference noise estimation unit, = 1, · · ·, /, / is the pilot corresponding to the data stream in the interference suppression region The number of carriers;
^;)为所述接收端在该数据流对应的第 I个导频子载波上的接收信号; (广为 ( )的共轭转置矩阵; 以及 βηι为计算 icMD时,赋予第 /个干扰噪声估计单元中该数据流对应的各导 频子载波对应的 ( 与 ( 的乘积的权值, | |为
Figure imgf000006_0001
^;) is the received signal of the receiving end on the first pilot subcarrier corresponding to the data stream; (a wide () conjugate transposed matrix; When β η ι is calculated as ic MD , it is assigned to each pilot subcarrier corresponding to the data stream in the third interference noise estimation unit (the weight of the product of (, | | is
Figure imgf000006_0001
包含的导频子载波的个数。 本发明的方法中,按式 (b)计算 ^釆用的权值 ,中, / = 1,2,···,Μ , β 大于等于其他的权值。 本发明的方法中, 所述对该数据流对应的每一数据子载波, 将对该数据 子载波位置的干扰噪声协方差矩阵估计值进行对角加载后得到的结果, 作为 该数据子载波位置的干扰噪声协方差矩阵的步骤中, 釆用式(c)计算数据子 载波位置的干扰噪声协方差矩阵: RM_D 0) = RM_D (j) + βΑ ( c ) 其中, RW_D /)为该干扰抑制区域中该数据流对应的第 _;个数据子载波位 置的干扰噪声协方差矩阵, · = 1,··· , J为该干扰抑制区域中该数据流对应的 数据子载波个数; ≥0; 为该干扰抑制区域中该数据流对应的第 _;个 数据子载波位置的干扰噪声协方差矩阵估计值; β>0 Λ表示 N xN¾的对角 矩阵, 以及 ¾表示所述接收端的接收天线的数目。 本发明的方法中, 所述第 个数据子载波位置的干扰噪声协方差矩阵为:
Figure imgf000006_0002
The number of pilot subcarriers included. In the method of the present invention, the weights used for the calculation are calculated according to the formula (b), where / = 1, 2, ···, Μ, β are greater than or equal to other weights. In the method of the present invention, the data interference result obtained by diagonally loading the interference noise covariance matrix estimation value of the data subcarrier position is used as the data subcarrier position for each data subcarrier corresponding to the data stream. In the step of the interference noise covariance matrix, the interference noise covariance matrix of the data subcarrier position is calculated by using equation (c): R M _ D 0) = R M _ D (j) + βΑ ( c ) where R W _ D /) is the interference noise covariance matrix of the _th data subcarrier position corresponding to the data stream in the interference suppression region, · = 1,···, J is the data stream corresponding to the interference suppression region The number of data subcarriers; ≥ 0; is the estimated value of the interference noise covariance matrix of the _; data subcarrier positions corresponding to the data stream in the interference suppression region; β>0 Λ represents the diagonal matrix of N xN 3⁄4 And 3⁄4 denote the number of receiving antennas at the receiving end. In the method of the present invention, the interference noise covariance matrix of the first data subcarrier position is:
Figure imgf000006_0002
其中, 0≤ ≤1; tr( MD(_/》表示对矩阵 求迹; /为该干扰抑制区 域内该数据流对应的导频子载波的个数, 1为^„^^¾的单位矩阵。 本发明的方法中, 所述在一干扰抑制区域内, 用该方法对其中承载的一 个数据流进行干扰噪声估计时, 按以下方式计算该数据流对应的每一导频子 载波位置的信道系数估计值: 将该导频子载波上的接收信号与发送端在该导频子载波上发送的导频信 号的共轭相乘, 得到该导频子载波位置的信道系数估计值。 相应地, 本发明还提供了一种宽带同频干扰噪声估计的系统, 应用于正 交频分复用(OFDM)或正交频分多址 (OFDMA)系统的接收端, 在一干扰抑制 区域内对其中承载的一个数据流进行干扰噪声估计, 该干扰抑制区域为接收 数据承载区域中的一时频二维资源块, 该系统包括: 第一装置, 其设置为: 对该数据流对应的每一数据子载波, 将该数据流 对应的各导频子载波上的接收信号与该接收信号的共轭转置的乘积的加权平 均, 作为该数据子载波位置的干扰噪声协方差矩阵估计值; 以及 第二装置, 其设置为: 对该数据流对应的每一数据子载波, 将对所述第 一装置得到的该数据子载波位置的干扰噪声协方差矩阵估计值进行对角加载 后得到的结果, 作为该数据子载波位置的干扰噪声协方差矩阵。 本发明的系统中, 所述第一装置是设置为釆用式( a )计算出该数据子载 波位置的干扰噪声协方差矩阵估计值:
Figure imgf000007_0001
Where 0 ≤ ≤1; tr( MD (_/′ represents the trace of the matrix; / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region, 1 is ^„^^ 3⁄4 In the method of the present invention, when the interference noise estimation is performed on a data stream carried by the method in the interference suppression region, each pilot subcarrier position corresponding to the data stream is calculated in the following manner. Estimated value of the channel coefficient: Multiply the received signal on the pilot subcarrier with the conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier to obtain an estimated channel coefficient of the pilot subcarrier position. Correspondingly, the present invention also provides a system for wideband co-channel interference noise estimation, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression region. Performing interference noise estimation on a data stream carried therein, the interference suppression area being a time-frequency two-dimensional resource block in the received data bearer area, the system comprising: a first device, configured to: each corresponding to the data stream a data subcarrier, a weighted average of a product of a received signal on each pilot subcarrier corresponding to the data stream and a conjugate transposed of the received signal, as an interference noise covariance matrix estimation value of the data subcarrier position; And a second device, configured to: diagonally load the estimated value of the interference noise covariance matrix of the data subcarrier position obtained by the first device for each data subcarrier corresponding to the data stream As a result, the interference noise covariance matrix is the position of the data subcarrier. In the system of the present invention, the first device is configured to calculate an interference noise covariance matrix estimate of the data subcarrier position using equation (a):
Figure imgf000007_0001
其中, 为该干扰抑制区域中该数据流对应的第 个数据子载波位 置的干扰噪声协方差矩阵估计值, · = 1,· · · , J为该干扰抑制区域中该数据流 对应的数据子载波个数; = 1, · · ·,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; ^为计算该干扰抑制区域中该数据流对应的第 个数据子载 波位置的干扰噪声协方差矩阵估计值时赋予 ( ( i)) (yp (i )f 的权值,Wherein, the estimated value of the interference noise covariance matrix of the first data subcarrier position corresponding to the data stream in the interference suppression region, · = 1, · · · , J is the data corresponding to the data stream in the interference suppression region The number of carriers; = 1, · · ·, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region; ^ is to calculate the first data subcarrier corresponding to the data stream in the interference suppression region imparting ((I)) weight (y p (i) f the noise covariance matrix estimate interference position,
∑βϋ = \ ; W为所述接收端在该干扰抑制区域中该数据流对应的第 ,个导 频子载波上的接收信号, (广为 ( 的共轭转置矩阵。 所述系统还包括第三装置, 其中: 所述第三装置设置为: 将该干扰抑制区域划分为一个或多个干扰噪声估 计单元, 每一干扰噪声估计单元为一个时域二维资源块且其中包含至少一个 导频子载波和一个数据子载波; 所述第一装置按式 (a)计算数据子载波位置的干扰噪声协方差矩阵估计 值时, 为同一干扰噪声估计单元中各个导频子载波上的接收信号与其共轭转 置的乘积, 赋予相同的权值。 所述系统还包括第四装置, 其中: 所述第四装置设置为: 将该干扰抑制区域划分为 个干扰噪声估计单 元, 每一干扰噪声估计单元为一个时域二维资源块且其中包含至少一个导频 子载波和一数据子载波, 为正整数; 对该数据流对应的每一数据子载波, 所述第一装置是设置为按如下公式 将该数据流对应的各导频子载波上的接收信号与该接收信号的共轭转置的乘 积的加权平均, 作为该数据子载波位置的干扰噪声协方差矩阵估计值:
Figure imgf000008_0001
其中,
Σβ ϋ = \; W to the receiving end in the area of the interference suppression corresponding to the first data stream, the received signal on pilot subcarriers, (widely (conjugate transpose matrix of the system further. A third device is included, where: the third device is configured to: divide the interference suppression region into one or more interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least one a pilot subcarrier and a data subcarrier; the first device calculates the interference noise covariance matrix estimation value of the data subcarrier position according to the equation (a), and is the reception on each pilot subcarrier in the same interference noise estimation unit Signal and its conjugate Set the product, giving the same weight. The system further includes a fourth device, where: the fourth device is configured to: divide the interference suppression region into interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least a pilot subcarrier and a data subcarrier are positive integers; for each data subcarrier corresponding to the data stream, the first device is set to be on each pilot subcarrier corresponding to the data stream according to the following formula A weighted average of the product of the received signal and the conjugate transpose of the received signal as the estimated value of the interference noise covariance matrix for the data subcarrier position:
Figure imgf000008_0001
among them,
O 第 m 个干扰噪声估计单元中该数据流对应的每一数据子载波位 置的干扰噪声协方差矩阵估计值, m = 1, 2, · · · ,Μ;  O The estimated value of the interference noise covariance matrix for each data subcarrier position corresponding to the data stream in the mth interference noise estimation unit, m = 1, 2, · · ·, Μ;
/为一循环变量, / = 1,2,···,Μ; / is a loop variable, / = 1,2,···,Μ;
为第 /个干扰噪声估计单元包含的该数据流对应的导频子载波的索引 I 的集合, = 1,···,/, /为该干扰抑制区域中该数据流对应的导频子载波的个 数; a set of index I of pilot subcarriers corresponding to the data stream included in the first interference noise estimation unit, = 1, . . . , /, / is the pilot subcarrier corresponding to the data stream in the interference suppression region Number of
^;)为所述接收端在该数据流对应的第 I个导频子载波上的接收信号; yp (if为 ()的共轭转置矩阵; fiml为计算 —D时,赋予第 /个干扰噪声估计单元中该数据流对应的各导 ^;) is the received signal of the receiving end on the first pilot subcarrier corresponding to the data stream; y p (if is the conjugate transposed matrix of (); fi ml is calculated - D , the first Each channel corresponding to the data stream in the interference noise estimation unit
M  M
频子载波对应的 w与 ( 的乘积的权值, ΣΙΩ=1,0≤ ≤1, | |为 The weight of the product of w and (the ΣΙ Ω =1 =1 , 0 ≤ ≤ 1 , | |
1=1  1=1
包含的导频子载波的个数, 大于等于其他的权值。 本发明的系统中, 对该数据流对应的每一数据子载波, 所述第二装置是 设置为按如下公式将对该数据子载波位置的干扰噪声协方差矩阵估计值进行 对角加载后得到的结果, 作为该数据子载波位置的干扰噪声协方差矩阵: NI-D The number of pilot subcarriers included is greater than or equal to other weights. In the system of the present invention, for each data subcarrier corresponding to the data stream, the second device is configured to diagonally load the estimated value of the interference noise covariance matrix of the data subcarrier position according to the following formula: As a result, the interference noise covariance matrix as the data subcarrier position: NI-D
其中, RW_D /)为该干扰抑制区域中该数据流对应的第 _;个数据子载波位 置的干扰噪声协方差矩阵, · = 1,··· , J为该干扰抑制区域中该数据流对应的 数据子载波个数; ≥0; 为该干扰抑制区域中该数据流对应的第 _;个 数据子载波位置的干扰噪声协方差矩阵估计值; β>0 Λ表示 N xN¾的对角 矩阵, ¾表示所述接收端的接收天线的数目。 本发明的系统中, Where R W _ D /) is the interference noise covariance matrix of the _th data subcarrier position corresponding to the data stream in the interference suppression region, · = 1,···, J is the interference suppression region The number of data subcarriers corresponding to the data stream; ≥ 0; is the estimated value of the interference noise covariance matrix of the _; data subcarrier positions corresponding to the data stream in the interference suppression region; β>0 Λ represents N xN 3⁄4 The diagonal matrix, 3⁄4, represents the number of receiving antennas at the receiving end. In the system of the present invention,
述第二装置计算出的第 个数据子载波位置的干扰噪声协方差矩阵为:
Figure imgf000009_0001
The interference noise covariance matrix of the first data subcarrier position calculated by the second device is:
Figure imgf000009_0001
其中, 0≤ ≤1; tr( MD(_/》表示对矩阵 求迹; /为该干扰抑制区 域内该数据流对应的导频子载波的个数, 1为^„^^¾的单位矩阵。 所述系统还包括第五装置, 其中: 所述第五装置设置为按以下方式得到发送端在该数据流对应的每一导频 子载波位置的信道系数估计值并输出到所述第一装置: 将该导频子载波上的接收信号与发送端在该导频子载波上发送的导频信 号的共轭相乘, 得到该导频子载波位置的信道系数估计值。 Where 0 ≤ ≤1; tr( MD (_/′ represents the trace of the matrix; / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region, 1 is ^„^^ 3⁄4 The unit matrix further includes: a fifth device, wherein: the fifth device is configured to obtain, at the following manner, a channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream at the transmitting end and output the The first device: multiplies the received signal on the pilot subcarrier by a conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier to obtain an estimated channel coefficient of the pilot subcarrier position.
上述宽带同频干扰的估计方法可以得到较为准确的干扰噪声特征, 有利 于提高干扰抑制的性能和数据检测的准确性。 The above-mentioned wideband co-channel interference estimation method can obtain more accurate interference noise characteristics, which is beneficial to improve the performance of interference suppression and the accuracy of data detection.
本发明的另一目的是提供一种宽带同频干扰的抑制方法及相应系统, 以 解决相邻小区存在同频干扰时干扰抑制性能较差的问题。 为了解决上述技术问题, 本发明提供了一种宽带同频干扰抑制的方法, 应用于正交频分复用 (OFDM)或正交频分多址(OFDMA) 系统的接收端, 在一干扰抑制区域内, 用该方法对其中承载的一个数据流进行干扰抑制时, 该方法包括: 按如上所述的干扰噪声估计方法, 得到该数据流对应的各导频子载波位 置的信道系数估计值和各数据子载波位置的干扰噪声协方差矩阵; 对该数据流对应的每一数据子载波, 将该数据流对应的各个导频子载波 位置的信道系数估计值的加权平均, 作为该数据子载波位置的信道系数估计 值; 以及 对该数据流对应的每一数据子载波, 根据该数据子载波上的接收信号, 及该数据子载波位置的信道系数估计值和干扰噪声协方差矩阵, 计算得到该 数据子载波上的数据信号估计; 其中, 该干扰抑制区域为接收数据承载区域中的一时频二维资源块。 所述方法还包括: 对该数据流对应的每一数据子载波, 在计算得到该数据载波上的数据信 号估计后, 对该估计出的数据信号进行解调, 并将解调后得到的该信号各比 特对应的软信息 L[,…'!F进行调整, 调整后各比特对应的软信息分别为 G{j)Li,-,G{j)UF; 其中, ' = l, ···, J, J为该干扰抑制区域中该数据流对应的数据子载波个数;Another object of the present invention is to provide a method for suppressing broadband co-channel interference and a corresponding system to solve the problem of poor interference suppression performance when adjacent cells have co-channel interference. In order to solve the above technical problem, the present invention provides a method for wideband co-channel interference suppression, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression manner. In this area, when this method is used to perform interference suppression on a data stream carried in it, The method includes: obtaining, according to the interference noise estimation method, the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream, and the interference noise covariance matrix of each data subcarrier position; corresponding to the data stream a weighted average of the channel coefficient estimates of the respective pilot subcarrier positions corresponding to the data stream, as the channel coefficient estimation value of the data subcarrier position; and each data substring corresponding to the data stream Transmitting, according to the received signal on the data subcarrier, and the channel coefficient estimation value and the interference noise covariance matrix of the data subcarrier position, calculating a data signal estimate on the data subcarrier; wherein the interference suppression region is receiving A time-frequency two-dimensional resource block in the data bearer area. The method further includes: after calculating the data signal estimation on the data carrier for each data subcarrier corresponding to the data stream, demodulating the estimated data signal, and demodulating the obtained data signal Soft information corresponding to each bit of the signal L[,...'! F is adjusted, and the soft information corresponding to each bit after adjustment is G{j)Li, -, G{j)U F ; wherein, ' = l, ···, J, J is the data in the interference suppression area The number of data subcarriers corresponding to the stream;
F为该干扰抑制区域中该数据流对应的第 J个数据载波上的数据信号中包含 的比特个数, 表示该干扰抑制区域中该数据流对应的第 _/·个数据载波对 应 的 信 道 增 益 4 ( )(RW— (7) 或 者 信 干 噪 比 值 F is the number of bits included in the data signal on the Jth data carrier corresponding to the data stream in the interference suppression area, and indicates the channel gain corresponding to the _/th data carrier corresponding to the data stream in the interference suppression area. 4 ( ) (R W — (7) or signal to interference and noise ratio
2
Figure imgf000010_0001
2
Figure imgf000010_0001
hd {j)为该干扰抑制区域中该数据流对应的第 J个数据子载波位置的信道 系数估计值, /)为 4( ·;)的共轭转置矩阵, RM_D /;>为第 ·个数据子载波位 置的干扰噪声协方差矩阵, (RMD ( ))'1为 RM_D (_/·)的逆矩阵; = 1, ···,/, /为该干扰抑制区域中该数据流对应的导频子载波个数; ft ή 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道系数估计 值, ( 为 的共轭转置矩阵, M_PW为第 ,个导频子载波位置的干扰 噪声协方差矩阵, 1为^ 的逆矩阵; ()为所述接收端在该干 扰抑制区域中该数据流对应的第 ,个导频子载波上的接收信号, 以及 为 发送端在该干扰抑制区域中该数据流对应的第 , 个导频子载波上的发射信 号。 本发明的方法中, 所述在一干扰抑制区域内, 用该方法对其中承载的一 个数据流进行干扰抑制时, 还将该干扰抑制区域划分为 f个信道估计单元, 每一信道估计单元为一个时域二维资源块且其中包含至少一个导频子载波和 一个数据子载波, f为正整数; 所述对该数据流对应的每一数据子载波, 将该数据流对应的各个导频子 载波位置的信道系数估计值的加权平均, 作为该数据子载波位置的信道系数 估计值的步骤, 釆用的计算公式如下:
Figure imgf000011_0001
其中,
h d {j) is a channel coefficient estimation value of the Jth data subcarrier position corresponding to the data stream in the interference suppression region, /) is a conjugate transposed matrix of 4 (·;), R M _ D /; > is the interference noise covariance matrix of the first data subcarrier position, (R MD ( )) ' 1 is the inverse matrix of R M _ D (_/·); = 1, ···, /, / The number of pilot subcarriers corresponding to the data stream in the interference suppression region; ft ή is the channel coefficient estimation of the first pilot subcarrier position corresponding to the data stream in the interference suppression region Value, ( for the conjugate transposed matrix, M _ P W is the first, the interference noise covariance matrix of the pilot subcarrier position, 1 is the inverse matrix of ^; () is the receiving end in the interference suppression region a received signal on the first pilot subcarrier corresponding to the data stream, and a transmit signal on the first pilot subcarrier corresponding to the data stream in the interference suppression region by the transmitting end. And performing interference suppression on a data stream carried by the method in the interference suppression region, and dividing the interference suppression region into f channel estimation units, where each channel estimation unit is a time domain two-dimensional a resource block and including at least one pilot subcarrier and one data subcarrier, where f is a positive integer; and for each data subcarrier corresponding to the data stream, a channel of each pilot subcarrier location corresponding to the data stream The weighted average of the coefficient estimates, as a step of estimating the channel coefficients for the data subcarrier position, is calculated as follows:
Figure imgf000011_0001
among them,
为第 个信道估计单元中该数据流对应的每一数据子载波位置的信道 系数估计值, k = \,2,---,K  The channel coefficient estimate for each data subcarrier location corresponding to the data stream in the first channel estimation unit, k = \, 2, ---, K
/为一循环变量, ! = 1,2,···,Κ / is a loop variable, ! = 1,2,···,Κ
Ω,为第 个信道估计单元包含的该干扰抑制区域中该数据流对应的导频 子载波的索引 ,的集合, = 1,···,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道 系数估计值;  Ω, which is a set of indices of pilot subcarriers corresponding to the data stream in the interference suppression region included in the first channel estimation unit, = 1, . . . , /, / is corresponding to the data stream in the interference suppression region The number of pilot subcarriers; the estimated channel coefficient of the first pilot subcarrier position corresponding to the data stream in the interference suppression region;
Okl为计算 时,赋予第 /个信道估计单元中各导频子载波位置的 ίιΡ ()的O kl is calculated, giving the first / channel estimation unit for each pilot sub-carrier position ίι Ρ () of
Κ Κ
权值, ∑|Ω=1,0≤ ≤1, | |表示 包含的导频子载波的个数, 且在权 Weight, ∑| Ω =1 , 0 ≤ ≤1 , | | indicates the number of pilot subcarriers included, and is in the right
1=1 值 中, l = \, 2 , K , 大于等于其他的权值。 1=1 In the value, l = \, 2, K, is greater than or equal to other weights.
相应地, 本发明还提供了一种宽带同频干扰抑制的系统, 应用于正交频 分复用 (OFDM )或正交频分多址(OFDMA )系统的接收端, 在一干扰抑制 区域内对其中承载的一个数据流进行干扰抑制, 该干扰抑制区域为接收数据 承载区域中的一时频二维资源块, 该系统包括: 第一子系统, 其设置为: 按如上所述的宽带同频干扰噪声估计的系统相 同的方式, 得到该数据流对应的各导频子载波位置的信道系数估计值和各数 据子载波位置的干扰噪声协方差矩阵; 第二子系统, 其设置为: 对该数据流对应的每一数据子载波, 将该数据 流对应的各个导频子载波位置的信道系数估计值的加权平均, 作为该数据子 载波位置的信道系数估计值; 以及 第三子系统, 其设置为: 对该数据流对应的每一数据子载波, 根据该数 据子载波上的接收信号, 及该数据子载波位置的信道系数估计值和干扰噪声 协方差矩阵, 计算得到该数据子载波上的数据信号估计。 所述系统还包括第四子系统, 其中: 所述第四子系统设置为: 对该数据流对应的每一数据子载波, 在所述第 三子系统计算得到该数据载波上的数据信号估计后, 对该估计出的数据信号 进行解调, 并将解调后得到的该信号各比特对应的软信息 进行调整, 调整后各比特对应的软信息分别为
Figure imgf000012_0001
; 其中, ' = l, · · · , J , J为该干扰抑制区域中该数据流对应的数据子载波个数;
Correspondingly, the present invention also provides a system for wideband co-channel interference suppression, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression region. Interference suppression is performed on a data stream carried in the received data bearer region, the system includes: a first subsystem, configured to: In the same manner as the system for interference noise estimation, the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream and the interference noise covariance matrix of each data subcarrier position are obtained; the second subsystem is set to: Each data subcarrier corresponding to the data stream, a weighted average of channel coefficient estimates of respective pilot subcarrier positions corresponding to the data stream, as a channel coefficient estimate of the data subcarrier position; and a third subsystem, Set to: each data subcarrier corresponding to the data stream, according to the received signal on the data subcarrier, and the location of the data subcarrier Channel coefficient estimates and the covariance matrix of interference noise, the data signal is calculated on the data subcarrier estimated. The system further includes a fourth subsystem, wherein: the fourth subsystem is configured to: calculate, for each data subcarrier corresponding to the data stream, a data signal estimate on the data carrier calculated by the third subsystem After that, the estimated data signal is demodulated, and the soft information corresponding to each bit of the signal obtained after demodulation is adjusted, and the soft information corresponding to each bit after adjustment is respectively
Figure imgf000012_0001
Where ' = l, · · · , J , J are the number of data subcarriers corresponding to the data stream in the interference suppression region;
F为该干扰抑制区域中该数据流对应的第 J个数据载波上的数据信号中包含 的比特个数, 表示该干扰抑制区域中该数据流对应的第 _/·个数据载波对F is the number of bits included in the data signal on the Jth data carrier corresponding to the data stream in the interference suppression area, and indicates the _/·th data carrier pair corresponding to the data stream in the interference suppression area.
( )(RW— (7) 或 者 信 干 噪 比 值
Figure imgf000012_0002
hd {j)为该干扰抑制区域中该数据流对应的第 j个数据子载波位置的信道 系数估计值, /)为 4( ·;)的共轭转置矩阵, RM_D /;>为第 ·个数据子载波位 置的干扰噪声协方差矩阵, (RMD ( ))'1为 RM_D (_/·)的逆矩阵; = 1, ···,/, /为该干扰抑制区域中该数据流对应的导频子载波个数; ft ή 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道系数估计 值, 为 的共轭转置矩阵, M_PW为第 ,个导频子载波位置的干扰 噪声协方差矩阵, 1为^ 的逆矩阵; ()为所述接收端在该干 扰抑制区域中该数据流对应的第 ,个导频子载波上的接收信号, 为发送 端在该干扰抑制区域中该数据流对应的第 I个导频子载波上的发射信号。 所述系统还包括第五子系统和第六子系统, 其中: 所述第五子系统设置为: 将该干扰抑制区域划分为 K个信道估计单元, 每一信道估计单元为一个时域二维资源块且其中包含至少一个导频子载波和 一个数据子载波, f为正整数; 所述第六子系统设置为: 对该数据流对应的每一数据子载波, 釆用如下 公式将该数据流对应的各个导频子载波位置的信道系数估计值的加权平均, 作为该数据子载波位置的信道系数估计值:
Figure imgf000013_0001
( ) (R W — (7) or signal to interference and noise ratio
Figure imgf000012_0002
h d {j) is an estimated channel coefficient of the jth data subcarrier position corresponding to the data stream in the interference suppression region, /) is a conjugate transposed matrix of 4 (·;), R M _ D /; > is the interference noise covariance matrix of the first data subcarrier position, (R MD ( )) ' 1 is the inverse matrix of R M _ D (_/·); = 1, ···, /, / The number of pilot subcarriers corresponding to the data stream in the interference suppression region; ft ή is the channel coefficient estimation value of the first pilot subcarrier position corresponding to the data stream in the interference suppression region, and is a conjugate rotation The matrix, M _ P W is the interference noise covariance matrix of the first pilot subcarrier position, where 1 is the inverse matrix of ^; () is the number corresponding to the data stream in the interference suppression region of the receiving end, The received signal on the pilot subcarrier is the transmit signal on the first pilot subcarrier corresponding to the data stream in the interference suppression region. The system further includes a fifth subsystem and a sixth subsystem, wherein: the fifth subsystem is configured to: divide the interference suppression region into K channel estimation units, each channel estimation unit is a time domain two-dimensional a resource block and including at least one pilot subcarrier and one data subcarrier, where f is a positive integer; the sixth subsystem is set to: for each data subcarrier corresponding to the data stream, using the following formula to the data The weighted average of the channel coefficient estimates for the respective pilot subcarrier locations corresponding to the stream, as the channel coefficient estimates for the data subcarrier locations:
Figure imgf000013_0001
/=1 ieQf 其中,  /=1 ieQf where,
为第 个信道估计单元中该数据流对应的每一数据子载波位置的信道 系数估计值, k = \,2,---,K  The channel coefficient estimate for each data subcarrier location corresponding to the data stream in the first channel estimation unit, k = \, 2, ---, K
/为一循环变量, ! = 1,2,···,Κ / is a loop variable, ! = 1,2,···,Κ
Ω,为第 个信道估计单元包含的该干扰抑制区域中该数据流对应的导频 子载波的索引 ,的集合, = 1,···,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道 系数估计值; akl为计算 时,赋予第 /个信道估计单元中各导频子载波位置的 ίΐΡ ()的Ω, which is a set of indices of pilot subcarriers corresponding to the data stream in the interference suppression region included in the first channel estimation unit, = 1, . . . , /, / is corresponding to the data stream in the interference suppression region The number of pilot subcarriers; the channel of the first pilot subcarrier position corresponding to the data stream in the interference suppression region Coefficient estimate; a kl is the ΐ Ρ () assigned to the position of each pilot subcarrier in the channel estimation unit when calculating
Κ Κ
权值, ∑|Ω=1,0≤ ≤1, | |表示 包含的导频子载波的个数, 且在权 Weight, ∑| Ω =1 , 0 ≤ ≤1 , | | indicates the number of pilot subcarriers included, and is in the right
1=1  1=1
值 0¾中, 1 = \,2,···,Κ , 大于等于其他的权值。 In value 03⁄4, 1 = \, 2,···, Κ , is greater than or equal to other weights.
上述宽带同频干扰的抑制方法基于较为准确的干扰噪声特征, 可以提高 干扰抑制的性能和数据检测的准确性。 The above-mentioned method for suppressing broadband co-channel interference is based on relatively accurate interference noise characteristics, which can improve the performance of interference suppression and the accuracy of data detection.
附图概述 图 1是现有技术中相邻多小区的示意图; 图 2是本发明实施例宽带同频干扰噪声估计和干扰抑制方法的流程图; 图 3至图 5分别是对干扰抑制区域样式一进行划分的 3种方式的示意图, 图加粗线圈表示导频子载波, 细线圈表示数据子载波, 后面的图 6和图 7同 此; 图 6是对干扰抑制区域样式二进行划分的一种方式的示意图; 图 7是对干扰抑制区域样式三进行划分的一种方式的示意图; 图 8是对干扰抑制区域样式四进行划分的一种方式的示意图, 图中粗线 圈表示数据流 1对应的导频子载波, 点划线圈表示数据流 2对应的导频子载 波, 细线圈表示数据子载波, 后面的图 9同此; 图 9是对干扰抑制区域样式五进行划分的一种方式的示意图。 BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is a schematic diagram of a neighboring multi-cell in the prior art; FIG. 2 is a flowchart of a method for estimating and suppressing interference of wideband co-channel interference according to an embodiment of the present invention; FIG. 3 to FIG. A schematic diagram of three ways of dividing, the bold coil represents the pilot subcarrier, the thin coil represents the data subcarrier, and the following FIG. 6 and FIG. 7 are the same; FIG. 6 is a division of the interference suppression region pattern 2 FIG. 7 is a schematic diagram of a manner of dividing the interference suppression region pattern three; FIG. 8 is a schematic diagram of a manner of dividing the interference suppression region pattern four, where the thick coil indicates the data stream 1 corresponding The pilot subcarrier, the dotted coil represents the pilot subcarrier corresponding to the data stream 2, and the thin coil represents the data subcarrier, which is the same as FIG. 9; FIG. 9 is a manner of dividing the interference suppression region pattern five. Schematic diagram.
本发明的较佳实施方式 为使本发明的目的、 技术方案和优点更加清楚明白, 下文中将结合附图 对本发明的实施例进行详细说明。 需要说明的是, 在不冲突的情况下, 本申 请中的实施例及实施例中的特征可以相互任意组合。 本实施例宽带同频干扰的估计和抑制方法应用于正交频分复用The embodiments of the present invention will be described in detail below with reference to the accompanying drawings. It should be noted that, in the case of no conflict, the features in the embodiments and the embodiments in the present application may be arbitrarily combined with each other. The method for estimating and suppressing wideband co-channel interference in this embodiment is applied to orthogonal frequency division multiplexing
(Orthogonal Frequency Division Multiplexing , OFDM)或正交频分多址 (Orthogonal Frequency Division Multiple Access , OFDMA)系统。 文中的发送 端可以是基站和中继站等控制设备, 也可以是手机、 笔记本电脑、 以及手持 电脑等终端设备。 类似地, 接收端用于接收发送端的数据信号, 接收端可以 是手机、 笔记本电脑、 以及手持电脑等终端设备, 也可以是基站和中继站等 控制设备。 接收端将接收数据承载区域划分为一个或多个干扰抑制区域, 每一干扰 抑制区域为帧或半帧结构中的一个时频二维资源块, 即每一个干扰抑制区域 在时间上包含多个连续的 OFDM/OFDMA符号,在频域上包含多个连续的子 载波。 接收数据承载区域可能包括一个时频二维资源块, 也可能包括多个分 离的时频二维资源块, 在本实施例中, 将其中的每一个独立的时频二维资源 块作为一个干扰抑制区域。 当然, 在其他实施例中, 接收数据承载区域中的 相对独立的各个时频二维护资源也可以被进一步划分为多个干扰抑制区域。 在 OFDM/OFDMA系统中,上述干扰抑制区域可以承载一个或多个数据 流, 每一数据流对应一个或多个数据子载波和导频子载波, 不同数据流对应 的导频子载波不同。 如图 2所示, 在每一个干扰抑制区域内, 按本实施例方法对其中承载的 一个数据流进行宽带同频干扰噪声估计和干扰抑制时, 该方法包括: 步骤 10、 对该数据流对应的每一数据子载波, 将该数据流对应的各导频 子载波上的接收信号与该接收信号的共轭转置的乘积的加权平均, 作为该数 据子载波位置的干扰噪声协方差矩阵估计值; 用 DsC( )表示该干扰抑制区域中该数据流对应的第 j 个数据子载波, 7 = 1,···, J , 则 DsC(/)对应的干扰噪声协方差矩阵估计值^ 按式(1 )得
Figure imgf000015_0001
(Orthogonal Frequency Division Multiplexing, OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system. The transmitting end in the text may be a control device such as a base station and a relay station, or may be a terminal device such as a mobile phone, a notebook computer, or a handheld computer. Similarly, the receiving end is configured to receive the data signal of the transmitting end, and the receiving end may be a terminal device such as a mobile phone, a notebook computer, and a handheld computer, or may be a control device such as a base station and a relay station. The receiving end divides the received data bearer area into one or more interference suppression areas, and each interference suppression area is a time-frequency two-dimensional resource block in a frame or a field structure, that is, each interference suppression area includes multiple times in time. A continuous OFDM/OFDMA symbol comprising a plurality of consecutive subcarriers in the frequency domain. The receiving data bearer area may include a time-frequency two-dimensional resource block, and may also include a plurality of separate time-frequency two-dimensional resource blocks. In this embodiment, each of the independent time-frequency two-dimensional resource blocks is used as an interference. Suppress area. Certainly, in other embodiments, the relatively independent time-frequency two maintenance resources in the received data bearer area may be further divided into multiple interference suppression areas. In an OFDM/OFDMA system, the interference suppression area may carry one or more data streams, and each data stream corresponds to one or more data subcarriers and pilot subcarriers, and different pilot streams corresponding to different data streams are different. As shown in FIG. 2, in each interference suppression region, when performing broadband wide-band interference noise estimation and interference suppression on a data stream carried by the method in this embodiment, the method includes: Step 10: Corresponding to the data stream Weighted average of the product of each data subcarrier, the received signal on each pilot subcarrier corresponding to the data stream and the conjugate transposed of the received signal, as the interference noise covariance matrix estimation of the data subcarrier position Value; DsC( ) is used to represent the jth data subcarrier corresponding to the data stream in the interference suppression region, 7 = 1,···, J , then the estimated value of the interference noise covariance matrix corresponding to DsC(/) ^ Formula (1)
Figure imgf000015_0001
其中, 为计算 DsC(/)对应的 时, 赋予 (_;)与 ( (7)广的乘积的 权值, ∑y¾ =1 , 部分权值可以为 0; W为接收端在该干扰抑制区域中该Where, in order to calculate the correspondence of DsC(/), the product of (_;) and ((7) is given Weight, ∑y3⁄4 =1 , part of the weight may be 0; W is the receiving end in the interference suppression area
7=1 7=1
数据流对应的第 ,个导频子载波上的接收信号, = ι,···,/, /为该干扰抑制区 域中该数据流对应的导频子载波的个数, J为该干扰抑制区域中该数据流对 应的数据子载波个数, 以及 ( ^表示矩阵 ( (· 的共轭转置。 步骤 20、 对该数据流对应的每一数据子载波, 将对该数据子载波位置的 干扰噪声协方差矩阵估计值进行对角加载后得到的结果, 作为该数据子载波 位置的干扰噪声协方差矩阵; The received signal on the first pilot subcarrier corresponding to the data stream, = ι,···, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region, and J is the interference suppression The number of data subcarriers corresponding to the data stream in the region, and (^ represents a matrix (conjugate transpose of (·). Step 20, for each data subcarrier corresponding to the data stream, the location of the data subcarrier will be The result obtained by diagonally loading the estimated value of the interference noise covariance matrix as the interference noise covariance matrix of the data subcarrier position;
即 DsC(/)位置的干扰噪声协方差矩阵 RW_D /)式(2)得到:
Figure imgf000016_0001
That is, the interference noise covariance matrix R W _ D /) of the DsC(/) position is obtained by the equation (2):
Figure imgf000016_0001
其中, ≥0, β≥0 , Λ表示 N xN¾的对角矩阵, 即除了对角线上元素 不为 0外, 其它位置的元素都为 0, N¾表示接收天线的数目。 较佳地, Wherein, ≥0, β≥0, Λ represents N xN ¾ diagonal matrix, i.e., in addition to the non-zero diagonal elements, the elements of the other positions are 0, N ¾ denotes the number of receive antennas. Preferably,
RwD O) = RM_D O) + 7„ ( 3 ) 其中, 。≤ ≤l, 且 由 /决定; tr( MD( 》表示对矩阵 求迹, 即 表示矩阵 MD (j)中所有对角线元素的累加和; /为该干扰抑制区域内该数据 流对应的导频子载波的个数, I为 N X N¾的单位矩阵。 通过以上两步, 接收端已经完成了对该干扰抑制区域的宽带同频干扰噪 声的估计。 对数据承载区域内的各干扰抑制区域均按上述方法计算后, 就完 成了对该数据承载区域的宽带同频干扰噪声估计。 步骤 30、 对该数据流对应的每一数据子载波, 根据该数据子载波上的接 收信号, 及该数据子载波位置的信道系数估计值和干扰噪声协方差矩阵, 计 算得到该数据子载波上的数据信号估计。 本步骤的运算是常规运算。 例如, 干扰抑制区域中该数据流对应的数据 子载波 DsC()对应的数据信号估计 ^ 通过以下方式计算得到: 当 ^( 表示为列向量时,
Figure imgf000017_0001
R wD O) = R M _ D O) + 7 „ ( 3 ) where ≤ ≤ l, and determined by /; tr( MD ( ′′ indicates the trace of the matrix, ie the matrix MD (j), an accumulated sum of all diagonal elements; / interference area for the data stream corresponding to the number of pilot subcarriers is suppressed, I is a unit matrix NXN ¾ above two steps, the receiving end has been completed. The estimation of the broadband co-channel interference noise in the interference suppression region is performed. After the interference suppression regions in the data bearer region are calculated according to the above method, the broadband co-channel interference noise estimation for the data bearer region is completed. And calculating, on each data subcarrier corresponding to the data stream, the data on the data subcarrier according to the received signal on the data subcarrier, and the channel coefficient estimation value and the interference noise covariance matrix of the data subcarrier position. Signal estimation The operation of this step is a conventional operation. For example, the data signal corresponding to the data subcarrier DsC() corresponding to the data stream in the interference suppression region is calculated by the following method: When ^ (represented as column direction) Time,
Figure imgf000017_0001
当 4( ·)表示为行向量时,
Figure imgf000017_0002
When 4 ( · ) is represented as a row vector,
Figure imgf000017_0002
其中, 4( ·)为数据子载波 DsC(/)对应的信道系数估计值, 为 4( ·) 的共辄转置, 》表示对 4( 的元素取共辄, (Rw—。(_/))— 1为 Rw— )的 逆矩阵, 以及 (_/·)为 DsC()上的接收信号。 本实施例中, (_/)表示为列向 量, 如 (7·)表示为行向量, 上述公式需做适应性变化, 不再赘述。 对于同一干扰抑制区域中的各数据流, 均可按上述步骤得到对应的数据 计算, 当然具体的权值选择可以是不同的。 如前所述, 上述步骤中该数据流 对应的导频子载波和数据子载波均指当前干扰抑制区域内的导频子载波和数 据子载波。 可以将上述得到的各数据子载波上的数据信号估计送到解调译码装置, 完成数据的检测。 Where 4(·) is the estimated channel coefficient corresponding to the data subcarrier DsC(/), which is a conjugated transposition of 4(·), ′′ indicates that the element of 4(取) is 辄, (R w —.(_ /)) - 1 is the inverse matrix of R w — ), and (_/·) is the received signal on DsC(). In this embodiment, (_/) is represented as a column vector, and if (7·) is expressed as a row vector, the above formula needs to be adaptively changed, and will not be described again. For each data stream in the same interference suppression area, the corresponding data calculation can be obtained according to the above steps. Of course, the specific weight selection can be different. As described above, the pilot subcarriers and the data subcarriers corresponding to the data stream in the foregoing steps all refer to the pilot subcarriers and the data subcarriers in the current interference suppression region. The data signal estimation on each of the data subcarriers obtained above can be sent to a demodulation decoding device to complete the detection of the data.
在对 ( ·)解调后, 还可以对该 ( ·)解调后得到的各比特对应的软信息度 量 L[,…'! ^进行调整, 调整后的各比特对应的软信息度量分别为 G(j)Li,--,G(j)UF O 其中, F为 ( ·)中包含的比特个数, 表示第 _/·个数据After (·) demodulation, the soft information metric L[,...'! corresponding to each bit obtained after demodulation can also be performed. ^Adjustment, the soft information metric corresponding to each bit after adjustment is G(j)Li, --,, G(j)U FO where F is the number of bits contained in (·), indicating the _/· Data
Η(υ ΐι—1^ ) 或者信干噪比值 在完成调整之后, 再将经过调整的所
Figure imgf000017_0003
Η(υ ΐι— 1 ^ ) or the signal to noise ratio after the adjustment is completed, then the adjusted place
Figure imgf000017_0003
有软信息送入译码器完成数据译码。 其中, 为 (·)的共轭转置矩阵, = 1, ···,/, /为该干扰抑制区域中该数据流对应的导频子载波个数; 为该 干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道系数估计值, hp H ()为 p ()的共轭转置矩阵, ()为接收端在该干扰抑制区域中该数据流 对应的第 ,个导频子载波上的接收信号, 为发送端在该干扰抑制区域中 该数据流对应的第 ,个导频子载波上的发射信号。 其他参数定义同上, 在此 不再进行赘述。 本实施例中, 上述宽带同频干扰噪声估计和干扰抑制方法的步骤中用到 的导频子载波和数据子载波位置的信道系数估计值^ W及 ^( ·;),可以通过以 下方式计算得到: Soft information is sent to the decoder to complete the data decoding. Wherein, the conjugate transposed matrix of (·), = 1, ···, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region; the data stream in the interference suppression region Corresponding channel coefficient estimated value of the first pilot subcarrier position, h p H () is a conjugate transposed matrix of p (), and () is the number corresponding to the data stream of the receiving end in the interference suppression region, The received signal on the pilot subcarriers is the transmit signal on the first pilot subcarrier corresponding to the data stream in the interference suppression region. Other parameters are defined as above, and will not be described here. In this embodiment, the estimated channel coefficient estimates of the pilot subcarriers and the data subcarrier positions used in the steps of the wideband co-channel interference noise estimation and interference suppression method are calculated by the following methods: get:
步骤一, 对该干扰抑制区域中该数据流对应的每一导频子载波, 接收端 将该导频子载波上的接收信号与发送端在该导频子载波上发送的导频信号的 共轭相乘, 得到该导频子载波位置的信道系数估计值; 该干扰抑制区域中该数据流对应的第 I个导频子载波 PsC(z)对应的信道 系数估计值hP ()由下式得到: Step 1: For each pilot subcarrier corresponding to the data stream in the interference suppression area, the receiving end shares the received signal on the pilot subcarrier with the pilot signal sent by the transmitting end on the pilot subcarrier. Multiplying the yoke to obtain a channel coefficient estimation value of the pilot subcarrier position; the channel coefficient estimation value h P () corresponding to the first pilot subcarrier PsC(z) corresponding to the data stream in the interference suppression region is Get:
^( =^( ( (6) 其中, = ι,···,/, /为该干扰抑制区域中该数据流对应的导频子载波的个 数, (Ζ)为接收端在该干扰抑制区域中该数据流对应的第 ,个导频子载波上 的接收信号, 为发送端在该干扰抑制区域中该数据流对应的第 ,个导频 子载波上发送的导频信号 (两端可以约定) , 表示对 取共轭。 因为相邻小区在同一导频子载波上的导频信号相关性比较低, 通过上述 运算, 可以滤除导频子载波上相邻小区导频带来的干扰信号, 得到较为准确 的信道系数估计值。 进而, 基于各导频子载波位置的信道系数估计值的加权 平均得到的数据子载波位置的信道系数估计值也较为准确。 步骤二, 对该数据流对应的每一数据子载波, 接收端将该干扰抑制区域 中该数据流对应的各导频子载波位置的信道系数估计值的加权平均, 作为该 数据子载波位置的信道系数估计值; 该干扰抑制区域中该数据流对应的第 ·个数据子载波记为 DsC(/),DsC() 位置的信道系数估计值 /)按式(7)得到:
Figure imgf000018_0001
其中, 为计算 DsC()对应的 /)时, 赋予^ ()的权值, 部 分 (0的权值《y可以为 0, 其他参数含义如上文所述。 接收端可以将该干扰抑制区域再划分为 f个时频二维资源块, =1,2,...; 每个时频二维资源块作为一个信道估计单元, 每一信道估计单元中包括至少 一个导频子载波和一个数据子载波。 在进行信道估计单元划分的一实施例中, 在按公式( 7 )计算某个数据子 载波位置的信道系数估计值时, 为同一信道估计单元中该数据流对应的各个 导频子载波位置的信道系数估计值赋予的权值相同。 在进行信道估计单元划分的另一实施例中, 在按公式(7 )计算同一信道 估计单元中该数据流对应的各个数据子载波位置的信道系数估计值时, 取一 组相同的权值《y., = ν··,/, 7 = 1,···, J, 得到的该数据流对应的各数据子载波 位置的信道系数估计值相同。 在进行信道估计单元划分的又一实施例, 可以结合上述两个实施例的方 式。 如下:
^( =^( ( 6 ) where = ι,···, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region, (Ζ) is the interference suppression at the receiving end The received signal on the first pilot subcarrier corresponding to the data stream in the area is a pilot signal sent by the transmitting end on the first pilot subcarrier corresponding to the data stream in the interference suppression area. Conventional), indicating that the pair is conjugated. Because the correlation of the pilot signals on the same pilot subcarrier is relatively low, the interference signal of the adjacent cell pilot band on the pilot subcarrier can be filtered out by the above operation. Obtaining a relatively accurate channel coefficient estimation value. Further, the channel coefficient estimation value of the data subcarrier position obtained based on the weighted average of the channel coefficient estimation values of the pilot subcarrier positions is also relatively accurate. Step 2: Corresponding to the data stream For each data subcarrier, the receiving end weights the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream in the interference suppression region as the channel coefficient estimation value of the data subcarrier position. The interference suppression area, the data stream corresponding to a first-data subcarriers denoted DsC (/), DsC () the position of the channel coefficient estimates /) according to equation (7) to give:
Figure imgf000018_0001
Wherein, in order to calculate /) corresponding to DsC(), the weight of ^() is given, and the part (the weight of 0 " y can be 0, and other parameters have the meanings as described above. The receiving end may further divide the interference suppression area into f time-frequency two-dimensional resource blocks, =1, 2, . . .; each time-frequency two-dimensional resource block is used as a channel estimation unit, and each channel estimation unit includes At least one pilot subcarrier and one data subcarrier. In an embodiment in which the channel estimation unit partitioning is performed, when the channel coefficient estimation value of a certain data subcarrier position is calculated according to the formula (7), the pilot subcarrier positions corresponding to the data stream in the same channel estimation unit are The channel coefficient estimates give the same weight. In another embodiment in which the channel estimation unit partitioning is performed, when calculating the channel coefficient estimation values of the respective data subcarrier positions corresponding to the data stream in the same channel estimation unit according to formula (7), a set of the same weight is taken. y ., = ν··, /, 7 = 1,···, J, The channel coefficient estimates of the respective data subcarrier positions corresponding to the obtained data stream are the same. In still another embodiment of performing channel estimation unit division, the manner of the above two embodiments may be combined. as follows:
定义第 k个信道估计单元包含的导频子载波的索引构成的集合为 A, k = \,2,---,K; % k个信道估计单元中该数据流对应的每一数据子载波位置的信道系数 估计值相等, 记为 , 接收端按下式来计算该 :  Defining a set of index constituents of the pilot subcarriers included in the kth channel estimation unit is A, k = \, 2, ---, K; % k channel estimation units corresponding to each data subcarrier of the data stream The estimated channel coefficient of the position is equal, and it is recorded that the receiving end calculates the following:
其中, /为一循环变量, 1 = \,2, ···, , 0¾为计算 时, 赋予第 /个信道估 计单元中该数据流对应的各导频子载波位置的信道系数估计值的权值, 因为 Where / is a loop variable, 1 = \, 2, ···, , 03⁄4 is the weight of the channel coefficient estimate assigned to each pilot subcarrier position corresponding to the data stream in the channel estimation unit when calculating Value because
K  K
是加权平均, 要满足条件∑| | =1,0≤ ≤1, 其中 | |表示导频索引 Is a weighted average, to satisfy the condition ∑| | =1,0≤ ≤1, where | | denotes the pilot index
1=1  1=1
集合 包含的导频子载波的个数。 在时频上, 与某个数据子载波位置越近的 导频子载波, 信道相关性就越强。 因此较佳地, 在计算 釆用的权值" w中,The number of pilot subcarriers included in the set. At the time-frequency, the closer the pilot subcarriers are to a certain data subcarrier, the stronger the channel correlation. Therefore, preferably, in calculating the weight used, " w ,
"tt大于等于其他的权值, ! = 、K。 可以看出, 本实施例在按公式( 7 )计算某个数据子载波位置的信道系数 估计值时, 对于同一信道估计单元中各导频子载波位置的信道系数估计值, 取相同的权值, 且计算同一信道估计单元中各数据子载波位置的信道系数估 计值时, 通过取相同的一套权值, 使得得到的各数据子载波位置的信道系数 估计值相同。 时频区域内, 与某个数据子载波位置越近的导频子载波, 信道相关性就 越强。 因此较佳地, 在计算 釆用的权值 ,中, ^½大于等于其他的权值, ί = 1,2,"',Κ。 釆用上述基于信道估计单元的方式可以简化计算。 上述宽带同频干扰噪声估计和干扰抑制方法中,步骤 20的加权平均可以 基于干扰噪声估计单元来进行。接收端将干扰抑制区域再划分为 个时频二 维资源块, Λ/=1,2,...; 每个时频二维资源块作为一个干扰噪声估计单元, 每 一干扰噪声估计单元中包括至少一个导频子载波。 同一干扰抑制区域中信道 估计单元和干扰噪声估计单元的划分可以相同, 也可以不同。 在进行干扰噪声估计单元划分的一实施例中, 按公式 ( 1 )计算某个数据 子载波位置的干扰噪声协方差矩阵估计值时, 为同一干扰噪声估计单元中各 个导频子载波位置的干扰噪声协方差矩阵赋予的权值相同。 在进行干扰噪声估计单元划分的另一实施例, 按公式 ( 1 )计算同一干扰 噪声估计单元中各个数据子载波位置的干扰噪声协方差矩阵估计值时, 取相 同的一组权值 7 = 1,···, J, 得到相同的干扰噪声协方差矩阵估计值。 在进行干扰噪声估计单元划分的又一实施例, 可以结合上述两个实施例 的方式。 ^下: 定义第 m 个干扰噪声估计单元包含的导频子载波的索引构成的集合为 , m = \,2 ,M。 第 w个干扰噪声估计单元中该数据流对应的每一数据子 载波位置的干扰噪声协方差矩阵估计值相等, 记为 M mD , 接收端按下式来计 算:
Figure imgf000020_0001
"tt is greater than or equal to other weights, ! = , K. It can be seen that in this embodiment, when calculating the channel coefficient estimation value of a certain data subcarrier position according to formula (7), for each pilot in the same channel estimation unit Estimating the channel coefficient of the subcarrier position, taking the same weight, and calculating the channel coefficient estimate of each data subcarrier position in the same channel estimation unit When calculating, the same set of weights is taken, so that the obtained channel coefficient estimates of the respective data subcarrier positions are the same. In the time-frequency region, the closer the pilot subcarriers are to a certain data subcarrier, the stronger the channel correlation. Therefore, preferably, in calculating the weight of the weight, ^1⁄2 is greater than or equal to the other weights, ί = 1, 2, "', Κ. 釆 The above method based on the channel estimation unit can simplify the calculation. In the same-frequency interference noise estimation and interference suppression method, the weighted average of step 20 may be performed based on the interference noise estimation unit. The receiving end further subdivides the interference suppression region into a time-frequency two-dimensional resource block, Λ/=1, 2,. Each time-frequency two-dimensional resource block is used as an interference noise estimation unit, and each interference noise estimation unit includes at least one pilot sub-carrier. The division of the channel estimation unit and the interference noise estimation unit in the same interference suppression region may be the same. In an embodiment in which interference noise estimation unit division is performed, when the interference noise covariance matrix estimation value of a certain data subcarrier position is calculated according to formula (1), each pilot in the same interference noise estimation unit is used. The interference noise covariance matrix of the subcarrier position is given the same weight. In another embodiment of performing interference noise estimation unit partitioning, calculation is performed according to formula (1) When an interference noise covariance matrix estimation value of each data subcarrier position in the interference noise estimation unit is obtained, the same set of weights 7 = 1, . . . , J is obtained, and the same interference noise covariance matrix estimation value is obtained. A further embodiment of performing interference noise estimation unit partitioning may be combined with the manner of the above two embodiments. ^Bottom: defining a set of index constituents of pilot subcarriers included in the mth interference noise estimating unit is m = \, 2, M. The estimated noise interference covariance matrix of each data subcarrier position corresponding to the data stream in the w interference noise estimation unit is equal, denoted as M mD , and the receiving end calculates according to the following formula:
Figure imgf000020_0001
其中, /为一循环变量, / = 1,2,···,Μ ; 为计算 ^— D时, 赋予第 /个干 扰噪声估计单元中各导频子载波位置的 ()与 {if乘积的权值, 因为是加 权平均, ,要满足条件∑|4| = 1,0≤ ≤1 , 其中 | |表示导频索引集 Where / is a loop variable, / = 1,2,···,Μ; for calculating ^ -D , assigning the () and {if product of each pilot subcarrier position in the first/interference noise estimation unit Weight, because it is plus The weight average, , must satisfy the condition ∑|4| = 1 , 0 ≤ ≤1 , where | | denotes the pilot index set
1=1  1=1
合 包含的导频子载波的个数。 可以看出, 本实施例在按公式( 1 )计算某个数据子载波位置的干扰噪声 协方差矩阵估计值时, 对于同一干扰噪声估计单元中各导频子载波位置的干 扰噪声协方差矩阵, 取相同的权值; 且在计算同一干扰噪声估计单元中各数 据子载波位置的干扰噪声协方差矩阵估计值时, 通过取相同的一套权值, 使 得各数据子载波位置的干扰噪声协方差矩阵估计值相同。 时频区域内, 与某个数据子载波位置越近的导频子载波, 信道相关性越 强。 因此较佳地, 计算^ _D釆用的权值 中, / = 1, 2,···,Μ , Α ^大于等于 其他的权值。 釆用上述基于干扰噪声估计单元的方式可以简化计算。 相应地, 本实施例还提供了一种宽带同频干扰噪声估计的系统, 应用于 正交频分复用(OFDM)或正交频分多址 (OFDMA)系统的接收端, 在一干扰抑 制区域内对其中承载的一个数据流进行干扰噪声估计, 该干扰抑制区域为接 收数据承载区域中的一时频二维资源块, 该系统包括: 第一装置, 其设置为: 对该数据流对应的每一数据子载波, 将该数据流 对应的各导频子载波上的接收信号与该接收信号的共轭转置的乘积的加权平 均, 作为该数据子载波位置的干扰噪声协方差矩阵估计值; 第二装置, 其设置为: 对该数据流对应的每一数据子载波, 将对第一装 置得到的该数据子载波位置的干扰噪声协方差矩阵估计值进行对角加载后得 到的结果, 作为该数据子载波位置的干扰噪声协方差矩阵。 上述第一装置是设置为按式 ( a )计算出该数据子载波位置的干扰噪声协 方差矩阵估计值:
Figure imgf000021_0001
The number of pilot subcarriers included. It can be seen that, in this embodiment, when calculating the interference noise covariance matrix estimation value of a data subcarrier position according to formula (1), the interference noise covariance matrix of each pilot subcarrier position in the same interference noise estimation unit, Taking the same weight; and calculating the interference noise covariance matrix estimation value of each data subcarrier position in the same interference noise estimation unit, by taking the same set of weights, making the interference noise covariance of each data subcarrier position The matrix estimates are the same. In the time-frequency region, the closer the pilot subcarriers are to a certain data subcarrier, the stronger the channel correlation. Therefore, preferably, in calculating the weight of ^ _ D釆, / = 1, 2,····, Μ, Α ^ is greater than or equal to other weights. The above calculation based on the interference noise estimation unit can simplify the calculation. Correspondingly, the embodiment further provides a system for wideband co-channel interference noise estimation, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression manner. Interference noise estimation is performed on a data stream carried in the area, where the interference suppression area is a time-frequency two-dimensional resource block in the received data bearer area, the system includes: a first device, configured to: correspond to the data stream a weighted average of the product of the received signal on each pilot subcarrier corresponding to the data stream and the conjugate transposed of the received signal, as the estimated value of the interference noise covariance matrix of the data subcarrier position And a second device, configured to: perform a diagonal loading on the estimated value of the interference noise covariance matrix of the data subcarrier position obtained by the first device for each data subcarrier corresponding to the data stream, As the interference noise covariance matrix of the data subcarrier position. The first device is configured to calculate an interference noise covariance matrix estimate of the data subcarrier position according to equation (a):
Figure imgf000021_0001
其中, 为该干扰抑制区域中该数据流对应的第 j个数据子载波位 置的干扰噪声协方差矩阵估计值, · = 1,··· , J为该干扰抑制区域中该数据流 对应的数据子载波个数; = 1, · · ·,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; ^为计算该干扰抑制区域中该数据流对应的第 _;个数据子载 波位置的干扰噪声协方差矩阵估计值时赋予 ( (i)) (yp (l)f 的权值,Wherein, the estimated value of the interference noise covariance matrix of the jth data subcarrier position corresponding to the data stream in the interference suppression region, · = 1,···, J is the data flow in the interference suppression region Corresponding number of data subcarriers; = 1, · · ·, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region; ^ is to calculate the corresponding data stream in the interference suppression region _; the interference noise covariance matrix estimate of the data subcarrier position is given ((i)) the weight of (y p (l)f,
∑βϋ = \ ·, W为所述接收端在该干扰抑制区域中该数据流对应的第 ,个导 频子载波上的接收信号, (广为 ( 的共轭转置矩阵。 该宽带同频干扰噪声估计的系统还可包括第三装置, 其中: 所述第三装置设置为: 将该干扰抑制区域划分为一个或多个干扰噪声估 计单元, 每一干扰噪声估计单元为一个时域二维资源块且其中包含至少一个 导频子载波和一个数据子载波; 相应地, 所述第一装置按式 (a)计算数据子载波位置的干扰噪声协方差矩 阵估计值时, 为同一干扰噪声估计单元中各个导频子载波上的接收信号与其 共轭转置的乘积, 赋予相同的权值。 该宽带同频干扰噪声估计的系统还可包括第四装置, 其中: 所述第四装置设置为: 将该干扰抑制区域划分为 个干扰噪声估计单 元, 每一干扰噪声估计单元为一个时域二维资源块且其中包含至少一个导频 子载波和一数据子载波, 为正整数; 相应地, 对该数据流对应的每一数据子载波, 所述第一装置是设置为按 如下公式将该数据流对应的各导频子载波上的接收信号与该接收信号的共轭 转置的乘积的加权平均, 作为该数据子载波位置的干扰噪声协方差矩阵估计 值:
Figure imgf000022_0001
其中,
Σβ ϋ = \ ·, W is the receiving end of, the interference suppression of guide area corresponding to the data stream received signal on pilot subcarriers, (widely (the conjugate transpose matrix of the same broadband The system for frequency interference noise estimation may further include a third device, wherein: the third device is configured to: divide the interference suppression region into one or more interference noise estimation units, each interference noise estimation unit is a time domain two a dimension resource block and including at least one pilot subcarrier and one data subcarrier; correspondingly, the first device calculates the interference noise covariance matrix estimation value of the data subcarrier position according to formula (a), and is the same interference noise The product of the received signal on each pilot subcarrier in the estimation unit and its conjugate transpose is given the same weight. The system for wideband co-channel interference noise estimation may further comprise a fourth device, wherein: the fourth device is set The interference suppression area is divided into interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least one pilot. And a data subcarrier, which is a positive integer; correspondingly, for each data subcarrier corresponding to the data stream, the first device is set to be on each pilot subcarrier corresponding to the data stream according to the following formula A weighted average of the product of the received signal and the conjugate transpose of the received signal as an estimate of the interference noise covariance matrix for the data subcarrier position:
Figure imgf000022_0001
among them,
O 第 m 个干扰噪声估计单元中该数据流对应的每一数据子载波位 置的干扰噪声协方差矩阵估计值, m = 1, 2, · · · ,Μ; /为一循环变量, / = 1, 2, · · ·, ; 为第 /个干扰噪声估计单元包含的该数据流对应的导频子载波的索引 Ζ· 的集合, = 1,···,/, /为该干扰抑制区域中该数据流对应的导频子载波的个 数;O The estimated value of the interference noise covariance matrix for each data subcarrier position corresponding to the data stream in the mth interference noise estimation unit, m = 1, 2, · · · , Μ; / is a cyclic variable, / = 1 , 2, · · ·, ; a set of indices Ζ· of pilot subcarriers corresponding to the data stream included in the first interference noise estimation unit, = 1,···, /, / is the pilot corresponding to the data stream in the interference suppression region The number of carriers;
^;)为所述接收端在该数据流对应的第 I个导频子载波上的接收信号; yp (if为 ()的共轭转置矩阵; fiml为计算 —D时,赋予第 /个干扰噪声估计单元中该数据流对应的各导 频子载波对应的 ( 与 ( 的乘积的权值,^;) is the received signal of the receiving end on the first pilot subcarrier corresponding to the data stream; y p (if is the conjugate transposed matrix of (); fi ml is calculated - D , the first / interference noise estimation unit corresponding to each pilot subcarrier corresponding to the data stream (the weight of the product of (,
Figure imgf000023_0001
, | |为 包含的导频子载波的个数, 大于等于其他的权值。 对该数据流对应的每一数据子载波, 所述第二装置是设置为按如下公式 将对该数据子载波位置的干扰噪声协方差矩阵估计值进行对角加载后得到的 结果, 作为该数据子载波位置的干扰噪声协方差矩阵:
Figure imgf000023_0001
, | | is the number of pilot subcarriers included, and is greater than or equal to other weights. For each data subcarrier corresponding to the data stream, the second device is configured to diagonally load the interference noise covariance matrix estimation value of the data subcarrier position as the data, as the data Interference noise covariance matrix for subcarrier position:
NI-D NI-D
其中, RW_D /)为该干扰抑制区域中该数据流对应的第 _;个数据子载波位 置的干扰噪声协方差矩阵, · = 1,··· , J为该干扰抑制区域中该数据流对应的 数据子载波个数; ≥0; 为该干扰抑制区域中该数据流对应的第 _;个 数据子载波位置的干扰噪声协方差矩阵估计值; β>0 Λ表示 N xN¾的对角 矩阵, ¾表示所述接收端的接收天线的数目。 该宽带同频干扰噪声估计的系统中, 所述第二装置计算出的第 个数据 子载波位置的干扰噪声协方差矩阵为: RW_D 0) = RM_D (j) + ; 其中, 0≤ ≤1; tr( MDC/》表示对矩阵 — β( ·)求迹; /为该干扰抑制区 域内该数据流对应的导频子载波的个数, 1为^„^^¾的单位矩阵。 该宽带同频干扰噪声估计的系统还可包括第五装置: 所述第五装置设置为: 按以下方式得到发送端在该数据流对应的每一导 频子载波位置的信道系数估计值并输出到所述第一装置: 将该导频子载波上的接收信号与发送端在该导频子载波上发送的导频信 号的共轭相乘, 得到该导频子载波位置的信道系数估计值。 Where R W _ D /) is the interference noise covariance matrix of the _th data subcarrier position corresponding to the data stream in the interference suppression region, · = 1,···, J is the interference suppression region The number of data subcarriers corresponding to the data stream; ≥ 0; is the estimated value of the interference noise covariance matrix of the _; data subcarrier positions corresponding to the data stream in the interference suppression region; β>0 Λ represents N xN 3⁄4 The diagonal matrix, 3⁄4, represents the number of receiving antennas at the receiving end. In the system for estimating the wideband co-channel interference noise, the interference noise covariance matrix of the first data subcarrier position calculated by the second device is: R W _ D 0) = R M _ D (j) + ; , 0≤ ≤1; tr (M - D C / " represents a matrix - β (·) track seeking; / interference area for the data stream corresponding to the number of pilot subcarriers is suppressed, 1 ^" ^ ^ ¾ of the matrix may be the same wideband frequency interference noise estimation system comprises fifth means: the fifth means is arranged to: obtain the data stream at the transmitting end corresponding to each pilot subcarrier position in the following manner The channel coefficient estimate is output to the first device: Multiplying the received signal on the pilot subcarrier with the conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier to obtain an estimated channel coefficient of the pilot subcarrier position.
相应地, 本实施例还提供了一种宽带同频干扰抑制的系统, 应用于正交 频分复用 (OFDM ) /正交频分多址(OFDMA ) 系统的接收端, 在一干扰抑 制区域内对其中承载的一个数据流进行干扰抑制, 该干扰抑制区域为接收数 据承载区域中的一时频二维资源块, 该系统包括: 第一子系统, 其设置为: 按上述宽带同频干扰噪声估计的系统相同的方 式, 得到该数据流对应的各导频子载波位置的信道系数估计值和各数据子载 波位置的干扰噪声协方差矩阵; 第二子系统, 其设置为: 对该数据流对应的每一数据子载波, 将该数据 流对应的各个导频子载波位置的信道系数估计值的加权平均, 作为该数据子 载波位置的信道系数估计值; 以及 第三子系统, 其设置为: 对该数据流对应的每一数据子载波, 根据该数 据子载波上的接收信号, 及该数据子载波位置的信道系数估计值和干扰噪声 协方差矩阵 , 计算得到该数据子载波上的数据信号估计。 进一步地, 上述宽带同频干扰抑制的系统中还可包括第四子系统: 所述第四子系统设置为: 对该数据流对应的每一数据子载波, 在所述第 三子系统计算得到该数据载波上的数据信号估计后, 对该估计出的数据信号 进行解调, 并将解调后得到的该信号各比特对应的软信息 进行调整, 调整后各比特对应的软信息分别为
Figure imgf000024_0001
; 其中, · = ι, ···, J , J为该干扰抑制区域中该数据流对应的数据子载波个数; F为该干扰抑制区域中该数据流对应的第 j个数据载波上的数据信号中包含 的比特个数, 表示该干扰抑制区域中该数据流对应的第 _/·个数据载波对
Correspondingly, the embodiment further provides a system for wideband co-channel interference suppression, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM)/Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression region. Performing interference suppression on a data stream carried therein, the interference suppression area being a time-frequency two-dimensional resource block in the received data bearer area, the system comprising: a first subsystem, configured to: according to the above-mentioned broadband co-channel interference noise The estimated system obtains the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream and the interference noise covariance matrix of each data subcarrier position; the second subsystem is set to: the data stream Corresponding to each data subcarrier, a weighted average of channel coefficient estimates of respective pilot subcarrier positions corresponding to the data stream is used as an estimated channel coefficient of the data subcarrier position; and a third subsystem is set to : for each data subcarrier corresponding to the data stream, according to the received signal on the data subcarrier, and the channel of the data subcarrier position Estimated value covariance matrix and the interference noise, the data signal estimate is calculated on the data subcarriers. Further, the system for the broadband inter-frequency interference suppression may further include a fourth subsystem: the fourth subsystem is configured to: calculate, for each data subcarrier corresponding to the data stream, the third subsystem After estimating the data signal on the data carrier, the estimated data signal is demodulated, and the soft information corresponding to each bit of the signal obtained after demodulation is adjusted, and the soft information corresponding to each bit after adjustment is respectively
Figure imgf000024_0001
Where = · ι, ···, J, J is the number of data subcarriers corresponding to the data stream in the interference suppression region; F is the jth data carrier corresponding to the data stream in the interference suppression region The number of bits included in the data signal indicates the _/·th data carrier pair corresponding to the data stream in the interference suppression region
FF ( ')(RW-D ( '))— ( ·) 或 者 信 干 噪 比 值
Figure imgf000024_0002
hd {j)为该干扰抑制区域中该数据流对应的第 j个数据子载波位置的信道 系数估计值, /)为 4( ·;)的共轭转置矩阵, RM_D /;>为第 ·个数据子载波位 置的干扰噪声协方差矩阵, (RMD ( ))'1为 RM_D (_/·)的逆矩阵; = 1, ···,/, /为该干扰抑制区域中该数据流对应的导频子载波个数; ft ή 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道系数估计 值, 为 的共轭转置矩阵, M_PW为第 ,个导频子载波位置的干扰 噪声协方差矩阵, 1为^ 的逆矩阵; ()为所述接收端在该干 扰抑制区域中该数据流对应的第 ,个导频子载波上的接收信号, 为发送 端在该干扰抑制区域中该数据流对应的第 I个导频子载波上的发射信号。 上述宽带同频干扰抑制的系统中还可包括第五子系统和第六子系统: 所述第五子系统设置为: 将该干扰抑制区域划分为 K个信道估计单元, 每一信道估计单元为一个时域二维资源块且其中包含至少一个导频子载波和 一个数据子载波, f为正整数; 所述第六子系统设置为: 对该数据流对应的每一数据子载波, 釆用如下 计算公式将该数据流对应的各个导频子载波位置的信道系数估计值的加权平 均, 作为该数据子载波位置的信道系数估计值的步骤:
Figure imgf000025_0001
F F ( ')(R W - D ( '))— ( )) or signal to interference and noise ratio
Figure imgf000024_0002
h d {j) is an estimated channel coefficient of the jth data subcarrier position corresponding to the data stream in the interference suppression region, /) is a conjugate transposed matrix of 4 (·;), R M _ D /; > is the interference noise covariance matrix of the first data subcarrier position, (R MD ( )) ' 1 is the inverse matrix of R M _ D (_/·); = 1, ···, /, / The number of pilot subcarriers corresponding to the data stream in the interference suppression region; ft ή is the channel coefficient estimation value of the first pilot subcarrier position corresponding to the data stream in the interference suppression region, and is a conjugate rotation The matrix, M _ P W is the interference noise covariance matrix of the first pilot subcarrier position, where 1 is the inverse matrix of ^; () is the number corresponding to the data stream in the interference suppression region of the receiving end, The received signal on the pilot subcarrier is the transmit signal on the first pilot subcarrier corresponding to the data stream in the interference suppression region. The system for wideband co-channel interference suppression may further include a fifth subsystem and a sixth subsystem: the fifth subsystem is configured to: divide the interference suppression region into K channel estimation units, where each channel estimation unit is a time domain two-dimensional resource block and including at least one pilot subcarrier and one data subcarrier, where f is a positive integer; the sixth subsystem is set to: for each data subcarrier corresponding to the data stream, The weighting average of the channel coefficient estimation values of the respective pilot subcarrier positions corresponding to the data stream is used as a step of estimating the channel coefficient of the data subcarrier position as follows:
Figure imgf000025_0001
/=1 ieQf 其中,  /=1 ieQf where,
为第 个信道估计单元中该数据流对应的每一数据子载波位置的信道 系数估计值, k = \,2,---,K  The channel coefficient estimate for each data subcarrier location corresponding to the data stream in the first channel estimation unit, k = \, 2, ---, K
/为一循环变量, ! = 1,2,···,Κ / is a loop variable, ! = 1,2,···,Κ
Ω,为第 个信道估计单元包含的该干扰抑制区域中该数据流对应的导频 子载波的索引 ,的集合, = 1,···,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道 系数估计值; 以及 akl为计算 时,赋予第 /个信道估计单元中各导频子载波位置的 ( )的Ω, which is a set of indices of pilot subcarriers corresponding to the data stream in the interference suppression region included in the first channel estimation unit, = 1, . . . , /, / is corresponding to the data stream in the interference suppression region The number of pilot subcarriers; the channel of the first pilot subcarrier position corresponding to the data stream in the interference suppression region Coefficient estimated value; and a kl is the ( ) given to the position of each pilot subcarrier in the channel estimation unit when calculating
Κ Κ
权值, ∑|Ω= 1,0≤ ≤1 , | |表示 包含的导频子载波的个数, 且在权 Weight, ∑| Ω= 1 , 0 ≤ ≤1 , | | indicates the number of pilot subcarriers included, and is in the right
1=1  1=1
值 0¾中, l = \, 2, - - -,K , 大于等于其他的权值 ( In the value 03⁄4, l = \, 2, - - -, K, is greater than or equal to other weights (
下面用一些应用示例对本发明所述方法进行进一步说明,在以下示例中, 各参数的含义与上述实施例方案相同, 且假定接收端已经获得每个数据子载 波上的接收信号。 示例中主要说明在不同的干扰抑制区域样式和干扰噪声估 计单元划分的情况下, 如何进一步计算得到数据子载波位置的干扰噪声协方 差矩阵, 对于数据信号估计见上文, 也不再重复。 应用示例一 The method of the present invention will be further described below with some application examples. In the following examples, the meaning of each parameter is the same as that of the above embodiment, and it is assumed that the receiving end has obtained the received signal on each data subcarrier. The example mainly shows how to further calculate the interference noise covariance matrix of the data subcarrier position in the case of different interference suppression region patterns and interference noise estimation unit partitioning. For the data signal estimation, see above, it will not be repeated. Application example one
上包含 15个连续的 OFDM/OFDMA符号, 在频域上包含 4个连续子载波, 其中承载一个数据流。 在本示例中, 该干扰抑制区域被划分为 1个干扰噪声估计单元, 该干扰 抑制区域中包含的 20个导频子载波的索引属于 1个导频索引集合, 1~20均 属于 。 在进行干扰噪声估计时, 该干扰噪声估计单元内的各个数据子载波位置 的干扰噪声协方差矩阵估计值均为 ^ , 有: It contains 15 consecutive OFDM/OFDMA symbols, and contains 4 consecutive subcarriers in the frequency domain, which carries one data stream. In this example, the interference suppression region is divided into one interference noise estimation unit, and the index of the 20 pilot subcarriers included in the interference suppression region belongs to one pilot index set, and 1 to 20 belong to . When performing interference noise estimation, the estimated value of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit is ^, with:
^ 20 H ^ 20 H
'=ι  '=ι
其中, A满足条件 A =丄。  Where A satisfies the condition A = 丄.
20  20
应用示例二  Application example two
抑制区域划分成 2个干扰噪声估计单元,该干扰抑制区域中包含的 20个导频 子载波分别属于 2个导频索引集合, 其中 1~8属于 Ω, , 9~20属于 Ω2。 在进行干扰噪声估计时: 第一个干扰噪声估计单元内的各个数据子载波位置的干扰噪声协方差矩 阵估计值均为 ^ , 有: The suppression region is divided into two interference noise estimation units, and the 20 pilot subcarriers included in the interference suppression region belong to two pilot index sets, wherein 1 to 8 belong to Ω, and 9 to 20 belong to Ω 2 . When performing interference noise estimation: The estimated value of the interference noise covariance matrix of each data subcarrier position in the first interference noise estimation unit is ^, with:
% H 20 H % H 20 H
i=\ i=9 第二个干扰噪声估计单元内的各个数据子载波位置的干扰噪声协方差矩 阵估计值均为 , 有:  i=\ i=9 The estimated values of the interference noise covariance matrix for each data subcarrier position in the second interference noise estimation unit are:
其中, k = i, 2 , | |表示导频索引集合
Figure imgf000027_0001
Where k = i, 2 , | | represent the set of pilot indices
Figure imgf000027_0001
中包含的导频子载波的个数, βα > ^~ The number of pilot subcarriers included, β α > ^~
1=1 1=1
应用示例三  Application example three
抑制区域划分成 2个干扰噪声估计单元,该干扰抑制区域中包含的 20个导频 子载波分别属于 2个导频索引集合, 其中: 1、 2、 5、 6、 9、 10、 13、 14、 17、 18属于 Ω, , 其余导频索引属于 Ω2。 在进行干扰噪声估计时: 第一个干扰噪声估计单元内的各个导频子载波位置的干扰噪声协方差矩 阵均为 ^ , 有:
Figure imgf000027_0002
第二个干扰噪声估计单元内的各个导频子载波位置的干扰噪声协方差矩 阵均为 , 有:
Figure imgf000027_0003
则该干扰抑制区域内第 j个数据子载波位置的干扰噪声协方差矩阵估计 值为:
The suppression region is divided into two interference noise estimation units, and the 20 pilot subcarriers included in the interference suppression region belong to two pilot index sets, where: 1, 2, 5, 6, 9, 10, 13, 14 17, 18, 18 belong to Ω, and the rest of the pilot index belongs to Ω 2 . When performing interference noise estimation: the interference noise covariance matrix of each pilot subcarrier position in the first interference noise estimation unit is ^, with:
Figure imgf000027_0002
The interference noise covariance matrix of each pilot subcarrier position in the second interference noise estimation unit is:
Figure imgf000027_0003
Then, the interference noise covariance matrix estimation of the jth data subcarrier position in the interference suppression region The value is:
其 , 0≤βα≤1 , k = l, 2 , | |表示导频索引集合
Figure imgf000028_0001
It, 0 ≤ β α ≤ 1, k = l, 2 , | | represents the pilot index set
Figure imgf000028_0001
中包含的导频子载波的个数, β > ^~ , η + η2! = 1。 The number of pilot subcarriers included, β > ^~ , η + η 2! = 1.
1=1  1=1
应用示例四 如图 6所示, 在本示例中的干扰抑制区域为干扰抑制区域样式二。 在时 域上包含 12个连续的 OFDM/OFDMA符号,在频域上包含 4个连续子载波, 其中承载一个数据流。 在本示例中, 该干扰抑制区域被划分为 1个干扰噪声估计单元, 该干扰 抑制区域中包含的 16个导频子载波的索引属于 1个导频索引集合, 1~16均 属于 。 在进行干扰噪声估计时, 该干扰噪声估计单元内的各个数据子载波位置 的干扰噪声协方差矩阵估计值均为 ^ , 有:
Figure imgf000028_0002
其中, A满足条件 A =^。
Application Example 4 As shown in FIG. 6, the interference suppression area in this example is the interference suppression area pattern 2. It contains 12 consecutive OFDM/OFDMA symbols in the time domain and 4 consecutive subcarriers in the frequency domain, which carries one data stream. In this example, the interference suppression region is divided into one interference noise estimation unit, and the index of the 16 pilot subcarriers included in the interference suppression region belongs to one pilot index set, and all of 1 to 16 belong. When performing interference noise estimation, the estimated values of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit are ^, with:
Figure imgf000028_0002
Where A satisfies the condition A = ^.
16  16
应用示例五 如图 7所示, 在本示例中的干扰抑制区域为干扰抑制区域样式三。 在时 域上包含 9个连续的 OFDM/OFDMA符号, 在频域上包含 4个连续子载波, 其中承载一个数据流。 在本示例中, 该干扰抑制区域被划分为 1个干扰噪声估计单元, 该干扰 抑制区域中包含的 12个导频子载波的索引属于 1个导频索引集合, 1~12均 属于 。 在进行干扰噪声估计时, 该干扰噪声估计单元内的各个数据子载波位置 的干扰噪声协方差矩阵估计值均为 ^ , 有: — = ∑ (Φ ' 其中, Α满足条件 = 。 下面用应用示例六和七说明一个干扰抑制区域内承载两个数据流时的干 扰噪声估计方法。 假设: Application Example 5 As shown in FIG. 7, the interference suppression area in this example is the interference suppression area pattern 3. It contains 9 consecutive OFDM/OFDMA symbols in the time domain and 4 consecutive subcarriers in the frequency domain, which carries one data stream. In this example, the interference suppression region is divided into one interference noise estimation unit, and the index of the 12 pilot subcarriers included in the interference suppression region belongs to one pilot index set, and 1 to 12 belong to each other. When performing interference noise estimation, the estimated values of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit are ^, with: — = ∑ (Φ ' where Α satisfies the condition = . The following uses the application examples 6 and 7 to illustrate the interference noise estimation method for carrying two data streams in an interference suppression region.
第一个数据流对应的第 个导频子载波上的信道系数估计值为 fip ( ); 第 二个数据流对应的第 个导频子载波上的信道系数估计值为 fip2 ^>; 发送端在 第一个数据流对应的第 个导频子载波上发送的导频信号为 发送端在 第二个数据流对应的第 个导频子载波上发送的导频信号为 ρ^ή 接收端在 第一个数据流对应的第 ,个导频子载波上接收的接收信号为 以及接收 端在第二个数据流对应的第 I个导频子载波上接收的接收信号为 2 ( )。 应用示例六 如图 8所示, 本应用示例中的干扰抑制区域为干扰抑制区域样式四, 在 时域上包含 15个连续的 OFDM/OFDMA符号, 在频域上包含 4个连续子载 波, 其中承载两个数据流。 在本示例中, 将干扰抑制区域分成一个干扰噪声估计单元, 该干扰抑制 区域中包含的每个数据流对应的 10个导频子载波的索引分别属于 2个导频索 引集合, 其中: 1~10属于 Ω, , 11~20属于 ί¾。 在进行干扰噪声估计时: 对于第一个数据流, 该干扰噪声估计单元内的各个数据子载波位置的干 扰噪声协方差矩阵估计值均为 ά11 , 有:
Figure imgf000029_0001
The channel coefficient estimate on the first pilot subcarrier corresponding to the first data stream is fi p ( ); the channel coefficient estimate on the first pilot subcarrier corresponding to the second data stream is fi p2 ^>; The pilot signal sent by the transmitting end on the first pilot subcarrier corresponding to the first data stream is the pilot signal sent by the transmitting end on the first pilot subcarrier corresponding to the second data stream is ρ^ή receiving The received signal received by the terminal on the first pilot subcarrier corresponding to the first data stream and the received signal received by the receiving end on the first pilot subcarrier corresponding to the second data stream are 2 ( ). Application Example 6 As shown in FIG. 8, the interference suppression region in this application example is an interference suppression region pattern 4, which includes 15 consecutive OFDM/OFDMA symbols in the time domain and 4 consecutive subcarriers in the frequency domain, where Carry two data streams. In this example, the interference suppression region is divided into an interference noise estimation unit, and the indexes of the 10 pilot subcarriers corresponding to each data stream included in the interference suppression region belong to two pilot index sets, where: 1~ 10 belongs to Ω, and 11~20 belongs to ί3⁄4. When performing interference noise estimation: For the first data stream, the estimated value of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit is ά 11 , having:
Figure imgf000029_0001
对于第二个数据流, 该干扰噪声估计单元内的各个数据子载波位置的干 扰噪声协方差矩阵估计值均为 ά12 , 有:
Figure imgf000029_0002
其中, =A = 。 上述对每个数据流进行干扰噪声估计的方法^^本相同的。 当然在另一 实施例中, 针对不同数据流的干扰噪声估计单元的划分也是可以不同的。 应用示例七 如图 9所示, 本应用示例中的干扰抑制区域为干扰抑制区域样式五, 在 时域上包含 6个连续的 OFDM/OFDMA符号,在频域上包含 6个连续子载波, 其中承载两个数据流。 在本示例中, 将干扰抑制区域分成 1个干扰噪声估计单元, 该干扰抑制 区域中包含的每个数据流对应的 4个导频子载波索引均属于 1个导频索引集 合 Ω, , 即 1~4属于 Ω,。 在进行干扰噪声估计时: 对于第一个数据流, 干扰噪声估计单元内的各个数据子载波位置的干扰 噪声协方差矩阵估计值均为 ά11 , 有:
Figure imgf000030_0001
对于第二个数据流, 干扰噪声估计单元内的各个数据子载波位置的干扰 噪声协方差矩阵估计值均为 ά12 , 有:
For the second data stream, the estimated value of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit is ά 12 , having:
Figure imgf000029_0002
Where =A = . The above method for performing interference noise estimation for each data stream is the same. Of course, in another embodiment, the division of the interference noise estimation unit for different data streams may also be different. Application Example 7 As shown in FIG. 9, the interference suppression region in this application example is an interference suppression region pattern five, which includes 6 consecutive OFDM/OFDMA symbols in the time domain and 6 consecutive subcarriers in the frequency domain, where Carry two data streams. In this example, the interference suppression region is divided into one interference noise estimation unit, and the four pilot subcarrier indexes corresponding to each data stream included in the interference suppression region belong to one pilot index set Ω, that is, 1 ~4 belongs to Ω,. When performing interference noise estimation: For the first data stream, the estimated value of the interference noise covariance matrix of each data subcarrier position in the interference noise estimation unit is ά 11 , with:
Figure imgf000030_0001
For the second data stream, the estimated values of the interference noise covariance matrix for each data subcarrier position within the interference noise estimation unit are ά 12 , with:
其中, A
Figure imgf000030_0002
以下应用示例主要说明对上述估计出的干扰噪声估计协方差矩阵估计值 进行对角加载。 需要说明的是, 下文以在一个干扰抑制区域内仅承载一个数 据流的情形进行描述, 对于同时承载多个数据流的情形, 对每个数据流, 利 用与下述应用示例相同的方法进行各数据流对应的干扰噪声估计协方差矩阵 估计值分别进行对角加载即可。 应用示例八
Among them, A
Figure imgf000030_0002
The following application examples mainly illustrate the diagonal loading of the estimated interference noise estimation covariance matrix estimates above. It should be noted that the following describes a case where only one data stream is carried in one interference suppression area. For the case where multiple data streams are simultaneously carried, for each data stream, each method is performed by the same method as the following application example. The estimated values of the interference noise estimation covariance matrix corresponding to the data stream are respectively diagonally loaded. Application example eight
在本示例中, 假设接收端的接收天 N = 8 第 个数据子载波位置的 干扰噪声协方差矩阵估计值为 UJ、 对于任意大小的干
Figure imgf000031_0001
In this example, assume that the receiver's receive day N = 8 The first data subcarrier position's interference noise covariance matrix is estimated to be UJ, for any size of the stem
Figure imgf000031_0001
扰抑制区域, 对第 _/·个数据子载波位置的干扰噪声协方差矩阵估计值完成对 角加载后得到该数据子载波位置的干扰噪声协方差矩阵 RMD /;), 有: The interference suppression region, the interference noise covariance matrix estimated value of the _/· data subcarrier positions is obtained after the diagonal loading, and the interference noise covariance matrix R MD /;) of the data subcarrier position is obtained,
Figure imgf000031_0002
Figure imgf000031_0002
其中, /表示当前数据流在该干扰抑制区域内包含的导频子载波数目。 γ 为一个常数, 满足条件 0≤ ≤1。 应用示例九 在本示例中, 假设接收端的接收 N = 8, 第 _/·个数据子载波位置的 干扰噪声协方差矩阵估计值为 NI_。ij) 在当前干扰抑制区
Figure imgf000031_0003
Where / represents the number of pilot subcarriers that the current data stream contains in the interference suppression region. γ is a constant and satisfies the condition 0 ≤ ≤1. Application Example 9 In this example, it is assumed that the receiving end receives N = 8, and the estimated noise covariance matrix of the _/· data subcarrier positions is estimated to be NI _. Ij) in the current interference suppression zone
Figure imgf000031_0003
域内, 对第 _/·个数据子载波位置的干扰噪声协方差矩阵估计值完成对角加载 后得到该数据子载波位置的干扰噪声协方差矩阵 RM_DC/;>, 有: In the domain, the interference noise covariance matrix estimated value of the _/· data subcarrier positions is subjected to diagonal loading to obtain the interference noise covariance matrix R M _ D C/;> of the data subcarrier position, and there are:
Figure imgf000031_0004
其中, /表示当前数据流在该干扰抑制区域内包含的导频子载波数目, 在如图 3所示干扰抑制区中, / = 20; 在如图 6所示干扰抑制区中, / = 16。 γ 为一个常数,满足条件 0≤^≤1,且其取值与 f有关,例如,当^:≤12时, =0.01; 当 f>12时, =0。 以下应用示例主要说明接收端在解调之后进行的软信息调整步骤。 应用示例十 在接收端完成对每个数据子载波上的发射信号解调之后, 假设发送端在 第 个数据子载波上发送的信号中包含 F个软信息, 分别为 A,...JF , 则可以 利用信道增益对软信息进行调整, 调整后的各比特对应的软信息度量分别为 Gm " ,其中 G(_7)的计算方法为 C/)(RM—。(_/))— 当然 软 信 息 进 行 调 整 , 即 :
Figure imgf000032_0001
Figure imgf000031_0004
Where / represents the number of pilot subcarriers included in the interference suppression region of the current data stream, in the interference suppression region shown in Figure 3, / = 20; in the interference suppression region shown in Figure 6, / = 16 . γ is a constant satisfying the condition 0 ≤ ^ ≤ 1, and its value is related to f, for example, when ^: ≤ 12, = 0.01; when f > 12, =0. The following application examples mainly describe the soft information adjustment steps performed by the receiving end after demodulation. Application Example 10 After the receiving end demodulates the transmitted signal on each data subcarrier, it is assumed that the transmitting end transmits F soft information in the signal transmitted on the first data subcarrier, respectively, A, ... J F , the soft information can be adjusted by using the channel gain, and the soft information metric corresponding to each bit after adjustment is Gm " , wherein G (_7) is calculated as C /) (R M - (_ /)) - Of course, the soft information is adjusted, namely:
Figure imgf000032_0001
以下应用示例主要说明利用本发明所述方法进行宽带同频干扰噪声估计 和干扰抑制。 需要说明的是, 下文以在一个干扰抑制区域内仅承载一个数据 流的情形进行描述, 对于同时承载多个数据流的情形, 对每个数据流, 利用 与下述应用示例相同的方法进行干扰抑制信号的接收即可。 应用示例十一 在本应用示例中, 以干扰抑制区域样式一为例进行说明, 在本示例中, 信道估计单元的划分方法和干扰噪声估计单元的划分方法相同, 即各信道估 计单元中包含的子载波与各干扰噪声估计单元中包含的子载波对应相同。 The following application examples primarily illustrate broadband inter-frequency interference noise estimation and interference suppression using the method of the present invention. It should be noted that the following describes the case where only one data stream is carried in one interference suppression area. For the case of carrying multiple data streams simultaneously, for each data stream, interference is performed by the same method as the following application example. Suppress the reception of the signal. Application Example 11 In this application example, the interference suppression region pattern 1 is taken as an example. In this example, the channel estimation unit division method and the interference noise estimation unit are divided in the same manner, that is, included in each channel estimation unit. The subcarriers are identical to the subcarriers included in each interference noise estimation unit.
首先, 接收端将干扰抑制区域内每个导频子载波上的接收信号 与发 送端在导频子载波上发送的导频信号的共轭对应相乘, 即: W = ^ »;) , 其中 表示发送端在导频子载波 I上发送的导频信号 的共轭; 然后, 将样式一的干扰抑制区域分成 1份, 如图 3所示, 每一份为一个 信道估计单元, 亦为一个干扰噪声估计单元, 该干扰抑制区域中包含的 20 个导频子载波属于 1个导频索引集合, 即: 1~20属于 在进行信道估计的时候:  First, the receiving end multiplies the received signal on each pilot subcarrier in the interference suppression region by the conjugate corresponding to the pilot signal transmitted by the transmitting end on the pilot subcarrier, that is, W = ^ »;) Representing the conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier I; then, dividing the interference suppression area of the pattern one into one, as shown in FIG. 3, each one is a channel estimating unit, which is also a The interference noise estimation unit, the 20 pilot subcarriers included in the interference suppression region belong to one pilot index set, that is, 1 to 20 belong to when performing channel estimation:
该信道估计单元内的所有数据子载波位置的信道系数估计值均为 , 有:
Figure imgf000033_0001
在进行干扰噪声估计时: 该干扰噪声估计单元内的所有数据子载波位置的干扰噪声协方差矩阵估 计值均为 ^ , 有:
The channel coefficient estimates for all data subcarrier locations within the channel estimation unit are:
Figure imgf000033_0001
When performing interference noise estimation: the estimated value of the interference noise covariance matrix of all data subcarrier positions in the interference noise estimation unit is ^, with:
1 20 H 1 20 H
RL=^∑^ (0(^(0) 在完成数据子载波 DsC(/)对应的信道估计和干扰噪声协方差矩阵估计值 R L=^∑^ (0(^(0) Completion of channel estimation and interference noise covariance matrix estimation corresponding to data subcarrier DsC(/)
D ;>估计之后,接着完成对该干扰噪声协方差矩阵估计值的对角加载,得 到干扰噪声协方差矩阵 RM_DC/;>, 有: D;> After estimation, then complete the loading of interference noise covariance diagonal covariance matrix estimates to obtain the interference noise covariance matrix R M _ D C /;> , are:
Figure imgf000033_0002
然后接收端进行数据检测, 有:
Figure imgf000033_0002
Then the receiver performs data detection, which has:
其中的 4( )为列向量; 当 4( ·)为行向量时, 相应的数据检测公式为:
Figure imgf000033_0003
Where 4() is the column vector; when 4(·) is the row vector, the corresponding data detection formula is:
Figure imgf000033_0003
其中 co«»表示对输入向量或者标量的每个元素求取共轭。 然后进行软信息调整(此步骤为可选, 即可以忽略)  Where co«» denotes the conjugate of each element of the input vector or scalar. Then perform soft information adjustment (this step is optional, ie can be ignored)
将检测得到的数据 ^ 送入解调器完成解调后, 得到每个数据符号 s 包含的每个比特对应的软信息: , 可选的, 对该软信息进行相应的调 整, 调整后的各比特对应的软信息度量分别为(?^/;^,… ^/;)^。 其中,
Figure imgf000033_0004
当然上述的软信息调整过程也可以忽略, 即直接将解调之后的信息送入 到译码模块完成整个数据接收过程。
After the detected data is sent to the demodulator to complete the demodulation, the soft information corresponding to each bit included in each data symbol s is obtained: , optionally, the soft information is adjusted accordingly, and the adjusted The soft information metrics corresponding to the bits are (?^/;^,... ^/;)^. among them,
Figure imgf000033_0004
Of course, the above soft information adjustment process can also be neglected, that is, the information after demodulation is directly sent. Go to the decoding module to complete the entire data receiving process.
本领域普通技术人员可以理解上述方法中的全部或部分步骤可通过程序 来指令相关硬件完成, 所述程序可以存储于计算机可读存储介质中, 如只读 存储器、 磁盘或光盘等。 可选地, 上述实施例的全部或部分步骤也可以使用 一个或多个集成电路来实现。 相应地, 上述实施例中的各模块 /单元可以釆用 硬件的形式实现, 也可以釆用软件功能模块的形式实现。 本发明不限制于任 何特定形式的硬件和软件的结合。 One of ordinary skill in the art will appreciate that all or a portion of the above steps may be accomplished by a program instructing the associated hardware, such as a read-only memory, a magnetic disk, or an optical disk. Alternatively, all or part of the steps of the above embodiments may also be implemented using one or more integrated circuits. Correspondingly, each module/unit in the above embodiment may be implemented in the form of hardware or in the form of a software function module. The invention is not limited to any specific form of combination of hardware and software.
工业实用性 本发明的宽带同频干扰的估计方法和系统可以得到较为准确的干扰噪声 特征, 有利于提高干扰抑制的性能和数据检测的准确性。 并且, 本发明的宽 带同频干扰的抑制方法和系统基于较为准确的干扰噪声特征, 可以提高干扰 抑制的性能和数据检测的准确性。 Industrial Applicability The method and system for estimating the wideband co-channel interference of the present invention can obtain more accurate interference noise characteristics, which is advantageous for improving the performance of interference suppression and the accuracy of data detection. Moreover, the wideband co-channel interference suppression method and system of the present invention is based on relatively accurate interference noise characteristics, and can improve the performance of interference suppression and the accuracy of data detection.

Claims

权 利 要 求 书 Claim
1、 一种宽带同频干扰噪声估计的方法, 应用于正交频分复用(OFDM)或 正交频分多址 (OFDMA)系统的接收端, 在一干扰抑制区域内, 用该方法对其 中承载的一个数据流进行干扰噪声估计时, 该方法包括: 对该数据流对应的每一数据子载波, 将该数据流对应的各导频子载波上 的接收信号与该接收信号的共轭转置的乘积的加权平均, 作为该数据子载波 位置的干扰噪声协方差矩阵估计值; 以及 对该数据流对应的每一数据子载波, 将对该数据子载波位置的干扰噪声 协方差矩阵估计值进行对角加载后得到的结果, 作为该数据子载波位置的干 扰噪声协方差矩阵; 其中, 该干扰抑制区域为接收数据承载区域中的一时频二维资源块。 A method for estimating a wideband co-channel interference noise, which is applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, and is used in an interference suppression region. When the data stream is carried in the interference noise estimation, the method includes: ???the conjugate of the received signal on each pilot subcarrier corresponding to the data stream and the received signal for each data subcarrier corresponding to the data stream a weighted average of the product of the transposition, as an estimate of the interference noise covariance matrix for the data subcarrier position; and for each data subcarrier corresponding to the data stream, an interference noise covariance matrix estimate for the data subcarrier position The value obtained after the diagonal loading is used as the interference noise covariance matrix of the data subcarrier position; wherein the interference suppression region is a time-frequency two-dimensional resource block in the received data bearer region.
2、 如权利要求 1所述的方法, 其中: 所述对该数据流对应的每一数据子载波, 将该数据流对应的各导频子载 波上的接收信号与该接收信号的共轭转置的乘积的加权平均, 作为该数据子 载波位置的干扰噪声协方差矩阵估计值的步骤中, 釆用式(a )计算数据子载 波位置的干扰噪声协方差矩阵估计值:
Figure imgf000035_0001
2. The method according to claim 1, wherein: the data conjugate corresponding to the data stream, the conjugate of the received signal on each pilot subcarrier corresponding to the data stream and the received signal The weighted average of the products, as the interference noise covariance matrix estimate for the data subcarrier position, uses equation (a) to calculate the interference noise covariance matrix estimate for the data subcarrier position:
Figure imgf000035_0001
其中, 为该干扰抑制区域中该数据流对应的第 个数据子载波位 置的干扰噪声协方差矩阵估计值, · = 1,· · · , J为该干扰抑制区域中该数据流 对应的数据子载波个数; = 1, · · ·,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; ^为计算该干扰抑制区域中该数据流对应的第 _;个数据子载 波位置的干扰噪声协方差矩阵估计值时赋予 ( ( i)) (yp (i )f 的权值,Wherein, the estimated value of the interference noise covariance matrix of the first data subcarrier position corresponding to the data stream in the interference suppression region, · = 1, · · · , J is the data corresponding to the data stream in the interference suppression region The number of carriers; = 1, · · ·, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region; ^ is to calculate the _th data corresponding to the data stream in the interference suppression region The interference noise covariance matrix estimate for the subcarrier position is given the weight of ( ( i ) ) ( y p (i )f ,
∑βϋ = \ ; W为所述接收端在该干扰抑制区域中该数据流对应的第 ,个导 频子载波上的接收信号, 以及 ( 为 ( 的共轭转置矩阵。 Σβ ϋ = \; received subcarrier signal W is at the receiving end of the interference suppression area, the data stream corresponding to a first, pilot, and (for the (conjugate transpose matrix.
3、 如权利要求 2所述的方法, 按式 (a)计算数据子载波位置的干扰噪声 协方差矩阵估计值的步骤之前, 所述方法还包括: 将该干扰抑制区域划分为一个或多个干扰噪声估计单元, 每一干扰噪声 估计单元为一个时域二维资源块且其中包含至少一个导频子载波和一个数据 子载波; 按式 (a)计算数据子载波位置的干扰噪声协方差矩阵估计值时, 为同一干 扰噪声估计单元中各个导频子载波上的接收信号与其共轭转置的乘积, 赋予 相同的权值。 3. The method of claim 2, before calculating the interference noise covariance matrix estimate of the data subcarrier position according to equation (a), the method further comprising: dividing the interference suppression region into one or more An interference noise estimating unit, each of the interference noise estimating units is a time domain two-dimensional resource block and includes at least one pilot subcarrier and one data subcarrier; and calculating an interference noise covariance matrix of the data subcarrier position according to formula (a) When estimating the value, the same weight is given to the product of the received signal on each pilot subcarrier in the same interference noise estimation unit and its conjugate transposed.
4、 如权利要求 1所述的方法, 所述方法还包括: 所述在一干扰抑制区域内, 用该方法对其中承载的一个数据流进行干扰 噪声估计时, 还将该干扰抑制区域划分为 个干扰噪声估计单元,每一干扰 噪声估计单元为一个时域二维资源块且其中包含至少一个导频子载波和一数 据子载波, 为正整数; 所述对该数据流对应的每一数据子载波, 将该数据流对应的各导频子载 波上的接收信号与该接收信号的共轭转置的乘积的加权平均, 作为该数据子 载波位置的干扰噪声协方差矩阵估计值的步骤中, 釆用式( b )计算数据子载 波位置的干扰噪声协方差矩阵估计值: 4. The method according to claim 1, further comprising: when the interference noise estimation is performed on a data stream carried by the method in the interference suppression region, the interference suppression region is further divided into Interference noise estimation unit, each interference noise estimation unit is a time domain two-dimensional resource block and includes at least one pilot subcarrier and one data subcarrier, which are positive integers; each data corresponding to the data stream Subcarrier, the weighted average of the product of the received signal on each pilot subcarrier corresponding to the data stream and the conjugate transposed of the received signal as a step of estimating the interference noise covariance matrix of the data subcarrier position , (b) Calculate the estimated value of the interference noise covariance matrix of the data subcarrier position:
RR
Figure imgf000036_0001
其中,
Figure imgf000036_0001
among them,
O 第 m 个干扰噪声估计单元中该数据流对应的每一数据子载波位 置的干扰噪声协方差矩阵估计值, m = 1, 2, · · · ,Μ;  O The estimated value of the interference noise covariance matrix for each data subcarrier position corresponding to the data stream in the mth interference noise estimation unit, m = 1, 2, · · ·, Μ;
/为一循环变量, / = 1, 2, · · ·,Μ ; / is a loop variable, / = 1, 2, · · ·, Μ;
为第 /个干扰噪声估计单元包含的该数据流对应的导频子载波的索引 ζ· 的集合, = 1, · · ·,/, /为该干扰抑制区域中该数据流对应的导频子载波的个 数; 为所述接收端在该数据流对应的第 I个导频子载波上的接收信号; (广为 ()的共轭转置矩阵; 以及 fiml为计算 —D时,赋予第 /个干扰噪声估计单元中该数据流对应的各导 频子载波对应的 ( 与 ( 的乘积的权值, | |为
Figure imgf000037_0001
a set of indices ζ· of the pilot subcarriers corresponding to the data stream included in the first interference noise estimation unit, = 1, · · ·, /, / is the pilot corresponding to the data stream in the interference suppression region The number of carriers; a receiving signal on the first pilot subcarrier corresponding to the data stream at the receiving end; (a conjugate transposed matrix of wide (); and fi ml being a calculation - D , giving the first/interference noise Estimating the corresponding pilot subcarrier corresponding to the data stream in the unit (the weight of the product of (, | | is
Figure imgf000037_0001
包含的导频子载波的个数。  The number of pilot subcarriers included.
5、如权利要求 4所述的方法,其中,按式 (b)计算^— D釆用的权值 An/中,The method according to claim 4, wherein the weight An/ in the calculation of ^ -D is calculated according to the formula (b),
/ = 1,2,···,Μ , 大于等于其他的权值。 / = 1,2,···,Μ , greater than or equal to other weights.
6、 如权利要求 2所述的方法, 其中: 所述对该数据流对应的每一数据子载波, 将对该数据子载波位置的干扰 噪声协方差矩阵估计值进行对角加载后得到的结果, 作为该数据子载波位置 的干扰噪声协方差矩阵的步骤中, 釆用式( c )计算数据子载波位置的干扰噪 声协方差矩阵: 6. The method according to claim 2, wherein: the data interference obtained by diagonally loading the estimated value of the interference noise covariance matrix of the data subcarrier position for each data subcarrier corresponding to the data stream In the step of the interference noise covariance matrix of the data subcarrier position, the interference noise covariance matrix of the data subcarrier position is calculated by using equation (c):
RM_D(j) = aRM_D(j) + (c) 其中, RW_D /)为该干扰抑制区域中该数据流对应的第 个数据子载波位 置的干扰噪声协方差矩阵, · = 1,··· , J为该干扰抑制区域中该数据流对应的 数据子载波个数; ≥0; 为该干扰抑制区域中该数据流对应的第 _;个 数据子载波位置的干扰噪声协方差矩阵估计值; β>0 Λ表示 N xN¾的对角 矩阵, 以及 ¾表示所述接收端的接收天线的数目。 R M _ D (j) = aR M _ D (j) + (c) where R W _ D /) is the interference noise covariance matrix of the first data subcarrier position corresponding to the data stream in the interference suppression region , · = 1,··· , J is the number of data subcarriers corresponding to the data stream in the interference suppression region; ≥ 0; is the _th data subcarrier position corresponding to the data stream in the interference suppression region Interference noise covariance matrix estimate; β>0 Λ denotes a diagonal matrix of N xN 3⁄4 , and 3⁄4 denotes the number of receive antennas at the receiving end.
7、 如权利要求 6所述的方法, 其中: 述第 j个数据子载波位置的干扰噪声协方差矩阵为:
Figure imgf000037_0002
7. The method of claim 6, wherein: the interference noise covariance matrix of the jth data subcarrier position is:
Figure imgf000037_0002
其中, 0≤ ≤1; tr( MD(_/》表示对矩阵 求迹; /为该干扰抑制区 域内该数据流对应的导频子载波的个数, 1为^„^^¾的单位矩阵。 Where 0 ≤ ≤1; tr( MD (_/′ represents the trace of the matrix; / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region, 1 is ^„^^ 3⁄4 Unit matrix.
8、 如权利要求 1至 7中任一项所述的方法, 其中: 所述在一干扰抑制区域内, 用该方法对其中承载的一个数据流进行干扰 噪声估计时, 按以下方式计算该数据流对应的每一导频子载波位置的信道系 数估计值: 将该导频子载波上的接收信号与发送端在该导频子载波上发送的导频信 号的共轭相乘, 得到该导频子载波位置的信道系数估计值。 The method according to any one of claims 1 to 7, wherein: in the interference suppression area, when the method performs interference noise estimation on a data stream carried therein, the data is calculated as follows The channel coefficient estimation value of each pilot subcarrier position corresponding to the stream: multiplying the received signal on the pilot subcarrier by the conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier, to obtain the guide Estimated channel coefficient for the frequency subcarrier position.
9、 一种宽带同频干 4尤抑制的方法, 应用于正交频分复用 (OFDM )或正 交频分多址( OFDMA ) 系统的接收端, 在一干扰抑制区域内, 用该方法对 其中承载的一个数据流进行干扰抑制时, 该方法包括: 按权利要求 8中所述的干扰噪声估计方法, 得到该数据流对应的各导频 子载波位置的信道系数估计值和各数据子载波位置的干扰噪声协方差矩阵; 对该数据流对应的每一数据子载波, 将该数据流对应的各个导频子载波 位置的信道系数估计值的加权平均, 作为该数据子载波位置的信道系数估计 值; 以及 对该数据流对应的每一数据子载波, 根据该数据子载波上的接收信号, 及该数据子载波位置的信道系数估计值和干扰噪声协方差矩阵, 计算得到该 数据子载波上的数据信号估计; 其中, 该干扰抑制区域为接收数据承载区域中的一时频二维资源块。 9. A method for wideband co-frequency interference suppression, applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, in an interference suppression region, using the method When performing interference suppression on a data stream carried in the method, the method includes: estimating, by using the interference noise estimation method according to claim 8, the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream, and each data sub An interference noise covariance matrix of a carrier position; a weighted average of channel coefficient estimates of respective pilot subcarrier positions corresponding to the data stream for each data subcarrier corresponding to the data stream, as a channel of the data subcarrier position An estimated value of the coefficient; and each data subcarrier corresponding to the data stream, the data is calculated according to the received signal on the data subcarrier, and the channel coefficient estimation value and the interference noise covariance matrix of the data subcarrier position. The data signal estimation on the carrier; wherein the interference suppression area is a time-frequency two-dimensional resource block in the received data bearer area.
10、 如权利要求 9所述的方法, 所述方法还包括: 对该数据流对应的每一数据子载波, 在计算得到该数据载波上的数据信 号估计后, 对该估计出的数据信号进行解调, 并将解调后得到的该信号各比 特对应的软信息 L[,…'!F进行调整, 调整后各比特对应的软信息分别为 10. The method according to claim 9, the method further comprising: performing, after calculating the data signal estimation on the data carrier for each data subcarrier corresponding to the data stream, performing the estimated data signal Demodulate, and demodulate the soft information corresponding to each bit of the signal L[,...'! F is adjusted, and the soft information corresponding to each bit after adjustment is
其中, ' = l, · · · , J , J为该干扰抑制区域中该数据流对应的数据子载波个数; F为该干扰抑制区域中该数据流对应的第 J个数据载波上的数据信号中包含 的比特个数, 表示该干扰抑制区域中该数据流对应的第 _/·个数据载波对 应 的 信 道 增 益 4 ( )(RW— ( ) 或 者 信 干 噪 比 值
Figure imgf000039_0001
Wherein, ' = l, · · · , J , J are the number of data subcarriers corresponding to the data stream in the interference suppression region; F is the data on the Jth data carrier corresponding to the data stream in the interference suppression region The number of bits included in the signal indicates the _/·th data carrier pair corresponding to the data stream in the interference suppression region Shoulder channel gain 4 ( ) (R W — ( ) or signal to interference and noise ratio
Figure imgf000039_0001
hd {j)为该干扰抑制区域中该数据流对应的第 j个数据子载波位置的信道 系数估计值, /)为 4( ·;)的共轭转置矩阵, RM_D /;>为第 ·个数据子载波位 置的干扰噪声协方差矩阵, (RwD(_/))— 1为 RM— 的逆矩阵; = 1, ···,/, /为该干扰抑制区域中该数据流对应的导频子载波个数; ft ή 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道系数估计 值, 为 的共轭转置矩阵, M_PW为第 ,个导频子载波位置的干扰 噪声协方差矩阵, 1为^ 的逆矩阵; ()为所述接收端在该干 扰抑制区域中该数据流对应的第 ,个导频子载波上的接收信号, 以及; 为 发送端在该干扰抑制区域中该数据流对应的第 , 个导频子载波上的发射信 号。 h d {j) is an estimated channel coefficient of the jth data subcarrier position corresponding to the data stream in the interference suppression region, /) is a conjugate transposed matrix of 4 (·;), R M _ D /; > is the interference noise covariance matrix of the first data subcarrier position, (R wD (_/)) — 1 is the inverse matrix of R M —; 1, 1, ···, /, / is the interference suppression The number of pilot subcarriers corresponding to the data stream in the region; ft ή is the channel coefficient estimation value of the first pilot subcarrier position corresponding to the data stream in the interference suppression region, and is a conjugate transposed matrix, M _ P W is the interference noise covariance matrix of the first pilot subcarrier position, where 1 is the inverse matrix of ^; () is the first pilot of the data stream corresponding to the receiving end in the interference suppression region a received signal on the carrier, and a transmit signal on the first pilot subcarrier corresponding to the data stream in the interference suppression region by the transmitting end.
11、 如权利要求 9或 10所述的方法, 其中: 所述在一干扰抑制区域内, 用该方法对其中承载的一个数据流进行干扰 抑制时, 还将该干扰抑制区域划分为 f个信道估计单元, 每一信道估计单元 为一个时域二维资源块且其中包含至少一个导频子载波和一个数据子载波, f为正整数; 所述对该数据流对应的每一数据子载波, 将该数据流对应的各个导频子 载波位置的信道系数估计值的加权平均, 作为该数据子载波位置的信道系数 估计值的步骤, 釆用的计算公式如下:
Figure imgf000039_0002
其中,
The method according to claim 9 or 10, wherein: in the interference suppression area, when the interference suppression is performed on a data stream carried by the method, the interference suppression area is further divided into f channels. An estimation unit, each channel estimation unit is a time domain two-dimensional resource block and includes at least one pilot subcarrier and one data subcarrier, where f is a positive integer; and each data subcarrier corresponding to the data stream is The weighted average of the channel coefficient estimation values of the respective pilot subcarrier positions corresponding to the data stream is used as the step of estimating the channel coefficient of the data subcarrier position, and the calculation formula is as follows:
Figure imgf000039_0002
among them,
为第 个信道估计单元中该数据流对应的每一数据子载波位置的信道 系数估计值, k = \,2,--、K Channel for each data subcarrier location corresponding to the data stream in the first channel estimation unit Coefficient estimate, k = \, 2, --, K
/为一循环变量, ! = 1,2,···, K / is a loop variable, ! = 1,2,···, K
Ω,为第 个信道估计单元包含的该干扰抑制区域中该数据流对应的导频 子载波的索引 ,的集合, = ι,···,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道 系数估计值;  Ω, which is a set of indices of pilot subcarriers corresponding to the data stream in the interference suppression region included in the first channel estimation unit, = ι,···, /, / is corresponding to the data stream in the interference suppression region The number of pilot subcarriers; the estimated channel coefficient of the first pilot subcarrier position corresponding to the data stream in the interference suppression region;
Okl为计算 时,赋予第 /个信道估计单元中各导频子载波位置的 ίιΡ ()的O kl is calculated, giving the first / channel estimation unit for each pilot sub-carrier position ίι Ρ () of
Κ Κ
权值, ∑|Ω=1,0≤ ≤1, | |表示 包含的导频子载波的个数, 且在权 Weight, ∑| Ω =1 , 0 ≤ ≤1 , | | indicates the number of pilot subcarriers included, and is in the right
1=1  1=1
值" w中, l = \,2 ,K , 大于等于其他的权值。 The value "w, l = \, 2, K, is greater than or equal to other weights.
12、 一种宽带同频干扰噪声估计的系统, 应用于正交频分复用(OFDM) 或正交频分多址 (OFDMA)系统的接收端,在一干扰抑制区域内对其中承载的 一个数据流进行干扰噪声估计, 该干扰抑制区域为接收数据承载区域中的一 时频二维资源块, 该系统包括: 第一装置, 其设置为: 对该数据流对应的每一数据子载波, 将该数据流 对应的各导频子载波上的接收信号与该接收信号的共轭转置的乘积的加权平 均, 作为该数据子载波位置的干扰噪声协方差矩阵估计值; 以及 第二装置, 其设置为: 对该数据流对应的每一数据子载波, 将对所述第 一装置得到的该数据子载波位置的干扰噪声协方差矩阵估计值进行对角加载 后得到的结果, 作为该数据子载波位置的干扰噪声协方差矩阵。 12. A system for wideband co-channel interference noise estimation applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, one of which is carried in an interference suppression region The data stream performs interference noise estimation, where the interference suppression area is a time-frequency two-dimensional resource block in the received data bearer area, the system includes: a first device, configured to: each data subcarrier corresponding to the data stream, a weighted average of a product of a received signal on each pilot subcarrier corresponding to the data stream and a conjugate transpose of the received signal as an estimated value of an interference noise covariance matrix of the data subcarrier position; and a second device And being configured to: perform, after each of the data subcarriers corresponding to the data stream, a result obtained by diagonally loading an interference noise covariance matrix estimation value of the data subcarrier position obtained by the first device, as the data sub Interference noise covariance matrix for carrier position.
13、 如权利要求 12所述的系统, 其中: 所述第一装置是设置为釆用式(a)计算出该数据子载波位置的干扰噪声 协方差矩 :
Figure imgf000040_0001
其中, 为该干扰抑制区域中该数据流对应的第 _;个数据子载波位 置的干扰噪声协方差矩阵估计值, · = 1,· · · , J为该干扰抑制区域中该数据流 对应的数据子载波个数; = 1, · · ·,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数; ^为计算该干扰抑制区域中该数据流对应的第 _;个数据子载 波位置的干扰噪声协方差矩阵估计值时赋予 ( ( i)) (yp (i )f 的权值, ∑βϋ = \ ; W为所述接收端在该干扰抑制区域中该数据流对应的第 ,个导 频子载波上的接收信号, (if为 ( 的共轭转置矩阵。
13. The system of claim 12, wherein: said first device is configured to calculate an interference noise covariance moment of said data subcarrier position using equation (a):
Figure imgf000040_0001
Wherein, the estimated value of the interference noise covariance matrix of the data subcarrier position corresponding to the data stream in the interference suppression region, · = 1, · · · , J is the corresponding data stream in the interference suppression region The number of data subcarriers; = 1, · · ·, /, / is the number of pilot subcarriers corresponding to the data stream in the interference suppression region; ^ is to calculate the _ corresponding to the data stream in the interference suppression region; when a given interference noise covariance matrix estimation value of the data sub-carrier position ((I)) weight (y p (i) f the, Σβ ϋ = \; W to the receiving end in the area of the interference suppression The received signal on the first pilot subcarrier corresponding to the data stream, (if is the conjugate transposed matrix of .
14、 如权利要求 13所述的系统, 所述系统还包括第三装置, 其中: 所述第三装置设置为: 将该干扰抑制区域划分为一个或多个干扰噪声估 计单元, 每一干扰噪声估计单元为一个时域二维资源块且其中包含至少一个 导频子载波和一个数据子载波; 所述第一装置按式 (a)计算数据子载波位置的干扰噪声协方差矩阵估计 值时, 为同一干扰噪声估计单元中各个导频子载波上的接收信号与其共轭转 置的乘积, 赋予相同的权值。 14. The system of claim 13, the system further comprising a third device, wherein: the third device is configured to: divide the interference suppression region into one or more interference noise estimation units, each interference noise The estimating unit is a time domain two-dimensional resource block and includes at least one pilot subcarrier and one data subcarrier; when the first device calculates the interference noise covariance matrix estimation value of the data subcarrier position according to formula (a), The same weight is assigned to the product of the received signal on each pilot subcarrier in the same interference noise estimation unit and its conjugate transpose.
15、 如权利要求 12所述的系统, 所述系统还包括第四装置, 其中: 所述第四装置设置为: 将该干扰抑制区域划分为 个干扰噪声估计单 元, 每一干扰噪声估计单元为一个时域二维资源块且其中包含至少一个导频 子载波和一数据子载波, 为正整数; 对该数据流对应的每一数据子载波, 所述第一装置是设置为按如下公式 将该数据流对应的各导频子载波上的接收信号与该接收信号的共轭转置的乘 积的加权平均, 作为该数据子载波位置的干扰噪声协方差矩阵估计值:
Figure imgf000041_0001
其中,
The system of claim 12, the system further comprising a fourth device, wherein: the fourth device is configured to: divide the interference suppression region into interference noise estimation units, each interference noise estimation unit is a time domain two-dimensional resource block and including at least one pilot subcarrier and one data subcarrier, which are positive integers; for each data subcarrier corresponding to the data stream, the first device is set to be according to the following formula A weighted average of the product of the received signal on each pilot subcarrier corresponding to the data stream and the conjugate transposed of the received signal as the estimated value of the interference noise covariance matrix of the data subcarrier position:
Figure imgf000041_0001
among them,
U m 个干扰噪声估计单元中该数据流对应的每一数据子载波位 置的干扰噪声协方差矩阵估计值, m = 1, 2, · · · ,Μ; /为一循环变量, / = 1,2,···,Μ; Each data subcarrier bit corresponding to the data stream in the U m interference noise estimation units The estimated value of the interference noise covariance matrix, m = 1, 2, · · · , Μ; / is a cyclic variable, / = 1, 2, ···, Μ;
为第 /个干扰噪声估计单元包含的该数据流对应的导频子载波的索引 ζ· 的集合, = 1,···,/, /为该干扰抑制区域中该数据流对应的导频子载波的个 数; a set of indices ζ· of pilot subcarriers corresponding to the data stream included in the first interference noise estimation unit, = 1,···, /, / is the pilot corresponding to the data stream in the interference suppression region The number of carriers;
^;)为所述接收端在该数据流对应的第 I个导频子载波上的接收信号; yp (if为 ()的共轭转置矩阵; fiml为计算 —D时,赋予第 /个干扰噪声估计单元中该数据流对应的各导 频子载波对应的 ( 与 ( 的乘积的权值, | |为
Figure imgf000042_0001
^;) is the received signal of the receiving end on the first pilot subcarrier corresponding to the data stream; y p (if is the conjugate transposed matrix of (); fi ml is calculated - D , the first Corresponding to each pilot subcarrier corresponding to the data stream in the interference noise estimation unit (with the weight of the product of (, | |
Figure imgf000042_0001
包含的导频子载波的个数, 大于等于其他的权值。  The number of pilot subcarriers included is greater than or equal to other weights.
16、 如权利要求 13所述的系统, 其中: 对该数据流对应的每一数据子载波, 所述第二装置是设置为按如下公式 将对该数据子载波位置的干扰噪声协方差矩阵估计值进行对角加载后得到的 结果, 作为该数据子载波位置的干扰噪声协方差矩阵: 16. The system of claim 13, wherein: for each data subcarrier corresponding to the data stream, the second device is configured to estimate an interference noise covariance matrix for the data subcarrier position according to the following formula The value obtained after diagonal loading is used as the interference noise covariance matrix of the data subcarrier position:
RN1_D(j) = aR NI-D R N1 _ D (j) = aR NI-D
其中, RW_D /)为该干扰抑制区域中该数据流对应的第 _;个数据子载波位 置的干扰噪声协方差矩阵, · = 1,··· , J为该干扰抑制区域中该数据流对应的 数据子载波个数; ≥0; 为该干扰抑制区域中该数据流对应的第 _;个 数据子载波位置的干扰噪声协方差矩阵估计值; β>0 Λ表示 N xN¾的对角 矩阵, ¾表示所述接收端的接收天线的数目。 Where R W _ D /) is the interference noise covariance matrix of the _th data subcarrier position corresponding to the data stream in the interference suppression region, · = 1,···, J is the interference suppression region The number of data subcarriers corresponding to the data stream; ≥ 0; is the estimated value of the interference noise covariance matrix of the _; data subcarrier positions corresponding to the data stream in the interference suppression region; β>0 Λ represents N xN 3⁄4 The diagonal matrix, 3⁄4, represents the number of receiving antennas at the receiving end.
17、 如权利要求 16所述的系统, 其中: 所述第二装置计算出的第 个数据子载波位置的干扰噪声协方差矩阵为:
Figure imgf000042_0002
17. The system of claim 16, wherein: the interference noise covariance matrix of the first data subcarrier position calculated by the second device is:
Figure imgf000042_0002
其中, 0≤ ≤1; tr( MD(_/》表示对矩阵 求迹; /为该干扰抑制区 域内该数据流对应的导频子载波的个数, 1为^„^ 的单位矩阵。 Where 0 ≤ ≤1; tr( MD (_/′ represents the trace of the matrix; / is the interference suppression zone The number of pilot subcarriers corresponding to the data stream in the domain, where 1 is the identity matrix of ^^^.
18、 如权利要求 12至 17中任一项所述的系统, 所述系统还包括第五装 置, 其中: 所述第五装置设置为按以下方式得到发送端在该数据流对应的每一导频 子载波位置的信道系数估计值并输出到所述第一装置: 将该导频子载波上的接收信号与发送端在该导频子载波上发送的导频信 号的共轭相乘, 得到该导频子载波位置的信道系数估计值。 The system according to any one of claims 12 to 17, the system further comprising a fifth device, wherein: the fifth device is configured to obtain, according to the following manner, each of the transmitting ends corresponding to the data stream Estimating the channel coefficient of the frequency subcarrier position and outputting to the first device: multiplying the received signal on the pilot subcarrier by a conjugate of the pilot signal transmitted by the transmitting end on the pilot subcarrier, The estimated channel coefficient of the pilot subcarrier position.
19、 一种宽带同频干 4尤抑制的系统, 应用于正交频分复用 (OFDM )或 正交频分多址(OFDMA ) 系统的接收端, 在一干扰抑制区域内对其中承载 的一个数据流进行干扰抑制, 该干扰抑制区域为接收数据承载区域中的一时 频二维资源块, 该系统包括: 第一子系统,其设置为:按权利要求 18所述的宽带同频干扰噪声估计的 系统相同的方式, 得到该数据流对应的各导频子载波位置的信道系数估计值 和各数据子载波位置的干扰噪声协方差矩阵; 第二子系统, 其设置为: 对该数据流对应的每一数据子载波, 将该数据 流对应的各个导频子载波位置的信道系数估计值的加权平均, 作为该数据子 载波位置的信道系数估计值; 以及 第三子系统, 其设置为: 对该数据流对应的每一数据子载波, 根据该数 据子载波上的接收信号, 及该数据子载波位置的信道系数估计值和干扰噪声 协方差矩阵, 计算得到该数据子载波上的数据信号估计。 19. A system for wideband co-frequency interference suppression, applied to a receiving end of an Orthogonal Frequency Division Multiplexing (OFDM) or Orthogonal Frequency Division Multiple Access (OFDMA) system, carried in an interference suppression region Interference suppression is performed on a data stream, the interference suppression area is a time-frequency two-dimensional resource block in the received data bearer area, and the system includes: a first subsystem configured to: the broadband co-channel interference noise according to claim 18. The estimated system obtains the channel coefficient estimation value of each pilot subcarrier position corresponding to the data stream and the interference noise covariance matrix of each data subcarrier position; the second subsystem is set to: the data stream Corresponding to each data subcarrier, a weighted average of channel coefficient estimates of respective pilot subcarrier positions corresponding to the data stream is used as an estimated channel coefficient of the data subcarrier position; and a third subsystem is set to : for each data subcarrier corresponding to the data stream, according to the received signal on the data subcarrier, and the channel system of the data subcarrier position And estimate interference noise covariance matrix, the calculated data signal on the data subcarrier estimated.
20、 如权利要求 19所述的系统, 所述系统还包括第四子系统, 其中: 所述第四子系统设置为: 对该数据流对应的每一数据子载波, 在所述第 三子系统计算得到该数据载波上的数据信号估计后, 对该估计出的数据信号 进行解调, 并将解调后得到的该信号各比特对应的软信息 进行调整, 调整后各比特对应的软信息分别为
Figure imgf000043_0001
; 其中, ' = l, ···, J, J为该干扰抑制区域中该数据流对应的数据子载波个数;
20. The system of claim 19, the system further comprising a fourth subsystem, wherein: the fourth subsystem is configured to: each of the data subcarriers corresponding to the data stream, in the third sub After the system calculates the data signal on the data carrier, the estimated data signal is demodulated, and the soft information corresponding to each bit of the signal obtained after the demodulation is adjusted, and the soft information corresponding to each bit is adjusted. Separately
Figure imgf000043_0001
; Wherein, ' = l, ···, J, J is the number of data subcarriers corresponding to the data stream in the interference suppression region;
F为该干扰抑制区域中该数据流对应的第 J个数据载波上的数据信号中包含 的比特个数, 表示该干扰抑制区域中该数据流对应的第 _/·个数据载波对F is the number of bits included in the data signal on the Jth data carrier corresponding to the data stream in the interference suppression area, and indicates the _/·th data carrier pair corresponding to the data stream in the interference suppression area.
( )(RW— (7) 或 者 信 干 噪 比 值
Figure imgf000044_0001
( ) (R W — (7) or signal to interference and noise ratio
Figure imgf000044_0001
hd {j)为该干扰抑制区域中该数据流对应的第 j个数据子载波位置的信道 系数估计值, /)为 4( ·;)的共轭转置矩阵, RM_D /;>为第 ·个数据子载波位 置的干扰噪声协方差矩阵, (RMD ( ))'1为 RM_D (_/·)的逆矩阵; = 1, ···,/, /为该干扰抑制区域中该数据流对应的导频子载波个数; ft ή 为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道系数估计 值, 为 的共轭转置矩阵, M_PW为第 ,个导频子载波位置的干扰 噪声协方差矩阵, 1为^ 的逆矩阵; ()为所述接收端在该干 扰抑制区域中该数据流对应的第 ,个导频子载波上的接收信号, 为发送 端在该干扰抑制区域中该数据流对应的第 I个导频子载波上的发射信号。 h d {j) is an estimated channel coefficient of the jth data subcarrier position corresponding to the data stream in the interference suppression region, /) is a conjugate transposed matrix of 4 (·;), R M _ D /; > is the interference noise covariance matrix of the first data subcarrier position, (R MD ( )) ' 1 is the inverse matrix of R M _ D (_/·); = 1, ···, /, / The number of pilot subcarriers corresponding to the data stream in the interference suppression region; ft ή is the channel coefficient estimation value of the first pilot subcarrier position corresponding to the data stream in the interference suppression region, and is a conjugate rotation The matrix, M _ P W is the interference noise covariance matrix of the first pilot subcarrier position, where 1 is the inverse matrix of ^; () is the number corresponding to the data stream in the interference suppression region of the receiving end, The received signal on the pilot subcarrier is the transmit signal on the first pilot subcarrier corresponding to the data stream in the interference suppression region.
21、 如权利要求 19或 20所述的系统, 所述系统还包括第五子系统和第 六子系统, 其中: 所述第五子系统设置为: 将该干扰抑制区域划分为 K个信道估计单元, 每一信道估计单元为一个时域二维资源块且其中包含至少一个导频子载波和 一个数据子载波, f为正整数; 所述第六子系统设置为: 对该数据流对应的每一数据子载波, 釆用如下 公式将该数据流对应的各个导频子载波位置的信道系数估计值的加权平均, 作为该数据子载波位置的信道系数估计值:
Figure imgf000044_0002
21. The system of claim 19 or 20, the system further comprising a fifth subsystem and a sixth subsystem, wherein: the fifth subsystem is configured to: divide the interference suppression region into K channel estimates a unit, each channel estimation unit is a time domain two-dimensional resource block and includes at least one pilot subcarrier and one data subcarrier, where f is a positive integer; the sixth subsystem is set to: corresponding to the data stream For each data subcarrier, the weighted average of the channel coefficient estimation values of the respective pilot subcarrier positions corresponding to the data stream is used as the channel coefficient estimation value of the data subcarrier position by using the following formula:
Figure imgf000044_0002
/=1 ie^ 为第 个信道估计单元中该数据流对应的每一数据子载波位置的信道 系数估计值, k = \,2,---,K /=1 ie^ The channel coefficient estimate for each data subcarrier location corresponding to the data stream in the first channel estimation unit, k = \, 2, ---, K
/为一循环变量, ! = 1,2,···, K / is a loop variable, ! = 1,2,···, K
Ω,为第 个信道估计单元包含的该干扰抑制区域中该数据流对应的导频 子载波的索引 ,的集合, = ι,···,/, /为该干扰抑制区域中该数据流对应的导 频子载波个数;  Ω, which is a set of indices of pilot subcarriers corresponding to the data stream in the interference suppression region included in the first channel estimation unit, = ι,···, /, / is corresponding to the data stream in the interference suppression region Number of pilot subcarriers;
为该干扰抑制区域中该数据流对应的第 I个导频子载波位置的信道 系数估计值;  An estimated channel coefficient of the first pilot subcarrier position corresponding to the data stream in the interference suppression region;
Okl为计算 时,赋予第 /个信道估计单元中各导频子载波位置的 ίιΡ ()的O kl is calculated, giving the first / channel estimation unit for each pilot sub-carrier position ίι Ρ () of
Κ Κ
权值, ∑|Ω= 1,0≤ ≤ 1, | |表示 包含的导频子载波的个数, 且在权 Weight, ∑| Ω ′′ = 1 , 0 ≤ ≤ 1 , | | indicates the number of pilot subcarriers included, and is in the right
1=1  1=1
值 0¾中, ί = 1,2,···,Κ , 大于等于其他的权值。 In value 03⁄4, ί = 1,2,···,Κ , is greater than or equal to other weights.
PCT/CN2011/075517 2010-12-15 2011-06-09 Methods for estimating broadband co-channel interference and noise and suppressing interference, and systems thereof WO2012079348A1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN201010590051.6A CN102546483B (en) 2010-12-15 2010-12-15 Method for broadband common-frequency interference noise estimation and interference suppression and corresponding system
CN201010590051.6 2010-12-15

Publications (1)

Publication Number Publication Date
WO2012079348A1 true WO2012079348A1 (en) 2012-06-21

Family

ID=46244034

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2011/075517 WO2012079348A1 (en) 2010-12-15 2011-06-09 Methods for estimating broadband co-channel interference and noise and suppressing interference, and systems thereof

Country Status (2)

Country Link
CN (1) CN102546483B (en)
WO (1) WO2012079348A1 (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103200138A (en) * 2013-04-03 2013-07-10 北京航空航天大学 Noise variance estimation method

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9231667B2 (en) 2013-02-22 2016-01-05 Futurewei Technologies, Inc. Systems and methods for a long-term evolution network
CN103227760B (en) * 2013-04-28 2016-03-16 中国铁路通信信号股份有限公司 Channel estimation methods under a kind of high-speed mobile environment
CN105338632B (en) * 2014-06-20 2020-09-29 索尼公司 Radio resource management system and radio resource management method
CN114257256B (en) * 2021-12-15 2023-08-15 哲库科技(北京)有限公司 Noise estimation method, apparatus, device, and readable storage medium
CN115865109B (en) * 2022-11-11 2024-05-14 北京智芯微电子科技有限公司 Interference suppression combining method, device and medium for multiple receiving antennas and receiving terminal

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101726730A (en) * 2009-12-07 2010-06-09 中国人民解放军空军雷达学院 Self-adaption anti-coherent interference technology based on characteristic component rejection
CN101753176A (en) * 2009-12-24 2010-06-23 北京北方烽火科技有限公司 Interference rejection combining method and system
WO2010119031A1 (en) * 2009-04-14 2010-10-21 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for direct estimation of impairment correlations for received signal processing

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101155156B (en) * 2006-09-25 2011-06-08 华为技术有限公司 Channel estimation method and device and pilot frequency sequence generation method and device

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2010119031A1 (en) * 2009-04-14 2010-10-21 Telefonaktiebolaget Lm Ericsson (Publ) Method and apparatus for direct estimation of impairment correlations for received signal processing
CN101726730A (en) * 2009-12-07 2010-06-09 中国人民解放军空军雷达学院 Self-adaption anti-coherent interference technology based on characteristic component rejection
CN101753176A (en) * 2009-12-24 2010-06-23 北京北方烽火科技有限公司 Interference rejection combining method and system

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103200138A (en) * 2013-04-03 2013-07-10 北京航空航天大学 Noise variance estimation method

Also Published As

Publication number Publication date
CN102546483B (en) 2014-12-31
CN102546483A (en) 2012-07-04

Similar Documents

Publication Publication Date Title
CN102447656B (en) Method for processing received OFDM data symbols and OFDM baseband receiver
CN102882575B (en) For determining the method and apparatus of channel condition information
CN101232474B (en) Apparatus and method for receiving signal in broadband wireless access system
CN102480444B (en) Method and corresponding system for broadband co-frequency interference noise estimation and interference suppression
CN101753176B (en) Interference rejection combining method and system
WO2011090028A1 (en) Communication apparatus and base station apparatus
CN102571659B (en) Interference noise estimation and interference inhibition method and corresponding system
WO2012079348A1 (en) Methods for estimating broadband co-channel interference and noise and suppressing interference, and systems thereof
WO2017097269A1 (en) Interference estimation method and device
WO2011035594A1 (en) Demodulation method and demodulator for orthogonal frequency division multiplexing - multiple input multiple output system
EP2713542B1 (en) Method and transmitter apparatus for composing a transmit signal and method and receiver apparatus for channel estimation
EP2608443B1 (en) Communications terminal, apparatus, and method for detecting rank indication
WO2005122447A1 (en) A method of realizing channel estimation in ofdm system
CN113660071A (en) Interference estimation method, device, equipment and medium
CN105099610A (en) Signal processing method and apparatus
WO2012092751A1 (en) Method and system for neighboring cell interference detection
WO2013174850A1 (en) Method and arrangement in a communications system
WO2012106963A1 (en) Method and device for eliminating interference and noise
JP5282042B2 (en) Method and apparatus for impairment correlation evaluation in multi-antenna receiver
US8687675B2 (en) Method and system for transmitting/receiving data in communication system
WO2012068871A1 (en) Method and corresponding system for broadband co-frequency interference noise estimation and interference suppression
CN102685060B (en) Multi-user multiple input multiple output (MIMO) receiving method and device for orthogonal frequency division multiplexing system
JP5904210B2 (en) Wireless communication apparatus and communication method
CN102447658B (en) Method and device for suppressing and merging interference
WO2012068865A1 (en) Channel estimation method and system in broadband co-frequency interference environment

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 11847929

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 11847929

Country of ref document: EP

Kind code of ref document: A1