CN100358261C - An array antenna MC-CDMA system up-link receiving method - Google Patents

An array antenna MC-CDMA system up-link receiving method Download PDF

Info

Publication number
CN100358261C
CN100358261C CNB2005100864156A CN200510086415A CN100358261C CN 100358261 C CN100358261 C CN 100358261C CN B2005100864156 A CNB2005100864156 A CN B2005100864156A CN 200510086415 A CN200510086415 A CN 200510086415A CN 100358261 C CN100358261 C CN 100358261C
Authority
CN
China
Prior art keywords
tau
signal
spatial domain
integral
prime
Prior art date
Legal status (The legal status 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 status listed.)
Expired - Fee Related
Application number
CNB2005100864156A
Other languages
Chinese (zh)
Other versions
CN1741414A (en
Inventor
杨维
刘俊英
陈俊仕
程时昕
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Jiaotong University
Original Assignee
Beijing Jiaotong University
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 Beijing Jiaotong University filed Critical Beijing Jiaotong University
Priority to CNB2005100864156A priority Critical patent/CN100358261C/en
Publication of CN1741414A publication Critical patent/CN1741414A/en
Application granted granted Critical
Publication of CN100358261C publication Critical patent/CN100358261C/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Noise Elimination (AREA)
  • Radio Transmission System (AREA)

Abstract

The present invention relates to an uplink receiving method of an array antenna MC-CDMA system, which comprises the steps that signals received by each array element of the array antenna are carried out with the separation of sub-carrier signals; the sub-carrier signals of each separated array element are treated by dispreading and matched filtering so as to realize the dispreading of subscriber signals and obtain the matched filtering output of subscriber sub-carrier signals; the signals belonging to the same subscriber and the same sub-carrier in the signals outputted by a matched filter on different array elements are carried out with spatial domain combination so as to obtain a spatial domain diversity gain; each obtained sub-carrier signals which are carried out with the spatial domain combination and belong to the same subscriber and the same sub-carrier are carried out with frequency domain combination so as to obtain a frequency domain diversity gain; a final decision variable of the subscriber signals after the spatial domain combination and the frequency domain combination is decided to obtain the decision result of the subscriber signals. The present invention has the advantages that processing is simple, the method is easy to realize, and the system has good receptivity because the combination of the spatial domain diversity gain and the frequency domain diversity gain is obtained.

Description

A kind of array antenna MC-CDMA system up-link receiving method
Technical field
The invention belongs to MC-CDMA cell mobile communication systems array antenna technical field.
Background technology
CDMA is a kind of topmost technology in the third generation (3G) mobile communication, and multi-carrier modulation will be the key technology of following wideband wireless transmission system.Merge CDMA technology with multi-transceiver technology, constituting multi-carrier CDMA system is one of important directions of future mobile communications development.The scheme that multi-transceiver technology combines with CDMA technology mainly contains CDMA multiple carrier (MC-CDMA), multi-carrier direct sequence spectrum CDMA (multicarrier DS-CDMA) and three kinds of principal modes of multitone modulation CDMA (MT-CDMA).Wherein, the MC-CDMA scheme also is one of the most competitive scheme of future mobile communication system owing to can adopt frequency diversity and good performance to be considered in three kinds of schemes the scheme of tool prospect.
The same with cdma system, use array antenna can improve capacity, spectrum efficiency, communication quality and the coverage of system significantly at the MC-CDMA system base-station.Existing array antenna MC-CDMA system is handled in the advanced line frequency of receiving terminal territory, and then carries out the spatial domain and handle, and frequency domain and spatial domain are handled and organically do not associated, and handles complicated.
Summary of the invention
The present invention has proposed a kind of array antenna MC-CDMA system up-link spatial domain and frequency domain merging method of reseptance for solving the problems of the technologies described above, compare with array antenna MC-CDMA system up-link receiving method in the past, the advanced line space of method of the present invention territory merges, and then carry out frequency domain and merge, the weight vector of spatial domain merging and the weight vector of frequency domain merging are organically associated, method of reseptance is handled simple, be easy to realize, and owing to obtained spatial domain, the frequency diversity gain of associating, system has good receptivity.
The technical solution adopted in the present invention is that the signal that at first each array element of array antenna is received carries out the separation of sub-carrier signal; Secondly, the sub-carrier signal of isolated each array element of institute is carried out despreading and matched filter processing, realize the despreading of subscriber signal, and the matched filtering that obtains user's sub-carrier signal is exported; Then, carry out the spatial domain merging, obtain the spatial domain diversity gain belonging to the same subcarrier signals of same user in the matched filter output signal on the different array elements; Afterwards, each sub-carrier signal that belongs to same user in the signal after more resultant spatial domain being merged carries out frequency domain and merges, and obtains the frequency diversity gain; At last, the conclusive judgement variable of resultant subscriber signal after spatial domain, frequency domain merge is adjudicated, obtain the court verdict of subscriber signal.
Below method of the present invention is discussed.
1. the separation of sub-carrier signal
Investigate array antenna MC-CDMA system up-link, array antenna is adopted in the base station, and travelling carriage adopts single antenna.Travelling carriage uplink scheme adopts traditional single antenna MC-CDMA transmission plan and QPSK modulation system.
K mobile subscriber arranged in certain cellular cell in the supposing the system.Each subcarrier of MC-CDMA transmission plan is carried identical user profile, and k the signal that the user launched can be expressed as:
s k ( t ) = Σ l = 1 L P k , l ( b k I ( t ) c k , l I ( t ) cos ( ω l t ) + b k Q ( t ) c k , l Q ( t ) sin ( ω l t ) ) [formula 1]
In the formula, P K, lRepresent k user l (signal power of individual subcarrier of 1≤l≤L), ω lBe the angular frequency of l subcarrier, bit period is T s, c k IAnd c k QBe respectively the frequency expansion sequence of I road and Q road correspondence, its expression formula is respectively c k I = [ c k , 1 I c k , 2 l . . . c k , L I ] With c k Q = [ c k , 1 Q c k , 2 Q . . . c k , L Q ] , Chip period T c=T sHomophase and quadrature signal component b k IAnd b k QBe expressed as respectively:
b k I ( t ) = Σ l = - ∞ ∞ d k I rect T ( t - iT s )
[formula 2]
b k Q ( t ) = Σ l = - ∞ ∞ d k Q rect T ( t - iT s )
In the formula, d k IAnd d k QBe corresponding digital signal, value is ± 1, rect T(t) being is T in the cycle sSquare wave.
To the MC-CDMA transmission plan, wireless channel is a frequency non-selective fading channel.The multiple low-pass impulse response of l sub-carrier channels of k user can be expressed as:
h K, l(t)=ρ K, le -j φ k, lδ (t-τ K, l) [formula 3]
Wherein, suppose complex channel coefficient ρ K, le -j φ k, lBe the multiple Gaussian random variable of zero-mean, ρ K, lAnd φ K, lBe respectively corresponding amplitude and phase gain, variance is σ K, l 2τ K, lBe the time delay of l subcarrier of user k,, suppose τ the different user different sub carrier K, lBe independent identically distributed stochastic variable, [0, T s) between evenly distribute.
Suppose that the base station adopts equidistantly linear battle array, like this base station array antenna n (=1 ..., N) signal that receives on the individual array element is:
r n ( t ) = Σ k = 1 K Σ l = 1 L P k , l ρ k , l [ b k I ( t - τ k , l ) c k , l I ( t - τ k , l ) cos ( ω l ( t - τ k , l ) + φ k , l )
[formula 4]
+ b k Q ( t - τ k , l ) c k , l Q ( t - τ k , l ) sin ( ω l ( t - τ k , l ) + φ k , l ) ] a k , l , n + n n ( t )
Wherein, a k , l , n = exp [ - j 2 πd λ k , l ( n - 1 ) sin ( θ k , l ) ] Be the array response of l subcarrier of user k n array element, λ K, lBe corresponding carrier wavelength, d is the distance between the adjacent array element, θ K, lIt is the angle of arrival of l subcarrier of user k.n n(t) be additive white Gaussian noise, bilateral power spectral density is N 0/ 2.
Like this array antenna received to resultant signal can be expressed as:
r ( t ) = [ r 1 ( t ) , r 2 ( t ) , · · · , r N ( t ) ] T
= Σ k = 1 K Σ l = 1 L P k , l ρ k , l [ b k I ( t - τ k , l ) c k , l I ( t - τ k , l ) cos ( ω l ( t - τ k , l ) + φ k , l ) [formula 5]
+ b k Q ( t - τ k , l ) c k , l Q ( t - τ k , l ) sin ( ω l ( t - τ k , l ) + φ k , l ) ] a k , l + n ( t )
Wherein, () TBe the transposition computing, n (t)=[n 1..., n N] TBe corresponding noise vector, a K, lBe the array vector of l subcarrier of k user,
a K, l=[a K, l, 1, a K, l, 2..., a K, l, N] T[formula 6]
Be without loss of generality, the signal of supposing the 1st user is a desired signal.The signal times that each array element is received is with the subcarrier cos (ω corresponding with the transmitting terminal modulation system l, t), sin (ω l, combination t), and corresponding time delay taken into account, just can realize the separation of each sub-carrier signal.To isolated the 1st the user l ' of QPSK modulation system base station array antenna (1≤l '≤L) individual subcarrier signals is:
x 1,l′(t)=[x 1,l′,1(t),x 1,l′,2(t),…,x 1,l′,N(t)] T=r(t)(cos(ω l′(t-τ 1,l′))+sin(ω l′(t-τ 1,l′)))
=[r 1(t)(cos(ω l′(t-τ 1,l′))+sin(ω l′(t-τ 1,l′))),r 2(t)(cos(ω l′(t-τ 1,l′))+sin(ωo l′(t-τ 1,l′))),
…,r N(t)(cos(ω l′(t-τ 1,l′))+sin(ω l′(t-τ 1,l′)))] T
[formula 7]
2. despreading and matched filter processing
Isolated sub-carrier signal is carried out despreading and matched filter processing, realize the despreading of subscriber signal and obtain the matched filtering output of user's sub-carrier signal.Base station array antenna is output as the matched filtering of i bit of the 1st the individual carrier signal of user l ':
y 1 , l ′ ( i ) = [ y 1 , l ′ , 1 , · · · , y 1 , l ′ , N ] T
= 1 T s ∫ iT s + τ 1 , l ′ ( i + 1 ) T s + τ 1 , l ′ x 1 , l ′ ( t ) c 1 , l ′ I ( t ) dt + j 1 T s ∫ iT s + τ 1 , l ′ ( i + 1 ) T s + τ 1 , l ′ x 1 , l ′ c 1 , l ′ Q ( t ) dt [formula 8]
= P 1 , l ′ ρ 1 , l ′ ( b 1 I ( i ) + jb 1 Q ( i ) ) a 1 , l ′ ( i ) + m 1 , l ′ ( i ) + n 1 , l ′ ( i )
Wherein, m 1, l '(i), n 1, l '(i) represent total interference and noise signal respectively.
The noise of system and interference characteristic have significant effects to the empty performance that merges method of reseptance frequently that is proposed.Below the kind of system noise and interference is analyzed earlier.
1) noise
Noise item can provide as follows:
n 1 , l ′ ( i ) = n 1 , l ′ I ( i ) + jn 1 , l ′ Q ( i ) [formula 9]
Homophase has identical variance with quadrature component
σ N I 2 = σ N Q 2 = σ N 2 / 2 = N 0 4 T s [formula 10]
2) disturb
To the MC-CDMA transmission plan, wireless channel is a frequency non-selective fading channel, still keeps orthogonality between the subcarrier, and therefore, total interference is the interference of other user's same carrier wave, comprises homophase and quadrature component, can be expressed as:
m 1 , l ′ ( i ) = m 1 , l ′ I ( i ) + jm 1 , l ′ Q ( i ) [formula 11]
m 1, l ' IAnd m 1, l ' QBe respectively:
m 1 , l ′ I ( i ) = Σ k = 2 K a k , l ′ P k , l ′ ρ k , l ′ cos ( ω l ′ ( τ k , l ′ - τ 1 , l ′ ) + ( φ k , l ′ - φ 1 , l ′ ) )
[formula 12]
c 1 , l ′ I [ ( b k I ( i - 1 ) c k , l ′ I + jb k Q ( i - 1 ) c k , l ′ Q ) R g ( τ k , l ′ ′ ) + ( b k I ( i ) c k , l ′ I + jb k Q ( i ) c k , l ′ Q ) R g ( T c - τ k , l ′ ′ ) ]
m 1 , l ′ Q ( i ) = Σ k = 2 K a k , l ′ P k , l ′ ρ k , l ′ cos ( ω l ′ ( τ k , l ′ - τ 1 , l ′ ) + ( φ k , l ′ - φ 1 , l ′ ) )
[formula 13]
c 1 , l ′ Q [ ( b k I ( i - 1 ) c k , l ′ I + jb k Q ( i - 1 ) c k , l ′ Q ) R g ( τ k , l ′ ′ ) + ( b k I ( i ) c k , l ′ I + jb k Q ( i ) c k , l ′ I ) R g ( T c - τ k , l ′ ′ ) ]
In the formula, R gBe the partial auto correlation of chip waveform (τ), be defined as follows:
R g ( τ ) = 1 T c ∫ 0 τ g ( t + T c - τ ) g ( t ) dt , 0 ≤ τ ≤ T c [formula 14]
Wherein, τ K, l '' be the relative time delay of l subcarrier of k user with respect to the 1st the individual subcarrier of user l '.
3. the spatial domain of signal merges
Below the correlation properties of disturbing are further analyzed, draw optimum and suboptimum merging weight that the received signal spatial domain merges.Disturb correlation properties to comprise auto-correlation and their cross correlation.
1) auto-correlation
Disturb vectorial m 1, l ' I(i) and m 1, l ' Q(i) element m 1, l ', n I(i) and m 1, l ', n Q(i) can be regarded as a series of independent Gaussian stochastic variables.Therefore, their auto-correlation function can approximate evaluation be:
E { m 1 , l ′ , n I ( i ) m 1 , l ′ m I * ( i + h ) } = E { m 1 , l ′ , n Q ( i ) m 1 , l ′ , n Q * ( i + h ) } ≅ σ m , 1 , l ′ 2 δ ( h ) , 1 ≤ n ≤ N [formula 15]
Wherein, () *Be conjugate operation, δ (h) is the delta function.
The variance that can get homophase and quadrature component is:
σ m , 1 , l ′ 2 = Σ k = 2 K P k , l ′ σ k , l ′ 2 G ( R g 2 ( τ k , l ′ ′ ) + R g 2 ( T s - τ k , l ′ ′ ) ) [formula 16]
Can suppose τ K, l '' be independent identically distributed stochastic variable, and [0, T c] go up and obey evenly distribution, g (t) is to be T in the cycle cSquare wave, can get:
σ m , 1 , l ′ 2 = Σ k = 2 K 2 3 σ k , l ′ 2 P k , l ′ [formula 17]
Therefore, in-phase component adds the population variance of quadrature component and is:
σ I 2 = 2 σ m , 1 , l ′ 2 [formula 18]
2) cross-correlation
m 1, l ' n I(i) and m 1, l ', n ' I(i) cross-correlation between is (perhaps at m 1, l ', n Q(i) and m 1, l ', n ' Q(i) be that the spatial domain of received signal on different array element n and n ' is relevant the cross-correlation between).It can be derived out by the cross-correlation function of linear array, and the cross-correlation function of linear array is expressed as follows:
E { a k , l ′ ( t ) a k , l ′ H ( t + τ ) } = J 0 ( 2 π f d τ ) · R k , l ′ [formula 19]
Wherein, () HThe expression conjugate transpose, f dIt is maximum Doppler frequency-shift.R K, l 'Be k the individual subcarrier array vector of user l ' a of expression K, l 'The N of correlation * N Hemitian Toeplitz matrix, its real part and imaginary part can be expressed as respectively:
Re { R k , l ′ ( n ′ , n ) } = J 0 ( r k , l ′ ( n ′ , n ) ) + 2 Σ q = 1 ∞ J 2 q ( r k , l ′ ( n ′ , n ) ) · cos ( 2 q θ k , l ′ ) [formula 20]
Im { R k , l ′ ( n ′ , n ) } = 2 Σ q = 1 ∞ J 2 q + 1 ( r k , l ′ ( n ′ , n ) ) · sin ( ( 2 q + 1 ) θ k , l ′ ) [formula 21]
In the formula, J qBe first kind q rank Bessel functions,
r K, l '(n ', n)=2 π d|n '-n|/λ K, l '[formula 22]
3) total interference plus noise correlation matrix
Definition interference correlation matrix M (h) is:
M ( h ) = E { m 1 , l ′ I ( i ) m 1 , l ′ IH ( i + h ) + m 1 , l ′ Q ( k ) m 1 , l ′ QH ( i + h ) } = Mδ ( h ) [formula 23]
Matrix M can be derived by formula 17 and 18:
M = 2 Σ k = 2 K 2 3 σ k , l ′ 2 P k , l ′ a k , l ′ a k , l ′ H [formula 24]
The distribution function f (θ) that supposes the angle of arrival is equally distributed within [0,2 π], and formula 24 can further be expressed as:
M = 2 σ m , 1 , l ′ 2 ∫ θ a ( θ ) a H ( θ ) · f ( θ ) dθ = σ I 2 ∫ θ R ( θ ) f ( θ ) dθ [formula 25]
R (θ) is the matrix in formula 20 and the formula 21, so the element of the row of the n ' row n in the matrix M is:
m n ′ , n = σ I 2 J 0 ( r ( n ′ , n ) ) [formula 26]
Therefore, total interference plus noise correlation matrix is:
R T = E { ( m 1 , l ′ + n 1 , l ′ ) ( m 1 , l ′ + n 1 , l ′ ) H } = M + σ N 2 I [formula 27]
In the formula, I is a unit matrix.
R TDetermining the optimum weight vector that merges in spatial domain of following formula, promptly the optimum weight vector that merges in the spatial domain of the 1st the individual carrier wave of user l ' is:
w 1 , l ′ , opt s = α R T - 1 a 1 , l ′ H [formula 28]
Wherein, a 1, l ' HBe the Hermite conversion of the individual subcarrier array response of desired user l ', α is a constant, can be taken as 1.
Work as R TWhen approaching unit matrix, interference plus noise and can be considered to empty white noise frequently, promptly have
R T = M + σ N 2 I = ( σ I 2 + σ N 2 ) I = σ 2 I [formula 29]
Wherein, σ I 2Be defined in formula 18, σ N 2It is the noise variance that is shown in formula 10.Like this, can obtain suboptimum spatial domain merging weight vector is:
w 1 , l ′ s ( i ) = [ w 1 , l ′ , 1 s ( i ) , w 1 , l ′ , 2 s ( i ) , . . . , w 1 , l ′ , N s ( i ) ] = α a 1 , l ′ H ( i ) [formula 30]
The present invention only uses the suboptimum spatial domain of having simplified to merge weight, compare with the optimum merging, systematic function descends little, but because the suboptimum spatial domain merges the knowledge that weight only need be known user's subcarrier array vector, do not merge the contrary of the necessary interference plus noise correlation matrix of weight and do not need to calculate optimum spatial domain, greatly reduce the complexity of weight calculation.The signal of the 1st the individual carrier wave of user l ' merges i the bit in back through the spatial domain and is output as:
z 1 , l ′ ( i ) = w 1 , l ′ s ( i ) y 1 , l ′ ( i ) = P 1 , l ′ ρ 1 , l ′ ( b 1 I ( i ) + jb 1 Q ( i ) ) | a 1 , l ′ ( i ) | + η 1 , l ′ ( i ) [formula 31]
Wherein, | a 1 , l ′ | = a 1 , l ′ H · a 1 , l ′ a 1 , l ′ H · a 1 , l ′ , η 1 , l ′ ( i ) For the spatial domain merge that the back is disturbed and noise item with.
Corresponding homophase and quadrature judgment variables are:
z 1 , l ′ I ( i ) = Re { z 1 , l ′ ( i ) } = P 1 , l ′ ρ 1 . l ′ b 1 I ( i ) | a 1 , l ′ | + η 1 , l ′ I ( i )
[formula 32]
z 1 , l ′ Q ( i ) = Im { z 1 , l ′ ( i ) } = P 1 , l ′ ρ 1 . l ′ b 1 Q ( i ) | a 1 , l ′ | + η 1 , l ′ Q ( i )
4. signal frequency domain combination
Each subcarrier to the same user of MC-CDMA scheme carries identical user profile.The output signal that the present invention merges each subcarrier spatial domain is carried out frequency domain again and is merged, and adopts high specific to merge (MRC), and the high specific on the individual subcarrier of l ' merges weight and is:
w 1 , l ′ f = | a 1 , l ′ | [formula 33]
| a 1, l '| will be estimated when merging, so the weight that frequency domain merges does not need to reappraise in the signal spatial domain.
Like this, through the conclusive judgement variable after the empty merging frequently be:
Z 1 ( i ) = Σ l ′ = 1 L w 1 , l ′ f z 1 , l ′ ( i ) [formula 34]
Corresponding homophase and quadrature judgment variables are:
Z 1 I ( i ) = Re { Z 1 ( i ) }
[formula 35]
Z 1 Q ( i ) = Im { Z 1 ( i ) }
5. the judgement of signal
The present invention utilizes simple polarity decision method to adjudicate to resulting conclusive judgement variable, obtains the court verdict of user's sub-carrier signal.To the QPSK modulation system, be:
b 1 I ( i ) = 1 , Z 1 I ( i ) > 0 - 1 , Z 1 I ( i ) < 0 [formula 36]
b 1 Q ( i ) = 1 , Z 1 Q ( i ) > 0 - 1 , Z 1 Q ( i ) < 0 [formula 37]
Beneficial effect of the present invention:
At first, because isolated sub-carrier signal is carried out despreading and matched filter processing earlier, suppressed the multiple access interference to a great extent, to make the signal processing of receiver based on disturbing repressed signal to carry out, stability and performance that this can improve receiver greatly make the method for reseptance that is proposed have practical outstanding advantage.
Secondly, the spatial domain of signal merges has adopted suboptimum to merge weight, does not need to calculate optimum spatial domain and merges the contrary of the necessary interference plus noise correlation matrix of weight, greatly reduce the complexity of weight calculation, and systematic function descends very little; The frequency domain that carries out afterwards merges the high specific merging method that adopted, and merging weight does not need to reappraise, and it is simple that method of reseptance is handled, and is easy to realize.
At last, the empty method of reseptance that frequently merges that is proposed has obtained spatial domain, frequency diversity gain by spatial domain, frequency domain merging, makes system have good receptivity.
Description of drawings
Fig. 1 is that array antenna MC-CDMA system up-link (arbitrary user k) receives structured flowchart;
Fig. 2 is the simulation result of error rate of system (BER) when adopting optimum and suboptimum merging weight vector;
Fig. 3 is the not error performance of simultaneity factor of array number.
Embodiment
Below in conjunction with accompanying drawing method of the present invention is described in detail.
Method of the present invention is applicable to the array antennas mobile communication system of any employing MC-CDMA transmission plan.
Receive structured flowchart with reference to the array antenna MC-CDMA system up-link of Fig. 1, a kind of concrete steps of array antenna MC-CDMA system up-link receiving method comprise, to arbitrary user (k):
Step 1, the signal r that array antenna the 1st array element is received 1(t) send into L carrier signal separation module (1-1-1) respectively, (1-2-1) ..., (1-L-1) in, obtain output signal x K, 1,1(t)=r 1(t) (cos (ω 1(t-τ K, l))+sin (ω 1(t-τ K, l))), x K, 2,1(t)=r 1(t) (cos (ω 2(t-τ K, 2))+sin (ω 2(t-τ K, 2))) ..., x K, L, 1(t)=r 1(t) (cos (ω L(t-τ K, L))+sin (ω L(t-τ K, L))); The signal r that the 2nd array element is received simultaneously 2(t) send into L carrier signal separation module (1-1-2) respectively, (1-2-2) ..., (1-L-2) in, obtain output signal x K, 1,2(t)=r 2(t) (cos (ω 1(t-τ K, l))+sin (ω 1(t-τ K, l))), x K, 2,2(t)=r 2(t) (cos (ω 2(t-τ K, 2))+sin (ω 2(t-τ K, 2))) ..., x K, L, 2(t)=r 2(t) (cos (ω L(t-τ K, L))+sin (ω L(t-τ K, L))); ...; Simultaneously with N the signal r that array element received N(t) send into L carrier signal separation module (1-1-N) respectively, (1-2-N) ..., (1-L-N) in, obtain output signal x K, 1, N(t)=r N(t) (cos (ω 1(t-τ K, 1))+sin (ω 1(t-τ K, 1))), x K, 2, N(t)=r N(t) (cos (ω 2(t-τ K, 2))+sin (ω 2(t-τ K, 2))) ..., x K, L, N(t)=r N(t) (cos (ω L(t-τ K, L))+sin (ω L(t-τ K, L)));
Step 2 is with signal x K, 1,1(t), x K, 2,1(t) ..., x K, L, 1(t) send into L despreading and matched filtering module (2-1-1) respectively, (2-2-1) ..., (2-L-1), carry out despreading and matched filter processing, obtain the i bit signal:
y k , 1,1 ( i ) = 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1,1 ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , 1 c k , 1 Q ( t ) dt ,
y k , 2,1 ( i ) = 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2,1 ( t ) c k , 2 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , 1 c k , 2 Q ( t ) dt , . . . ,
y k , L , 1 ( i ) = 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , 1 ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , 1 c k , 1 Q ( t ) dt ;
With signal x K, 1,2(t), x K, 2,2(t) ..., x K, L, 2(t) send into L despreading and matched filtering module (2-1-2) respectively, (2-2-2) ..., (2-L-2) in, obtain output signal:
y k , 1,2 ( i ) = 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , 2 ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , 2 c k , 1 Q ( t ) dt ,
y k , 2,2 ( i ) = 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2,2 ( t ) c k , 2 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , 2 c k , 2 Q ( t ) dt , . . . ,
y k , L , 2 ( i ) = 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , 2 ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , 2 c k , 1 Q ( t ) dt ;
With signal x K, 1, N(t), x K, 2, N(t) ..., x K, L, N(t) send into despreading and matched filtering module (2-1-N) respectively, (2-2-N) ..., (2-L-N) in, obtain output signal:
y k , 1 , N ( i ) = 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , N ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , N c k , 1 Q ( t ) dt ,
y k , 2 , N ( i ) = 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , N ( t ) c k , 2 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , N c k , 2 Q ( t ) dt , . . . ,
y k , L , N ( i ) = 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , N ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , N c k , 1 Q ( t ) dt ;
Step 3 is with N matched filtering output signal y K, 1,1(i), y K, 1,2(i) ..., y K, 1, N(i) send into the spatial domain and merge module (3-1), obtain vector y K, 1(i)=[y K, 1,1(i), y K, 1,2(i) ..., y K, 1, N(i)] T, merge weight vector by the spatial domain w k , 1 s ( i ) = [ w k , 1 , 1 s ( i ) , w k , 1 , 2 s ( i ) , . . . , w k , 1 , N s ( i ) ] , The spatial domain of finishing the 1st sub-carrier signal of user k merges, and obtains the signal after the spatial domain merges z k , 1 ( i ) = w k , 1 s ( i ) &CenterDot; y k , 1 ( i ) ; Simultaneously with N matched filtering output signal y K, 2,1(i), y K, 2,2(i) ..., y K, 2, N(i) send into the spatial domain and merge module (3-2), obtain vector y K, 2(i)=[y K, 2,1(i), y K, 2,2(i) ..., y K, 2, N(i)] T, merge weight vector by the spatial domain w k , 2 s ( i ) = [ w k , 2,1 s ( i ) , w k , 2,2 s ( i ) , . . . , w k , 2 , N s ( i ) ] , The spatial domain of finishing the 2nd sub-carrier signal of user k merges, and obtains the signal after the spatial domain merges z k , 2 ( i ) = w k , 2 s ( i ) &CenterDot; y k , 2 ( i ) , . . . , Simultaneously with N matched filtering output signal y K, L, 1(i), y K, L, 2(i) ..., y K, L, N(i) send into the spatial domain and merge module (3-L), obtain vector y K, L(i)=[y K, L, 1(i), y K, L, 2(i) ..., y K, L, N(i)] T, merge weight vector by the spatial domain w k , L s ( i ) = [ w k , L , 1 s ( i ) , w k , L , 2 s ( i ) , . . . , w k , L , N s ( i ) ] , The spatial domain of finishing L sub-carrier signal of user k merges, and obtains the signal after the spatial domain merges z k , L ( i ) = w k , L s ( i ) &CenterDot; y k , L ( i ) ;
Step 4 is with L sub-carrier signal z after the merging of spatial domain K, 1(i), z K, 2(i) ..., z K, L(i) send into carrier wave and merge module (4), obtain vector z k(i)=[z K, 1(i), z K, 2(i) ..., z K, L(i)] T, merge weight vector by frequency domain w k f ( i ) = [ w k , 1 f ( i ) , w k , 2 f ( i ) , &CenterDot; &CenterDot; &CenterDot; , w k , L f ( i ) ] Finish the merging of multi-carrier signal, obtain the final judgment variables of user k Z k ( i ) = w k f ( i ) &CenterDot; z k ( i ) ;
Step 5 is with final judgment variables Z k(i) send into signal decision module (5) and carry out polarity decision, obtain the court verdict of arbitrary user (k) signal.
A kind of array antenna MC-CDMA system up-link receiving method, the spatial domain that it is characterized in that step 3 merge and to be based on suboptimum and to merge that weight vector carries out, and the weight vector of l carrier signal spatial domain merging of k user is w k , l s ( i ) = &alpha; a k , l H ( i ) .
A kind of array antenna MC-CDMA system up-link receiving method is characterized in that the carrier wave merging in the step 4 is undertaken by the high specific merging criterion, and k user's multi-carrier signal merges weight vector and is w k f ( i ) = [ | a k , 1 ( i ) | , | a k , 2 ( i ) | , . . . , | a k , L ( i ) | ] .
Fig. 2 and Fig. 3 have provided the performance simulation result of a kind of array antenna MC-CDMA system up-link receiving method that adopts the present invention's proposition.Supposing the system has reached synchronously in emulation is crossed, and the array element distance of the even linear battle array of employing is λ/2, if there is not channel fading, all transmit a signal to and have identical power when reaching the base station so, and Normalized Signal/Noise Ratio is E b/ N 0=10dB, the spread processing gain is 128 with sub-carrier number, and channel is the Rayleigh slow fading channel, and the user is distributed in the sub-district randomly, and its ripple reaches the angle and evenly distributes between [0,2 π].
Fig. 2 has provided the simulation result when number of users BER performance when the inventive method does not adopt optimum and suboptimum merging weight vector respectively simultaneously.As can be seen from Figure 2 adopt suboptimum to merge weight vector and compare with adopting the optimum weight vector that merges, the BER decreased performance of system is very little, but can reduce to calculate the amount of calculation of weight greatly, and it is suitable therefore adopting suboptimum merging method.
Fig. 3 has provided the BER performance that adopts the inventive method system when the array number of array antenna is respectively 1,3 and 5.As can be seen from Figure 3, along with the increase of array number,, improved the performance of system significantly owing to obtained bigger space-frequency diversity gain.

Claims (3)

1. the big line MC-CDMA of an array system up-link method of reseptance is characterized in that comprising following receiving step, to arbitrary user k:
Step 1, the signal r that array antenna the 1st array element is received 1(t) send into the pairing L of the 1st array element carrier signal separation module (1-1-1) respectively, (1-2-1) ..., (1-L-1) in, obtain output signal
x k,1,1(t)=r 1(t)(cos(ω 1(t-τ k,1))+sin(ω 1(t-τ k,1))),
x k,2,1(t)=r 1(t)(cos(ω 2(t-τ k,2))+sin(ω 2(t-τ k,2))),...,
x K, L, 1(t)=r 1(t) (cos (ω L(t-τ K, L))+sin (ω L(t-τ K, L))); The signal r that the 2nd array element is received simultaneously 2(t) send into the pairing L of the 2nd array element carrier signal separation module (1-1-2) respectively, (1-2-2) ..., (1-L-2) in, obtain output signal x K, 1,2(t)=r 2(t) (cos (ω 1(t-τ K, 1))+sin (ω 1(t-τ K, 1))), x K, 2,2(t)=r 2(t) (cos (ω 2(t-τ K, 2))+sin (ω 2(t-τ K, 2))) ..., x K, L, 2(t)=r 2(t) (cos (ω L(t-τ K, L))+sin (ω L(t-τ K, L))); ...; Simultaneously with N the signal r that array element received N(t) send into the pairing L of N array element carrier signal separation module (1-1-N) respectively, (1-2-N) ..., (1-L-N) in, obtain output signal x K, 1, N(t)=r N(t) (cos (ω 1(t-τ K, 1))+sin (ω 1(t-τ K, 1))),
x k,2,N(t)=r N(t)(cos(ω 2(t-τ k,2))+sin(ω 2(t-τ k,2))),...,
x k,L,N(t)=r N(t)(cos(ω L(t-τ k,L))+sin(ω L(t-τ k,L)));
Step 2 is with signal x K, 1,1(t), x K, 2,1(t) ..., x K, L, 1(t) send into the pairing L of the 1st an array element despreading and matched filtering module (2-1-1) respectively, (2-2-1) ..., (2-L-1), carry out despreading and matched filter processing, obtain the i bit signal:
y k , 1,1 ( i ) = 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1,1 ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1,1 ( t ) c k , 1 Q ( t ) dt ,
y k , 2,1 ( i ) = 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2,1 ( t ) c k , 2 I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , , 21 ( t ) c k , 2 Q ( t ) dt , . . . ,
y k , L , l ( i ) = 1 T s &Integral; i T s + &tau; k , L ( i + 1 ) T s + &tau; k , 1 x k , L , l ( t ) c k , 1 l ( t ) dt + j 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , L , 1 ( t ) c k , 1 Q ( t ) dt ;
With signal x K, 1,2(t), x K, 2,2(t) ..., x K, L, 2(t) send into the pairing L of the 2nd an array element despreading and matched filtering module (2-1-2) respectively, (2-2-2) ..., (2-L-2) in, obtain output signal:
y k , 1 , 2 ( i ) = 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , 2 ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , 2 ( t ) c k , 1 Q ( t ) dt ,
y k , 2 , 2 ( i ) = 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , 2 ( t ) c k , 2 I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , 2 ( t ) c k , 2 Q ( t ) dt , . . . ,
y k , L , 2 ( i ) = 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , 2 ( t ) c k , L I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , 2 ( t ) c k , L Q ( t ) dt ,
With signal x K, 1, N(t), x K, 2, N(t) ..., x K, L, N(t) send into the pairing L of a N array element despreading and matched filtering module (2-1-N) respectively, (2-2-N) ..., (2-L-N) in, obtain output signal:
y k , 1 , N ( i ) = 1 T s &Integral; iT s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , N ( t ) c k , 1 I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , 1 ( i + 1 ) T s + &tau; k , 1 x k , 1 , N ( t ) c k , 1 Q ( t ) dt ,
y k , 2 , N ( i ) = 1 T s &Integral; iT s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , N ( t ) c k , 2 I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , 2 ( i + 1 ) T s + &tau; k , 2 x k , 2 , N ( t ) c k , 2 Q ( t ) dt , . . . ,
y k , L , N ( i ) = 1 T s &Integral; iT s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , N ( t ) c k , L I ( t ) dt + j 1 T s &Integral; i T s + &tau; k , L ( i + 1 ) T s + &tau; k , L x k , L , N ( t ) c k , L Q ( t ) dt ,
Step 3 is with N matched filtering output signal y K, 1,1(i), y K, 1,2(i) ..., y K, 1, N(i) send into the 1st the pairing spatial domain of subcarrier and merge module (3-1), obtain vector y K, 1(i)=[y K, 1,1(i), y K, 1,2(i) ..., y K, 1, N(i)] T, merge weight vector by the spatial domain w k , 1 s ( i ) = [ w k , l , 1 s ( i ) , w k , l , 2 s ( i ) , . . . , w k , l , N s ( i ) ] , The spatial domain of finishing the 1st sub-carrier signal of user k merges, and obtains the signal after the spatial domain merges z k , l ( i ) = w k , l s ( i ) &CenterDot; y k , l ( i ) ; Simultaneously with N matched filtering output signal y K, 2,1(i), y K, 2,2(i) ..., y K, 2, N(i) send into the 2nd the pairing spatial domain of subcarrier and merge module (3-2), obtain vector y K, 2(i)=[y K, 2,1(i), y K, 2,2(i) ..., y K, 2, N(i)] T, merge weight vector by the spatial domain w k , 2 s ( i ) = [ w k , 2,1 s ( i ) , w k , 2,2 s ( i ) , . . . , w k , 2 , N s ( i ) ] , The spatial domain of finishing the 2nd sub-carrier signal of user k merges, and obtains the signal after the spatial domain merges z k , 2 ( i ) = w k , 2 s ( i ) &CenterDot; y k , 2 ( i ) , . . . , Simultaneously with N matched filtering output signal y K, L, 1(i), y K, L, 2(i) ..., y K, L, N(i) send into L the pairing spatial domain of subcarrier and merge module (3-L), obtain vector y K, L(i)=[y K, L, 1(i), y K, L, 2(i) ..., y K, L, N, (i)] T, merge weight vector by the spatial domain w k , L s ( i ) = [ w k , L , 1 s ( i ) , w k , L , 2 s ( i ) , . . . , w k , L , N s ( i ) ] , The spatial domain of finishing L sub-carrier signal of user k merges, and obtains the signal after the spatial domain merges z k , L ( i ) = w k , L s ( i ) &CenterDot; y k , L ( i ) ;
Step 4 is with L sub-carrier signal z after the merging of spatial domain K, 1(i), z K, 2(i) ..., z K, L(i) send into carrier wave and merge module (4), obtain vector z k(i)=[z K, 1(i), z K, 2(i) ..., z K, L(i)] T, merge weight vector by frequency domain w k f ( i ) = [ w k , 1 f ( i ) , w k , 2 f ( i ) , . . . , w k , L f ( i ) ] Finish the merging of multi-carrier signal, obtain the final judgment variables of user k z k ( i ) = w k f ( i ) &CenterDot; z k ( i ) ;
Step 5 is with final judgment variables Z k(i) send into signal decision module (5) and carry out polarity decision, obtain the court verdict of arbitrary user k signal;
ω in the step in the above lBe angular frequency (1≤l≤L), the τ of l subcarrier of user K, lBe the time delay of l subcarrier of user k, T sBe bit period, [c K, 1 Ic K, 2 I... c K, L I] and [c K, 1 Qc K, 2 Q... c K, L Q] be respectively the I road of user k and the frequency expansion sequence of Q road correspondence, c wherein K, l I, c K, l QSpreading code for user k l subcarrier I road and Q road.
2. according to the described a kind of array antenna MC-CDMA system up-link receiving method of claim 1, it is characterized in that spatial domain in the step 3 merges to be based on suboptimum and to merge that weight vector carries out that the weight vector of l carrier signal spatial domain merging of user k is w k , l s ( i ) = &alpha; a k , l H ( i ) ;
Wherein, α K, l HBe the Hermite conversion of l subcarrier array response vector of user k, α is a constant, is taken as 1.
3. according to the described a kind of array antenna MC-CDMA system up-link receiving method of claim 1, it is characterized in that the carrier wave merging in the step 4 is undertaken by the high specific merging criterion, k user's multi-carrier signal merges weight vector and is w k f ( i ) = [ | a k , l ( i ) | , | a k , 2 ( i ) | , . . . , | a k , L ( i ) | ] ;
Wherein, α K, 1(i), α K, 2(i) ..., α K, L(i) be user k the 1st, 2 respectively ..., the array response vector of L subcarrier.
CNB2005100864156A 2005-09-13 2005-09-13 An array antenna MC-CDMA system up-link receiving method Expired - Fee Related CN100358261C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CNB2005100864156A CN100358261C (en) 2005-09-13 2005-09-13 An array antenna MC-CDMA system up-link receiving method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CNB2005100864156A CN100358261C (en) 2005-09-13 2005-09-13 An array antenna MC-CDMA system up-link receiving method

Publications (2)

Publication Number Publication Date
CN1741414A CN1741414A (en) 2006-03-01
CN100358261C true CN100358261C (en) 2007-12-26

Family

ID=36093651

Family Applications (1)

Application Number Title Priority Date Filing Date
CNB2005100864156A Expired - Fee Related CN100358261C (en) 2005-09-13 2005-09-13 An array antenna MC-CDMA system up-link receiving method

Country Status (1)

Country Link
CN (1) CN100358261C (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101106412B (en) * 2007-07-23 2010-09-08 北京交通大学 MC-DS-CDMA uplink receiving method and device for array antenna
CN103684687B (en) * 2013-11-18 2017-02-01 北京交通大学 Mine roadway uplink MC-CDMA wireless transmission cooperation method
US9900036B2 (en) * 2014-11-26 2018-02-20 Mediatek Singapore Pte. Ltd. Dual receive processing in wireless communications

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003169036A (en) * 2001-11-30 2003-06-13 Japan Telecom Co Ltd Orthogonal frequency-division multiplex system and transmitter/receiver
CN1489404A (en) * 2002-09-05 2004-04-14 ������������ʽ���� MC-CDMA down link transmitting method
WO2005011172A1 (en) * 2003-07-28 2005-02-03 Matsushita Electric Industrial Co., Ltd. Radio communication apparatus and radio communication method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2003169036A (en) * 2001-11-30 2003-06-13 Japan Telecom Co Ltd Orthogonal frequency-division multiplex system and transmitter/receiver
CN1489404A (en) * 2002-09-05 2004-04-14 ������������ʽ���� MC-CDMA down link transmitting method
WO2005011172A1 (en) * 2003-07-28 2005-02-03 Matsushita Electric Industrial Co., Ltd. Radio communication apparatus and radio communication method

Also Published As

Publication number Publication date
CN1741414A (en) 2006-03-01

Similar Documents

Publication Publication Date Title
CN101232484B (en) Signal transmission method, apparatus and communication system
JP5230779B2 (en) Wireless communication apparatus and wireless communication method
CN101385249B (en) Pattern diversity to support a mimo communications system and associated methods
JP5230766B2 (en) Arrival direction estimation apparatus and arrival direction estimation method
US8379620B2 (en) Transmission apparatus and transmission method using a plurality of divided frequency bands in a communication band
CN101083647B (en) Method for realizing synchronization in multi-input multi-output OFDM system
CN1941663B (en) Multi-antenna channel duplicating method and wavebeam shaping method
EP1985125A2 (en) Method and apparatus for performing spatial-division multiple access
CN101296011B (en) Self-adapting accidental beam mode selection method of wireless cognition network
CN101765123B (en) Method for distributing uplink resource of multi-user MIMO-OFDM
CN100358257C (en) Multi-antenna transmitting receiving system and method using random beam forming
CN100358261C (en) An array antenna MC-CDMA system up-link receiving method
CN105939541B (en) The equivalent subcarrier distribution system and method for multi-user&#39;s time frequency coding diversity MT-CDMA
CN100512083C (en) Uplink receiving method for array antenna MT-CDMA system
CN101170340B (en) Sub-channel distribution method for low-complexity, multi-user and multi-antenna OFDM system
Sun et al. A framework to construct three-dimensional complementary codes for multiuser MIMO systems
CN101106412B (en) MC-DS-CDMA uplink receiving method and device for array antenna
CN108768471A (en) The MIMO-OFDM-CDMA spectrum spreading methods of Combination selection mapping
CN103209496B (en) Based on the user scheduling method of beam selection in multi-base station cooperative system
CN101989867A (en) Cooperative communication method and system, base station and mobile terminal device
Alimi et al. Performance analysis of bit-error-rate and channel capacity of MIMO communication systems over multipath fading channels
CN101309249A (en) Uplink transmitting and receiving method for space multiplexed MIMO MT-CDMA system
Shrivastava et al. Iterative random beamforming for MIMO-OFDM systems
Ng et al. Multicarrier spread space-spectrum multiple access for the MIMO forward link transmission
Tian et al. Code Index Modulation CSK Spread Spectrum Underwater Acoustic Communications

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
C17 Cessation of patent right
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20071226

Termination date: 20100913