CN101834631B - Method and device for testing uplink synchronous maintenance time characteristic - Google Patents

Method and device for testing uplink synchronous maintenance time characteristic Download PDF

Info

Publication number
CN101834631B
CN101834631B CN 201010173937 CN201010173937A CN101834631B CN 101834631 B CN101834631 B CN 101834631B CN 201010173937 CN201010173937 CN 201010173937 CN 201010173937 A CN201010173937 A CN 201010173937A CN 101834631 B CN101834631 B CN 101834631B
Authority
CN
China
Prior art keywords
speed
data sequence
mobile terminal
sampling
correlation
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.)
Active
Application number
CN 201010173937
Other languages
Chinese (zh)
Other versions
CN101834631A (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 STARPOINT TECHNOLOGY Co Ltd
Original Assignee
BEIJING STARPOINT TECHNOLOGY Co Ltd
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 STARPOINT TECHNOLOGY Co Ltd filed Critical BEIJING STARPOINT TECHNOLOGY Co Ltd
Priority to CN 201010173937 priority Critical patent/CN101834631B/en
Publication of CN101834631A publication Critical patent/CN101834631A/en
Application granted granted Critical
Publication of CN101834631B publication Critical patent/CN101834631B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Synchronisation In Digital Transmission Systems (AREA)

Abstract

The invention discloses a method and a device for testing a TD-SCDMA (Time Division-Synchronization Code Division Multiple Access) uplink synchronous maintenance time characteristic, wherein the method comprises the following steps of: receiving a signal sent by a mobile terminal at the nth frame and sampling at a speed; carrying out interpolation operation on a sampled data sequence so as to obtain a data sequence of b speed; dividing the data sequence of the b speed into b groups, and carrying out relevance on b groups of data sequences respectively by utilizing local complex intermediate codes to obtain b correlation peak values and the positions of the b correlation peak values; fitting the secondary correlation curve of the data sequence of the b speed according to the b correlation peak values and the positions of the b correlation peak values, and working out the peak value position P1 of the secondary correlation curve; carrying out sampling, interpolation, relevance and fitting on a signal sent by the mobile terminal at the n+mth frame by a receiving end so as to work out the peak value position P2 of the secondary correlation curve, and estimating position offset tau according to the peak value position P1 and the peak value position P2; and comparing the position offset tau with real sending time delay tau' of the mobile terminal so as to work out the error between the position offset tau and the real sending time delay tau'.

Description

The method and apparatus of testing uplink synchronous maintenance time characteristic
Technical field
The present invention relates to TD-SCDMA mobile communication technology field, more specifically, relate to a kind of method and apparatus of the TD-SCDMA of test uplink synchronous maintenance time characteristic.
Background technology
The current epoch that entered 3G technology and application fast development of mobile communication.Time Division-Synchronous Code Division Multiple Access access (the Time Division-Synchronization Code DivisionMultiple Access of China's independent development, TD-SCDMA) the extensive commercialization of standard has critical strategic importance for the communications industry of China, wherein, the measuring technology of mobile terminal is having a strong impact on again the large-scale production of TD-SCDMA terminal, therefore becomes the focus in industrial chain.
TD-SCDMA is time division duplex (Time Division Duplex, TDD) pattern, it is the form transmission data with burst at a plurality of time slots, the multiple access that not exclusively synchronously can produce between different user disturbs, therefore Timing Synchronization is exactly basis and the prerequisite of whole system normal operation, has very important effect.In the TD-SCDMA system, uplink synchronous is divided into uplink synchronous foundation (Establishment of Uplink Synchronization) and uplink synchronous is safeguarded (Maintenance of Uplink Synchronization) two processes (specifically, referring to the introduction of 5.2.1 chapters and sections in " 2GHzTD-SCDMA digital mobile cellular telecommunication net height speed upstream packet access (HSUPA) Uu interface physical-layer techniques requires the 5th part: physical layer procedure ", the 3GPP agreement of correspondence is the joint of 5.1 in " 25224-840 ").On the basis that mobile terminal is set up in uplink synchronous, accurately adjust the synchronization accuracy of upward signal by the uplink synchronous maintenance process, so the time precision that uplink synchronous is safeguarded has determined the highest synchronization accuracy of whole system upward signal to a great extent.
Along with TD-SCDMA being strengthened to technology high speed downlink grouped access (High SpeedDownlink Package Access, HSDPA), high speed uplink packet access (High SpeedUplink Packet Access, HSUPA), multimedia broadcast multi-broadcasting business (MultimediaBroadcast and Multicast Services, MBMS), the research of HSPA+ (HSPA+ is the downward evolution version of HSPA (3GPPR6)) and application gos deep into, transmission rate and channel utilization are more and more higher, also more and more stricter to the requirement of the time response of mobile terminal synchronization.Traditional uplink synchronous maintenance time characteristic method of testing is by directly data sampling rate being improved and carries out the relevant lifting that obtains measuring accuracy to intermediate code (Midamble) again, the occupancy of its operand and memory space is all very large, and precision improvement is higher, the cost spent and the resource expended are also just larger.Simultaneously, the method time delay is obvious, to test, brings larger error.Therefore, how effectively to improve the study hotspot that uplink synchronous maintenance test precision and efficiency become TD-SCDMA terminal test in situation cheaply.
Summary of the invention
The technical problem that the present invention will solve is to provide a kind of method of the TD-SCDMA of test uplink synchronous maintenance time characteristic, can effectively improve measuring accuracy with lower cost.
The invention provides a kind of method of the TD-SCDMA of test uplink synchronous maintenance time characteristic, comprise: between mobile terminal and receiving terminal in the uplink synchronous maintenance process, the signal that the receiving terminal mobile terminal receive sends at the n frame, and signal mobile terminal sent with a speed is sampled, wherein, 0≤n<MAX_SFN; By interpolation method, the data sequence after sampling is carried out to interpolative operation to obtain the data sequence with respect to the b speed of mobile terminal transmitted signal, wherein, b>a; The data sequence of b speed is divided into to the b group, and the multiple intermediate code in this locality of utilizing receiving terminal is carried out relevantly to b group data sequence respectively, obtains b correlation peak and b the relevant position of correlation peak in the data sequence of b speed; The secondary correlation curve of the data sequence of the relevant position matching b speed in the data sequence of b speed according to b correlation peak and b correlation peak, and calculate the peak P1 of secondary correlation curve; The signal that receiving terminal sends at the n+m frame mobile terminal is sampled, interpolation, relevant and matching to be to calculate the peak P2 of secondary correlation curve, and according to the peak P1 of secondary correlation curve and the peak P2 estimated position side-play amount τ of secondary correlation curve, wherein, m>0; The true transmission delay τ ' of the position offset τ that receiving terminal is estimated and mobile terminal compares the error with the true transmission delay τ ' of calculating location side-play amount τ and mobile terminal.
According to an embodiment of the inventive method, the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code that the data sequence of carrying out interpolative operation is the sampling of a speed and the sampling of a speed, wherein, the length that L is associated window.
According to another embodiment of the inventive method, interpolation method is: in time domain, by the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code to the sampling of a speed and the sampling of a speed, carry out zero padding to realize interpolative operation; Or before and after the multiple intermediate code of sampling by multiple intermediate code and a speed of the sampling of a speed to after Fourier transform in frequency domain, the data sequence of each L/2 chip is carried out zero padding to realize interpolative operation.
According to the another embodiment of the inventive method, according to b correlation peak and b correlation peak, the step of the secondary correlation curve of the data sequence of the relevant position matching b speed in the data sequence of b speed comprises: structure secondary correlative curve model; Utilize the coefficient of least square method relevant position calculating secondary correlative curve model in the data sequence of b speed according to b correlation peak and b correlation peak.
According to an embodiment again of the inventive method, obtain the true transmission delay τ ' of mobile terminal according to following step: extract the uplink synchronous control command word that receiving terminal sends to mobile terminal between n frame and n+m frame; True transmission delay τ ' according to uplink synchronous control command word and uplink synchronous step size computation mobile terminal.
The method of testing uplink synchronous maintenance time characteristic provided by the invention, overcome operand and the large problem of memory space occupancy while reaching the conventional method equal accuracy by the packet of b speed being carried out to relevant method, and effectively dwindled time delay.The method can also further improve measuring accuracy by the matching computing simultaneously.
Another technical problem that the present invention will solve is to provide a kind of device of the TD-SCDMA of test uplink synchronous maintenance time characteristic, can effectively improve measuring accuracy with lower cost.
The invention provides a kind of device of the TD-SCDMA of test uplink synchronous maintenance time characteristic, comprise: receive and sampling module, for the signal sent at uplink synchronous maintenance process mobile terminal receive, and signal mobile terminal sent with a speed is sampled; Interpose module, be connected with sampling module with reception, for by interpolation method, the data sequence after sampling being carried out to interpolative operation to obtain the data sequence with respect to the b speed of mobile terminal transmitted signal, wherein, b>a; Correlation module, with interpose module, be connected, for the data sequence of b speed is divided into to the b group, the multiple intermediate code in this locality of use device is carried out relevantly to b group data sequence respectively, obtains b correlation peak and b the relevant position of correlation peak in the data sequence of b speed; Matching and detection module, be connected with correlation module, for the secondary correlation curve in the data sequence of the relevant position of the data sequence of b speed matching b speed according to b correlation peak and b correlation peak, and calculates the peak of secondary correlation curve; The position offset estimation module, with matching, with detection module, be connected, the peak P2 of the secondary correlation curve of the signal that the peak P1 of the secondary correlation curve of the signal sent at the n frame for the mobile terminal according to simulating and the mobile terminal simulated send at the n+m frame estimates position offset τ, wherein, 0≤n<MAX_SFN, m>0; Comparison module, be connected with the position offset estimation module, for the true transmission delay τ ' by position offset τ and mobile terminal, compares the error with the true transmission delay τ ' of calculating location side-play amount τ and mobile terminal.
An embodiment of the apparatus according to the invention, the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code that the data sequence that interpose module carries out interpolative operation is the sampling of a speed and the sampling of a speed, wherein, the length that L is associated window.
According to another embodiment of apparatus of the present invention, the interpolation method that interpose module adopts is: in time domain, by the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code to the sampling of a speed and the sampling of a speed, carry out zero padding to realize interpolative operation; Or before and after the multiple intermediate code of sampling by multiple intermediate code and a speed of the sampling of a speed to after Fourier transform in frequency domain, the data sequence of each L/2 chip is carried out zero padding to realize interpolative operation.
According to the another embodiment of apparatus of the present invention, matching and detection module comprise: modeling unit, for constructing the secondary correlative curve model; The matching unit, be connected with modeling unit, for utilizing least square method, according to b correlation peak and b correlation peak, calculates the coefficient of secondary correlative curve model in the relevant position of the data sequence of b speed; Detecting unit, be connected with the matching unit, the peak of the secondary correlation curve gone out for digital simulation.
An embodiment again according to apparatus of the present invention, this device also comprises true transmission delay computing module, true transmission delay computing module comprises: the command word extraction unit, for being extracted in the uplink synchronous control command word sent to mobile terminal between n frame and n+m frame; The time-delay calculation unit, be connected with the command word extraction unit, for the true transmission delay τ ' according to uplink synchronous control command word and uplink synchronous step size computation mobile terminal.
The device of testing uplink synchronous maintenance time characteristic provided by the invention, by correlation module, the packet of b speed is carried out to relevant method and overcome operand and the large problem of memory space occupancy while reaching the conventional method equal accuracy, and effectively dwindled time delay.The method can also be carried out the matching operation further to improve measuring accuracy by matching and detection module simultaneously.
The accompanying drawing explanation
Accompanying drawing described herein is used to provide a further understanding of the present invention, forms the application's a part.In the accompanying drawings:
Fig. 1 is the schematic flow sheet of an embodiment of the inventive method.
Fig. 2 is the schematic flow sheet of the another embodiment of the inventive method.
Fig. 3 is the shock response curve schematic diagram of raised cosine filter.
Fig. 4 is the structure of time slot schematic diagram of TD-SCDMA upward signal.
Fig. 5 is the divide into groups distribution schematic diagram of correlation curve sampled point of the present invention.
Fig. 6 is the divide into groups position distribution schematic diagram of relevant peaks of the present invention.
Fig. 7 is the matching correlation curve schematic diagram of the present invention's 16 speed data.
Fig. 8 is the structural representation of an embodiment of apparatus of the present invention.
Fig. 9 is the structural representation of another embodiment of apparatus of the present invention.
Figure 10 is the structural representation of the another embodiment of apparatus of the present invention.
Embodiment
With reference to the accompanying drawings the present invention is described more fully, exemplary embodiment of the present invention wherein is described.Exemplary embodiment of the present invention and explanation thereof are used for explaining the present invention, but do not form inappropriate limitation of the present invention.
High-accuracy high-efficiency rate TD-SCDMA uplink synchronous maintenance time characteristic method of testing provided by the invention has very important significance in the uniformity test of terminal.For example, the present invention can be used for to the response test of mobile terminal to uplink synchronization command word, can also apply the present invention in the robustness test of uplink synchronous.It should be noted that the present invention is not limited in above-mentioned two kinds of application, above-mentioned application is exemplary illustration.
The present invention is after the signal to the mobile terminal transmission is sampled, and mainly process interpolation of data, grouping slip are correlated with, secondary is correlated with curve and four steps of time difference calculating.Wherein, interpolation of data is to be reduced to higher speed by receiving signal, thereby provides favourable guarantee for the raising of conic fitting precision.Grouping is relevant is after the data after interpolation are divided into groups, and utilizes local multiple intermediate code He Ge road grouped data to be correlated with respectively and tries to achieve the relevant Ge of grouping road correlation peak location.The relevant curve of secondary is the relevant peaks of the relevant Ge of grouping road signal to be fitted to the correlation curve of whole signal.Time difference is calculated the side-play amount of the correlation peak location that at first records twice signal and the actual delay of mobile terminal transmitted signal, utilizes the corresponding relation of spreading rate and time to obtain uplink synchronous maintenance time characteristic.
Fig. 1 is the schematic flow sheet of an embodiment of the inventive method.
As shown in Figure 1, the method can comprise the following steps:
S102, between mobile terminal and receiving terminal in the uplink synchronous maintenance process, the signal that the receiving terminal mobile terminal receive sends at the n frame, and with a speed (wherein, the a speed refers to chip sample a time) signal that mobile terminal is sent sampled, wherein, 0≤n<MAX_SFN, the maximum subframe numbers that MAX_SFN is the TD-SCDMA system.
S104, by interpolation method, the data sequence after sampling is carried out to interpolative operation to obtain the data sequence with respect to the b speed of mobile terminal transmitted signal, wherein, b>a, for example, can be through an interpolation filter in time domain, by the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code to the sampling of a speed and the sampling of a speed, to carry out zero padding to realize interpolative operation; Perhaps can moving by frequency spectrum,, before and after the multiple intermediate code of sampling by multiple intermediate code and a speed of the sampling of a speed to after Fourier transform in frequency domain, the data sequence of each L/2 chip is carried out zero padding to realize interpolative operation, the length that L is associated window, preferably, it can be 8 or 16 chips.
S106, be divided into the b group by the data sequence of b speed, and the multiple intermediate code in this locality of utilizing receiving terminal is carried out relevantly to b group data sequence respectively, obtains b correlation peak and b the relevant position of correlation peak in the data sequence of b speed.
S108, the secondary correlation curve of the data sequence of the relevant position matching b speed in the data sequence of b speed according to b correlation peak and b correlation peak, and calculate the peak P1 of secondary correlation curve, for example, at first can construct the secondary correlative curve model, the recycling least square method is calculated the coefficient of secondary correlative curve model according to b correlation peak and b the relevant position of correlation peak in the data sequence of b speed, thereby can obtain the secondary correlation curve after matching.
S110, the signal that receiving terminal sends at the n+m frame mobile terminal repeats the peak P2 that above-mentioned sampling, interpolation, relevant and matching operate to calculate another secondary correlation curve, and according to the peak P1 of secondary correlation curve and the peak P2 estimated position side-play amount τ of secondary correlation curve,, τ=| P1-P2|, wherein, m>0.
S112, the true transmission delay τ ' of the position offset τ that receiving terminal is estimated and mobile terminal compares the error that calculates the true transmission delay τ ' of position offset τ and mobile terminal with the corresponding relation according to spreading rate and time, for example, can obtain according to following step the true transmission delay τ ' of mobile terminal: extract receiving terminal uplink synchronous control command word to the mobile terminal transmission between n frame and n+m frame, then according to the true transmission delay τ ' of uplink synchronous control command word and uplink synchronous step size computation mobile terminal.
This embodiment overcomes operand and the large problem of memory space occupancy while reaching the conventional method equal accuracy by the packet of b speed being carried out to relevant method, and has effectively dwindled time delay.The method can also further improve measuring accuracy by the matching computing simultaneously.
Owing to needing at first to carry out the foundation of uplink synchronous in the TD-SCDMA system, carry out interpolation and related operation so can only choose the data sequence of multiple intermediate code and correlation window length sum, can guarantee relevant peaks can be detected in correlated process.
In another embodiment of the inventive method, the data sequence of carrying out interpolative operation can be the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code of a speed sampling and the sampling of a speed, wherein, and the length that L is associated window.Owing to only above-mentioned very short data sequence being carried out to interpolative operation, compared to existing technology, without the data to whole time slot, carry out interpolative operation, improve to a great extent arithmetic speed, and significantly saved memory space, reduced and realized cost.
Fig. 2 is the schematic flow sheet of the another embodiment of the inventive method.
As shown in Figure 2, this embodiment can comprise the following steps:
S202, the signal that the receiving terminal mobile terminal receive sends at the n frame, wherein, 0≤n<MAX_SFN;
S204, the signal that receiving terminal sends mobile terminal with lower a speed is sampled, and a is less than the inverse of the resolution of system requirements usually, for example, the uplink synchronous precision that supposing the system requires is 1/16, and a is less than 16 usually, can reach by following interpolative operation 1/16 precision;
S206, use interpolation method to obtain the b speed data sequence of higher corresponding mobile terminal transmitted signal to the data sequence after sampling, wherein, and b>a;
S208, the data sequence of b speed is divided into to the b group, utilize the multiple intermediate code (Midamble code) in this locality of receiving terminal to slide respectively relevant in associated window to above-mentioned b group data sequence, obtain the b group data relevant sequence of correlation values of sliding, find respectively organize sequence of correlation values peak value (, and the correspondence position of b relevant peaks in b speed data sequence maximum in sequence of correlation values);
S210, the curve formed due to position and the peak value of b relevant peaks says the impulse response of raised cosine filter from the strict sense, its time-domain expression is
Figure GSA00000122947400081
wherein, t is the time, T sfor symbol period, α is roll-off factor, 0≤α≤1, its curve as shown in Figure 3, as can be seen from Figure 3, near the curve impulse response peak value of this raised cosine filter can be approximated to be conic section, thereby, can utilize least square method b grouped data correlation peak and position to be simulated to the secondary correlation curve of b speed data, for approaching the shock response curve of above-mentioned raised cosine filter, and obtain the vertex position P1 of the secondary correlation curve of b speed data, turn S214;
S212, the signal that the receiving terminal mobile terminal receive sends at the n+m frame, repeat above-mentioned S204~S210 step, obtains the vertex position P2 of the secondary correlation curve of another b speed data, turns to S214, wherein, m>0;
S214, compare the position P1 of two relevant peaks and the position offset of P2, mapping relations according to spreading rate and time, try to achieve the time delay of receiving terminal, and contrasted with the true transmission delay of mobile terminal, calculate the error of the true transmission delay of position offset and mobile terminal, thereby obtain uplink synchronous maintenance time characteristic.
Preferably, a can get 4, b and can get 16.
Preferably, in order to test time migration in a big way, the length that can choose associated window is 16 chips.
In addition, can adopt least square method to carry out the matching of secondary correlation curve.For example, can construct a ternary linear function group, utilize the displace analysis equation group of b group correlation peak and correlation peak, obtain the parameter of the quadratic equation with one unknown of secondary correlation curve.
This embodiment has reached the requirement of high measuring accuracy by the matching computing, conventional method operand and the larger problem of memory space occupancy when reaching high accuracy have also been overcome simultaneously, and time delay is little, it is a kind of TD-SCDMA uplink synchronous maintenance time characteristic method of testing of high-accuracy high-efficiency rate.
The instantiation of below take is example, describes above-mentioned each processing procedure in detail, wherein, suppose, and a=4, b=16, L=16:
(1) interpolation of data
Fig. 4 is the structure of time slot schematic diagram of TD-SCDMA upward signal.
As shown in Figure 4, the uplink service time slot comprises data field 1, data field 2 and Midamble code territory, totally 864 chips (chip).If the length of associated window is 16 chip, the data sequence of carrying out interpolative operation is: front 8 chip of last 8 chip of data field 1, Midamble code territory (144 chip) and data field 2, be total to 144+8+8=160 chip, as shown in the empty frame in Fig. 4.
Define symbol is as follows:
Data sequence after 4 speed samplings is: { x (n), 0≤n≤160 * 4-1};
The data sequence that data sequence after 4 speed samplings is carried out after Fourier transform is: { X (k), 0≤k≤160 * 4-1};
Data sequence after 16 speed fast Fourier transforms is: { Y (k), 0≤k≤160 * 16-1};
The time domain data sequence of 16 speeds is: { y (n), 0≤n≤160 * 16-1};
At first the data sequence x (n) of 4 speeds carried out to fast Fourier transform, obtain X (k), X (k) is exactly the frequency spectrum reflection of x (n), inserts 160 * 4 * 30 to realize the time domain interpolation to x (n) in X (k) sequence,
Y ( k ) = X ( k ) ( 0 &le; k &le; 80 &times; 4 - 1 ) Y ( k ) = 0 ( 80 &times; 4 &le; k &le; 560 &times; 4 - 1 ) Y ( k ) = X ( k - 560 &times; 4 ) ( 560 &times; 4 &le; k &le; 640 &times; 4 - 1 )
Again Y (k) is carried out to inverse Fourier transform, can obtain the data sequence y (n) of 16 speeds.
(2) grouping is slided relevant
Define symbol is as follows:
Every group of data sequence after 16 speed data sequence groupings is: y ' and (n), 0≤n≤160-1};
Local multiple intermediate code is: { M (n), 0≤n≤144-1};
Sequence after grouping is relevant is: { p j(i), 0≤i≤15,0≤j≤15};
The position of grouping relevant peaks is: { A j(m j, n j), 0≤j≤15};
At first the data sequence y of 16 speeds (n) is divided into to 16 groups, , at every turn respectively from y (0)~y (15), take 16 as the interval extracted data, form new data sequence y ' (n), by the conjugate form of the multiple intermediate code M (n) of this locality, with the data sequence sampled is corresponding, multiply each other, and each product addition is obtained to correlation, and then sampled data is shifted, carry out new related operation once, altogether slide 16 times (, the length of associated window), carry out respectively 16 correlations computings (, at first choose sampled data 1~144, to M (n) do 144 relevant, and then choose 2~145 sampled datas and carry out relevant to M (n), by that analogy, carry out altogether sliding for 16 times relevant), thereby obtain the grouping correlation p of j group j(i), as shown in Figure 5.
Next, in these 16 correlations, choose maximum one (that is, the relevant peaks of every group), the position in its data sequence y in 16 speeds (n) and corresponding correlation are stored in to the position A of grouping relevant peaks j(m j, n j) in (wherein, m jrepresent position, n jfor correlation peak), as shown in Figure 6.For example, this peak m jcan be m j=b*pos j+ index j, wherein, pos jbe the position that j organizes the relevant peaks of single speed, 0≤pos j<L, index jbe the index that j organizes the b speed that single speed data extract jthe road signal, 0≤index j<b.
(3) matching of secondary correlation curve
Through the preparation of front two steps, we have obtained the position A of 16 grouping relevant peaks j(m j, n j), following problem is exactly how to utilize these 16 points to obtain the correlation peak location of 16 haplotype data sequences with maximal accuracy, the present invention can adopt least square method to carry out the matching of secondary correlation curve, it should be noted that, the present invention can adopt multiple curve-fitting method, be not limited in least square method, the following least square method provided is only for exemplary illustration.Curve, the most frequently usedly simulate the curve that can represent these data characteristics by one group of data.The correlation curve of multiple intermediate code is approximately a conic section, if the data point that first two steps provide fails accurately to drop on this parabola, is the error due to experimental data.
Definition secondary correlative curve model, and Coefficient of determination a 0, a 1and a 2:
y=a 0+a 1x+a 2x 2
From background geometry, considering, is exactly will be with a 0, a 1and a 2for undetermined coefficient, determine that conic section makes 16 corresponding 16 points of grouping relevant peaks be positioned as close to this secondary correlation curve.In general, data point can all not drop on this curve, if the data of j point drop on curve just, the coordinate of this point meets quadratic curve equation, that is,
n j=a 0+a 1m j+a 2m j 2
If this point does not drop on curve, its coordinate does not meet above-mentioned curvilinear equation, can have certain error (that is, residual error).So all the overall error at some place can be expressed as by residual sum of squares (RSS):
F ( a 0 , a 1 , a 2 ) = &Sigma; j = 1 16 [ ( a 0 + a 1 m j + a 2 m j 2 ) - n j &prime; ] 2 ;
When j point by chance drops on this curve, n j=n j'; When j point do not drop on this curve, n j' the measured value that goes out for Practical Calculation, and
Figure GSA00000122947400112
Above-mentioned F (a0, a1, a2) is about a 0, a 1and a 2three meta-functions, can choose a 0, a 1and a 2, make this function minimalization.In order to ask the minimum point of this function, can make:
Figure GSA00000122947400113
Figure GSA00000122947400114
Figure GSA00000122947400115
Obtain:
&Sigma; j = 1 16 2 [ ( a 0 + a 1 m j + a 2 m j 2 ) - n j &prime; ] = 0 &Sigma; j = 1 16 2 [ ( a 0 + a 1 m j + a 2 m j 2 ) - n j &prime; ] m j = 0 &Sigma; j = 1 16 2 [ ( a 0 + a 1 m j + a 2 m j 2 ) - n j &prime; ] m j 2 = 0
This is about undetermined coefficient a 0, a 1and a 2system of linear equations, solve this equation group and can obtain three undetermined coefficient a in the quadratic fit function 0, a 1and a 2thereby the curve after matching as shown in Figure 7.
The conic section simulated according to said method is very large with respect to discrete point set precision improvement, and the precision promoted is equivalent to analog-and digital-difference.
(4) calculating of time difference
The abscissa position at the top of conic section is the correlation peak location of 16 speed data sequences, and recording its side-play amount is P1 (only recording abscissa gets final product).Mobile terminal transmitted signal again after the m frame, repeat above-mentioned processing procedure, obtain relevant peaks side-play amount P2, try to achieve the difference Δ P of twice side-play amount, according to the corresponding relation of spreading rate and time (for example, spreading rate is 1.28Mcps), utilize m, Δ P, uplink synchronous control command word and synchronous step-length to try to achieve time error, be TD-SCDMA uplink synchronous maintenance time characteristic value.
Fig. 8 is the structural representation of an embodiment of apparatus of the present invention.
As shown in Figure 8, the device of this embodiment comprises reception and sampling module 11, and for the signal sent at uplink synchronous maintenance process mobile terminal receive, and signal mobile terminal sent with a speed is sampled; Interpose module 12, be connected with sampling module 11 with reception, for by interpolation method, the data sequence after sampling being carried out to interpolative operation to obtain the data sequence with respect to the b speed of mobile terminal transmitted signal, wherein, b>a; Correlation module 13, with interpose module 12, be connected, for the data sequence of b speed is divided into to the b group, the multiple intermediate code in this locality of use device is carried out relevantly to b group data sequence respectively, obtains b correlation peak and b the relevant position of correlation peak in the data sequence of b speed; Matching and detection module 14, with correlation module 13, be connected, for the secondary correlation curve in the data sequence of the relevant position of the data sequence of b speed matching b speed according to b correlation peak and b correlation peak, and calculate the peak of secondary correlation curve; Position offset estimation module 15, with matching, with detection module 14, be connected, the peak P2 of the secondary correlation curve of the signal that the peak P1 of the secondary correlation curve of the signal sent at the n frame for the mobile terminal according to simulating and the mobile terminal simulated send at the n+m frame estimates position offset τ, wherein, 0≤n<MAX_SFN, m>0; Comparison module 16, be connected with position offset estimation module 15, for the true transmission delay τ ' by position offset τ and mobile terminal, compares the error with the true transmission delay τ ' of calculating location side-play amount τ and mobile terminal.
This embodiment carries out relevant method by correlation module by the packet of b speed and overcomes operand and the large problem of memory space occupancy while reaching the conventional method equal accuracy, and has effectively dwindled time delay.The method can also be carried out the matching operation further to improve measuring accuracy by matching and detection module simultaneously.
Owing to needing at first to carry out the foundation of uplink synchronous in the TD-SCDMA system, carry out interpolation and related operation so can only choose the data sequence of multiple intermediate code and correlation window length sum, can guarantee relevant peaks can be detected in correlated process, so it can be the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code of a speed sampling and the sampling of a speed that interpose module carries out the data sequence of interpolative operation, wherein, the length that L is associated window.Owing to only above-mentioned data sequence being carried out to interpolative operation, compared to existing technology, without the data to whole time slot, carry out interpolative operation, improve to a great extent arithmetic speed, and significantly saved memory space.
For example, if the length of associated window is 16 chip, front 8 chip of last 8 chip, Midamble code territory (144 chip) and data field 2 that the data sequence of carrying out interpolative operation is data field 1 as shown in Figure 4,144+8+8=160 chip altogether.
Particularly, interpose module can adopt following method to carry out interpolative operation: in time domain, by the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code to the sampling of a speed and the sampling of a speed, carry out zero padding to realize interpolative operation; Or before and after the multiple intermediate code of sampling by multiple intermediate code and a speed of the sampling of a speed to after Fourier transform in frequency domain, the data sequence of each L/2 chip is carried out zero padding to realize interpolative operation.
Operand and the large problem of memory space occupancy when this embodiment overcomes conventional method arrival equal accuracy by the packet of b speed being carried out to relevant method, and effectively dwindled time delay.This device can also further improve measuring accuracy by matching and detection module simultaneously.
Fig. 9 is the structural representation of another embodiment of apparatus of the present invention.
As shown in Figure 9, with the embodiment in Fig. 8, compare, the matching in this embodiment and detection module 21 comprise: modeling unit 211, for constructing the secondary correlative curve model; Matching unit 212, be connected with modeling unit 211, for utilizing least square method, according to b correlation peak and b correlation peak, calculates the coefficient of secondary correlative curve model in the relevant position of the data sequence of b speed; Detecting unit 213, be connected with matching unit 212, the peak of the secondary correlation curve gone out for digital simulation.
For example, can be by modeling unit 211 definition secondary correlative curve model y=a 0+ a 1x+a 2x 2, by matching unit 212 according to b correlation peak and b the relevant position Coefficient of determination a of correlation peak in the data sequence of b speed 0, a 1and a 2.When matching unit 212 utilizes least square method to carry out curve fitting, function can be set
Figure GSA00000122947400131
and order
Figure GSA00000122947400132
Figure GSA00000122947400133
Figure GSA00000122947400134
thereby obtain coefficient a 0, a 1and a 2, can determine the secondary correlation curve.
The matching operation that this embodiment carries out by matching and detection module has further improved the measuring accuracy of uplink synchronous.
Figure 10 is the structural representation of the another embodiment of apparatus of the present invention.
As shown in figure 10, with the embodiment in Fig. 8, compare, device in this embodiment also comprises true transmission delay computing module 31, true transmission delay computing module 31 comprises: command word extraction unit 311, for being extracted in the uplink synchronous control command word sent to mobile terminal between n frame and n+m frame; Time-delay calculation unit 312, be connected with command word extraction unit 311, for the true transmission delay τ ' according to uplink synchronous control command word and uplink synchronous step size computation mobile terminal.
Description of the invention provides for example with for the purpose of describing, and is not exhaustively or limit the invention to disclosed form.Many modifications and variations are apparent for the ordinary skill in the art.Selecting and describing embodiment is for better explanation principle of the present invention and practical application, thereby and makes those of ordinary skill in the art can understand the various embodiment with various modifications that the present invention's design is suitable for special-purpose.

Claims (10)

1. a method of testing the TD-SCDMA uplink synchronous maintenance time characteristic, is characterized in that, described method comprises:
Between mobile terminal and receiving terminal in the uplink synchronous maintenance process, described receiving terminal receives the signal that described mobile terminal sends at the n frame, and the signal described mobile terminal sent with a speed is sampled, wherein, 0≤n<MAX_SFN, the maximum subframe numbers that MAX_SFN is the TD-SCDMA system;
By interpolation method, the data sequence after sampling is carried out to interpolative operation to obtain the data sequence with respect to the b speed of described mobile terminal transmitted signal, wherein, b>a;
The data sequence of described b speed is divided into to the b group, and the multiple intermediate code in this locality of utilizing described receiving terminal is carried out relevantly to b group data sequence respectively, obtains b correlation peak and the relevant position of a described b correlation peak in the data sequence of described b speed;
The secondary correlation curve of the data sequence of the described b speed of relevant position matching in the data sequence of described b speed according to a described b correlation peak and a described b correlation peak, and calculate the peak P1 of secondary correlation curve;
The signal that described receiving terminal sends at the n+m frame described mobile terminal carries out described sampling, interpolation, relevant and matching to calculate the peak P2 of secondary correlation curve, and according to the peak P1 of described secondary correlation curve and the peak P2 estimated position side-play amount τ of described secondary correlation curve, wherein, τ=| P1-P2|, m>0;
The true transmission delay τ ' of the position offset τ that described receiving terminal is estimated and described mobile terminal compares to calculate the error of the true transmission delay τ ' of described position offset τ and described mobile terminal.
2. method according to claim 1, is characterized in that, the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code that described data sequence of carrying out interpolative operation is the sampling of a speed and the sampling of described a speed, wherein, the length that L is associated window.
3. method according to claim 2, is characterized in that, described interpolation method is:
Carry out zero padding to realize interpolative operation by the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code to described a speed sampling and the sampling of described a speed in time domain; Or
The data sequence of each L/2 chip of multiple intermediate code front and back of sampling by multiple intermediate code and the described a speed of the described a speed sampling to after Fourier transform in frequency domain is carried out zero padding to realize interpolative operation.
4. method according to claim 1, is characterized in that, described according to b correlation peak and a described b correlation peak, the step of the secondary correlation curve of the data sequence of the described b speed of relevant position matching in the data sequence of described b speed comprises:
Structure secondary correlative curve model;
Utilize least square method to calculate the coefficient of described secondary correlative curve model in the relevant position in the data sequence of described b speed according to a described b correlation peak and a described b correlation peak.
5. method according to claim 1, is characterized in that, obtains the true transmission delay τ ' of described mobile terminal according to following step:
Extract the uplink synchronous control command word that described receiving terminal sends to described mobile terminal between n frame and n+m frame;
True transmission delay τ ' according to described uplink synchronous control command word and the described mobile terminal of uplink synchronous step size computation.
6. a device of testing the TD-SCDMA uplink synchronous maintenance time characteristic, is characterized in that, described device comprises:
Receive and sampling module, in the uplink synchronous maintenance process, receiving the signal that described mobile terminal sends, and the signal described mobile terminal sent with a speed is sampled;
Interpose module, be connected with sampling module with described reception, for by interpolation method, the data sequence after sampling being carried out to interpolative operation to obtain the data sequence with respect to the b speed of described mobile terminal transmitted signal, wherein, b>a;
Correlation module, with described interpose module, be connected, for the data sequence of described b speed is divided into to the b group, the multiple intermediate code in this locality of utilizing described device is carried out relevantly to b group data sequence respectively, obtains b correlation peak and the relevant position of a described b correlation peak in the data sequence of described b speed;
Matching and detection module, with described correlation module, be connected, for the secondary correlation curve in the data sequence of the described b speed of the relevant position of the data sequence of described b speed matching according to a described b correlation peak and a described b correlation peak, and calculate the peak of secondary correlation curve;
The position offset estimation module, with described matching, with detection module, be connected, the peak P2 of the secondary correlation curve of the signal that the peak P1 of the secondary correlation curve of the signal sent at the n frame for the described mobile terminal according to simulating and the described mobile terminal simulated send at the n+m frame estimates position offset τ, wherein, τ=| P1-P2|, 0≤n<MAX_SFN, m>0,0≤n<MAX_SFN, the maximum subframe numbers that MAX_SFN is the TD-SCDMA system;
Comparison module, be connected with described position offset estimation module, compare to calculate the error of the true transmission delay τ ' of described position offset τ and described mobile terminal for the true transmission delay τ ' by described position offset τ and described mobile terminal.
7. device according to claim 6, it is characterized in that, the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code that the data sequence that described interpose module carries out interpolative operation is the sampling of a speed and the sampling of described a speed, wherein, the length that L is associated window.
8. device according to claim 7, is characterized in that, the interpolation method that described interpose module adopts is:
Carry out zero padding to realize interpolative operation by the data sequence of each L/2 chip before and after the multiple intermediate code of the multiple intermediate code to described a speed sampling and the sampling of described a speed in time domain; Or
The data sequence of each L/2 chip of multiple intermediate code front and back of sampling by multiple intermediate code and the described a speed of the described a speed sampling to after Fourier transform in frequency domain is carried out zero padding to realize interpolative operation.
9. device according to claim 6, is characterized in that, described matching and detection module comprise:
Modeling unit, for constructing the secondary correlative curve model;
The matching unit, be connected with described modeling unit, for utilizing least square method, according to a described b correlation peak and a described b correlation peak, calculates the coefficient of described secondary correlative curve model in the relevant position of the data sequence of described b speed;
Detecting unit, be connected with described matching unit, the peak of the secondary correlation curve gone out for digital simulation.
10. device according to claim 6, is characterized in that, described device also comprises true transmission delay computing module, and described true transmission delay computing module comprises:
The command word extraction unit, for being extracted in the uplink synchronous control command word sent to described mobile terminal between n frame and n+m frame;
The time-delay calculation unit, be connected with described command word extraction unit, for the true transmission delay τ ' according to described uplink synchronous control command word and the described mobile terminal of uplink synchronous step size computation.
CN 201010173937 2010-05-10 2010-05-10 Method and device for testing uplink synchronous maintenance time characteristic Active CN101834631B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201010173937 CN101834631B (en) 2010-05-10 2010-05-10 Method and device for testing uplink synchronous maintenance time characteristic

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201010173937 CN101834631B (en) 2010-05-10 2010-05-10 Method and device for testing uplink synchronous maintenance time characteristic

Publications (2)

Publication Number Publication Date
CN101834631A CN101834631A (en) 2010-09-15
CN101834631B true CN101834631B (en) 2013-06-05

Family

ID=42718583

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010173937 Active CN101834631B (en) 2010-05-10 2010-05-10 Method and device for testing uplink synchronous maintenance time characteristic

Country Status (1)

Country Link
CN (1) CN101834631B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102647735B (en) * 2012-04-13 2014-11-12 上海创远仪器技术股份有限公司 Fast uplink synchronization method and device of TD-SCDMA (Time Division-Synchronous Code Division Multiple Access)
CN108937957B (en) * 2018-06-05 2021-11-09 武汉久乐科技有限公司 Detection method, device and detection equipment

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1968049A (en) * 2005-11-16 2007-05-23 大唐移动通信设备有限公司 Synchronization method and device after terminal dormancy awaken of TD_SCDMA mobile phone
CN101436900A (en) * 2008-12-25 2009-05-20 重庆重邮信科通信技术有限公司 Method and apparatus for scanning whole frequency band frequency point of advance system

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1968049A (en) * 2005-11-16 2007-05-23 大唐移动通信设备有限公司 Synchronization method and device after terminal dormancy awaken of TD_SCDMA mobile phone
CN101436900A (en) * 2008-12-25 2009-05-20 重庆重邮信科通信技术有限公司 Method and apparatus for scanning whole frequency band frequency point of advance system

Also Published As

Publication number Publication date
CN101834631A (en) 2010-09-15

Similar Documents

Publication Publication Date Title
CN102457870B (en) Method and device for detecting primary synchronization signal and method and system for searching neighborhoods
CN102413079B (en) Initial fraction frequency offset estimation method used in downlink of 3 generation partnership project long term evolution (3GPP-LTE) system
CN102298097B (en) Method for estimating thunder impulse signal Time Difference of Arrival (TDOA)
CN101489238B (en) Time difference measuring method, system and apparatus
CN101541078A (en) Method, system and device for estimating TDOA
CN105467361A (en) Radio signal complex envelope and carrier phase information combination BVR (beyond visual range) target direct positioning method
CN101814160A (en) RBF neural network modeling method based on feature clustering
CN101982953B (en) Frequency domain multi-dimensional parameterized model of broadband wireless communication channel and modeling method
CN101425855B (en) Wideband channel parameter extraction method, wideband channel simulation method and apparatus thereof
CN101834631B (en) Method and device for testing uplink synchronous maintenance time characteristic
EP3554138A1 (en) Method for determining candidate line-of-sight path, and wireless local area network device
CN102356671A (en) Method and device of communication system frequency search
CN101426213B (en) Wideband channel simulation method and apparatus thereof
CN103391182B (en) Frame detection method and device
CN101447823A (en) TD-SCDMA system down link synchronisation method
CN101611606B (en) Method for symbol synchronization of received digital signal and digital signal receiver using the same method
CN102664835B (en) Multicell channel estimation method based on interference elimination route by route
CN102983956A (en) Moving speed estimation method for terminal and base station
CN102083176B (en) Time synchronization method and device in cell search process
CN102111228B (en) Cognitive radio frequency spectrum sensing method based on circulation symmetry
CN106411793B (en) DMRS channel parameter estimation method, device and user terminal
CN100459446C (en) An estimation and search method for channel path of multi-path fading channel
CN110552682A (en) Descending analysis method for vertical well yield of carbonate rock triple medium gas reservoir
CN102664656B (en) Direct sequence spread spectrum signal compressed domain united synchronization and demodulation method based on structured sparsity
CN1863408B (en) Method and system for implementing common frequency measurement in TD-SCDMA system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB02 Change of applicant information

Address after: 100102 Beijing City, Chaoyang District Lize two Road No. 2 three floor, four floor of the South

Applicant after: BEIJING STARPOINT TECHNOLOGY CO., LTD.

Address before: 100083, Haidian District, Xueyuan Road, No. 40 research, 1-1 floor, two floor, Beijing

Applicant before: Xingheliangdian Communication Software Co., Ltd., Beijing

COR Change of bibliographic data

Free format text: CORRECT: APPLICANT; FROM: XINGHELIANGDIAN COMMUNICATION SOFTWARE CO LTD, BEIJING TO: BEIJING STARPOINT TECHNOLOGY CO., LTD.

C14 Grant of patent or utility model
GR01 Patent grant