WO2012117788A1 - 送信装置、受信装置、通信システム、通信方法、および集積回路 - Google Patents

送信装置、受信装置、通信システム、通信方法、および集積回路 Download PDF

Info

Publication number
WO2012117788A1
WO2012117788A1 PCT/JP2012/051923 JP2012051923W WO2012117788A1 WO 2012117788 A1 WO2012117788 A1 WO 2012117788A1 JP 2012051923 W JP2012051923 W JP 2012051923W WO 2012117788 A1 WO2012117788 A1 WO 2012117788A1
Authority
WO
WIPO (PCT)
Prior art keywords
signal
unit
demodulation reference
reference signal
dmrs
Prior art date
Application number
PCT/JP2012/051923
Other languages
English (en)
French (fr)
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 シャープ株式会社
Priority to US14/002,142 priority Critical patent/US9362995B2/en
Publication of WO2012117788A1 publication Critical patent/WO2012117788A1/ja

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0452—Multi-user MIMO systems
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/0204—Channel estimation of multiple channels
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/0224—Channel estimation using sounding signals
    • H04L25/0226—Channel estimation using sounding signals sounding signals per se
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202—Channel estimation
    • H04L25/024—Channel estimation channel estimation algorithms
    • H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00—Baseband systems
    • H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/06—Dc level restoring means; Bias distortion correction ; Decision circuits providing symbol by symbol detection
    • H04L25/067—Dc level restoring means; Bias distortion correction ; Decision circuits providing symbol by symbol detection providing soft decisions, i.e. decisions together with an estimate of reliability
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L5/00—Arrangements affording multiple use of the transmission path
    • H04L5/003—Arrangements for allocating sub-channels of the transmission path
    • H04L5/0048—Allocation of pilot signals, i.e. of signals known to the receiver
    • H04L5/005—Allocation of pilot signals, i.e. of signals known to the receiver of common pilots, i.e. pilots destined for multiple users or terminals
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04B—TRANSMISSION
    • H04B7/00—Radio transmission systems, i.e. using radiation field
    • H04B7/02—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04—Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413—MIMO systems
    • H04B7/0456—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • H04B7/046—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking physical layer constraints into account
    • H04B7/0465—Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting taking physical layer constraints into account taking power constraints at power amplifier or emission constraints, e.g. constant modulus, into account
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04J—MULTIPLEX COMMUNICATION
    • H04J11/00—Orthogonal multiplex systems, e.g. using WALSH codes
    • H04J11/0023—Interference mitigation or co-ordination
    • H04J11/0026—Interference mitigation or co-ordination of multi-user interference
    • H04J11/003—Interference mitigation or co-ordination of multi-user interference at the transmitter

Definitions

  • the present invention relates to a transmission device, a reception device, a communication system, a communication method, and an integrated circuit.
  • MIMO Multi-Input Multi-Output
  • ⁇ LP MU-MIMO> In MIMO, a terminal different from Single-User MIMO (SU-MIMO) in which a base station device (Base Station: BS) transmits a plurality of signals to one terminal device (Mobile Station: MS) at the same time and the same frequency. There is Multi-User MIMO (MU-MIMO) that transmits signals to devices at the same time and the same frequency.
  • SU-MIMO Single-User MIMO
  • MS Mobile Station: MS
  • MU-MIMO Multi-User MIMO
  • Non-Patent Document 1 LTE (Long Term Evolution) and LTE-Advanced, downlink (DL) MU-MIMO using linear precoding (LP) has already been specified (see Non-Patent Document 1 below).
  • LTE Long Term Evolution
  • DL downlink
  • LP linear precoding
  • ⁇ NLP MU-MIMO> a base station device multiplies a linear filter so that transmission signals are orthogonalized and interference between terminal devices (Multi-User Interference: MUI) is reduced.
  • MUI Multi-User Interference
  • Nonlinear Precoding (NLP) MU-MIMO has been proposed as another method for realizing spatial multiplexing.
  • the terminal device receives signals in an in-phase component (I-ch) and quadrature component (Quadrature channel: Q-ch) directions in an integral multiple of a constant width (Modulo width).
  • I-ch in-phase component
  • Q-ch quadrature channel
  • a modulo operation is performed in which the points that are translated in parallel are regarded as the same point.
  • the base station apparatus can add a signal (perturbation vector) that is an integer multiple of the modulo width to the modulation signal, appropriately select the perturbation vector, add it to the signal addressed to each terminal apparatus, and increase the transmission power.
  • VP MU-MIMO When the terminal apparatus performs a modulo operation on the received signal, the base station apparatus has a degree of freedom to add a signal having an arbitrary integral multiple of the modulo width to each modulated signal. This signal that can be added is called a perturbation vector.
  • a VP (Vector Perturbation) MU-MIMO method is a method in which all perturbation vectors that improve the power efficiency are searched in consideration of the propagation path states of all the terminal devices that are spatially multiplexed.
  • VP MU-MIMO is an NLP MU-MIMO system that can obtain a full transmission diversity gain and exhibits very good characteristics although the amount of computation of the base station apparatus is large (see Non-Patent Document 2 below).
  • THP Tomlinson-Harashima precoding
  • MU-MIMO is a method for calculating a perturbation vector to be sequentially added to a signal addressed to each terminal device in consideration of inter-user interference received by each terminal device. Call.
  • THP MU-MIMO has low complexity of the transmission processing of the base station apparatus, full transmission diversity cannot be obtained in all terminal apparatuses (see Non-Patent Document 3 below).
  • LR-THP is a method capable of obtaining full transmission diversity gain with a smaller amount of computation than VP MU-MIMO by adding processing called lattice reduction (LR) to THP MU-MIMO ( Non-patent document 3 below).
  • DMRS Demodulation Reference Signal
  • DMRS is a signal for the base station apparatus to notify each terminal apparatus in advance of the amplitude and phase of the data signal precoded by NLP MU-MIMO.
  • the base station apparatus adds a perturbation vector to the DMRS (or performs a modulo operation) and transmits the DMRS. Therefore, the terminal apparatus also performs a modulo operation on the DMRS. There is a need. Therefore, the terminal device needs to know in advance the modulo width necessary for the modulo calculation. Since the modulo width is proportional to the amplitude of the modulation signal in the reception signal, the terminal device needs to know the reception gain (complex gain of the propagation path) of the data signal after nonlinear precoding.
  • the terminal apparatus cannot know the reception gain (complex gain of the propagation path) without performing propagation path estimation using DMRS. That is, the terminal apparatus is in a state that “a reception gain (complex gain of a propagation path) cannot be obtained unless a propagation path is estimated from DMRS, but a DMRS cannot be estimated without knowing the reception gain”. Therefore, the above problem occurs.
  • each terminal apparatus is configured with orthogonal radio resources (areas divided in the time direction and the frequency direction and do not interfere with each other even if different data signals and reference signals are allocated).
  • Send DMRS since the terminal device does not need to perform a modulo operation on the DMRS, a reception gain can be obtained.
  • Patent Document 1 the method of arranging DMRSs addressed to each terminal device in orthogonal radio resources requires as many DMRS dedicated radio resources as the number of terminal devices. There was a problem that the increase became large.
  • the present invention has been made in view of such circumstances, and in the NLP MU-MIMO system, the DMRS is spatially multiplexed and the increase in overhead due to DMRS insertion is minimized, while the data signal amplitude and
  • An object is to provide a technique capable of estimating a phase (complex gain / reception gain).
  • a transmitting device that transmits data signals to a plurality of receiving devices at the same time and the same frequency, the demodulation reference signal addressed to the first receiving device, and the first
  • a transmission device comprising: a transmission unit that transmits a demodulation reference signal addressed to a second reception device different from the reception device at the same time and the same frequency.
  • the spatially multiplexed demodulation reference signal can transmit the demodulation reference signal to all the mobile station devices with one radio resource, so that sufficient radio resources are allocated to arrange data signals.
  • a plurality of demodulation reference signals can be transmitted to all mobile station apparatuses.
  • Each mobile station apparatus (the following receiving apparatus) can estimate a propagation path using a plurality of demodulation reference signals. Therefore, if the present invention is used, the insertion loss of the demodulation reference signal can be greatly reduced.
  • non-linear precoding unit that adds a signal that is an integral multiple of a predetermined signal to the demodulation reference signal. Further, a non-linear precoding unit that performs non-linear precoding on the demodulation reference signal may be provided.
  • the transmission unit transmits a plurality of data signals at the same time and the same frequency
  • the nonlinear precoding unit performs the same nonlinear precoding processing as the data signal on the demodulation reference signal.
  • nonlinear precoding can be performed on the same principle using the same filter as the data signal.
  • a DMRS correction unit that corrects the demodulation reference signal.
  • the DMRS correction unit preferably includes a two-dimensional Euclidean algorithm unit that applies a two-dimensional Euclidean algorithm to the demodulation reference signal.
  • the two-dimensional Euclidean algorithm section includes a signal obtained by rotating the phase of the second demodulation reference signal and the second demodulation reference signal by 90 degrees from the first demodulation reference signal, The signal obtained by rotating the phase of the demodulation reference signal by 180 degrees and the signal obtained by rotating the phase of the second demodulation reference signal by 270 degrees may be subtracted to have a difference vector calculation unit.
  • the present invention is a receiving device characterized by having a perturbation vector addition unit for adding a signal of an integral multiple of a predetermined width to a demodulation reference signal. It is preferable to have a provisional channel estimation unit that estimates a channel using a demodulation reference signal obtained by adding the signals.
  • a perturbation vector candidate selection unit that selects a plurality of different signals, and a perturbation vector estimation that selects one of the signals based on a channel estimation result obtained by channel estimation using the temporary channel estimation unit using each of the signals. It is good to have a part.
  • a demodulator that softly estimates a data signal based on each of the propagation path estimation results corresponding to a plurality of different signals and calculates a log likelihood ratio, respectively, and a variance of each log likelihood ratio It is preferable that a perturbation vector evaluation value calculation unit to calculate is selected, and the perturbation vector estimation unit selects the signal corresponding to the largest of the variances. It is preferable that a two-dimensional Euclidean algorithm for calculating an irreducible vector is provided by applying a two-dimensional Euclidean algorithm to the plurality of demodulation reference signals. It is preferable to have a complex gain calculation unit that calculates a complex gain of the propagation path using the irreducible vector.
  • the present invention also provides a nonlinear precoding unit that adds a signal that is an integer multiple of a predetermined signal to a demodulation reference signal, and the demodulation reference signal and other signals at the same time and the same frequency.
  • a communication system comprising: a transmitter having a transmitter for transmitting; and a receiver having a perturbation vector adder for adding a signal having an integer multiple of a predetermined width to a demodulation reference signal. is there.
  • the present invention also includes a step of adding a signal that is an integral multiple of a predetermined signal to the demodulation reference signal, and a step of transmitting the demodulation reference signal and another signal at the same time and the same frequency.
  • a reception method including a step of adding a signal having an integer multiple of a predetermined width to a demodulation reference signal.
  • the present invention also provides a nonlinear precoding unit that adds a signal that is an integer multiple of a predetermined signal to a demodulation reference signal, and the demodulation reference signal and other signals at the same time and the same frequency.
  • An integrated circuit having a transmitting step.
  • the present invention is an integrated circuit characterized by having a perturbation vector addition unit for adding a signal having an integer multiple of a predetermined width to a demodulation reference signal.
  • the present invention may be a program for causing a content to execute the communication method described above, or a computer-readable recording medium for recording the program.
  • FIG. 1 is a conceptual diagram showing a configuration example of a communication system 1 according to the first embodiment of the present invention.
  • the base station apparatus A1 transmits a common reference signal (Common Reference Signal: CRS).
  • CRS is a signal in which the base station apparatus A1 and the first to fourth terminal apparatuses B11 to B14 store their reference signals in advance.
  • Each of the first to fourth terminal apparatuses B11 to B14 estimates the propagation path state based on the CRS transmitted by the base station apparatus A1, and transmits the propagation path state information based on the estimated propagation path state. Notify A1.
  • the base station apparatus A1 transmits DMRS and data signals to the first to fourth terminal apparatuses B11 to B14.
  • base station apparatus A1 performs precoding on the DMRS and data signal, and transmits the multiplied DMRS and data signal.
  • the multiplexed first to fourth terminal apparatuses B11 to B14 are based on the DMRS received from the base station apparatus A1, and the equivalent propagation path (hereinafter simply referred to as equivalent propagation) in which the precoding process is regarded as a part of the propagation path. And a data signal is obtained based on equivalent channel state information indicating the channel state of the estimated equivalent channel.
  • FIG. 2 is a sequence diagram showing an example of the operation of the communication system 1 according to the present embodiment. This figure is a diagram showing an example of the operation of the communication system 1 in the case of FIG.
  • Step S101 The base station apparatus A1 transmits the CRS to the first to fourth terminal apparatuses B11 to B14. Thereafter, the process proceeds to step S102.
  • Step S102 The first to fourth terminal apparatuses B11 to B14 estimate the propagation path state based on the CRS transmitted in step S101. Thereafter, the process proceeds to step S103.
  • Step S103 The first to fourth terminal apparatuses B11 to B14 calculate channel state information based on the channel state estimated in step S102. Thereafter, the process proceeds to step S104.
  • Step S104 The first to fourth terminal apparatuses B11 to B14 notify the base station apparatus A1 of the propagation path state information calculated in step S103. Thereafter, the process proceeds to step S105.
  • Step S105 Base station apparatus A1 calculates the filter used for nonlinear precoding based on the propagation path state information notified by step S104.
  • the base station apparatus A1 performs nonlinear precoding on the generated DMRS and data signal using the filter, and generates the DMRS and data signal. Thereafter, the process proceeds to step S106.
  • Step S106 The base station apparatus A1 transmits the DMRS signal generated in step S105 to the first to fourth terminal apparatuses B11, B12, B13, and B14. Thereafter, the process proceeds to step S107.
  • Step S107 The first to fourth terminal apparatuses B11 to B14 estimate the propagation path state of the equivalent propagation path based on the DMRS signal transmitted in step S106. Thereafter, the process proceeds to step S108.
  • Step S108 The base station apparatus A1 transmits the data signal generated in step S105 to the terminal apparatuses B11 to B14. Thereafter, the process proceeds to step S109.
  • Step S109 The first to fourth terminal apparatuses B11 to B14 detect and acquire a data signal based on the equivalent channel state information indicating the channel state of the equivalent channel estimated in step S108.
  • terminal apparatuses B11 to B14 are multiplexed.
  • N first to Nth terminal apparatuses B11 to B1N are multiplexed thereafter.
  • the terminal devices B11 to B1N are collectively referred to as a terminal device B1n.
  • FIG. 3 is a functional block diagram illustrating a schematic configuration example of the base station apparatus A1 according to the present embodiment.
  • the base station apparatus A1 includes first to Nth antennas a101-1 to a101-N, first to Nth receiving units a102-1 to a102-N, first to Nth antennas.
  • GI Guard Interval removal units a103-1 to a103-N, first to Nth FFT (Fast Fourier Transform) units a104-1 to a104-N, propagation path state information acquisition units a105, filter calculation unit a11, first to Nth code units a121-1 to a121-N, first to Nth modulation units a122-1 to a122-N, DMRS generation unit a124, eigensignal configuration unit a125, nonlinear precoding unit a13, CRS generation unit a141, frame configuration unit a142, first to Nth I FT (Inverse Fast Fourier Transform) sections a143-1 to a143-N, first to Nth GI insertion sections a144-1 to a144-N, and first to Nth transmission sections a145-1 to a145-N.
  • I FT Inverse Fast Fourier Transform
  • N 4
  • OFDM Orthogonal Frequency Division Multiplexing
  • the device A1 may use a time division multiplexing (TDM) method or a frequency division multiplexing (FDM) method in one or both of the uplink and the downlink.
  • a signal (carrier frequency signal) is received. This signal includes propagation path state information.
  • the first to Nth receivers a102-n downconvert the received signal and perform A / D (analog / digital) conversion to generate a baseband digital signal.
  • the first to Nth receiving units a102-n output the generated digital signal to the first to Nth GI removing units a103-n.
  • the first to Nth GI removal units a103-n remove the GI from the digital signals input from the first to Nth reception units a102-n, and the removed signals from the first to the Nth To the FFT unit a104-n.
  • the first to N-th FFT units a104-n perform frequency-domain signals by performing FFT on the signals input from the first to N-th GI removal units a103-n.
  • the first to N-th FFT units a104-n output the generated frequency domain signals to the propagation path state information acquisition unit a105.
  • the propagation path state information acquisition unit a105 demodulates signals input from the first to N-th FFT units a104-n, and extracts propagation path state information from the demodulated information.
  • the propagation path state information acquisition unit a105 outputs the extracted propagation path state information to the filter calculation unit a11.
  • signals other than the signal of the propagation path state information are demodulated by a control unit (not shown).
  • control information is used to control the base station apparatus A1, and data other than the control information is transmitted to other base station apparatuses, server apparatuses, and the like.
  • the filter calculation unit a11 calculates a filter used for nonlinear precoding based on the propagation path state information input from the propagation path state information acquisition unit a105. Details of the filter calculation processing performed by the filter calculation unit a11 will be described later.
  • the filter calculation unit a11 inputs the calculated filter to the nonlinear precoding unit a13.
  • the first to N-th code units a121-n perform error correction coding on the input information bits and output the encoded bits to the first to N-th modulation units a122-n.
  • the first to Nth modulation units a122-n modulate the encoded bits input from the first to Nth code units a121-n to generate a data signal addressed to the terminal device B1n.
  • the first to Nth modulation units a122-n output the generated data signal to the unique signal configuration unit a125.
  • the base station apparatus A1 determines modulation schemes to be used in the first to N-th modulation units a122-n based on the propagation path state information, and transmits the modulation information indicating the determined modulation schemes to the modulation units a122-n. To the first to N-th terminal devices B1n from the base station device A1.
  • the DMRS generator a124 generates DMRSs addressed to the first to Nth terminal devices B1n.
  • the DMRS generation unit a124 outputs the generated DMRS to the specific signal configuration unit a125.
  • the unique signal configuration unit a125 receives the data signals addressed to the first to Nth terminal apparatuses B1n input from the first to Nth modulation units 122-n and the first to Nth data signals input from the DMRS generation unit a124.
  • the DMRSs addressed to the terminal devices B1n up to the Nth are related.
  • Each of the information addressed to the first to Nth terminal apparatuses B1n related by the unique signal configuration unit a125 is referred to as a unique signal of the first to Nth terminal apparatuses B1n.
  • the unique signal configuration unit a125 outputs each of the unique signals of the first to Nth terminal apparatuses B1n generated by the association to the nonlinear precoding unit a13.
  • the non-linear precoding unit a13 performs non-linear precoding on the eigensignals (data signal and DMRS) of the first to Nth terminal apparatuses B1n input from the eigensignal configuration unit a125. Details of the nonlinear precoding performed by the nonlinear precoding unit a13 will be described later.
  • the nonlinear precoding unit a13 outputs the unique signal subjected to the nonlinear precoding to the frame configuration unit a142.
  • the CRS generation unit a141 generates a CRS having a known reference signal in the base station apparatus A1 and the first to Nth terminal apparatuses B1n, and outputs the generated CRS to the frame configuration unit a142.
  • the frame configuration unit a142 maps the unique signal input from the nonlinear precoding unit a13 and the CRS input from the CRS generation unit a141. Note that the frame configuration unit a142 may map the unique signal and the CRS to different frames, or may map to the same frame. For example, only the CRS may be mapped to a certain frame, and the CRS and unique signal may be mapped to another frame. It is assumed that the base station apparatus A1 maps the CRS and the unique signal to the frame according to a predetermined mapping, and the first to Nth terminal apparatuses B1n know the mapping in advance.
  • the frame configuration unit a142 outputs the signal transmitted by the antenna a101-n among the mapped signals to the first to Nth IFFT units a143-n in units of frames.
  • the 1st to N-th IFFT units a143-n generate time-domain signals by performing IFFT on the signals input from the frame configuration unit a142.
  • the first to Nth IFFT units a143-n output the generated time domain signals to the first to Nth GI insertion units a144-n.
  • the first to Nth GI insertion sections a144-n add a guard interval to the signals input from the first to Nth IFFT sections a143-n,
  • the data is output to the Nth transmitters a145-n.
  • the first to Nth transmitters a145-n D / A convert the signals (baseband digital signals) input from the first to Nth GI inserters a144-n. .
  • the first to Nth transmitters a145-n upconvert the converted signal to generate a carrier frequency signal.
  • the first to Nth transmitters a145-n transmit the generated signals via the first to Nth antennas a101-n.
  • FIG. 4 is a schematic diagram illustrating an example of the configuration of the unique signal according to the present embodiment.
  • the horizontal axis represents time.
  • This figure represents the unique signal output by the unique signal configuration unit a125. Further, this figure shows the arrangement of unique signals of the first to Nth terminal apparatuses B1n transmitted at the same frequency, with the time axis aligned.
  • DMRSs addressed to the first to Nth terminal apparatuses B1n are represented by “DMRS-MSn”.
  • the unique signal is a signal for each of the first to Nth terminal apparatuses B1n and is composed of a DMRS and a data signal.
  • the signal with the code S11 indicates DMRS (DMRS-MS1) addressed to the terminal device B11.
  • FIG. 4 shows that the DMRSs addressed to the first to Nth terminal apparatuses B1n are all output from the unique signal configuration unit a125 at the same time.
  • This figure also shows that data signals addressed to the first to N-th terminal devices B1n are also output from the unique signal configuration unit a125 at the same time.
  • both data signals and DMRS assigned at the same time are spatially multiplexed by nonlinear precoding a13, and are transmitted from base station apparatus A1 at the same time and the same frequency (that is, the same subcarrier of the same OFDM symbol).
  • Sent by FIG. 4 shows that the DMRS and the data signal are output from the unique signal configuration unit a125 at different times. For example, DMRS addressed to the mobile station B1n The first through N are output to the time t 1 and time t 2, the data signal is output to t 3 ⁇ t L.
  • the configuration of the unique signal in FIG. 4 is an example, and the present invention is not limited to this.
  • specific signal configuration unit a125 after the DMRS addressed to the mobile station B1n is output (t 3 or later) on it data signal is output, specific signal configuration unit a125 also outputs the DMRS after outputting the data signals Good.
  • the unique signal configuration unit a125 may alternately output the data signal and the DMRS on the time axis, or may output in another order.
  • FIG. 5A is a schematic diagram showing an example of a frame configuration according to the present embodiment. This figure shows the structure of a frame in which signals are mapped by the frame configuration unit a142. This figure also shows the structure of frames transmitted from the first to N-th antennas a101-n with the time axis aligned.
  • the CRS transmitted from the first to Nth antennas a101-n is represented by “CRS-Txn”.
  • CRS-Txn the CRS transmitted from the first to Nth antennas a101-n is represented by “CRS-Txn”.
  • each of the first to Nth antennas a101-n is composed of a CRS and a unique signal (a unique signal after nonlinear precoding composed of a data signal and DMRS).
  • the signal denoted by reference numeral S21 indicates CRS (CRS-Tx1) transmitted by the antenna a101-1.
  • a signal denoted by reference sign S22 is a unique signal transmitted by the antenna a 101-1, and indicates a unique signal after nonlinear precoding.
  • FIG. 5A shows that CRSs transmitted by antenna a101-n are arranged in different time bands by frame configuration unit a142 and transmitted from base station apparatus A1. For example, CRS sent by the antenna a101-1 is sent to the time t 1, CRS sent by the antenna a101-2 is transmitted to the time t 2.
  • FIG. 5A shows that the CRS and the specific signal are arranged in different time bands in the frame configuration unit a142.
  • the CRS transmitted by the antenna a101-n is arranged in a time band transmitted from time t 1 to t N
  • the unique signal is arranged in a band after t N + 1 .
  • This figure also shows that all the unique signals transmitted by the first to Nth antennas a101-n are arranged in the same time band by the frame configuration unit a142.
  • the frame configuration in FIG. 5A is an example, and the present invention is not limited to this.
  • all CRS transmitted by the antenna a101-n are arranged ahead of time bands (before t N + 1), after the time band (t N + 1 and later) are disposed unique signal
  • the frame configuration unit a142 may arrange the unique signal in the previous time band and arrange the CRS in the later time band.
  • the frame configuration unit a142 may be arranged in the time band in which the unique signal and the CRS alternate, or may be arranged in another order.
  • the frame configuration unit a142 may arrange the CRS transmitted by the antenna a 101-1 in the time band after the time tN + 1 and arrange the unique signal in the time band before the time tN + 1 .
  • base station apparatus A1 transmits the frame which has arrange
  • FIG. 5B and FIG. 5C show a configuration example of a frame transmitted by each of the first to N-th antennas a101-n.
  • 4 and 5A show examples in which DMRS and data signals, or CRS and unique signals (DMRS and data signals are included) are arranged in the time direction, but as shown in FIGS. 5B and 5C,
  • the CRS, DMRS, and data signal may be arranged on a two-dimensional matrix in the time direction (t) and the frequency direction (f).
  • FIG. 6 is a functional block diagram illustrating a configuration example of the terminal device B1n according to the present embodiment.
  • the terminal device B1n includes an antenna b101, a receiving unit b102, a GI removing unit b103, an FFT unit b104, a signal separating unit b105, a CRS channel estimating unit b107, a DMRS channel estimating unit b12, and a channel compensating unit.
  • b106 a modulo calculation unit b109, a demodulation unit b110, a decoding unit b111, a propagation path state information generation unit b108, an IFFT unit b131, a GI insertion unit b132, and a transmission unit b133.
  • the receiving unit b102 receives a signal (carrier frequency signal) transmitted from each terminal device B1n via the antenna b101.
  • the receiving unit b102 down-converts the received signal and performs A / D (analog / digital) conversion to generate a baseband digital signal.
  • the receiving unit b102 outputs the generated digital signal to the GI removing unit b103.
  • the GI removal unit b103 removes the GI from the digital signal input from the reception unit b102, and outputs the signal after the removal to the FFT unit b104.
  • the FFT unit b104 generates a frequency domain signal by performing a fast Fourier transform on the signal input from the GI removing unit b103.
  • the FFT unit b104 outputs the generated frequency domain signal to the signal separation unit b105.
  • the signal separation unit b105 demaps the signal input from the FFT unit b104 based on the mapping information notified from the base station apparatus A1. Of the demapped signals, the signal separator b105 outputs CRS to the CRS channel estimator b107, and outputs DMRS to the DMRS channel estimator b12. The signal separator b105 outputs the data signal to the propagation path compensator b106. Further, the signal separator b105 inputs the data signal to the DMRS propagation path estimator b12.
  • the CRS channel estimator b107 estimates the channel state based on the CRS input from the signal separator b105, and outputs information indicating the estimated channel state to the channel state information generator b108.
  • the DMRS propagation path estimation unit b12 estimates a propagation path state of an equivalent propagation path in which a filter used for nonlinear precoding is regarded as a part of the propagation path based on the DMRS input from the signal separation unit b105. Details of the DMRS propagation path estimation unit b12 will be described later.
  • the DMRS channel estimator b12 outputs equivalent channel state information indicating the estimated channel state of the equivalent channel to the channel compensator b106.
  • the propagation path compensation unit b106 performs propagation path compensation on the signal input from the signal separation unit b105 using the equivalent propagation path state information input from the DMRS propagation path estimation unit b12.
  • the propagation path compensation unit b106 outputs the signal after propagation path compensation to the modulo calculation unit b109.
  • the modulo calculator b109 performs a modulo calculation on the signal input from the propagation path compensator b106 based on the modulation information notified from the base station apparatus A1.
  • the modulo operation can be expressed by the following equation.
  • Equation (1-1) ⁇ -floor ((Re [ ⁇ ] + ⁇ / 2) / ⁇ ) ⁇ -i ⁇ floor ((Im [ ⁇ ] + ⁇ / 2) / ⁇ ) ⁇ (1-1)
  • floor (x) represents a maximum integer not exceeding x
  • Re [ ⁇ ] and Im [ ⁇ ] represent a real part and an imaginary part of the complex number ⁇ , respectively.
  • i is an imaginary unit.
  • ⁇ represents the modulo width.
  • ⁇ may be 2 ⁇ 2 times the average amplitude of the QPSK signal, 8 / ⁇ 10 times that of the 16QAM signal if 16QAM, and 16 / ⁇ 42 times that of the 64QAM signal if 64QAM. desirable.
  • the modulo width may be another value.
  • the modulo calculation unit b109 outputs the signal after the modulo calculation to the demodulation unit b110.
  • the demodulator b110 demodulates the signal input from the modulo calculator b109 in the modulation scheme indicated by the modulation information notified from the base station apparatus A1.
  • the demodulator b110 outputs the demodulated information (the hard-decided encoded bit or the soft estimated value of the encoded bit) to the decoder b111.
  • the decoding unit b111 acquires information bits by decoding the information input from the demodulation unit b110, and outputs the acquired information bits.
  • the propagation path state information generation unit b108 generates propagation path state information from the propagation path state input from the CRS propagation path estimation unit b107 (referred to as propagation path state information generation processing).
  • h n1 , h n2 ,..., h nN are propagation path states between the antenna a 101-1, the antenna a 101-2,.
  • the propagation path state information generation unit b108 does not necessarily need to notify the base station apparatus A1 using the row vector
  • may be notified to the base station apparatus A1 as the propagation path state information of the row vector C ⁇ h n /
  • Information approximate to a predetermined value may be notified to the base station apparatus A1 as propagation path state information.
  • the propagation path state information generation unit b108 modulates the generated propagation path state information and outputs the modulated propagation path state information signal to the IFFT unit b131.
  • the IFFT unit b131 generates a time domain signal by performing inverse fast Fourier transform on the signal input from the propagation path state information generation unit b108.
  • the IFFT unit b131 outputs the generated time domain signal to the GI insertion unit b132.
  • the GI insertion unit b132 adds a guard interval to the signal input from the IFFT unit b131, and outputs the signal after the addition to the transmission unit b133.
  • the transmission unit b133 D / A converts the signal (baseband digital signal) input from the GI insertion unit b132.
  • the transmitter b133 generates a carrier frequency signal by up-converting the converted signal.
  • the transmitter b133 transmits the generated signal via the antenna b101.
  • the filter calculation unit a11 illustrated in FIG. 3 configures the propagation path matrix H from the propagation path state information input from the propagation path state information acquisition unit a105.
  • H is a matrix of N rows and N columns, and the p row and q column components indicate the complex gain of the propagation path between the p th terminal device B1p and the q th antenna a101-q of the base station device (here, p and q are arbitrary integers from 1 to N).
  • the propagation path state acquired from the propagation path state information notified from each terminal apparatus B1p becomes a row vector, and a propagation path matrix H can be generated by generating a matrix having a row vector corresponding to all terminal apparatuses in each line. . Even if the norm of the propagation path state information is normalized in the terminal apparatus B1n, the base station apparatus A1 keeps the propagation path state information notified from the terminal apparatus B1n as it is, that is, the normalized propagation path state in each row. To generate a channel matrix H.
  • each unique signal input from the unique signal component a125 is denoted by sn, and a vertical vector having all s1 to sN as components is denoted by s.
  • the unique signal indicates a data signal and the DMRS generated by the DMRS generator a124.
  • the nonlinear precoding unit a13 searches for a combination of the N-dimensional integer vertical vectors z1 and z2 that minimizes the norm of the transmission signal multiplied by the filter W. This is expressed as follows.
  • (Z1, Z2) argmin (z1, z2)
  • (Z1, Z2) be a combination of (z1, z2) that minimizes the norm
  • L BS an integer: [- L BS, -L BS +1, ..., - 2, -1,0 , 1, 2,..., L BS -1, L BS ]
  • z1 ⁇ + iz2 ⁇ is called a perturbation vector.
  • the base station apparatus A1 must normalize the total transmission power of data signals in a certain number of subcarriers and a certain number of OFDM symbols (referred to as “power normalization units”) in order to keep the transmission power constant. Don't be.
  • the power normalization unit indicates, for example, the entire frame unit shown in FIGS. 5B and 5C.
  • the nonlinear precoding unit a13 multiplies the eigensignal x (data signal and DMRS) by the reciprocal of the power normalization coefficient g, and inputs the multiplied data signal to the frame configuration unit a142.
  • each component of the signal x indicates a transmission signal transmitted from each of the first to N-th antennas a101-1 to a101-N in order.
  • the propagation path estimation unit b12 estimates the propagation path, and the propagation path compensation unit b106 multiplies the amplitude with respect to the data signal by g. Can be compensated. That is, a signal can be correctly detected by multiplying the received data signal by g.
  • the nonlinear precoding unit a13 applies the same nonlinear precoding to the DMRS as the data signal, including normalization of power.
  • FIG. 7 shows a detailed configuration of the DMRS propagation path estimation unit b12.
  • the DMRS propagation channel estimation unit b12 includes a perturbation vector candidate selection unit b121, a perturbation vector addition unit b122, a temporary propagation channel estimation unit b123, a propagation channel compensation unit b124, a modulo operation unit b125, a demodulation unit b126, and a perturbation vector evaluation value calculation unit b127. , And a perturbation vector estimation unit b128.
  • the perturbation vector candidate selection unit b121 selects a candidate perturbation vector candidate, and inputs the selected perturbation vector to the perturbation vector addition unit b122.
  • a modulo width for DMRS common to the base station apparatus A1 and the terminal apparatus B1n is used.
  • L BS L MS
  • L BS L MS
  • the temporary propagation path estimation unit b123 inputs the amplitude
  • the propagation path compensation unit b124 uses the amplitude and phase input from the temporary propagation path estimation unit b123 to perform propagation path compensation on the data signal in the frame in which the DMRS is arranged. That is, the received signal is divided by the complex gain h.
  • the data signal after propagation path compensation is input to the modulo arithmetic unit b125.
  • Step S25 The modulo operation unit b125 performs a modulo operation on the data signal subjected to the propagation path compensation by the propagation path compensation unit b124.
  • a modulo width for data signals common to the base station apparatus A1 and the terminal apparatus B1n is used.
  • the modulo operation unit b125 inputs the data signal after the modulo operation to the demodulation unit b126.
  • Step S26 The demodulator b126 performs soft estimation on the input data signal after the modulo calculation, and inputs the soft estimation value to the perturbation vector evaluation value calculation unit b127.
  • LLR log likelihood ratio
  • the demodulator b126 inputs the LLR corresponding to each bit assigned to each data signal to the perturbation vector evaluation value calculator b127 based on the equations (1-4) and (1-5). That is, the number of data signals ⁇ 2 for QPSK and the number of data signals ⁇ 4 for 16QAM are calculated and input to the perturbation vector evaluation value calculation unit b127.
  • the perturbation vector evaluation value calculation unit b127 calculates the variance of the input LLR, and inputs the calculated LLR variance and the corresponding LLR to the perturbation vector estimation unit b128.
  • the LLR variance is calculated by the following equation.
  • V is the number of data signals soft-estimated for each perturbation vector candidate.
  • M is the number of bits allocated to each modulation system, and is 2 for QPSK and 4 for 16QAM.
  • L v m represents the LLR of the mth bit assigned to the vth data signal.
  • Step S28 The perturbation vector estimation unit b128 selects the largest of the LLR variances corresponding to each perturbation vector candidate input from the perturbation vector evaluation value calculation unit b127, and sets the perturbation vector having the maximum variance.
  • the amplitude and phase of the corresponding data signal (referred to as “equivalent channel state information”) are output and input to the channel compensator b106 (FIG. 6) outside the DMRS channel estimator b12.
  • a large variance of the LLR indicates that the mutual information amount transmitted from the base station apparatus A1 is the largest when each perturbation vector candidate is assumed.
  • the LLR has a larger absolute value as the “probability / probability” of whether each bit is 1 or 0 is higher. Therefore, assuming each perturbation vector candidate, by calculating the variance of each LLR and selecting the perturbation vector having the largest LLR variance, the perturbation vector is the perturbation vector having the highest “probability / probability”. It can be estimated that there is.
  • the propagation path compensation unit b124, the modulo calculation unit b125, and the demodulation unit b126 in the DMRS propagation path estimation unit b12 are the propagation path compensation unit b106, the modulo calculation unit b109, and the demodulation unit b110 outside the DMRS propagation path estimation unit b12. Therefore, both components may be made a common circuit in order to reduce the circuit scale.
  • the soft estimation value corresponding to the perturbation vector estimated by the perturbation vector estimation unit b128 among the soft estimation values once calculated by the demodulation unit b126 is used by the decoding unit b111, so that the soft estimation value is calculated twice. The calculation amount may be reduced by avoiding the above.
  • the base station apparatus spatially multiplexes and transmits DMRS using NLP MU-MIMO, so that overhead due to DMRS insertion can be reduced with high power efficiency.
  • the data signal that is softly estimated for estimating the perturbation vector may be a data signal that is a part of the frame.
  • the perturbation vector may be estimated using an index other than the LLR variance.
  • Second Embodiment Two-dimensional Euclidean Mutation Method
  • the terminal apparatus B1n estimates the most likely perturbation vector from the perturbation vector candidates.
  • DMRS can be spatially multiplexed by NLP MU-MIMO.
  • both the base station apparatus A2 according to the present embodiment and the terminal apparatus B2n according to the present embodiment introduce a process called “two-dimensional Euclidean mutual division method”, so that the computation amount of the terminal apparatus B2n is increased. Is reduced as compared with the first embodiment.
  • the base station apparatus A2 according to the present embodiment has a DMRS correction unit a226.
  • the base station apparatus A1 of the first embodiment except for the operation of the DMRS propagation path estimation unit b22, the base station apparatus A1 of the first embodiment
  • the terminal device B1n has the same configuration as the terminal device B1n.
  • FIG. 15 shows the configuration of the base station apparatus A2. As described above, a DMRS correction unit a226 is newly added to FIG.
  • the DMRS correction unit a226 of the base station apparatus A2 and the DMRS propagation path estimation unit b22 of the terminal apparatus B2n will be described in detail.
  • Both the base station apparatus A2 and the terminal apparatus B2n according to the present embodiment have a two-dimensional Euclidean algorithm division unit that performs a process called "two-dimensional Euclidean algorithm". Since the two-dimensional Euclidean algorithm section is one of the characteristic parts of this embodiment, the principle will be described first.
  • the 2D Euclidean algorithm is an algorithm that extends the normal Euclidean algorithm to 2D.
  • the normal Euclidean algorithm is called a one-dimensional Euclidean algorithm.
  • An extension to complex numbers is the two-dimensional Euclidean algorithm.
  • the method is to calculate a “base vector (or irreducible vector)” corresponding to the greatest common divisor of both Gaussian integers when there are two complex numbers (Gaussian integers) whose real part and imaginary part are integers. It is.
  • a base vector or irreducible vector
  • a Gaussian integer having a large norm (absolute value) of 3 + i and -1 + i and a small Gaussian integer multiplied by +1, ⁇ 1, + i, and ⁇ i are added (Gauss integer is multiplied by ⁇ 1, + i, and ⁇ (Multiplying by i corresponds to rotating the phase of a Gaussian integer by +180 degrees, +90 degrees, and +270 degrees, respectively).
  • This embodiment is characterized in that DMRS propagation path estimation is performed using the two-dimensional Euclidean mutual division method as described above. The principle of this propagation path estimation will be described. In the present embodiment, it is assumed that the reference DMRS is always one of the points of the QPSK signal.
  • FIG. 8 shows an example of the positions of two DMRSs (referred to as DMRS1 and DMRS2) on the signal point plane.
  • the base station apparatus transmits DMRS1 at a position of 3 + i and DMRS2 transmits at a position of -1 + i.
  • the unit of this signal point plane is set so that the amplitude of the QPSK signal is ⁇ 2. That is, each signal point of QPSK is a point of ( ⁇ 1 ⁇ i), and the modulo width is 4 (for convenience, the unit of amplitude is changed from that of the first embodiment).
  • the DMRS1 and DMRS2 and ⁇ other points are all QPSK signal points and points obtained by adding arbitrary perturbation vectors to the QPSK signal points. As shown in FIG. The dots are aligned on the top.
  • DMRS1 is a signal obtained by adding the perturbation vector +4 to -1 + i of the QPSK signal.
  • the complex number 1 + i of “the greatest common divisor” can be obtained.
  • This irreducible vector corresponds to one of the four points of QPSK.
  • the irreducible vector is also multiplied by a to be a (1 + i) An irreducible vector can be obtained.
  • DMRS1 and DMRS2 are h (3 + i) and h ( ⁇ 1 + i) multiplied by the complex gain h of the propagation path. Noise is ignored here.
  • an irreducible vector h (1 + i) is obtained. This is a signal obtained by multiplying one of the four points of the QPSK signal to which the perturbation vector is not added by the complex gain of the propagation path. Since it is known that each point of QPSK has a norm of ⁇ 2, the absolute value
  • the terminal device cannot determine whether the signal h (1 + i) obtained by the two-dimensional Euclidean mutual division method at the time of DMRS channel estimation is multiplied by h of 4 points ( ⁇ 1 ⁇ i) of QPSK. Assuming (1 + i), the phase arg (h) is obtained. Even if the other three points are ( ⁇ 1 + i), (1-i), and ( ⁇ 1 ⁇ i), the phase difference of the actual complex gain is always +90 degrees, +180 degrees, or +270 degrees. Become.
  • a modulo operation is performed on DMRS1 and DMRS2.
  • the modulo operation is the same processing for both I-ch and Q-ch, and is symmetric about each axis, and is symmetric about rotation of 90 degrees. Therefore, even if arg (h) is rotated by 90 degrees ⁇ integer multiple, there is no problem and modulo arithmetic can be performed on DMRS1 and DMRS2.
  • DMRS1 after the modulo calculation is h (1 + i)
  • DMRS2 is h (-1 + i).
  • DMRS1 is h (1 + i) + ⁇ 1
  • DMRS2 is h ( ⁇ 1 + i) + ⁇ 2.
  • h ( ⁇ 1 ⁇ i) corresponding to the irreducible vector may not be calculated by the two-dimensional Euclidean algorithm. is there.
  • the irreducible vector is 3 + 3i and does not coincide with the QPSK signal point.
  • the base station apparatus also uses the two-dimensional Euclidean algorithm to determine that the case shown in FIG. 9 occurs, and transmits two points with no problem by changing the signal point. To do.
  • 10 includes a vector storage unit 301, a difference vector calculation unit 303, a difference vector norm calculation unit 305, a convergence determination unit 307, and a norm calculation unit 309.
  • Step S31 First, the vector storage unit 301 inputs the two input DMRSs into the norm calculation unit 309.
  • the DMRS transmitted from the base station apparatus is multiplied by the complex gain h of the propagation path, and noise is added to the two DMRSs.
  • the norm calculation unit 309 calculates the norms of the two input DMRSs, and inputs each norm to the vector storage unit 301.
  • Step S33 The vector storage unit 301 inputs the norm corresponding to the two input DMRSs to the difference vector calculation unit 303.
  • the difference vector calculation unit 303 has a norm of two input vectors (hereinafter, a signal recursively calculated by the difference vector calculation unit 303 with respect to two DMRSs is referred to as a “vector”).
  • a vector a signal recursively calculated by the difference vector calculation unit 303 with respect to two DMRSs.
  • Step S35 The difference vector norm calculation unit 305 calculates the norm of the four input vectors c1 to c4, inputs the vector having the smallest norm and the norm corresponding to the vector to the vector storage unit 301, and sets the most norm The norm of the vector having a small value is input to the convergence determination unit 307.
  • Step S ⁇ b> 36 The convergence determination unit 307 compares the magnitude relationship between the input norm and a predetermined positive number T. When the norm is larger than T, the input norm and vector are input to the vector storage unit 301. When the norm is smaller than T, information indicating that the two-dimensional Euclidean algorithm is completed is input to the vector storage unit 301.
  • the convergence condition is not set when the norm becomes 0, because the DMRS includes noise on the terminal device side. This is because the norm of the portion 307 does not become zero. However, since noise is usually smaller than DMRS, it is determined that convergence has occurred when the norm becomes smaller than T using a predetermined constant T.
  • T The value of T is determined in advance by computer simulation in advance by a trade-off between the following two factors.
  • Step S37 The vector storage unit 301 1) a norm corresponding to the vector input from the difference vector calculation unit 303, 2) Two norms corresponding to the vector input from the difference vector norm calculation unit 305 are input to the difference vector calculation unit 303, and the process is repeated again from step S34.
  • Step S ⁇ b> 38 When information indicating that the two-dimensional Euclidean algorithm has been completed is input from the convergence determination unit 307, the vector storage unit 301 outputs a vector b (irreducible vector).
  • FIG. 11 shows the configuration of the DMRS correction unit a226.
  • the DMRS correction unit a ⁇ b> 226 includes a two-dimensional Euclidean algorithm unit 300, an irreducible vector confirmation unit 320, and a perturbation vector addition unit 340.
  • This DMRS correction unit a226 aims to avoid in advance the above-described problem that the irreducible vector does not match the minimum lattice vector by the base station apparatus A2.
  • the DMRS correction unit a226 first acquires the DMRS to which the perturbation vector is added from the nonlinear precoding unit a13.
  • the perturbation vector is added to the DMRS acquired by the DMRS correction unit a226.
  • a signal (q + z1 ⁇ + iz2 ⁇ ) before adding the filter W is used.
  • two DMRSs are acquired at a time, and these two are paired in advance, and which two are paired is known by the base station apparatus A2 and the terminal apparatus B2n. is there.
  • the two-dimensional Euclidean algorithm unit 300 obtains the irreducible vector (denoted as “d”) obtained by applying the two-dimensional Euclidean algorithm to the two DMRSs and the nonlinear precoding unit a13. These two DMRSs are input to the irreducible vector confirmation unit 320.
  • the irreducible vector confirmation unit 320 uses the two DMRSs acquired from the nonlinear precoding unit a13. Is again input to the nonlinear precoding unit a13.
  • the irreducible vector confirmation unit 320 uses two DMRSs acquired from the nonlinear precoding unit a13. Both are input to the perturbation vector adder 340.
  • the perturbation vector addition unit 340 adds the perturbation vector ⁇ (or - ⁇ , i ⁇ , -i ⁇ ) to either one of the DMRSs acquired from the nonlinear precoding unit a13.
  • which perturbation vector ⁇ is added to which DMRS is randomly selected using a random number or the like.
  • the perturbation vector addition unit 340 inputs two DMRSs including the DMRS obtained by adding the perturbation vectors to the two-dimensional Euclidean algorithm unit 300 again.
  • the 2D Euclidean algorithm section 300 performs the 2D Euclidean algorithm again using the new DMRS. Thereafter, the irreducible vector confirmation unit 320 determines that the irreducible vector d is one ( ⁇ 1 ⁇ i) of DMRS before adding the perturbation vector, that is, QPSK signal ( ⁇ 1 ⁇ i). The processing of the irreducible vector confirmation unit 320 and the perturbation vector addition unit 340 is continued.
  • the irreducible vector confirmation unit 320 outputs two corrected DMRSs and inputs them to the nonlinear precoding unit a13.
  • the nonlinear precoding unit a13 multiplies the DMRS by the filter W without searching for a new perturbation vector, and reciprocal g ⁇ 1 of the power normalization coefficient g. Multiply
  • FIG. 12A shows the detailed configuration of the DMRS propagation path estimation unit b12.
  • the DMRS propagation path estimation unit b12 includes a two-dimensional Euclidean division unit 300 and a complex gain calculation unit 350.
  • the complex gain calculation unit 350 includes a provisional complex gain calculation unit 351, a DMRS propagation path compensation unit 352, a DMRS-Modulo unit 353, and a vector division unit 354, as shown in FIG.
  • the DMRS propagation path estimation unit b12 acquires two DMRSs received by the terminal device B2n. These two are the DMRSs (referred to here as p 1 and p 2 ) that are the pair described in the DMRS correction unit a 226 in the base station apparatus.
  • the two-dimensional Euclidean algorithm unit 300 performs a two-dimensional Euclidean algorithm on the two DMRSs to calculate an irreducible vector pred.
  • irreducible vector p red because the complex gain h is multiplied is different from the irreducible vector d computed by the DMRS corrector A226.
  • the irreducible vector pred is a signal obtained by multiplying one of the four points of the QPSK signal to which the perturbation vector is not added by the complex gain h of the propagation path and adding an error due to noise. is there.
  • ( pred / ⁇ 2) of the complex gain of the propagation path can also be obtained.
  • the provisional complex gain calculation unit 351 inputs the absolute value
  • of the complex gain of the propagation path and the phase ⁇ ( arg ( pred / (1 + i))) calculated here to the DMRS propagation path compensation unit 352.
  • the DMRS propagation path compensation unit 352 performs propagation path compensation for each DMRSp 1 and p 2 using the input absolute value
  • DMRS propagation path compensation section 352 inputs DMRS q 1 and q 2 after propagation path compensation, the absolute value
  • the DMRS-Modulo unit 353 uses the absolute value
  • of the complex gain input from the DMRS propagation path compensation unit 352 and the phase ⁇ for the signals e 1 and e 2 on which the modulo operation has been performed, to generate the signal p mod1 (
  • ⁇ e 1 ⁇ exp (2 ⁇ i ⁇ )) and p mod2 (
  • the vector division unit 354 converts the signals p mod1 and p mod2 into the DMRS reference signal q (that is, the DMRS before the perturbation vector is added in the base station apparatus) for which the terminal apparatus previously knows the DMRS received by the terminal apparatus.
  • the vector division unit 354 inputs the estimated complex gain h of the propagation path to the propagation path compensation unit b106.
  • the base station apparatus A2 and the terminal apparatus B2n perform processing based on the two-dimensional Euclidean mutual division method, so that DMRS is spatially multiplexed by NLP MU-MIMO with a relatively small amount of calculation of the terminal apparatus B2n. It becomes possible to do. Accordingly, it is possible to reduce overhead due to insertion of DMRS.
  • the base station apparatuses B1n and B2n according to the first embodiment and the second embodiment both perform VP using the nonlinear precoding unit a13 and the filter calculation unit a11.
  • This embodiment is a method using THP having a calculation amount lower than that of VP.
  • the base station apparatus according to the present embodiment is referred to as a base station apparatus A3
  • the terminal apparatuses according to the present embodiment are referred to as terminal apparatuses B31 to B3N
  • an arbitrary one is referred to as a terminal apparatus B3n.
  • the non-linear precoding unit according to the present embodiment is referred to as a non-linear precoding unit a33, and the filter calculation unit according to the present embodiment is referred to as a filter calculation unit a31.
  • the filter calculation unit a31 calculates the same propagation path matrix H as in the first embodiment, and then QR-decomposes H.
  • H H QR (3-1)
  • R is an upper triangular matrix and Q is a unitary matrix.
  • H be the conjugate transpose of the matrix.
  • A be a diagonal matrix obtained by extracting only the diagonal component of RH .
  • the linear filter P 0 is calculated by the following equation using A and Q.
  • the interference coefficient filter F is calculated by the following equation using A and R.
  • the filter calculation unit a31 inputs the linear filter P 0 the interference coefficient filter F to the non-linear precoding unit a33.
  • g trace ⁇ ((P 0 C v P 0 H ) / P tr ) 1/2 ⁇ (3-4)
  • g ⁇ trace ((P 0 C v P 0 H ) / P tr ) ⁇ 1/2
  • P s is the average power of the modulation signal
  • P v is the average power of the data signal after the modulo calculation
  • C v has diagonal components [P s , P v , P v ,. ] (3-5) Is a diagonal matrix.
  • the filter calculation unit a31 inputs a new matrix obtained by multiplying the linear filter P 0 calculated by the equation (3-2) by g ⁇ 1 as the linear filter P to the nonlinear precoding unit a33. At this time, power normalization as shown in the first embodiment is not performed.
  • the configuration of the non-linear precoding unit a33 is shown in FIG. 13, and a flowchart showing the operation thereof is shown in FIG.
  • operations in the nonlinear precoding unit a33 will be described in order.
  • Step S1 The interference calculation unit a131 and the linear filter multiplication unit a134 obtain the interference coefficient filter F and the linear filter P from the filter calculation unit a31, respectively.
  • Step S2 1 is substituted into the number n indicating the number of the terminal device that is calculating the transmission signal.
  • Step S3 The unique signal s1 addressed to the terminal device B31 is set to v1.
  • Step S5 The interference calculation unit a131 calculates the interference signal f2 received by the terminal device B32 using v1 using the following equation.
  • F (p, q) a 2 ⁇ 1 component of the matrix F.
  • Step S6 The interference subtraction unit a133-2 subtracts f2 from the unique signal s2 addressed to the terminal device B32 to calculate the signal s2-f2.
  • Step S7 The first modulo operation unit a132-2 applies the modulo operation to s2-f2 to calculate the signal v2.
  • Step S4 The value of n is increased by 1.
  • Step S5 The interference calculation unit a131 calculates the interference signal fn received by the terminal apparatus n using v1 to v (n-1) by the following equation.
  • F (n, 1: n ⁇ 1) * [v1, v2,. . . , V (n ⁇ 1)] t (3-7)
  • F (n, 1: n ⁇ 1) represents a horizontal vector indicating the components in the 1st to n ⁇ 1th columns of the nth row of the matrix F.
  • Step S6 The n-th interference subtraction unit a133-n subtracts fn from the specific signal sn addressed to the n-th terminal apparatus B3n to calculate a signal sn-fn.
  • Step S7 The modulo operation unit a132-n applies the modulo operation to sn-fn to calculate the signal vn. By this modulo calculation, the transmission signal power addressed to each terminal apparatus is reduced.
  • x be a signal obtained by multiplying t by the linear filter P.
  • each component of the signal x is a transmission signal transmitted in order from each antenna a101-1 to a101-N. This signal x is input to the frame configuration unit a142.
  • THP can be considered as a perturbation vector search algorithm in the same way as VP.
  • it is not an algorithm that searches for a truly optimal perturbation vector, but an algorithm that searches for a suboptimal perturbation vector with a small amount of computation compared to VP.
  • the nonlinear precoding unit a33 x P 0 (I + F) ⁇ 1 s (3-8) The operation is performed.
  • the signal x calculated by the non-linear precoding unit a33 is not necessarily a signal obtained by adding the optimum perturbation vector expressed by the equation (1-2), but becomes a sub-optimal signal x because power suppression is performed by modulo calculation. When the modulo operation is not performed, the norm is always smaller than Ws.
  • THP using ordering or LR-THP is applicable similarly.
  • a program for realizing the functions described in the present embodiment is recorded on a computer-readable recording medium, and the program recorded on the recording medium is read into a computer system and executed to execute processing of each unit. May be performed.
  • the “computer system” here includes an OS and hardware such as peripheral devices.
  • the “computer system” includes a homepage providing environment (or display environment) if a WWW system is used.
  • the “computer-readable recording medium” means a storage device such as a flexible disk, a magneto-optical disk, a portable medium such as a ROM and a CD-ROM, and a hard disk incorporated in a computer system. Furthermore, the “computer-readable recording medium” dynamically holds a program for a short time like a communication line when transmitting a program via a network such as the Internet or a communication line such as a telephone line. In this case, a volatile memory in a computer system serving as a server or a client in that case, and a program that holds a program for a certain period of time are also included.
  • the program may be a program for realizing a part of the above-described functions, or may be a program that can realize the above-described functions in combination with a program already recorded in a computer system.
  • the present invention can be used for communication devices.
  • CRS generation unit a142 ... frame configuration unit, a143 ... IFFT unit, a144 ... GI insertion unit, a145 ... transmission unit, a226 ... CRS correction unit, b12 ... propagation for DMRS Route estimation unit, b101 ... antenna, b102 ... receiving unit, b103 ... GI removal B104 ... FFT unit, b105 ... signal separation unit, b106 ... propagation path compensation unit, b107 ... CRS propagation path estimation unit, b108 ... propagation path state information generation unit, b109 ... modulo operation unit, b110 ... demodulation unit, b111 Decoding unit, b121 ... Perturbation vector candidate selection unit, b122 ...
  • Perturbation vector addition unit b123 ... Temporary propagation channel estimation unit, b124 ... Propagation channel compensation unit, b125 ... Modulo operation unit, b126 ... Demodulation unit, b127 ... Perturbation vector evaluation Value calculation unit, b128 ... perturbation vector estimation unit, b131 ... IFFT unit, b132 ... GI insertion unit, b133 ... transmission unit, 300 ... two-dimensional Euclidean mutual removal unit, 301 ... vector storage unit, 303 ... difference vector calculation unit, 305 ... Difference vector norm calculation unit, 307 ... Convergence determination unit, 309 ... Norm calculation unit, 320 ... About vector checking unit, 340 ... perturbation vector adding section, 350 ... complex gain calculation unit, 351 ... temporary complex gain calculator, 352 ... DMRS channel compensation unit, 353 ... DMRS-Modulo unit, 354 ... Vector divider.

Landscapes

  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Power Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Transmission System (AREA)

Abstract

複数の受信装置に対して同一時刻・同一周波数でデータ信号を送信する送信装置であって、第１の受信装置宛の復調用参照信号と、前記第１の受信装置とは異なる第２の受信装置宛の復調用参照信号とを、同一時刻・同一周波数で送信する送信部、を有することを特徴とする送信装置。

Description

送信装置、受信装置、通信システム、通信方法、および集積回路
　本発明は、送信装置、受信装置、通信システム、通信方法、および集積回路に関する。
＜ＭＩＭＯ＞
　近年、無線データ通信の高速化を限られた周波数帯域で実現するため、周波数利用効率向上のための研究が多くなされてきた。その中でも、複数のアンテナを同時に利用することにより、単位周波数当たりの伝送容量を増やすＭＩＭＯ（Ｍｕｌｔｉ－Ｉｎｐｕｔ　Ｍｕｌｔｉ－Ｏｕｔｐｕｔ；多入力多出力）技術は注目されている。
＜ＬＰ　ＭＵ－ＭＩＭＯ＞
　ＭＩＭＯには、基地局装置（Ｂａｓｅ　Ｓｔａｔｉｏｎ：　ＢＳ）が１つの端末装置（Ｍｏｂｉｌｅ　Ｓｔａｔｉｏｎ：ＭＳ）に同一時刻・同一周波数で複数の信号を送信するＳｉｎｇｌｅ－Ｕｓｅｒ　ＭＩＭＯ（ＳＵ－ＭＩＭＯ）と、異なる端末装置に同一時刻・同一周波数で信号を送信するＭｕｌｔｉ－Ｕｓｅｒ　ＭＩＭＯ（ＭＵ－ＭＩＭＯ）がある。
　ＳＵ－ＭＩＭＯは、端末装置が持つアンテナ数より多くのストリームを多重できないため、最大ストリーム数は端末装置の物理的なアンテナ数により制限される。一方で、基地局装置は端末装置よりも多くのアンテナを持つことができるため、基地局装置の余ったアンテナを最大限活かすためには、ＭＵ－ＭＩＭＯも必須となる。既にＬＴＥ（Ｌｏｎｇ　Ｔｅｒｍ　Ｅｖｏｌｕｔｉｏｎ）やＬＴＥ－Ａｄｖａｎｃｅｄでは、線形プレコーディング（Ｌｉｎｅａｒ　Ｐｒｅｃｏｄｉｎｇ：ＬＰ）を用いたダウンリンク（Ｄｏｗｎｌｉｎｋ：ＤＬ）ＭＵ－ＭＩＭＯが仕様化されている（下記非特許文献１参照）。
＜ＮＬＰ　ＭＵ－ＭＩＭＯ＞
　しかし、ＬＰを用いたＭＵ－ＭＩＭＯ（ＬＰ　ＭＵ－ＭＩＭＯ）は、基地局装置が線形フィルタを乗算することで、送信信号を直交させて、端末装置間の干渉（Ｍｕｌｔｉ－Ｕｓｅｒ　Ｉｎｔｅｒｆｅｒｅｎｃｅ：ＭＵＩ）を除去しなければならず、空間多重できる端末装置の組み合わせの柔軟性が低下してしまう。
　一方で、空間多重を実現する別の方法として、非線形プレコーディング（Ｎｏｎｌｉｎｅａｒ　Ｐｒｅｃｏｄｉｎｇ：ＮＬＰ）ＭＵ－ＭＩＭＯが提案されている。ＮＬＰ　ＭＵ－ＭＩＭＯにおいて、端末装置は、受信信号を同相成分（Ｉｎ－ｐｈａｓｅ　ｃｈａｎｎｅｌ：Ｉ－ｃｈ）と直交成分（Ｑｕａｄｒａｔｕｒｅ　ｃｈａｎｎｅｌ：Ｑ－ｃｈ）の方向に、一定の幅（Ｍｏｄｕｌｏ幅）の整数倍だけ平行移動した点を同一の点とみなすＭｏｄｕｌｏ演算を行う。これにより、基地局装置は、Ｍｏｄｕｌｏ幅の任意の整数倍の信号（摂動ベクトル）を変調信号に加算可能となり、摂動ベクトルを適切に選択して各端末装置宛の信号に加算して送信電力を低減する（下記非特許文献２参照）。
＜ＶＰ　ＭＵ－ＭＩＭＯ＞
　端末装置が受信信号に対してＭｏｄｕｌｏ演算を施すことによって、基地局装置は、各変調信号に対してＭｏｄｕｌｏ幅の任意の整数倍の信号を加算する自由度を得る。この加算可能な信号を摂動ベクトル（Ｐｅｒｔｕｒｂａｔｉｏｎ　Ｖｅｃｔｏｒ）と呼ぶ。そして、この摂動ベクトルのうち最も電力効率を改善するものを、空間多重する全ての端末装置の伝搬路状態を考慮して、全探索する方法がＶＰ（Ｖｅｃｔｏｒ　Ｐｅｒｔｕｒｂａｔｉｏｎ）ＭＵ－ＭＩＭＯ方式である。ＶＰ　ＭＵ－ＭＩＭＯは基地局装置の演算量が大きいものの、フル送信ダイバーシチ利得を得ることができ、非常に良好な特性を示すＮＬＰ　ＭＵ－ＭＩＭＯ方式である（下記非特許文献２参照）。
＜ＴＨＰ　ＭＵ－ＭＩＭＯ＞
　ＶＰ　ＭＵ－ＭＩＭＯと異なり、各端末装置が受けるユーザ間干渉を考慮して、逐次的に各端末装置宛の信号に加算する摂動ベクトルを算出する方法をＴＨＰ（Ｔｏｍｌｉｎｓｏｎ－Ｈａｒａｓｈｉｍａ　ｐｒｅｃｏｄｉｎｇ）ＭＵ－ＭＩＭＯと呼ぶ。ＴＨＰ　ＭＵ－ＭＩＭＯは、基地局装置の送信処理の複雑度が低いものの、全端末装置でフル送信ダイバーシチを得ることができない（下記非特許文献３参照）。
＜ＬＲ－ＴＨＰ＞
　また、ＴＨＰ　ＭＵ－ＭＩＭＯに、格子基底縮小（Ｌａｔｔｉｃｅ　Ｒｅｄｕｃｔｉｏｎ：ＬＲ）という処理を加えることで、ＶＰ　ＭＵ－ＭＩＭＯより少ない演算量でフル送信ダイバーシチ利得を得ることができる方法がＬＲ－ＴＨＰである（下記非特許文献３参照）。
＜ＤＭＲＳ＞
　ＮＬＰ　ＭＵ－ＭＩＭＯシステムにおいては、基地局装置が、ＤＭＲＳ（ＤｅＭｏｄｕｌａｔｉｏｎ　Ｒｅｆｅｒｅｎｃｅ　Ｓｉｇｎａｌ：復調用参照信号）を各端末装置に対して送信する必要がある。しかし、基地局装置がデータ信号と同じ非線形プレコーディングをＤＭＲＳに施して送信しても、端末装置は、伝搬路推定することができないという問題がある。
　ＤＭＲＳは、基地局装置が、ＮＬＰ　ＭＵ－ＭＩＭＯによってプレコーディングを施したデータ信号の振幅と位相を予め各端末装置に通知するための信号である。ＤＭＲＳに非線形プレコーディングを施したとすると基地局装置は、ＤＭＲＳに対して摂動ベクトルを加算して（またはＭｏｄｕｌｏ演算を施して）から送信するため、端末装置においてもＤＭＲＳに対してＭｏｄｕｌｏ演算を施す必要がある。そのため、端末装置は、Ｍｏｄｕｌｏ演算に必要なＭｏｄｕｌｏ幅をあらかじめ知っておく必要がある。Ｍｏｄｕｌｏ幅は、受信信号における変調信号の振幅に比例するため、端末装置は、非線形プレコーディング後のデータ信号の受信ゲイン（伝搬路の複素利得）を知る必要がある。しかし、端末装置は、ＤＭＲＳを用いて伝搬路推定をしなければ、受信ゲイン（伝搬路の複素利得）を知ることができない。つまり、端末装置は、「ＤＭＲＳから伝搬路推定しなければ受信ゲイン（伝搬路の複素利得）が得られないが、その受信ゲインを知らなければ、ＤＭＲＳを伝搬路推定できない。」という状態となるため、上記問題が起こる。
　そこで、特許文献１に記載の技術では、直交した無線リソース（時間方向および周波数方向に分割した領域であって、異なるデータ信号や参照信号を割り当てても互いに干渉しない領域）で、各端末装置にＤＭＲＳを送信する。この場合、端末装置は、ＤＭＲＳにＭｏｄｕｌｏ演算を施す必要が無いため、受信ゲインを得ることが出来る。
特開２００９－１８２８９４号公報
３ＧＰＰ　Ｔｅｃｈｎｉｃａｌ　Ｓｐｅｃｉｆｉｃａｔｉｏｎ　３６．２１１　ｖ８．９．０ Ｂ．Ｍ．Ｈｏｃｈｗａｌｄ，Ｃ．Ｂ．Ｐｅｅｌ，Ａ．Ｌ．Ｓｗｉｎｄｌｅｈｕｒｓｔ，"Ａ　Ｖｅｃｔｏｒ－Ｐｅｒｔｕｒｂａｔｉｏｎ　Ｔｅｃｈｎｉｑｕｅ　ｆｏｒ　Ｎｅａｒ－Ｃａｐａｃｉｔｙ　Ｍｕｌｔｉａｎｔｅｎｎａ　Ｍｕｌｔｉｕｓｅｒ　Ｃｏｍｍｕｎｉｃａｔｉｏｎ－Ｐａｒｔ　ＩＩ　Ｐｅｒｔｕｒｂａｔｉｏｎ，"　ＩＥＥＥ　ＴＲＡＮＳＡＣＴＩＯＮＳ　ＯＮ　ＣＯＭＭＵＮＩＣＡＴＩＯＮＳ，　ＶＯＬ．５３，ＮＯ．３，ｐｐ．５３７－５４４，Ｍａｒｃｈ　２００５ Ｆ. Ｌｉｕ, Ｌ. Ｊｉａｎｇ, Ｃ. Ｈｅ, "Ｌｏｗ ｃｏｍｐｌｅｘｉｔｙ ｌａｔｔｉｃｅ ｒｅｄｕｃｔｉｏｎ ａｉｄｅｄ ＭＭＳＥ ｐｒｅｃｏｄｉｎｇ ｄｅｓｉｇｎ ｆｏｒ ＭＩＭＯ ｓｙｓｔｅｍｓ," Ｐｒｏｃｅｅｄｉｎｇｓ ｏｆ ＩＣＣ ２００７, ｐｐ. ２５９８ － ２６０３, Ｊｕｎｅ ２００７.
　しかしながら、上記特許文献１に示したように、直交した無線リソースに各端末装置宛のＤＭＲＳを配置する方法は、ＤＭＲＳ専用の無線リソースが端末装置数分だけ必要であるため、ＤＭＲＳ挿入によるオーバーヘッドの増加が大きくなるという問題があった。
　本発明は、係る事情に鑑みてなされたものであり、ＮＬＰ　ＭＵ－ＭＩＭＯシステムにおいてＤＭＲＳを空間多重してＤＭＲＳ挿入によるオーバーヘッドの増加を最小限に抑えながら、端末装置において正常にデータ信号の振幅および位相（複素利得・受信ゲイン）を推定できる技術を提供することを目的とする。
　本発明の一観点によれば、複数の受信装置に対して同一時刻・同一周波数でデータ信号を送信する送信装置であって、第１の受信装置宛の復調用参照信号と、前記第１の受信装置とは異なる第２の受信装置宛の復調用参照信号とを、同一時刻・同一周波数で送信する送信部、を有することを特徴とする送信装置が提供される。
　送信装置において、空間多重された復調用参照信号（DMRS）は、一つの無線リソースで全移動局装置に対して復調用参照信号を送信できるので、データ信号を配置するための無線リソースを十分確保しながら、複数の復調用参照信号を全移動局装置に送信することができる。また各移動局装置（下記受信装置）は、複数の復調用参照信号を用いて伝搬路推定することができる。そのため本発明を用いれば復調用参照信号の挿入損を大きく低減できる。
　前記復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算する非線形プレコーディング部を有することが好ましい。また、前記復調用参照信号に対して、非線形プレコーディングを施す非線形プレコーディング部を有するようにしても良い。
　前記送信部は、同一時刻・同一周波数で複数のデータ信号を送信し、前記非線形プレコーディング部は、前記復調用参照信号を、前記データ信号と同じ非線形プレコーディング処理することが好ましい。これにより、データ信号と全く同じフィルタを用いて同じ原理で非線形プレコーディングを行うことができる。
　また、復調用参照信号を補正するＤＭＲＳ補正部を有するようにしても良い。前記ＤＭＲＳ補正部は、前記復調用参照信号に２次元ユークリッド互除法を適用する２次元ユークリッド互除法部を有することが好ましい。また、前記２次元ユークリッド互除法部は、第１の前記復調用参照信号から、第２の前記復調用参照信号、第２の前記復調用参照信号の位相を９０度回転させた信号、第２の前記復調用参照信号の位相を１８０度回転させた信号、および第２の前記復調用参照信号の位相を２７０度回転させた信号を減算して、差分ベクトル算出部を有するようにすると良い。
　また、本発明は、復調用参照信号にあらかじめ決めた所定の幅の整数倍の信号を加算する摂動ベクトル加算部を有することを特徴とする受信装置である。前記信号を加算した復調用参照信号を用いて伝搬路推定する仮伝搬路推定部を有することが好ましい。複数の異なる前記信号を選択する摂動ベクトル候補選択部と、前記信号それぞれを用いて前記仮伝搬路推定部で伝搬路推定した伝搬路推定結果に基づいて、１つの前記信号を選択する摂動ベクトル推定部を有するようにすると良い。また、複数の異なる前記信号に対応する前記伝搬路推定結果それぞれに基づいてデータ信号を軟推定して対数尤度比をそれぞれ算出する復調部と、前記対数尤度比ぞれぞれの分散を算出する摂動ベクトル評価値算出部と、を有し、前記摂動ベクトル推定部は、前記分散のうち最も大きいものに対応する前記信号を選択することが好ましい。複数の前記復調用参照信号に対して２次元ユークリッドの互除法を適用し、既約ベクトルを算出する２次元ユークリッドの互除法部を有することが好ましい。前記既約ベクトルを用いて伝搬路の複素利得を算出する複素利得算出部を有することが好ましい。
　また、本発明は、復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算する非線形プレコーディング部と、復調用参照信号と他の信号を、同一時刻・同一周波数で送信する送信部とを有する送信装置と、復調用参照信号にあらかじめ決めた所定の幅の整数倍の信号を加算する摂動ベクトル加算部を有する受信装置と、を有することを特徴とする通信システムである。
　また、本発明は、復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算するステップと、復調用参照信号と他の信号を、同一時刻・同一周波数で送信するステップを有する送信方法と、復調用参照信号にあらかじめ決めた所定の幅の整数倍の信号を加算するステップを有する受信方法と、を有することを特徴とする通信方法である。
　また、本発明は、復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算する非線形プレコーディング部と、復調用参照信号と他の信号を、同一時刻・同一周波数で送信するステップを有する送信部と、を有することを特徴とする集積回路である。
　また、本発明は、復調用参照信号にあらかじめ決めた所定の幅の整数倍の信号を加算する摂動ベクトル加算部を有することを特徴とする集積回路である。
　まあ、本発明は、上記に記載の通信方法を、コンテンツに実行させるためのプログラムであっても良く、当該プログラムを記録するコンピュータ読み取り可能な記録媒体であっても良い。
　本明細書は本願の優先権の基礎である日本国特許出願２０１１－０４３７５５号の明細書および／または図面に記載される内容を包含する。
本発明の第１の実施形態に係る通信システムの一構成例を示す概念図である。 本実施形態に係る通信システムの動作の一例を示すシーケンス図である。 本実施形態に係る基地局装置の概略構成例を示す機能ブロック図である。 本実施形態に係る固有信号の構成の一例を示す概略図である。 本実施形態に係るフレームの構成の一例を示す概略図である。 第１から第Ｎまでの各アンテナａ１０１－ｎで送信するフレームの第１の構成例を示す図である。 第１から第Ｎまでの各アンテナａ１０１－ｎで送信するフレームの第２の構成例を示す図である。 本実施形態に係る端末装置の一構成例を示す機能ブロック図である。 ＤＭＲＳ伝搬路推定部の詳細な構成を示す機能ブロック図である。 ２つのＤＭＲＳ（ＤＭＲＳ１およびＤＭＲＳ２と呼ぶ。）の位置の一例を信号点平面上に示す図である。 既約ベクトルが、ＱＰＳＫの信号点と一致しない位置の一例を信号点平面上に示す図である。 ２次元ユークリッドの互除法部の一構成例を示す図である。 ＤＭＲＳ補正部の一構成例を示す図である。 ＤＭＲＳ伝搬路推定部（ａ）と、複素利得算出部（ｂ）、とを示す図である。 非線形プレコーディング部の一構成例を示す図である。 非線形プレコーディング部の動作を示すフローチャート図である。 本発明の第２の実施の形態による基地局装置の概略構成例を示す機能ブロック図である。
　以下、本発明を実施するための形態について、図面を参照しながら詳細に説明する。
（第１の実施形態：　対数尤度比の分散ベース）
＜通信システム１について＞
　図１は、本発明の第１の実施形態に係る通信システム１の一構成例を示す概念図である。通信システム１は、基地局装置Ａ１、および、第１から第Ｎまでの端末装置Ｂ１１～Ｂ１Ｎを具備する（図１は、基地局装置Ａ１が、第１から第４までの端末装置Ｂ１１、Ｂ１２、Ｂ１３、およびＢ１４を選択した場合の一例（Ｎ＝４）を示す図である。）。
　基地局装置Ａ１は、共通参照信号（Ｃｏｍｍｏｎ　Ｒｅｆｅｒｅｎｃｅ　Ｓｉｇｎａｌ：ＣＲＳ）を送信する。なお、ＣＲＳは、基地局装置Ａ１と第１から第４までの端末装置Ｂ１１～Ｂ１４とが、その基準信号を予め記憶する信号である。第１から第４までの端末装置Ｂ１１～Ｂ１４各々は、基地局装置Ａ１が送信したＣＲＳに基づいて伝搬路状態を推定し、推定した伝搬路状態に基づいて伝搬路状態情報を、基地局装置Ａ１に通知する。
　基地局装置Ａ１は、第１から第４までの端末装置Ｂ１１～Ｂ１４に対して、ＤＭＲＳおよびデータ信号を送信する。ここで、基地局装置Ａ１は、ＤＭＲＳおよびデータ信号に対して、プレコーディングを施し、乗算後のＤＭＲＳおよびデータ信号を送信する。
　多重された第１から第４までの端末装置Ｂ１１～Ｂ１４は、基地局装置Ａ１から受信したＤＭＲＳに基づいて、プレコーディング処理を伝搬路の一部とみなした等価伝搬路（以下、単に等価伝搬路と称する。）の伝搬路状態を推定し、推定した等価伝搬路の伝搬路状態を示す等価伝搬路状態情報に基づいてデータ信号を取得する。
　図２は、本実施形態に係る通信システム１の動作の一例を示すシーケンス図である。この図は、図１の場合の通信システム１の動作の例を示す図である。
（ステップＳ１０１）　基地局装置Ａ１は、第１から第４までの端末装置Ｂ１１～Ｂ１４に、ＣＲＳを送信する。その後、ステップＳ１０２に進む。
（ステップＳ１０２）　第１から第４までの端末装置Ｂ１１～Ｂ１４は、ステップＳ１０１で送信されたＣＲＳに基づいて伝搬路状態を推定する。その後、ステップＳ１０３に進む。
（ステップＳ１０３）　第１から第４までの端末装置Ｂ１１～Ｂ１４は、ステップＳ１０２で推定した伝搬路状態に基づいて伝搬路状態情報を算出する。その後、ステップＳ１０４に進む。
（ステップＳ１０４）　第１から第４までの端末装置Ｂ１１～Ｂ１４は、ステップＳ１０３で算出した伝搬路状態情報を、基地局装置Ａ１へ通知する。その後、ステップＳ１０５に進む。
（ステップＳ１０５）　基地局装置Ａ１は、ステップＳ１０４で通知された伝搬路状態情報に基づいて、非線形プレコーディングに用いるフィルタを算出する。基地局装置Ａ１は、生成したＤＭＲＳおよびデータ信号に対して、当該フィルタを用いて非線形プレコーディングを行い、ＤＭＲＳおよびデータ信号を生成する。その後、ステップＳ１０６に進む。
（ステップＳ１０６）　基地局装置Ａ１は、ステップＳ１０５で生成したＤＭＲＳの信号を、第１から第４までの端末装置Ｂ１１、Ｂ１２、Ｂ１３、およびＢ１４へ送信する。その後、ステップＳ１０７に進む。
（ステップＳ１０７）　第１から第４までの端末装置Ｂ１１～Ｂ１４は、ステップＳ１０６で送信されたＤＭＲＳの信号に基づいて、等価伝搬路の伝搬路状態を推定する。その後、ステップＳ１０８に進む。
（ステップＳ１０８）　基地局装置Ａ１は、ステップＳ１０５で生成したデータ信号を、各端末装置Ｂ１１～Ｂ１４へ送信する。その後、ステップＳ１０９に進む。
（ステップＳ１０９）　第１から第４までの端末装置Ｂ１１～Ｂ１４は、ステップＳ１０８で推定した等価伝搬路の伝搬路状態を示す等価伝搬路状態情報に基づいて、データ信号を検出して取得する。
　以下、図１・２では端末装置Ｂ１１～Ｂ１４まで４個多重したが、以降Ｎ個の第１から第Ｎまでの端末装置Ｂ１１～Ｂ１Ｎを多重するとする。また、端末装置Ｂ１１～Ｂ１Ｎを総称して端末装置Ｂ１ｎと呼ぶ。
＜基地局装置Ａ１について＞
　図３は、本実施形態に係る基地局装置Ａ１の概略構成例を示す機能ブロック図である。この図において、基地局装置Ａ１は、第１から第Ｎまでのアンテナａ１０１－１～ａ１０１－Ｎ、第１から第Ｎまでの受信部ａ１０２－１～ａ１０２－Ｎ、第１から第ＮまでのＧＩ（Ｇｕａｒｄ Ｉｎｔｅｒｖａｌ；ガードインターバル）除去部ａ１０３－１～ａ１０３－Ｎ、第１から第ＮまでのＦＦＴ（Ｆａｓｔ　Ｆｏｕｒｉｅｒ　Ｔｒａｎｓｆｏｒｍ；高速フーリエ変換）部ａ１０４－１～ａ１０４－Ｎ、伝搬路状態情報取得部ａ１０５、フィルタ算出部ａ１１、第１から第Ｎまでの符号部ａ１２１－１～ａ１２１－Ｎ、第１から第Ｎまでの変調部ａ１２２－１～ａ１２２－Ｎ、ＤＭＲＳ生成部ａ１２４、固有信号構成部ａ１２５、非線形プレコーディング部ａ１３、ＣＲＳ生成部ａ１４１、フレーム構成部ａ１４２、第１から第ＮまでのＩＦＦＴ（Ｉｎｖｅｒｓｅ　Ｆａｓｔ　Ｆｏｕｒｉｅｒ　Ｔｒａｎｓｆｏｒｍ；逆高速フーリエ変換）部ａ１４３－１～ａ１４３－Ｎ、第１から第ＮまでのＧＩ挿入部ａ１４４－１～ａ１４４－Ｎ、および、第１から第Ｎまでの送信部ａ１４５－１～ａ１４５－Ｎを含んで構成される。
　なお、図３の基地局装置Ａ１は、Ｎ本のアンテナａ１０１－１～ａ１０１－Ｎを備え、Ｎ個の端末装置を多重する場合の基地局装置である（例えば、図１、２の例では、Ｎ＝４）。また、図３の基地局装置Ａ１では、一例として上りリンクおよび下りリンクともに直交周波数分割多重（Ｏｒｔｈｏｇｏｎａｌ　Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ：ＯＦＤＭ）方式を用いる場合について説明するが、本発明はこれに限らず、基地局装置Ａ１は、上りリンクおよび下りリンクの一方又は両方で、時間分割多重（Ｔｉｍｅ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ：ＴＤＭ）方式や周波数分割多重（Ｆｒｅｑｕｅｎｃｙ　Ｄｉｖｉｓｉｏｎ　Ｍｕｌｔｉｐｌｅｘｉｎｇ：ＦＤＭ）方式を用いてもよい。
　第１から第Ｎまでの受信部ａ１０２－ｎ（ｎ＝１、２、・・・、Ｎ）は、第１から第Ｎまでのアンテナａ１０１－ｎを介して、各端末装置Ｂ１ｎから送信された信号（搬送波周波数の信号）を受信する。この信号には、伝搬路状態情報が含まれる。第１から第Ｎまでの受信部ａ１０２－ｎは、受信した信号をダウンコンバージョンし、Ａ／Ｄ（アナログ／デジタル）変換することで、ベースバンドのデジタル信号を生成する。第１から第Ｎまでの受信部ａ１０２－ｎは、生成したデジタル信号を第１から第ＮまでのＧＩ除去部ａ１０３－ｎに出力する。
　第１から第ＮまでのＧＩ除去部ａ１０３－ｎは、第１から第Ｎまでの受信部ａ１０２－ｎから入力されたデジタル信号からＧＩを除去し、除去後の信号を第１から第ＮまでのＦＦＴ部ａ１０４－ｎに出力する。
　第１から第ＮまでのＦＦＴ部ａ１０４－ｎは、第１から第ＮまでのＧＩ除去部ａ１０３－ｎから入力された信号に対して、ＦＦＴを行うことで、周波数領域の信号を生成する。第１から第ＮまでのＦＦＴ部ａ１０４－ｎは、生成した周波数領域の信号を伝搬路状態情報取得部ａ１０５に出力する。
　伝搬路状態情報取得部ａ１０５は、第１から第ＮまでのＦＦＴ部ａ１０４－ｎから入力された信号を復調し、復調した情報から伝搬路状態情報を抽出する。伝搬路状態情報取得部ａ１０５は、抽出した伝搬路状態情報をフィルタ算出部ａ１１に出力する。なお、第１から第ＮまでのＦＦＴ部ａ１０４－ｎが出力する信号のうち伝搬路状態情報の信号以外の信号は、制御部（図示せず）で復調される。復調された情報のうち制御情報は基地局装置Ａ１の制御に用いられ、また、制御情報以外のデータは他の基地局装置やサーバ装置等へ送信される。
　フィルタ算出部ａ１１は、伝搬路状態情報取得部ａ１０５から入力された伝搬路状態情報に基づいて、非線形プレコーディングに用いるフィルタを算出する。フィルタ算出部ａ１１が行うフィルタ算出処理の詳細については、後述する。フィルタ算出部ａ１１は、算出したフィルタを非線形プレコーディング部ａ１３に入力する。
　第１から第Ｎまでの符号部ａ１２１－ｎには、各端末装置Ｂ１ｎ宛（例えば、図１、２の例では、Ｎ＝４）の情報ビット（データ）が入力される。第１から第Ｎまでの符号部ａ１２１－ｎは、入力された情報ビットを誤り訂正符号化し、符号化後の符号化ビットを第１から第Ｎまでの変調部ａ１２２－ｎに出力する。
　第１から第Ｎまでの変調部ａ１２２－ｎは、第１から第Ｎまでの符号部ａ１２１－ｎから入力された符号化ビットを変調することで、端末装置Ｂ１ｎ宛のデータ信号を生成する。第１から第Ｎまでの変調部ａ１２２－ｎは、生成したデータ信号を固有信号構成部ａ１２５に出力する。なお、基地局装置Ａ１は、伝搬路状態情報に基づいて第１から第Ｎまでの変調部ａ１２２－ｎで用いる変調方式を決定し、決定した変調方式を示す変調情報を、変調部ａ１２２－ｎに出力し、基地局装置Ａ１から第１から第Ｎまでの端末装置Ｂ１ｎへ通知する。
　ＤＭＲＳ生成部ａ１２４は、第１から第Ｎまでの各端末装置Ｂ１ｎ宛のＤＭＲＳを生成する。ＤＭＲＳ生成部ａ１２４は、生成したＤＭＲＳを、固有信号構成部ａ１２５に出力する。
　固有信号構成部ａ１２５は、第１から第Ｎまでの各変調部１２２－ｎから入力された第１から第Ｎまでの端末装置Ｂ１ｎ宛のデータ信号とＤＭＲＳ生成部ａ１２４から入力された第１から第Ｎまでの端末装置Ｂ１ｎ宛のＤＭＲＳとを関係付ける。固有信号構成部ａ１２５が関係付けた第１から第Ｎまでの端末装置Ｂ１ｎ宛の情報各々を、第１から第Ｎまでの端末装置Ｂ１ｎの固有信号という。固有信号構成部ａ１２５は、関係付けることで生成した第１から第Ｎまでの端末装置Ｂ１ｎの固有信号各々を、非線形プレコーディング部ａ１３に出力する。
　非線形プレコーディング部ａ１３は、固有信号構成部ａ１２５から入力された第１から第Ｎまでの端末装置Ｂ１ｎの固有信号（データ信号およびＤＭＲＳ）に対して非線形プレコーディングを行う。非線形プレコーディング部ａ１３が行う非線形プレコーディングの詳細については、後述する。非線形プレコーディング部ａ１３は、非線形プレコーディングを行った固有信号を、フレーム構成部ａ１４２に出力する。
　ＣＲＳ生成部ａ１４１は、基地局装置Ａ１と第１から第Ｎまでの端末装置Ｂ１ｎで既知の基準信号を有するＣＲＳを生成し、生成したＣＲＳをフレーム構成部ａ１４２に出力する。
　フレーム構成部ａ１４２は、非線形プレコーディング部ａ１３から入力された固有信号、および、ＣＲＳ生成部ａ１４１から入力されたＣＲＳをマッピングする。なお、フレーム構成部ａ１４２は、固有信号とＣＲＳとを別のフレームにマッピングしてもよいし、同じフレームにマッピングしてもよい。例えば、ＣＲＳのみをあるフレームにマッピングし、ＣＲＳおよび固有信号を他のフレームにマッピングしてもよい。なお、基地局装置Ａ１は、あらかじめ決められたマッピングに従ってＣＲＳと固有信号をフレームにマッピングし、第１から第Ｎまでの端末装置Ｂ１ｎは、マッピングをあらかじめ把握しているものとする。
　フレーム構成部ａ１４２は、マッピング後の信号のうちアンテナａ１０１－ｎで送信する信号を、第１から第ＮまでのＩＦＦＴ部ａ１４３－ｎに、フレーム単位で出力する。
　第１から第ＮまでのＩＦＦＴ部ａ１４３－ｎは、フレーム構成部ａ１４２から入力された信号に対して、ＩＦＦＴを行うことで、時間領域の信号を生成する。第１から第ＮまでのＩＦＦＴ部ａ１４３－ｎは、生成した時間領域の信号を第１から第ＮまでのＧＩ挿入部ａ１４４－ｎに出力する。
　第１から第ＮまでのＧＩ挿入部ａ１４４－ｎは、第１から第ＮまでのＩＦＦＴ部ａ１４３－ｎから入力された信号に対して、ガードインターバルを付与し、付与後の信号を第１から第Ｎまでの送信部ａ１４５－ｎに出力する。
　第１から第Ｎまでの送信部ａ１４５－ｎは、第１から第ＮまでのＧＩ挿入部ａ１４４－ｎから入力された信号（ベースバンドのデジタル信号）をＤ／Ａ（デジタル／アナログ）変換する。第１から第Ｎまでの送信部ａ１４５－ｎは、変換後の信号をアップコンバージョンすることで搬送波周波数の信号を生成する。第１から第Ｎまでの送信部ａ１４５－ｎは、生成した信号を第１から第Ｎまでのアンテナａ１０１－ｎを介して送信する。
　図４は、本実施形態に係る固有信号の構成の一例を示す概略図である。この図において、横軸は時間を表す。この図は、固有信号構成部ａ１２５が出力した固有信号を表す。また、この図は、同一の周波数で送信する、第１から第Ｎまでの端末装置Ｂ１ｎの固有信号の配置を時間軸を揃えて表したものである。
　なお、図４では、第１から第Ｎまでの端末装置Ｂ１ｎ宛のＤＭＲＳを「ＤＭＲＳ－ＭＳｎ」との表記で表している。この図は、固有信号が、第１から第Ｎまでの端末装置Ｂ１ｎ毎の信号であってＤＭＲＳとデータ信号とから構成されることを示す。例えば、符号Ｓ１１を付した信号は、端末装置Ｂ１１宛のＤＭＲＳ（ＤＭＲＳ‐ＭＳ１）を示す。
　図４は、第１から第Ｎまでの端末装置Ｂ１ｎ宛のＤＭＲＳ各々が、固有信号構成部ａ１２５から全て同じ時間に出力されることを示す。また、この図は、第１から第Ｎまでの各端末装置Ｂ１ｎ宛のデータ信号も、固有信号構成部ａ１２５から同じ時間に出力されることを示す。なお、上述のように、同時刻に割り当てられたデータ信号同士、ＤＭＲＳ同士はともに非線形プレコーディングａ１３で空間多重され、基地局装置Ａ１から同一時刻・同一周波数（つまり同一ＯＦＤＭシンボルの同一サブキャリア）で送信される。また、図４は、ＤＭＲＳとデータ信号とが、固有信号構成部ａ１２５から異なる時間に出力されることを示す。例えば、第１から第Ｎまでの端末装置Ｂ１ｎ宛のＤＭＲＳは時間ｔ１と時間ｔ２に出力され、データ信号はｔ３～ｔＬに出力される。
　なお、図４における固有信号の構成は一例であり、本発明はこれに限られない。例えば、図４では、端末装置Ｂ１ｎ宛のＤＭＲＳが出力された後（ｔ３以降）にデータ信号が出力されているが、固有信号構成部ａ１２５は、データ信号を出力後にＤＭＲＳを出力してもよい。また、固有信号構成部ａ１２５は、データ信号とＤＭＲＳを時間軸上で交互に出力してもよいし、その他の順序で出力してもよい。
　また、図５Ａは、本実施形態に係るフレームの構成の一例を示す概略図である。この図は、フレーム構成部ａ１４２が信号をマッピングしたフレームの構成を示す。また、この図は、第１から第Ｎまでの各アンテナａ１０１－ｎで送信するフレームの構成を、時間軸を揃えて表したものである。
　なお、図５Ａでは、第１から第Ｎまでのアンテナａ１０１－ｎから送信するＣＲＳを「ＣＲＳ－Ｔｘｎ」との表記で表している。この図は、第１から第Ｎまでのアンテナａ１０１－ｎ毎のフレーム各々が、ＣＲＳと固有信号（データ信号とＤＭＲＳから構成される非線形プレコーディング後の固有信号）とから構成されることを示す。例えば、符号Ｓ２１を付した信号は、アンテナａ１０１－１で送信するＣＲＳ（ＣＲＳ－Ｔｘ１）を示す。また、符号Ｓ２２を付した信号は、アンテナａ１０１－１で送信する固有信号であって、非線形プレコーディング後の固有信号を示す。
　図５Ａは、アンテナａ１０１－ｎで送信するＣＲＳ各々が、フレーム構成部ａ１４２で互いに異なる時間帯域に配置され、基地局装置Ａ１から送信されることを示す。例えば、アンテナａ１０１－１で送信されるＣＲＳは時間ｔ１に送信され、アンテナａ１０１－２で送信されるＣＲＳは時間ｔ２に送信される。
　また、図５Ａは、ＣＲＳと固有信号とが、フレーム構成部ａ１４２で異なる時間帯域に配置されることを示す。例えば、アンテナａ１０１－ｎで送信するＣＲＳは時間ｔ１～ｔＮに送信される時間帯域に配置され、固有信号はｔＮ＋１以降の帯域に配置される。また、この図は、第１から第Ｎまでのアンテナａ１０１－ｎで送信する固有信号が、全てフレーム構成部ａ１４２で同じ時間帯域に配置されることを示す。
　なお、図５Ａにおけるフレームの構成は一例であり、本発明はこれに限られない。例えば、図５Ａでは、アンテナａ１０１－ｎで送信するＣＲＳの全てが先の時間帯域（ｔＮ＋１より前）に配置され、後の時間帯域（ｔＮ＋１以降）に固有信号が配置されているが、フレーム構成部ａ１４２は、固有信号を先の時間帯域に配置してＣＲＳを後の時間帯域に配置してもよい。また、フレーム構成部ａ１４２は、固有信号とＣＲＳとを交互する時間帯域に配置してもよいし、他の順序で配置してもよい。例えば、フレーム構成部ａ１４２は、アンテナａ１０１－１で送信するＣＲＳを時間ｔＮ＋１以降の時間帯域に配置し、固有信号を時間ｔＮ＋１より前の時間帯域に配置してもよい。なお、基地局装置Ａ１は、固有信号の送信を開始する前には、ＣＲＳのみを配置したフレームを送信する。
　さらに図５Ｂおよび図５Ｃに、第１から第Ｎまでの各アンテナａ１０１－ｎで送信するフレームの一構成例を示す。図４および図５Ａでは、ＤＭＲＳとデータ信号、または、ＣＲＳと固有信号（ＤＭＲＳとデータ信号が含まれる。）を時間方向に並べた例を図示したが、図５Ｂおよび図５Ｃに示すように、時間方向（ｔ）および周波数方向（ｆ）の２次元マトリックス上にＣＲＳ、ＤＭＲＳ、およびデータ信号を配置してもよい。
＜端末装置Ｂ１ｎについて＞
　図６は、本実施形態に係る端末装置Ｂ１ｎの一構成例を示す機能ブロック図である。この図において、端末装置Ｂ１ｎは、アンテナｂ１０１、受信部ｂ１０２、ＧＩ除去部ｂ１０３、ＦＦＴ部ｂ１０４、信号分離部ｂ１０５、ＣＲＳ用伝搬路推定部ｂ１０７、ＤＭＲＳ用伝搬路推定部ｂ１２、伝搬路補償部ｂ１０６、Ｍｏｄｕｌｏ演算部ｂ１０９、復調部ｂ１１０、復号部ｂ１１１、伝搬路状態情報生成部ｂ１０８、ＩＦＦＴ部ｂ１３１、ＧＩ挿入部ｂ１３２、および、送信部ｂ１３３を含んで構成される。
　受信部ｂ１０２は、アンテナｂ１０１を介して、各端末装置Ｂ１ｎから送信された信号（搬送波周波数の信号）を受信する。受信部ｂ１０２は、受信した信号をダウンコンバージョンし、Ａ／Ｄ（アナログ／デジタル）変換することで、ベースバンドのデジタル信号を生成する。受信部ｂ１０２は、生成したデジタル信号をＧＩ除去部ｂ１０３に出力する。
　ＧＩ除去部ｂ１０３は、受信部ｂ１０２から入力されたデジタル信号からＧＩを除去し、除去後の信号をＦＦＴ部ｂ１０４に出力する。
　ＦＦＴ部ｂ１０４は、ＧＩ除去部ｂ１０３から入力された信号に対して、高速フーリエ変換を行うことで、周波数領域の信号を生成する。ＦＦＴ部ｂ１０４は、生成した周波数領域の信号を信号分離部ｂ１０５に出力する。
　信号分離部ｂ１０５は、基地局装置Ａ１から通知されたマッピング情報に基づいて、ＦＦＴ部ｂ１０４から入力された信号をデマッピングする。信号分離部ｂ１０５は、デマッピングした信号のうち、ＣＲＳをＣＲＳ用伝搬路推定部ｂ１０７に出力し、ＤＭＲＳをＤＭＲＳ用伝搬路推定部ｂ１２に出力する。信号分離部ｂ１０５は、データ信号を伝搬路補償部ｂ１０６に出力する。さらに信号分離部ｂ１０５は、データ信号をＤＭＲＳ用伝搬路推定部ｂ１２に入力する。
　ＣＲＳ用伝搬路推定部ｂ１０７は、信号分離部ｂ１０５から入力されたＣＲＳに基づいて伝搬路状態を推定し、推定した伝搬路状態を示す情報を伝搬路状態情報生成部ｂ１０８に出力する。
　ＤＭＲＳ用伝搬路推定部ｂ１２は、信号分離部ｂ１０５から入力されたＤＭＲＳに基づいて非線形プレコーディングに用いるフィルタを伝搬路の一部とみなした等価伝搬路の伝搬路状態を推定する。なお、ＤＭＲＳ用伝搬路推定部ｂ１２の詳細については後述する。ＤＭＲＳ用伝搬路推定部ｂ１２は、推定した等価伝搬路の伝搬路状態を示す等価伝搬路状態情報を伝搬路補償部ｂ１０６に出力する。伝搬路補償部ｂ１０６は、ＤＭＲＳ用伝搬路推定部ｂ１２から入力された等価伝搬路状態情報を用いて、信号分離部ｂ１０５から入力された信号に対して伝搬路補償を行う。伝搬路補償部ｂ１０６は、伝搬路補償後の信号をＭｏｄｕｌｏ演算部ｂ１０９に出力する。
　Ｍｏｄｕｌｏ演算部ｂ１０９は、基地局装置Ａ１から通知された変調情報に基づいて、伝搬路補償部ｂ１０６から入力された信号に対して、Ｍｏｄｕｌｏ演算を行う。Ｍｏｄｕｌｏ演算は、下式で表すことが出来る。
mod(α)=α-floor((Re[α]+τ/2)/τ)τ-i×floor((Ｉｍ［α］+τ/2)/τ)τ
　　　　　　　　　　　　　　　　　　　　　　　 （１－１）
　ここで、式（１－１）は信号αにＭｏｄｕｌｏ幅τのＭｏｄｕｌｏ演算を施した場合を示している。また、ｆｌｏｏｒ（ｘ）はｘを超えない最大の整数を示し、Ｒｅ［α］およびＩｍ［α］は複素数αの実部および虚部をそれぞれ表す。なお、ｉは虚数単位である。また、τはＭｏｄｕｌｏ幅を表す。データ信号の変調方式がＱＰＳＫならばτは、ＱＰＳＫ信号の平均振幅の２√２倍、１６ＱＡＭならば１６ＱＡＭ信号の８／√１０倍、６４ＱＡＭならば６４ＱＡＭ信号の１６／√４２倍にすることが望ましい。ただし、基地局装置と端末装置とで共通の値を用いればＭｏｄｕｌｏ幅は他の値でもよい。
　Ｍｏｄｕｌｏ演算部ｂ１０９は、Ｍｏｄｕｌｏ演算後の信号を復調部ｂ１１０に出力する。
　復調部ｂ１１０は、基地局装置Ａ１から通知された変調情報が示す変調方式で、Ｍｏｄｕｌｏ演算部ｂ１０９から入力された信号を復調する。復調部ｂ１１０は、復調後の情報（硬判定した符号化ビット又は符号化ビットの軟推定値）を復号部ｂ１１１に出力する。
　復号部ｂ１１１は、復調部ｂ１１０から入力された情報を復号することで、情報ビットを取得し、取得した情報ビットを出力する。
　伝搬路状態情報生成部ｂ１０８は、ＣＲＳ用伝搬路推定部ｂ１０７から入力された伝搬路状態から伝搬路状態情報を生成する（伝搬路状態情報生成処理という）。ここで、推定した伝搬路状態を行ベクトルｈｎ＝［ｈｎ１，ｈｎ２，…，ｈｎＮ］とおく。ｈｎ１，ｈｎ２，…ｈｎＮはそれぞれアンテナａ１０１－１、アンテナａ１０１－２、…、アンテナａ１０１－Ｎと端末装置Ｂ１ｎの間の伝搬路状態である。伝搬路状態情報生成部ｂ１０８は、必ずしも行ベクトル｜ｈｎ｜をそのまま伝搬路状態情報にして基地局装置Ａ１に通知する必要は無い。例えば、行ベクトルのノルム｜ｈｎ｜を、所定の数Ｃに正規化した行ベクトルＣ×ｈｎ／｜ｈｎ｜を伝搬路状態情報として、基地局装置Ａ１に通知しても良いし、あらかじめ決められた値に近似したものを伝搬路状態情報として基地局装置Ａ１に通知してもよい。
　伝搬路状態情報生成部ｂ１０８は、生成した伝搬路状態情報を変調し、変調後の伝搬路状態情報の信号をＩＦＦＴ部ｂ１３１に出力する。
　ＩＦＦＴ部ｂ１３１は、伝搬路状態情報生成部ｂ１０８から入力された信号に対して、逆高速フーリエ変換を行うことで、時間領域の信号を生成する。ＩＦＦＴ部ｂ１３１は、生成した時間領域の信号をＧＩ挿入部ｂ１３２に出力する。
　ＧＩ挿入部ｂ１３２は、ＩＦＦＴ部ｂ１３１から入力された信号に対して、ガードインターバルを付与し、付与後の信号を送信部ｂ１３３に出力する。
　送信部ｂ１３３は、ＧＩ挿入部ｂ１３２から入力された信号（ベースバンドのデジタル信号）をＤ／Ａ（デジタル／アナログ）変換する。送信部ｂ１３３は、変換後の信号をアップコンバージョンすることで搬送波周波数の信号を生成する。送信部ｂ１３３は、生成した信号を、アンテナｂ１０１を介して送信する。
＜フィルタ算出部ａ１１＞
　図３に示すフィルタ算出部ａ１１は、伝搬路状態情報取得部ａ１０５から入力された伝搬路状態情報から伝搬路行列Ｈを構成する。Ｈは、Ｎ行Ｎ列の行列であり、ｐ行ｑ列成分が、ｐ番目の端末装置Ｂ１ｐと基地局装置のｑ番目のアンテナａ１０１－ｑの間の伝搬路の複素利得を示す（ここでｐとｑは１からＮまでの任意の整数である。）。つまり、各端末装置Ｂ１ｐから通知された伝搬路状態情報から取得した伝搬路状態が行ベクトルとなり、全端末装置に対応する行ベクトルを各行に持つ行列を生成することで伝搬路行列Ｈを生成できる。また端末装置Ｂ１ｎで伝搬路状態情報のノルムが正規化されていても基地局装置Ａ１は、端末装置Ｂ１ｎから通知された伝搬路状態情報をそのまま、つまり、正規化後の伝搬路状態を各行に用いて伝搬路行列Ｈを生成する。
　フィルタ算出部ａ１１は、伝搬路行列Ｈの逆行列Ｗ（＝Ｈ－１）を算出し、非線形プレコーディング部ａ１３に入力する。
＜非線形プレコーディング部ａ１３詳細＞
　固有信号構成部ａ１２５から入力された各固有信号をｓｎとおき、全ｓ１～ｓＮを各成分に持つ縦ベクトルをｓとする。ここで固有信号とは、テータ信号およびＤＭＲＳ生成部ａ１２４で生成されたＤＭＲＳを示す。ここでＤＭＲＳ信号はデータ信号とは独立にτの値を決めることが出来るが、基地局装置と端末装置とで共通の値を用いる必要がある。例えばＤＭＲＳ信号がＱＰＳＫ信号の１点に一致するとすれば、データ信号が１６ＱＡＭで変調されていてもτ＝２√２とする。非線形プレコーディング部ａ１３は、フィルタＷを乗算後の送信信号のノルムが最小となるようなＮ次元整数縦ベクトルｚ１とｚ２の組み合わせを探索する。これを式で表すと以下のようになる。
（Ｚ１，Ｚ２）＝ａｒｇｍｉｎ（ｚ１，ｚ２）｜Ｗ（ｓ＋ｚ１τ＋ｉｚ２τ）｜
　　　　　　　　　　　　　　　　　　　　　　　 （１－２）
　右辺のノルム｜Ｗ（ｓ＋ｚ１τ＋ｉｚ２τ）｜を最小にする（ｚ１，ｚ２）の組み合わせを（Ｚ１，Ｚ２）とおく。ｚ１とｚ２の各成分はいかなる整数をも取り得るので、総当たりで探索するのは不可能である。そこで、ｚ１とｚ２の各成分の探索範囲を、あらかじめ決められた範囲（たとえば絶対値がＬＢＳ以下の整数：［－ＬＢＳ，－ＬＢＳ＋１，．．．，－２，－１，０，１，２，．．．，ＬＢＳ－１，ＬＢＳ］）に限定する。このように信号点Ｗｓを中心とした点を候補にして探索を行うのが望ましいが、これ以外の方法で、ノルムを最小とする点を探索してもよい。いまｚ１τ＋ｉｚ２τを摂動ベクトルと呼ぶ。算出した（Ｚ１，Ｚ２）を用いた信号ｘ＝Ｗ（ｓ＋Ｚ１τ＋ｉＺ２τ）の、電力正規化を行う。
　基地局装置Ａ１は、送信電力を一定にするために、一定数のサブキャリアおよび一定数のＯＦＤＭシンボル（「電力正規化単位」と呼ぶ。）内のデータ信号の総送信電力を正規化しなければならない。電力正規化単位は、例えば図５Ｂや図５Ｃに示したフレーム単位全体を示す。まず、非線形プレコーディングにより算出したデータ信号ｘの電力の電力正規化単位に亘る総和Ｐｘを算出する。１つの電力正規化単位のデータ信号の送信に基地局装置Ａ１が割り当てられる総電力がＰｔｒであるとすると、電力正規化係数ｇ（　＝（Ｐｘ／Ｐｔｒ）１／２　）を算出する。非線形プレコーディング部ａ１３は電力正規化係数ｇの逆数を固有信号ｘ（データ信号およびＤＭＲＳ）に乗算し、乗算後のデータ信号をフレーム構成部ａ１４２に入力する。ここで、信号ｘの各成分は、順番に第１から第Ｎまでの各アンテナａ１０１－１～ａ１０１－Ｎで送信する送信信号を示す。
　なお、端末装置Ｂ１ｎでは、ＤＭＲＳ信号自体がｇ－１倍されているので、ＤＭＲＳ用伝搬路推定部ｂ１２で伝搬路推定し、伝搬路補償部ｂ１０６で、データ信号に対して振幅をｇ倍して補償することが出来る。すなわち受信したデータ信号をｇ倍することで正しく信号を検出できる。
　以上のように、本発明に係る非線形プレコーディング部ａ１３は電力の正規化も含めてデータ信号と同じ非線形プレコーディングをＤＭＲＳに対しても施す。
＜ＤＭＲＳ伝搬路推定部ｂ１２詳細＞
　基地局装置Ｂ１でＷ（＝Ｈ－１）をデータ信号に乗算しているので理想的な環境では受信信号がＨＷｓ＝ｓとなり、データ信号を伝搬路補償する必要が無い。しかし、１）基地局装置Ａ１で電力の正規化を行っており、２）端末装置Ｂ１ｎが伝搬路状態情報をフィードバックした時と基地局装置Ａ１がデータ信号を送信する時とで伝搬路状態が変動するため、あらためて、データ信号受信時に、データ信号と同じプレコーディング処理が施されたＤＭＲＳを使ってデータ信号の受信ゲイン（プレコーディングも含めた等価的な伝搬路の複素利得）を推定する必要がある。
　ＤＭＲＳ伝搬路推定部ｂ１２の詳細な構成を図７に示す。ＤＭＲＳ伝搬路推定部ｂ１２は、摂動ベクトル候補選択部ｂ１２１、摂動ベクトル加算部ｂ１２２、仮伝搬路推定部ｂ１２３、伝搬路補償部ｂ１２４、Ｍｏｄｕｌｏ演算部ｂ１２５、復調部ｂ１２６、摂動ベクトル評価値算出部ｂ１２７、および摂動ベクトル推定部ｂ１２８から構成される。
（ステップＳ２１）　まず、摂動ベクトル候補選択部ｂ１２１が、候補となる摂動ベクトル候補を選択して、選択した摂動ベクトルを摂動ベクトル加算部ｂ１２２に入力する。ここで候補となる摂動ベクトルＺａは、
　　　 Ｚａ＝（ｚ１ａ＋ｉｚ２ａ）τ
　　　　　　　　　　　　　　　　　　　　　　　 （１－３）
で表される。ここでｚ１ａとｚ２ａは、Ｎ次元縦ベクトルであり、各成分があらかじめ決められた範囲（たとえば絶対値がＬＭＳ以下の整数：［－ＬＭＳ，－ＬＭＳ＋１，．．．，－２，－１，０，１，２，．．．，ＬＭＳ－１，ＬＭＳ］）をとる。そのため式（１－３）の摂動ベクトルの候補は、（２ＬＭＳ＋１）２Ｎ個存在する。ここでは、基地局装置Ａ１と端末装置Ｂ１ｎで共通のＤＭＲＳ用のＭｏｄｕｌｏ幅を用いる。
　なお、式（１－３）において、基地局装置でも摂動ベクトルの候補を同様に制限したが、ここではＬＢＳ＝ＬＭＳであることが望ましいが、端末装置Ｂ１ｎの演算量低減のためＬＢＳ＞ＬＭＳとしても良い。
　また、第３の実施形態で後述するＴＨＰやＬＲ－ＴＨＰの場合においては、ＬＢＳを決めることが出来ないので、端末装置Ｂ１ｎの演算量によってＬＭＳを決めても良い。
（ステップＳ２２）　次に、摂動ベクトル加算部ｂ１２２は、基準信号（基地局装置から送信される信号であって、非線形プレコーディング前のＤＭＲＳ信号）ｑに摂動ベクトル候補Ｚａ＝（Ｚ１ａ＋ｉＺ２ａ）τを加算して信号ｑ＋Ｚａを算出する。摂動ベクトル加算部ｂ１２２は、信号ｑ＋Ｚａを仮伝搬路推定部ｂ１２３に入力する。
（ステップＳ２３）　次に、仮伝搬路推定部ｂ１２３は、信号ｑ＋Ｚａで受信したＤＭＲＳｐを除算する。すると複素利得は、ｈ＝ｐ／（ｑ＋Ｚａ）と推定できる。
　仮伝搬路推定部ｂ１２３は、複素利得ｈの振幅｜ｈ｜と位相ａｒｇ（ｈ）とを伝搬路補償部ｂ１２４、および、摂動ベクトル推定部ｂ１２８に入力する。
（ステップＳ２４）　伝搬路補償部ｂ１２４は、仮伝搬路推定部ｂ１２３から入力された振幅および位相を用いて、当該ＤＭＲＳが配置されたフレーム内のデータ信号を伝搬路補償する。つまり受信信号を複素利得ｈで除算する。伝搬路補償後のデータ信号をＭｏｄｕｌｏ演算部ｂ１２５に入力する。
（ステップＳ２５）　Ｍｏｄｕｌｏ演算部ｂ１２５は、伝搬路補償部ｂ１２４が伝搬路補償したデータ信号に対してＭｏｄｕｌｏ演算を行う。基地局装置Ａ１と端末装置Ｂ１ｎで共通のデータ信号用のＭｏｄｕｌｏ幅を用いる。Ｍｏｄｕｌｏ演算部ｂ１２５はＭｏｄｕｌｏ演算後のデータ信号を復調部ｂ１２６に入力する。
（ステップＳ２６）　復調部ｂ１２６は、入力されたＭｏｄｕｌｏ演算後のデータ信号に対して、軟推定を行い、軟推定値を摂動ベクトル評価値算出部ｂ１２７に入力する。
　なお、復調部ｂ１２６は軟推定値を対数尤度比（Ｌｏｇ　Ｌｉｋｅｌｉｈｏｏｄ　Ｒａｔｉｏ：ＬＬＲ）で算出するとする。ＬＬＲは次式で算出する。
Figure JPOXMLDOC01-appb-M000001
Figure JPOXMLDOC01-appb-M000002
Figure JPOXMLDOC01-appb-I000003
　復調部ｂ１２６は式（１－４）と式（１－５）に基づいて、各データ信号に割り当てられた各ビットに対応するＬＬＲを摂動ベクトル評価値算出部ｂ１２７に入力する。すなわち、ＱＰＳＫの場合はデータ信号の数×２個、１６ＱＡＭの場合はデータ信号の数×４個のＬＬＲを算出して、摂動ベクトル評価値算出部ｂ１２７に入力する。
（ステップＳ２７）　摂動ベクトル評価値算出部ｂ１２７は、入力されたＬＬＲの分散を算出し、算出したＬＬＲの分散および対応するＬＬＲを摂動ベクトル推定部ｂ１２８に入力する。なおＬＬＲの分散は下式で算出する。
Figure JPOXMLDOC01-appb-M000004
　ここで、Ｖは摂動ベクトルの候補毎に軟推定したデータ信号の個数である。またＭは各変調方式に割り当てられたビット数であり、ＱＰＳＫで２、１６ＱＡＭで４となる。またＬｖ ｍは第ｖ番目のデータ信号に割り当てられたｍ番目のビットのＬＬＲを示す。
（ステップＳ２８）　摂動ベクトル推定部ｂ１２８は、摂動ベクトル評価値算出部ｂ１２７から入力された各摂動ベクトル候補に対応するＬＬＲの分散の中で最も大きいものを選択し、最大の分散を有する摂動ベクトルに対応するデータ信号の振幅と位相（「等価伝搬路状態情報」と呼ぶ。）を出力し、ＤＭＲＳ伝搬路推定部ｂ１２の外の伝搬路補償部ｂ１０６（図６）に入力する。
　ＬＬＲの分散が大きいということは、各摂動ベクトル候補を仮定した時に、基地局装置Ａ１から伝わった相互情報量が最も大きいということを示す。ＬＬＲは、各ビットが１か０かという「確からしさ・確率」が高いほど、絶対値が大きくなる。そのため、各摂動ベクトル候補を仮定した中で、それぞれのＬＬＲの分散を算出し、ＬＬＲの分散が最も大きい摂動ベクトルを選ぶことで、その摂動ベクトルが最も「確からしさ・確率」が高い摂動ベクトルであると推定できる。
　尚、ＤＭＲＳ伝搬路推定部ｂ１２内の伝搬路補償部ｂ１２４、Ｍｏｄｕｌｏ演算部ｂ１２５、および復調部ｂ１２６は、ＤＭＲＳ伝搬路推定部ｂ１２外の伝搬路補償部ｂ１０６、Ｍｏｄｕｌｏ演算部ｂ１０９、および復調部ｂ１１０とそれぞれ同じ動作を行うので、回路規模低減のため、両構成部分を共通の回路としても良い。また、一度復調部ｂ１２６で算出した軟推定値のうち摂動ベクトル推定部ｂ１２８で推定された摂動ベクトルに対応する軟推定値を復号部ｂ１１１で用いることで、二重に軟推定値を算出することを避けて演算量を低減しても良い。
＜効果＞
　本実施形態で示したように、基地局装置がＤＭＲＳをＮＬＰ　ＭＵ－ＭＩＭＯにより、空間多重して送信することで、電力効率良く、ＤＭＲＳ挿入によるオーバーヘッドを低減できる。
　尚、摂動ベクトルの推定のために軟推定するデータ信号は、フレームの一部のデータ信号でもよい。また、ＬＬＲの分散以外の指標を用いて摂動ベクトルを推定してもよい。
（第２の実施形態：　２次元ユークリッドの互除法）
　第１の実施形態では、端末装置Ｂ１ｎが摂動ベクトルの候補から最尤な摂動ベクトルを推定した。これによりＤＭＲＳをＮＬＰ　ＭＵ－ＭＩＭＯによって空間多重することが可能となった。本実施形態では、本実施形態に係る基地局装置Ａ２と本実施形態に係る端末装置Ｂ２ｎの両装置が、「２次元ユークリッドの互除法」という処理を導入することで、端末装置Ｂ２ｎの演算量を第１の実施形態よりも低減する。
　また本実施形態に係る基地局装置Ａ２はＤＭＲＳ補正部ａ２２６を有すること以外、本実施形態に係る端末装置Ｂ２ｎにおいてはＤＭＲＳ伝搬路推定部ｂ２２の動作以外、第１の実施形態の基地局装置Ａ１、端末装置Ｂ１ｎと全く同じ構成を有するので、相違する部分を除き詳細な説明は省略する。図１５に基地局装置Ａ２の構成を示す。上述の通り、図３に対して、ＤＭＲＳ補正部ａ２２６が新たに追加されている。
　以降、基地局装置Ａ２のＤＭＲＳ補正部ａ２２６と端末装置Ｂ２ｎのＤＭＲＳ伝搬路推定部ｂ２２について詳細に説明する。
＜２次元ユークリッドの互除法＞
　本実施形態に係る基地局装置Ａ２および端末装置Ｂ２ｎは、両方とも「２次元ユークリッドの互除法」という処理を行う２次元ユークリッドの互除法部を有する。２次元ユークリッドの互除法部は、本実施形態の特徴部分の１つであるため、原理を先に説明する。
　２次元ユークリッドの互除法は、通常のユークリッドの互除法を２次元に拡張したアルゴリズムである。ここでは通常のユークリッドの互除法を１次元ユークリッドの互除法と呼ぶ。１次元ユークリッドの互除法は、異なる２つの整数の最大公約数を求めるアルゴリズムである。一例として１５と３６の最大公約数を求める。まず、１５と３６のうち大きい方から小さい方を減算する。すると、３６－１５＝２１となり、２１が得られる。次は１５と３６のうち小さい方である１５と２１を対象として、大きい方から小さい方を減算する。すると２１－１５＝６となる。この処理を順に繰り返す。すると、以下のように算出できる。
（３６，１５）⇒（２１，１５）⇒（１５，６）⇒（９，６）⇒（６，３）⇒（３，３）⇒（３，０）
　最後に０が出てきたところでアルゴリズムは終了する。残った数の０以外の数「３」が最大公約数となる。以上が１次元ユークリッドの互除法である。
　次に、これを複素数に拡張する。複素数に拡張したものが、２次元ユークリッドの互除法である。対象とするのは、実部および虚部が整数の複素数（ガウス整数）が２つ有る場合に、両ガウス整数の最大公約数に対応する「基底ベクトル（もしくは既約ベクトル）」を算出する方法である。一例として、ガウス整数３＋ｉと－１＋ｉに２次元ユークリッドの互除法を適用する例について説明する。
　まず、３＋ｉと－１＋ｉのノルム（絶対値）が大きいガウス整数と、＋１倍、－１倍、＋ｉ倍および－ｉ倍した小さいガウス整数を加算する（ガウス整数を－１倍、＋ｉ倍および－ｉ倍することは、それぞれガウス整数の位相を＋１８０度、＋９０度、＋２７０度回転させることに対応する。）。
　すると、
　　　 （３＋ｉ）＋（＋１）（－１＋ｉ）＝２＋２ｉ
　　　 （３＋ｉ）＋（－１）（－１＋ｉ）＝４
　　　 （３＋ｉ）＋（＋ｉ）（－１＋ｉ）＝２
　　　 （３＋ｉ）＋（－ｉ）（－１＋ｉ）＝４＋２ｉ
となる。この４つの値（２＋２ｉ、４、２、４＋２ｉ）の中で最もノルムが小さいガウス整数は２となる。次に、３＋ｉと－１＋ｉのノルムが小さい方と２を用いて同様の処理を繰り返す。
　　　 ２＋（＋１）（－１＋ｉ）＝１＋ｉ
　　　 ２＋（－１）（－１＋ｉ）＝３－ｉ
　　　 ２＋（＋ｉ）（－１＋ｉ）＝１－ｉ
　　　 ２＋（－ｉ）（－１＋ｉ）＝３＋ｉ
　この４つの値の中で最もノルムが小さいのは１＋ｉと１－ｉである。ノルムが同じものが２つ以上ある時は、どちらを選択しても良い。例えば１＋ｉを選択したとして、－１＋ｉと１＋ｉで同様の処理を行うと、
　　　 （－１＋ｉ）＋（＋１）（１＋ｉ）＝２ｉ
　　　 （－１＋ｉ）＋（－１）（１＋ｉ）＝－２
　　　 （－１＋ｉ）＋（＋ｉ）（１＋ｉ）＝－２＋２ｉ
　　　 （－１＋ｉ）＋（－ｉ）（１＋ｉ）＝０
となる（２つのガウス整数のノルムが等しい時はどちらを、＋１倍、－１倍、＋ｉ倍および－ｉ倍してもよい。）。
　よって、０が出てきたので、１＋ｉまたは－１＋ｉが１次元ユークリッドの互除法における最大公約数に相当するもの（「既約ベクトル」と呼ぶ。）であり、以上で２次元ユークリッドの互除法は終了となる。
　本実施形態は、前に述べたように２次元ユークリッドの互除法を使ってＤＭＲＳの伝搬路推定を行うことが特徴であるが、この伝搬路推定の原理を説明する。本実施形態では基準となるＤＭＲＳは、必ずＱＰＳＫ信号の各点のうち１つであるとする。
　図８に、２つのＤＭＲＳ（ＤＭＲＳ１およびＤＭＲＳ２と呼ぶ。）の位置の一例を信号点平面上に示す。ここで、基地局装置は、ＤＭＲＳ１を３＋ｉの位置で送信し、ＤＭＲＳ２が－１＋ｉの位置で送信したとする。また本信号点平面の単位は、ＱＰＳＫ信号の振幅が√２となるようにする。すなわち、ＱＰＳＫの各信号点は、（±１±ｉ）の点であり、Ｍｏｄｕｌｏ幅は４となる（便宜上、実施形態１と振幅の単位を変えている。）。また、ＤＭＲＳ１およびＤＭＲＳ２を示す●とそれ以外の点○はすべてＱＰＳＫの信号点、およびＱＰＳＫの信号点に対して任意の摂動ベクトルが加算された点であり、図８に示したように格子点上に整列した点となっている。ＤＭＲＳ１は、ＱＰＳＫ信号の－１＋ｉに摂動ベクトル＋４が加算された信号である。
　このような信号点平面において、３＋ｉと－１＋ｉに対して、先ほどの２次元ユークリッドの互除法を実行すれば、「最大公約数」の複素数１＋ｉを得ることが出来る。これは、図８に示した格子点群（○および●）の最小の格子ベクトル（既約ベクトル）を算出したことを意味する。そして、この既約ベクトルルはＱＰＳＫの４点のうちいずれかに相当する。また、最初の信号がａ（３＋ｉ）とａ（－１＋ｉ）というように、「ａ」という複素数の共通因数が乗算されている場合は、既約ベクトルもａ倍されて、ａ（１＋ｉ）という既約ベクトルを得ること出来る。
　ここで、端末装置においては、ＤＭＲＳ１およびＤＭＲＳ２は伝搬路の複素利得ｈが乗算されたｈ（３＋ｉ）とｈ（－１＋ｉ）である。ここで雑音は無視している。この２つの信号に対して、前述の２次元ユークリッドの互除法を適用すると既約ベクトルｈ（１＋ｉ）が得られる。これは、摂動ベクトルが加算されていないＱＰＳＫ信号の４点のうち１点に伝搬路の複素利得が乗算された信号である。ＱＰＳＫの各点はノルムが√２であることが分かっているので、伝搬路の複素利得の絶対値｜ｈ｜を求めることが出来る。またａｒｇ（ｈ）も求めることが出来るので、伝搬路の複素利得ｈの位相を求めることが出来る。
　ただし、端末装置はＤＭＲＳ伝搬路推定時に２次元ユークリッドの互除法で得た信号ｈ（１＋ｉ）は、ＱＰＳＫの４点（±１±ｉ）のいずれにｈが乗算されたものか判断できないので、仮に（１＋ｉ）とおいて位相ａｒｇ（ｈ）を求める。もし、他の３点（－１＋ｉ）、（１－ｉ）および（－１－ｉ）であったとしても、実際の複素利得の位相差は必ず＋９０度、＋１８０度または＋２７０度のうちいずれかとなる。
　次に、ＤＭＲＳ１とＤＭＲＳ２に対してＭｏｄｕｌｏ演算を行う。ここでＭｏｄｕｌｏ演算は、Ｉ－ｃｈ、Ｑ－ｃｈともに同様の処理を行い、各軸に線対称な処理であって、９０度の回転に対称な演算である。そのためａｒｇ（ｈ）が９０度×整数倍だけ回転しても問題無く、ＤＭＲＳ１とＤＭＲＳ２に対してＭｏｄｕｌｏ演算を施すことが出来る。
　この場合、Ｍｏｄｕｌｏ演算後のＤＭＲＳ１はｈ（１＋ｉ）に、ＤＭＲＳ２はｈ（－１＋ｉ）となる。雑音μ１、μ２を考慮すれば、ＤＭＲＳ１はｈ（１＋ｉ）＋μ１、ＤＭＲＳ２はｈ（－１＋ｉ）＋μ２となる。ここで端末装置は、非線形プレコーディング前、ＤＭＲＳ１は（１＋ｉ）、ＤＭＲＳ２は（－１＋ｉ）という形となっていることを知っているので、各ＤＭＲＳにより、
　　　 ｛ｈ（１＋ｉ）＋μ１｝／（１＋ｉ）＝ｈ＋μ１／（１＋ｉ）
　　　 ｛ｈ（－１＋ｉ）＋μ２｝／（－１＋ｉ）＝ｈ＋μ２／（－１＋ｉ）
という形で伝搬路推定をすることができる。μ１／（１＋ｉ）とμ２／（－１＋ｉ）は伝搬路推定誤差であり、ＤＭＲＳ１とＤＭＲＳ２の伝搬路推定結果を最大比合成することで伝搬路推定誤差を低減できる。
　ここで、注意すべきなのが、一部のＤＭＲＳの組み合わせにおいては、２次元ユークリッドの互除法によって、既約ベクトルに対応するｈ（±１±ｉ）を算出することができないことがある点である。ＤＭＲＳが図９のような配置になったときは、既約ベクトルが３＋３ｉとなり、ＱＰＳＫの信号点と一致しない。このような状態を回避すべく、基地局装置においても、２次元ユークリッドの互除法を用いて、図９のような場合が起こることを判断し、信号点を変えて問題のない２点を送信する。
　以上が、本実施形態の２次元ユークリッド互除法部の処理の概要である。これらの処理を、２次元ユークリッドの互除法部の構成を示した図１０を用いて再度説明する。前述のように、２次元ユークリッドの互除法部３００は、基地局装置と端末装置との両方に存在するが、まず、端末装置の２次元ユークリッドの互除法部について説明する。
　図１０に示した２次元ユークリッドの互除法部３００は、ベクトル保存部３０１、差分ベクトル算出部３０３、差分ベクトルノルム算出部３０５、収束判定部３０７、ノルム算出部３０９からなる。
（ステップＳ３１）　まず、ベクトル保存部３０１は、入力された２つのＤＭＲＳをノルム算出部３０９に入力する。ここで２つのＤＭＲＳには、基地局装置から送信されたＤＭＲＳに、伝搬路の複素利得ｈが乗算され、さらに雑音が加算されている。
（ステップＳ３２）　ノルム算出部３０９は、入力された２つのＤＭＲＳのノルムをそれぞれ算出し、各ノルムをベクトル保存部３０１に入力する。
（ステップＳ３３）　ベクトル保存部３０１は、入力された２つのＤＭＲＳと対応するノルムを差分ベクトル算出部３０３に入力する。
（ステップＳ３４）　差分ベクトル算出部３０３は、入力された２つのベクトル（以降、２つのＤＭＲＳに対して再帰的に差分ベクトル算出部３０３で算出した信号を「ベクトル」と呼ぶ。）のうちノルムが大きいベクトルをａとし、小さいベクトルをｂとして、ｃ１=ａ＋ｂ, ｃ２=ａ－ｂ, ｃ３＝ａ＋ｉｂ，ｃ４＝ａ－ｉｂの４つのベクトルを算出し、差分ベクトルノルム算出部３０５に入力する。またベクトルｂとベクトルｂのノルムをベクトル保存部３０１に入力する。
（ステップＳ３５）　差分ベクトルノルム算出部３０５は、入力された４つのベクトルｃ１～ｃ４のノルムを算出し、最もノルムが小さいベクトルとそのベクトルに対応するノルムをベクトル保存部３０１に入力し、最もノルムが小さいベクトルのノルムを収束判定部３０７に入力する。
（ステップＳ３６）　収束判定部３０７は、入力されたノルムと所定の正の数Ｔとの大小関係を比較する。Ｔよりノルムが大きい場合、入力されたノルムとベクトルをベクトル保存部３０１に入力する。Ｔよりノルムが小さい場合、ベクトル保存部３０１に、２次元のユークリッドの互除法が終了したことを示す情報を入力する。
　ここで、前述の２次元ユークリッド互除法の原理のように、ノルムが０になったときを収束条件としないのは、端末装置側では、ＤＭＲＳが雑音を含むことから、最終的に、収束判定部３０７のノルムが０にならないからである。しかし、雑音は、ＤＭＲＳよりも通常小さいので、所定の定数Ｔを用いて、ノルムがＴより小さくなったときに収束したと判定する。
　Ｔの値は、次の２つの要素のトレードオフで、事前にコンピュータシミュレーションによってあらかじめ決定しておく。
１）Ｔが大きくなると、収束判定条件が緩くなり、雑音による誤差が無い時に０でないのに収束していると判定してしまう確率が上昇する。
２）Ｔが小さくなると、収束判定条件が厳しくなり、雑音による誤差のために、収束したと判定すべきであるにもかかわらず、まだ収束していないと判定してしまうことがある。
　そのため、この２つの要素のトレードオフで事前に決めておく。
（ステップＳ３７）　ベクトル保存部３０１は、
１）差分ベクトル算出部３０３から入力されたベクトルと対応するノルム、
２）差分ベクトルノルム算出部３０５から入力されたベクトルと対応するノルム
の２つを差分ベクトル算出部３０３に入力し、再びステップＳ３４から処理を繰り返す。
（ステップＳ３８）　ベクトル保存部３０１は、収束判定部３０７から、２次元のユークリッドの互除法が終了したことを示す情報が入力された場合、ベクトルｂ（既約ベクトル）を出力する。
　また、基地局装置における２次元ユークリッドの互除法部３００は、ＤＭＲＳに対して伝搬路の複素利得ｈが乗算されておらず、雑音も加算されていないので収束判定部３０７は判定条件としているＴ＝０とする。
＜ＤＭＲＳ補正部ａ２２６詳細＞
　次に、基地局装置Ａ２において２次元ユークリッドの互除法を行うＤＭＲＳ補正部ａ２２６の詳細な構成と動作を説明する。図１１は、ＤＭＲＳ補正部ａ２２６の構成を示したものである。ＤＭＲＳ補正部ａ２２６は、２次元ユークリッドの互除法部３００、既約ベクトル確認部３２０、および摂動ベクトル加算部３４０、からなる。
　このＤＭＲＳ補正部ａ２２６は、前述の、既約ベクトルが最小の格子ベクトルと一致しない問題を、基地局装置Ａ２で処理することで、あらかじめ回避することを目的とする。
　ＤＭＲＳ補正部ａ２２６は、まず非線形プレコーディング部ａ１３から摂動ベクトルが加算されたＤＭＲＳを取得する。ここで、ＤＭＲＳ補正部ａ２２６が取得するＤＭＲＳは、摂動ベクトルが加算されている。フィルタＷを加算する前の信号（ｑ＋ｚ１τ＋ｉｚ２τ）とする。また、取得するＤＭＲＳは１度に２つであり、この２つは、あらかじめ２つずつ対になっており、どの２つを対にするかは基地局装置Ａ２と端末装置Ｂ２ｎとで既知である。
　次に２次元ユークリッドの互除法部３００が、この２つのＤＭＲＳに対して２次元ユークリッドの互除法を適用して得た既約ベクトル（「ｄ」とおく。）と非線形プレコーディング部ａ１３から取得した２つのＤＭＲＳを既約ベクトル確認部３２０に入力する。
　既約ベクトル確認部３２０は、既約ベクトルｄが摂動ベクトル加算前のＤＭＲＳ信号ｑすなわちＱＰＳＫ信号のうちの１つ（±１±ｉ）であれば、非線形プレコーディング部ａ１３から取得した２つのＤＭＲＳを再び非線形プレコーディング部ａ１３に入力する。
　一方、既約ベクトル確認部３２０は、既約ベクトルｄが摂動ベクトル加算前のＤＭＲＳすなわちＱＰＳＫ信号のうちの１つ（±１±ｉ）でない場合、非線形プレコーディング部ａ１３から取得したＤＭＲＳを２つとも摂動ベクトル加算部３４０に入力する。
　摂動ベクトル加算部３４０は、非線形プレコーディング部ａ１３から取得したＤＭＲＳのうちどちらか一方に対して摂動ベクトルτ（または－τ、ｉτ、－ｉτ）を加算する。ここで、どの摂動ベクトルτをどちらのＤＭＲＳに加算するかは、乱数などを用いてランダムに選択する。摂動ベクトル加算部３４０は、摂動ベクトルを加算したＤＭＲＳを含む２つのＤＭＲＳを再び２次元ユークリッド互除法部３００に入力する。
　２次元ユークリッド互除法部３００は、新たなＤＭＲＳを用いて再び２次元ユークリッドの互除法を行う。以降、既約ベクトル確認部３２０が、既約ベクトルｄが摂動ベクトル加算前のＤＭＲＳすなわちＱＰＳＫ信号のうちの１つ（±１±ｉ）であると判定するまで、２次元ユークリッドの互除法部３００、既約ベクトル確認部３２０、摂動ベクトル加算部３４０の処理を続ける。
　最終的に、既約ベクトル確認部３２０が、補正後のＤＭＲＳを２つ出力し、非線形プレコーディング部ａ１３に入力する。
　また、非線形プレコーディング部ａ１３は、ＤＭＲＳ補正部ａ２２６でＤＭＲＳが補正された場合、新たに摂動ベクトルを探索することなく、当該ＤＭＲＳにフィルタＷを乗算し、電力正規化係数ｇの逆数ｇ－１を乗算する。
＜ＤＭＲＳ伝搬路推定部ｂ１０７詳細＞
　次に、端末装置ｂ２に係るＤＭＲＳ伝搬路推定部ｂ１２の詳細な構成と動作とを説明する。図１２（ａ）は、ＤＭＲＳ伝搬路推定部ｂ１２の詳細な構成を示したものである。ＤＭＲＳ伝搬路推定部ｂ１２は、２次元ユークリッドの互除法部３００と複素利得算出部３５０から構成される。
　また複素利得算出部３５０は，図１２（ｂ）に示したように、内部に仮複素利得算出部３５１、ＤＭＲＳ伝搬路補償部３５２、ＤＭＲＳ－Ｍｏｄｕｌｏ部３５３、およびベクトル除算部３５４を有する。
　まず、ＤＭＲＳ伝搬路推定部ｂ１２は、端末装置Ｂ２ｎが受信した２つのＤＭＲＳを取得する。この２つは、基地局装置におけるＤＭＲＳ補正部ａ２２６でも説明した対となるＤＭＲＳ（ここではｐ１とｐ２とおく。）である。
　この２つのＤＭＲＳに対して２次元ユークリッドの互除法部３００で２次元ユークリッドの互除法を行い、既約ベクトルｐｒｅｄを算出する。なお、既約ベクトルｐｒｅｄは、複素利得ｈが乗算されているため、ＤＭＲＳ補正部ａ２２６で算出した既約ベクトルｄとは異なる。ここで前述のように、既約ベクトルｐｒｅｄは、摂動ベクトルが加算されていないＱＰＳＫ信号の４点のうち１点に伝搬路の複素利得ｈが乗算され、さらに雑音による誤差を加算した信号である。
　仮複素利得算出部３５１は、この既約ベクトルｐｒｅｄが、ＱＰＳＫ信号の４点のうち（１＋ｉ）に複素利得が乗算されたものと仮定して、位相θ（　＝ａｒｇ（ｐｒｅｄ／（１＋ｉ））　）を求める（前述のように、もし、他の（－１＋ｉ）、（１－ｉ）および（－１－ｉ）であったとしても、実際の複素利得の位相差は必ず＋９０度、－９０度または＋１８０度となる。）。ＱＰＳＫの各点はノルムが√２であることが分かっているので、伝搬路の複素利得の絶対値｜ｈ｜（＝ｐｒｅｄ／√２）も求めることが出来る。仮複素利得算出部３５１は、ここで算出した伝搬路の複素利得の絶対値｜ｈ｜と位相θ（　＝ａｒｇ（ｐｒｅｄ／（１＋ｉ））　）をＤＭＲＳ伝搬路補償部３５２に入力する。
　ＤＭＲＳ伝搬路補償部３５２は、入力された伝搬路の複素利得の絶対値｜ｈ｜と位相θとを用いて各ＤＭＲＳｐ１とｐ２の伝搬路補償を行う。具体的には、ＤＭＲＳの振幅を１／｜ｈ｜倍し、位相を－θだけ回転する。つまり、
ｑ１＝ｐ１／｜ｈ｜×ｅｘｐ（－２πｉθ）
ｑ２＝ｐ２／｜ｈ｜×ｅｘｐ（－２πｉθ）
とする。ＤＭＲＳ伝搬路補償部３５２は、伝搬路補償後のＤＭＲＳｑ１とｑ２、伝搬路の複素利得の絶対値｜ｈ｜、および位相θをＤＭＲＳ－Ｍｏｄｕｌｏ部３５３に入力する。
ＤＭＲＳ－Ｍｏｄｕｌｏ部３５３は、入力された伝搬路補償後の各ＤＭＲＳｑ１とｑ２に対してＭｏｄｕｌｏ演算を行って信号ｅ１、ｅ２を算出する。ＤＭＲＳは振幅・位相ともに伝搬路補償されているため、通常のＱＰＳＫ信号に対応するＭｏｄｕｌｏ幅τ（＝２√２）を用いたＭｏｄｕｌｏ演算を行うことができる。ＤＭＲＳ－Ｍｏｄｕｌｏ部３５３はＭｏｄｕｌｏ演算を行った信号ｅ１、ｅ２に対して、ＤＭＲＳ伝搬路補償部３５２から入力された複素利得の絶対値｜ｈ｜と位相θを用いて、信号ｐｍｏｄ１（　＝｜ｈ｜×ｅ１×ｅｘｐ（２πｉθ）　）およびｐｍｏｄ２（　＝｜ｈ｜×ｅ２×ｅｘｐ（２πｉθ）　）を算出する。その後信号ｐｍｏｄ１とｐｍｏｄ２をベクトル除算部３５４に入力する。
　ベクトル除算部３５４は、信号ｐｍｏｄ１とｐｍｏｄ２を、端末装置が受信したＤＭＲＳをあらかじめ端末装置が把握している各ＤＭＲＳの基準信号ｑ（すなわち基地局装置で摂動ベクトルが加算される前のＤＭＲＳの信号点であり、ここではＱＰＳＫの信号点４点のうちいずれか１つ。）で除算することで伝搬路の複素利得ｈ（＝ｐｍｏｄ１／ｑまたはｐｍｏｄ２／ｑ）を推定する。ここでＤＭＲＳが２つあるため、同じ伝搬路の複素利得ｈに対する推定値が２つ得られるが、最大比合成することで雑音の影響による誤差を抑えてもよい。ベクトル除算部３５４は、推定した伝搬路の複素利得ｈを、伝搬路補償部ｂ１０６に入力する。
＜効果＞
　本実施形態に示したように基地局装置Ａ２と端末装置Ｂ２ｎでそれぞれ２次元ユークリッド互除法に基づく処理を行うことにより、比較的少ない端末装置Ｂ２ｎの演算量によってＤＭＲＳをＮＬＰ　ＭＵ－ＭＩＭＯによる空間多重することが可能となる。これによりＤＭＲＳを挿入することによるオーバーヘッドを低減することが可能となる。
（第３の実施形態：　ＶＰ⇒ＴＨＰなど）
　第１の実施形態、第２実施形態に係る基地局装置Ｂ１ｎ、Ｂ２ｎは、いずれも非線形プレコーディング部ａ１３およびフィルタ算出部ａ１１でＶＰを行っていた。本実施形態は、ＶＰより演算量の低いＴＨＰを用いた方式である。本実施形態に係る基地局装置を基地局装置Ａ３とし、本実施形態に係る端末装置を端末装置Ｂ３１～Ｂ３Ｎと呼び、任意の１つを端末装置Ｂ３ｎと呼ぶ。
　本実施形態に係る非線形プレコーディング部を非線形プレコーディング部ａ３３と呼び、本実施形態に係るフィルタ算出部をフィルタ算出部ａ３１と呼ぶ。
　まず、フィルタ算出部ａ３１では、実施形態１と同じ伝搬路行列Ｈを算出した後、ＨをＱＲ分解する。
　　　 ＨＨ＝ＱＲ
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－１）
　ここでＲは上三角行列、Ｑはユニタリ行列である。またＨを行列の共役転置とする。ＲＨの対角成分だけ取り出した対角行列をＡとおく。ＡとＱを用いて線形フィルタＰ０を下式により算出する。
　　　 Ｐ０＝ＱＡ－１
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－２）
　また、ＡとＲを用いて干渉係数フィルタＦを下式により算出する。
　　　 Ｆ＝ＲＨＡ－１－Ｉ
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－３）
　ここで、ＩはＮ行Ｎ列の単位行列とする。最後に、フィルタ算出部ａ３１は線形フィルタＰ０と干渉係数フィルタＦを非線形プレコーディング部ａ３３に入力する。
　なお、電力正規化を行うために、フィルタ算出部ａ３１で、
　　　　　　　ｇ＝ｔｒａｃｅ｛（（Ｐ０ＣｖＰ０ Ｈ）／Ｐｔｒ）１／２｝
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－４）
　　　　　　　ｇ＝｛ｔｒａｃｅ（（Ｐ０ＣｖＰ０ Ｈ）／Ｐｔｒ）｝１／２
により電力正規化係数ｇを算出する。ここで、Ｐｓを変調信号の平均電力、ＰｖをＭｏｄｕｌｏ演算後のデータ信号の平均電力とすると、Ｃｖは対角成分が左上から順に
　　　　　　　［Ｐｓ，Ｐｖ，Ｐｖ，…Ｐｖ］
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－５）
となる対角行列である。なお、ＰｖはＭｏｄｕｌｏ幅τによって変化し、Ｐｓ＝１とすると、ＱＰＳＫ（τ＝２√２）では４／３、１６ＱＡＭ（τ＝８／√１０）では１６／１５、６４ＱＡＭ（τ＝１６／√４２）では６４／６３となる。これは、Ｍｏｄｕｌｏ演算後のデータ信号が原点を中心とするＭｏｄｕｌｏ幅の範囲に、一様分布することを利用したものである。フィルタ算出部ａ３１は、式（３－２）で算出した線形フィルタＰ０にｇ－１を乗算した行列を新たに線形フィルタＰとして非線形プレコーディング部ａ３３に入力する。このときは第１の実施形態に示したような電力の正規化は行わない。
　非線形プレコーディング部ａ３３の構成を、図１３に、その動作を示すフローチャートを図１４に示した。以下、順番に非線形プレコーディング部ａ３３内の動作を説明する。
（ステップＳ１）　干渉算出部ａ１３１と線形フィルタ乗算部ａ１３４は、それぞれ干渉係数フィルタＦと線形フィルタＰとをフィルタ算出部ａ３１から取得する。　
（ステップＳ２）　送信信号を計算中の端末装置の番号を示す番号ｎに１を代入する。
（ステップＳ３）　端末装置Ｂ３１宛の固有信号ｓ１をｖ１とおく。
（ステップＳ４）　ｎにｎ＋１を代入する。すなわちｎ＝２とする。
（ステップＳ５）　干渉算出部ａ１３１は、ｖ１を用いて端末装置Ｂ３２が受ける干渉信号ｆ２を下式により算出する。
　　　 ｆ２＝Ｆ（２，１）＊ｖ１
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－６）
　ここで、Ｆ（ｐ，ｑ）は行列Ｆの２行１列成分を表す。
（ステップＳ６）　干渉減算部ａ１３３－２は端末装置Ｂ３２宛の固有信号ｓ２からｆ２を減算し、信号ｓ２－ｆ２を算出する。
（ステップＳ７）　第１のＭｏｄｕｌｏ演算部ａ１３２－２がｓ２－ｆ２に対してＭｏｄｕｌｏ演算を適用し、信号ｖ２を算出する。
（ステップＳ８）　ｎ＝２なので、ステップＳ４から再び次の端末装置Ｂ３３宛の信号の演算（ステップＳ４～ステップＳ７）を行う。
　以降、ステップＳ４～Ｓ８の処理をｎ＝Ｎになるまで繰り返すので、一例としてｎ番目の端末装置Ｂ３ｎ宛の信号を算出する処理を説明する。
（ステップＳ４）　ｎの値を１増やす。
（ステップＳ５）　干渉算出部ａ１３１はｖ１～ｖ（ｎ－１）を用いて端末装置ｎが受ける干渉信号ｆｎを下式により算出する。
　　　 ｆｎ＝Ｆ（ｎ，１：ｎ－１）＊［ｖ１，ｖ２，．．．，ｖ（ｎ－１）］ｔ
　　　　　　　　　　　　　　　　　　　　　　　　　　　（３－７）
　ここで、Ｆ（ｎ，１：ｎ－１）は、行列Ｆのｎ行目の１～ｎ－１列目の成分を示す横ベクトルを示す。
（ステップＳ６）　第ｎの干渉減算部ａ１３３－ｎは、第ｎの端末装置Ｂ３ｎ宛の固有信号ｓｎからｆｎを減算し、信号ｓｎ－ｆｎを算出する。
（ステップＳ７）　Ｍｏｄｕｌｏ演算部ａ１３２－ｎが、ｓｎ－ｆｎに対してＭｏｄｕｌｏ演算を適用し、信号ｖｎを算出する。このＭｏｄｕｌｏ演算によって各端末装置宛の送信信号電力を低減する。
（ステップＳ８）　ｎ＜Ｎのときは、再びステップＳ４を行う。またｎ＝ＮのときはステップＳ９に進む。
（ステップＳ９）　信号ｖ＝（ｖ１，ｖ２，．．．，ｖＮ）ｔに線形フィルタＰを乗算して得た信号をｘとおく。ここで第１の実施形態と同様に信号ｘの各成分は、順番に各アンテナａ１０１－１～ａ１０１－Ｎで送信する送信信号である。この信号ｘをフレーム構成部ａ１４２に入力する。
　ＴＨＰは、ＶＰと同様に摂動ベクトル探索アルゴリズムと考えることが出来ることを説明する。ただし、ＴＨＰの場合は、真に最適な摂動ベクトルを探索するアルゴリズムではなく，ＶＰと比較して少量の演算量で準最適な摂動べクトルを探索するアルゴリズムである。
　まず、非線形プレコーディング部ａ３３のＭｏｄｕｌｏ演算部ａ１３２－２～ａ１３２－Ｎが無いものと仮定して信号を演算すると、電力の正規化係数が１であるとすれば、非線形プレコーディング部ａ３３は、
　　　 ｘ＝Ｐ０（Ｉ＋Ｆ）－１ｓ
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－８）
という演算を行っている。ここでＰ０とＦに式（３－２）と（３－３）を代入すると、
　　　 ｘ＝ＱＡ－１｛Ｉ＋（ＲＨＡ－１－Ｉ）｝－１ｓ
　　　　 ＝ＱＡ－１（ＲＨＡ－１）－１ｓ
　　　　 ＝Ｑ（ＲＨ）－１ｓ＝（ＲＨＱＨ）－１ｓ｛（ＱＲ）Ｈ｝－１ｓ 
　　　　 ＝｛（ＨＨ）Ｈ｝－１ｓ　＝Ｈ－１ｓ＝Ｗｓ
　　　　　　　　　　　　　　　　　　　　　　　　　　　 （３－９）
ｘ＝ＱＡ－１｛Ｉ＋（ＲＨＡ－１－Ｉ）｝－１ｓ
　　　　 ＝ＱＡ－１（ＲＨＡ－１）－１ｓ
　　　　 ＝Ｑ（ＲＨ）－１ｓ＝（ＲＨＱＨ）－１ｓ＝｛（ＱＲ）Ｈ｝－１ｓ 
　　　　 ＝｛（ＨＨ）Ｈ｝－１ｓ＝Ｈ－１ｓ＝Ｗｓ
となり、式（１－２）において（ｚ１，ｚ２）＝（０，０）とした場合に対応する。
　さらに、Ｍｏｄｕｌｏ演算部ｓ１３２－２～ｓ１３２－Ｎは、固有信号ｓにＭｏｄｕｌｏ幅の整数倍の信号をＩ－ｃｈ又はＱ－ｃｈに加算する処理である。そのため、本実施形態に係る非線形プレコーディング部ａ３３で算出した信号ｘもｘ＝Ｗ（ｓ＋ｚ１τ＋ｉｚ２τ）という式で表すことが出来る。非線形プレコーディング部ａ３３で算出した信号ｘは必ずしも式（１－２）で表される最適な摂動ベクトルを加算した信号ではないが、Ｍｏｄｕｌｏ演算による電力抑圧を行っているため準最適な信号ｘとなり、Ｍｏｄｕｌｏ演算を行わない場合Ｗｓと比較して必ずノルムが小さくなる。
　以上、説明したようなアルゴリズムで信号ｘを算出することで、第１の実施形態のＶＰのように、候補となる摂動ベクトルを全探索する必要が無くなるため、演算量を低減することが出来る。これによりＮＬＰ　ＭＵ－ＭＩＭＯシステムにおいてＤＭＲＳを空間多重してＤＭＲＳ挿入によるオーバーヘッドの増加を最小限に抑えながら、端末装置において正常にデータ信号の複素利得を推定することが可能になる。
　なお、ここではＴＨＰを用いた準最適な摂動ベクトル探索アルゴリズムを用いる場合について説明したが、オーダリングを用いたＴＨＰ、又はＬＲ－ＴＨＰでも同様に適用可能である。
　上記の実施の形態において、添付図面に図示されている構成等については、これらに限定されるものではなく、本発明の効果を発揮する範囲内で適宜変更することが可能である。その他、本発明の目的の範囲を逸脱しない限りにおいて適宜変更して実施することが可能である。
　また、本実施の形態で説明した機能を実現するためのプログラムをコンピュータ読み取り可能な記録媒体に記録して、この記録媒体に記録されたプログラムをコンピュータシステムに読み込ませ、実行することにより各部の処理を行ってもよい。尚、ここでいう「コンピュータシステム」とは、ＯＳや周辺機器等のハードウェアを含むものとする。
　また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
　また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ－ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また前記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
　本発明は、通信装置に利用可能である。
Ａ１…基地局装置、Ｂ１１～Ｂ１４…端末装置、ａ１１…フィルタ算出部、ａ１３…非線形プレコーディング部、ａ３１…フィルタ算出部、ａ３３…非線形プリコーディング部、ａ１０１…アンテナ、ａ１０２…受信部、ａ１０３…ＧＩ除去部、ａ１０４…ＦＦＴ部、ａ１０５…伝搬路状態情報取得部、ａ１２１…符号化部、ａ１２２…変調部、ａ１２４…ＤＲＭＳ生成部、ａ１２５…固有信号構成部、ａ１３３…干渉減算部、ａ１３２…Ｍｏｄｕｌｏ演算部、ａ１３４…線形フィルタ乗算部、ａ１４１…ＣＲＳ生成部、ａ１４２…フレーム構成部、ａ１４３…ＩＦＦＴ部、ａ１４４…ＧＩ挿入部、ａ１４５…送信部、ａ２２６…ＣＲＳ補正部、ｂ１２…ＤＭＲＳ用伝搬路推定部、ｂ１０１…アンテナ、ｂ１０２…受信部、ｂ１０３…ＧＩ除去部、ｂ１０４…ＦＦＴ部、ｂ１０５…信号分離部、ｂ１０６…伝搬路補償部、ｂ１０７…ＣＲＳ用伝搬路推定部、ｂ１０８…伝搬路状態情報生成部、ｂ１０９…Ｍｏｄｕｌｏ演算部、ｂ１１０…復調部、ｂ１１１…復号部、ｂ１２１…摂動ベクトル候補選択部、ｂ１２２…摂動ベクトル加算部、ｂ１２３…仮伝搬路推定部、ｂ１２４…伝搬路補償部、ｂ１２５…Ｍｏｄｕｌｏ演算部、ｂ１２６…復調部、ｂ１２７…摂動ベクトル評価値算出部、ｂ１２８…摂動ベクトル推定部、ｂ１３１…ＩＦＦＴ部、ｂ１３２…ＧＩ挿入部、ｂ１３３…送信部、３００…２次元ユークリッドの互除去部、３０１…ベクトル保存部、３０３…差分ベクトル算出部、３０５…差分ベクトルノルム算出部、３０７…収束判定部、３０９…ノルム算出部、３２０…既約ベクトル確認部、３４０…摂動ベクトル加算部、３５０…複素利得算出部、３５１…仮複素利得算出部、３５２…ＤＭＲＳ伝搬路補償部、３５３…ＤＭＲＳ－Ｍｏｄｕｌｏ部、３５４…ベクトル除算部。
　本明細書で引用した全ての刊行物、特許および特許出願をそのまま参考として本明細書にとり入れるものとする。

Claims (18)

  1. 　複数の受信装置に対して同一時刻・同一周波数でデータ信号を送信する送信装置であって、
    　第１の受信装置宛の復調用参照信号と、前記第１の受信装置とは異なる第２の受信装置宛の復調用参照信号とを、同一時刻・同一周波数で送信する送信部、
    を有することを特徴とする送信装置。
  2. 　前記復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算する非線形プレコーディング部を有することを特徴とする請求項１に記載の送信装置。
  3. 　前記復調用参照信号に対して、非線形プレコーディングを施す非線形プレコーディング部を有することを特徴とする請求項１または請求項２に記載の送信装置。
  4. 　前記送信部は、
    　同一時刻・同一周波数で複数のデータ信号を送信し、
    　前記非線形プレコーディング部は、
    　前記復調用参照信号と、前記データ信号とに、同じ非線形プレコーディング処理することを特徴とする請求項３に記載の送信装置。
  5. 　前記非線形プレコーディング部は、
    　前記復調用参照信号と、前記データ信号とに、同じ電力の正規化を施すことを特徴とする請求項４に記載の送信装置。
  6. 　復調用参照信号を補正するＤＭＲＳ補正部を有することを特徴とする請求項１から５までのいずれか１項に記載の送信装置。
  7. 　前記ＤＭＲＳ補正部は、
    　前記復調用参照信号に２次元ユークリッド互除法を適用する２次元ユークリッド互除法部を有することを特徴とする請求項６に記載の送信装置。
  8. 　前記２次元ユークリッド互除法部は、
    第１の前記復調用参照信号から、
    第２の前記復調用参照信号、第２の前記復調用参照信号の位相を９０度回転させた信号、第２の前記復調用参照信号の位相を１８０度回転させた信号、および第２の前記復調用参照信号の位相を２７０度回転させた信号を減算する差分ベクトル算出部を有することを特徴とする請求項７に記載の送信装置。
  9. 　復調用参照信号にあらかじめ決めた所定の幅の整数倍の信号を加算する摂動ベクトル加算部を有することを特徴とする受信装置。
  10. 　前記信号を加算した復調用参照信号を用いて伝搬路推定する仮伝搬路推定部を有することを特徴とする請求項９に記載の受信装置。
  11. 　複数の異なる前記信号を選択する摂動ベクトル候補選択部と、
    　前記信号それぞれを用いて前記仮伝搬路推定部で伝搬路推定した伝搬路推定結果に基づいて、１つの前記信号を選択する摂動ベクトル推定部を有することを特徴とする請求項１０に記載の受信装置。
  12. 　複数の異なる前記信号に対応する前記伝搬路推定結果それぞれに基づいてデータ信号を軟推定して対数尤度比をそれぞれ算出する復調部と、
    　前記対数尤度比それぞれの分散を算出する摂動ベクトル評価値算出部と、
    を有し、
    　前記摂動ベクトル推定部は、前記分散のうち最も大きいものに対応する前記信号を選択することを特徴とする請求項１０に記載の受信装置。
  13. （追加）
    　　複数の前記復調用参照信号に対して２次元ユークリッドの互除法を適用し、既約ベクトルを算出する２次元ユークリッドの互除法部を有することを特徴とする請求項９に記載の受信装置
  14. 　前記既約ベクトルを用いて伝搬路の複素利得を算出する複素利得算出部を有する請求項１３に記載の受信装置。
  15. 　復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算する非線形プレコーディング部と、復調用参照信号と他の信号を、同一時刻・同一周波数で送信する送信部とを有する送信装置と、
    　復調用参照信号にあらかじめ決めた所定の幅の整数倍の信号を加算する摂動ベクトル加算部を有する受信装置と、
    を有することを特徴とする通信システム。
  16. 　復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算するステップと、復調用参照信号と他の信号を、同一時刻・同一周波数で送信するステップを有する送信方法と、
    　復調用参照信号にあらかじめ決めた所定の幅の整数倍の信号を加算するステップを有する受信方法と、
    を有することを特徴とする通信方法。
  17. 　復調用参照信号に対して、あらかじめ決めた所定の信号の整数倍の信号を加算する非線形プレコーディング部と、復調用参照信号と他の信号を、同一時刻・同一周波数で送信する送信部と、
    を有することを特徴とする集積回路。
  18. 　復調用参照信号に、あらかじめ決めた所定の幅の整数倍の信号を加算する摂動ベクトル加算部を有することを特徴とする集積回路。
PCT/JP2012/051923 2011-03-01 2012-01-30 送信装置、受信装置、通信システム、通信方法、および集積回路 WO2012117788A1 (ja)

Priority Applications (1)

Application Number Priority Date Filing Date Title
US14/002,142 US9362995B2 (en) 2011-03-01 2012-01-30 Transmitter apparatus, receiver apparatus, communication system, communication method, and integrated circuit

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2011-043755 2011-03-01
JP2011043755A JP5687524B2 (ja) 2011-03-01 2011-03-01 送信装置、受信装置、通信システム、通信方法、および集積回路

Publications (1)

Publication Number Publication Date
WO2012117788A1 true WO2012117788A1 (ja) 2012-09-07

Family

ID=46757729

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2012/051923 WO2012117788A1 (ja) 2011-03-01 2012-01-30 送信装置、受信装置、通信システム、通信方法、および集積回路

Country Status (3)

Country Link
US (1) US9362995B2 (ja)
JP (1) JP5687524B2 (ja)
WO (1) WO2012117788A1 (ja)

Families Citing this family (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103229532B (zh) * 2011-02-25 2016-04-20 夏普株式会社 无线通信系统、基站装置以及无线通信路径选择方法
JP5804594B2 (ja) * 2011-08-05 2015-11-04 シャープ株式会社 プリコーディング装置、プリコーディング用プログラムおよび集積回路
US20160127058A1 (en) * 2013-06-07 2016-05-05 Hewlett-Packard Development Company, L.P. Target channel identification for a wireless communication
JP6405155B2 (ja) * 2014-08-27 2018-10-17 三星電子株式会社Ｓａｍｓｕｎｇ Ｅｌｅｃｔｒｏｎｉｃｓ Ｃｏ．，Ｌｔｄ． 信号処理装置、信号処理方法、及びプログラム
EP3214782B1 (en) * 2014-11-25 2019-11-13 Huawei Technologies Co., Ltd. Method and apparatus for generating pilot frequency signal
CN111884957B (zh) * 2015-01-28 2023-05-16 索尼公司 无线通信设备和无线通信方法
CN107836089B (zh) * 2015-04-15 2021-05-11 梁平 混合波束成形多天线无线系统
WO2017034105A1 (ko) * 2015-08-24 2017-03-02 엘지전자 주식회사 Dmrs를 전송하는 방법 및 이를 위한 장치
US10778500B2 (en) * 2016-06-27 2020-09-15 Lg Electronics Inc. Method for transmitting and receiving reference signal in wireless communication system and device therefor
JP2020010072A (ja) 2016-11-02 2020-01-16 シャープ株式会社 基地局装置、端末装置および通信方法
WO2018128399A1 (ko) 2017-01-09 2018-07-12 엘지전자 주식회사 무선 통신 시스템에서, 참조 신호를 전송하는 방법 및 이를 위한 장치
WO2018171783A1 (zh) * 2017-03-24 2018-09-27 华为技术有限公司 信号传输方法、装置及系统
CN108632005B (zh) 2017-03-24 2023-12-15 华为技术有限公司 一种参考信号传输方法、装置及系统
WO2019215874A1 (ja) * 2018-05-10 2019-11-14 株式会社Ｎｔｔドコモ 受信装置及び送信装置
CN110784292B (zh) * 2018-07-31 2021-10-01 华为技术有限公司 参考信号的处理方法和装置
US11956006B2 (en) * 2022-07-01 2024-04-09 Qualcomm Incorporated Precoding perturbation to allow large array digital post distortion

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009182894A (ja) * 2008-01-31 2009-08-13 Toshiba Corp 無線送信方法及び装置
WO2010053059A1 (ja) * 2008-11-05 2010-05-14 株式会社 東芝 無線通信装置及び方法
WO2010150313A1 (ja) * 2009-06-25 2010-12-29 株式会社 東芝 通信装置
WO2011108319A1 (ja) * 2010-03-01 2011-09-09 シャープ株式会社 通信システム、送信装置、受信装置

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8139685B2 (en) * 2005-05-10 2012-03-20 Qualcomm Incorporated Systems, methods, and apparatus for frequency control
JP4695003B2 (ja) * 2005-09-05 2011-06-08 Ｏｋｉセミコンダクタ株式会社 等化器及び等化方法
US9094167B2 (en) * 2009-02-02 2015-07-28 Samsung Electronics Co., Ltd. System and method for multi-user and multi-cell MIMO transmissions
EP2453593B1 (en) * 2009-06-19 2018-08-15 BlackBerry Limited Downlink reference signal for type II relay
JP2011254171A (ja) * 2010-05-31 2011-12-15 Sharp Corp 通信システム、送信装置、送信制御方法、送信制御プログラム、及びプロセッサ

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2009182894A (ja) * 2008-01-31 2009-08-13 Toshiba Corp 無線送信方法及び装置
WO2010053059A1 (ja) * 2008-11-05 2010-05-14 株式会社 東芝 無線通信装置及び方法
WO2010150313A1 (ja) * 2009-06-25 2010-12-29 株式会社 東芝 通信装置
WO2011108319A1 (ja) * 2010-03-01 2011-09-09 シャープ株式会社 通信システム、送信装置、受信装置

Also Published As

Publication number Publication date
US20130336282A1 (en) 2013-12-19
JP5687524B2 (ja) 2015-03-18
JP2012182627A (ja) 2012-09-20
US9362995B2 (en) 2016-06-07

Similar Documents

Publication Publication Date Title
JP5687524B2 (ja) 送信装置、受信装置、通信システム、通信方法、および集積回路
JP5686427B2 (ja) 送信装置、受信装置、無線通信システム、送信制御方法、受信制御方法、及び、プロセッサ
JP4702883B2 (ja) 送信装置、受信装置、ｍｉｍｏ−ｏｆｄｍ通信システム及びｍｉｍｏ−ｏｆｄｍ通信システムにおけるｉｑインバランス補償方法
RU2742820C1 (ru) Передающее устройство, приемное устройство, способ передачи и способ приема
JP5221285B2 (ja) 無線通信装置及び方法
JP5689353B2 (ja) フィルタ算出装置、送信装置、受信装置、プロセッサおよびフィルタ算出方法
JP5908307B2 (ja) プリコーディング装置、無線送信装置、無線受信装置、無線通信システムおよび集積回路
WO2011040258A1 (ja) 無線通信システム、通信装置、無線通信方法および端末装置
WO2012060177A1 (ja) 基地局装置、移動局装置及びそれらを用いた無線通信システム
JP2012070366A (ja) 送信装置、受信装置、中継装置、通信システム、送信方法、受信方法、中継方法、通信方法、コンピュータプログラムおよび半導体チップ
US20130003641A1 (en) Communication system, transmitter and receiver
CN103997473B (zh) 一种信号干扰的滤波方法及相关装置
JP2010502140A (ja) 等化構造及び等化方法
US20130343320A1 (en) Terminal device, base station device, and wireless communication system
JP2013030940A (ja) 送信装置、受信装置、および通信システム
JP6159839B1 (ja) 基地局装置
CN102377719A (zh) 数据符号正交处理方法及装置
JP5557543B2 (ja) 無線通信システム、送信装置、受信装置、受信制御方法、及び、受信制御プログラム
WO2013018466A1 (ja) 受信装置、プログラムおよび集積回路
JP6484205B2 (ja) 基地局装置、摂動ベクトル探索方法及びプログラム
JP2012105146A (ja) 端末装置、基地局装置、通信システム、通信方法、及びプロセッサ
JP2010246095A (ja) 送信ウェイト決定方法
JP5676221B2 (ja) 送信装置、通信システム、通信方法、プロセッサ
RU2777352C2 (ru) Передающее устройство, приемное устройство, способ передачи и способ приема

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: 12752992

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 14002142

Country of ref document: US

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 12752992

Country of ref document: EP

Kind code of ref document: A1