Summary of the invention
For solving the deficiencies in the prior art, the object of the present invention is to provide a kind of MIMO blind Channel Estimation ambiguity removal method, the method is based on Orthogonal Space-Time Block Code (Orthogonal Space-Time Block Code, OSTBC), can efficient solution the ambiguity problem of blind Channel Estimation in cooperation MIMO communication by no means.
In order to realize above-mentioned target, the present invention adopts following technical scheme:
A kind of MIMO blind Channel Estimation ambiguity removal method based on Orthogonal Space-Time Block Code, is characterized in that, comprises the following steps:
(1), set up, analyze and process MIMO space-time system model:
(1) set up one and there is N
tindividual transmitting antenna and N
rthe MIMO space-time system of individual reception antenna, each frame signal of transmission antennas transmit is made up of training sequence, Cyclic Prefix and data three parts, and aforementioned training sequence is time domain orthogonal training sequence, and channel is smooth slow change rayleigh fading channel, supposes continuous N
fthe channel status of frame signal experience is identical;
(2) analyze time domain orthogonal training sequence signal:
Suppose the time domain orthogonal training sequence signal N of k arbitrary antenna of moment
t× 1 dimension complex vector s
tr(k) represent, receive signal indication and be:
y
tr(k)=Hs
tr(k)+v(k)
In formula, H is N
r× N
tdimension rayleigh fading channel response matrix, y
tr(k) be N
r× 1 dimension received signal vector, v (k) is N
r× 1 dimension noise vector, aforementioned noise obeys that average is 0, variance is
gaussian Profile;
(3) before transmitting, data-signal is divided into groups:
Make s (k)=[s
1(k) s
2(k) ... s
n(k)]
tfor armed k the packet being formed by N symbol, and each symbol independent same distribution;
(4) process k packet s (k):
S (k) is mapped as to N through Space Time Coding
t× L dimension encoder matrix C (k):
In formula, A
nand B
nbe respectively corresponding to n symbol s
n(k) real part (s
rn) and imaginary part (s (k)
in(k) encoder matrix), L is the number of encoder matrix time slot, receives signal indication to be:
Y(k)=HC(k)+V(k)
In formula, Y (k) is N
r× L dimension receives signal matrix, and V (k) is N
r× L ties up noise matrix, and Y (k), V (k) obey that average is 0, variance is
gaussian Profile;
(2), channel estimating and decoding:
(1) export the estimated signal of first data block:
First, combine and utilize channel estimation methods blind and based on time domain orthogonal training sequence, first data block is carried out to blind Channel Estimation; Then, adopt signal and channel associating two-dimensional search to do maximum likelihood decoding to first data block, export the estimated signal of first data block;
(2) estimated signal of output subsequent data blocks:
First, utilize the time domain orthogonal training sequence estimated value that zero forcing equalization obtains to carry out the channel estimating based on time domain orthogonal training sequence to subsequent data blocks, make the estimation channel of subsequent data blocks there is the phase ambiguity identical with first data block; Then, adopt transmitted signal linear search to do maximum likelihood decoding to subsequent data blocks, the estimated signal of output subsequent data blocks.
The aforesaid MIMO blind Channel Estimation ambiguity removal method based on Orthogonal Space-Time Block Code, is characterized in that, in the time that first data block is carried out to blind Channel Estimation, utilizes chu sequence to do time domain orthogonal training sequence, adopts Alamouti delivery plan.
The aforesaid MIMO blind Channel Estimation ambiguity removal method based on Orthogonal Space-Time Block Code, is characterized in that, in the time that first data block is carried out to blind Channel Estimation, utilizes QPSK modulation sequence to do time domain orthogonal training sequence, adopts Alamouti delivery plan.
Usefulness of the present invention is: by adopting signal and channel associating two-dimensional search, utilize the correlation between signal in orthogonal Space Time Coding, on the basis of blind Channel Estimation, eliminate the intrinsic phase ambiguity of channel estimation results, dwindle the scope of channel feasible solution, thereby be that the phase place of determining decoded signal is provided convenience; In addition, this scheme is by utilizing respectively channel estimation methods blind and based on time domain orthogonal training sequence at different pieces of information piece, by phase deviation unification to fixed value of the different pieces of information piece decoded data once sending, solve the problem that causes decoding error in non-cooperative communication due to the ambiguity problem of blind Channel; Method of the present invention can be applied to various multi-antenna signals, the blind identification of collaboration communication signal, blind Detecting.
Embodiment
Below in conjunction with the drawings and specific embodiments, the present invention is done to concrete introduction.
One, set up, analyze and process MIMO space-time system model
1, set up MIMO space-time system model
Set up one and there is N
tindividual transmitting antenna and N
rthe MIMO space-time system of individual reception antenna, each frame signal of transmission antennas transmit is made up of training sequence, Cyclic Prefix and data three parts, and channel is smooth slow change rayleigh fading channel, supposes continuous N
fthe channel status of frame signal experience is identical, and a data block is by N
fframe composition, as shown in Figure 1.
Training sequence in wireless communication system is divided into time domain orthogonal training sequence, frequency domain quadrature training sequence and code territory quadrature training sequence, and in the space-time system model of setting up in the present invention, training sequence adopts time domain orthogonal training sequence.
2, analyze time domain orthogonal training sequence signal
Suppose the time domain orthogonal training sequence signal N of k arbitrary antenna of moment
t× 1 dimension complex vector s
tr(k) represent, receive signal indication and be:
y
tr(k)=Hs
tr(k)+v(k)
In formula, H is N
r× N
tdimension rayleigh fading channel response matrix, y
tr(k) be N
r× 1 dimension received signal vector, v (k) is N
r× 1 dimension noise vector, wherein, noise obeys that average is 0, variance is
gaussian Profile.
3, process MIMO space-time system model
Data-signal, because needs carry out Space Time Coding, therefore first divided into groups before transmitting, makes s (k)=[s
1(k) s
2(k) ... s
n(k)]
tfor armed k the packet being formed by N symbol, and each symbol independent same distribution.
Next process k packet s (k), concrete, s (k) is mapped as to N through Space Time Coding
t× L dimension encoder matrix C (k):
In formula, A
nand B
nbe respectively corresponding to n symbol s
n(k) real part (s
rn) and imaginary part (s (k)
in(k) encoder matrix), L is the number of encoder matrix time slot.Receiving signal can be expressed as:
Y(k)=HC(k)+V(k)
In formula, Y (k) is N
r× L dimension receives signal matrix, and V (k) is N
r× L ties up noise matrix, and Y (k), V (k) obey that average is 0, variance is
gaussian Profile.
In addition,, because channel becomes slowly, the channel of adjacent data block n+1 and data block n is followed following relation:
H
nobey multiple Gaussian Profile with △ H, β is constant coefficient, and the less channel variation of β value is less.
Two, channel estimating and decoding
In Blind channel estimation algorithm in the past, the phase ambiguity that each channel estimating produces is not identical and be random, and this brings very large difficulty to being correctly decoded.
Next,, with reference to Fig. 2, introduce in detail channel estimating involved in the present invention and decoding scheme.
1, first data block is carried out to channel estimating and decoding
Channel estimating part: combine and utilize channel estimation methods blind and based on time domain orthogonal training sequence, first data block is carried out to blind Channel Estimation, by phase ambiguity unification to fixed value of first data block, this fixed value is determined by semanteme judgement.
Decoded portion: due to the intrinsic ambiguity of blind Channel Estimation, channel estimating exists many groups to separate.The present invention combines two-dimensional search in decoded portion by employing signal and channel, utilize the correlation between signal in Space Time Coding, first data block is done to maximum likelihood decoding, on the basis of channel estimating, further dwindle the scope of channel feasible solution, thereby provide convenience for the definite of final phase ambiguity.
Unified for the decoded signal phase place deflection that all data blocks are caused by blind Channel Estimation ambiguity, we only carry out the two-dimensional search of signal and channel in the maximum likelihood decoding of first data block, and detailed process is shown in Fig. 3.
When first data block is carried out to blind Channel Estimation, because this blind Channel Estimation is utilized time domain orthogonal training sequence and completed, and for different time domain orthogonal training sequences, the performance of scheme exists difference.
Processing scheme while having there is in recent years numerous multi-antenna space for different communication environment.Wherein, Space-Time Block Coding (STBC) scheme, by providing diversity gain improve connective stability and improve data transmission rate, is scheme while having representational sky.But its performance depends on the accuracy of channel estimating to a great extent.
Taking classical Alamouti scheme as example, its encoder matrix can be expressed as
Therefore, below we will be taking two transmitting antennas as example, the channel estimating and the decoding problem that are done the Alamouti delivery plan of time domain orthogonal training sequence by chu sequence and QPSK modulation sequence are discussed respectively.
(1), utilize chu sequence as time domain orthogonal training sequence
In the time not considering noise, there is the phase ambiguity of 180 ° or 0 ° in whole mimo system channel matrix.Concrete derivation is as follows:
Make h
1, h
2represent respectively the true subchannel of two transmit antennas, s
1, s
2represent actual transmission signal, receive signal and be:
Order
represent channel estimation value,
represent final decoded result, if the subchannel estimated value of two antennas all exists 180 ° of phase place deflections, while decoding with such channel, we can obtain:
Order
α
1=α
2=-1, result is identical with reception signal, and error vector is 0, therefore in the case, 180 ° of phase place deflections can occur decoded signal.
From deriving, for Alamouti delivery plan, decoded portion has certain selection effect to phase ambiguity, makes final decoded result only may have two kinds, is respectively (s
1s
2) and (s
1-s
2), wherein, (s
1s
2) be correct decoding.
(2), utilize QPSK modulation sequence to do time domain orthogonal training sequence
In the time not considering noise, there is the phase ambiguity of 0 °, 180 °, 90 ° or-90 ° in whole mimo system channel matrix.Concrete derivation is as follows:
Make h
1, h
2represent respectively the true subchannel of two transmit antennas, s
1, s
2represent actual transmission signal, receive signal and be:
Order
represent channel estimation value,
represent final decoded result, while decoding with such channel, α
1with α
2respectively desirable 1 ,-1, i or-i, remove and be correctly decoded (α
1=α
2=1) situation, we need to consider to also have three kinds of situations:
If 1. the subchannel estimated value of two antennas all exists 180 ° of phase place deflections, we can obtain:
Order
α
1=α
2=-1, result is identical with reception signal, and error vector is 0, therefore in the case, 180 ° of phase place deflections can occur decoded signal.
If 2. first subchannel estimated value exists 90 ° of phase place deflections, there are-90 ° of phase place deflections in the second sub-channels estimated value, and we can obtain:
Order
α
1=-i, α
2=i, result is identical with reception signal, and error vector is 0, and therefore in the case, there are-90 ° of phase place deflections in the decoded signal of antenna 1, and there are 90 ° of phase place deflections in the decoded signal of antenna 2.
If 3. first subchannel estimated value exists-90 ° of phase place deflections, there are 90 ° of phase place deflections in the second sub-channels estimated value, and we can obtain:
Order
α
1=i, α
2=-i, result is identical with reception signal, and error vector is 0, and therefore in the case, there are 90 ° of phase place deflections in the decoded signal of antenna 1, and there are-90 ° of phase place deflections in the decoded signal of antenna 2.
To sum up aforementioned, in the time using QPSK modulation sequence, as time domain orthogonal training sequence, Alamouti delivery plan is carried out to channel estimating and decoding, decoded result has following four kinds of possibilities: (s
1s
2), (s
1-s
2), (is
1is
2), (is
1-is
2).
2, subsequent data blocks is carried out to channel estimating and decoding
Channel estimating part: utilize the time domain orthogonal training sequence estimated value that zero forcing equalization obtains to carry out the channel estimating based on time domain orthogonal training sequence to subsequent data blocks, make the estimation channel of subsequent data blocks there is the phase ambiguity identical with first data block, thereby by phase ambiguity unification to fixed value of the different pieces of information piece once sending.
Decoded portion: adopt transmitted signal linear search to do maximum likelihood decoding to subsequent data blocks, the estimated signal of output subsequent data blocks.
Visible, further determining (after first data block being carried out to channel estimating and decode) after channel, what to the decoding of subsequent data blocks, we adopted is general maximum likelihood algorithm, adopt transmitted signal linear search to do maximum likelihood decoding, thereby ensured that the phase place deflection that different pieces of information piece decoded signal carries is identical, this detailed process is shown in Fig. 4.
Three, scheme emulation and performance evaluation
In order to obtain estimated performance, we use respectively the time domain orthogonal training sequence being made up of chu sequence and the time domain orthogonal training sequence being made up of QPSK modulation sequence, estimate channel by direct SVD decomposition method, and phase ambiguity is eliminated with real channel, the statistics error rate, result is specific as follows:
1, chu sequence is as time domain orthogonal training sequence, employing Alamouti delivery plan
Simulated conditions: set up one and there are two the eight MIMO space-time systems of receiving, adopt Alamouti delivery plan, send and receive two data blocks, each data block 4 frames.Send time domain orthogonal training sequence in data and use chu sequence, be made up of 32 symbols, data division is made up of 64 symbols.Carry out altogether Monte Carlo simulation 1000 times.
The planisphere that receives signal is shown in Fig. 5.
The planisphere of equalizing signal is shown in Fig. 6.
Decoding performance curve chart is shown in Fig. 7.
As can be seen here: use chu sequence as time domain orthogonal training sequence, under Alamouti space time coding scheme, the algorithm of carrying equilibrium be preferably subject to channel effect to be aliasing in reception signal together, thereby correctly recover system transmitted signal, and the error rate reduces with the increase of signal to noise ratio, reaches 10 at 12dB
-2.
2, QPSK modulation sequence is as time domain orthogonal training sequence, employing Alamouti delivery plan
Simulated conditions: set up one and there are two the eight MIMO space-time systems of receiving, adopt Alamouti delivery plan, send and receive two data blocks, each data block 4 frames.Send time domain orthogonal training sequence in data and use QPSK modulation sequence, be made up of 32 symbols, data division is made up of 64 symbols.Carry out altogether Monte Carlo simulation 1000 times, the statistics error rate.
The planisphere that receives signal is shown in Fig. 8.
The planisphere of equalizing signal is shown in Fig. 9.
Decoding performance curve chart is shown in Figure 10.
As can be seen here: use QPSK modulation sequence as time domain orthogonal training sequence, under Alamouti space time coding scheme, the algorithm of carrying equilibrium be preferably subject to channel effect to be aliasing in reception signal together, thereby correctly recover system transmitted signal, and the error rate reduces with the increase of signal to noise ratio, reaches 10 at 15dB
-2.
While utilizing QPSK modulation sequence as time domain orthogonal training sequence, the decoding performance of the decoding performance of Alamouti delivery plan and the complete known system of channel information is compared, see Figure 11.
As can be seen here: under the scene of considering at this patent, receiving terminal cannot be known time domain orthogonal training sequence particular content and channel condition information, therefore under the bit error rate performance of blind Channel Estimation and decoding and cooperative communication scene, the bit error rate performance that known real channel is decoded has certain gap, but also can realize to a certain extent system proper communication.
To sum up aforementioned, method of the present invention is first by adopting signal and channel associating two-dimensional search, utilize the correlation between signal in Space Time Coding, on the basis of blind Channel Estimation, eliminate the intrinsic phase ambiguity of channel estimation results, dwindle the scope of channel feasible solution, thereby be that the phase place of determining decoded signal is provided convenience.
In addition, the solution of the present invention is by utilize respectively channel estimation methods blind and based on time domain orthogonal training sequence at different pieces of information piece, by phase deviation unification to fixed value of the different pieces of information piece decoded data once sending.
Generally speaking, method of the present invention has solved the impact that the intrinsic phase ambiguity of blind Channel Estimation causes system decodes.
Method of the present invention can be applied to various multi-antenna signals, the blind identification of collaboration communication signal, blind Detecting.
It should be noted that, above-described embodiment does not limit the present invention in any form, and all employings are equal to replaces or technical scheme that the mode of equivalent transformation obtains, all drops in protection scope of the present invention.