Summary of the invention
Technology of the present invention is dealt with problems and is: a kind of subspace ZF sign indicating number householder method that can suppress digital narrow band interference in the cdma system is provided, can reduces existing direct ZF sign indicating number householder method computation complexity, keep system's high output Signal to Interference plus Noise Ratio performance simultaneously.
Technical solution of the present invention is: the subspace of high-performance low computational complexity compels zero, and (ZF, ZeroForcing) the sign indicating number householder method comprises the following steps:
(1) receives wireless communication signals, comprise CDMA signal, white noise and digital narrow band interference.
(2) said wireless communication signals is down-converted to intermediate-freuqncy signal.
(3) the said intermediate-freuqncy signal of digitlization is to provide digital signal.
(4) with said digital demodulation signal, so that baseband signal to be provided.
(5) with said baseband signal envelope r (t) through cutting general matched filtering so that sampled signal r to be provided (m), wherein m is sampled signal label and m=0,1,2 ...
Sampled signal r (m) can be modeled as:
r(m)=y(m)+i(m)+ε(m)
=y
0(m)+z(m)+i(m)+ε(m)
Y in the formula (m) representes the CDMA sampled signal, comprises the signal component y of desired user 0
0(m) and multiple access disturb the signal component z (m) of (MAI, Multiple Access Interference); I (m) representes digital narrow band interference sampled signal; ε (m) expression power spectral density does
The white Gaussian noise sampled signal,
And system synchronization.
The CDMA sampled signal can be modeled as:
Wherein the sampled signal model of desired user 0 and MAI is respectively
K in the formula
zExpression MAI number of users (K
z>=1), A
kRepresent that k user receives signal amplitude, A
0For desired user receives signal amplitude, b
k(n) k subscriber signal stream of expression (1 or-1), b
0(n) be expectation subscriber signal stream, s
kRepresent k user's direct sequence spread spectrum codes sampled value, s
0Be desired user direct sequence spread spectrum codes sampled value, the time domain frequency expansion sequence does before the sampling
{ s in the formula
K, j: j=1 ..., N} is k user's a spreading code (1 or-1), and k=0,1 ... K
Z, ψ
c() is duration T
cReturn-change waveform, N=T
b/ T
cBe spreading gain, T
bBe the signal period, T
cBe the spreading code cycle, and spreading code is independent of signal.
Numeral narrow band interference sampled signal can be modeled as:
K in the formula
iRepresent digital narrow band interference number, A
IkRepresent that k digital narrow band interference receives signal amplitude, b
Ik(n) k digital narrow-band interference signal stream of expression (1 or-1), the duration is T
Ik, the promptly digital narrow band interference cycle.Here T
Ik>>T
cAnd N
Vk=T
b/ T
IkBe integer, υ () representation unit height rectangular pulse.
(6) with said sampled signal r (m) through the windowing memory, send signal for n, handling [nT at interval
b, (n+1) T
b] in to sampled signal r (m) windowing obtain (N * 1) dimension windowing vector r (n)=[r (nN+N-1), r (nN+N-2) ..., r (nN)]
T, wherein n is for sending signal label and n=0,1,2 ..., T
bBe the signal period, N is a spreading gain.
Windowing vector r (n) can be modeled as:
r(n)=y(n)+i(n)+ε(n)
=y
0(n)+z(n)+i(n)+ε(n)
Y in the formula (n) comprises the signal component y of desired user 0 for CDMA adds window signal
0(n) and the signal component z of MAI (n), i (n) adds window signal for digital narrow band interference, and ε (n) adds window signal for white Gaussian noise.
CDMA windowing vector can be modeled as:
The windowing signal model that wherein comprises desired user 0 and MAI
y
0(n)=A
0b
0(n)s
0
In the formula
Numeral narrow band interference windowing vector can be modeled as:
In the formula
And wherein each symbol repeats N
Ik=T
Ik/ T
cIndividual element, N
Ik=T
Ik/ T
cBe digital narrow band interference gain.
(7) said windowing vector r (n) is carried out the auxiliary filtering of subspace ZF sign indicating number of high-performance low computational complexity, comprise the steps:
Step 1: estimate the statistics autocorrelation matrix
R
rr(n)=E{r(n)r
T(n)|b}
=R
yy(n)+R
ii(n)+R
εε(n)
=R
yy0(n)+R
zz(n)+R
ii(n)+R
εε(n)
R in the formula
Yy(n)=E{y (n) y
T(n) } be CDMA windowing signal autocorrelation matrix, comprise the autocorrelation matrix of desired
user 0
Autocorrelation matrix R with MAI
Zz(n)=E{z (n) z
T(n) }, R
Ii(n)=E{i (n) i
T(n) } be digital narrow band interference windowing signal autocorrelation matrix, R
ε ε(n)=E{ ε (n) ε
T(n) } be white Gaussian noise windowing autocorrelation matrix.
FI estimates the CDMA autocorrelation matrix
The autocorrelation matrix that comprises desired user 0 and MAI
Estimative figure narrow band interference autocorrelation matrix
Estimate the white Gaussian noise autocorrelation matrix
I in the formula
NExpression N dimension diagonal matrix.
Step 2: autocorrelation matrix characteristic value decomposition
Λ in the formula
s=diag{ λ
1, λ
2..., λ
KComprise K=1+K
z+ K
iIndividual greater than characteristic value
Diagonal matrix, U
s={ u
1, u
2..., u
K] be the matrix that a corresponding K orthogonal vectors constitute, satisfy
U
ε=[u
K+1, u
K+2..., u
N] be N-K characteristic value
The matrix that corresponding orthogonal vectors constitute,
Be the virtual CDMA subspace that CDMA user and digital narrow band interference are opened, its quadrature component is U
εThe noise subspace opened of row.
Step 3: estimator space ZF filter vector
Setting the auxiliary filter vector of subspace ZF sign indicating number is w=[w
1, w
2..., w
N]
T(N * 1 dimension) is with w=U
sX substitution ZF sign indicating number aided algorithm cost function, x is an intermediate variable
Obtain
Make that the following formula gradient is zero, obtain
Limited Lagrangian factor generation time following formula with following formula substitution constraints
obtains can get
According to
U wherein
sWith U
εThe character of quadrature can obtain
The x that above two formulas are obtained
ZFSubstitution w=U
sX can estimator space ZF sign indicating number aided algorithm filter vector expression formula
Step 4: directly detected symbol judgement
According to directly detecting rule
Carry out symbol judgement after the detection and promptly obtain useful information.
The present invention's beneficial effect compared with prior art is:
The present invention has realized the high-performance low computational complexity filtering of cdma system numeral narrow band interference.The subspace ZF sign indicating number householder method matrix inversion part that provides is improved to K dimension diagonal matrix (K<N) by N * N dimension non-singular matrix of existing direct ZF sign indicating number householder method; Greatly reduce the system-computed complexity, kept the height output Signal to Interference plus Noise Ratio performance close simultaneously with direct ZF sign indicating number householder method.
Embodiment
In the detailed description of the present invention, with reference to appended drawing, the specific exemplary embodiment of these accompanying drawing explainations invention can be implemented in these exemplary embodiments below.These embodiment describe with sufficient details, allowing those skilled in the art's embodiment of the present invention, but can utilize other embodiment, and can make logic, machinery, electrical equipment with other change, and do not depart from standard of the present invention.Therefore, following detailed should not regarded as restrictive, and scope of the present invention is limited by appended claims only.
The embodiment of subspace ZF sign indicating number householder method comprises the following steps:
(1) receives wireless communication signals, comprise CDMA signal, white noise and digital narrow band interference.
(2) said wireless communication signals is down-converted to intermediate-freuqncy signal.
(3) the said intermediate-freuqncy signal of digitlization is to provide digital signal.
(4) with said digital demodulation signal, so that baseband signal to be provided.
(5) with said baseband signal envelope r (t) through cutting general matched filtering so that sampled signal r to be provided (m), wherein m is sampled signal label and m=0,1,2 ...
Sampled signal r (m) can be modeled as:
r(m)=y(m)+i(m)+ε(m)
=y
0(m)+z(m)+i(m)+ε(m)
Y in the formula (m) representes the CDMA sampled signal, comprises the signal component y of desired user 0
0(m) and multiple access disturb the signal component z (m) of (MAI, Multiple Access Interference); I (m) representes digital narrow band interference sampled signal; ε (m) expression power spectral density does
The white Gaussian noise sampled signal,
And system synchronization.
The CDMA sampled signal can be modeled as:
Wherein the sampled signal model of desired user 0 and MAI is respectively
K in the formula
zExpression MAI number of users (K
z>=1), A
kRepresent that k user receives signal amplitude, A
0For desired user receives signal amplitude, b
k(n) k subscriber signal stream of expression (1 or-1), b
0(n) be expectation subscriber signal stream, s
kRepresent k user's direct sequence spread spectrum codes sampled value, s
0Be desired user direct sequence spread spectrum codes sampled value, the time domain frequency expansion sequence does before the sampling
{ s in the formula
K, j: j=1 ..., N} is k user's a spreading code (1 or-1), and k=0,1 ... K
Z, ψ
c() is duration T
cThe normalization waveform, N=T
b/ T
cBe spreading gain, T
bBe the signal period, T
cBe the spreading code cycle, and spreading code is independent of signal.
Numeral narrow band interference sampled signal can be modeled as:
K in the formula
iRepresent digital narrow band interference number, A
IkRepresent that k digital narrow band interference receives signal amplitude, b
Ik(n) k digital narrow-band interference signal stream of expression (1 or-1), the duration is T
Ik, the promptly digital narrow band interference cycle.Here T
Ik>>T
cAnd N
Vk=T
b/ T
IkBe integer, υ () representation unit height rectangular pulse.
(6) with said sampled signal r (m) through the windowing memory, send signal for n, handling [nT at interval
b, (n+1) T
b] in to sampled signal r (m) windowing obtain (N * 1) dimension windowing vector r (n)=[r (nN+N-1), r (nN+N-2) ..., r (nN)]
T, wherein n is for sending signal label and n=0,1,2 ..., T
bBe the signal period, N is a spreading gain.
Windowing vector r (n) can be modeled as:
r(n)=y(n)+i(n)+ε(n)
=y
0(n)+z(n)+i(n)+ε(n)
Y in the formula (n) comprises the signal component y of desired user 0 for CDMA adds window signal
0(n) and the signal component z of MAI (n), i (n) adds window signal for digital narrow band interference, and ε (n) adds window signal for white Gaussian noise.
CDMA windowing vector can be modeled as:
The windowing signal model that wherein comprises desired user 0 and MAI
y
0(n)=A
0b
0(n)s
0
In the formula
Numeral narrow band interference windowing vector can be modeled as:
In the formula
And wherein each symbol repeats N
Ik=T
Ik/ T
cIndividual element, N
Ik=T
Ik/ T
cBe digital narrow band interference gain.
(7) said windowing vector r (n) is carried out the auxiliary filtering of subspace ZF sign indicating number of high-performance low computational complexity, comprising:
At first: estimate the statistics autocorrelation matrix
R
rr(n)=E{r(n)r
T(n)|b}
=R
yy(n)+R
ii(n)+R
εε(n)
=R
yy0(n)+R
zz(n)+R
ii(n)+R
εε(n)
R in the formula
Yy(n)=E{y (n) y
T(n) } be CDMA windowing signal autocorrelation matrix, comprise the autocorrelation matrix of desired
user 0
Autocorrelation matrix R with MAI
Zz(n)=E{z (n) z
T(n) }, R
Ii(n)=E{i (n) i
T(n) } be digital narrow band interference windowing signal autocorrelation matrix, R
ε ε(n)=E{ ε (n) ε
T(n) } be white Gaussian noise windowing autocorrelation matrix.
Estimate the CDMA autocorrelation matrix
The autocorrelation matrix that comprises desired user 0 and MAI
Estimative figure narrow band interference autocorrelation matrix
Estimate the white Gaussian noise autocorrelation matrix
I in the formula
NExpression N dimension diagonal matrix.
Secondly: the autocorrelation matrix characteristic value decomposition
Λ in the formula
s=diag{ λ
1, λ
2..., λ
KComprise K=1+K
z+ K
iIndividual greater than characteristic value
Diagonal matrix, U
s=[u
1, u
2..., u
K] be the matrix that a corresponding K orthogonal vectors constitute, satisfy
U
ε=[u
K+1, u
K+2..., u
N] be N-K characteristic value
The matrix that corresponding orthogonal vectors constitute,
Be the virtual CDMA subspace that CDMA user and digital narrow band interference are opened, its quadrature component is U
εThe noise subspace opened of row.
Once more: estimator space ZF filter vector
Setting the auxiliary filter vector of subspace ZF sign indicating number is w=[w
1, w
2..., w
N]
T(N * 1 dimension) is with w=U
sX substitution ZF sign indicating number aided algorithm cost function, x is an intermediate variable
Obtain
Make that the following formula gradient is zero, obtain
Limited Lagrangian factor generation time following formula with following formula substitution constraints
obtains can get
According to
U wherein
sWith U
εThe character of quadrature can obtain
The x that above two formulas are obtained
ZFSubstitution w=U
sX can estimator space ZF sign indicating number aided algorithm filter vector expression formula
At last: directly detected symbol judgement
According to directly detecting rule
Carry out symbol judgement after the detection and promptly obtain useful information.
Fig. 1 describes the cdma system reception block diagram that existing direct ZF sign indicating number householder method suppresses digital narrow band interference, is made up of down conversion module 101, the auxiliary filtration module 102 of direct ZF sign indicating number and the auxiliary filter vector computing module 103 of direct ZF sign indicating number.
As shown in the figure, the wireless communication signals 105 that antenna 104 receives comprises useful signal CDMA signal, digital narrow band interference and white Gaussian noise.Antenna 104 is coupled to down conversion module 101.In down conversion module 101, at first handle wireless communication signals 105 through band pass filter 106, this filter is optimally selected the frequency wanted, the frequency that for example joins with CDMA signal correction.Band pass filter 106 is coupled to amplifier 107, and this amplifier amplifies the signal from band pass filter 106.Blender 108 mixes the output of amplifier 107 with oscillator signal from local oscillator 109.Like this, the output signal of blender 108 down-conversion amplifiers 107 is to provide intermediate-freuqncy signal 110.After initial down-conversion, intermediate-freuqncy signal 110 is transformed into digital signal 112 through modulus a/d transducer 111.Because the QPSK signal can be regarded as the stack of two quadrature 2PSK signals, so the coherent carrier that intersects with two-way goes demodulation.Wherein one road signal gets into blender 115 and mixes with signal from digital controlled oscillator 113, another road signal entering blender 116 and with mix from the signal of digital controlled oscillator 113 after pi/2 phase shift 114.The signal of blender 115 and blender 116 outputs is connected low pass filter 117 and low pass filter 118 respectively.Thereafter; The output signal of low pass filter 117 and the output signal of low pass filter 118 are connected sampling decision device 119 and sampling decision device 120 respectively; And the output signal of will sample decision device 119 and sampling decision device 120 becomes baseband signal 122 outputs after parallel/serial device 121 conversion.
Like this; Band pass filter 106, amplifier 107, blender 108, local oscillator 109 have been accomplished the process that is down-converted to intermediate-freuqncy signal; Modulus a/d transducer 111, digital controlled oscillator 113, pi/2 phase shift 114, blender 115; 116, low pass filter 117,118, sampling decision device 119,120, parallel/serial device 121 have been accomplished the process that is down-converted to baseband signal by intermediate-freuqncy signal.Above process has been accomplished down conversion module 101.
Existing direct ZF sign indicating number householder method is carried out filtering through the
auxiliary filtration module 102 of direct ZF sign indicating number.Baseband signal 122 envelope r (t) are cut general matched
filtering sampling 123, sampled signal r (m) 124 is provided.With sampled signal r (m) 124 through
windowing memory 125, to sampled signal r (m) 124 windowings obtain (N * 1) dimension windowing vector r (n)=[r (nN+N-1), r (nN+N-2) ..., r (nN)]
T126.The filter vector element w that the auxiliary filter
vector computing module 103 of direct ZF sign indicating number is obtained
1131,
w
2132...
w
N13N sends into linear combination estimator 1:w with element r (nN+N-1) 141, r (nN) 142...r (nN+N-2) 14N of windowing vector r (n) 126 respectively
1R (nN+N-1) 151, linear combination estimator 2:w
2R (nN+N-2) 152... linear combination estimator N:w
NR (nN) 15N obtains the linear combination signal through
adder 161
At last, bit estimated signal
Send into
symbol judgement device 162, obtain
useful signal 163.
Like this, and general filtering sampling 123; Sampling windowing storage 125; Linear combination estimates 1 151, linear combination estimates that 2 152... linear combinations estimate N 15N; Adder 161; Symbol judgement device 162 has constituted the auxiliary filtration module 102 of direct ZF sign indicating number jointly.
The direct required filter vector element w of the auxiliary filtration module of ZF
sign indicating number 102
1131,
w
2132...
w
N13N is provided by the auxiliary filter
vector computing module 103 of direct ZF sign indicating number.According to the channel estimating principle, windowing vector r (n) 126 is sent into estimation statistics autocorrelation matrix computing module 127R
Rr(n)=E{r (n) r
T(n) | b}=R
Yy(n)+R
Ii(n)+R
ε ε(n)=R
Yy0(n)+R
Zz(n)+R
Ii(n)+R
ε ε(n).Send estimating the relevant parameter that statistics autocorrelation
matrix computing module 127 obtains into the direct ZF filter
vector computing module 128 of estimation
To estimate the filter vector element w that direct ZF filter
vector computing module 128 obtains at last
1131,
w
2132...
w
N13N sends into the
auxiliary filtration module 102 of direct ZF sign indicating number.
Like this, estimate statistics autocorrelation matrix computing module 127, add up direct ZF filter vector computing module 128 and constituted the auxiliary filter vector computing module 103 of direct ZF sign indicating number.
The cdma system that the subspace ZF sign indicating number householder method that Fig. 2 describes the high-performance low computational complexity suppresses digital narrow band interference receives block diagram, is made up of down conversion module 201, the auxiliary filtration module 202 of subspace ZF sign indicating number and the auxiliary filter vector computing module 203 of subspace ZF sign indicating number.
As shown in the figure, the wireless communication signals 205 that antenna 204 receives comprises useful signal CDMA signal, digital narrow band interference and white Gaussian noise.Antenna 204 is coupled to down conversion module 201.In down conversion module 201, at first handle wireless communication signals 205 through band pass filter 206, this filter is optimally selected the frequency wanted, the frequency that for example joins with CDMA signal correction.Band pass filter 206 is coupled to amplifier 207, and this amplifier amplifies the signal from band pass filter 206.Blender 208 mixes the output of amplifier 207 with oscillator signal from local oscillator 209.Like this, the output signal of blender 208 down-conversion amplifiers 207 is to provide intermediate-freuqncy signal 210.After initial down-conversion, intermediate-freuqncy signal 210 is transformed into digital signal 212 through modulus a/d transducer 211.Because the QPSK signal can be regarded as the stack of two quadrature 2PSK signals, so the coherent carrier that intersects with two-way goes demodulation.Wherein one road signal gets into blender 215 and mixes with signal from digital controlled oscillator 213, another road signal entering blender 216 and with mix from the signal of digital controlled oscillator 213 after pi/2 phase shift 214.The signal of blender 215 and blender 216 outputs is connected low pass filter 217 and low pass filter 218 respectively.Thereafter; The output signal of low pass filter 217 and the output signal of low pass filter 218 are connected sampling decision device 219 and sampling decision device 220 respectively; And the output signal of will sample decision device 219 and sampling decision device 220 becomes baseband signal 222 outputs after parallel/serial device 221 conversion.
Like this; Band pass filter 206, amplifier 207, blender 208, local oscillator 209 have been accomplished the process that is down-converted to intermediate-freuqncy signal; Modulus a/d transducer 211, digital controlled oscillator 213, pi/2 phase shift 214, blender 215; 216, low pass filter 217,218, sampling decision device 219,220, parallel/serial device 221 have been accomplished the process that is down-converted to baseband signal by intermediate-freuqncy signal.Above process has been accomplished down conversion module 201.
Subspace ZF sign indicating number householder method is carried out filtering through the
auxiliary filtration module 202 of subspace ZF sign indicating number.Baseband signal 222 envelope r (t) are cut general matched
filtering sampling 223, sampled signal r (m) 224 is provided.With sampled signal r (m) 224 through
windowing memory 225, to sampled signal r (m) 224 windowings obtain (N * 1) dimension windowing vector r (n)=[r (nN+N-1), r (nN+N-2) ..., r (nN)]
T226.The filter vector element w that ZF sign indicating number auxiliary filter
vector computing module 203 in subspace is obtained
1231,
w
2232...w
N23N sends into linear combination estimator 1:w with element r (nN+N-1) 241, r (nN) 242...r (nN+N-2) 24N of windowing vector r (n) 226 respectively
1R (nN+N-1) 251, linear combination estimator 2:w
2R (nN+N-2) 252... linear combination estimator N:w
NR (nN) 25N obtains the linear combination signal through
adder 261
At last, with the bit estimated signal
Send into
symbol judgement device 262, obtain
useful signal 263.
Like this, and general filtering sampling 223; Sampling windowing storage 225; Linear combination estimates 1 251, linear combination estimates that 2 252... linear combinations estimate N 25N; Adder 261; Symbol judgement device 262 has constituted the auxiliary filtration module 202 of subspace ZF sign indicating number jointly.
The filter vector element w that the auxiliary filtration module of subspace ZF
sign indicating number 202 is required
1231,
w
2232...w
N23N is provided by the auxiliary filter
vector computing module 203 of subspace ZF sign indicating number.According to the channel estimating principle, windowing vector r (n) 226 is sent into estimation statistics autocorrelation matrix computing module 227R
Rr(n)=E{r (n) r
T(n) | b}=R
Yy(n)+R
Ii(n)+R
Zz(n)=R
Yy0(n)+R
Zz(n)+R
Ii(n)+R
ε ε(n).Send estimating the autocorrelation matrix that statistics autocorrelation
matrix computing module 227 obtains into autocorrelation matrix characteristic
value decomposition module 228
The relevant parameter that again autocorrelation matrix characteristic
value decomposition module 228 is obtained is sent into estimator space ZF filter
vector computing module 229
The filter vector element w that at last estimator space ZF filter
vector computing module 229 is obtained
1231,
w
2232...w
N23N sends into the
auxiliary filtration module 202 of subspace ZF sign indicating number.
Like this, estimate statistics autocorrelation matrix computing module 227, autocorrelation matrix characteristic value decomposition module 228, statistics subspace ZF filter vector computing module 229 has constituted the auxiliary filter vector computing module 203 of subspace ZF sign indicating number.
Fig. 3 describes the subspace ZF sign indicating number householder method of high-performance low computational complexity and the performance simulation correlation curve that existing direct ZF sign indicating number householder method suppresses digital narrow band interference, wherein comprises the emulation correlation curve to single-tone disturbs and multitone disturbs respectively.
The cdma system that simulated conditions is set spreading gain N=63 comprises 4 users; Wherein
user 0 is a desired user; Signal power, i.e. unit signal energy
other MAI and desired user constant power; The CDMA spreading code is all chosen the Gold sequence of
coefficient correlation 1/N.Channel circumstance white Gaussian noise power spectral density
is that the relative ambient noise power of signal power is 20dB (after the despreading).Interference power is with " J " expression, and unit is dB.Set and only have a digital narrow band interference, N
v=4 and-10dB≤J≤30dB.Output Signal to Interference plus Noise Ratio (SINR, Signal to Interference and Noise Rate) has reflected the ratio of system's output mixed signal desired signal and interference and noise power, the dB of unit.
As shown in Figure 3, the output SINR of SINR 301 expression systems, J 302 expression interference powers.When disturbing to digital narrow band interference, the performance 303 of subspace ZF sign indicating number householder method is close with the performance 304 of direct ZF sign indicating number householder method.This has also proved for digital narrow band interference; The subspace ZF sign indicating number householder method that the present invention proposes is with directly ZF sign indicating number householder method performance is close; But its matrix inversion part is improved to K dimension diagonal matrix (K<N), greatly reduce the system-computed complexity by the N * N dimension non-singular matrix of existing direct ZF sign indicating number householder method.Explain that ZF sign indicating number householder method more existing direct ZF sign indicating number householder method in subspace has significant advantage.
Fig. 4 is the flow chart of method.In step 401, wireless communication signals 205 is down-converted to intermediate-freuqncy signal 210.In step 402, at first use modulus a/d transducer 211 that intermediate-freuqncy signal 210 is digitized as digital signal, carry out QPSK demodulation output baseband signal 222 thereafter.In step 403, baseband signal 222 is obtained sampled signal r (m) 224 through cutting general matched filtering sampler 223.In step 404, sampled signal r (m) 224 is obtained windowing vector r (n) 226 through windowing memory 225.In step 405, through estimating statistics autocorrelation matrix computing module 227 counting statistics autocorrelation matrix R
Rr(n)=E{r (n) r
T(n) | b}=R
Yy(n)+R
Ii(n)+R
ε ε(n)=R
Yy0(n)+R
Zz(n)+R
Ii(n)+R
ε ε(n).In step 406, decompose through 228 pairs of autocorrelation matrix characteristic values of autocorrelation matrix characteristic value decomposition module
In step 407, through estimator space ZF filter vector computing module 229 calculating subspace ZF filter vectors
In step 408, the filter vector element w that ZF sign indicating number auxiliary filter vector computing module 203 in subspace is obtained
1231, w
2232...w
N23N sends into linear combination estimator 1:w with element r (nN+N-1) 241, r (nN) 242...r (nN+N-2) 24N of windowing vector r (n) 226 respectively
1R (nN+N-1) 251, linear combination estimator 2:w
2R (nN+N-2) 252... linear combination estimator N:w
NR (nN) 25N obtains the linear combination signal through adder 261
In step 409, with the bit estimated signal
Send into symbol judgement device 262, obtain useful signal 263.In step 410, symbolic label increases progressively n=n+1.In step 411, judge n whether greater than sending the signal message sum, if, process ends, if not, as from 411 send return arrow indicated.