CN104253771A - Multi-parameter joint estimation method and apparatus - Google Patents

Abstract

The invention provides a multi-parameter joint estimation method and apparatus. The apparatus comprises a channel estimation unit which performs channel estimation on received signals through utilizing locally-generated channel state information reference signals (CSI-RS) so as to obtain the channel state information of the received signals, a fast Fourier inverse transformation unit transforms the channel state information from frequency domain to time domain, a signal power computing unit which calculates the signal power of the channel state information in the time domain, and a joint estimation unit which estimates the time offset, delay spread and noise power of the received signals through utilizing the signal power of the channel state information in the time domain. With the multi-parameter joint estimation method and apparatus of the invention adopted, multiple parameters including the time offset, delay spread and noise power can be estimated jointly, and therefore, the PDSCH demodulation performance of R11UE which is configured as QCL B type can be improved, and the performance of data throughput can be improved.

Description

Multi-parameter combined estimation method and device

Technical field

The present invention relates to radio communication, particularly relate to a kind of multi-parameter combined estimation method and device.

Background technology

At LTE-A(Long Term Evolution-Advanced, the Long Term Evolution strengthened) R11(Release11 of system, version 11) in, in order to support TM10(Transmission Mode10, transmission mode 10) in multipoint cooperative (CoMP, Coordinate Multipoint) transmission, define the accurate colocated (quasi co-location, QCL) between various antenna port.For Serving cell, the UE(User Equipment configured in TM10, subscriber equipment, referred to as user) be configured to one of following two kinds of QCL types.

Type A, for PDSCH(Physical Downlink Shared CHannel, Physical Downlink Shared Channel) CRS(Common Reference Signal, public reference signal), CSI-RS(Channel State Information Reference Signals, channel state information reference signals) and DM-RS(DeModulation Reference Signal, demodulated reference signal) accurate colocated about " Doppler frequency shift, doppler spread, average delay and delay spread " can be assumed to be.

Type B, the DM-RS for PDSCH indicated by high level and specific CSI-RS are assumed to be the accurate colocated about " Doppler frequency shift, doppler spread, average delay and delay spread ".For each CSI-RS, network can pass through RRC(Radio Resource Control, wireless heterogeneous networks) signaling indicates CSI-RS port and the CRS port of a community, can be assumed to be the accurate colocated about " Doppler frequency shift and doppler spread ".

In some CoMP scenes, as shown in Figure 1, a UE is synchronous with Serving cell TP#1, and meanwhile, its other transfer points (TP#2) concentrated from CoMP receive PDSCH.But, (timing offset) Δ t may be offset if having time between Serving cell (TP#1) and PDSCH community (TP#2).Current, if this UE is configured to type B, then the CRS for PDSCH demodulation can not be used to estimate average delay (average delay) and delay spread (delay spread).Therefore, the new time migration estimation method of research and delay spread method of estimation is needed.

In addition, in descending CoMP, the Noise and Interference that the Noise and Interference that CRS-RS is subject to and CRS are subject to is different.CSI-RS unit is only subject to the Noise and Interference outside CoMP collection.But the estimation based on CRS provides the summation disturbed in CoMP collection external noise and CoMP collection.Therefore, in TM10, based on CSI-RS Noise and Interference power estimate for CQI(Channel Quality Indicator, CQI) calculate be required.

Above it should be noted that, just conveniently to technical scheme of the present invention, clear, complete explanation is carried out to the introduction of technical background, and facilitate the understanding of those skilled in the art to set forth.Only can not think that technique scheme is conventionally known to one of skill in the art because these schemes have carried out setting forth in background technology part of the present invention.

Summary of the invention

Embodiments provide a kind of multi-parameter combined estimation method and device, to solve the problem pointed by background technology.

According to the first aspect of the embodiment of the present invention, provide a kind of multi-parameter Combined estimator device, wherein, described device comprises:

Channel estimating unit, it utilizes the local channel state information reference signals (CSI-RS) generated to carry out channel estimating to received signal, obtains the channel condition information of described Received signal strength;

Inverse fast Fourier transform unit, described channel condition information is transformed to time domain from frequency domain by it;

Signal power computing unit, it calculates the signal power of the channel condition information in time domain;

Combined estimator unit, it utilizes, and the signal power of the channel condition information in described time domain estimates the time migration of described Received signal strength, delay is expanded and noise power.

According to the second aspect of the embodiment of the present invention, provide a kind of multi-parameter combined estimation method, wherein, described method comprises:

Utilize the local CSI-RS generated to carry out channel estimating to received signal, obtain the channel condition information of described Received signal strength;

Inverse fast Fourier transform is carried out to described channel condition information, described channel condition information is transformed to time domain from frequency domain;

Calculate the signal power of the channel condition information in time domain;

Utilize the signal power of the channel condition information in described time domain estimate described Received signal strength time migration, postpone expansion and noise power.

The beneficial effect of the embodiment of the present invention is, by the multi-parameter combined estimation method of the embodiment of the present invention and device to comprising time migration, multiple parameters of delay spread and noise power carry out Combined estimator, improve the demodulation performance of the PDSCH of the R11UE being configured to QCL type B, improve the performance of data throughout.

With reference to explanation hereinafter and accompanying drawing, disclose in detail particular implementation of the present invention, specifying principle of the present invention can adopted mode.Should be appreciated that, thus embodiments of the present invention are not restricted in scope.In the spirit of claims and the scope of clause, embodiments of the present invention comprise many changes, amendment and are equal to.

The feature described for a kind of execution mode and/or illustrate can use in one or more other execution mode in same or similar mode, combined with the feature in other execution mode, or substitutes the feature in other execution mode.

Should emphasize, term " comprises/comprises " existence referring to feature, one integral piece, step or assembly when using herein, but does not get rid of the existence or additional of one or more further feature, one integral piece, step or assembly.

Accompanying drawing explanation

With reference to following accompanying drawing a lot of aspects that the present invention may be better understood.Parts in accompanying drawing are not proportional draftings, and just in order to illustrate principle of the present invention.For the ease of illustrating and describing some parts of the present invention, in accompanying drawing, corresponding part may be exaggerated or reduce.The element described in an accompanying drawing of the present invention or a kind of execution mode and feature can combine with the element shown in one or more other accompanying drawing or execution mode and feature.In addition, in the accompanying drawings, similar label represents parts corresponding in several accompanying drawing, and can be used for instruction more than the corresponding component used in a kind of execution mode.

In the accompanying drawings:

Fig. 1 is the system pattern schematic diagram of CoMP scene;

Fig. 2 is the composition schematic diagram of an execution mode of the multi-parameter Combined estimator device of the embodiment of the present invention;

Fig. 3 is the composition schematic diagram of another execution mode of the multi-parameter Combined estimator device of the embodiment of the present invention;

Fig. 4 is the composition schematic diagram of Combined estimator unit in the execution mode of Fig. 2 or Fig. 3;

Fig. 5 is the composition schematic diagram of embodiment of the present invention application case in terminal equipment;

Fig. 6 is the schematic diagram of the multi-parameter Combined estimator scheme in time domain;

Fig. 7 is the process chart of the multi-parameter Combined estimator of Combined estimator unit;

Fig. 8 is the schematic diagram of the execution mode that time migration is estimated;

Fig. 9 is the schematic diagram of another execution mode that time migration is estimated;

Figure 10 is the demodulation performance schematic diagram adopting multi-parameter Combined estimator;

Figure 11 is the flow chart of the multi-parameter combined estimation method of the embodiment of the present invention.

Embodiment

With reference to accompanying drawing, by specification below, the aforementioned and further feature of the embodiment of the present invention will become obvious.These execution modes are exemplary, are not limitations of the present invention.In order to enable those skilled in the art to easily to understand principle of the present invention and execution mode, embodiments of the present invention are described for the multi-parameter combined estimation method in lte-a system, but be appreciated that, the embodiment of the present invention is not limited to said system, all applicable for the other system relating to multi-parameter Combined estimator.

Embodiment 1

Embodiments provide a kind of multi-parameter Combined estimator device.Fig. 2 is the composition schematic diagram of this device, please refer to Fig. 2, this device comprises: channel estimating unit 21, inverse fast Fourier transform (IFFT, Inverse Fast Fourier Transform) unit 22, signal power computing unit 23 and Combined estimator unit 24, wherein

Channel estimating unit 21 utilizes the local CSI-RS generated to carry out channel estimating to received signal, obtains the channel condition information of described Received signal strength.

Wherein, the device of the present embodiment can also comprise a CSI-RS copy generation unit to generate CSI-RS, so that the channel estimating unit 21 being supplied to the present embodiment carries out channel estimating accordingly.Wherein, the mode that CSI-RS copy generation unit generates CSI-RS copy can be realized by existing means, does not repeat them here.

Wherein, the present embodiment does not limit concrete channel estimation methods, typically, can adopt ZF(Zero Forcing, ZF) channel estimation methods.

By the channel estimating unit 21 of the present embodiment, the channel condition information of Received signal strength can be obtained.

Described channel condition information is transformed to time domain from frequency domain by inverse fast Fourier transform (IFFT) unit 22.

Wherein, in order to carry out the Combined estimator of multi-parameter, channel condition information is transformed to time domain from frequency domain by IFFT unit 22 by the present embodiment.

Wherein, the present embodiment is transformed to example with IFFT, but not in this, as restriction, any protection range that channel condition information all can be contained in the present embodiment from the transform method that frequency domain transforms to time domain.

Signal power computing unit 23 calculates the signal power of the channel condition information in time domain.

Wherein, utilize the output of IFFT unit 22, signal power computing unit 23 can calculate the signal power of the channel condition information in time domain, to utilize, the signal power calculated carries out time migration, delay is expanded and the estimation of noise power.

Combined estimator unit 24 utilize the signal power of the channel condition information in described time domain estimate described Received signal strength time migration, postpone expansion and noise power.

Wherein, utilize the output of signal power computing unit 23, Combined estimator unit 24 can to comprising time migration, postpone expansion and multiple parameters of noise power are estimated, to carry out to received signal correcting and the solution doing data is further in harmonious proportion the calculating of feedback information.

Wherein, the mode that concrete estimated time offsets, the mode of delay spread and the mode of estimating noise power can adopt existing means to realize, the method that other embodiments of the present invention also can be adopted to provide realizes, and specifically will be described in detail in the examples below.

By the multi-parameter Combined estimator device of the embodiment of the present invention to comprising time migration, multiple parameters of delay spread and noise power carry out Combined estimator, improve the demodulation performance of PDSCH, improve the performance of data throughout.

Embodiment 2

The embodiment of the present invention additionally provides a kind of multi-parameter Combined estimator device.Fig. 3 is the composition schematic diagram of this device, please refer to Fig. 3, with the multi-parameter Combined estimator device of embodiment 1 unlike, in the present embodiment, this multi-parameter Combined estimator device is except comprising channel estimating unit 31, inverse fast Fourier transform (IFFT) unit 32, signal power computing unit 33 and Combined estimator unit 34, also comprise: window adding device 35 and zero padding element 36, wherein

The channel condition information of the Received signal strength that window adding device 35 pairs of channel estimating unit 31 estimate adds window function.

Wherein, add by empty subcarrier the power leakage caused in order to reduce, in the present embodiment, the channel estimation results of channel estimating unit 31 can be multiplied by certain window function.The present embodiment is not limited window function here, typically, can adopt Hanning window.

Zero padding element 36 carries out zero padding to the channel condition information that with the addition of window function, to obtain the input of inverse fast Fourier transform unit 32.

Wherein, zero padding is exactly the power side IFFT of formation 2, to reduce implementation complexity.In the present embodiment, the function of channel estimating unit 31, inverse fast Fourier transform (IFFT) unit 32, signal power computing unit 33 and Combined estimator unit 34 is identical respectively with the channel estimating unit 21 of embodiment 1, inverse fast Fourier transform (IFFT) unit 22, signal power computing unit 23 and Combined estimator unit 24, its content is incorporated in this, does not repeat them here.

By the multi-parameter Combined estimator device of the embodiment of the present invention to comprising time migration, multiple parameters of delay spread and noise power carry out Combined estimator, improve the demodulation performance of PDSCH, improve the performance of data throughout.Relative to embodiment 1, the use of window function can improve the estimated accuracy of each parameter.

Embodiment 3

The embodiment of the present invention additionally provides a kind of multi-parameter Combined estimator device.In the present embodiment, the composition of this multi-parameter Combined estimator device is identical with embodiment 1 or embodiment 2.In addition, in the present embodiment, what Combined estimator unit 24/34 was concrete utilizes signal power sampled value skew estimated time in the time window in time domain and postpones expansion, and utilizes the noise power sampled value estimating noise power in the noise power estimation window in time domain.It can comprise: time migration estimation module 41, delay spread estimation module 42 and noise power estimation module 43.Fig. 4 is the composition schematic diagram of this Combined estimator unit 24/34, please refer to Fig. 4.

In an execution mode of time migration estimation module 41, signal power sampled value moves together by this time migration estimation module 41 first application cycle displacement, then according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, the signal power sampled value of each the signal power sampled value selected and follow-up second predetermined quantity thereof is added, obtain the energy value of each signal power sampled value of corresponding selection, maximum energy value is selected from energy value obtained above, and determine time offset according to the time of advent of its correspondence.

Wherein, the second predetermined quantity is less than the first predetermined quantity, but the present embodiment does not limit the value of the first predetermined quantity and the second predetermined quantity, and it can preset, also can according to other strategy decisions, and the present embodiment is not in this, as restriction.

Wherein, if the quantity of its follow-up signal power sampled value is less than the second predetermined quantity, then using its follow-up all signal power sampled values as the signal power sampled value of described second predetermined quantity, computational methods are the same, do not repeat them here.

Wherein, determine the time of advent, can determine time offset accordingly, the mode of side-play amount of specifically really fixing time can adopt existing means, and the present embodiment is not in this, as restriction.

In another execution mode of time migration estimation module 41, signal power sampled value moves together by this time migration estimation module 41 first application cycle displacement, then maximum signal power sampled value is chosen, using the signal power of the signal power sampled value little pre-determined threshold more maximum than this as signal power threshold value, the signal power sampled value alternatively signal power sampled value of preset signals power threshold will be greater than, the time of advent according to the nearest candidate signal power samples value of described candidate signal power samples value middle distance initial synchronization time determines time offset.

Wherein, the present embodiment does not limit the size of predetermined threshold value, this in advance threshold value can preset, also can according to other strategy decisions, the present embodiment is not in this, as restriction.

Wherein, similar with last execution mode, determine the time of advent, can determine time offset accordingly, the mode of side-play amount of specifically really fixing time can adopt existing means, and the present embodiment is not in this, as restriction.

In an execution mode of delay spread estimation module 42, signal power sampled value moves together by this delay spread estimation module 42 first application cycle displacement, then according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, to the energy of each signal power sampled value selected and be weighted on average to the distance of up-to-date synchronization timing, obtain delay spread amount.

In the another one execution mode of delay spread estimation module 42, signal power sampled value moves together by this delay spread estimation module 42 first application cycle displacement, then maximum signal power sampled value is chosen, using the signal power of signal power sampled value little certain pre-determined threshold more maximum than this as signal power threshold value, to be greater than in the signal power sampled value of preset signals power threshold, apart from the time migration distance corresponding to up-to-date synchronization timing signal power sampled value farthest as delay spread amount.

In an execution mode of noise power estimation module 43, this noise power estimation module 43 is according to the original position of the noise power estimation window preset and final position, noise power sampled value within the scope of the original position of described noise power estimation window and final position is averaged, obtains noise power estimation result.

Wherein, noise power estimation window, referred to as noise window, in the present embodiment, pre-sets original position and the final position of the noise window of corresponding different system bandwidth.Noise power estimation module 43 is when estimating noise power, original position and the final position of the noise window of its correspondence is determined according to system bandwidth, again within the scope of this original position and final position, noise power sampled value is averaged, obtains noise power estimation result.

In the another one execution mode of noise power estimation module 43, this noise power estimation module 43 is according to the original position of delay spread value determination noise power estimation window, and determine the final position of noise power estimation window according to counting of inverse fast Fourier transform, noise power sampled value within the scope of the original position of described noise power estimation window and final position is averaged, obtains noise power estimation result.

Wherein, counting of inverse fast Fourier transform is different according to system bandwidth difference, and also can pre-set, the present embodiment is not in this, as restriction.

Wherein, the original position of noise window and final position are that the mode of counting represents, this point is same as the prior art, does not repeat them here.

In the present embodiment, be divided into by Combined estimator unit three modules so that the function of its multi-parameter Combined estimator to be described, but, the present embodiment is not in this, as restriction, such as, time migration estimation module 41 and delay spread estimation module 42 can merge execution, corresponding respective previous execution mode, can find out the signal power sampled value (peak value) of the first predetermined quantity, then estimated time offsets and delay spread accordingly; A corresponding respective rear execution mode, can find out signal power threshold value, and then estimated time offsets and delay spread accordingly.In addition, two execution modes of noise power estimation module 43 go for time migration estimation module 41 and delay spread estimation module 42 separates situation about performing, and also go for the situation of time migration estimation module 41 and delay spread estimation module 42 merging execution.

In the present embodiment; the function of time migration estimation module 41, delay spread estimation module 42 and noise power estimation module 43 is described respectively with two execution modes; but; the present embodiment is not in this, as restriction; any signal power sampled value that utilizes offsets and delay spread estimated time, and utilizes the mode of noise power sampled value estimating noise power to be all contained in the protection range of the present embodiment.

By the multi-parameter Combined estimator device of the present embodiment to comprising time migration, multiple parameters of delay spread and noise power carry out Combined estimator, improve the demodulation performance of PDSCH, especially improve the demodulation performance of the PDSCH of the R11UE being configured to QCL type B, and improve the performance of data throughout.

In order to the constitute and function of the multi-parameter Combined estimator device making embodiment 1-embodiment 3 is clearly understandable, 4 be described in detail in conjunction with the multi-parameter Combined estimator device of an instantiation to the present embodiment by the following examples.

Embodiment 4

Embodiments provide a kind of terminal equipment.Fig. 5 is the composition schematic diagram of this terminal equipment, and this terminal equipment contains multi-parameter Combined estimator device, and it can be realized by the multi-parameter Combined estimator device of embodiment 1-embodiment 3.Please refer to Fig. 5, this terminal equipment comprises: ZF channel estimating unit 51, window adding device 52, zero padding element 53, IFFT unit 54, Norm unit 55, TO estimate and compensating unit 56, delay spread estimation unit 57, noise power estimation unit 58, DM-RS channel estimating unit 59, NZP CSI-RS channel estimating unit 510 and CSI-RS copy generation unit 511.

Wherein, ZF channel estimating unit 51, window adding device 52, zero padding element 53, IFFT unit 54, Norm unit 55, TO estimates and compensating unit 56, delay spread estimation unit 57, noise power estimation unit 58 can respectively by the channel estimating unit 21/31 of previous embodiment, window adding device 35, zero padding element 36, IFFT unit 22/32, signal power computing unit 23/33, time migration estimation module 41, delay spread estimation module 42, noise power estimation module 43 realizes, concrete constitute and function is described in previous embodiment 1-embodiment 3, the no longer repeat specification of content something in common.

Please refer to Fig. 5, α ^(p)represent CSI-RS copy (p=15 ~ 22).Y is the Received signal strength on frequency domain, and it can be expressed as formula (1):

$Y(l,k)=H(l,k)\×S(l,k)\×e^{{-j2\mathrm{\π\Δtk}/N}_{\mathrm{fft}}}+n(l,k)---\left(1\right)$

For ZF channel estimating unit 51, as described in Example 1, it utilizes the local CSI-RS generated to carry out channel estimating to received signal, obtains the channel condition information of described Received signal strength.ZF channel estimating can be represented as formula (2):

$H_{\mathrm{ZF}}(a,b,n_s,m)=\frac{1}{2}\{Y(a,b,n_s,l,k)*{\mathrm{\α}}^{\left(b\right)}(n_s,l,k)+Y(a,b,n_s,l+1,k)*{\mathrm{\α}}^{\left(b\right)}(n_s,l+1,k)\}---\left(2\right)$

Wherein, a receives branch index (0≤a≤N _rx-1), the index of also i.e. reception antenna number; B sends CSI-RS port index (0≤b≤N _tx-1), n _sbe subframe index, l is OFDM(Orthogonal Frequency Division Multiplexing, OFDM) index of symbol, k is the index of subcarrier, and m is positioned at RB(Resource Block, Resource Block) CSI-RS unit (l in index, k), 0≤m≤N _rB-1.

For window adding device 52, as described in Example 2, its channel condition information to the Received signal strength that ZF channel estimating unit 51 estimates adds window function.Wherein, add window function and be also multiplied by certain window function by ZF estimated result, can reduce thus and add by empty subcarrier the power leakage caused.As formula (3):

$H_{\mathrm{ZF}}^{\mathrm{win}}(a,b,n_s,m)=H_{\mathrm{ZF}}(a,b,n_s,m)*\mathrm{win}\left(m\right)---\left(3\right)$

Wherein, window function can be Hanning window, also can be other window functions, the present embodiment not in this, as restriction, wherein, the expression formula of Hanning window as shown in Equation (4):

$\mathrm{win}\left(m\right)=0.5+0.5\cos \left(\frac{N_{\mathrm{CSIRS}}\mathrm{\&pi;}-2\mathrm{\&pi;m}}{N_{\mathrm{CSIRS}}}\right),m<\frac{N_{\mathrm{CSIRS}}}{2}$ （4）

$\mathrm{win}\left(m\right)=0.5+0.5\cos \left(\frac{N_{\mathrm{CSIRS}}\mathrm{\&pi;}+2\mathrm{\&pi;m}}{N_{\mathrm{CSIRS}}}\right),m\&GreaterEqual;\frac{N_{\mathrm{CSIRS}}}{2}$

For zero padding element 53, as described in Example 2, it carries out zero padding to the channel condition information that with the addition of window function, to obtain the input of inverse fast Fourier transform unit 54.Wherein, filling zero is exactly formation 2 power IFFT.The input signal of IFFT can be expressed as formula (5):

Wherein, the size N of IFFT _iFFTas shown in table 1, according to different system bandwidths, the size of IFFT is different.

Table 1

System bandwidth (System Bandwidth)	IFFT size (N IFFT）
		1.4M	64
3M	64
		5M	64
10M	64
		15M	128
20M	128

For IFFT unit 54, as described in Example 1, the signal after zero padding element 53 is filled is transformed to time domain from frequency domain by it.The output signal of IFFT unit can be passed through formula (6) and represent.

$H_{\mathrm{ifft}\_\mathrm{out}}(a,b,n_s,m)=\underset{k=0}{\overset{N_{\mathrm{IFFT}}-1}{\mathrm{\&Sigma;}}}{H}_{\mathrm{ifft}\_\mathrm{in}}(a,b,n_s,m)e^{{j2\mathrm{\&pi;mk}/N}_{\mathrm{IFFT}}}---\left(6\right)$

For signal power computing unit 55, as previously mentioned, it can calculate the signal power of the channel condition information in time domain according to the output of IFFT unit 54.Wherein, the signal power of the IFFT of output can be used to estimate timing offset, delay spread and noise power.This signal power can be passed through formula (7) and represent.

$H_{\mathrm{norm}}^2(a,n_s,m)=\frac{1}{N_{\mathrm{tx}}}\underset{a=0}{\overset{N_{\mathrm{tx}}-1}{\mathrm{\&Sigma;}}}{\left|H_{\mathrm{ifft}\_\mathrm{out}}(a,b,n_s,m)\right|}^2---\left(7\right)$

In the present embodiment, calculated the signal power of the channel condition information in time domain by signal power computing unit 55 after, the Combined estimator of multi-parameter can be implemented.

Fig. 6 illustrates the main thought of the multi-parameter Combined estimator that the present embodiment proposes, and it can be realized by the Combined estimator unit 24/34 described in embodiment 1 and embodiment 2.Concrete, in time domain, the signal power sampled value in T window is used to skew estimated time and delay spread.Then, based on the noise power sampled value in N window, calculate the noise power of CSI-RS channel estimating.

Fig. 7 illustrates the handling process of multi-parameter Combined estimator, please refer to Fig. 7, and this flow process comprises:

Step 701: signaling rate sampled value moves together by application cycle displacement.T ₁represent cyclic shift window length.

Step 702: utilize signal power sampled value skew estimated time in T window and delay spread.

Step 703: utilize the noise power sampled value estimating noise power in N window.

In a step 702, time migration is estimated, can be realized by time migration estimation and compensating unit 56, as described in Example 3, it utilizes signal power sampled value skew estimated time in the time window in time domain, in the present embodiment, propose two kinds of execution modes, but the present embodiment is not in this, as restriction.In addition, the time migration realized for this time migration estimation and compensating unit 56 compensates, and can adopt existing compensation method to realize, not repeat them here.

Fig. 8 illustrates the execution mode that time migration is estimated.As shown in Figure 8, first K(first predetermined quantity is detected) individual peak value path.Corresponding path indexing and signal power sampled value are expressed as t _peak(k) and V _peak(k), wherein, 0≤k≤K-1.Processing procedure is as follows, wherein, and T _maskrepresent mask length of window.

By above process, can determine the path at the peak that energy value is maximum, the time of advent of its correspondence can be expressed as formula (10).Wherein, each energy value corresponding to peak value path can be expressed as formula (9).Wherein, K _sumrepresent total length of window.

$t_d=\arg \underset{0\&le;k\&le;K}{\max }\left(G\left(k\right)\right)---\left(10\right)$

$\mathrm{\&Delta;t}=(t_d-T_1)*\frac{N_{\mathrm{FFT}\_\mathrm{sys}}}{12*N_{\mathrm{IFFT}}}---\left(14\right)$

The time of advent in the path at the peak having had energy maximum, time migration can be determined accordingly, concrete, formula (14) can be passed through and represent.

$\mathrm{\&Delta;t}=(t_d-T_1)*\frac{N_{\mathrm{FFT}\_\mathrm{sys}}}{12*N_{\mathrm{IFFT}}}---\left(14\right)$

Wherein, N _{fFT_sys}for the FFT size of whole system bandwidth.

Fig. 9 illustrates another execution mode that time migration is estimated.As shown in Figure 9, in this embodiment, a threshold value V is defined _tHfind the first footpath (the peak value path that energy is maximum).Specifically be described in embodiment 3.Can be represented by following formula (11)-(13).

${\mathrm{\&alpha;}}_{\mathrm{TH}}={10}^{V_{\mathrm{TH}}/10}---\left(11\right)$ $V_{\max }=\arg \underset{{0\&le;t\&le;T}_{\mathrm{window}}-1}{\max }\left(V(a,n_s,t)\right)---\left(12\right)$

t _d=the_first_path(V(a,n _s,t)≥α _TH*V _max) （13）

Same, utilize the time of advent in this first footpath, time migration can be determined according to formula (14).

After estimating time migration Δ t, can compensate accordingly, the data after compensation can be expressed as formula (15).

$Y_{\mathrm{comp}}(a,b,n_s,l,k)=\left\{\begin{array}{c}{Y}_{\mathrm{comp}}(a,b,n_s,l,k)*e^{j\frac{2\mathrm{\&pi;\&Delta;t}*(k-N_{\mathrm{FFT}\_\mathrm{sys}}/2)}{N_{\mathrm{FFT}}}},(0\&le;k\&le;\frac{N_{\mathrm{FFT}\_\mathrm{sys}}}{2}-1)\\ Y_{\mathrm{comp}}(a,b,n_s,l,k)*e^{j\frac{2\mathrm{\&pi;\&Delta;r}*(k-N_{\mathrm{FFT}\_\mathrm{sys}}/2+1)}{N_{\mathrm{FFT}\_\mathrm{sys}}}},(\frac{N_{\mathrm{FFT}\_\mathrm{sys}}}{2}\&le;k\&le;N_{\mathrm{FFT}\_\mathrm{sys}}-1)\end{array}\right.---\left(15\right)$

In a step 702, delay spread is estimated, can be realized by delay spread estimation unit 57, it is mainly used in DM-RS CHEST(Channel estimation, channel estimating) and CSI-RS CHEST execution MMSE(Minimum MSE, minimum MSE(Minimum Squared Error, Minimum Mean Square Error)) channel estimating.As described in Example 3, it is also utilize the signal power sampled value estimation time delay in the time window in time domain to expand, and in the present embodiment, propose two kinds of execution modes equally, but the present embodiment is not in this, as restriction.

Wherein execution mode can with reference to an explanation for Fig. 8 and embodiment 3, and Output rusults is the mean value in K peak value path, the estimated value of the delay spread also namely finally calculated, and is exactly the distance of K peak value and the weighted average of energy.By which, if use the time migration estimation method shown in Fig. 9, K peak value path also can be obtained based on this threshold value.Concrete, the estimated result of time delay expansion can pass through formula (16)-(18) and represent.

${\mathrm{\&tau;}}_{\mathrm{tmpl}}(a,n_s)=\frac{\underset{k=0}{\overset{K-1}{\mathrm{\&Sigma;}}}{V}_{\mathrm{peak}}\left(k\right)\mathrm{\&omega;}\left(k\right)*|t_{\mathrm{peak}}\left(k\right)-t_d|}{\underset{k=0}{\overset{K-1}{\mathrm{\&Sigma;}}}{V}_{\mathrm{peak}}\left(k\right)\mathrm{\&omega;}\left(k\right)}---\left(16\right)$

${\mathrm{\&tau;}}_{\mathrm{tmp}2}(a,n_s)=\frac{\underset{k=0}{\overset{K-1}{\mathrm{\&Sigma;}}}{V}_{\mathrm{peak}}\left(k\right)\mathrm{\&omega;}\left(k\right)*{|t_{\mathrm{peak}}\left(k\right)-t_d|}^2}{\underset{k=0}{\overset{K-1}{\mathrm{\&Sigma;}}}{V}_{\mathrm{peak}}\left(k\right)\mathrm{\&omega;}\left(k\right)}---\left(17\right)$

${\mathrm{\&tau;}}_{\mathrm{rms}}(a,n_s)={\mathrm{\&tau;}}_{\mathrm{tmp}2}(a,n_s)-{\mathrm{\&tau;}}_{\mathrm{tmp}1}^2(a,n_s)---\left(18\right)$

Wherein, ω (k) is the power adjustments factor.

Another kind of execution mode can with reference to the explanation of Fig. 9 and embodiment 3.Distance between this path, path and the 1st footpath are maximum, are considered to delay spread and export.Also namely, relative to up-to-date synchronization timing t _d, farthest that distance corresponding to energy sampled point is just considered to the estimated value of delay spread.Formula (19) can be expressed as.

${\mathrm{\&tau;}}_{\mathrm{rms}}(a,n_s)=\underset{0\&le;k\&le;K}{\max }\left\{\right|t_{\mathrm{peak}}\left(k\right)-t_d\left|\right\}---\left(19\right)$

In step 703, for noise power estimation, can be realized by noise power estimation unit 58, as described in Example 3, it utilizes the noise power sampled value estimating noise power in the noise power estimation window in time domain.

In the present embodiment, if obtain noise power estimation window N _window, noise power exports and can be represented as:

V _noise(a,n _s,k)=V _noise(a,n _s,(k+t _left)),(0≤k≤N _window-1) （20）

${\mathrm{\&sigma;}}^2(a,n_s)=\frac{\mathrm{CDM}\_\mathrm{GAIN}*\mathrm{\&beta;}}{N_{\mathrm{window}}}\underset{k=0}{\overset{N_{\mathrm{window}}-1}{\mathrm{\&Sigma;}}}{V}_{\mathrm{noise}}(a,n_s,k)---\left(21\right)$

N _window=t _right-t _left+1 （22）

$\mathrm{\&beta;}=\frac{N_{\mathrm{IFFT}}}{\underset{m=0}{\overset{N_{\mathrm{CSIRS}}}{\mathrm{\&Sigma;}}}{\left|\mathrm{win}\left(m\right)\right|}^2}---\left(23\right)$

${\mathrm{\&sigma;}}_{\mathrm{ave}}^2(a,n_s)=(1-{\mathrm{\&lambda;}}_{\mathrm{noise}}){\mathrm{\&sigma;}}_{\mathrm{ave}}^2(a,n_s-1)+{\mathrm{\&lambda;}}_{\mathrm{noise}}{\mathrm{\&sigma;}}^2(a,n_s)---\left(24\right)$

Wherein, t _leftand t _rightrepresent window edge.

Wherein, window edge can be determined by two kinds of execution modes, but the present embodiment is not in this, as restriction.

Wherein in a kind of execution mode, the edge of noise window can pass through formula (25) and represent.

$\left\{\begin{array}{c}{t}_{\mathrm{left}}=T_{\mathrm{left}}+t_d\\ t_{\mathrm{right}}=T_{\mathrm{right}}+t_d\end{array}\right.---\left(25\right)$

Wherein, T _leftand T _rightcan be as shown in table 2 by offline parameter file default settings.T _dit is the time of advent in the first footpath, as previously mentioned.

Table 2

In another execution mode, the edge of the estimated result determination noise window of delay spread can be utilized, specifically can pass through formula (26) and represent.

$\left\{\begin{array}{c}{t}_{\mathrm{left}}={\mathrm{\&tau;}}_{\mathrm{rms}}+T_{\mathrm{gap}}\\ t_{\mathrm{right}}=N_{\mathrm{IFFT}}-T_{\mathrm{gap}}\end{array}\right.---\left(26\right)$

Wherein, T _gapfor the protection interval of some, concrete value is determined on a case-by-case basis.

The terminal equipment provided by the embodiment of the present invention to comprising time migration, multiple parameters of delay spread and noise power carry out Combined estimator, improve the demodulation performance of PDSCH, improve the performance of data throughout.

Based on RAN4 test, the performance of the multi-parameter Combined estimator device that the embodiment of the present invention provides is significantly increased, and Figure 10 illustrates the demodulation performance of the multi-parameter Combined estimator device utilizing the embodiment of the present invention.As shown in Figure 10, in demodulation part, time migration to received signal is first had to compensate.Then, the delay spread estimated is used to carry out the PDSCH channel estimating based on DM-RS.Further, the noise power estimation result of multi-parameter Combined estimator is sent to CSI feedback part to calculate CQI.Figure 10 illustrates under CoMP scene, when the time migration of presence service community and PDSCH minizone, use the reception solution of the multi-parameter Combined estimator device of the embodiment of the present invention, as can be seen from Figure, the multi-parameter Combined estimator device that the embodiment of the present invention provides significantly improves throughput performance.

The embodiment of the present invention additionally provides a kind of multi-parameter combined estimation method, as as described in the following examples 5, the multi-parameter Combined estimator device of the principle of dealing with problems due to the method and embodiment 1-4 is similar, therefore its concrete enforcement can with reference to the enforcement of the device of embodiment 1-4, the no longer repeat specification of content something in common.

Embodiment 4

Embodiments provide a kind of multi-parameter combined estimation method.Figure 11 is the flow chart of the method, please refer to Figure 11, and the method comprises:

Step 1101: utilize the local CSI-RS generated to carry out channel estimating to received signal, obtain the channel condition information of described Received signal strength;

Step 1102: carry out inverse fast Fourier transform to described channel condition information, transforms to time domain by described channel condition information from frequency domain;

Step 1103: the signal power calculating the channel condition information in time domain;

Step 1104: utilize the signal power of the channel condition information in described time domain estimate described Received signal strength time migration, postpone expansion and noise power.

In one embodiment, before step 1102, described method can also comprise:

Step 1101a: window function is added to the channel condition information of described Received signal strength;

Step 1101b: carry out zero padding to the channel condition information that with the addition of window function, to carry out inverse fast Fourier transform to described channel condition information.

In step 1104, signal power sampled value skew estimated time in the time window in time domain can be utilized and postpone expansion, and utilizing the noise power sampled value estimating noise power in the noise power estimation window in time domain.Concrete, can be realized by the method described in Fig. 7, its content is incorporated in this, does not repeat them here.

Wherein, the step of signal power sampled value skew estimated time in the time window in time domain is utilized to have following two kinds of execution modes:

Execution mode one: signal power sampled value moves together by application cycle displacement, according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, the signal power sampled value of each the signal power sampled value selected and follow-up second predetermined quantity thereof is added, obtain the energy value of each signal power sampled value of corresponding selection, from energy value obtained above, select maximum energy value, and determine time offset according to the time of advent of its correspondence.

Execution mode two: signal power sampled value moves together by application cycle displacement, choose maximum signal power sampled value, using the signal power of the signal power sampled value little pre-determined threshold more maximum than this as signal power threshold value, to be greater than the signal power sampled value alternatively signal power sampled value of preset signals power threshold, the time of advent according to the nearest candidate signal power samples value of described candidate signal power samples value middle distance initial synchronization time determines time offset.

Wherein, the step utilizing the signal power sampled value estimation time delay in the time window in time domain to expand has following two kinds of execution modes:

Execution mode one: signal power sampled value moves together by application cycle displacement, according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, to the energy of each signal power sampled value selected and the time of advent side-play amount distance be weighted on average, obtain delay spread amount.

Execution mode two: signal power sampled value moves together by application cycle displacement, choose maximum signal power sampled value, using the signal power of signal power sampled value little certain pre-determined threshold more maximum than this as signal power threshold value, to be greater than in the signal power sampled value of preset signals power threshold, apart from the distance of up-to-date synchronization timing signal power sampled value farthest and described time offset as delay spread amount.

Wherein, the step of the noise power sampled value estimating noise power in the noise power estimation window in time domain is utilized to have following two kinds of execution modes:

Execution mode one: according to the original position of noise power estimation window preset and final position, is averaged to the noise power sampled value within the scope of the original position of described noise power estimation window and final position, obtains noise power estimation result.

Execution mode two: according to the original position of delay spread value determination noise power estimation window, and determine the final position of noise power estimation window according to counting of inverse fast Fourier transform, noise power sampled value within the scope of the original position of described noise power estimation window and final position is averaged, obtains noise power estimation result.

The method provided by the embodiment of the present invention to comprising time migration, multiple parameters of delay spread and noise power carry out Combined estimator, improve the demodulation performance of PDSCH, improve the performance of data throughout.

Apparatus and method more than the present invention can by hardware implementing, also can by combination of hardware software simulating.The present invention relates to such computer-readable program, when this program is performed by logical block, this logical block can be made to realize device mentioned above or component parts, or make this logical block realize various method mentioned above or step.The processor etc. used in logical block such as field programmable logic parts, microprocessor, computer.The invention still further relates to the storage medium for storing above program, as hard disk, disk, CD, DVD, flash memory etc.

More than in conjunction with concrete execution mode, invention has been described, but it will be apparent to those skilled in the art that these descriptions are all exemplary, is not limiting the scope of the invention.Those skilled in the art can make various variants and modifications according to spirit of the present invention and principle to the present invention, and these variants and modifications also within the scope of the invention.

Claims (10)

1. a multi-parameter Combined estimator device, wherein, described device comprises:

Channel estimating unit, it utilizes the local channel state information reference signals (CSI-RS) generated to carry out channel estimating to received signal, obtains the channel condition information of described Received signal strength;

Inverse fast Fourier transform unit, described channel condition information is transformed to time domain from frequency domain by it;

Signal power computing unit, it calculates the signal power of the channel condition information in time domain;

Combined estimator unit, it utilizes, and the signal power of the channel condition information in described time domain estimates the time migration of described Received signal strength, delay is expanded and noise power.

2. device according to claim 1, wherein, described device also comprises:

Window adding device, its channel condition information to the described Received signal strength that described channel estimating unit estimates adds window function;

Zero padding element, it carries out zero padding to the channel condition information that with the addition of window function, to obtain the input of described inverse fast Fourier transform unit.

3. device according to claim 1, wherein, described Combined estimator unit utilizes signal power sampled value skew estimated time in the time window in time domain and postpones expansion, and utilizes the noise power sampled value estimating noise power in the noise power estimation window in time domain.

4. device according to claim 3, wherein, described Combined estimator unit comprises:

Time migration estimation module, signal power sampled value moves together by the displacement of its application cycle, according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, the signal power sampled value of each the signal power sampled value selected and follow-up second predetermined quantity thereof is added, obtain the energy value of each signal power sampled value of corresponding selection, from energy value obtained above, select maximum energy value, and determine time offset according to the time of advent of its correspondence; Or, signal power sampled value moves together by the displacement of its application cycle, choose maximum signal power sampled value, using the signal power of the signal power sampled value little pre-determined threshold more maximum than this as signal power threshold value, to be greater than the signal power sampled value alternatively signal power sampled value of preset signals power threshold, the time of advent according to the nearest candidate signal power samples value of described candidate signal power samples value middle distance initial synchronization time determines time offset.

5. device according to claim 3, wherein, described Combined estimator unit comprises:

Delay spread estimation module, signal power sampled value moves together by the displacement of its application cycle, according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, to the energy of each signal power sampled value selected and the time of advent side-play amount distance be weighted on average, obtain delay spread amount; Or, signal power sampled value moves together by the displacement of its application cycle, choose maximum signal power sampled value, using the signal power of signal power sampled value little certain pre-determined threshold more maximum than this as signal power threshold value, to be greater than in the signal power sampled value of preset signals power threshold, apart from the distance of up-to-date synchronization timing signal power sampled value farthest and described time offset as delay spread amount.

6. device according to claim 3, wherein, described Combined estimator unit comprises:

Noise power estimation module, it is according to the original position of the noise power estimation window preset and final position, noise power sampled value within the scope of the original position of described noise power estimation window and final position is averaged, obtains noise power estimation result; Or, it is according to the original position of delay spread value determination noise power estimation window, and determine the final position of noise power estimation window according to counting of inverse fast Fourier transform, noise power sampled value within the scope of the original position of described noise power estimation window and final position is averaged, obtains noise power estimation result.

7. a multi-parameter combined estimation method, wherein, described method comprises:

Utilize the local CSI-RS generated to carry out channel estimating to received signal, obtain the channel condition information of described Received signal strength;

Inverse fast Fourier transform is carried out to described channel condition information, described channel condition information is transformed to time domain from frequency domain;

Calculate the signal power of the channel condition information in time domain;

Utilize the signal power of the channel condition information in described time domain estimate described Received signal strength time migration, postpone expansion and noise power.

8. method according to claim 7, wherein, described method also comprises:

Window function is added to the channel condition information of described Received signal strength;

Zero padding is carried out to the channel condition information that with the addition of window function, to carry out inverse fast Fourier transform to described channel condition information.

9. method according to claim 7, wherein, utilize the signal power of the channel condition information in described time domain estimate described Received signal strength time migration, postpone expansion and the step of noise power comprises:

Utilize signal power sampled value skew estimated time in the time window in time domain and postpone expansion, and utilizing the noise power sampled value estimating noise power in the noise power estimation window in time domain.

10. method according to claim 9, wherein,

The step of signal power sampled value skew estimated time in the time window in time domain is utilized to comprise: signal power sampled value moves together by application cycle displacement, according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, the signal power sampled value of each the signal power sampled value selected and follow-up second predetermined quantity thereof is added, obtain the energy value of each signal power sampled value of corresponding selection, maximum energy value is selected from energy value obtained above, and determine time offset according to the time of advent of its correspondence, or, signal power sampled value moves together by application cycle displacement, choose maximum signal power sampled value, using the signal power of the signal power sampled value little pre-determined threshold more maximum than this as signal power threshold value, to be greater than the signal power sampled value alternatively signal power sampled value of preset signals power threshold, the time of advent according to the nearest candidate signal power samples value of described candidate signal power samples value middle distance initial synchronization time determines time offset,

The step utilizing the signal power sampled value estimation time delay in the time window in time domain to expand comprises: signal power sampled value moves together by application cycle displacement, according to the signal power sampled value of signal power sampled value selective sequential first predetermined quantity from big to small, to the energy of each signal power sampled value selected and the time of advent side-play amount distance be weighted on average, obtain delay spread amount; Or, signal power sampled value moves together by application cycle displacement, choose maximum signal power sampled value, using the signal power of signal power sampled value little certain pre-determined threshold more maximum than this as signal power threshold value, to be greater than in the signal power sampled value of preset signals power threshold, apart from the distance of up-to-date synchronization timing signal power sampled value farthest and described time offset as delay spread amount;

The step of the noise power sampled value estimating noise power in the noise power estimation window in time domain is utilized to comprise: according to original position and the final position of the noise power estimation window preset, noise power sampled value within the scope of the original position of described noise power estimation window and final position is averaged, obtains noise power estimation result; Or, according to the original position of delay spread value determination noise power estimation window, and determine the final position of noise power estimation window according to counting of inverse fast Fourier transform, noise power sampled value within the scope of the original position of described noise power estimation window and final position is averaged, obtains noise power estimation result.