CN101312378A - Estimation method for time-delay spreading parameter of receiving machine and time-delay spreading parameter estimation unit - Google Patents

Abstract

The invention discloses an estimation method and a delay spreading parameter estimation unit for the delay spreading parameter of a receiver. The method comprises cross-correlating a received signal and a training sequence which are deviation-compensated by a carrier frequency to obtain a correlation sequence, and through the judgment of the threshold, finding out a location index value I<first> of one foremost path and a location index value I<last> of one rearmost path in the correlation sequence, thereby obtaining the delay spreading parameter, a remainder of subtracting the I<first> from the I<last>. The delay spreading parameter estimation unit comprises a correlation sequence generating module, a threshold judging module and a delay spreading parameter computing module, wherein the correlation sequence generating module is used to cross-correlate the received signal and the training sequence which are deviation-compensated by the carrier frequency to obtain the correlation sequence, the threshold judging module is used to find out the location index value I<first> of one foremost path and the location index value I<last> of one rearmost path in the correlation sequence through the judgment of the threshold, and the delay spreading parameter computing module is used to receive the I<first> and the I<last>, and to obtain the delay spreading parameter by subtracting the I<first> from the I<last>. The estimation method and the delay spreading parameter estimation unit can enable the receiver to obtain comparatively accurate delay spreading.

Description

The method of estimation of receiver delay spreading parameter and delay spreading parameter estimation unit

Technical field

The present invention relates to the communications field, particularly a kind of method of estimation of receiver delay spreading parameter and delay spreading parameter estimation unit.

Background technology

In communication system,, just need carry out rational channel estimating for the data of being transmitted in the correct demodulation of receiving terminal.At present, be used to estimate that the method for wireless channel has two classes greatly, a class is based on training sequence, and a class is based on pilot tone.

Be the channel estimating of example explanation with DTMB (Digital Terrestrial/Television Multimedia Broadcasting, digital TV ground multimedia broadcasting) based on training sequence.

At first specifically introduce frame structure and frame head sequence in the DTMB system, the frame head PN sequence as training sequence in the DTMB standard has three kinds: PN420, PN595 and PN945.

Fig. 1 has provided four layers of frame structure based on multi-frame of DTMB system.Signal frame is the most basic transmission unit, comprises frame head and frame.Frame head is filled the PN sequence, and frame can be the data of single carrier mode, also can be the data of multi-carrier mode.Under multi-carrier mode, the size of IFFT (Inverse Fast FourierTransform, fast adverse Fourier transform) is 3780 points.The length of frame head PN sequence has three kinds: 420 data symbols, 595 data symbols and 945 data symbols, thus there are three kinds of signal frame length accordingly.The baseband sampling rate of DTMB system is 7.56MSPS (Mega-Samples-Per-Second, per second sample 1,000,000 times).Therefore, the time span of Dui Ying three kinds of PN frame head sequences is respectively 420/7.56=55.56 microsecond, 595/7.56=78.703 microsecond and 945/7.56=125 microsecond.And the time span of frame is 500 microseconds.Thereby the time span of each signal frame is respectively 555.56 microseconds, 578.703 microseconds and 625 microseconds accordingly.Be superframe on signal frame, the time span of each superframe is unified to be 125ms.Accordingly, a superframe comprises 225,216 and 200 signal frames respectively.For dividing frame, time span is one minute on the superframe.A branch frame comprises 480 superframes.Up be a day frame again, corresponding one day 24 hours.One day, frame comprised 1440 branch frames.

The building method of three kinds of PN frame head sequences is also inconsistent.Fig. 2 has provided the structure chart of frame head PN420 sequence.The centre is the PN255 sequence, is the m sequence on one 8 rank, and reaching " 1 " through " 0 " to+1 value is the binary character of non-return-to-zero to the mapping transformation of-1 value.The PN420 sequence comprises synchronizing sequence (length is 83 symbols) behind a preamble sequence (length is 82 symbols), PN255 sequence and.Preamble sequence is the Cyclic Prefix sequence of PN255 sequence, and back synchronizing sequence is the cyclic suffix sequence of PN255 sequence.The PN255 sequence adopts LFSR (Linear Feedback Shift Register, linear feedback shift register) to generate, and the initial phase of different LFSR will generate different PN255 sequences.

Fig. 3 has provided the structure of frame head PN945 sequence.The structure of PN945 and the similar of PN420 all have preamble sequence and back synchronizing sequence.The intermediate sequence of PN945 sequence is a PN511 sequence, adopts the m sequence on one 9 rank to generate, and the mapping transformation to-1 value is the binary symbol sequence of non-return-to-zero to+1 value and " 1 " through " 0 " again.The concrete value of PN945 sequence also is variable, and is relevant with the initial phase of LFSR.And in order to reduce adjacent PN420 sequence or the correlation between the PN945 sequence, the DTMB standard has been arranged initial phase once meticulously by Computer Simulation specially.

The building method of frame head PN420 and PN945 is the same, but the building method of the building method of frame head PN595 and PN420 and PN945 differs widely.The PN595 sequence is to adopt 10 rank maximum length pseudo-random binary sequence brachymemmas to form, and is that length is preceding 595 chips of 1023 m sequence.And the initial phase of LFSR that generates the m sequence of this 1023 length is also fixed, and is 0000000001, and promptly the PN595 sequence that adopts of each signal frame is identical.Preceding 595 chips of pseudo random sequence, the mapping transformation to-1 value is the binary symbol sequence of non-return-to-zero to+1 value and " 1 " through " 0 ", is PN595 frame head sequence.

Receiver is after received signal, earlier the supposition known transmitter is launched which kind of frame head sequence one of (three kinds frame head sequence PN420, PN595 and PN945), on receiver, utilize frame head sequence that is stored in local correspondence in advance and the signal that receives to carry out cross-correlation operation then, seek the output result's of cross-correlation operation amplitude peak value, and then the original position of PN frame head in definite received signal, know the frame head PN pattern (one of PN420, PN595 or PN945) and the frame head PN sequence of transmitter emission, and estimate carrier frequency offset; After this on these results' basis, also need to carry out more meticulous channel estimating.

Time delay expansion is estimated follow-up receiver channel and the training sequence elimination all can be influential; Especially for DTMB etc. because channel quality is relatively poor relatively, the time delay expansion of multipath channel can relatively large system, time delay expansion can exert an influence to estimated result especially.

Summary of the invention

The technical problem to be solved in the present invention provides a kind of method of estimation and delay spreading parameter estimation unit of receiver delay spreading parameter, be applicable to that receiver can determine the communication system of training sequence, can carry out the estimation of delay spreading parameter, thereby finish more precise channels estimation.

In order to address the above problem, the invention provides a kind of method of estimation of receiver delay spreading parameter, comprising:

(A) received signal and training sequence through the carrier frequency offset compensation are carried out cross-correlation operation, obtain correlated series;

(B) judge by thresholding, find out in this correlated series the location index value I in last footpath _FirstLocation index value I with last footpath _Last

(C) delay spreading parameter is I _LastDeduct I _FirstAfter poor.

Further, steps A specifically comprises:

Intercepting one piece of data sequence is as pending data sequence r from the received signal that compensates through carrier frequency offset _Pn(n); With training sequence as secondary processing sequence s _Pn(m); Pending data sequence and secondary processing sequence are carried out cross-correlation operation, obtain correlated series.

Further, in the steps A:

Pending data sequence r _Pn(n) length is 2 * L _{Ds, max}+ P, n are the index of sample of signal, and the value of n is from I _Sync-L _{Ds, max}To I _Sync+ L _{Ds, max}Each integer of+P-1 comprises I _Sync-L _{Ds, max}And I _Sync+ L _{Ds, max}+ P-1;

Wherein P is the length of training sequence, and frame synchronization position index is I _Sync, the requirement of the maximum delay expansion of receiver is L _{Ds, max}

Further, when this method is applied to the DTMB system, in the steps A training sequence is meant as secondary processing sequence: choose local corresponding PN sequence of being preserved as secondary processing sequence s according to frame head mode that estimates and PN sequence phase _Pn(m), this sequence length is P, and wherein the value of index m is each integer from 0 to P-1, comprises 0 and P-1.

Further, when this method is applied to the DTMB system, when frame head PN pattern is PN420, P=420; When frame head PN pattern is PN595, P=595; When frame head PN pattern is PN945, P=945.

Further, to sequence r _Pn(n) and sequence s _Pn(m) do cross-correlation operation according to following formula, obtain correlated series s _Corr(l):

$s_{\mathrm{corr}}\left(l\right)=\underset{m=0}{\overset{P-1}{\mathrm{\&Sigma;}}}\mathrm{conj}\left(s_{\mathrm{pn}}\left(m\right)\right)\&CenterDot;r_{\mathrm{pn}}(I_{\mathrm{sync}}-L_{\mathrm{ds},\max }+l+m)$

Wherein conj () expression complex conjugate operation, the value of l is to 2 * L from 0 _{Ds, max}Each integer, comprise 0 and 2 * L _{Ds, max}

Further, step B specifically comprises:

(B1) the power factor threshold T is set _Power

(B2) ask the power sequence of correlated series;

(B3) ask the average power P of power sequence _{Corr, avg}

(B4) power sequence P is judged on sample-by-sample ground from front to back _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}Wherein, greater than T _Power* P _{Corr, avg}The index value of first sample be I _First, greater than T _Powet* P _{Corr, avg}The index value of last sample be I _Last

Further, when this method is applied to the DTMB system, among the step B1, T _PowerBe 16.

Further, step B4 specifically comprises:

(91) S set of structure, sample-by-sample ground judges that whether each sample in the power sequence is greater than T from front to back _Power* P _{Corr, avg}If greater than, then the index value of this sample is the element of S set just; Otherwise, be not the element of S set;

(92) if S set for null set, then first element value of S set is the location index value in last footpath, is designated as I _FirstLast element value of S set is the location index value in last footpath, is designated as I _Last

The present invention also provides a kind of receiver delay spreading parameter estimation unit, comprising:

The correlated series generation module is used to receive the received signal through the carrier frequency offset compensation, and itself and training sequence are carried out cross-correlation operation, obtains correlated series;

The thresholding judge module is used for judging by thresholding, finds out in the described correlated series the location index value I in last footpath _FirstLocation index value I with last footpath _Last

The delay spreading parameter computing module is used to receive I _FirstAnd I _Last, and use I _LastDeduct I _FirstObtain delay spreading parameter L _Ds

Further, described correlated series generation module comprises:

Pending sequence generates submodule, is used to receive the received signal through the carrier frequency offset compensation, and therefrom intercepts the one piece of data sequence as pending data sequence r _Pn(n);

The secondary processing sequence generation module is used to provide training sequence as secondary processing sequence s _Pn(m);

The cross-correlation operation submodule is used for pending data sequence r _Pn(n) and secondary processing sequence s _Pn(m) carry out cross-correlation operation, obtain correlated series s _Corr(l).

Further, described pending sequence generates the pending data sequence r that submodule generates _Pn(n) length is 2 * L _{Ds, max}+ P, wherein n is the index of sample of signal, the value of n is from I _Sync-L _{Ds, max}To I _Sync+ L _{Ds, max}Each integer of+P-1 comprises I _Sync-L _{Ds, max}And I _Sync+ L _{Ds, max}+ P-1;

Wherein P is the length of frame head PN sequence, and frame synchronization position index is I _Sync, the requirement of the maximum delay expansion of receiver is L _{Ds, max}

Further, when this estimation unit was applied to the DTMB system, described secondary processing sequence generates submodule provided training sequence to be meant as secondary processing sequence: secondary processing sequence generates submodule and chooses local corresponding PN sequence of being preserved as secondary processing sequence s according to frame head mode that estimates and PN sequence phase _Pn(m), this sequence length is P, and wherein the value of index m is each integer from 0 to P-1, comprises 0 and P-1.

Further, when this estimation unit is applied to the DTMB system, when frame head PN pattern is PN420, P=420; When frame head PN pattern is PN595, P=595; When frame head PN pattern is PN945, P=945.

Further, described cross-correlation operation submodule carries out cross-correlation operation to two sequences according to following formula, obtains correlated series s _Corr(l):

$s_{\mathrm{corr}}\left(l\right)=\underset{m=0}{\overset{P-1}{\mathrm{\&Sigma;}}}\mathrm{conj}\left(s_{\mathrm{pn}}\left(m\right)\right)\&CenterDot;r_{\mathrm{pn}}(I_{\mathrm{sync}}-L_{\mathrm{ds},\max }+l+m)$

Wherein conj () expression complex conjugate operation, the value of l is to 2 * L from 0 _{Ds, max}Each integer, comprise 0 and 2 * L _{Ds, max}

Further, described thresholding judge module comprises:

The threshold value sub module stored is used to store set power factor threshold T _Power

Power sequence generates submodule, is used to ask the power sequence of correlated series;

Average power generates submodule, is used to ask the average power P of power sequence _{Orr, avg}

Index value is judged submodule, is used for from front to back sample-by-sample ground and judges that whether each sample of power sequence is greater than T _Power* P _{Corr, avg}, and will be greater than T _Power* P _{Orr, avg}The index value of first sample be designated as I _First, will be greater than T _Power* P _{Corr, avg}The index value of last sample be designated as I _Last

Further, when this method was applied to the DTMB system, described threshold value was provided with the set power factor threshold T of submodule _PowerBe 16.

Further, described index value is judged submodule sample-by-sample ground judgement from front to back power sequence P _CorrWhether each sample (l) is greater than T _Powet* P _{Corr, avg}, and will be greater than T _Powet* P _{Corr, avg}The index value of first sample be designated as I _First, will be greater than T _Power* P _{Corr, avg}The index value of last sample be designated as I _LastSpecifically be meant:

Index value is judged S set of submodule structure, and power sequence P is judged on sample-by-sample ground from front to back _CorrWhether each sample (l) is greater than T _Powet* P _{Corr, avg}If greater than, the index value of this sample is designated as the element of S set; Otherwise, be not designated as the element of S set; If S set is not null set, then index value judges that submodule is designated as I with first element value of S set _First, this is the location index value in last footpath, and last element value is designated as I _Last, this is the location index value in last footpath.

The present invention is directed to receiver and can determine the communication system of training sequence, proposed a kind of receiver and carried out the scheme that delay spreading parameter is estimated, can make receiver obtain comparatively accurate time delay expansion.Further, the present invention has provided the processing procedure of detailed cross-correlation operation process, thresholding judgement, and the processing details of estimation scheme has been carried out the refinement description.Further, the present invention has also provided concrete parameter value at the similar system of DTMB system and other and DTMB.In addition, the present invention gives preferred power factor threshold value, makes the reliability of estimation scheme stronger.

Description of drawings

Fig. 1 is based on the schematic diagram of four layers of frame structure of multi-frame in the existing DTMB system;

Fig. 2 is the structural representation of the frame head PN420 sequence in the existing DTMB system;

Fig. 3 is the structural representation of the frame head PN945 sequence in the existing DTMB system;

Fig. 4 is the concrete implementing procedure figure of the method for estimation of receiver delay spreading parameter of the present invention;

Fig. 5 is in the method for estimation of receiver delay spreading parameter of the present invention, carries out the concrete implementing procedure figure that thresholding is judged;

Fig. 6 is the concrete enforcement block diagram of the estimation unit of receiver delay spreading parameter of the present invention;

Fig. 7 is in the estimation unit of receiver delay spreading parameter of the present invention, the concrete enforcement block diagram of correlated series generation module;

Fig. 8 is in the estimation unit of receiver delay spreading parameter of the present invention, the concrete enforcement block diagram of thresholding judge module.

Embodiment

Below in conjunction with drawings and Examples technical scheme of the present invention is described in detail.

When receiver obtains the frame synchronization position, knows training sequence, and estimate carrier frequency offset, carried out carrying out the estimation of delay spreading parameter after the compensation of carrier frequency offset to received signal; The length of note training sequence is P baseband sample point; Frame synchronization position (this frame synchronization position has identified the starting point of the frame head PN sequence) index that the Synchronous Processing step obtains before the note is I _SyncThe requirement that the note system docking is received the maximum delay expansion of machine is L _{Ds, max}Remember under the current multipath channel that the delay spreading parameter that will estimate is L _Ds

The method of estimation of the receiver delay spreading parameter that the present invention proposes as shown in Figure 4, comprising:

A, received signal and training sequence that the process carrier frequency offset is compensated carry out cross-correlation operation, obtain correlated series;

B, judge, find out in this correlated series the location index value I in last footpath by thresholding _FirstLocation index value I with last footpath _Last

C, delay spreading parameter are I _LastDeduct I _FirstAfter poor, i.e. L _Ds=I _Last-I _First

Wherein, steps A specifically comprises:

Intercepting one piece of data sequence is as pending data sequence r from the received signal that compensates through carrier frequency offset _Pn(n); With training sequence as secondary processing sequence s _Pn(m); Pending data sequence and secondary processing sequence are carried out cross-correlation operation, obtain correlated series.

Wherein, in the steps A, training sequence can but be not limited to frame head PN sequence, also can be other sequence, as complete 0 or complete 1 sequence etc.

Wherein, in the steps A, the pending data sequence r that is intercepted _Pn(n) length is 2 * L _{Ds, max}+ P, wherein n is the index of sample of signal, the value of n is from I _Sync-L _{Ds, max}To I _Sync+ L _{Ds, max}Each integer of+P-1 comprises I _Sync-L _{Ds, max}And I _Sync+ L _{Ds, max}+ P-1; That is:

n＝I _sync-L _ds，max，I _sync-L _ds，max+1，......，I _sync+L _ds，max+P-1。

Secondary processing sequence length is P as can be known, and wherein the value of index m is each integer from 0 to P-1, comprises 0 and P-1, that is: m=0, and 1 ..., P-1.

When this method is applied to similar system of DTMB system or other and DTMB, when frame head PN pattern is PN420, P=420; When frame head PN pattern is PN595, P=595; When frame head PN pattern is PN945, P=945.

When this method is applied to similar system of DTMB system or other and DTMB, in the steps A training sequence is meant as secondary processing sequence: choose local corresponding PN sequence of being preserved as secondary processing sequence s according to frame head mode that estimates and PN sequence phase _Pn(m).

Wherein, in the steps A to sequence r _Pn(n) and sequence s _Pn(m) do cross-correlation operation, obtain correlated series s _Corr(l), the sample point s of correlated series _Corr(l) obtain according to following formula:

$s_{\mathrm{corr}}\left(l\right)=\underset{m=0}{\overset{P-1}{\mathrm{\&Sigma;}}}\mathrm{conj}\left(s_{\mathrm{pn}}\left(m\right)\right)\&CenterDot;r_{\mathrm{pn}}(I_{\mathrm{sync}}-L_{\mathrm{ds},\max }+l+m)$

Wherein conj () expression complex conjugate operation, the value of l is to 2 * L from 0 _{Ds, max}Each integer, comprise 0 and 2 * L _{Ds, max}, that is: l=0,1 ..., 2 * L _{Ds, max}

Wherein, step B specifically can comprise as shown in Figure 5:

B1, the power factor threshold T is set _Power, this threshold value obtains according to emulation;

B2, ask correlated series S _Corr(l) power sequence P _Corr(l)=| S _Corr(l) | ², l=0,1 ..., 2L _{Ds, max}

B3, ask power sequence P _Corr(l) average power P _{Corr, avg}

B4, the power sequence of sample-by-sample ground judgement from front to back P _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}Wherein, greater than T _Power* P _{Corr, avg}The index value l of first sample be I _First, greater than T _Power* P _{Crr, avg}The index value l of last sample be I _Last

When practical application, thresholding is judged the judgement that is not limited to power, also can be the judgement to range value, such as the amplitude factor threshold T can be set _Amp, wherein can get $T_{\mathrm{amp}}=\sqrt{T_{\mathrm{power}}},$ Other processing copies B2 to carry out to B4.

When this method is applied to similar system of DTMB system or other and DTMB, in the step B1, T _PowerCan but be not limited to 16.

Wherein, in the step B3, specifically ask average power P _{Corr, avg}Method can but to be not limited to be P _{Corr, avg}=mean (P _Corr(l)), wherein mean () represents to average operation.

Wherein, step B4 can specifically comprise:

(1) S set of structure, the element of this set is the relative indexing value of every path position in the multipath channel; The method that obtains the element of this set is: power sequence P is judged on sample-by-sample ground from front to back _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}If greater than, then the index value l of this sample is the element of S set just; Otherwise, be not the element of S set;

(2) if S set for null set, then first element value of S set is the location index value in last footpath, is designated as I _FirstLast element value of S set is the location index value in last footpath, is designated as I _Last

This method can also comprise: if thresholding is judged failure (is null set such as S set), represent that then the delay spreading parameter estimation is unsuccessful, adjust T _PowerAnd/or L _{Ds, max}(such as turning T down _Power, transfer big L _{Ds, max}), if only adjusted T _Power, then return above-mentioned steps B, if adjusted L _{Ds, max}, return above-mentioned steps A after then generating new pending sequence, promptly carry out the estimation of delay spreading parameter once more, until success.

The receiver delay spreading parameter estimation unit that the present invention proposes as shown in Figure 6, comprising: correlated series generation module, thresholding judge module and delay spreading parameter computing module; Wherein:

Described correlated series generation module is used to receive the received signal through the carrier frequency offset compensation, and itself and training sequence are carried out cross-correlation operation, obtains correlated series and sends to the thresholding judge module;

Described thresholding judge module is used for judging by thresholding, finds out in the described correlated series the location index value I in last footpath _FirstLocation index value I with last footpath _Last, and send to the delay spreading parameter computing module;

The delay spreading parameter computing module is used to receive I _FirstAnd I _Last, and use I _LastDeduct I _FirstObtain delay spreading parameter L _Ds

Wherein, described correlated series generation module comprises that pending sequence generates submodule, secondary processing sequence generates submodule and cross-correlation operation submodule as shown in Figure 7;

Described pending sequence generates submodule and is used to receive the received signal that compensates through carrier frequency offset, and therefrom intercepts the one piece of data sequence as pending data sequence r _Pn(n), send to the cross-correlation operation submodule;

Described secondary processing sequence generates submodule and is used to provide training sequence as secondary processing sequence s _PnAnd send to the cross-correlation operation submodule (m);

Described cross-correlation operation submodule is used for pending data sequence r _Pn(n) and secondary processing sequence s _Pn(m) carry out cross-correlation operation, obtain correlated series S _Corr(l).

Wherein, described secondary processing sequence generate training sequence that submodule provided can but be not limited to frame head PN sequence, also can be other sequence, as complete 0 or complete 1 sequence etc.

Wherein, described pending sequence generates the pending data sequence r that submodule generates _Pn(n) length is 2 * L _{Ds, max}+ P, wherein n is the index of sample of signal, the value of n is from I _Sync-L _{Ds, max}To I _Sync+ L _{Ds, max}Each integer of+P-1 comprises I _Sync-L _{Ds, max}And I _Sync+ L _{Ds, max}+ P-1; That is:

n＝I _sync-L _ds，max，I _sync-L _ds，max+1，......，I _sync+L _ds，max+P-1。

Wherein, described secondary processing sequence generates the secondary processing sequence s that submodule sends _Pn(m) length is P, and wherein the value of index m is each integer from 0 to P-1, comprises 0 and P-1, that is: m=0, and 1 ..., P-1.

When this estimation unit is applied to similar system of DTMB system or other and DTMB, when frame head PN pattern is PN420, P=420; When frame head PN pattern is PN595, P=595; When frame head PN pattern is PN945, P=945.

When this estimation unit was applied to similar system of DTMB system or other and DTMB, described secondary processing sequence generates submodule provided training sequence to be meant as secondary processing sequence: secondary processing sequence generates submodule and chooses local corresponding PN sequence of being preserved as secondary processing sequence s according to frame head mode that estimates and PN sequence phase _Pn(m).

Described cross-correlation operation submodule carries out cross-correlation operation to two sequences according to following formula, obtains correlated series s _Corr(l):

$s_{\mathrm{corr}}\left(l\right)=\underset{m=0}{\overset{P-1}{\mathrm{\&Sigma;}}}\mathrm{conj}\left(s_{\mathrm{pn}}\left(m\right)\right)\&CenterDot;r_{\mathrm{pn}}(I_{\mathrm{sync}}-L_{\mathrm{ds},\max }+l+m)$

Wherein conj () expression complex conjugate operation, the value of l is to 2 * L from 0 _{Ds, max}Each integer, comprise 0 and 2 * L _{Ds, max}, that is: l=0,1 ..., 2 * L _{Ds, max}

Wherein, described thresholding judge module can comprise that threshold value is provided with submodule, power sequence generates submodule, average power generation submodule and index value and judges submodule as shown in Figure 8;

Described threshold value sub module stored is used to store set power factor threshold T _Power, this threshold value obtains according to emulation;

Described power sequence generates submodule and is used to ask correlated series S _Corr(l) power sequence P _Corr(l)=| S _Corr(l) | ², l=0,1 ..., 2L _{Ds, max}

Described average power generates submodule and is used to ask power sequence P _Corr(l) average power P _{Corr, avg}

Described index value judges that submodule is used for the power sequence of sample-by-sample ground judgement from front to back P _Corr(l) each sample in carries out P _Corr(l) whether greater than T _Power* P _{Corr, avg}, and will be greater than T _Power* P _{Corr, avg}The index value l of first sample be designated as I _First, will be greater than T _Power* P _{Corr, avg}The index value l of last sample be designated as I _Last

When practical application, the thresholding judge module is not limited to the judgement to power, also can be the judgement to range value, such as threshold value submodule is set and can be used to be provided with the amplitude factor threshold T _Amp, wherein can get $T_{\mathrm{amp}}=\sqrt{T_{\mathrm{power}}},$ Other each submodule is copied above-mentioned submodule design.

When this method was applied to similar system of DTMB system or other and DTMB, described threshold value was provided with the set power factor threshold T of submodule _PowerCan but be not limited to 16.

Wherein, described average power generate submodule can but be not limited to ask average power P by following formula _{Corr, avg}: P _{Corr, avg}=mean (P _Corr(l)), wherein mean () represents to average operation.

Wherein, described index value is judged submodule sample-by-sample ground judgement from front to back power sequence P _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}, and will be greater than T _Power* P _{Corr, avg}The index value l of first sample be designated as I _First, will be greater than T _Power* P _{Corr, avg}The index value l of last sample be designated as I _LastSpecifically can be meant:

Index value is judged S set of submodule structure, and the element of this set is the relative indexing value of every path position in the multipath channel; Index value is judged submodule sample-by-sample ground judgement from front to back power sequence P _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}If greater than, this index value l is designated as the element of S set; Otherwise, be not designated as the element of S set; If S set is not null set, then index value judges that submodule is designated as I with first element value of S set _First, this is the location index value in last footpath, and last element value is designated as I _Last, this is the location index value in last footpath.

Described thresholding judge module can also be used for: if thresholding is judged failure, then T is adjusted in indication _PowerAnd/or L _{Ds, max}If only adjusted T _Power, then the thresholding judge module rejudges; If adjusted L _{Ds, max}, then indicate pending sequence to generate submodule and generate new pending sequence; Accordingly, generate new pending sequence after, cross-correlation operation submodule and thresholding judge module carry out respective handling according to new pending sequence successively, if also adjusted simultaneously T _Power, the new T of foundation when then the thresholding judge module is judged _Power

Certainly; the present invention also can have other various embodiments; under the situation that does not deviate from spirit of the present invention and essence thereof; those of ordinary skill in the art work as can make various corresponding changes and distortion according to the present invention, but these corresponding changes and distortion all should belong to the protection range of claim of the present invention.

Claims (18)

1, a kind of method of estimation of receiver delay spreading parameter comprises:

(A) received signal and training sequence through the carrier frequency offset compensation are carried out cross-correlation operation, obtain correlated series;

(B) judge by thresholding, find out in this correlated series the location index value I in last footpath _FirstLocation index value I with last footpath _Last

(C) delay spreading parameter is I _LastDeduct I _FirstAfter poor.

2, method of estimation as claimed in claim 1 is characterized in that, steps A specifically comprises:

Intercepting one piece of data sequence is as pending data sequence r from the received signal that compensates through carrier frequency offset _Pn(n); With training sequence as secondary processing sequence s _Pn(m); Pending data sequence and secondary processing sequence are carried out cross-correlation operation, obtain correlated series.

3, method of estimation as claimed in claim 2 is characterized in that, in the steps A:

Pending data sequence r _Pn(n) length is 2 * L _{Ds, max}+ P, n are the index of sample of signal, and the value of n is from I _Sync-L _{Ds, max}To I _Sync+ L _{Ds, max}Each integer of+P-1 comprises I _Sync-L _{Ds, max}And I _Sync+ L _{Ds, max}+ P-1;

Wherein P is the length of training sequence, and frame synchronization position index is I _Sync, the requirement of the maximum delay expansion of receiver is L _{Ds, max}

4, method of estimation as claimed in claim 2 is characterized in that:

When this method is applied to the DTMB system, in the steps A training sequence is meant as secondary processing sequence: choose local corresponding PN sequence of being preserved as secondary processing sequence s according to frame head mode that estimates and PN sequence phase _Pn(m), this sequence length is P, and wherein the value of index m is each integer from 0 to P-1, comprises 0 and P-1.

5, as claim 3 or 4 described methods of estimation, it is characterized in that:

When this method is applied to the DTMB system, when frame head PN pattern is PN420, P=420; When frame head PN pattern is PN595, P=595; When frame head PN pattern is PN945, P=945.

6, method of estimation as claimed in claim 3 is characterized in that, to sequence r _Pn(n) and sequence s _Pn(m) do cross-correlation operation according to following formula, obtain correlated series s _Corr(l):

$s_{\mathrm{corr}}\left(l\right)=\underset{m=0}{\overset{P-1}{\mathrm{\&Sigma;}}}\mathrm{conj}\left(s_{\mathrm{pn}}\left(m\right)\right)\&CenterDot;r_{\mathrm{pn}}(I_{\mathrm{sync}}-L_{\mathrm{ds},\max }+l+m)$

Wherein conj () expression complex conjugate operation, the value of l is to 2 * L from 0 _{Ds, max}Each integer, comprise 0 and 2 * L _{Ds, max}

7, as each described method of estimation in the claim 1 to 4,6, it is characterized in that step B specifically comprises:

(B1) the power factor threshold T is set _Power

(B2) ask the power sequence of correlated series;

(B3) ask the average power P of power sequence _{Corr, avg}

(B4) power sequence P is judged on sample-by-sample ground from front to back _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}Wherein, greater than T _Power* P _{Corr, avg}The index value of first sample be I _First, greater than T _Power* P _{Corr, avg}The index value of last sample be I _Last

8, method of estimation as claimed in claim 7 is characterized in that:

When this method is applied to the DTMB system, among the step B1, T _PowerBe 16.

9, method of estimation as claimed in claim 7 is characterized in that, step B4 specifically comprises:

(91) S set of structure, sample-by-sample ground judges that whether each sample in the power sequence is greater than T from front to back _Power* P _{Corr, avg}If greater than, then the index value of this sample is the element of S set just; Otherwise, be not the element of S set;

(92) if S set for null set, then first element value of S set is the location index value in last footpath, is designated as I _FirstLast element value of S set is the location index value in last footpath, is designated as I _Last

10, a kind of receiver delay spreading parameter estimation unit is characterized in that, comprising:

The correlated series generation module is used to receive the received signal through the carrier frequency offset compensation, and itself and training sequence are carried out cross-correlation operation, obtains correlated series;

The thresholding judge module is used for judging by thresholding, finds out in the described correlated series the location index value I in last footpath _FirstLocation index value I with last footpath _Last

The delay spreading parameter computing module is used to receive I _FirstAnd I _Last, and use I _LastDeduct I _FirstObtain delay spreading parameter L _Ds

11, estimation unit as claimed in claim 10 is characterized in that, described correlated series generation module comprises:

Pending sequence generates submodule, is used to receive the received signal through the carrier frequency offset compensation, and therefrom intercepts the one piece of data sequence as pending data sequence r _Pn(n);

The secondary processing sequence generation module is used to provide training sequence as secondary processing sequence s _Pn(m);

The cross-correlation operation submodule is used for pending data sequence r _Pn(n) and secondary processing sequence s _Pn(m) carry out cross-correlation operation, obtain correlated series s _Corr(l).

12, estimation unit as claimed in claim 11 is characterized in that:

Described pending sequence generates the pending data sequence r that submodule generates _Pn(n) length is 2 * L _{Ds, max}+ P, wherein n is the index of sample of signal, the value of n is from I _Sync-L _{Ds, max}To I _Sync+ L _{Ds, max}Each integer of+P-1 comprises I _Sync-L _{Ds, max}And I _Sync+ L _{Ds, max}+ P-1;

Wherein P is the length of frame head PN sequence, and frame synchronization position index is I _Sync, the requirement of the maximum delay expansion of receiver is L _{Ds, max}

13, estimation unit as claimed in claim 11 is characterized in that:

When this estimation unit was applied to the DTMB system, described secondary processing sequence generates submodule provided training sequence to be meant as secondary processing sequence: secondary processing sequence generates submodule and chooses local corresponding PN sequence of being preserved as secondary processing sequence s according to frame head mode that estimates and PN sequence phase _Pn(m), this sequence length is P, and wherein the value of index m is each integer from 0 to P-1, comprises 0 and P-1.

14, as claim 12 or 13 described estimation units, it is characterized in that:

When this estimation unit is applied to the DTMB system, when frame head PN pattern is PN420, P=420; When frame head PN pattern is PN595, P=595; When frame head PN pattern is PN945, P=945.

15, estimation unit as claimed in claim 12 is characterized in that:

Described cross-correlation operation submodule carries out cross-correlation operation to two sequences according to following formula, obtains correlated series s _Corr(l):

$s_{\mathrm{corr}}\left(l\right)=\underset{m=0}{\overset{P-1}{\mathrm{\&Sigma;}}}\mathrm{conj}\left(s_{\mathrm{pn}}\left(m\right)\right)\&CenterDot;r_{\mathrm{pn}}(I_{\mathrm{sync}}-L_{\mathrm{ds},\max }+l+m)$

Wherein conj () expression complex conjugate operation, the value of l is to 2 * L from 0 _{Ds, max}Each integer, comprise 0 and 2 * L _{Ds, max}

16, as claim 10 to 13,15 each described estimation units, it is characterized in that described thresholding judge module comprises:

The threshold value sub module stored is used to store set power factor threshold T _Power

Power sequence generates submodule, is used to ask the power sequence of correlated series;

Average power generates submodule, is used to ask the average power P of power sequence _{Corr, avg}

Index value is judged submodule, is used for from front to back sample-by-sample ground and judges that whether each sample of power sequence is greater than T _Power* P _{Corr, avg}, and will be greater than T _Power* P _{Corr, avg}The index value of first sample be designated as I _First, will be greater than T _Power* P _{Corr, avg}The index value of last sample be designated as I _Last

17, estimation unit as claimed in claim 16 is characterized in that:

When this method was applied to the DTMB system, described threshold value was provided with the set power factor threshold T of submodule _PowerBe 16.

18, estimation unit as claimed in claim 16 is characterized in that, described index value is judged submodule sample-by-sample ground judgement from front to back power sequence P _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}, and will be greater than T _Power* P _{Corr, avg}The index value of first sample be designated as I _First, will be greater than T _Power* P _{Corr, avg}The index value of last sample be designated as I _LastSpecifically be meant:

Index value is judged S set of submodule structure, and power sequence P is judged on sample-by-sample ground from front to back _CorrWhether each sample (l) is greater than T _Power* P _{Corr, avg}If greater than, the index value of this sample is designated as the element of S set; Otherwise, be not designated as the element of S set; If S set is not null set, then index value judges that submodule is designated as I with first element value of S set _First, this is the location index value in last footpath, and last element value is designated as I _Last, this is the location index value in last footpath.

Images