CN104184524B - A kind of transmitting frequency calibration method and device - Google Patents

Abstract

The invention discloses a kind of transmitting frequency calibration method and device, relate to communication technology, after carrying out coarse frequency offset and the compensation of thick frequency deviation, data after thick frequency deviation compensates intercept two segment datas, determine the correlation of two intercepted segment datas, thin frequency deviation deviant can be determined according to correlation, obtain cell information without information alternately and be used for determining that accurate time slot original position carries out thin frequency deviation estimation, reduce the complexity of frequency calibration, improve testing efficiency.

Description

A kind of transmitting frequency calibration method and device

Technical field

The present invention relates to communication technology, particularly relate to a kind of transmitting frequency calibration method and device.

Background technology

Along with the fast development of mobile communication system, support that the terminal of multiple types is also developed rapidly, Manufacturer terminal has been developed for supporting GSM(Global System for Mobile communication at present, Global system for mobile communications), WCDMA(Wireless Code Division Multiple Access wireless code Point multiple access accesses), TD-SCDMA(time division synchronous CDMA), LTE(Long Term Evolution, Long Term Evolution) etc. the Related product of multiple types.Before these end products network, need it is carried out Uniformity test, and RF consistency therein test is most basic test therein, test case is also More.

At present terminal to these different systems carries out usually combining of RF consistency test Test equipment, is referred to as comprehensive test instrument.Comprehensive test instrument is typically except possessing the relevant association according to RF consistency test View carries out outside the test of relevant measurement item, in addition it is also necessary to possess quickly calibrated function, terminal is carried out frequency, Power and the calibration of gain.Calibration to frequency is commonly referred to as AFC(Automatic Frequency Control, Automatic frequency controls) calibration, generally require and large range of frequency error is estimated fast and accurately.

The AFC calibration function that comprehensive test instrument realizes at present typically requires knows that (these information include cell information Community ID(identifies), frequency pilot sign configuration etc.), reference marks, after knowing cell information, is entered by comprehensive test instrument Line reconstruction, utilizes local reference marks can dock the collection of letters number with the dependency received between reference marks and carries out essence When determining, reduce by the timing synchronization errors impact on Frequency Estimation, carry out frequency according to Timing Synchronization result again Inclined estimation.

Conventional AFC calibration need comprehensive test instrument after knowing community configured information, according to community configured information Restructural this locality reference marks, utilizes the local reference marks of reconstruct to carry out essence with the dependency receiving reference signal True Timing Synchronization, can improve measuring accuracy, but this mode adds information interaction time and test process Complexity, testing efficiency is relatively low.

Summary of the invention

The embodiment of the present invention provides a kind of transmitting frequency calibration method and device, to reduce the complexity of frequency calibration, Improve testing efficiency.

A kind of transmitting frequency calibration method, including:

Receive the phase-shift keying (PSK) psk modulation signal that terminal sends, and determine the leading edge position of described signal；

Sufficiently long data are intercepted according to leading edge position；

Thick frequency deviation deviant is determined according to the data intercepted, and by this thick frequency deviation deviant to being intercepted Data carry out thick frequency deviation compensation；

Data after thick frequency deviation compensates intercept two segment datas, determines the relevant of two intercepted segment datas Value, and determine thin frequency deviation deviant according to described correlation；

According to described thick frequency deviation deviant and the sum of described thin frequency deviation deviant, carry out frequency calibration.

A kind of frequency calibration device, including:

Rising edge determines unit, for receiving the psk modulation signal that terminal sends, and determines described signal Leading edge position；

Data cutout unit, for intercepting sufficiently long data according to leading edge position；

Coarse frequency offset unit is for determining thick frequency deviation deviant according to the data intercepted and thick by this The frequency deviation deviant data to being intercepted carry out thick frequency deviation compensation；

Thin frequency deviation estimating unit, intercepts two segment datas the data after compensating from thick frequency deviation, determines and cut The correlation of two segment datas taken, and determine thin frequency deviation deviant according to described correlation；

Frequency calibration unit, for according to described thick frequency deviation deviant and the sum of described thin frequency deviation deviant, enters Line frequency is calibrated.

The embodiment of the present invention provides a kind of transmitting frequency calibration method and device, is carrying out coarse frequency offset and thick frequency deviation After compensation, the data after thick frequency deviation compensates intercept two segment datas, determines the phase of two intercepted segment datas Pass is worth, and can determine thin frequency deviation deviant according to correlation, obtains cell information without information alternately and uses In determining that accurate time slot original position carries out thin frequency deviation estimation, reduce the complexity of frequency calibration, improve Testing efficiency.

Accompanying drawing explanation

The transmitting frequency calibration method flow chart that Fig. 1 provides for the embodiment of the present invention；

The more specific transmitting frequency calibration method flow chart that Fig. 2 provides for the embodiment of the present invention；

The transmitting frequency calibration method flow chart of the corresponding embodiment one that Fig. 3 provides for the embodiment of the present invention；

The transmitting frequency calibration method flow chart of the corresponding embodiment two that Fig. 4 provides for the embodiment of the present invention；

The frequency calibration apparatus structure schematic diagram that Fig. 5 provides for the embodiment of the present invention.

Detailed description of the invention

The embodiment of the present invention provides a kind of transmitting frequency calibration method and device, is carrying out coarse frequency offset and thick frequency deviation After compensation, the data after thick frequency deviation compensates intercept two segment datas, determines the phase of two intercepted segment datas Pass is worth, and can determine thin frequency deviation deviant according to correlation, obtains cell information without information alternately and uses In determining that accurate time slot original position carries out thin frequency deviation estimation, reduce the complexity of frequency calibration, improve Testing efficiency.

As it is shown in figure 1, the transmitting frequency calibration method that the embodiment of the present invention provides, including:

The PSK(Phase Shift Keying that step S101, reception terminal send, phase-shift keying (PSK)) modulation letter Number, and determine the leading edge position of signal；

Step S102, according to leading edge position intercept sufficiently long data；

Step S103, determine thick frequency deviation deviant according to the data intercepted, and by this thick frequency deviation deviant The data intercepted are carried out thick frequency deviation compensation；

Step S104, from thick frequency deviation compensate after data intercept two segment datas, determine two intercepted hop counts According to correlation, and determine thin frequency deviation deviant according to correlation；

Step S105, according to thick frequency deviation deviant and the sum of thin frequency deviation deviant, carry out frequency calibration.

In step S102, the data intercepted are the longest, and the result precision carrying out channel estimation is the highest, is carried The complexity come is the highest, and when the data that intercept are shorter, the result precision carrying out channel estimation is relatively low, multiple Miscellaneous degree is relatively low, and generally, in step S102, the data intercepted are at least the data of a time slot.

In step S104, it is not necessary to obtained pilot tone by information alternately or local training sequence determines standard True time slot original position, directly intercepts two segment datas the data after thick frequency deviation compensates, determines and intercepted The correlation of two segment datas, thin frequency deviation deviant can be determined according to correlation, reduce frequency calibration Complexity, improves testing efficiency.

Concrete, as in figure 2 it is shown, the method includes:

Step S201, host computer configurating terminal send M unit PSK(such as BPSK, QPSK, 8PSK Modulation etc.) modulated signal, the purpose being configured so that is to obtain the modulation of standard M unit PSK planisphere Signal, to facilitate removal modulation intelligence, thus carries out frequency deviation estimation；

Step S202, comprehensive test instrument carry out power searching, determine the leading edge position receiving signal；

Step S203, according to leading edge position, by arranging sufficiently large surplus, comprehensive test instrument intercept long enough The effective data of determination (determining the data after real rising edge)；

Step S204, seek thick frequency deviation deviant: (M is that modulation refers to first the data intercepted to be asked M power Number: for BPSK, M=2；For QPSK, M=4, for 8PSK, M=8；Behaviour by that analogy) Make to eliminate the modulation intelligence of M unit psk modulation signal；Again to the data padding obtained, carry out FFT change Change, and the data after FFT are carried out FFT shifting function so that zero-frequency is in center；So After seek the peak of absolute value of the data after FFT shifting function；Finally according to peak and centre bit The difference put and the value of FFT resolution and M obtain thick frequency deviation deviant；

Step S205, to intercept data carry out thick frequency deviation compensation；

Step S206, seek thin frequency deviation deviant: the data after compensating thick frequency deviation seek M power, and (M is for adjusting Index processed: for BPSK, M=2；For QPSK, M=4, for 8PSK, M=8；By that analogy) Thus remove modulation intelligence, then from the data obtained, intercept two segment datas, determine two intercepted segment datas Correlation, and determine thin frequency deviation deviant according to correlation；

Step S207, thick frequency deviation deviant is added with thin frequency deviation deviant, obtains final frequency deviation and estimate knot Really.

Concrete, in a wcdma system, terminal can be configured so that terminal only sends DPCCH(Dedicated Physical Control Channel, Dedicated Physical Control Channel) signal, do not send out Send DPDCH(Dedicated Physical Data Channel, Dedicated Physical Data Channel) signal, thus Obtain the QPSK(Quadrature Phase Shift Keying of standard, QPSK) modulated signal, When described signal is TD-SCDMA signal, it is also possible to terminal is configured so that terminal does not sends Control channel signal, only send traffic channel signal, and then be easy to remove modulation intelligence and carry out frequency deviation estimation.

In step S101, the psk modulation signal that receipts terminal sends, and determine the leading edge position of signal, Can be determined by the performance number of psk modulation signal sampling point sliding window, the sampling point further according to each sampled point is sliding The performance number of dynamic window determines that the leading edge position of signal realizes.

In step S102, intercept the data of at least one slot length according to leading edge position, particularly as follows:

The data of a slot length are intercepted according to leading edge position.

Concrete, intercept the data of a slot length according to leading edge position, specifically include:

Determine data intercept be r=[r (0), r (1) ..., r (N-1)], N=N_FOE, wherein r (i)=r (i+I_start+N_Δ), N_Δ For the surplus pre-set, N_FOEIt is the symbolic number of a time slot, I_startLeading edge position for signal.

In step S103, determine thick frequency deviation deviant according to the data intercepted, and inclined by this thick frequency deviation The shifting value data to being intercepted carry out thick frequency deviation compensation, specifically include:

The data intercepted carry out M power operation, and wherein, M is modulation index；

To the data padding after the M power operation obtained, carry out FFT(fast Fourier change) conversion, Wherein, FFT count in set point so that FFT resolution thin frequency deviation estimate in the range of；

Data after FFT are carried out FFT shifting function so that zero-frequency is in center；

Determine the peak of the absolute value of the data after FFT shifting function；

It is inclined that value according to peak and the difference of center and FFT resolution and M obtains thick frequency deviation Shifting value；

Thick frequency deviation compensation is carried out by these thick frequency deviation deviant data to being intercepted.

When the data intercepted according to leading edge position are more than a slot length, in basis in step S103 The data intercepted determine thick frequency deviation deviant, and are carried out by these thick frequency deviation deviant data to being intercepted Before thick frequency deviation compensates, also include:

Determine the sampling location that sampling deviation is minimum；

From the sampling location that sampling deviation is minimum, again intercept the data for carrying out FFT computing.

In step S104, the data after thick frequency deviation compensates intercept two segment datas, determines intercepted two The correlation of segment data, and determine thin frequency deviation deviant according to correlation, specifically include:

Take result r after thick frequency deviation compensates_comp=[r_comp(0),r_comp(1),…,r_comp(N-1)], N is data length；

According to data intercept original position start_pos set in advance, data length corr_length and two sections Spacing gap of data, intercepts two segment datas d_b, wherein b=1,2, gap is nonnegative integer, and

d₁=r_comp(start_pos:corr_length+start_pos-1)

d₂=r_comp(corr_length+start_pos+gap:2corr_length+gap+start_pos-1)；

By to d_bAsk biquadratic to remove modulation intelligence, obtain:

d_b,sqr=[d_b,sqr(0),d_b,sqr(1),...,d_b,sqr(corr_length-1)],

Wherein d_b,sqr(i)=(d_b(i))²·(d_b(i))², it is the number of symbols of 1 data block；

Determine d_1,sqrAnd d_2,sqrThe correlation R of middle element:

Determine that thin frequency offset estimation result is: $f_{e2}=\frac{1}{2\mathrm{\π}{T}_c(\mathrm{corr}\_\mathrm{length}+\mathrm{gap})}\arg \left(R\right)/4,$ Wherein T_cFor chip Cycle；

Below by specific embodiment, the transmitting frequency calibration method of the present invention is described in detail:

Embodiment one,

In TD-SCDMA system, terminal transmission up DPCH signal, comprehensive test instrument utilization receives DPCH signal carries out frequency deviation estimation.As it is shown on figure 3, the specifically comprising the following steps that of frequency calibration

Step S301, detection signal rising edge:

The performance number of calculating reception signal sampling point sliding window:

$P_{\mathrm{win}}\left(i\right)=\underset{n=i}{\overset{i+L_{\mathrm{win}}-1}{\mathrm{\Σ}}}{\left|r\left(n\right)\right|}^2,i=\mathrm{0,1},...,N-L_{\mathrm{win}};$

Wherein, L_winThe sampling point number comprised for power calculation window length；

Search rising edge:

Fori=0:N-L_win

if $\frac{P_{\mathrm{win}}(i+L_{\mathrm{win}})}{P_{\mathrm{win}}\left(i\right)}>P_{\lim }$

I_start=i+2L_win

end

break;

end

Wherein P_limFor rising edge difference power threshold value (linear value).

Step S302, the data of one slot length of taking-up:

R=[r (1), r (2) ..., r (864)]；

Wherein r (i)=r (i+I_start+N_Δ), N_ΔAfter guaranteeing to get the position that authentic data starts Data and the surplus that arranges；

Step S303, carry out coarse frequency offset:

R is asked 4 power operations, to eliminate the phase hit that modulation brings, obtains:

r_pow4=[r_pow4(0),r_pow4(1),...,r_pow4(N-1)], wherein r_pow4(i)=(r_FFT(i))⁴=(r_FFT(i))²(r_FFT(i))²；

To r_pow4Zero padding, obtains N_FFTPoint carries out the data of FFT, wherein N_FFTIt it is the integral number power of 2；

Carry out FFT computing, and FFT shift(FFT displacement) operation:

FFT computing obtains: r_{FFT_OUT}=FFT (r_{FFT_IN})；

FFT shift obtains: r_{FFT_Shift}=[r_{FFT_Shift}(0),r_{FFT_Shift}(1),...,r_{FFT_Shift}(N_FFT-1)]；

Wherein $r_{\mathrm{FFT}\_\mathrm{Shift}}\left(i\right)=\left\{\begin{array}{c}{r}_{\mathrm{FFT}\_\mathrm{OUT}}(i+\frac{N_{\mathrm{FFT}}}{2}),0\&le;i<\frac{N_{\mathrm{FFT}}}{2}\\ r_{\mathrm{FFT}\_\mathrm{OUT}}(i-\frac{N_{\mathrm{FFT}}}{2}),\frac{N_{\mathrm{FFT}}}{2}\&le;i<N_{\mathrm{FFT}}\end{array}\right.;$

Determine the position of maximum absolute value: $A_{\mathrm{index}}=\arg \underset{i=\mathrm{0,1},...,N_{\mathrm{FFT}}}{\max }\left(\right|r_{\mathrm{FFT}\_\mathrm{Shift}}\left(i\right)\left|\right);$

Obtain coarse frequency offset result: $f_{e1}=(A_{\mathrm{index}}-\frac{N_{\mathrm{FFT}}}{2})*\mathrm{\&Delta;f}/4;$

WhereinFrequency resolution for FFT.

Step S304, carry out thick frequency deviation compensation, obtain: r_comp=[r_comp(0),r_comp(1),…,r_comp(N-1)]；

Wherein, N=864 is the number of chips of 1 time slot,

$r_{\mathrm{comp}}\left(i\right)=r\left(i\right)\&CenterDot;e^{-j2\mathrm{\&pi;}{f}_{e1}i{T}_c},$

$r_{\mathrm{comp}}\left(i\right)=r\left(i\right)\&CenterDot;e^{-j2\mathrm{\&pi;}{f}_{e1}i{T}_c}=r_{\mathrm{comp}}\left(i\right)(\cos (2\mathrm{\&pi;}{f}_{e1}i/f_c)-j\sin (2\mathrm{\&pi;}{f}_{e1}i/f_c)),$

I=0,1 ..., N_c-1；

Step S305, carry out thin frequency deviation estimation:

Take result r after the frequency deviation compensation of a time slot_comp=[r_comp(0),r_comp(1),…,r_comp(N-1)]；

Intercepting two segment datas and do relevant, the original position of first paragraph data is start_pos, data a length of Corr_length, the spacing of two segment datas is gap, and these three parameter is the parameter that can be pre-configured with；

Determine that two data blocks carrying out correlation computations are d_b, wherein b=1,2；

d₁=r_comp(start_pos:corr_length+start_pos-1)

d₂=r_comp(corr_length+start_pos+gap:2corr_length+gap+start_pos-1)

Removal modulation intelligence:

By to d_bAsk biquadratic to remove modulation intelligence, obtain:

d_b,sqr=[d_b,sqr(0),d_b,sqr(1),...,d_b,sqr(corr_length-1)],

Wherein d_b,sqr(i)=(d_b(i))²·(d_b(i))², it is the number of symbols of 1 data block；

Ask relevant:

Calculate d_1,sqrAnd d_2,sqrThe correlation R of middle element:

Then residual frequency deviation is calculated: $f_{e2}=\frac{1}{2\mathrm{\&pi;}{T}_c(\mathrm{corr}\_\mathrm{length}+\mathrm{gap})}\arg \left(R\right)/4;$

Step S306, obtain total frequency deviation: f_e=f_e1+f_e2；

Step S307, carry out frequency calibration according to total frequency deviation.

Embodiment two,

In a wcdma system, the reception data of a radio frames are:

R=[r (0), r (1) ..., r (N-1)], N=OSR*N_c*N_slot；Wherein, OSR is over-sampling rate, N_c=2560 For the number of chips of one time slot of DPCCH channel, N_slotThe number of time slot that=15 is input data.

Now, as shown in Figure 4, the specifically comprising the following steps that of frequency calibration

Step S401, detection signal rising edge:

The performance number of calculating reception signal sampling point sliding window:

$P_{\mathrm{win}}\left(i\right)=\underset{n=i}{\overset{i+L_{\mathrm{win}}-1}{\mathrm{\&Sigma;}}}{\left|r\left(n\right)\right|}^2,i=\mathrm{0,1},...,N-L_{\mathrm{win}}$

Wherein L_winThe sampling point number comprised for power calculation window length；

Search rising edge:

Fori=0:N-L_win

if $\frac{P_{\mathrm{win}}(i+L_{\mathrm{win}})}{P_{\mathrm{win}}\left(i\right)}>P_{\lim }$

I_start=i+2L_win

end

break;

end

Wherein P_limFor rising edge difference power threshold value (linear value)；

Take out in order to the data carrying out frequency deviation estimation:

R '=[r ' (0), r ' (1) ..., r ' (N '-1)], N '=OSR*N_FOE；

Wherein r ' (i)=r (i+I_start+N_Δ* OSR), N_ΔIn order to guarantee to get position that authentic data starts it After data and the surplus that arranges, N_FOEFor carrying out effective number of chips of coarse frequency offset；

Step S402, carry out coarse frequency offset:

Ask the sampling location that sampling deviation is minimum:

Front N to r '_var* OSR data are grouped, and obtain OSR group 1 sampling data:

R ＇₀, r ＇₁..., r ＇_OSR-1；

The jth sampling point value of i-th group of data is:

r′_i(j)=r ' (j*OSR+i), i=0,1 ..., OSR-1, j=0,1 ..., N_var-1；

OSR group data are asked respectively amplitude square, obtain: r ＇_{S, 0}, r ＇_{S, 1}..., r ＇_{S, OSR-1}；

Wherein r ＇_s,i(j)=| r ＇_i(j)|², i=0,1 ..., OSR-1, j=0,1 ..., N_var-1；

The OSR group data obtained are averaged respectively, obtain:

Wherein $\stackrel{\&OverBar;}{r_{s,i}^{\&prime;}}=\frac{1}{N_{\mathrm{var}}}\underset{j=0}{\overset{N_{\mathrm{var}}-1}{\mathrm{\&Sigma;}}}{r}_{s,i}^{\&prime;}\left(j\right),$ i=0,1,…,OSR-1；

The OSR group data obtained are determined respectively the variance of mould side's vector, obtains σ₀,σ₁,...,σ_OSR-1；

Wherein ${\mathrm{\&sigma;}}_i=\frac{1}{N_{\mathrm{var}}}\underset{j=0}{\overset{N_{\mathrm{var}}-1}{\mathrm{\&Sigma;}}}{(r_{s,i}^{\&prime;}\left(j\right)-\stackrel{\&OverBar;}{r_{s,i}^{\&prime;}})}^2,$ i=0,1,…,OSR-1；

The sampling deviation of the road grouped data that variance is minimum is minimum:

FFT computing is used to carry out coarse frequency offset:

Intercept in order to the sampling point without over-sampling carrying out FFT computing:

R ＇_FFT=[r ＇_FFT(0), r ＇_FFT(1) ..., r ＇_FFT(N_FOE-1)]

Wherein, r ＇_FFT(i)=r(i*OSR+I_start+N_Δ*OSR+I_opt),i=0,1,...,N_FOE-1,N_FOECarry out FFT Effective number of samples of computing；

To r ＇_FFTAsk 4 power operations, to eliminate the phase hit that modulation brings, obtain:

R ＇_pow4=[r ＇_pow4(0), r ＇_pow4(1) ..., r ＇_pow4(N_FOE-1)];

To r ＇_pow4Zero padding obtains carrying out the data of FFT:

Wherein, N_FFTCount for FFT；

Carry out FFT computing, and FFT shift operate:

FFT computing obtains: r_{FFT_OUT}=FFT (r_{FFT_IN})；

FFT shift obtains: r_{FFT_Shift}=[r_{FFT_Shift}(0),r_{FFT_Shift}(1),...,r_{FFT_Shift}(N_FFT-1)]；

Wherein $r_{\mathrm{FFT}\_\mathrm{Shift}}\left(i\right)=\left\{\begin{array}{c}{r}_{\mathrm{FFT}\_\mathrm{OUT}}(i+\frac{N_{\mathrm{FFT}}}{2}),0\&le;i<\frac{N_{\mathrm{FFT}}}{2}\\ r_{\mathrm{FFT}\_\mathrm{OUT}}(i-\frac{N_{\mathrm{FFT}}}{2}),\frac{N_{\mathrm{FFT}}}{2}\&le;i<N_{\mathrm{FFT}}\end{array}\right..$

Determine the position of maximum absolute value: $A_{\mathrm{index}}=\arg \underset{i=\mathrm{0,1},...,N_{\mathrm{FFT}}}{\max }\left(\right|r_{\mathrm{FFT}\_\mathrm{Shift}}\left(i\right)\left|\right);$

Obtain coarse frequency offset result: $f_{e1}=(A_{\mathrm{index}}-\frac{N_{\mathrm{FFT}}}{2})*\mathrm{\&Delta;f}/4;$

WhereinFrequency resolution for FFT.

Step S403, to data r ＇_pow4Carry out thick frequency deviation compensation, obtain:

r_comp=[r_comp(0),r_comp(1),…,r_comp(N_FOE-1)],

Wherein, $r_{\mathrm{comp}}\left(i\right)=r_{\mathrm{pow}4}^{\&prime;}\left(0\right)e^{-j2\mathrm{\&pi;i}4f_{e1}{T}_c},i=\mathrm{0,1},...,N_{\mathrm{FOE}}-1;$

Step S404, carry out thin frequency deviation estimation:

To signal r_compCarry out before and after's related operation and obtain e_corr:

e_corr=[e_corr(0),e_corr(1),...,e_corr(corr_size-1)], whereincorr_size For distance relevant front and back；

Will relevant rear signal e_corrI () asks the phase place of its sum to obtain after being added:

Wherein " arg () " is for asking angle computing；

Phase place converts and obtains essence frequency offset estimation result:

Step S405, obtain total frequency deviation: f_e=f_e1+f_e2；

Step S406, carry out frequency calibration according to total frequency deviation.

The embodiment of the present invention also provides for a kind of frequency calibration device, and this device can be specially comprehensive test instrument, such as figure Shown in 5, including:

Rising edge determines unit 501, for receiving the psk modulation signal that terminal sends, and determines signal Leading edge position；

Data cutout unit 502, for intercepting sufficiently long data according to leading edge position；

Coarse frequency offset unit 503, for determining thick frequency deviation deviant according to the data intercepted, and passes through These thick frequency deviation deviant data to being intercepted carry out thick frequency deviation compensation；

Thin frequency deviation estimating unit 504, intercepts two segment datas the data after compensating from thick frequency deviation, determines The correlation of two segment datas intercepted, and determine thin frequency deviation deviant according to correlation；

Frequency calibration unit 505, for according to thick frequency deviation deviant and the sum of thin frequency deviation deviant, carries out frequency Rate is calibrated.

When signal is WCDMA signal, rising edge determines that unit 501 receives the PSK tune that terminal sends Signal processed, and before determining the leading edge position of signal, be additionally operable to:

Terminal is configured so that terminal only sends DPCCH signal, does not send DPDCH signal；

When signal is TD-SCDMA signal, rising edge determines that unit 501 receives the PSK that terminal sends Modulated signal, and before determining the leading edge position of described signal, also include:

Terminal is configured so that terminal does not send control channel signal, only sends traffic channel signal.

Wherein, rising edge determine unit 501 specifically for:

Determine the performance number of psk modulation signal sampling point sliding window；

The performance number of the sampling point sliding window according to each sampled point determines the leading edge position I of signal_start。

Data cutout unit 502 specifically for:

The data of a slot length are intercepted according to leading edge position.

Data cutout unit 502 specifically for:

Determine data intercept be r=[r (0), r (1) ..., r (N_FOE-1)], wherein r (i)=r (i+I_start+N_Δ), N_ΔFor in advance The surplus first arranged, N_FOEIt is the symbolic number of a time slot, I_startLeading edge position for signal.

Coarse frequency offset unit 503 specifically for:

The data intercepted carry out M power operation, and wherein, M is modulation index；

To the data padding after the M power operation obtained, carry out FFT, wherein, FFT Count in set point so that FFT resolution is in the range of thin frequency deviation is estimated；

Data after FFT are carried out FFT shifting function so that zero-frequency is in center；

Determine the peak of the absolute value of the data after FFT shifting function；

It is inclined that value according to peak and the difference of center and FFT resolution and M obtains thick frequency deviation Shifting value；

Thick frequency deviation compensation is carried out by these thick frequency deviation deviant data to being intercepted.

When the data intercepted according to leading edge position are more than a slot length, coarse frequency offset unit 503 Thick frequency deviation deviant is being determined according to the data intercepted, and by this thick frequency deviation deviant number to being intercepted Before carrying out thick frequency deviation compensation, it is additionally operable to:

Determine the sampling location that sampling deviation is minimum；

From the sampling location that sampling deviation is minimum, again intercept the data for carrying out FFT computing.

Thin frequency deviation estimating unit 504 specifically for:

Take result r after thick frequency deviation compensates_comp=[r_comp(0),r_comp(1),…,r_comp(N-1)], N is data length；

According to data intercept original position start_pos set in advance, data length corr_length and two sections Spacing gap of data, intercepts two segment datas d_b, wherein b=1,2, gap is nonnegative integer, and

d₁=r_comp(start_pos:corr_length+start_pos-¹)

d₂=r_comp(corr_length+start_pos+gap:2corr_length+gap+start_pos-1)；

By to d_bAsk biquadratic to remove modulation intelligence, obtain:

d_b,sqr=[d_b,sqr(0),d_b,sqr(1),...,d_b,sqr(corr_length-1)],

Wherein d_b,sqr(i)=(d_b(i))²·(d_b(i))², it is the number of symbols of 1 data block；

Determine d_1,sqrAnd d_2,sqrThe correlation R of middle element:

Determine that thin frequency offset estimation result is: $f_{e2}=\frac{1}{2\mathrm{\&pi;}{T}_c(\mathrm{corr}\_\mathrm{length}+\mathrm{gap})}\arg \left(R\right)/4,$ Wherein T_cFor chip Cycle；

The frequency calibration device that the embodiment of the present invention provides can be specially comprehensive test instrument, and this comprehensive test instrument includes can Receive the antenna of the psk modulation signal that terminal sends, be used for determining thick frequency deviation deviant and thin frequency deviation skew Value is gone forward side by side the CPU of line frequency calibration process, and for when CPU carries out respective handling to related data Carry out the memory element cached.

The transmitting frequency calibration method provided by the embodiment of the present invention, comprehensive test instrument can not know that terminal transmission is believed In the case of number cell information, the big frequency deviation of terminal is estimated accurately and effectively.This frequency deviation estimating method Need not synchronization timing accurately, as long as receiving valid data, being divided into thick frequency deviation to estimate on this basis Meter and thin frequency deviation are estimated, coarse frequency offset therein can ensure that the scope that frequency deviation is estimated, thin frequency deviation is estimated to protect The precision that card is estimated.Reduce the requirement to timing, and process is simple, eliminate comprehensive test instrument acquisition community and join It is mutual that confidence ceases, and is a kind of highly efficient transmitting frequency calibration method.

The embodiment of the present invention provides a kind of transmitting frequency calibration method and device, is carrying out coarse frequency offset and thick frequency deviation After compensation, the data after thick frequency deviation compensates intercept two segment datas, determines the phase of two intercepted segment datas Pass is worth, and can determine thin frequency deviation deviant according to correlation, without determining accurate time slot original position Carry out thin frequency deviation estimation, reduce the complexity of frequency calibration, improve testing efficiency.

Those skilled in the art are it should be appreciated that embodiments of the invention can be provided as method, system or meter Calculation machine program product.Therefore, the present invention can use complete hardware embodiment, complete software implementation or knot The form of the embodiment in terms of conjunction software and hardware.And, the present invention can use and wherein wrap one or more Computer-usable storage medium containing computer usable program code (include but not limited to disk memory, CD-ROM, optical memory etc.) form of the upper computer program implemented.

The present invention is with reference to method, equipment (system) and computer program product according to embodiments of the present invention The flow chart of product and/or block diagram describe.It should be understood that can by computer program instructions flowchart and / or block diagram in each flow process and/or flow process in square frame and flow chart and/or block diagram and/ Or the combination of square frame.These computer program instructions can be provided to general purpose computer, special-purpose computer, embedding The processor of formula datatron or other programmable data processing device is to produce a machine so that by calculating The instruction that the processor of machine or other programmable data processing device performs produces for realizing at flow chart one The device of the function specified in individual flow process or multiple flow process and/or one square frame of block diagram or multiple square frame.

These computer program instructions may be alternatively stored in and computer or the process of other programmable datas can be guided to set In the standby computer-readable memory worked in a specific way so that be stored in this computer-readable memory Instruction produce and include the manufacture of command device, this command device realizes in one flow process or multiple of flow chart The function specified in flow process and/or one square frame of block diagram or multiple square frame.

These computer program instructions also can be loaded in computer or other programmable data processing device, makes Sequence of operations step must be performed to produce computer implemented place on computer or other programmable devices Reason, thus the instruction performed on computer or other programmable devices provides for realizing flow chart one The step of the function specified in flow process or multiple flow process and/or one square frame of block diagram or multiple square frame.

Although preferred embodiments of the present invention have been described, but those skilled in the art once know base This creativeness concept, then can make other change and amendment to these embodiments.So, appended right is wanted Ask and be intended to be construed to include preferred embodiment and fall into all changes and the amendment of the scope of the invention.

Obviously, those skilled in the art can carry out various change and modification without deviating from this to the present invention Bright spirit and scope.So, if the present invention these amendment and modification belong to the claims in the present invention and Within the scope of its equivalent technologies, then the present invention is also intended to comprise these change and modification.

Claims (12)

1. a transmitting frequency calibration method, it is characterised in that including:

Terminal is configured so that terminal only sends special used for physical control channel DPCCH signal, does not sends out Sending Dedicated Physical Data Channel DPDCH signal, wherein, described signal is WCDMA WCDMA Signal；Or, terminal is configured so that terminal does not send control channel signal, only sends business letter Road signal, wherein, described signal is TD SDMA TD-SCDMA signal；

Receive the phase-shift keying (PSK) psk modulation signal that terminal sends, and determine the leading edge position of described signal；

Sufficiently long data are intercepted according to leading edge position；

Thick frequency deviation deviant is determined according to the data intercepted, and by this thick frequency deviation deviant to being intercepted Data carry out thick frequency deviation compensation；

Data after thick frequency deviation compensates intercept two segment datas, determines the relevant of two intercepted segment datas Value, and determine thin frequency deviation deviant according to described correlation；

According to described thick frequency deviation deviant and the sum of described thin frequency deviation deviant, carry out frequency calibration.

2. the method for claim 1, it is characterised in that the PSK that described reception terminal sends adjusts Signal processed, and determine the leading edge position of described signal, specifically include:

Determine the performance number of described psk modulation signal sampling point sliding window；

The performance number of the sampling point sliding window according to each sampled point determines the leading edge position I of described signal_start。

3. the method for claim 1, it is characterised in that described according to leading edge position intercept extremely The data of a few slot length, particularly as follows:

The data of a slot length are intercepted according to leading edge position.

4. method as claimed in claim 3, it is characterised in that described according to leading edge position intercepting one The data of individual slot length, specifically include:

Determine data intercept be r=[r (0), r (1) ..., r (N-1)], N=N_FOE, wherein r (i)=r (i+I_start+N_Δ), N_Δ For the surplus pre-set, N_FOEIt is the symbolic number of a time slot, I_startLeading edge position for described signal.

5. the method for claim 1, it is characterised in that determine thick frequency according to the data intercepted Deviant partially, and carry out thick frequency deviation compensation by these thick frequency deviation deviant data to being intercepted, specifically include:

The data intercepted carry out M power operation, and wherein, M is modulation index；

To the data padding after the M power operation obtained, carry out FFT, wherein, FFT Count in set point so that FFT resolution is in the range of thin frequency deviation is estimated；

Data after FFT are carried out FFT shifting function so that zero-frequency is in center；

Determine the peak of the absolute value of the data after FFT shifting function；

It is inclined that value according to peak and the difference of center and FFT resolution and M obtains thick frequency deviation Shifting value；

Thick frequency deviation compensation is carried out by these thick frequency deviation deviant data to being intercepted.

6. the method for claim 1, it is characterised in that described data after thick frequency deviation compensates Middle intercepting two segment data, determines the correlation of two intercepted segment datas, and determines carefully according to described correlation Frequency deviation deviant, specifically includes:

Take result r after thick frequency deviation compensates_comp=[r_comp(0),r_comp(1),…,r_comp(N-1)], N is data length；

According to data intercept original position start_pos set in advance, data length corr_length and two sections Spacing gap of data, intercepts two segment datas d_b, wherein b=1,2, gap is nonnegative integer, and

d₁=r_comp(start_pos:corr_length+start_pos-1)

d₂=r_comp(corr_length+start_pos+gap:2corr_length+gap+start_pos-1)；

By to d_bAsk biquadratic to remove modulation intelligence, obtain:

d_b,sqr=[d_b,sqr(0),d_b,sqr(1),...,d_b,sqr(corr_length-1)],

Wherein d_b,sqr(i)=(d_b(i))²·(d_b(i))², it is the number of symbols of 1 data block；

Determine d_1,sqrAnd d_2,sqrThe correlation R of middle element:

Determine that thin frequency offset estimation result is:Wherein T_cFor chip Cycle.

7. a frequency calibration device, it is characterised in that including:

Rising edge determines unit, for configuring terminal so that terminal only sends special physical control letter Road DPCCH signal, does not send Dedicated Physical Data Channel DPDCH signal, and wherein, described signal is WCDMA WCDMA signal；Or, terminal is configured so that terminal does not send control Channel signal, only sends traffic channel signal, and wherein, described signal is TD SDMA TD-SCDMA signal；And, for receiving the psk modulation signal that terminal sends, and determine described letter Number leading edge position；

Data cutout unit, for intercepting sufficiently long data according to leading edge position；

Coarse frequency offset unit is for determining thick frequency deviation deviant according to the data intercepted and thick by this The frequency deviation deviant data to being intercepted carry out thick frequency deviation compensation；

Thin frequency deviation estimating unit, intercepts two segment datas the data after compensating from thick frequency deviation, determines and cut The correlation of two segment datas taken, and determine thin frequency deviation deviant according to described correlation；

Frequency calibration unit, for according to described thick frequency deviation deviant and the sum of described thin frequency deviation deviant, enters Line frequency is calibrated.

8. device as claimed in claim 7, it is characterised in that described rising edge determines that unit is specifically used In:

Determine the performance number of described psk modulation signal sampling point sliding window；

The performance number of the sampling point sliding window according to each sampled point determines the leading edge position I of described signal_start。

9. device as claimed in claim 7, it is characterised in that described data cutout unit specifically for:

The data of a slot length are intercepted according to leading edge position.

10. device as claimed in claim 9, it is characterised in that described data cutout unit specifically for:

Determine data intercept be r=[r (0), r (1) ..., r (N_FOE-1)], wherein r (i)=r (i+I_start+N_Δ), N_ΔFor The surplus pre-set, N_FOEIt is the symbolic number of a time slot, I_startLeading edge position for described signal.

11. devices as claimed in claim 7, it is characterised in that described coarse frequency offset unit is specifically used In:

The data intercepted carry out M power operation, and wherein, M is modulation index；

To the data padding after the M power operation obtained, carry out FFT, wherein, FFT Count in set point so that FFT resolution is in the range of thin frequency deviation is estimated；

Data after FFT are carried out FFT shifting function so that zero-frequency is in center；

Determine the peak of the absolute value of the data after FFT shifting function；

It is inclined that value according to peak and the difference of center and FFT resolution and M obtains thick frequency deviation Shifting value；

Thick frequency deviation compensation is carried out by these thick frequency deviation deviant data to being intercepted.

12. devices as claimed in claim 7, it is characterised in that described thin frequency deviation estimating unit is specifically used In:

Take result r after thick frequency deviation compensates_comp=[r_comp(0),r_comp(1),…,r_comp(N-1)], N is data length；

According to data intercept original position start_pos set in advance, data length corr_length and two sections Spacing gap of data, intercepts two segment datas d_b, wherein b=1,2, gap is nonnegative integer, and

d₁=r_comp(start_pos:corr_length+start_pos-1)

d₂=r_comp(corr_length+start_pos+gap:2corr_length+gap+start_pos-1)；

By to d_bAsk biquadratic to remove modulation intelligence, obtain:

d_b,sqr=[d_b,sqr(0),d_b,sqr(1),...,d_b,sqr(corr_length-1)],

Wherein d_b,sqr(i)=(d_b(i))²·(d_b(i))², it is the number of symbols of 1 data block；

Determine d_1,sqrAnd d_2,sqrThe correlation R of middle element:

Determine that thin frequency offset estimation result is:Wherein T_cFor chip Cycle.