CN105306115A - MIMO system transmitting terminal pre-encoding method based on mercury water injection - Google Patents

Abstract

The invention discloses an MIMO system transmitting terminal pre-encoding method based on mercury water injection. The method comprises the following steps: S1, using an SVD technology to decompose a channel matrix to obtain N parallel sub-channels; S2, carrying out power allocation for the signal-to-noise ratio gamma i of the ith sub-channel according to a mercury water injection method to obtain a power allocation matrix sigma MWF; S3, constituting a modulation diversity unitary matrix VM according to a modulation diversity planisphere rotation method to reduce overlapped constellation points; and S4, using a product of an original pre-encoding matrix VH in S1, the matrix sigma MWF in S2 and the matrix VM in S3 as a transmitting terminal pre-encoding matrix, namely, G = VH sigma MWF VM, to obtain the pre-encoding method disclosed by the invention, and using the matrix to transmit signals. By adopting the MIMO system transmitting terminal pre-encoding method disclosed by the invention, a pre-encoding computation time cost is reduced, mutual information of the channel is improved, and the requirements of practical application are satisfied.

Description

A kind of mimo system transmitting terminal method for precoding based on mercury water filling

Technical field

The present invention relates to wireless communication technology field, be specifically related to a kind of mimo system transmitting terminal method for precoding based on mercury water filling.

Background technology

Multi-antenna technology refers to and adopts multiple antenna at the transmitting terminal of radio communication or receiving terminal, combines a kind of complex art that advanced signal processing technology realizes simultaneously.If adopt multiple antenna at transmitting terminal and receiving terminal simultaneously, be then called MIMO technique (MIMO, MultipleInputMultipleOutput).MIMO is the important breakthrough of wireless communication field intelligent antenna technology, and it can, under the prerequisite not increasing bandwidth, utilize Space-Time Block Coding and spatial reuse code technology to improve capacity and the availability of frequency spectrum of communication system exponentially.

Mimo system transmitting terminal precoding technique needs transmitting terminal accurately to obtain channel condition information (CSI, ChannelSideInformation).For becoming channel slowly, the transmitting terminal as WLAN (wireless local area network) (WLAN, WirelessLocalAreaNetwork) obtains CSI by feedback or time division duplex (TDD, TimeDivisionDuplex) symmetry principle.In fact, in 3G standard WCDMA, transmitting end channel state information (CSIT, ChanelSideInformationAtTheTransmitter) all can obtain at closed loop diversity or send in adaptive array techniques with receiving terminal channel condition information (CSIR, ChanelSideInformationAtTheReceiver).Based on this basis, transmitting terminal precoding technique obtains develop rapidly.Based on different design standards, precoding technique is also different.Such as, minimum mean square error criterion (MMSE), maximum signal noise ratio principle (SNR), maximum mutual information criterion (MI), or minimum bit-error rate criterion (BER) etc.Precoding technique, according to the pretreatment mode of transmitting terminal, can be divided into linear predictive coding and nonlinear precoding.It is well known that singular value decomposition (SVD, SingularValueDecomposition) precoding technique, it is by resolving into equivalent parallel subchannel by mimo channel, and by the power distributing technique of transmitting terminal to reach the maximum channel capacity of system.

Current SVD precoding technique also exists several method, channel diagonalization precoding (CD respectively, ChannelDiagonalizationPrecoder), classical water filling precoding (CWF, ClassicWater-fillingPrecoder), mercury water filling precoding (MWF, MurcuryWater-fillingPrecoder), maximum diversity precoding (MD, MaxDiversityPrecoder), based on the exhaustive precoding (MDE of modulation diversity, odulationDiversity-ExhaustivePrecoder), also has optimum precoding (OP, OptimalPrecoder) etc.But this several method for precoding has weak point, although the little channel mutual information reached of such as channel diagonalization, classical water filling, mercury water filling and maximum diversity precoding expense computing time is also little, be difficult to meet practical application request; Although the mutual information reached based on exhaustive precoding and the optimum precoding of modulation diversity expense large but computing time is also large, be difficult in actual applications realize.

Summary of the invention

In order to overcome the shortcoming of prior art existence with not enough, the invention provides a kind of MIMO transmitting terminal method for precoding based on mercury water filling.

The present invention adopts following technical scheme:

Based on a mimo system transmitting terminal method for precoding for mercury water filling, comprise the steps:

S1 adopts SVD technology to decompose channel matrix H at transmitting terminal, and obtain N number of parallel subchannel, described H resolves into v _hassociate matrix, the two can be changed mutually, Σ _hrepresent the singular value matrix of H;

S2 is according to the signal to noise ratio γ of mercury water flood for the i-th sub-channels _icarry out power division, obtain power division matrix Σ _mWF, wherein σ _idiagonal matrix Σ _hi-th singular value, P is total transmitting power, σ ²it is noise power;

S3 adopts modulation diversity planisphere spinning solution, constructs a modulation diversity unitary matrice V _m, be specially:

$V_M=\frac{1}{\sqrt{N_t}}\left[\begin{array}{cccc}1& {\mathrm{\α}}_1& \mathrm{...}& {\mathrm{\α}}_1^{N_t-1}\\ 1& {\mathrm{\α}}_2& \mathrm{...}& {\mathrm{\α}}_2^{N_t-1}\\ & & .& \\ & & .& \\ & & .& \\ 1& {\mathrm{\α}}_{N_t}& \mathrm{...}& {\mathrm{\α}}_{N_t}^{N_t-1}\end{array}\right]$

${\mathrm{\α}}_k=\exp \[j\mathrm{\π}\frac{2(k-1)+q_M}{N_t}\],k=1,2,...,N_t$

Wherein, N _tnumber of transmit antennas, N _rreception antenna number, q _mthen determine according to different modulation systems, j represents imaginary unit, i.e. j × j=-1;

S4 is by the original precoding matrices V in S1 _hwith the power division matrix Σ obtained in S2 _mWFand the modulation diversity unitary matrice V obtained in S3 _mproduct as transmitting terminal pre-coding matrix, i.e. G=V _hΣ _mWFv _m, and utilize this matrix emissive signal.

S1 Chinese style namely given channel matrix H application singular value decomposition is decomposed into three matrix multiples on the right of equal sign.

In described S2, mercury water flood specific implementation is according to following formula:

p _i＝0,γ _i≤γ ₀

γ _iMMSE(p _iγ _i)＝γ ₀,γ _i>γ ₀

Wherein p _irepresent the power that i-th parallel subchannel is assigned to, γ ₀be a threshold value, specifically determine according to channel matrix H, N is number of subchannels, and MMSE is Minimum Mean Square Error.

When signal to noise ratio is high, large subchannel can distribute to less power on the contrary.

When two sub-channels situations are identical, the channel that modulation radix M is larger can be assigned to more power.

Described q _mits selection principle be the compound planisphere that emission signal vector is formed constellation point between Euclidean distance as far as possible large.

The present invention is based on the mimo system transmitting terminal method for precoding of mercury water filling, the power distribution method of mercury water filling and modulation diversity planisphere spinning solution combine by the method.First adopt SVD technology, mimo channel is resolved into equivalent parallel subchannel, then according to the channel conditions of each parallel subchannel, carry out power division according to mercury water filling strategy, then avoid constellation point overlapping in conjunction with planisphere spinning solution.Thus design the method for precoding meeting practical application request.

Beneficial effect of the present invention:

The present invention, with maximum mutual information criterion, uses mercury water-filling method to carry out power division, reduces expense computing time; Modulation diversity planisphere spinning solution is used to improve the error performance of receiving terminal.Both combine, and reduce the amount of calculation of method for precoding on the one hand, improve channel mutual information on the other hand, thus meet practical application request.

Accompanying drawing explanation

Fig. 1 is workflow diagram of the present invention;

Fig. 2 is mimo system multi-antenna channel model;

Fig. 3 is mimo system equivalent parallel subchannel model;

Fig. 4 is mercury water injection power allocation strategy of the present invention;

Fig. 5 is the channel mutual information curve that in the present embodiment, 2 × 2MIMO-BPSK modulation adopts different method for precoding;

Fig. 6 is the channel mutual information curve that in the present embodiment, 2 × 2MIMO-QPSK modulation adopts different method for precoding;

Fig. 7 is the channel mutual information curve that in the present embodiment, 4 × 4MIMO-BPSK modulation adopts different method for precoding;

Fig. 8 is the channel mutual information curve that in the present embodiment, 4 × 4MIMO-QPSK modulation adopts different method for precoding.

Embodiment

Below in conjunction with embodiment and accompanying drawing, the present invention is described in further detail, but embodiments of the present invention are not limited thereto.

Embodiment

As shown in Figure 1, a kind of mimo system transmitting terminal method for precoding based on mercury water filling, core is that two key steps obtain matrix Σ respectively _mWFand V _m, final goal obtains pre-coding matrix G.

Fig. 2 is a N _ttransmit antennas and N _rthe mimo system channel model of root reception antenna, as seen from the figure, always co-exists in h ₁₁, h ₁₂..., n altogether _r× N _tthe coupling channel that bar is different.

Consider that the baseband equivalence model of mimo system is: y=HGx+n.

Wherein x is N _tthe transmitting complex vector of × 1, y is N _rthe reception complex vector of × 1, H is N _r× N _tcomplex channel matrix, G is N _t× N _tpre-coding matrix.Carry out SVD to pre-coding matrix G equally to decompose and can obtain: $G=U_G{\mathrm{\Σ}}_G{V}_G^h.$

In method for precoding of the present invention, the left side unitary matrice U of pre-coding matrix G _gget V _h, diagonal matrix Σ _gthen get Σ according to mercury water injection power distribution method _mWF, the right unitary matrice get V _mso method for precoding of the present invention is: G=V _hΣ _mWFv _m.

Fig. 3 is the MIMO equivalent parallel subchannel model that a SVD decomposes, and altogether resolves into N number of subchannel, wherein N=min (N _t, N _r).Can be obtained by above-mentioned formula: distortion can obtain: $\tilde{y}={\mathrm{\Σ}}_H{\mathrm{\Σ}}_{\mathrm{MWF}}\tilde{x}+\tilde{n},$ Wherein $\tilde{y}=U_H^{h}y,\tilde{x}=V_{M}x,\tilde{n}=U_H^{h}n.$ In uncoupled mimo system model, Σ _hand Σ _mWFit is all diagonal matrix. with be i-th singular value of these two matrixes respectively, the former reflects channel status, and the latter represents power division weights.

Fig. 4 is mercury water injection power allocation strategy.In figure, each container is an equivalent parallel subchannel.The process of mercury water filling is, first toward each container influx height be 1/ γ _isolid matrix, next, toward each container perfusion mercury, the height comprising substrate is reached finally toward priming petock in all containers, until water surface elevation rises to 1/ γ ₀.Here suppose that chamber wall is transparent to water, and can not penetrate mercury and solid matrix.So the water-column p of each container _ibe exactly the power division weights that we will obtain, p _i=P _i/ P.

First this example sets up the mimo system of two transmit antennas and two reception antennas, and adopts binary phase shift keying (BPSK) to modulate, and channel matrix is $H=\left[\begin{array}{cc}2& 1\\ 1& 1\end{array}\right].$ Below that two key steps of core obtain power division matrix Σ respectively _mWFwith modulation diversity matrix V _m.

The first step, first adopts SVD technology to decompose channel matrix H, obtains two equivalent parallel subchannels and determine power division matrix Σ thus _mWF.The power division obtaining two sub-channels according to mercury water-filling algorithm is as shown in table 1 below: (supposing that total transmitting power is 2)

Table 1

Signal to noise ratio (dB)	Subchannel 1	Subchannel 2
			-20	2	0
-17.5	2	0
			-15	2	0
-12.5	2	0
			-10	2	0
-7.5	2	0
			-5	2	0
-2.5	1.58	0.42
			0	0.92	1.08
2.5	0.54	1.46
			5	0.33	1.67
7.5	0.21	1.79
			10	0.14	1.86
12.5	0.1	1.9
			15	0.08	1.92
17.5	0.06	1.94
			20	0.06	1.94

Then when signal to noise ratio is-20dB, ${\mathrm{\Σ}}_{MWF}=\left[\begin{array}{cc}\sqrt{2}& 0\\ 0& 0\end{array}\right].$

Second step, determines modulation diversity matrix V _m.Owing to adopting BPSK modulation, q _m=1; If adopt Quadrature Phase Shift Keying (QPSK) modulation, q _m=1/2.So the modulation diversity matrix of this example is: $V_M=\frac{1}{\sqrt{2}}\left[\begin{array}{cc}1& e^{j(\mathrm{\π}/2)}\\ 1& e^{j(3\mathrm{\π}/2)}\end{array}\right],$ So we just obtain pre-coding matrix of the present invention: G=V _hΣ _mWFv _m.

The simulation result of this example uses software Matlab to obtain.The calculating of channel mutual information adopts Lower Bound Formula, thus simplifies calculating.

Simulation result as shown in Figure 5, can be seen, lists the channel mutual information of six kinds of method for precoding under different state of signal-to-noise in figure respectively.These six kinds of methods are a kind of mimo system transmitting terminal method for precoding (ModifiedMWF) based on mercury water filling of the present invention, the exhaustive precoding (ExhaustiveMethod) based on modulation diversity, classical water filling precoding (ClassicWF), mercury water filling precoding (MercuryWF), channel diagonalization precoding (ChannelDiagonalization) and maximum diversity precoding (MaxDiversity) respectively.As seen from the figure, the mutual information that ModifiedMWF with ExhaustiveMethod reaches is substantially identical, be optimum in six kinds of methods in figure, and MaxDiversity is suboptimum when being greater than-7.5dB; When being greater than-5dB, ClassicWF is the poorest, and the performance can seeing MercuryWF is better than ClassicWF; When being greater than 0dB, ChannelDiagonalization is secondary poor.And ModifiedMWF and ExhaustiveMethod to converge to the limit be all the fastest, but the former expense computing time is less.

Adopt the simulation result of QPSK modulation as shown in Figure 6.As seen from the figure, in-2.5dB to 7.5dB interval, the performance of ModifiedMWF with ExhaustiveMethod is no longer consistent, the performance of ModifiedMWF and ExhaustiveMethod is basically identical, the mutual information that the former reaches comparatively the latter almost, but the former also reaches the mutual information of the latter 97.7% when 2.5dB difference is maximum.

The mimo system of four transmit antennas and four reception antennas adopts the simulation result of BPSK modulation as shown in Figure 7.The performance of same ModifiedMWF with ExhaustiveMethod is no longer consistent, at-7.5dB in 2.5dB interval, the mutual information that the former reaches comparatively the latter almost, when-5dB difference is maximum, the former reaches the mutual information of the latter 99.3%.

Adopt the simulation result of QPSK modulation as shown in Figure 8.At-2.5dB in 10dB interval, the mutual information that the former reaches comparatively the latter almost, when 5dB difference is maximum, the former reaches the mutual information of the latter 96.4%.The simulation result of this situation runs with the CPU of IntelCorei5-45903.30GHZ to obtain, and the running time adopting ModifiedMWF is 3.4525min, and adopts be 16.5995h the running time of ExhaustiveMethod, is the former about 288 times.The former expense computing time is much little as can be seen here, and the latter exponentially increases along with its time overhead that increases of order of modulation and number of antennas, be difficult to realize in practical engineering application so that in the practical application of extensive multiple antennas (MassiveMIMO).

The software program that the method described in conjunction with embodiment disclosed herein or step can perform as DSP and FPGA, processor with hardware, or the combination of the two is implemented.Software program can be placed in the storage medium of other form any known in random asccess memory (RAM), internal memory, read-only memory (ROM), electrically programmable ROM, register, hard disk, moveable magnetic disc, CD-ROM or technical field.

Above-described embodiment is the present invention's preferably execution mode; but embodiments of the present invention are not limited by the examples; change, the modification done under other any does not deviate from Spirit Essence of the present invention and principle, substitute, combine, simplify; all should be the substitute mode of equivalence, be included within protection scope of the present invention.

Claims (6)

1., based on a mimo system transmitting terminal method for precoding for mercury water filling, it is characterized in that, comprise the steps:

S1 adopts SVD technology to decompose channel matrix H at transmitting terminal, and obtain N number of parallel subchannel, described H resolves into v _hassociate matrix, the two can be changed mutually, Σ _hrepresent the singular value matrix of H;

S2 is according to the signal to noise ratio γ of mercury water flood for the i-th sub-channels _icarry out power division, obtain power division matrix Σ _mWF, wherein σ _idiagonal matrix Σ _hi-th singular value, P is total transmitting power, σ ²it is noise power;

S3 adopts modulation diversity planisphere spinning solution, constructs a modulation diversity unitary matrice V _m, be specially:

$V_M=\frac{1}{\sqrt{N_t}}\left[\begin{array}{cccc}1& {\mathrm{\α}}_1& \mathrm{...}& {\mathrm{\α}}_1^{N_t-1}\\ 1& {\mathrm{\α}}_2& \mathrm{...}& {\mathrm{\α}}_2^{N_t-1}\\ & & .& \\ & & .& \\ & & .& \\ 1& {\mathrm{\α}}_{N_t}& \mathrm{...}& {\mathrm{\α}}_{N_t}^{N_t-1}\end{array}\right]$

${\mathrm{\α}}_k=\exp \[j\mathrm{\π}\frac{2(k-1)+q_M}{N_t}\],k=1,2,\mathrm{...},N_t$

Wherein, N _tnumber of transmit antennas, N _rreception antenna number, q _mthen determine according to different modulation systems, j represents imaginary unit, i.e. j × j=-1;

S4 is by the original precoding matrices V in S1 _hwith the power division matrix Σ obtained in S2 _mWFand the modulation diversity unitary matrice V obtained in S3 _mproduct as transmitting terminal pre-coding matrix, i.e. G=V _hΣ _mWFv _m, and utilize this matrix emissive signal.

2. method for precoding according to claim 1, is characterized in that, S1 Chinese style namely given channel matrix H application singular value decomposition is decomposed into three matrix multiples on the right of equal sign.

3. method for precoding according to claim 1, is characterized in that, in described S2, mercury water flood specific implementation is according to following formula:

p _i＝0,γ _i≤γ ₀

γ _iMMSE(p _iγ _i)＝γ ₀,γ _i>γ ₀

Wherein p _irepresent the power that i-th parallel subchannel is assigned to, γ ₀be a threshold value, specifically determine according to channel matrix H, N is number of subchannels, and MMSE is Minimum Mean Square Error.

4. method for precoding according to claim 1, is characterized in that, when signal to noise ratio is high, large subchannel can distribute to less power on the contrary.

5. method for precoding according to claim 1, is characterized in that, when the channel conditions of two sub-channels is identical, the subchannel that modulation radix M is larger can be assigned to more power.

6. method for precoding according to claim 1, is characterized in that, described q _mits selection principle be the compound planisphere that emission signal vector is formed constellation point between Euclidean distance as far as possible large.