CN101848183B - Channel Estimation Method and Device in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing System - Google Patents
Channel Estimation Method and Device in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing System Download PDFInfo
- Publication number
- CN101848183B CN101848183B CN 201010139474 CN201010139474A CN101848183B CN 101848183 B CN101848183 B CN 101848183B CN 201010139474 CN201010139474 CN 201010139474 CN 201010139474 A CN201010139474 A CN 201010139474A CN 101848183 B CN101848183 B CN 101848183B
- Authority
- CN
- China
- Prior art keywords
- antenna
- data
- frequency
- transmitting antenna
- time domain
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 29
- 230000004044 response Effects 0.000 claims abstract description 93
- 238000009792 diffusion process Methods 0.000 claims abstract description 38
- 230000009466 transformation Effects 0.000 claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 24
- 239000000284 extract Substances 0.000 claims abstract description 17
- 238000004364 calculation method Methods 0.000 claims abstract description 9
- 230000001629 suppression Effects 0.000 claims description 22
- 238000006243 chemical reaction Methods 0.000 claims description 19
- 125000004122 cyclic group Chemical group 0.000 claims description 16
- 238000013507 mapping Methods 0.000 claims description 13
- 239000000969 carrier Substances 0.000 claims description 7
- 238000003780 insertion Methods 0.000 claims description 6
- 230000037431 insertion Effects 0.000 claims description 6
- 238000005070 sampling Methods 0.000 claims description 6
- 238000010586 diagram Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000005540 biological transmission Effects 0.000 description 3
- 239000011159 matrix material Substances 0.000 description 3
- 238000005562 fading Methods 0.000 description 2
- 238000004891 communication Methods 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
Images
Landscapes
- Radio Transmission System (AREA)
Abstract
Description
技术领域 technical field
本发明涉及无线通信领域,特别地涉及一种多入多出正交频分复用系统中信道估计方法及装置。The invention relates to the field of wireless communication, in particular to a channel estimation method and device in a multiple-input multiple-output orthogonal frequency division multiplexing system.
背景技术 Background technique
针对存在虚子载波的MIMO-OFDM(多入多出正交频分复用)系统,目前现有的信道估计方法主要是频域插值方法,即线性插值、高斯插值、维纳滤波插值以及变换域插值,但是现有的变化域插值只能有效抑制部分噪声,但却无法抑制由于虚子载波存在所导致的时域径扩散效应,从而使得在较高信噪比时,信道估计出现“地板效应”,时域径扩散所带来的性能恶化现象,使得信道估计的精度偏差大。For MIMO-OFDM (Multi-Input Multiple-Output Orthogonal Frequency Division Multiplexing) systems with virtual subcarriers, the existing channel estimation methods are mainly frequency domain interpolation methods, namely linear interpolation, Gaussian interpolation, Wiener filter interpolation and transform domain However, the existing variable domain interpolation can only effectively suppress part of the noise, but it cannot suppress the time-domain path diffusion effect caused by the existence of virtual subcarriers, so that when the signal-to-noise ratio is high, the channel estimation has a "floor effect" , the performance deterioration caused by the time-domain path diffusion makes the accuracy of channel estimation deviate greatly.
由此可见,当前需要一种新的信道估计的技术方案,来解决当前存在虚子载波的MIMO-OFDM系统中采用的频域插值方法存在的上述问题。It can be seen that a new channel estimation technical solution is currently needed to solve the above-mentioned problems existing in the frequency domain interpolation method used in the current MIMO-OFDM system with virtual subcarriers.
发明内容 Contents of the invention
本发明所要解决的技术问题是提供一种多入多出正交频分复用系统中信道估计方法及装置,解决了抑制进行时域信道估计时径扩散现象的问题,同时也可有效抑制部分噪声,从而最终提高了信道估计的精度。The technical problem to be solved by the present invention is to provide a channel estimation method and device in a multiple-input multiple-output OFDM system, which solves the problem of suppressing the phenomenon of time-path diffusion in time-domain channel estimation, and can effectively suppress some noise, which ultimately improves the accuracy of channel estimation.
为了解决上述问题,本发明提供了一种多入多出正交频分复用系统中信道估计方法,包括:In order to solve the above problems, the present invention provides a channel estimation method in a multiple-input multiple-output OFDM system, comprising:
系统对各发送天线的导频子载波数据与相应发送天线所分配的沃尔什Walsh正交码进行编码后,将编码后的导频子载波块和数据子载波组成频域的正交频分复用OFDM信号;通过对每根发送天线上所述频域的OFDM信号进行反向快速傅立叶变换IFFT变换后,生成时域的OFDM信号传输到接收天线;After the system encodes the pilot subcarrier data of each transmitting antenna and the Walsh orthogonal code assigned by the corresponding transmitting antenna, the coded pilot subcarrier block and data subcarrier form an orthogonal frequency division in the frequency domain Multiplexing the OFDM signal; after performing inverse fast Fourier transform IFFT to the OFDM signal in the frequency domain on each transmitting antenna, generating an OFDM signal in the time domain and transmitting it to the receiving antenna;
系统对每根接收天线上接收的时域OFDM信号进行快速傅立叶变换FFT变换处理,提取每根接收天线上的连续导频子载波块数据,与发送端相匹配的每根发送天线的Walsh正交码进行相乘并累加,获得每根发送天线上的导频接收数据;The system performs fast Fourier transform (FFT) transformation processing on the time-domain OFDM signal received on each receiving antenna, extracts the continuous pilot subcarrier block data on each receiving antenna, and is orthogonal to the Walsh of each transmitting antenna matching the transmitting end. The codes are multiplied and accumulated to obtain the pilot reception data on each transmit antenna;
系统通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,根据变换域方式获得时域信道响应值;对所获得的时域信道响应值通过峰值搜索、峰值点扩散系数计算以及抵消迭代,输出的有效的时域多径信道响应值,获得各发送天线上各数据子载波的信道频率响应值,完成信道估计。The system performs frequency domain LS channel estimation on the pilot received data on each transmit antenna, and obtains the frequency response of all pilot subcarriers on each transmit antenna, and then obtains the time domain channel response value according to the transform domain method; The time-domain channel response value of the time-domain channel response value is obtained through peak search, peak point spread coefficient calculation and offset iteration, and the effective time-domain multipath channel response value is output, and the channel frequency response value of each data subcarrier on each transmitting antenna is obtained to complete the channel estimation.
进一步地,上述方法还可包括,所述系统对各发送天线的导频子载波数据与相应发送天线所分配的Walsh正交码进行编码后,将编码后的导频子载波块和数据子载波组成频域的OFDM信号;通过对每根发送天线上所述频域的OFDM信号进行IFFT变换后,生成时域的OFDM信号传输到接收天线,具体分为以下步骤:Further, the above method may further include, after the system encodes the pilot subcarrier data of each transmitting antenna and the Walsh orthogonal code allocated by the corresponding transmitting antenna, and encodes the coded pilot subcarrier block and data subcarrier block Form the OFDM signal in the frequency domain; after performing IFFT transformation on the OFDM signal in the frequency domain on each transmitting antenna, generate the OFDM signal in the time domain and transmit it to the receiving antenna, which is specifically divided into the following steps:
所述系统将发送天线上的连续导频子载波块数据与相应各发送天线所分配的Walsh正交码进行编码,对连续导频块内的非导频子载波和非导频数据块内数据子载波上的数据由单独各发送天线上的比特流通过星座映射完成,将编码后的导频子载波块和星座映射完成后的数据子载波组成频域的OFDM信号;The system encodes the continuous pilot subcarrier block data on the transmitting antennas with the Walsh orthogonal codes assigned to the corresponding transmitting antennas, and encodes the non-pilot subcarriers in the continuous pilot blocks and the data in the non-pilot data blocks The data on the subcarriers is completed by constellation mapping of the bit streams on the individual transmit antennas, and the coded pilot subcarrier blocks and the data subcarriers after constellation mapping are completed to form OFDM signals in the frequency domain;
所述系统对每根发送天线上所述频域的OFDM信号进行IFFT变换后,进行并串转换同时添加循环前缀后,生成时域的OFDM信号,通过无线多入多出MIMO信道传输到接收天线。The system performs IFFT transformation on the OFDM signal in the frequency domain on each transmitting antenna, performs parallel-to-serial conversion and adds a cyclic prefix at the same time, generates an OFDM signal in the time domain, and transmits it to the receiving antenna through a wireless multiple-input multiple-output MIMO channel .
进一步地,上述方法还可包括,所述系统对每根接收天线上接收的时域OFDM信号进行FFT变换处理,提取每根接收天线上的连续导频子载波块数据,具体是指:Further, the above method may also include that the system performs FFT transformation processing on the time-domain OFDM signal received on each receiving antenna, and extracts continuous pilot subcarrier block data on each receiving antenna, specifically referring to:
所述系统对每根接收天线上接收的时域OFDM信号进行去除循环前缀和串并转换后,通过FFT变换处理,生成每根接收天线上频域的OFDM信号,提取每根接收天线上的连续导频子载波块数据。After the system removes the cyclic prefix and serial-to-parallel conversion of the time-domain OFDM signal received on each receiving antenna, it generates an OFDM signal in the frequency domain on each receiving antenna through FFT transformation processing, and extracts the continuous signal on each receiving antenna. Pilot subcarrier block data.
进一步地,上述方法还可包括,所述系统通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,根据变换域方式获得时域信道响应值,具体是指:Further, the above method may also include that the system performs frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, and after obtaining the frequency responses of all pilot subcarriers on each transmitting antenna, according to the transform domain Way to obtain the time domain channel response value, specifically refers to:
所述系统通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,插入若干个零,进而再进行IFFT变换,获得时域信道响应值,其中,插入零的个数是由导频子载波个数、FFT点数以及导频子载波的插入间隔所确定。The system performs frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, obtains the frequency response of all pilot subcarriers on each transmitting antenna, inserts several zeros, and then performs IFFT transformation to obtain The time-domain channel response value, wherein the number of inserted zeros is determined by the number of pilot subcarriers, the number of FFT points, and the insertion interval of pilot subcarriers.
进一步地,上述方法还可包括,所述系统对所获得的时域信道响应值通过峰值搜索、峰值点扩散系数计算以及抵消迭代,输出的有效的时域多径信道响应值,获得各发送天线上各数据子载波的信道频率响应值,完成信道估计,具体分为以下步骤:Further, the above method may further include that the system outputs an effective time-domain multipath channel response value through peak search, peak point spread coefficient calculation, and offset iteration for the obtained time-domain channel response value, and obtains each transmitting antenna The channel frequency response value of each data subcarrier is used to complete the channel estimation, which is specifically divided into the following steps:
所述系统对所获得的时域信道响应值,通过进行峰值搜索确定峰值及其对应位置后,计算该位置采样点在所有时域样点上的扩散系数,并将其从最初所获得的时域信道响应值中抵消,再进行新一轮的峰值搜索迭代,直至搜索出的该峰值降到设定门限以下,停止迭代,对输出的有效的时域多径信道响应值进行FFT变换处理,获得各发送天线上各数据子载波的信道频率响应值,完成信道估计。For the obtained time-domain channel response value, the system determines the peak value and its corresponding position by performing a peak search, then calculates the diffusion coefficient of the sample point at this position on all time-domain sample points, and converts it from the initially obtained time-domain channel response value to Offset in the domain channel response value, and then perform a new round of peak search iterations until the searched peak value drops below the set threshold, stop the iteration, and perform FFT transformation processing on the output effective time domain multipath channel response value, The channel frequency response value of each data subcarrier on each transmitting antenna is obtained to complete channel estimation.
本发明还提供了一种多入多出正交频分复用系统中信道估计装置,包括天线编码模块,天线解码模块和径扩散抑制模块,其中,The present invention also provides a channel estimation device in a multiple-input multiple-output OFDM system, including an antenna encoding module, an antenna decoding module and a radial diffusion suppression module, wherein,
所述天线编码模块,用于对各发送天线的导频子载波数据与相应发送天线所分配的沃尔什Walsh正交码进行编码后,将编码后的导频子载波块和数据子载波组成频域的正交频分复用OFDM信号;通过对每根发送天线上所述频域的OFDM信号进行反向快速傅立叶变换IFFT变换后,生成时域的OFDM信号传输到所述天线解码模块;The antenna encoding module is configured to encode the pilot subcarrier data of each transmitting antenna and the Walsh orthogonal code assigned by the corresponding transmitting antenna, and then form the encoded pilot subcarrier blocks and data subcarriers Orthogonal frequency division multiplexing OFDM signal in the frequency domain; after performing inverse fast Fourier transform IFFT to the OFDM signal in the frequency domain on each transmitting antenna, generate an OFDM signal in the time domain and transmit it to the antenna decoding module;
所述天线解码模块,用于对接收的每根接收天线上时域OFDM信号进行快速傅立叶变换FFT变换处理,提取每根接收天线上的连续导频子载波块数据,与发送端相匹配的每根发送天线的Walsh正交码进行相乘并累加,获得每根发送天线上的导频接收数据,并发送到所述径扩散抑制模块;The antenna decoding module is used to perform fast Fourier transform (FFT) transform processing on the received time-domain OFDM signal on each receiving antenna, extract continuous pilot subcarrier block data on each receiving antenna, and match each The Walsh orthogonal codes of the root transmit antennas are multiplied and accumulated to obtain the pilot reception data on each transmit antenna, and sent to the radial diffusion suppression module;
所述径扩散抑制模块,用于通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,根据变换域方式获得时域信道响应值;对所获得的时域信道响应值通过峰值搜索、峰值点扩散系数计算以及抵消迭代,输出的有效的时域多径信道响应值,获得各发送天线上各数据子载波的信道频率响应值,完成信道估计。The diameter diffusion suppression module is used to perform frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, and obtain the frequency responses of all pilot subcarriers on each transmitting antenna, and then obtain the time Domain channel response value; through peak search, peak point spread coefficient calculation and offset iteration for the obtained time domain channel response value, the effective time domain multipath channel response value output is obtained to obtain the channel of each data subcarrier on each transmitting antenna The frequency response value completes the channel estimation.
进一步地,上述装置还可包括,所述天线编码模块对各发送天线的导频子载波数据与相应发送天线所分配的Walsh正交码进行编码后,将编码后的导频子载波块和数据子载波组成频域的OFDM信号;通过对每根发送天线上所述频域的OFDM信号进行IFFT变换后,生成时域的OFDM信号,具体是指:Further, the above device may further include, after the antenna encoding module encodes the pilot subcarrier data of each transmitting antenna and the Walsh orthogonal code allocated by the corresponding transmitting antenna, and encodes the encoded pilot subcarrier block and data The subcarriers form the OFDM signal in the frequency domain; after performing IFFT transformation on the OFDM signal in the frequency domain on each transmitting antenna, the OFDM signal in the time domain is generated, specifically referring to:
所述天线编码模块将发送天线上的连续导频子载波块数据与相应各发送天线所分配的Walsh正交码进行编码,对连续导频块内的非导频子载波和非导频数据块内数据子载波上的数据由单独各发送天线上的比特流通过星座映射完成,将编码后的导频子载波块和星座映射完成后的数据子载波组成频域的OFDM信号;对每根发送天线上所述频域的OFDM信号进行IFFT变换后,进行并串转换同时添加循环前缀后,生成时域的OFDM信号。The antenna encoding module encodes the continuous pilot subcarrier block data on the transmitting antenna and the Walsh orthogonal codes allocated by the corresponding transmitting antennas, and encodes the non-pilot subcarriers and non-pilot data blocks in the continuous pilot block The data on the inner data subcarrier is completed by the bit stream on each transmit antenna through constellation mapping, and the coded pilot subcarrier block and the data subcarrier after the constellation mapping are completed to form an OFDM signal in the frequency domain; The OFDM signal in the frequency domain on the antenna is subjected to IFFT transformation, parallel-to-serial conversion is performed and a cyclic prefix is added at the same time, and an OFDM signal in the time domain is generated.
进一步地,上述装置还可包括,所述天线解码模块对每根接收天线上接收的时域OFDM信号进行FFT变换处理,提取每根接收天线上的连续导频子载波块数据,具体是指:Further, the above-mentioned device may further include that the antenna decoding module performs FFT transformation processing on the time-domain OFDM signal received on each receiving antenna, and extracts continuous pilot subcarrier block data on each receiving antenna, specifically referring to:
所述天线解码模块对每根接收天线上接收的时域OFDM信号进行去除循环前缀和串并转换后,通过FFT变换处理,生成每根接收天线上频域的OFDM信号,提取每根接收天线上的连续导频子载波块数据。After the antenna decoding module removes the cyclic prefix and serial-to-parallel conversion of the time-domain OFDM signal received on each receiving antenna, it generates an OFDM signal in the frequency domain on each receiving antenna through FFT transformation processing, and extracts the OFDM signal on each receiving antenna. Consecutive pilot subcarrier block data of .
进一步地,上述装置还可包括,所述径扩散抑制模块通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,根据变换域方式获得时域信道响应值,具体是指:Further, the above-mentioned device may further include that the channel diffusion suppression module performs frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, and after obtaining the frequency responses of all pilot subcarriers on each transmitting antenna, Obtaining the channel response value in the time domain according to the transformation domain method, specifically refers to:
所述径扩散抑制模块通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,插入若干个零,进而再进行IFFT变换,获得时域信道响应值,其中,插入零的个数是由导频子载波个数、FFT点数以及导频子载波的插入间隔所确定。The radial diffusion suppression module performs frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, and after obtaining the frequency responses of all pilot subcarriers on each transmitting antenna, inserts several zeros, and then performs IFFT Transform to obtain the time-domain channel response value, wherein the number of inserted zeros is determined by the number of pilot subcarriers, the number of FFT points, and the insertion interval of pilot subcarriers.
进一步地,上述装置还可包括,所述径扩散抑制模块对所获得的时域信道响应值通过峰值搜索、峰值点扩散系数计算以及抵消迭代,输出的有效的时域多径信道响应值,获得各发送天线上各数据子载波的信道频率响应值,具体是指:Further, the above-mentioned device may further include that the effective time-domain multi-path channel response value outputted by the path diffusion suppression module through peak search, calculation of peak point spread coefficient and cancellation iteration for the obtained time-domain channel response value is obtained. The channel frequency response value of each data subcarrier on each transmit antenna specifically refers to:
所述径扩散抑制模块对所获得的时域信道响应值,通过进行峰值搜索确定峰值及其对应位置后,计算该位置采样点在所有时域样点上的扩散系数,并将其从最初所获得的时域信道响应值中抵消,再进行新一轮的峰值搜索迭代,直至搜索出的该峰值降到设定门限以下,停止迭代,对输出的有效的时域多径信道响应值进行FFT变换处理,获得各发送天线上各数据子载波的信道频率响应值。For the obtained time-domain channel response value, the said diameter diffusion suppression module determines the peak value and its corresponding position by performing a peak search, calculates the diffusion coefficient of the sampling point at this position on all time-domain sample points, and converts it from the initially obtained Offset the obtained time-domain channel response value, and then perform a new round of peak search iterations until the searched peak value drops below the set threshold, stop the iteration, and perform FFT on the output effective time-domain multipath channel response value Transform processing to obtain the channel frequency response value of each data subcarrier on each transmit antenna.
与现有技术相比,应用本发明,对存在虚子载波的MIMO-OFDM系统在信道估计时出现的径扩散效应进行了一定的抑制,在很大程度上削弱了“地板效应”,同时也抑制了部分噪声,有效的改善了信道估计的精度,且其复杂度较低,实现相对简单。Compared with the prior art, the application of the present invention suppresses the channel diffusion effect of the MIMO-OFDM system with virtual subcarriers during channel estimation, weakens the "floor effect" to a large extent, and also suppresses Part of the noise is eliminated, the accuracy of channel estimation is effectively improved, and its complexity is low, and its implementation is relatively simple.
附图说明 Description of drawings
图1是本发明的多入多出正交频分复用系统中信道估计方法的流程图;Fig. 1 is the flow chart of channel estimation method in the multiple-input multiple-output OFDM system of the present invention;
图2是实例中入多出正交频分复用系统中信道估计的示意图;Fig. 2 is the schematic diagram of channel estimation in the input multiple output OFDM system in the example;
图3是实例中天线编码模块中结构示意图;Fig. 3 is a structural schematic diagram in the antenna coding module in the example;
图4是实例中天线解码模块中结构示意图;Fig. 4 is a structural schematic diagram in the antenna decoding module in the example;
图5是实例中径扩散抑制模块中结构示意图。Fig. 5 is a schematic diagram of the structure of the median diameter diffusion suppression module of the example.
具体实施方式 Detailed ways
下面结合附图和具体实施方式对本发明作进一步说明。The present invention will be further described below in conjunction with the accompanying drawings and specific embodiments.
本发明的实施条件是针对存在虚子载波的MIMO-OFDM系统。(其中,MIMO信道条件:块衰落(即在多个OFDM符号内的相同子载波上频率信道响应不变)+白噪声)。本发明的实施可以通过存储器和多个计算器(如乘法器、加法器等)实现。The implementation condition of the present invention is aimed at the MIMO-OFDM system with virtual sub-carriers. (Wherein, MIMO channel condition: block fading (that is, the frequency channel response is constant on the same subcarrier within multiple OFDM symbols) + white noise). The implementation of the present invention can be realized through memory and multiple calculators (such as multipliers, adders, etc.).
下面以NT×NR的MIMO-OFDM系统为例,(其中,NT为发送天线,NR为接收天线)本发明的多入多出正交频分复用系统中信道估计方法,包括以下步骤:The MIMO-OFDM system of NT × NR is taken as an example below, (wherein, NT is a transmitting antenna, and NR is a receiving antenna) the channel estimation method in the multiple-input multiple-output OFDM system of the present invention includes The following steps:
步骤110、系统将发送天线上的连续导频子载波块数据与相应各发送天线所分配的Walsh(沃尔什)正交码进行编码,对连续导频块内的非导频子载波和非导频数据块内数据子载波上的数据由单独各发送天线上的比特流通过星座映射完成,将编码后的导频子载波块和星座映射完成后的数据子载波组成频域的OFDM信号;
系统将NT根发送天线上的连续导频子载波块(其中,导频块长度为NC,各导频子载波上的数据全为“1”)与相应各发送天线所分配的Walsh正交码(码长也为NC)进行相乘(即进行编码),对于连续导频块内的非导频子载波和非导频数据块内数据子载波上的数据由单独各发送天线上的比特流通过星座映射完成。The system connects the continuous pilot subcarrier blocks on NT transmitting antennas (the length of the pilot block is N C , and the data on each pilot subcarrier is all "1") with the corresponding Walsh normal Cross-code (the code length is also N C ) is multiplied (i.e., encoded), and the data on the non-pilot subcarriers in the continuous pilot block and the data subcarriers in the non-pilot data block are transmitted by each individual transmit antenna The bitstream is done through constellation mapping.
步骤120、系统对每根发送天线上所述频域的OFDM信号进行NFFT点的IFFT变换后,进行并串转换(P/S)同时添加循环前缀(CP)后生成时域的OFDM信号,通过无线MIMO信道传输到接收天线;
系统对每根发送天线上进行NFFT点的IFFT(反向快速傅立叶变换)变换,再经并串转换(P/S)并添加循环前缀(CP)后成为时域的OFDM信号,最后通过无线衰落信道传输到接收端。The system performs IFFT (inverse fast Fourier transform) transformation of NFFT points on each transmitting antenna, and then undergoes parallel-to-serial conversion (P/S) and adds cyclic prefix (CP) to become an OFDM signal in the time domain. The fading channel is transmitted to the receiver.
步骤130、系统对每根接收天线上接收的时域OFDM信号进行去除循环前缀和串并转换(S/P)后,通过NFFT点的FFT(快速傅立叶变换)变换处理,生成每根接收天线上频域的OFDM信号,提取每根接收天线上的连续导频子载波块数据;
系统对每根接收天线上的接收数据首先去除循环前缀(CP)、串并转换(S/P),再通过NFFT点的FFT变换处理,最后提取每根接收天线上的连续导频子载波块数据。The system first removes the cyclic prefix (CP) and serial-to-parallel conversion (S/P) of the received data on each receiving antenna, and then processes through the FFT transformation of NFFT points, and finally extracts the continuous pilot subcarriers on each receiving antenna block data.
步骤140、系统将与发送端相匹配的每根发送天线的Walsh正交码,和提取后的每根接收天线上的导频子载波块数据进行相乘并累加,获得每根发送天线上的导频接收数据;
系统将与发送端相匹配的每根发送天线的Walsh正交码,将其与步骤130中最后所提取的每根接收天线上的导频子载波块数据进行相乘并累加,从而获得每根发送天线上的导频接收数据,即同时去除了其它发送天线带来的干扰。The system multiplies and accumulates the Walsh orthogonal code of each transmitting antenna matched with the transmitting end with the pilot subcarrier block data on each receiving antenna extracted in
步骤150、系统通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,根据变换域方式(即插入若干个零),进而再进行MFFT点的IFFT变换,获得MFFT点时域信道响应值;
系统通过频域LS信道估计方法,可首先获得每根发送天线上所有导频子载波的频率响应,并插入若干个零(个数是由导频子载波个数Np、FFT点数NFFT以及导频子载波的插入间隔ΔNf所决定),进而再进行MFFT点的IFFT变换。Through the frequency domain LS channel estimation method, the system can first obtain the frequency response of all pilot subcarriers on each transmitting antenna, and insert several zeros (the number is determined by the number of pilot subcarriers N p , the number of FFT points N FFT and The insertion interval ΔN f of the pilot subcarrier is determined), and then the IFFT transformation of the MFFT point is performed.
步骤160、系统对所获得的MFFT点时域信道响应值,首先进行峰值搜索,确定峰值及其对应位置,进而计算该位置采样点在所有时域样点上的扩散系数,并将其从最初所获得的MFFT点时域信道响应值中抵消,再进行新一轮的峰值搜索迭代,直至搜索出的该峰值降到设定门限以下,停止迭代,输出有效的MFFT点时域多径信道响应值,将其余样点全部置为0;
步骤170、系统对获得的MFFT点时域多径信道响应值,进行NFFT点FFT变换处理,获得各发送天线上各数据子载波的信道频率响应值,完成信道估计。
下面结合具体实例对本发明作进一步说明。The present invention will be further described below in conjunction with specific examples.
图2中符号解释如下:The symbols in Figure 2 are explained as follows:
第nt(nt=1,2,…NT)根发送天线上由导频子载波与数据子载波组成的第n个OFDM符号频域数据块向量; The nth OFDM symbol frequency-domain data block vector composed of pilot subcarriers and data subcarriers on the n t (n t =1, 2, ... N T ) root transmitting antenna;
第nt(nt=1,2,…NT)根发送天线上由导频子载波与数据子载波组成 The n t (n t = 1, 2, ... N T ) transmitting antenna is composed of pilot subcarriers and data subcarriers
的第n个OFDM符号时域数据块向量;The nth OFDM symbol time-domain data block vector of ;
h(n):无线MIMO信道的时域冲激响应,通常用其离散抽头延迟线模型表示为 h(n): The time-domain impulse response of a wireless MIMO channel, usually expressed by its discrete tapped delay line model as
第nr(nr=1,2,…NR)根接收天线上的时域接收信号; The time-domain received signal on the n rth (n r =1, 2, ... N R ) receiving antenna;
第nr(nr=1,2,…NR)根接收天线上去循环前缀(CP)后的频域接收信号; The frequency-domain received signal after removing the cyclic prefix (CP) on the n r (n r =1, 2, ... N R ) receiving antenna;
第nt根发送天线到第nr根接收天线之间时域信道冲激响应; The time-domain channel impulse response between the n t transmitting antenna and the n r receiving antenna;
第nt根发送天线到第nr根接收天线之间信道的频域响应; The frequency domain response of the channel between the n t transmitting antenna and the n r receiving antenna;
A:天线编码模块,完成各发送天线的上导频子载波数据的正交编码;A: The antenna encoding module completes the orthogonal encoding of the upper pilot subcarrier data of each transmitting antenna;
B:天线解码模块,即完成模块A的逆处理,主要由乘法器与累加器组成;B: Antenna decoding module, which completes the inverse processing of module A, mainly composed of a multiplier and an accumulator;
C:径扩散抑制模块,完成时域径扩散、噪声的抑制。C: Radial diffusion suppression module, which completes the time-domain radial diffusion and noise suppression.
图2中描述的信号流的过程为:The process of the signal flow described in Figure 2 is:
(1)模块A输出各发送天线上的导频子载波数据;(1) Module A outputs pilot subcarrier data on each transmit antenna;
(2)各发送天线上时域发送信号的产生:(2) Generation of time-domain transmission signals on each transmission antenna:
其中,TCP为循环前缀添加矩阵,为归一化的傅里叶变换矩阵。Among them, T CP adds the matrix for the cyclic prefix, is the normalized Fourier transform matrix.
(3)上述发送信号经过无线MIMO信道后,在接收端的接收形式为:(3) After the above transmitted signal passes through the wireless MIMO channel, the receiving form at the receiving end is:
其中,是由产生的循环矩阵,为高斯噪声。in, By The resulting circulant matrix, is Gaussian noise.
(4)进一步对上述每根接收天线上的接收数据进行FFT变化处理可得到其频域接收向量为:(4) Further perform FFT change processing on the received data on each of the above receiving antennas to obtain its frequency domain receiving vector as:
(5)模块B对所抽取的各天线上导频子载波数据进行正交分离;(5) Module B performs orthogonal separation on the extracted pilot subcarrier data on each antenna;
(6)模块C输出为有效的时域多径信道响应值,并通过NFFT点FFT变换得到最终各发送天线上各数据子载波的信道频率响应:(6) The output of module C is an effective time-domain multipath channel response value, and the channel frequency response of each data subcarrier on each transmitting antenna is finally obtained by NFFT point FFT transformation:
模块A中内部结构如图3所示。The internal structure of module A is shown in Figure 3.
图3中符号解释如下:The symbols in Figure 3 are explained as follows:
第nt(nt=1,2,…NT)根发送天线所分配的Walsh正交码,且有其中,各发送天线所分配不同的Walsh正交码具有如下正交特性,即 The Walsh orthogonal code assigned by the n t (n t =1, 2, ... N T )th transmit antenna, and has Among them, the different Walsh orthogonal codes assigned to each transmitting antenna have the following orthogonal characteristics, namely
第nt(nt=1,2,…NT)根发送天线上经编码后第kp个导频子载波的时域数据块,且有 The coded time-domain data block of the k p -th pilot subcarrier on the n t (n t =1, 2, ... N T )th transmitting antenna, and has
图3中描述的信号流过程为:即 The signal flow process described in Figure 3 is: Right now
模块B中内部结构如图4所示。The internal structure of module B is shown in Figure 4.
图4中符号解释如下:The symbols in Figure 4 are explained as follows:
第nr(nr=1,2,…NR)根接收信号中所提取的第kp个导频子载波处的连续导频块向量,且有或者 The continuous pilot block vector at the k pth pilot subcarrier extracted from the n r (n r =1, 2, ... NR ) root received signal, and has or
第nt(nt=1,2,…NT)根发送天线到nr(nr=1,2,…NR)根接收天线之间的第kp(p=0,1,…NP-1)个导频子载波处的接收数据。 The k p ( p=0 , 1 , . . . N P -1) received data at pilot subcarriers.
图4中描述的信号流过程为: The signal flow process described in Figure 4 is:
模块C中内部结构如图5所示。The internal structure of module C is shown in Figure 5.
图5中符号解释如下:The symbols in Figure 5 are explained as follows:
第nt(nt=1,2,…NT)根发送天线到nr(nr=1,2,…NR)根接收天线之间初始的时域信道冲激响应估计; The initial time-domain channel impulse response estimation between the n t (n t =1, 2, ... N T ) transmitting antenna and the n r (n r = 1, 2, ... NR ) receiving antenna;
nmax (m):第m次迭代后的峰值点位置, n max (m) : the position of the peak point after the mth iteration,
第m次迭代后的峰值点位置在所有样点上的扩散系数; The diffusion coefficient of the peak point position on all samples after the mth iteration;
经M次迭代抵消后,第nt(nt=1,2,…NT)根发送天线到nr(nr=1,2,…NR)根接收天线之间时域信道冲激响应估计。 After M times of iterative cancellation, the time-domain channel impulse between the n t (n t =1, 2, ... NT )th transmitting antenna and the n r (n r =1, 2, ... NR ) receiving antenna Response estimates.
图5中描述的信号流过程为:The signal flow process described in Figure 5 is:
(1)更新输入信号, (1) update the input signal,
(2)进行峰值搜索, (2) Perform peak search,
(3)计算扩散系数,(3) Calculate the diffusion coefficient,
本发明的一种多入多出正交频分复用系统中信道估计装置,包括天线编码模块,天线解码模块和径扩散抑制模块,其中,A channel estimation device in a multiple-input multiple-output OFDM system of the present invention includes an antenna encoding module, an antenna decoding module and a radial diffusion suppression module, wherein,
所述天线编码模块,用于对各发送天线的导频子载波数据与相应发送天线所分配的沃尔什Walsh正交码进行编码后,将编码后的导频子载波块和数据子载波组成频域的正交频分复用OFDM信号;通过对每根发送天线上所述频域的OFDM信号进行反向快速傅立叶变换IFFT变换后,生成时域的OFDM信号传输到所述天线解码模块;The antenna encoding module is configured to encode the pilot subcarrier data of each transmitting antenna and the Walsh orthogonal code assigned by the corresponding transmitting antenna, and then form the encoded pilot subcarrier blocks and data subcarriers Orthogonal frequency division multiplexing OFDM signal in the frequency domain; after performing inverse fast Fourier transform IFFT to the OFDM signal in the frequency domain on each transmitting antenna, generate an OFDM signal in the time domain and transmit it to the antenna decoding module;
所述天线解码模块,用于对接收的每根接收天线上时域OFDM信号进行快速傅立叶变换FFT变换处理,提取每根接收天线上的连续导频子载波块数据,与发送端相匹配的每根发送天线的Walsh正交码进行相乘并累加,获得每根发送天线上的导频接收数据,并发送到所述径扩散抑制模块;The antenna decoding module is used to perform fast Fourier transform (FFT) transform processing on the received time-domain OFDM signal on each receiving antenna, extract continuous pilot subcarrier block data on each receiving antenna, and match each The Walsh orthogonal codes of the root transmit antennas are multiplied and accumulated to obtain the pilot reception data on each transmit antenna, and sent to the radial diffusion suppression module;
所述径扩散抑制模块,用于通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,根据变换域方式获得时域信道响应值;对所获得的时域信道响应值通过峰值搜索、峰值点扩散系数计算以及抵消迭代,输出的有效的时域多径信道响应值,获得各发送天线上各数据子载波的信道频率响应值,完成信道估计。The diameter diffusion suppression module is used to perform frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, and obtain the frequency responses of all pilot subcarriers on each transmitting antenna, and then obtain the time Domain channel response value; through peak search, peak point spread coefficient calculation and offset iteration for the obtained time domain channel response value, the effective time domain multipath channel response value output is obtained to obtain the channel of each data subcarrier on each transmitting antenna The frequency response value completes the channel estimation.
所述天线编码模块对各发送天线的导频子载波数据与相应发送天线所分配的Walsh正交码进行编码后,将编码后的导频子载波块和数据子载波组成频域的OFDM信号;通过对每根发送天线上所述频域的OFDM信号进行IFFT变换后,生成时域的OFDM信号,具体是指:After the antenna encoding module encodes the pilot subcarrier data of each transmitting antenna and the Walsh orthogonal code assigned by the corresponding transmitting antenna, the encoded pilot subcarrier block and data subcarrier form an OFDM signal in the frequency domain; After performing IFFT transformation on the OFDM signal in the frequency domain on each transmitting antenna, an OFDM signal in the time domain is generated, specifically referring to:
所述天线编码模块将发送天线上的连续导频子载波块数据与相应各发送天线所分配的Walsh正交码进行编码,对连续导频块内的非导频子载波和非导频数据块内数据子载波上的数据由单独各发送天线上的比特流通过星座映射完成,将编码后的导频子载波块和星座映射完成后的数据子载波组成频域的OFDM信号;对每根发送天线上所述频域的OFDM信号进行IFFT变换后,进行并串转换同时添加循环前缀后,生成时域的OFDM信号。The antenna encoding module encodes the continuous pilot subcarrier block data on the transmitting antenna and the Walsh orthogonal codes allocated by the corresponding transmitting antennas, and encodes the non-pilot subcarriers and non-pilot data blocks in the continuous pilot block The data on the inner data subcarrier is completed by the bit stream on each transmit antenna through constellation mapping, and the coded pilot subcarrier block and the data subcarrier after the constellation mapping are completed to form an OFDM signal in the frequency domain; The OFDM signal in the frequency domain on the antenna is subjected to IFFT transformation, parallel-to-serial conversion is performed and a cyclic prefix is added at the same time, and an OFDM signal in the time domain is generated.
所述天线解码模块对每根接收天线上接收的时域OFDM信号进行FFT变换处理,提取每根接收天线上的连续导频子载波块数据,具体是指:The antenna decoding module performs FFT transformation processing on the time-domain OFDM signal received on each receiving antenna, and extracts the continuous pilot subcarrier block data on each receiving antenna, specifically referring to:
所述天线解码模块对每根接收天线上接收的时域OFDM信号进行去除循环前缀和串并转换后,通过FFT变换处理,生成每根接收天线上频域的OFDM信号,提取每根接收天线上的连续导频子载波块数据。After the antenna decoding module removes the cyclic prefix and serial-to-parallel conversion of the time-domain OFDM signal received on each receiving antenna, it generates an OFDM signal in the frequency domain on each receiving antenna through FFT transformation processing, and extracts the OFDM signal on each receiving antenna. Consecutive pilot subcarrier block data of .
所述径扩散抑制模块通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,根据变换域方式获得时域信道响应值,具体是指:The radial spread suppression module performs frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, and obtains the time-domain channel response according to the transform domain method after obtaining the frequency responses of all pilot subcarriers on each transmitting antenna value, specifically:
所述径扩散抑制模块通过对每根发送天线上的导频接收数据进行频域LS信道估计,获得每根发送天线上所有导频子载波的频率响应后,插入若干个零,进而再进行IFFT变换,获得时域信道响应值,其中,插入零的个数是由导频子载波个数、FFT点数以及导频子载波的插入间隔所确定。The radial diffusion suppression module performs frequency-domain LS channel estimation on the pilot received data on each transmitting antenna, and after obtaining the frequency responses of all pilot subcarriers on each transmitting antenna, inserts several zeros, and then performs IFFT Transform to obtain the time-domain channel response value, wherein the number of inserted zeros is determined by the number of pilot subcarriers, the number of FFT points, and the insertion interval of pilot subcarriers.
所述径扩散抑制模块对所获得的时域信道响应值通过峰值搜索、峰值点扩散系数计算以及抵消迭代,输出的有效的时域多径信道响应值,获得各发送天线上各数据子载波的信道频率响应值,具体是指:The path diffusion suppression module performs peak search, peak point spread coefficient calculation and offset iteration on the obtained time-domain channel response value, and the effective time-domain multi-path channel response value outputted obtains each data subcarrier on each transmitting antenna. Channel frequency response value, specifically refers to:
所述径扩散抑制模块对所获得的时域信道响应值,通过进行峰值搜索确定峰值及其对应位置后,计算该位置采样点在所有时域样点上的扩散系数,并将其从最初所获得的时域信道响应值中抵消,再进行新一轮的峰值搜索迭代,直至搜索出的该峰值降到设定门限以下,停止迭代,对输出的有效的时域多径信道响应值进行FFT变换处理,获得各发送天线上各数据子载波的信道频率响应值。For the obtained time-domain channel response value, the said diameter diffusion suppression module determines the peak value and its corresponding position by performing a peak search, calculates the diffusion coefficient of the sampling point at this position on all time-domain sample points, and converts it from the initially obtained Offset the obtained time-domain channel response value, and then perform a new round of peak search iterations until the searched peak value drops below the set threshold, stop the iteration, and perform FFT on the output effective time-domain multipath channel response value Transform processing to obtain the channel frequency response value of each data subcarrier on each transmit antenna.
以上所述,仅为本发明较佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉该技术的人在本发明所揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应该以权利要求的保护范围为准。The above is only a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any person familiar with the technology can easily think of changes or replacements within the technical scope disclosed in the present invention. , should be covered within the protection scope of the present invention. Therefore, the protection scope of the present invention should be determined by the protection scope of the claims.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010139474 CN101848183B (en) | 2010-04-06 | 2010-04-06 | Channel Estimation Method and Device in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing System |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010139474 CN101848183B (en) | 2010-04-06 | 2010-04-06 | Channel Estimation Method and Device in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing System |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101848183A CN101848183A (en) | 2010-09-29 |
CN101848183B true CN101848183B (en) | 2013-02-06 |
Family
ID=42772646
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010139474 Expired - Fee Related CN101848183B (en) | 2010-04-06 | 2010-04-06 | Channel Estimation Method and Device in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing System |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101848183B (en) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103428117B (en) * | 2012-05-14 | 2016-09-07 | 中兴通讯股份有限公司 | The determination method and apparatus of channel estimation methods |
CN103491049B (en) * | 2013-09-30 | 2016-08-31 | 电子科技大学 | Channel estimation methods based on decision-feedback in ofdm system |
CN104243370B (en) * | 2014-09-01 | 2017-06-30 | 中国科学院信息工程研究所 | A kind of time-domain channel estimating method applied in multiaerial system |
CN110166401B (en) * | 2019-07-12 | 2021-07-02 | 电子科技大学 | Phase Noise Suppression Method for Massive MIMO Orthogonal Frequency Division Multiplexing System |
CN112564874B (en) * | 2020-11-04 | 2022-09-13 | 展讯通信(上海)有限公司 | Pilot frequency sequence generation method and device |
CN115883294B (en) * | 2022-11-28 | 2024-11-05 | 中国电子科技集团公司第五十四研究所 | Channel estimation method and device for low-orbit satellite OFDM communication system |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1658528A (en) * | 2004-02-20 | 2005-08-24 | 电子科技大学 | An Adaptive Channel Estimation Method for MIMO-OFDM System |
CN101083646A (en) * | 2006-06-01 | 2007-12-05 | 电子科技大学 | Channel estimation optimizing method for amplitude-limiting OFDM system |
CN101136884A (en) * | 2007-09-28 | 2008-03-05 | 清华大学 | Channel Estimation Method for TDS-OFDM System |
-
2010
- 2010-04-06 CN CN 201010139474 patent/CN101848183B/en not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1658528A (en) * | 2004-02-20 | 2005-08-24 | 电子科技大学 | An Adaptive Channel Estimation Method for MIMO-OFDM System |
CN101083646A (en) * | 2006-06-01 | 2007-12-05 | 电子科技大学 | Channel estimation optimizing method for amplitude-limiting OFDM system |
CN101136884A (en) * | 2007-09-28 | 2008-03-05 | 清华大学 | Channel Estimation Method for TDS-OFDM System |
Also Published As
Publication number | Publication date |
---|---|
CN101848183A (en) | 2010-09-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107332797B (en) | Channel estimation method in power line OFDM communication system | |
CN101494528B (en) | Training sequence design and channel estimation method of transmission diversity block transmission system | |
CN103391268B (en) | Multi-I/O OFDM shallow-sea underwater acoustic communication pattern selection peak-to-average ratio suppressing method | |
CN100385824C (en) | An Adaptive Channel Estimation Method for MIMO-OFDM System | |
CN102291351B (en) | Timing synchronization method of receiver in OFDM wireless communication system | |
CN106302298B (en) | A Method for Eliminating Clipping Noise in OFDM Underwater Acoustic Communication System | |
CN107222438B (en) | Simplified BEM Channel Estimation Method for High Speed Mobile SC-FDMA System | |
CN101848183B (en) | Channel Estimation Method and Device in Multiple-Input Multiple-Output Orthogonal Frequency Division Multiplexing System | |
CN103428154B (en) | The transform domain multiplexing method of the dual-selection channel based on Vector OFDM | |
CN101242388A (en) | Channel Estimation Method for High Speed Single Carrier Frequency Domain Equalized UWB System | |
Kim et al. | Subspace-based channel estimation for MIMO-OFDM systems with few received blocks | |
CN108650005B (en) | Pilot structure and channel estimation method for MIMO-FBMC/OQAM system | |
CN103441980A (en) | Sideband-information-free shallow sea underwater acoustic communication pattern selection peak-to-average ratio restraining algorithm based on frequency reversal mirror technology | |
CN106341362B (en) | Pilot frequency transmission method, pilot frequency reception method and device thereof | |
CN102045285B (en) | Channel estimation method and device and communication system | |
CN103873406A (en) | Inter-frame interference elimination method used for underwater sound orthogonal frequency-division multiplexing communication system | |
CN102130860B (en) | Estimation method of two-dimensional discrete Fourier transform channel with phase compensation | |
CN104539562A (en) | MIMO-OFDM wideband HF channel estimation method | |
CN102025459B (en) | Nonparametric estimation ICA-based MIMO-OFDM system blind deconvolution method | |
CN102790746B (en) | Channel estimation method for OFDM (orthogonal frequency division multiplexing) system | |
CN106487738B (en) | A Selective Mapping Peak-to-Average Ratio Suppression Algorithm Based on Orthogonal Pilot Sequences for Underwater Acoustic OFDM Communication System | |
CN102487364B (en) | Channel estimation method and apparatus thereof | |
CN107592280A (en) | A kind of method that combination channel estimation avoids PTS transmitted sideband information | |
CN105119856B (en) | Single carrier frequency domain equalization method based on sub-block analysis | |
CN104967581B (en) | The channel estimation methods of single carrier frequency domain equalization under low energy disperse multipath channel |
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 | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20130206 Termination date: 20130406 |