[go: up one dir, main page]

CN1269325C - Channel estimation method for WCDMA communication system - Google Patents

Channel estimation method for WCDMA communication system Download PDF

Info

Publication number
CN1269325C
CN1269325C CN 200310103987 CN200310103987A CN1269325C CN 1269325 C CN1269325 C CN 1269325C CN 200310103987 CN200310103987 CN 200310103987 CN 200310103987 A CN200310103987 A CN 200310103987A CN 1269325 C CN1269325 C CN 1269325C
Authority
CN
China
Prior art keywords
channel estimation
symbol
channel
time slot
estimation value
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Expired - Fee Related
Application number
CN 200310103987
Other languages
Chinese (zh)
Other versions
CN1543106A (en
Inventor
万里龙
丁杰伟
任震
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Haimen Science And Technology Development General Corp
Original Assignee
ZTE Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by ZTE Corp filed Critical ZTE Corp
Priority to CN 200310103987 priority Critical patent/CN1269325C/en
Publication of CN1543106A publication Critical patent/CN1543106A/en
Application granted granted Critical
Publication of CN1269325C publication Critical patent/CN1269325C/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Landscapes

  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种用于宽带码分多址系统的信道估计方法,首先对输入RAKE接收机的宽带扩频信号进行解扰解扩,再用扩频码对其解扩,分别得到专用物理控制信道和专用物理数据信道;然后计算初始信道估计值,对于连续的经过解相关运算的专用物理控制信道符号串,用算数平均的算法求出窗口对应时刻的信道估计值Cem;最后将Cem补偿窗口对应时刻m时刻的专用物理数据信道符号,各个多径补偿后的结果输出给最大比合并模块。本发明所述方法克服了现有技术存在的在信道恶化的情况下信道估计值无法反映信道特性,从而导致的误判率增加、性能变差的缺点,能够有效提高在移动通信环境尤其是高速移动通信环境下的信道估计准确性,从而保证系统高性能。

The invention discloses a channel estimation method for a wideband code division multiple access system. First, descramble and despread the wideband spread spectrum signal input to a RAKE receiver, and then despread it with a spread spectrum code to obtain special physical Control channel and dedicated physical data channel; then calculate the initial channel estimation value, for the continuous dedicated physical control channel symbol string after decorrelation operation, use the algorithm of arithmetic mean to find the channel estimation value Cem at the corresponding moment of the window; finally, Cem compensation The window corresponds to the dedicated physical data channel symbol at time m, and each multipath compensated result is output to the maximum ratio combining module. The method of the present invention overcomes the disadvantages in the prior art that the channel estimation value cannot reflect the channel characteristics when the channel deteriorates, which leads to an increase in the misjudgment rate and poor performance, and can effectively improve the performance in the mobile communication environment, especially at high speeds. The channel estimation accuracy in the mobile communication environment, thus ensuring the high performance of the system.

Description

用于宽带码分多址通信系统的信道估计方法Channel Estimation Method for Wideband Code Division Multiple Access Communication System

技术领域technical field

本发明涉及无线通信系统的信道估计方法,特别是涉及一种宽带码分多址(WCDMA,Wide CDMA)通信系统的信道估计方法。The present invention relates to a channel estimation method of a wireless communication system, in particular to a channel estimation method of a Wideband Code Division Multiple Access (WCDMA, Wide CDMA) communication system.

背景技术Background technique

随着移动通信的迅猛发展,人们对移动通信的质量及其提供的业务类型要求也越来越高,第三代移动通信的提出及发展正是符合了人们的这种需求,因此备受瞩目。宽带码分多址(WCDMA,Wide CDMA)是第三代移动通信系统标准化组织(3rd Generation PartnershipProject,3GPP)提出的无线传输技术(Radio TransmissionTechnology,RTT)方案,它不但能提供高质量的话音服务,而且能够提供与固定网络类似的多媒体业务。With the rapid development of mobile communication, people have higher and higher requirements for the quality of mobile communication and the types of services it provides. The proposal and development of the third generation of mobile communication is exactly in line with this demand of people, so it has attracted much attention. . Wideband Code Division Multiple Access (WCDMA, Wide CDMA) is a Radio Transmission Technology (RTT) solution proposed by the 3rd Generation Partnership Project (3GPP), which can not only provide high-quality voice services , and can provide multimedia services similar to fixed networks.

如图1所示,在WCDMA系统中,为了提高上行链路质量,在上行专用物理控制信道(Dedicated Physical Control Channel,DPCCH)中加入了部分导频符号,用于信道估计以进行相干检测。导频符号的发射是不连续的,根据不同的时隙格式,每个时隙内的导频符号数介于3到8个之间,因此,采用何种方法来精确估计时隙内导频符号段和数据符号段对应的信道特性,对提高接收机在高速移动环境下的接收性能显得极其重要。As shown in Figure 1, in a WCDMA system, in order to improve the uplink quality, some pilot symbols are added to the uplink Dedicated Physical Control Channel (DPCCH) for channel estimation and coherent detection. The transmission of pilot symbols is discontinuous. According to different slot formats, the number of pilot symbols in each slot is between 3 and 8. Therefore, what method is used to accurately estimate the pilot in a slot The channel characteristics corresponding to the symbol segment and the data symbol segment are extremely important to improve the receiving performance of the receiver in a high-speed mobile environment.

目前较为常用的信道估计方法主要有以下几种:单时隙平均的信道估计方法、基于简化卡尔曼跟踪滤波的信道估计方法、多时隙加权平均的信道估计方法(Weighted Multi-Slot Averaging,WMSA)和迭代信道估计法(Iterative Channel Estimation,ICE)。At present, the more commonly used channel estimation methods mainly include the following: single-slot average channel estimation method, channel estimation method based on simplified Kalman tracking filter, and multi-slot weighted average channel estimation method (Weighted Multi-Slot Averaging, WMSA) And Iterative Channel Estimation (ICE).

单时隙平均的信道估计方法由于导频符号较少,无法有效滤除噪声的影响,因此不论是在低速还是高速移动环境下性能都较差。基于简化卡尔曼跟踪滤波的信道估计方法是前者的改进,在低速移动环境下性能有所改善,但在高速移动环境下性能仍然很差。The single-slot average channel estimation method has less pilot symbols and cannot effectively filter out the influence of noise, so its performance is poor no matter in low-speed or high-speed mobile environments. The channel estimation method based on simplified Kalman tracking filter is an improvement of the former, and its performance is improved in low-speed mobile environment, but its performance is still poor in high-speed mobile environment.

多时隙加权平均的信道估计方法(WMSA,Weighted Multi-SlotAveraging)的估计过程如下:1)利用单时隙平均的信道估计方法得到各时隙的信道估计值;2)将多个时隙的估计值加权组合,得到当前时隙的信道估计值。由于在多个时隙间进行加权平均,有效滤除了噪声的影响,因而在信号慢衰落的情况下,该方法能得到良好的信道估计性能,但随着信道的恶化,如移动速度加快,由于估计值不能很好地反映时隙内不同时段的信道特性,因而性能变差。The estimation process of the multi-slot weighted average channel estimation method (WMSA, Weighted Multi-SlotAveraging) is as follows: 1) use the single-slot average channel estimation method to obtain the channel estimation value of each time slot; Values are weighted and combined to obtain the channel estimation value of the current time slot. Since the weighted average is performed between multiple time slots, the influence of noise is effectively filtered, so in the case of slow signal fading, this method can obtain good channel estimation performance, but as the channel deteriorates, such as the speed of movement increases, due to The estimated value cannot well reflect the channel characteristics at different periods in the time slot, so the performance deteriorates.

迭代信道估计(ICE)方法在目前也有着较为普遍的应用,相应的专利有公开号为1384627的《一种用于宽带码分多址系统的信道估计方法》,其思路是:先通过多时隙加权平均(WMSA)得到各径当前时隙的信道估计值并对专用物理控制信道(DPCCH,Dedicated PhysicalControl Channel)进行补偿、合并,对合并后的结果进行硬判决,取出判决后的数据段符号作为虚拟导频并反馈回去对数据段符号重新进行信道估计,为提高精度,该过程可迭代多次。因为它以多时隙加权平均(WMSA)方法为基础,对数据段的信道估计、补偿仍然仅仅用到导频部分的信息,随着信道的恶化,多时隙加权平均(WMSA)方法的性能变差,合并后硬判决的误判率增加,这必将导致该方法的性能下降。The iterative channel estimation (ICE) method is also widely used at present. The corresponding patent has the publication number 1384627 "A Channel Estimation Method for Wideband Code Division Multiple Access System". The weighted average (WMSA) obtains the channel estimation value of the current time slot of each path, compensates and combines the dedicated physical control channel (DPCCH, Dedicated Physical Control Channel), performs hard judgment on the combined result, and takes out the data segment symbol after the judgment as The virtual pilot is fed back to re-estimate the channel for the data segment symbols. In order to improve the accuracy, this process can be iterated multiple times. Because it is based on the multi-slot weighted average (WMSA) method, the channel estimation and compensation of the data segment still only use the information of the pilot part. As the channel deteriorates, the performance of the multi-slot weighted average (WMSA) method becomes worse. , the misjudgment rate of hard decisions increases after merging, which will definitely lead to a decrease in the performance of the method.

发明内容Contents of the invention

本发明的目的就是克服现有技术存在的在信道恶化的情况下信道估计值无法反映信道特性,从而导致的误判率增加、性能变差的缺点,以期提出一种有效提高在移动通信环境尤其是高速移动通信环境下的信道估计准确性,从而保证高性能的信道估计方法。The purpose of the present invention is to overcome the shortcomings of the prior art that the channel estimation value cannot reflect the channel characteristics in the case of channel deterioration, resulting in an increase in misjudgment rate and poor performance, in order to propose a method that can effectively improve performance in mobile communication environments, especially in mobile communication environments. It is the channel estimation accuracy in the high-speed mobile communication environment, thereby ensuring a high-performance channel estimation method.

为实现上述目的,本发明提出一种用于宽带码分多址系统的信道估计方法,其特征在于,包括以下步骤:To achieve the above object, the present invention proposes a channel estimation method for wideband code division multiple access system, it is characterized in that, comprises the following steps:

第一步:对输入RAKE接收机的宽带扩频信号进行解扰解扩;Step 1: Descrambling and despreading the broadband spread spectrum signal input to the RAKE receiver;

第二步:根据已知的时隙格式从专用物理控制信道(DPCCH,Dedicated Physical Control Channel)的当前和未来一个时隙(以下简称“下一时隙”)内取出导频符号;对提取出的导频符号,每个导频符号与已知导频符号进行解相关运算,得到各导频符号的初始信道估计值;对每个时隙内所有导频符号的初始信道估计值进行算术平均,得到各时隙的初始信道估计值;The second step: extract the pilot symbols from the current and future time slot (hereinafter referred to as "the next time slot") of the dedicated physical control channel (DPCCH, Dedicated Physical Control Channel) according to the known time slot format; Pilot symbols, each pilot symbol is decorrelated with known pilot symbols to obtain the initial channel estimation value of each pilot symbol; the initial channel estimation value of all pilot symbols in each time slot is arithmetically averaged, Obtain the initial channel estimation value of each time slot;

第三步:对于导频部分,采用一个移动窗口,设窗口长度为L,窗口覆盖当前时隙中的当前导频符号的信道估计值以及历史L-1个符号的信道估计值,其中历史符号的信道估计值中既包括用理想导频解相关得到的导频符号的信道估计值,也包括非导频符号的信道估计值(求非导频符号的信道估计值的方法见步骤五)。窗口中的平均值为窗口中间位置对应时刻或者某个合适的时刻的信道估计值,该值也是信道估计模块最终输出并用来补偿DPDCH(Dedicated PhysicalData Channel)数据的信道估计值。这里的某个合适的时刻是指窗口中间偏向当前时刻的某个位置,主要是移动速度较慢时会设置一个较大的L值,这样会需要更多的数据存储器,为兼顾硬件资源,可以这种考虑一个合适时刻。这里和现有技术最大的不同的就是窗口中既包含了导频信息也包含非导频信息。有效地提高了信道估计的性能。对于非导频符号,首先用前一个窗口的结果和下一时隙的信道估计值的平均值作为当前符号的信道估计值,对当前符号进行补偿,然后对所有多径当前符号补偿后的值做最大比合并(MRC,Maximal RatioCombining),对合并后的结果再做硬判决,这样就得到当前非导频符号的虚拟的理想非导频(所谓虚拟是指本来我们并不预先知道,但我们根据其它已知信息,可信度较大地估计出它的发射端发送时的值)。用该虚拟值对当前接收非导频符号做解相关运算就得到当前非导频符号的信道估计值。最后用步骤四相同的方法,求出当前窗口中间位置对应时刻或者某个合适时刻位置的信道估计值。在这里,同样和现有技术有重要的不同,本方法不仅用上了未来时隙、当前时隙甚至历史时隙的导频部分信息,还用到了当前窗口中包含的已经估计出来的非导频部分的信息,进一步有效地提高了信道估计的性能。Step 3: For the pilot part, use a moving window, set the window length to L, and the window covers the channel estimation value of the current pilot symbol in the current time slot and the channel estimation value of the historical L-1 symbols, where the historical symbols The channel estimation value includes not only the channel estimation value of pilot symbols obtained by ideal pilot decorrelation, but also the channel estimation value of non-pilot symbols (see step 5 for the method of calculating the channel estimation value of non-pilot symbols). The average value in the window is the channel estimation value at the corresponding moment in the middle of the window or at a certain suitable moment, and this value is also the channel estimation value finally output by the channel estimation module and used to compensate DPDCH (Dedicated PhysicalData Channel) data. A suitable moment here refers to a certain position in the middle of the window that is biased towards the current moment. The main reason is that when the moving speed is slow, a larger L value will be set, which will require more data memory. In order to take into account hardware resources, you can This consideration is an opportune moment. The biggest difference between this and the prior art is that the window contains both pilot information and non-pilot information. Effectively improve the performance of channel estimation. For non-pilot symbols, first use the result of the previous window and the average value of the channel estimation value of the next time slot as the channel estimation value of the current symbol to compensate the current symbol, and then do Maximal RatioCombining (MRC, Maximal RatioCombining), make a hard decision on the combined result, so that the virtual ideal non-pilot of the current non-pilot symbol is obtained (the so-called virtual means that we do not know it in advance, but we based on For other known information, the value when it is sent by the transmitting end is estimated with a high degree of reliability). The channel estimation value of the current non-pilot symbol is obtained by performing decorrelation operation on the currently received non-pilot symbol with the dummy value. Finally, use the same method as Step 4 to obtain the channel estimation value at the time corresponding to the middle position of the current window or at an appropriate time. Here, there is also an important difference from the prior art. This method not only uses the pilot part information of future time slots, current time slots, and even historical time slots, but also uses the estimated non-guided frequency information contained in the current window. The information of the frequency part further effectively improves the performance of channel estimation.

第四步:用窗口输出的信道估计值对窗口中间位置时刻或者某个合适时刻的DPDCH(Dedicated Physical Data Channel)做补偿运算,补偿后的数据按照最大比合并的算法合并所有多径,得到的结果送给译码模块。Step 4: Use the channel estimation value output by the window to perform compensation calculation on the DPDCH (Dedicated Physical Data Channel) at the middle position of the window or at a certain suitable time. The compensated data is combined with all multipaths according to the maximum ratio combination algorithm, and the obtained The result is sent to the decoding module.

本发明所述的信道估计方法,与现有技术相比,多次充分利用了已经估计出来了的非导频的信息,在有限的移动窗口中充分地利用了所有可以利用的信息,有效地、明显地提高了信道估计模块的性能。采用发明阐述的信道估计方法,可以比现有技术更精确地得到一个时隙内每个符号位置的信道估计值,能在信道衰落慢和快时都能反映出一个时隙内不同时刻的信道特性变化情况,因而能够大大提高通信系统的接收性能。仿真结果表明,不管在移动速度慢还是快的时候,本发明阐述的信道估计方法的性能都要比现有技术的信道估计性能有明显的提高。Compared with the prior art, the channel estimation method of the present invention makes full use of the estimated non-pilot information many times, fully utilizes all available information in a limited moving window, and effectively , significantly improving the performance of the channel estimation module. Using the channel estimation method described in the invention, the channel estimation value of each symbol position in a time slot can be obtained more accurately than the prior art, and can reflect the channel at different times in a time slot when the channel fading is slow or fast Therefore, the receiving performance of the communication system can be greatly improved. Simulation results show that no matter when the moving speed is slow or fast, the performance of the channel estimation method described in the present invention is obviously improved compared with the channel estimation performance of the prior art.

附图说明Description of drawings

图1是WCDMA系统上行专用物理信道(DPCH,Dedicated PhysicalChannel)的帧结构;Fig. 1 is the frame structure of WCDMA system uplink dedicated physical channel (DPCH, Dedicated PhysicalChannel);

图2是现有的多时隙加权平均信道估计方法示意图。Fig. 2 is a schematic diagram of an existing multi-slot weighted average channel estimation method.

图3是本发明所述信道估计方法实施框图。Fig. 3 is a block diagram of implementing the channel estimation method of the present invention.

图4是本发明所述信道估计方法流程图。Fig. 4 is a flow chart of the channel estimation method of the present invention.

图5是本发明所述信道估计方法中迭代算法流程图。Fig. 5 is a flowchart of an iterative algorithm in the channel estimation method of the present invention.

图6是本发明所述信道估计方法中移动窗口算数平均方法示意图。Fig. 6 is a schematic diagram of the moving window arithmetic mean method in the channel estimation method of the present invention.

图7是本发明所述信道估计算法中求非导频符号信道估计值方法示意图。Fig. 7 is a schematic diagram of a method for calculating channel estimation values of non-pilot symbols in the channel estimation algorithm of the present invention.

具体实施方式Detailed ways

下面结合附图对技术方案的实施作进一步的详细描述,根据这些附图,同一领域的技术人员可以很容易实现本发明。The implementation of the technical solution will be further described in detail below in conjunction with the accompanying drawings. According to these drawings, those skilled in the art can easily realize the present invention.

本发明主要是提出了一种用于WCDMA系统上行信道的、基于充分利用导频和非导频信息的移动窗口算术平均的信道估计方法,以提高通信系统在移动环境下的接收性能。The present invention mainly proposes a channel estimation method for WCDMA system uplink channel based on moving window arithmetic mean which fully utilizes pilot frequency and non-pilot frequency information, so as to improve the receiving performance of the communication system in the mobile environment.

图1是WCDMA系统上行专用物理信道(DPCH)的帧结构;为了区别DPDCH中的数据域和DPCCH中的数据域,特别标识DPCCH中的数据域为非导频域。Fig. 1 is the frame structure of the WCDMA system uplink dedicated physical channel (DPCH); in order to distinguish the data field in the DPDCH and the data field in the DPCCH, the data field in the DPCCH is specially identified as a non-pilot field.

图2是现有技术中的多时隙加权平均信道估计算法的示意图。在该具体方法中,要采用多时隙加权平均信道估计方法来计算时隙2的信道估计值,具体过程是:先通过对时隙内各导频符号信道估计值进行算术平均分别得到时隙1、2、3和4的信道估计值,然后按图中各时隙对应的权值在相邻四个时隙间进行加权平均,得到时隙2最终的信道估计值,即Fig. 2 is a schematic diagram of a multi-slot weighted average channel estimation algorithm in the prior art. In this specific method, the multi-slot weighted average channel estimation method is used to calculate the channel estimation value of time slot 2. The specific process is: firstly, the channel estimation value of each pilot symbol in the time slot is arithmetically averaged to obtain the channel estimation value of time slot 1 respectively. , 2, 3 and 4 channel estimates, and then perform weighted average among adjacent four time slots according to the weights corresponding to each time slot in the figure, to obtain the final channel estimate of time slot 2, namely

EE. ^^ 22 == (( αα 22 EE. 11 ++ αα 11 EE. 22 ++ αα 11 EE. 33 ++ αα 22 EE. 44 )) // (( 22 αα 11 ++ 22 αα 22 ))

式中

Figure C20031010398700102
是时隙2的最终信道估计值,Ei,αi分别是各时隙内算术平均得到的信道估计值和进行多时隙加权平均时各时隙对应的权值。现有技术中的迭代信道估计法,就是用
Figure C20031010398700103
去补偿时隙2中的非导频数据,然后经过最大比合并、硬判决得到虚拟的已知非导频数据,再与解扩后的DPCCH中的非导频符号做解相关运算,得到非导频符号的信道估计值,从而得到整个时隙2的所有符号的信道估计值,再经过多时隙加权平均得到各个时隙的信道估计,该方法比普通的多时隙加权平均方法性能略好。In the formula
Figure C20031010398700102
is the final channel estimation value of time slot 2, E i , α i are the channel estimation value obtained by the arithmetic mean in each time slot and the weight value corresponding to each time slot when multi-slot weighted averaging is performed, respectively. The iterative channel estimation method in the prior art is to use
Figure C20031010398700103
To compensate the non-pilot data in time slot 2, and then obtain the virtual known non-pilot data through maximum ratio combination and hard decision, and then perform decorrelation operation with the non-pilot symbols in the despread DPCCH to obtain the non-pilot data The channel estimation value of the pilot symbols is used to obtain the channel estimation value of all symbols in the entire time slot 2, and then the channel estimation of each time slot is obtained through multi-slot weighted average. This method has slightly better performance than the ordinary multi-slot weighted average method.

图3是本发明所述信道估计方法实施框图,该图说明了要实施的本发明所述信道估计方法的模块在整个WCDMA RAKE接收机中的位置。在该实例中,AD转换器301将输入到接收机的模拟信号转变为数字信号送入解扰解扩模块304,解扰解扩模块304利用码跟踪器302和码生成器303产生的各条多径扰码对输入信号进行解扰解扩,并将解扰解扩后的各多径信号送入其对应的信道估计和补偿模块305,信道估计和补偿模块305分别计算出DPCCH每个符号对应时刻的信道估计值并对专用信道DPDCH进行补偿,补偿后的各径信号输出到最大比合并模块306进行最大比合并,最后将合并后的信号输入到解交织/译码器307内进行解交织和译码,得到发射端发射的原始信号。本发明提出的信道估计方法在实施框图中的位置是信道估计和补偿模块305。Fig. 3 is the implementation block diagram of the channel estimation method of the present invention, which has illustrated the position of the modules of the channel estimation method of the present invention to be implemented in the whole WCDMA RAKE receiver. In this example, the AD converter 301 converts the analog signal input to the receiver into a digital signal and sends it to the descrambling and despreading module 304. The multipath scrambling code descrambles and despreads the input signal, and sends each multipath signal after descrambling and despreading to its corresponding channel estimation and compensation module 305, and the channel estimation and compensation module 305 calculates each symbol of DPCCH respectively Compensate the dedicated channel DPDCH for the channel estimation value at the corresponding time, output the compensated signals of each path to the maximum ratio combination module 306 for maximum ratio combination, and finally input the combined signal into the deinterleaving/decoder 307 for decoding Interleave and decode to obtain the original signal transmitted by the transmitter. The position of the channel estimation method proposed in the present invention in the implementation block diagram is the channel estimation and compensation module 305 .

图4是本发明所述信道估计方法流程图,该图说明了要实施的本发明所述信道估计方法之前应该所做的准备工作,以及做完后处理结果的去向。Fig. 4 is a flow chart of the channel estimation method of the present invention, which illustrates the preparatory work that should be done before implementing the channel estimation method of the present invention, and the whereabouts of the processing results after completion.

这个逻辑流程始于步骤401,终止于步骤404。This logic flow starts at step 401 and ends at step 404 .

步骤401对输入RAKE接收机的宽带扩频信号进行解扰解扩。因为RAKE接收机接收到的信号都是经过加扰加扩的,因此需要用响应的扰码先对其解扰,即扰码和接收信号共轭相乘,在用扩频码对其解扩,分别得到DPCCH(Dedicated Physical Control Channel,专用物理控制信道)和DPDCH(Dedicated Physical Data Channel专用物理数据信道),这是后面步骤的准备工作,也只有解扰解扩后的信号才能够进一步作信道估计处理。Step 401 performs descrambling and despreading on the broadband spread spectrum signal input to the RAKE receiver. Because the signal received by the RAKE receiver is scrambled and spread, it needs to be descrambled first with the corresponding scrambling code, that is, the scrambling code is multiplied by the conjugate of the received signal, and then it is despread with the spreading code , to obtain DPCCH (Dedicated Physical Control Channel, dedicated physical control channel) and DPDCH (Dedicated Physical Data Channel dedicated physical data channel), which is the preparation work for the following steps, and only the signal after descrambling and despreading can be used for further channeling estimated processing.

顺序进入步骤402,计算初始信道估计值。根据已知的时隙格式从专用物理控制信道(DPCCH,Dedicated Physical Control Channel)的当前和未来一个时隙(以下简称“下一时隙”)内取出导频符号;对提取出的导频符号,每个导频符号与已知导频符号进行解相关运算,得到各导频符号的初始信道估计值;对每个时隙内所有导频符号的初始信道估计值进行算术平均,得到各时隙的初始信道估计值。Enter step 402 in order to calculate the initial channel estimation value. Extract the pilot symbols from the current and future time slot (hereinafter referred to as "the next time slot") of the dedicated physical control channel (DPCCH, Dedicated Physical Control Channel) according to the known time slot format; for the extracted pilot symbols, Each pilot symbol is decorrelated with the known pilot symbols to obtain the initial channel estimation value of each pilot symbol; the arithmetic average of the initial channel estimation values of all pilot symbols in each time slot is obtained to obtain the channel estimation value of each time slot The initial channel estimate of .

顺序进入步骤403,对于连续的经过解相关运算的DPCCH符号串,逐个逐个符号地移动窗口,用算数平均的算法求出窗口对应时刻的信道估计值CEm,这一步是本方法的核心部分,在图5中将做非常具体详尽的描述。Enter step 403 sequentially, for the continuous DPCCH symbol strings that have undergone decorrelation operations, move the window symbol by symbol one by one, and use the algorithm of arithmetic mean to obtain the channel estimation value CEm at the corresponding moment of the window. This step is the core part of this method. Figure 5 will do a very specific and detailed description.

顺序进入步骤404,用403得到的CEm补偿窗口对应时刻m时刻的DPDCH符号。各个多径补偿后的结果输出给最大比合并模块。Enter step 404 in sequence, use the CEm compensation window obtained in step 403 to correspond to the DPDCH symbol at time m. The results after each multipath compensation are output to the maximum ratio combining module.

图5是本发明所述信道估计方法中迭代算法流程图,该图详细解说了本发明所述的信道估计方法的计算流程。该图是本发明涉及的信道估计方法的核心部分,即如何既利用导频又利用非导频信息通过移动窗口求最精确的信道估计方法,图5是该方法的详细流程图,是步骤403的细化。该逻辑流程始于步骤501,终止于步骤512。Fig. 5 is a flowchart of an iterative algorithm in the channel estimation method of the present invention, which illustrates in detail the calculation flow of the channel estimation method of the present invention. This figure is the core part of the channel estimation method involved in the present invention, that is, how to use both pilot and non-pilot information to obtain the most accurate channel estimation method through a moving window. Figure 5 is a detailed flow chart of the method, which is step 403 refinement. The logic flow starts at step 501 and ends at step 512 .

步骤501,从步骤402得到的解扰解扩后的当前时隙和下一时隙的所有符号中,从当前时隙第一个符号开始,到当前时隙最后一个符号为止,一个符号一个符号地送给步骤502,一个经历图5的整个流程后再输入下一个符号。Step 501, from all the symbols of the descrambled and despreaded current slot and the next slot obtained in step 402, from the first symbol of the current slot to the last symbol of the current slot, symbol by symbol Sent to step 502, one enters the next symbol after going through the whole process of Fig. 5 .

步骤502,根据已知时隙格式(Slot Format)判别步骤501输送来的当前符号Sn是否为当前时隙中的第一个非导频符号,如果是步进到步骤506,如果否步进到503。Step 502, judge whether the current symbol Sn delivered by step 501 is the first non-pilot symbol in the current time slot according to the known time slot format (Slot Format), if it is stepping to step 506, if not stepping to 503.

步骤503,根据已知时隙格式判别当前符号Sn是否为导频符号,如果是步进到步骤505.1,如果否步进到步骤507。Step 503, judge whether the current symbol Sn is a pilot symbol according to the known time slot format, if yes, go to step 505.1, if not, go to step 507.

步骤504,根据已知时隙格式生成导频序列,供步骤505.1使用。Step 504, generate a pilot sequence according to the known time slot format for use in step 505.1.

步骤505,将接收到的符号与已知符号(或者估计出来的符号)做解相关运算,解相关运算实际就相乘,即去掉符号。其中步骤505.1的已知符号为步骤504根据已知时隙格式本地产生的导频符号序列中对应当前位置的符号。而步骤505.2中的已知符号,其实是经过步骤507、508、509、510估计出来,因为一般情况下,DPCCH的物理信道BER比较下,估计出来的符号可以近似当做已知符号应用于解相关运算。步骤505运算完后步进到步骤511。In step 505, a decorrelation operation is performed on the received symbol and the known symbol (or the estimated symbol), and the decorrelation operation is actually multiplication, that is, the symbol is removed. The known symbol in step 505.1 is the symbol corresponding to the current position in the pilot symbol sequence locally generated in step 504 according to the known time slot format. The known symbols in step 505.2 are actually estimated through steps 507, 508, 509, and 510, because in general, under the comparison of the physical channel BER of DPCCH, the estimated symbols can be approximately used as known symbols for decorrelation operation. Step 511 after step 505 operation.

步骤506,当步骤502接收到第一个非导频符号的时候就步进到本步骤,本步骤从402的输出结果中获取下一时隙的所有导频符号,再用根据已知时隙格式得到已知导频序列,对下一时隙的导频符号一一做解相关运算,然后将解相关运算后的结果求平均值就得到CEnext-slot。最后将CEnext-slot保存起来,因为估计所有的非导频符号时都将用到它。Step 506, just step into this step when step 502 receives the first non-pilot symbol, this step obtains all pilot symbols of the next time slot from the output result of 402, and then uses according to known time slot format The known pilot sequence is obtained, and the decorrelation operation is performed on the pilot symbols of the next time slot one by one, and then the CE next-slot is obtained by averaging the results after the decorrelation operation. Finally, save the CEnext-slot, because it will be used when estimating all non-pilot symbols.

步骤507,如图7所示,由CEm-1和CEnext-slot简单求平均值得到当前非导频符号的预信道估计值C’n,即C’n=(CEn-1+CEnext-slot)/2。其中CEm-1是由前一个符号Sn-1经过图5整个流程得到的结果,而CEnext-slot是在步骤506求得的结果。图7是求解当前非导频符号的预信道估计值的示意图,是步骤507的图形阐述表达方式。C’n将用来估计Sn在发射端的原符号。Step 507, as shown in Figure 7, simply calculate the average value of CE m-1 and CE next-slot to obtain the pre-channel estimated value C'n of the current non-pilot symbol, that is, C'n=(CE n-1 +CE next-slot )/2. Among them, CE m-1 is the result obtained from the previous symbol S n-1 through the entire process in FIG. 5 , and CE next-slot is the result obtained in step 506 . FIG. 7 is a schematic diagram of solving the pre-channel estimation value of the current non-pilot symbol, which is a graphical illustration of step 507. C'n will be used to estimate the original symbol of Sn at the transmitter.

步骤508,用C’n补偿Sn可以到补偿后的当前非导频符号S’n。Step 508, use C'n to compensate Sn to obtain the compensated current non-pilot symbol S'n.

步骤509,将所有多径中的S’n以最大比合并的方法合并得到合并后的符号S’n-source。Step 509, combining the S'n in all multipaths by the maximum ratio combination method to obtain the combined symbol S'n-source.

步骤510,对S’n-source做硬判决, S ^ n - source = ( S ' n - soucrce ≥ 0 ) ? 1 : - 1 , 那么

Figure C20031010398700133
就是发射端原符号Sn-source的估计值。Step 510, making a hard decision on S'n-source, S ^ no - source = ( S ' no - source &Greater Equal; 0 ) ? 1 : - 1 , So
Figure C20031010398700133
is the estimated value of the original symbol Sn-source at the transmitter.

步骤511,根据已知时隙格式判别当前符号Sn是否为导频符号,如果是选择开关指向左边接收505.1的输出结果,如果否选择开关指向右边接收505.2的输出结果。Step 511, judge whether the current symbol Sn is a pilot symbol according to the known time slot format, if the selection switch points to the left to receive the output result of 505.1, if not, select the switch to point to the right to receive the output result of 505.2.

步骤512,如图6所示,移动窗口,使窗口覆盖Cn-L+i.....Cn,对窗口内的信道估计值求算数平均值,得到窗口对应的m时刻的最终信道估计值CEm,即 CEm = 1 L Σ i = 1 L C n - L + i , 其中m>=n-int(L-1)/2,int为向下取整数,L为窗口长度,L根据移动速度不同取值6-60,速度越慢L取值可以越大,L可以根据多普勒频偏估计得到的移动速度来决定。m=n-int(L-1)/2时,补偿时时间对齐,性能会最好,但对于慢速衰落,即L比较大时,将需要较大的硬件存储空间,因此一般会取(n-int(L-1)/2,n)中间的一个折衷的值。CEm是m时刻的最终信道估计值,它将用于对DPDCH(Dedicated Physical Data Channel专用物理数据信道)做最终补偿,补偿结果即可经MRC(Maximal RatioCombining,最大比合并)后送给译码器译码。Step 512, as shown in Figure 6, move the window so that the window covers Cn-L+i.....Cn, calculate the arithmetic mean value of the channel estimation value in the window, and obtain the final channel estimation value at time m corresponding to the window CEm, namely CEm = 1 L Σ i = 1 L C no - L + i , Among them, m>=n-int(L-1)/2, int is an integer rounded down, L is the window length, L takes a value of 6-60 according to the moving speed, the slower the speed, the larger the value of L, and L can be It is determined based on the moving speed estimated by Doppler frequency offset. When m=n-int(L-1)/2, time alignment during compensation will have the best performance, but for slow fading, that is, when L is relatively large, a large hardware storage space will be required, so generally ( A compromise value between n-int(L-1)/2, n). CEm is the final channel estimation value at time m, which will be used to make final compensation for DPDCH (Dedicated Physical Data Channel), and the compensation result can be sent to the decoder after MRC (Maximal RatioCombining) decoding.

Claims (3)

1, a kind of channel estimation methods that is used for broadband CDMA system is characterized in that, may further comprise the steps:
Step 1: the wideband spread-spectrum signal to the input RAKE receiver carries out descrambling and de-spreading, obtains Dedicated Physical Control Channel and Dedicated Physical Data Channel respectively;
Step 2: calculate channel estimation value: in the current of Dedicated Physical Control Channel and next time slot, take out frequency pilot sign according to known time slot format; To the frequency pilot sign that extracts, each frequency pilot sign and known pilot symbols are carried out the decorrelation computing, obtain the channel estimation value of each frequency pilot sign;
Step 3: in all symbol strings for the current time slots of continuous Dedicated Physical Control Channel through the decorrelation computing and next time slot, from first symbol of current time slots S1, till last symbol SL of current time slots, be total to L symbol, symbol ground moving window one by one, L is length of window, for the pilot tone in the current time slots, the pilot frequency sequence of learning according to time slot format carries out participating in the window arithmetic average after the decorrelation and obtains channel estimation value CEm; For the non-pilot tone of current time slots, first pilot average value CE with next time slot Next-slotWith nearest channel estimation results CE M-1Ask the pre-channel estimation value C ' n of average this non-frequency pilot sign of conduct, C ' n is carried out channel compensation and maximum ratio merging, according to the hard decision result after merging this non-pilot tone is carried out decorrelation then, the result after the decorrelation participates in the window arithmetic average and obtains channel estimation value CEm; Wherein subscript m is a channel estimation value index in time, and n is a Dedicated Physical Control Channel symbol index in time;
Step 4: with the corresponding m constantly of the channel estimation value CEm backoff window Dedicated Physical Data Channel symbol constantly that step 3 obtains, the result after each multipath compensation exports to high specific and merges module.
2, the channel estimation methods that is used for broadband CDMA system according to claim 1 is characterized in that, described step 3 specifically realizes by following steps:
In all symbols for the current time slots behind the descrambling and de-spreading and next time slot, from first symbol of current time slots S1, till last symbol SL of current time slots, enter following flow process a symbol of a symbol, previous symbol is imported next symbol after finishing whole flow process again;
Whether, be first non-frequency pilot sign in current time slots, subscript n is a Dedicated Physical Control Channel symbol index in time, if step to step (5), steps to step (2) if not according to known time slot format if (1) differentiating current sign Sn;
(2) differentiate according to known time slot format whether current sign Sn is frequency pilot sign,, step to step (6) if not if step to step (4);
(3) generate pilot frequency sequence according to known time slot format;
(4) symbol and the known symbol that receive are done the decorrelation computing, wherein known symbol is the symbol of step (3) according to corresponding current location in the local sequence of pilot symbols that produces of known time slot format;
(5) obtain all frequency pilot signs of next time slot, use according to known time slot format again and obtain known pilot sequence, the frequency pilot sign of next time slot is done the decorrelation computing one by one, the decorrelation calculated result is averaged obtains CE then Next-slot, at last with CE Next-slotPreservation is got up;
(6) by CE M-1And CE Next-slotSimply average and obtain the pre-channel estimation value C ' n of current non-frequency pilot sign, i.e. C ' n=(CE M-1+ CE Next-slot)/2; CE wherein M-1Be the result who obtains through whole flow process by previous symbol Sn-1, CE Next-slotBe the result who tries to achieve in step (5);
(7) use current non-frequency pilot sign S ' n after C ' n compensation Sn is compensated;
(8) method that the S ' n in all multipaths is merged with high specific merges the symbol S ' n-source after obtaining merging;
(9) S ' n-source is done hard decision, S ^ n - soursce = ( S ' n - souscrce ≥ 0 ) 1:-1,
Figure C2003101039870004C2
It is exactly the estimated value of the former symbol Sn-source of transmitting terminal;
(10) symbol that the non-frequency pilot sign Sn-source of the Dedicated Physical Control Channel that receives and step (9) are estimated Do the decorrelation computing;
(11) differentiate according to known time slot format whether current sign Sn is frequency pilot sign, if selector switch points to the output result of left side receiving step (4), selector switch points to the output result of the right receiving step (10) if not;
(12) ask the arithmetic mean value, obtain the m final channel estimation value CEm constantly of window correspondence, promptly CEm = 1 L Σ i = 1 L C n - L + i , M>=n-int (L-1)/2 wherein, int is downward round numbers, i is an interim cyclic variable of asking average computation to use.
3, the channel estimation methods that is used for broadband CDMA system according to claim 2, it is characterized in that, in the described step (12), L is different and different according to translational speed, span is 6-60, the slow more L value of speed is big more, estimates that according to Doppler frequency deviation the translational speed that obtains decides.
CN 200310103987 2003-11-04 2003-11-04 Channel estimation method for WCDMA communication system Expired - Fee Related CN1269325C (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 200310103987 CN1269325C (en) 2003-11-04 2003-11-04 Channel estimation method for WCDMA communication system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 200310103987 CN1269325C (en) 2003-11-04 2003-11-04 Channel estimation method for WCDMA communication system

Publications (2)

Publication Number Publication Date
CN1543106A CN1543106A (en) 2004-11-03
CN1269325C true CN1269325C (en) 2006-08-09

Family

ID=34333377

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 200310103987 Expired - Fee Related CN1269325C (en) 2003-11-04 2003-11-04 Channel estimation method for WCDMA communication system

Country Status (1)

Country Link
CN (1) CN1269325C (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1848827B (en) * 2005-04-15 2011-04-20 中兴通讯股份有限公司 Channel estimating method and system for code division multiple access communication system
CN101179293B (en) * 2007-12-13 2011-09-14 华为技术有限公司 Channel estimation method, device and receiver in wireless communication system
CN101217305B (en) 2008-01-02 2011-05-04 华为技术有限公司 A processing method and device of data in special physical data channel
CN102611647B (en) * 2011-01-20 2016-04-13 深圳市中兴微电子技术有限公司 Channel estimation methods and device in Wideband Code Division Multiple Access (WCDMA) RAKE receiver
CN103220238A (en) * 2012-01-20 2013-07-24 电信科学技术研究院 Channel estimation method and device
CN103731378B (en) * 2012-10-10 2017-05-10 京信通信系统(中国)有限公司 Method and device for managing frequency deviation
CN107310227B (en) * 2017-06-08 2019-04-05 东华大学 A kind of anti-stab material based on friction self-locking principle and its making method
CN110651454B (en) * 2018-04-27 2022-09-06 深圳市汇顶科技股份有限公司 Channel estimation method, device and receiver
CN115801502B (en) * 2022-11-09 2025-08-29 北京邮电大学 A multi-slot channel estimation method and device based on neural network

Also Published As

Publication number Publication date
CN1543106A (en) 2004-11-03

Similar Documents

Publication Publication Date Title
US6452917B1 (en) Channel estimation in a CDMA wireless communication system
US6414988B1 (en) Amplitude and phase estimation method in a wireless communication system
CN1087527C (en) A method and apparatus for identifying a coded communication signal
CN1377530A (en) Adaptive Channel Estimation in Wireless Communication Systems
CN1346547A (en) Approximated MMSE-based channel estimator in a mobile communication system
CN1849793A (en) Unbiased signal to interference ratio in wireless communications devices and methods therefor
CN1430829A (en) Method and system for channel estimation using multi-slot average interpolation
CN1483250A (en) Parameter estimation for adaptive antenna system
CN1770650A (en) Maximum ratio combining of channel estimation for joint detection in TD-SCDMA systems
CN1269325C (en) Channel estimation method for WCDMA communication system
CN1114294C (en) Speed adaptive channel estimation method and its equipment
CN100486127C (en) Channel evaluation method in wide band CDMA communication system
CN1630997A (en) A method and an apparatus for Eb/Nt estimation for forward power control in spread spectrum communications systems
CN1151622C (en) Pilot channel tracking method based on multipath channel energy window gravity center tracking loop
CN1286534A (en) Channel estimating method used in WCDMA receiver
CN1355626A (en) Improved method for joint test spectrum-spread multi-user communicaltion system
CN1890891A (en) Method and apparatus of noise variance estimation for use in wireless communication systems
CN1312625A (en) A method and device for dynamically adjusting average length for channel estimation
CN1131614C (en) Ruike and equalization cascade receiving method under the code division multiple address low band-spreading ratio and its equipment
CN1260900C (en) Multi-path combining method for WCDMA communication system
CN1688144A (en) Receiver and method for ID SCDMA system
CN1269324C (en) SIR measuring method and device for WCDMA system up-link
CN102611647B (en) Channel estimation methods and device in Wideband Code Division Multiple Access (WCDMA) RAKE receiver
CN1848827B (en) Channel estimating method and system for code division multiple access communication system
CN1131615C (en) Pilot interference cancellation method of code division multiple address system

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
ASS Succession or assignment of patent right

Owner name: HAIMEN TECHNOLOGY DEVELOPMENT CORP.

Free format text: FORMER OWNER: ZTE CORPORATION

Effective date: 20130722

C41 Transfer of patent application or patent right or utility model
COR Change of bibliographic data

Free format text: CORRECT: ADDRESS; FROM: 518057 SHENZHEN, GUANGDONG PROVINCE TO: 226144 NANTONG, JIANGSU PROVINCE

TR01 Transfer of patent right

Effective date of registration: 20130722

Address after: 226144, No. 600, Beijing Road, Haimen, Jiangsu, Nantong province (room 0212 of administrative center)

Patentee after: Haimen science and Technology Development General Corporation

Address before: 518057 Nanshan District, Guangdong high tech Industrial Park, science and Technology Industrial Park, ZTE building, block A, layer 6, layer

Patentee before: ZTE Corporation

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20060809

Termination date: 20161104