CN1599364A - Data equalization method for burst communication - Google Patents

Abstract

This invention discloses a data balance method for burst communication based on the following: theory: An initial estimate vector of a channel parameter carries out a balance process to the received baseband deformation signals to output decisive data to recover the original delivery signal, the steps as estimating channel parameter vectors to the received signals based on the initial vector of the channel parameter and received signals to get estimated vectors to figure out the path branch, measured increment then to get measures of all necessary states and find out the state with the minimum branch measure to backtrack its path and transmit the decisive signal at the time scale.

Description

The data balancing method that is used for burst communication

Technical field

The invention belongs to communication technical field, particularly relate to a kind of data balancing method that is used for burst communication, promptly follow the tracks of and maximum likelihood equalization methods D-PSP based on the combined channel of depression of order.

Background technology

Burst communication is that a kind of efficient communication when being subjected to physics or electronic jamming ensures means, it is a kind of necessary complement of general communication, when the communication mode of routine is subjected to losing efficacy such as influence such as electronic jamming, earthquake and nuclear explosion, burst communication just can be brought into play the communication support effect at this crucial moment, therefore has special critical role in the communications field.Burst communication adopts the mode of burst to transmit information, transmission time is at random, and the duration is short, can signal be transmitted in moment to finish, therefore have advantages such as antijamming capability is strong, reliability is high, good concealment, be suitable for the civil and military communication under the particular surroundings.For a long time, burst communication theory and correlation technique thereof are special communications fields always, main in the world developed country adheres to the research work of this respect all the time, and has formed perfect burst communication system, and China also progressively carries out the key technology research of this respect.

Typical burst communication system operates mainly in shortwave, ultrashort involving under the channel circumstance such as meteoric trail.With respect to general communication system, burst communication has following characteristics.

1. the transmission channel environment is abominable.Because burst communication mainly utilizes ionospheric reflection or scattering process to communicate, transmission signals can be subjected to such as multipath transmisstion, decline and background noise etc. are serious fast disturbs, and therefore must carry out equilibrium treatment to data at receiving terminal;

2. communications distance.Utilize common omnidirectional antenna or directional antenna, communication distance requires to reach 1000～2000 kilometers;

3. communicate by letter effective duration is short.Under the burst communication situation,, require system in the time of 1-20ms, to enter stable operating state usually in order to guarantee the valid data percent of pass.

In view of above characteristics, in order to improve data pass rate, burst communication system is very high for the requirement of key technologies such as channel parameter tracking and channel equalization.

The present invention is exactly under such background, launches with the characteristics of communication in conjunction with the burst communication condition.The burst communication channel has mutability and randomness, generally is that unit transmits with the packet data frame, utilizes short transmission frame head to come the bid parameter as training sequence.On this basis, start follow-up data balancing according to the channel parameter first guess that captures.From burst communication mechanism as can be seen, what the used Burst Channel of transferring data frames may be with consecutive frame is different, so the equilibrium treatment of data should be that unit carries out with the Frame.

According to above communication mechanism, the structure of burst communication system as shown in Figure 1, wherein the data processing operation of receiver is mainly finished by high speed digital signal processor DSP.As can be seen, the transmission course of burst communication signal is: send signal x _k(k=1,2 ...) at first become s through coding _k(k=1,2 ...) to strengthen antijamming capability, carry out then serial/parallelly being for conversion into I road, Q road two paths of signals, and become baseband modulation signal according to constellation mapping.Modulated signal is carried out waveform by raised cosine filter and is formed, and is radiofrequency signal through upward frequency conversion again, by antenna transmission to channel.Transmission signals can produce signal distortion through the pollution of multipath effect and noise, so the receiving terminal at burst communication will carry out equilibrium treatment to the distorted signal that receives, recover original transmission signal, concrete receiving course is: the radiofrequency signal that antenna receives is through down-conversion and Filtering Processing, the baseband receiving signals d that obtains distorting _k(k=1,2 ...), at first according to the training sequence that receives, what promptly transmission frame head estimated channel parameter just estimates vector Utilize initial channel parameter Estimation vector then

The log-on data balanced unit carries out equilibrium treatment to recover original transmission signal to the effective information data after the transmission frame head.The judgement data of equalizer output are s ' _k(k=1,2 ...), through the decoded x ' that finally is output as _k(k=1,2 ...).Accurately estimating to have perfect performance as if the equalization methods that adopts on the basis of channel parameter, then the judgement data s ' of equalizer output _k(k=1,2 ...) should equal s _k(k=1,2 ...), corresponding decoding output x ' _k(k=1,2 ...) also equal original transmission signal x _k(k=1,2 ...), can realize the reliable reception of burst communication data, as shown in Figure 1.

As seen, the performance of data balancing is directly determining the systematic function of burst communication.Effective duration of burst communication is short, and channel acquisition, reaching frequency difference synchronously, to proofread and correct the desired data expense very big to the influence of payload data percent of pass.Therefore, the research processing speed fast, be suitable for all little data balancing method of Burst Channel characteristic, computing and storage overhead and become the key technology that must solve in the burst communication.

At present, existing all kinds of general balancing technique is divided into linear equalization and nonlinear equalization two classes.Linear equalization mainly is meant the linear equalizer based on transverse filter structure; Nonlinear equalization mainly comprises decision feedback equalization DFE, adaptive equalization and maximum likelihood equilibrium.

For linear equalization, simple in structure based on the linear equalizer of transverse filter structure, be easy to realize, but can't offset serious channel distortions, do not reach the systematic function requirement.

In nonlinear equalization, also there are problems.For example: 1. decision feedback equalization DFE is by eliminating that part of intersymbol interference that decision signal causes from current estimated value, reach portfolio effect preferably, but because the adjustment of adaptive algorithm is depended in the renewal of equalizer coefficients, therefore can't follow the tracks of the quick variation of channel, even the situation that do not restrain of coefficient appears, referring to U.S. John G Prokies. " DigitalCommunication " [M] .Beijing:Publishing House of Electronics Industry, 1998 and Wang Junfeng, Zhang Bin, Song Guo township " based on the nonlinear channel decision feedback equalization algorithm of wavelet transformation " system engineering and electronic technology, 2002 12 phases etc.2. adaptive equalization is utilized the contrary principle of simulation of sef-adapting filter, offsets the intersymbol interference that multipath channel produces, and its convergence rate is slow, and the work prerequisite of adaptive algorithm is that the channel parameter statistical property is constant.Under the burst communication environment, it is non-stationary that the sudden change of the characteristic of channel has signal, it is followed the tracks of and convergence capabilities has no way of adapting to, can't reach the performance requirement of burst communication system, referring to U.S. Simon Haykin.Adaptive Filter Theory (Fourth Edition) [M] .Beijing:Publishing House of Electronics Industry, 2002. and Xu Mingyuan, Lin Huafang, Qiu Gongan " based on the simulation study of the Adaptive Equalization System of LMS algorithm " system emulation journal, 2003 02 phases etc.3. based on Viterbi (VITERBI) though the balanced MLSE of the maximum likelihood of algorithm has good receptivity, but because the channel estimating of front end has introduced decision delay and computation complexity is big, therefore, also can't satisfy the real time signal processing requirement of burst communication, referring to U.S. John G Prokies. " DigitalCommunication " [M] .Beijing:Publishing House of Electronics Industry, 1998 and Song Liang, Hu Bo, insult mediate booth " based on a kind of NEW ADAPTIVE MLSE receiver of fast fading channel " electronic letters, vol, 2002,30 (5): 723-726.But on the balanced basis of maximum likelihood, follow the tracks of and maximum likelihood equalization methods PSP based on the combined channel of handling by survivor path, real-time channel tracking is incorporated in the survivor path processing of viterbi algorithm (VA), overcome the decision delay problem of traditional maximum likelihood equilibrium, has optimum in theory receptivity, referring to Italian Riccardo Raheli, Andreas Polydoros, Ching-Kae Tzou.Per-survivor processing:A general approach to MLSE inuncertain environments [J] .IEEE Trans on Comm, 1995, yet 43:354-364., owing to handle the PSP method when carrying out channel tracking and Data Detection by survivor path, the computing and the memory space that all need ML, wherein M is modulation system number, L is an equivalent channel length, is difficult on the engineering realize.At this situation, prior art mainly adopts following three class methods to simplify both at home and abroad:

The first kind, be that Canadian J.Omidi and Rollins propose respectively being used for the adaptive mode block method of channel estimating in associating channel tracking and the maximum likelihood equalization methods, carry out designs simplification, to reduce overall complexity, referring to Canadian J.Omidi, P.G.Gulak, and S.Pasupathy, Parallel Structuresfor Joint Channel Estimation and Data Detection Over Fading Channels[J], IEEETransactions on Selected Areas of Communications, 1998,16 (5): 1616-1629. and Canadian Rollins, M.E., Simmons, S.J.Simplified per-survivor Kalman processingin fast frequency-selective fading channels[J] .IEEE Trans on Comm., 1997,45 (5): 544-553.

Second class, be that Italian G Castellini proposes to reduce each number of times of channel tracking constantly, under the little situation of performance loss, only the state with minimum branch metric is carried out the channel tracking of a survivor path, and be the MSP method as the common signal channel parameter of these all states of moment with the channel parameter that estimates, referring to Italian G Castellini, F Conti, E Del, Re, L Pierucci.A continuouslyadaptive MLSE receiver for mobile communications:algorithm and performance[J] .IEEE Trans on Comm, 1997,45 (1): 80-89., or be used in combination the PSP method on this basis and the MSP method reduces operand, referring to Korea S Jung Suk Joo, Seung Chul Hong, YongHoon Lee.Adaptive MLSE receiver:hybrid of per-survivor processing andtentative decision MLSE[J] .Electronics Letters, 2000,36 (7): 678-680..

The 3rd class, be with channel tracking and Maximum Likelihood Detection associating and the method that combines with decision feedback equalization according to Finland Zhenhong Li etc., referring to Zhenhong Li, Piirainen, O., Mammela, A.Anew r educed-complexity a daptive P SP-MLSE receiver for GSM/EDGE systems[J] .IEEE International Symposium on Personal, Indoor and Mobile RadioCommunications, 2001 (1): 124-128., utilize decision feedback equalization processing section multipath, equivalence reduces the multipath channel length that associated detecting method is handled.

Above the whole bag of tricks all is to follow the tracks of with local module or the outside of maximum likelihood equalization methods PSP from combined channel to start with, and tries hard to by simplifying local algorithm or carrying out the operand that external treatment reduces this method.Yet the self structure of PSP method does not change, and still needs a large amount of computings and memory space, remains to be difficult to realize for the burst communication system that requires rapid data to handle.

The technology contents of invention

The objective of the invention is to solve above-mentioned existing combined channel follow the tracks of the computation complexity that exists with maximum likelihood equalization methods PSP too big, in burst communication, be difficult to the problem that realizes, variable Rate communication mechanism according to time dependent characteristics of burst communication channel and adaptive coding and modulating, a kind of data balancing method D-PSP based on channel tracking that is used for burst communication is provided, to satisfy the requirement of burst communication system fast processing, realize reliable reception to burst communication transmission data.

Technical scheme of the present invention is achieved in that

As one of key technology of burst communication, the equalization data method of employing is closely connected with the Burst Channel characteristic.As previously mentioned, follow the tracks of channel estimating delay and the error propagation problem that has overcome traditional equalization algorithm with equalization methods based on the combined channel of handling by survivor path, the real-time tracking and the data balancing of channel parameter vector are combined, have best receptivity.Yet its algorithm complex and googol make it be difficult to realize in the burst communication system of reality according to the computing memory space.Therefore, the present invention is from the burst communication principle, the data balancing method D-PSP that is used for burst communication of proposition based on channel tracking, and promptly a kind of combined channel based on depression of order is followed the tracks of and the maximum likelihood equalization methods.Effectively to reduce on the basis of computational complexity and memory space, fast processing and the quality data of realizing the burst communication signal receive, thereby improve the valid data percent of pass, open up the frontier of burst communication system signal processing key technology of new generation, for solid technical foundation is established in the exploitation and the realization of China's burst communication system.

As previously mentioned, in the burst communication system, the main effect of balanced part is to estimate vector according to channel parameter first To the base band distorted signal d that receives _k(k=1,2 ...) carrying out equilibrium treatment, data s ' is adjudicated in output _k(k=1,2 ...) to recover original transmission signal.The present invention utilizes the Burst Channel characteristics, the data balancing method D-PSP method that is used for burst communication that is provided based on channel tracking, and its concrete performing step can be divided into following four steps:

The first step is estimated vector according to channel parameter first With received signal d _k(k=1,2 ...), at each k constantly, adopt self adaptation RLS method sequence { d to received signal _k, d _K-1D _K-(L-1)Institute might state μ _kThe channel parameter vector estimate, obtain the constantly corresponding channel parameter estimation vector of k Wherein L is a Burst Channel length;

Second step is according to the constantly corresponding channel parameter estimation vector of the k that obtains Pass through formula ${\mathrm{\λ}}_N({\mathrm{\μ}}_k\→{\mathrm{\μ}}_{k+1})={|d_{k+1}-\hat{h}\left({\mathrm{\μ}}_k\right){S^{\′}}_N^T\left({\mathrm{\μ}}_{k+1}\right)|}^2$ Calculate received signal and when k, be carved into k+1 constantly from the possible state μ of institute _kTransfer to state μ _K+1Path branches tolerance increment λ _N(μ _k→ μ _K+1), in the formula: symbol () ^TRepresent the transposition computing, S ' _N(μ _K+1) representative for the judgement sequence that receives data s ' _K+1-(L-1)S ' _k, s ' _K+1, obtain by following formula

S′ _N(μ _k+1)＝S′ _N，Map(μ _k+1)+S′ _N，Trace(μ _k+1)

Wherein, N (N=1,2 ... L-1) combined channel is followed the tracks of and the exponent number of maximum likelihood equilibrium in order to reduce, S ' _{N, Map}(μ _K+1) represent state μ _K+1Directly the judgement sequence of Ying She a L-N information code element s ' _K+1-(L-N-1)... s ' _k, s ' _K+1, and S ' _{N, Trace}(μ _K+1) representative by individual path recall all the other N information code element that obtains the judgement sequence s ' _K+1-(L-1)... s ' _K+1-(L-N);

The 3rd step is according to the path branches tolerance increment λ that obtains _N(μ _k→ μ _K+1), received signal institute might state μ constantly at existing k _kPath branches tolerance Γ (μ _k) the basis on, utilize formula $\mathrm{\Γ}\left({\mathrm{\μ}}_{k+1}\right)=\underset{{\mathrm{\μ}}_k}{\min }[\mathrm{\Γ}\left({\mathrm{\μ}}_k\right)+{\mathrm{\λ}}_N({\mathrm{\μ}}_k\→{\mathrm{\μ}}_{k+1})]$ Obtain k+1 received signal sequence { d constantly _K+1, d _kD _K-(L-2)Institute might state μ _K+1Path branches tolerance Γ (μ _K+1);

The 4th step is according to the path branches tolerance Γ (μ of k+1 all states of the moment _K+1), find out state with minimum branch metric, carry out the path by state and recall minimum branch metric, export corresponding decision signal s ' of this moment _(k+1)-δ, wherein δ is the system determination degree of depth.

The first step in the said method, described employing self adaptation RLS method is to k moment received signal sequence { d _k, d _K-1D _K-(L-1)Institute might state μ _kThe channel parameter vector estimate, be to be undertaken by following four formulas:

$e\left({\mathrm{\μ}}_k\right)=d_k-{\hat{h}}^{*}\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)----\left(1\right)$

$g\left({\mathrm{\μ}}_k\right)=\frac{p\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)}{w+{S^{\′}}_N^{*}\left({\mathrm{\μ}}_k\right)p\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)}----\left(2\right)$

$p\left({\mathrm{\μ}}_k\right)=\frac{1}{w}[p\left({\mathrm{\μ}}_{k-1}\right)-g\left({\mathrm{\μ}}_k\right){S^{\′}}_N^{*}\left({\mathrm{\μ}}_k\right)p\left({\mathrm{\μ}}_{k-1}\right)]----\left(3\right)$

$\hat{h}\left({\mathrm{\μ}}_k\right)=\hat{h}\left({\mathrm{\μ}}_{k-1}\right)+g\left({\mathrm{\μ}}_k\right)e^{*}\left({\mathrm{\μ}}_k\right)----\left(4\right)$

Draw the constantly corresponding channel parameter estimation vector of k by above (1)～(4)

Symbol () in the formula ^*Represent conjugate operation, e (μ _k) be error signal, i.e. received signal d _kWith the desired output signal ${\hat{h}}^{*}\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)$ Difference, w is a forgetting factor, g (μ _k) be Kalman Kalman gain vector, p (μ _k) be correlation matrix, S ' _N(μ _k) represent k state μ constantly _kFor receiving data sequence { d _k, d _K-1D _K-(L-1)The judgement sequence s ' _K-(L-1)S ' _K-1, s ' _k.In self adaptation RLS method, when k=1 begins to start computing for the first time, p (μ _k) initial value p (μ ₀) get unit matrix,

Initial value be that channel parameter first estimated vector, promptly $\hat{h}\left({\mathrm{\μ}}_0\right)={\hat{h}}_0.$

In the said method second step, the judgement sequence S ' of described all the other N information code element _{N, Trace}(μ _K+1) recall by individual path and to obtain, concrete trace-back process is: at first find out arrival state μ _K+1Previous moment (be k constantly) state μ _k, draw corresponding decision signal s ' _K+1-(L-1)Find out arrival state μ more forward _kPrevious moment state (be k-1 constantly) state μ _K-1, draw corresponding decision signal s ' _K+1-(L-2)Similarly, find the state μ of (k+1)-n-hour always _(k+1)-N, draw corresponding decision signal s ' _K+1-(L-N)Thereby, the judgement sequence of N information code element of acquisition s ' _K+1-(L-1)... s ' _K+1-(L-N).

In the said method the 4th step, described state μ to k+1 minimum branch metric of the moment _K+1Carry out the path and recall, its detailed process is: if the system determination degree of depth is δ, at first find out arrival state μ _K+1Previous moment (be k constantly) state μ _kAnd then find out arrival state μ _kPrevious moment state (be k-1 constantly) state μ _K-1And the like, find (k+1)-δ state μ constantly always _(k+1)-δ, it is carried out bit decision, export corresponding decision signal s ' of this moment _(k+1)-δ

By above final output of four steps decision signal s ' _(k+1)-δThereby, finish the Data Detection of once uniting of the present invention.Return step (1) then, begin next associating Data Detection constantly.The rest may be inferred, can finish for received signal sequence d _k(k=1,2 ...) equilibrium, finally export system shown in Figure 1 and adjudicate data sequence s ' accordingly _k(k=1,2 ...), thereby the maximum likelihood reliable reception of realization burst communication data.

The present invention has following effect:

(1) computation complexity is low

Each is Path (D-PSP)=MM to the total path branches metric calculation number of times of all states constantly because of the present invention ^L-N=M ^L+1-N, and the total calculation times of each all state path branch metric of the moment of existing P SP method is Path (PSP)=MM ^L=M ^L+1, the path branches metric calculation that N rank method falls in visible the present invention has reduced M with respect to existing P SP method ^NDoubly.Table 1 is the computation complexity table of comparisons that algorithms of different is handled for each code element, and L is a channel length, and M is modulation system number.C wherein _CEThe computation complexity of a channel guess value of expression is determined by self adaptation RLS algorithm.As can be seen from Table 1, the complexity maximum of PSP method, the complexity of another prior art MSP is obviously greater than method of the present invention when M is big.As channel guess value complexity C _CEWhen determining with channel length L, it is big more that the modulation system is counted M, and the amount of calculation of the inventive method is littler.Table 2 has provided under the 16 system modulation systems, and channel length L is 4, when it is typical burst communication situation the complexity of the whole bag of tricks relatively, as seen, method of the present invention is all significantly decreasing aspect memory space and the computation complexity two.

The computation complexity of table 1 the whole bag of tricks relatively

Addition multiplication memory space channel guess value

PSP method (L+1) M ^L+1LM ^L+1(2+L) M ^LM ^LC _CE

MSP method (L+M) M ^LLM ^L2M ^L+ L C _CE

N rank D-PSP method

(L+1)M ^L+1-N?LM ^L+1-N (2+L)M ^L-N M ^L-N·C _CE

N＝1，2…(L-1)

The computation complexity of the whole bag of tricks relatively when table 2 L=4, M=16

Addition multiplication memory space channel guess value

PSP method 516 ⁵416 ⁵616 ⁴16 ⁴C _CE

MSP method 2016 ⁴416 ⁴216 ⁴+ 4 C _CE

1 rank D-PSP method 516 ⁴416 ⁴616 ³16 ³C _CE

2 rank D-PSP methods 516 ³416 ³616 ²16 ²C _CE

3 rank D-PSP methods 516 ²416 ²616 16C _CE

In practical communication, Frame for burst communication different modulating mode, count N (N=1 with the inventive method by selecting suitable depression of order, 2 ... L-1), can effectively control the computational complexity and the memory space of joint maximum likelihood equilibrium, each frame data is handled had roughly the same computing expense, thereby the limited instruction resource of the system that can make full use of is easy to realize the best quality of reception.

(2) systematic function is good

Not only computational complexity is low, method simple in the present invention, be easy to realization, and the systematic function of systematic function and existing combined channel tracking of employing and maximum likelihood equalization methods is approaching.Emulation experiment shows that under Quadrature Phase Shift Keying (QPSK) modulation system, the present invention gets different depression of orders and counts N, its signal to noise ratio E _s/ N ₀When 10.5dB was between 11.3dB, the error rate reached 1 * 10 respectively ^-3, as shown in Figure 4; Under 16 system phase shift keyings (16QAM) modulation system, the present invention gets different depression of orders and counts N, its signal to noise ratio E _s/ N ₀When 17.5dB was between 19dB, the error rate also reached 1 * 10 respectively ^-3, as shown in Figure 5.

Description of drawings

Fig. 1 is the burst communication system structure chart

Fig. 2 a is the schematic diagram of existing combined channel tracking and maximum likelihood equalization methods,

Fig. 2 b is the schematic diagram of depression of order combined channel tracking of the present invention and maximum likelihood equalization methods

Fig. 3 is the realization flow figure of the inventive method

Fig. 4 is the present invention under the QPSK modulation, the error rate of system curve chart of Computer Simulation

Fig. 5 is the present invention under the 16QAM modulation, the error rate of system curve chart of Computer Simulation

Embodiment

With reference to Fig. 1, the principle of existing burst communication system as stated in the Background Art.Wherein, the main effect of balanced part is to estimate vector according to channel parameter first

To the base band distorted signal d that receives _k(k=1,2 ...) carrying out equilibrium treatment, data s ' is adjudicated in output _k(k=1,2 ...) to recover original transmission signal, realize the reliable reception of data.The present invention utilizes the Burst Channel characteristics, associating channel tracking and maximum likelihood equalization methods PSP are simplified, provide a kind of combined channel to follow the tracks of and the maximum likelihood equalization methods based on depression of order, be data balancing method (D-PSP) method, to obtain optimum relatively receptivity and minimum computation complexity.

With reference to Fig. 2 a, existing combined channel is followed the tracks of with the principle of maximum likelihood equalization methods can represent that wherein CE represents the channel tracking computing with Viterbi (Viterbi) grid map, and channel length is L, each k M constantly ^LIndividual state μ _kComprise μ _{K, 1}, μ _{K, 2}... μ _{K, ML}, represent constantly corresponding L the receiving data sequence { d of k respectively _k, d _K-1D _K-(L-1)All possible state justify sequence s ' _k, s ' _K-1S ' _K-(L-1).Just estimate vector what capture channel parameter

Condition under, the judgement of code element has very strong correlation before and after the burst communication received signal, can think the k μ that gets the hang of constantly _kThe channel parameter estimation vector of path representative With go out μ _kThe channel parameter estimation vector of representative Almost constant.If channel length is L, then μ _K-1In have L-1 component to be retained in μ _kIn.

According to above principle, say intuitively, entering μ _kThe channel guess value basis on, if reduce the dimension of state, also can more correctly estimate out μ _kChannel, the advantage of doing like this is to reduce computation complexity and computing memory space.Because the bit decision value that front and back reservation component constantly promptly keeps is L-1, is L-1 so the depression of order of the inventive method is counted the maximum occurrences of N, i.e. N=1,2 ... L-1.

With reference to Fig. 2 b, with respect to Fig. 2 a the present invention at each k state μ constantly _kNumber be reduced to M ^L-NIndividual, i.e. μ _{K, 1}, μ _{K, 2}... μ _{K, ML-N}, they only represent L-N the judgement sequence of symhols s ' _k, s ' _K-1S ' _K-(L-N-1), and all the other N the judgement sequence of symhols s ' _K-(L-N)S ' _K-(L-1)Used the retention in a lasted moment, recall by the path of state and obtain.

With reference to Fig. 3, the present invention has provided based on the combined channel of depression of order and has followed the tracks of and the maximum likelihood equalization methods, i.e. the specific implementation flow process of data balancing method (D-PSP), and its process is:

At first, the burst communication receiver is carried out the initialization of system, realize catching of signal, and obtain initial channel parameter Estimation vector with synchronously

Then, on this basis to the payload data d after the frame head _k(k=1,2 ...) carry out following the tracks of and equilibrium treatment the corresponding judgement of output data s ' based on the combined channel of depression of order _k(k=1,2 ...).If the dicision depth of communication system Viterbi structure is δ, concrete as follows realization:

(1) obtains the channel parameter estimation vector

Realize catching of signal and synchronous at burst communication system, and obtain initial channel parameter Estimation vector Afterwards, at first estimate vector according to channel parameter first With received signal d _k(k=1,2 ...), at each k constantly, adopt adaptive recursive least-squares RLS method, by following four formulas sequence { d to received signal _k, d _K-1D _K-(L-1)Institute might state μ _kThe channel parameter vector estimate, obtain the constantly corresponding channel parameter estimation vector of k That is:

$e\left({\mathrm{\μ}}_k\right)=d_k-{\hat{h}}^{*}\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)$

$g\left({\mathrm{\μ}}_k\right)=\frac{p\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)}{w+{S^{\′}}_N^{*}\left({\mathrm{\μ}}_k\right)p\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)}$

$p\left({\mathrm{\μ}}_k\right)=\frac{1}{w}[p\left({\mathrm{\μ}}_{k-1}\right)-g\left({\mathrm{\μ}}_k\right){S^{\′}}_N^{*}\left({\mathrm{\μ}}_k\right)p\left({\mathrm{\μ}}_{k-1}\right)]$

$\hat{h}\left({\mathrm{\μ}}_k\right)=\hat{h}\left({\mathrm{\μ}}_{k-1}\right)+g\left({\mathrm{\μ}}_k\right)e^{*}\left({\mathrm{\μ}}_k\right)$

In the formula: symbol () ^*Represent conjugate operation, e (μ _k) be error signal, i.e. received signal d _kWith the desired output signal ${\hat{h}}^{*}\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)$ Difference, w is a forgetting factor, g (μ _k) be Kalman (Kalman) gain vector, p (μ _k) be correlation matrix, S ' _N(μ _k) represent k state μ constantly _kFor receiving data sequence { d _k, d _K-1D _K-(L-1)The judgement sequence s ' _K-(L-1)S ' _K-1, s ' _k.In self adaptation RLS method, when k=1 begins to start computing for the first time, p (μ _k) initial value p (μ ₀) get unit matrix, Initial value be that channel parameter first estimated vector, promptly $\hat{h}\left({\mathrm{\μ}}_0\right)={\hat{h}}_0.$

(2) obtain path branches tolerance increment λ (μ _k→ μ _K+1)

The k that obtains according to step (1) is all state μ constantly _kCorresponding channel parameter estimation vector Further calculate received signal and when k, be carved into k+1 constantly from the possible state μ of institute _kTransfer to state μ _K+1Path branches tolerance increment λ _N(μ _k→ μ _K+1).For falling N rank D-PSP method, the judgement sequence of wherein calculating used information code element is to be that the individual path of N is recalled and obtained by the degree of depth, and concrete execution mode is:

Estimating

The basis on, be carved into k+1 when calculating k and fall N rank D-PSP method constantly from state μ by following formula _kTransfer to state μ _K+1Path branches tolerance increment λ _N(μ _k→ μ _K+1):

${\mathrm{\λ}}_N({\mathrm{\μ}}_k\→{\mathrm{\μ}}_{k+1})={|d_{k+1}-\hat{h}\left({\mathrm{\μ}}_k\right){S^{\′}}_N^T\left({\mathrm{\μ}}_{k+1}\right)|}^2$

Symbol () in the formula ^TRepresent the transposition computing, S ' _N(μ _K+1) representative for the judgement sequence that receives data s ' _K+1-(L-1)S ' _k, s ' _K+1, obtain by following formula

S′ _N(μ _k+1)＝S′ _N，Map(μ _k+1)+S′ _N，Trace(μ _k+1)

For the D-PSP method of falling the N rank, wherein S ' _{N, Map}(μ _K+1) represent state μ _K+1Directly the judgement sequence of Ying She a L-N information code element s ' _K+1-(L-N-1)... s ' _k, s ' _K+1, and S ' _{N, Trace}(μ _K+1) then represent by individual path recall all the other N information code element that obtains the judgement sequence s ' _K+1-(L-1)... s ' _K+1-(L-N), the traceback depth of its individual path is depression of order exponent number N, promptly

S ' _{N, Trace}(μ _K+1)=Path_trace[μ _K-N→ μ _k] Path_trace[in the formula] represent the trace-back operation of individual path, concrete trace-back process is: at first find out arrival state μ _K+1Previous moment (be k constantly) state μ _k, draw corresponding decision signal s ' _K+1-(L-1)Find out arrival state μ more forward _kPrevious moment state (be k-1 constantly) state μ _K-1, draw corresponding decision signal S ' _K+1-(L-2)Similarly, find the state μ of (k+1)-n-hour always _(k+1)-N, draw corresponding decision signal s ' _K+1-(L-N)Thereby, the judgement sequence of N information code element of acquisition s ' _K+1-(L-1)... s ' _K+1-(L-N).

(3) obtain path branches tolerance Γ (μ _K+1)

Obtain state μ according to step (2) _kTransfer to state μ _K+1Path branches tolerance increment λ _N(μ _k→ μ _K+1), received signal institute might state μ constantly at existing k _kPath branches tolerance Γ (μ _k) the basis on, obtain k+1 received signal sequence { d constantly _K+1, d _kD _K-(L-2)Institute might state μ _K+1Path branches tolerance Γ (μ _K+1), concrete execution mode is:

Calculating path branches tolerance increment λ _N(μ _k→ μ _K+1) the basis on, further obtain state μ _K+1Survivor path branch metric function gamma (μ _K+1), it is by the state of arrival μ _K+1The decision of minimum survivor path state, promptly

$\mathrm{\Γ}\left({\mathrm{\μ}}_{k+1}\right)=\underset{{\mathrm{\μ}}_k}{\min }[\mathrm{\Γ}\left({\mathrm{\μ}}_k\right)+{\mathrm{\λ}}_N({\mathrm{\μ}}_k\→{\mathrm{\μ}}_{k+1})]$

$=\underset{{\mathrm{\μ}}_k}{\min }[\mathrm{\Γ}\left({\mathrm{\μ}}_k\right)+{|d_{k+1}-\hat{h}\left({\mathrm{\μ}}_k\right){S^{\′}}_N^T\left({\mathrm{\μ}}_{k+1}\right)|}^2]$

(4) output decision signal s ' _(k+1)-δ

The k+1 that obtains according to step (3) is all state μ constantly _K+1Path branches tolerance Γ (μ _K+1), find out state with minimum branch metric, carry out the path by state and recall minimum branch metric, export corresponding decision signal s ' of this moment _(k+1)-δConcrete trace-back process is: if the system determination degree of depth is δ, at first find out arrival state μ _K+1Previous moment (be k constantly) state μ _k, and then find out arrival state μ _kPrevious moment state (be k-1 constantly) state μ _K-1And the like, find (k+1)-δ state μ constantly always _(k+1)-δ, it is carried out bit decision, export corresponding decision signal s ' of this moment _(k+1)-δThereby finish the once associating Data Detection of falling N rank D-PSP method.

Can finish k+1 according to above-mentioned steps and follow the tracks of and balanced D-PSP processing based on the combined channel of depression of order constantly, export corresponding decision signal s ' _(k+1)-δReturn step (1) then, begin next constantly new D-PSP associating Data Detection ... and the like, finish for received signal sequence d _k(k=1,2 ...) equilibrium, finally export corresponding judgement data sequence s ' shown in Figure 1 _k(k=1,2 ...), thereby effectively reducing under the condition of computational complexity and memory space, realize the maximum likelihood reliable reception of burst communication data.

After burst communication system adopted the inventive method, its performance can provide by following computer simulation experiment.

Experiment 1:

The Watterson model that simulated channel adopts CCIR to recommend, Burst Channel fading rate f _rRandom sequence generation rate f with equivalent Rayleigh multidiameter fading channel _cBetween following relational expression is arranged.

$f_r=1.47*\frac{4.4f_c}{2\mathrm{\π}(2m+1)\sqrt{2}}$

Be under 4 the channel condition at the typical 1Hz fading rate of burst communication, channel length L, the message transmission rate that emulation is adopted is 32kbps, each symbol sample 5 point, carrier frequency is 40MHz, noise is an additive white Gaussian noise, the system transmissions data do not add and interweave and error correction coding is handled, and when modulation system was QPSK, the ber curve that emulation draws system as shown in Figure 4.For the ease of relatively, also simulated the Viterbi equalization methods under the known channel, PSP method and MSP method among Fig. 4, therefore can regard the ber curve under the known channel as the performance upper bound of the whole bag of tricks.As can be seen from Fig. 4, single order D-PSP method performance and PSP method under the QPSK modulation system are approaching, and the bit error rate performance of second order D-PSP method only descends about 0.5dB than the PSP method.The bit error rate performance of MSP method under Burst Channel is between the second order and three rank D-PSP methods of D-PSP, but its computation complexity is counted the increase of M along with the modulation system and increased sharply.Signal to noise ratio E _s/ N ₀When 10.5dB was between 11.3dB, the error rate of D-PSP method reached 1 * 10 respectively ^-3, be that error rate of system can reach 10 after 1/2 convolution code or the Turbo code coding through code efficiency ^-6Thereby, can realize the reliable reception of burst communication data.

Experiment 2:

Be similar to experiment 1, the Watterson model that simulated channel also adopts CCIR to recommend, and system emulation conditions such as channel fading rate, channel length, symbol sample point, carrier frequency are all identical.Under the modulation system of 16QAM, the message transmission rate that emulation is adopted is 64kbps, and the system transmissions data do not add and interweave and the error correction coding processing, and the ber curve that obtains as shown in Figure 5.For the ease of relatively, also simulated the Viterbi equalization methods under the known channel, PSP method and MSP method among Fig. 5, therefore can regard the ber curve under the known channel as the performance upper bound of the whole bag of tricks.Be similar to the result of Fig. 4, as can be seen from Fig. 5, single order D-PSP method performance and PSP method under the 16QAM modulation system are approaching, and the bit error rate performance of second order D-PSP method is than the PSP method 0.8dB that only descends.The bit error rate performance of MSP method under Burst Channel is between the second order and three rank D-PSP methods of D-PSP, but its computation complexity is counted the increase of M along with the modulation system and increased sharply.As signal to noise ratio E _s/ N ₀When 17.5dB was between 19dB, the error rate of D-PSP method also reached 1 * 10 respectively ^-3, be that error rate of system can reach 10 after 1/2 convolution code or the Turbo code coding through code efficiency ^-6

Above-mentioned two experimental results all prove, adopt method of the present invention can guarantee that the quality data of burst communication receives.

Claims (4)

1. data balancing method D-PSP who is used for burst communication, carry out as follows:

1. estimate vector by channel parameter first With received signal d _k(k=1,2 ...), at each k constantly, adopt self adaptation RLS method sequence { d to received signal _k, d _K-1D _K-(L-1)Institute might state μ _kThe channel parameter vector estimate, obtain the constantly corresponding channel parameter estimation vector of k Wherein L is a Burst Channel length;

2. by the constantly corresponding channel parameter estimation vector of the k that obtains Pass through formula ${\mathrm{\λ}}_N({\mathrm{\μ}}_k\→{\mathrm{\μ}}_{k+1})={|d_{k+1}-\hat{h}\left({\mathrm{\μ}}_k\right){S^{\′}}_N^T\left({\mathrm{\μ}}_{k+1}\right)|}^2$ Calculate received signal and when k, be carved into k+1 constantly from the possible state μ of institute _kTransfer to state μ _K+1Path branches tolerance increment λ _N(μ _k→ μ _K+1), in the formula: symbol () ^TRepresent the transposition computing, S ' _N(μ _K+1) representative for the judgement sequence that receives data s ' _K+1-(L-1)S ' _k, s ' _K+1, obtain by following formula

S′ _N(μ _k+1)＝S′ _N，Map(μ _k+1)+S′ _N，Trae(μ _k+1)

Wherein, N (N=1,2 ... L-1) combined channel is followed the tracks of and the exponent number of maximum likelihood equilibrium in order to reduce, S ' _{N, Map}(μ _K+1) represent state μ _K+1Directly the judgement sequence of Ying She a L-N information code element s ' _K+1-(L-N-1)... s ' _k, s ' _K+1, and S ' _{N, Trace}(μ _K+1) then represent by individual path recall all the other N information code element that obtains the judgement sequence s ' _K+1-(L-1)... s ' _K+1-(L-N);

3. by the path branches tolerance increment λ that obtains _N(μ _k→ μ _K+1), received signal institute might state μ constantly at existing k _kPath branches tolerance Γ (μ _k) the basis on, utilize formula $\mathrm{\Γ}\left({\mathrm{\μ}}_{k+1}\right)=\underset{{\mathrm{\μ}}_k}{\min }[\mathrm{\Γ}\left({\mathrm{\μ}}_k\right)+{\mathrm{\λ}}_N({\mathrm{\μ}}_k\→{\mathrm{\μ}}_{k+1})]$ Obtain k+1 received signal sequence { d constantly _K+1, d _kD _K-(L-2)Institute might state μ _K+1Path branches tolerance Γ (μ _K+1);

4. measure Γ (μ by the path branches of k+1 all states of the moment _K+1), find out state with minimum branch metric, carry out the path by state and recall minimum branch metric, export corresponding decision signal s ' of this moment _(k+1)-δ, wherein δ is the system determination degree of depth.

2. method according to claim 1, wherein the 1. described employing self adaptation of step RLS method is to k moment received signal sequence { d _k, d _K-1D _K-(L-1)Institute might state μ _kThe channel parameter vector estimate, be to draw the constantly corresponding channel parameter estimation vector of k by following four formulas

$e\left({\mathrm{\μ}}_k\right)=d_k-{\hat{h}}^{*}\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)$

$g\left({\mathrm{\μ}}_k\right)=\frac{p\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)}{w+{S^{\′}}_N^{*}\left({\mathrm{\μ}}_k\right)p\left({\mathrm{\μ}}_{k-1}\right){S^{\′}}_N^T\left({\mathrm{\μ}}_k\right)}$

$p\left({\mathrm{\μ}}_k\right)=\frac{1}{w}[p\left({\mathrm{\μ}}_{k-1}\right)-g\left({\mathrm{\μ}}_k\right){S^{\′}}_N^{*}\left({\mathrm{\μ}}_k\right)p\left({\mathrm{\μ}}_{k-1}\right)]$

$\hat{h}\left({\mathrm{\μ}}_k\right)=\hat{h}\left({\mathrm{\μ}}_{k-1}\right)+g\left({\mathrm{\μ}}_k\right)e^{*}\left({\mathrm{\μ}}_k\right)$

In the formula: symbol () ^*Represent conjugate operation, e (μ _k) be error signal, i.e. received signal d _kWith the desired output signal Difference, w is a forgetting factor, g (μ _k) be the kalman gain vector, p (μ _k) be correlation matrix, S ' _N(μ _k) represent k state μ constantly _kFor receiving data sequence { d _k, d _K-1D _K-(L-1)The judgement sequence s ' _K-(L-1)S ' _K-1, s ' _k;

When k=1 begins to start computing for the first time, p (μ _k) initial value p (μ ₀) get unit matrix, Initial value be that channel parameter first estimated vector, promptly $\hat{h}\left({\mathrm{\μ}}_0\right)={\hat{h}}_0.$

3. method according to claim 1, wherein the judgement sequence S ' of 2. described all the other N information code element of step _{N, Trace}(μ _K+1) recall the detailed process that obtains by individual path and be: at first find out arrival state μ _K+1Previous moment, i.e. k state μ constantly _k, draw corresponding decision signal s ' _K+1-(L-1)Find out arrival state μ more forward _kThe previous moment state, i.e. k-1 state μ constantly _K-1, draw corresponding decision signal s ' _K+1-(L-2)Similarly, find the state μ of (k+1)-n-hour always _(k+1)-N, draw corresponding decision signal s ' _K+1-(L-N)Thereby, the judgement sequence of N information code element of acquisition s ' _K+1-(L-1)... s ' _K+1-(L-N).

4. method according to claim 1, the wherein 4. described state μ of step to k+1 minimum branch metric of the moment _K+1The detailed process of carrying out recalling in the path is: if the system determination degree of depth is δ, at first find out arrival state μ _K+1Previous moment, i.e. k state μ constantly _k, and then find out arrival state μ _kThe previous moment state, i.e. k-1 state μ constantly _K-1And the like, find (k+1)-δ state μ constantly always _(k-1)-δ, it is carried out bit decision, export corresponding decision signal s ' of this moment _(k+1)-δ

Images