CN112702093A - Channel estimation method in FDD downlink multi-user large-scale MIMO system - Google Patents

Channel estimation method in FDD downlink multi-user large-scale MIMO system Download PDF

Info

Publication number
CN112702093A
CN112702093A CN202011526747.2A CN202011526747A CN112702093A CN 112702093 A CN112702093 A CN 112702093A CN 202011526747 A CN202011526747 A CN 202011526747A CN 112702093 A CN112702093 A CN 112702093A
Authority
CN
China
Prior art keywords
channel
user
estimation method
steering vector
mimo system
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.)
Granted
Application number
CN202011526747.2A
Other languages
Chinese (zh)
Other versions
CN112702093B (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.)
Hangzhou Dianzi University
Original Assignee
Hangzhou Dianzi University
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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN202011526747.2A priority Critical patent/CN112702093B/en
Publication of CN112702093A publication Critical patent/CN112702093A/en
Application granted granted Critical
Publication of CN112702093B publication Critical patent/CN112702093B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B7/00Radio transmission systems, i.e. using radiation field
    • H04B7/02Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas
    • H04B7/04Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas
    • H04B7/0413MIMO systems
    • H04B7/0456Selection of precoding matrices or codebooks, e.g. using matrices antenna weighting
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L25/00Baseband systems
    • H04L25/02Details ; arrangements for supplying electrical power along data transmission lines
    • H04L25/0202Channel estimation
    • H04L25/024Channel estimation channel estimation algorithms
    • H04L25/0242Channel estimation channel estimation algorithms using matrix methods

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Power Engineering (AREA)
  • Radio Transmission System (AREA)

Abstract

本发明提供FDD下行多用户大规模MIMO系统中的信道估计方法,包括步骤:S1.针对多个用户设计训练序列以及处理接收信号;S2.根据接收信号进行相应的角度估计;S3.通过角度估计值构造信道导向矢量;S4.根据构造的信道导向矢量重构下行信道状态信息。本发明基于少量的、统一的训练序列的信道估计方法,多个用户可同时估计信道状态信息,有效减少导频开销,为实现更高的系统吞吐量打下基础。

Figure 202011526747

The present invention provides a channel estimation method in an FDD downlink multi-user massive MIMO system, which includes steps: S1. Designing training sequences for multiple users and processing received signals; S2. Performing corresponding angle estimation according to the received signals; S3. Using angle estimation value to construct a channel steering vector; S4. Reconstruct downlink channel state information according to the constructed channel steering vector. The present invention is based on the channel estimation method of a small number of uniform training sequences, multiple users can estimate the channel state information at the same time, effectively reduce the overhead of pilot frequency, and lay a foundation for realizing higher system throughput.

Figure 202011526747

Description

FDD下行多用户大规模MIMO系统中的信道估计方法Channel Estimation Method in FDD Downlink Multi-User Massive MIMO System

技术领域technical field

本发明涉及无线通信技术领域,尤其涉及FDD下行多用户大规模MIMO系统中的信道估计方法。The present invention relates to the technical field of wireless communication, in particular to a channel estimation method in an FDD downlink multi-user massive MIMO system.

背景技术Background technique

在大规模多输入多输出(MIMO)系统中,利用基站(BS)的大规模天线阵,系统能够获得极高的空间分辨率和空间分割复用增益,在没有受到严重干扰的情况下,还能同时为多个用户提供服务,极大的提高了系统的频谱效率和传送速率。但是这些优良性能都取决于信道状态信息(CSI)的可用性,因此在实际系统中如何获取完整的CSI对系统性能的影响至关重要。In a massive multiple-input multiple-output (MIMO) system, using the massive antenna array of the base station (BS), the system can obtain extremely high spatial resolution and spatial division multiplexing gain, and also can achieve high spatial resolution without severe interference. It can provide services for multiple users at the same time, which greatly improves the spectral efficiency and transmission rate of the system. However, these excellent performances all depend on the availability of channel state information (CSI), so how to obtain a complete CSI in an actual system is very important to the system performance.

针对频分双工(FDD)大规模MIMO系统,要获取完整的下行CSI需要花费大量的导频开销。这是由于在实际的大规模MIMO系统中,训练次数和反馈开销与BS天线的数量成正比,采用传统的线性信道估计方法(如最小二乘算法(LS)和线性最小均方误差(LMMSE)算法)获取完整的CSI的资源需要发送和基站天线数相同的训练序列,从而造成近年来发展起来的FDD下行信道估计方法,几乎都是针对不同用户的信道信息设计不同的训练序列获得CSI,这样,对于多用户系统,会消耗更多的时间和资源,不切实际。For a frequency division duplex (FDD) massive MIMO system, a large amount of pilot overhead is required to obtain complete downlink CSI. This is because in practical massive MIMO systems, the training times and feedback overhead are proportional to the number of BS antennas, and traditional linear channel estimation methods such as least squares algorithm (LS) and linear least mean square error (LMMSE) are used. Algorithm) To obtain complete CSI resources, it is necessary to send training sequences with the same number of base station antennas, resulting in FDD downlink channel estimation methods developed in recent years, almost all of which are designed for different users’ channel information. Different training sequences to obtain CSI, so that , for a multi-user system, it will consume more time and resources, which is impractical.

针对以上技术问题,故需对其进行改进。In view of the above technical problems, it is necessary to improve it.

发明内容SUMMARY OF THE INVENTION

本发明的目的是针对现有技术的缺陷,提供了FDD下行多用户大规模MIMO系统中的信道估计方法。The purpose of the present invention is to provide a channel estimation method in an FDD downlink multi-user massive MIMO system in view of the defects of the prior art.

为了实现以上目的,本发明采用以下技术方案:In order to achieve the above purpose, the present invention adopts the following technical solutions:

FDD下行多用户大规模MIMO系统中的信道估计方法,包括步骤:The channel estimation method in the FDD downlink multi-user massive MIMO system includes the steps:

S1.针对多个用户设计训练序列以及处理接收信号;S1. Design training sequences for multiple users and process received signals;

S2.根据接收信号进行相应的角度估计;S2. Perform corresponding angle estimation according to the received signal;

S3.通过角度估计值构造信道导向矢量;S3. Construct a channel steering vector through the angle estimation value;

S4.根据构造的信道导向矢量重构下行信道状态信息。S4. Reconstruct downlink channel state information according to the constructed channel steering vector.

进一步的,所述步骤S1中多个用户设计训练序列以及处理接收信号中包括K个相互独立的单天线用户、一个基站、基站处配备M根天线。Further, in the step S1, multiple users design training sequences and process received signals, including K mutually independent single-antenna users, one base station, and M antennas at the base station.

进一步的,所述基站处配备的M根天线采用平面阵排列方式,且M=MhMv,其中每行有Mh根天线,每列有Mv根天线,则第k(k=1,2,3,…,K)个用户从基站到达接收端的信道具体表示为:Further, the M antennas equipped at the base station are arranged in a plane array, and M=M h M v , wherein each row has M h antennas, and each column has M v antennas, then the kth (k=1 , 2, 3, ..., K) users' channels from the base station to the receiver are specifically expressed as:

Figure BDA0002850812930000021
Figure BDA0002850812930000021

其中:Pk表示基站与第k个用户的接收端之间的传播路径数;gk,l表示第k个用户的第l条路径的信道复增益;α(θ,φ)表示信道的导向矢量;α(θk,lk,l)表示第k个用户的第l条路径的信道的导向矢量;θk,l和φk,l分别表示第k个用户的第l条路径的仰角和方位角。Where: P k represents the number of propagation paths between the base station and the receiver of the kth user; gk ,l represents the channel complex gain of the lth path of the kth user; α(θ,φ) represents the channel steering vector; α(θ k,l ,φk ,l ) represents the steering vector of the channel of the lth path of the kth user; θk ,l and φk ,l respectively represent the lth path of the kth user of elevation and azimuth.

进一步的,所述信道的导向矢量α(θ,φ)表示为:Further, the steering vector α(θ,φ) of the channel is expressed as:

Figure BDA0002850812930000022
Figure BDA0002850812930000022

Figure BDA0002850812930000023
Figure BDA0002850812930000023

Figure BDA0002850812930000024
Figure BDA0002850812930000024

其中:d表示天线阵元之间的距离,假设d=λ/2,λ为载波波长;αv(θ)中包含平面阵列的纵向角度信息,αh(θ,φ)包含平面阵列的横向角度信息。Among them: d represents the distance between the antenna elements, assuming d=λ/2, λ is the carrier wavelength; α v (θ) contains the longitudinal angle information of the planar array, and α h (θ, φ) contains the horizontal direction of the planar array. angle information.

进一步的,所述步骤S1中设计训练序列具体为:Further, the training sequence designed in the step S1 is specifically:

假设基站对所有的K个用户的发送信号均为SN,为M×N维的矩阵,具体表示为:It is assumed that the transmitted signals of the base station to all K users are S N , which is an M×N-dimensional matrix, which is specifically expressed as:

Figure BDA0002850812930000031
Figure BDA0002850812930000031

其中:Fn表示n维的离散傅里叶DFT矩阵;0n为n维的零矩阵;0(M-7n)×n为(M-7n)×n维的零矩阵,n满足4n=N,N为采样总数。Where: F n represents an n-dimensional discrete Fourier DFT matrix; 0 n is an n-dimensional zero matrix; 0 (M-7n)×n is a (M-7n)×n-dimensional zero matrix, and n satisfies 4n=N , and N is the total number of samples.

进一步的,所述步骤S1中处理接收信号具体为:Further, the processing of the received signal in the step S1 is specifically:

混入噪声后,第k(k=1,2,3,..,K)个用户端的接收信号表示为:After mixing with noise, the received signal of the kth (k=1, 2, 3, .., K) user terminal is expressed as:

Figure BDA0002850812930000032
Figure BDA0002850812930000032

其中:ρ为信噪比,zk为第k个用户处混入的噪声,假设是服从零均值,单位方差的高斯噪声;Among them: ρ is the signal-to-noise ratio, z k is the noise mixed in the kth user, assuming that it is Gaussian noise that obeys zero mean and unit variance;

对发送信号

Figure BDA0002850812930000035
进行反DFT变换,则在第k个用户处的实际接收信号表示为:to send a signal
Figure BDA0002850812930000035
Perform inverse DFT transform, then the actual received signal at the kth user is expressed as:

Figure BDA0002850812930000033
Figure BDA0002850812930000033

其中,yk表示接收信号。where y k represents the received signal.

进一步的,所述步骤S2具体包括:Further, the step S2 specifically includes:

S21.设定码本,将信道导向矢量α(θ,φ)中横向矢量αh(θ,φ)中的cosθsinφ和纵向矢量αv(θ)中的sinθ分别看成区间在(-1,1)中的两个整体;则在(-1,1)区间内分别创建两个m份的码本,且m=M/2;S21. Set the codebook, and regard the cosθsinφ in the horizontal vector α h (θ, φ) in the channel steering vector α (θ, φ) and the sinθ in the longitudinal vector α v (θ) as the interval in (-1, 1) two wholes; then create two m codebooks respectively in the (-1,1) interval, and m=M/2;

Figure BDA0002850812930000034
Figure BDA0002850812930000034

Figure BDA0002850812930000041
Figure BDA0002850812930000041

S22.第k个用户处的i(i=1,2,…,P)条路径的接收信号为yk(i),推出第k个用户的第l条路径的噪声,表示为:S22. The received signal of the i (i=1, 2, .

Figure BDA0002850812930000042
Figure BDA0002850812930000042

其中,

Figure BDA0002850812930000043
Figure BDA0002850812930000044
表示每次遍历设定的角度码本的取值;
Figure BDA0002850812930000045
表示每次遍历设定的角度码本后的取值;in,
Figure BDA0002850812930000043
and
Figure BDA0002850812930000044
Indicates the value of the angle codebook set for each traversal;
Figure BDA0002850812930000045
Indicates the value after each traversal of the set angle codebook;

S23.结合预先设定的码本,将码字设定如下:S23. In combination with the preset codebook, the codeword is set as follows:

Figure BDA0002850812930000046
Figure BDA0002850812930000046

结合设定的码字

Figure BDA0002850812930000047
遍历预先设定的码本,从yk,r(i)中挑选出投影功率最大的码字,表示为:Combining set codewords
Figure BDA0002850812930000047
Traverse the preset codebook, and select the codeword with the largest projection power from y k,r (i), which is expressed as:

Figure BDA0002850812930000048
Figure BDA0002850812930000048

码字中的

Figure BDA0002850812930000049
Figure BDA00028508129300000410
为第k个用户的第i条路径中的信道导向矢量中的角度信息的估计值,利用得到的角度信息的估计值,重新弄构造信道的导向矢量
Figure BDA00028508129300000411
in the codeword
Figure BDA0002850812930000049
and
Figure BDA00028508129300000410
is the estimated value of the angle information in the channel steering vector in the ith path of the kth user, and reconstructs the channel steering vector using the obtained estimated value of the angle information
Figure BDA00028508129300000411

S24.结合预先设定的码本,遍历角度码本,挑选得到角度信息的估计值,从而进行对第k个用户的第i条路径增益的估计,表示为:S24. Combine the preset codebook, traverse the angle codebook, select and obtain the estimated value of the angle information, so as to estimate the gain of the i-th path of the k-th user, which is expressed as:

Figure BDA00028508129300000412
Figure BDA00028508129300000412

则信道复增益的估计值为:

Figure BDA00028508129300000413
Then the estimated value of the channel complex gain is:
Figure BDA00028508129300000413

进一步的,所述步骤S3具体为:Further, the step S3 is specifically:

结合角度信息估计值

Figure BDA0002850812930000051
通过信道导向矢量α(θ,φ)中的αv(θ)和αh(θ,φ)的形式构造出导向矩阵横向矢量
Figure BDA0002850812930000052
和纵向矢量
Figure BDA0002850812930000053
二者做张量积
Figure BDA0002850812930000054
构造出信道导向矢量
Figure BDA0002850812930000055
Combined with angle information estimates
Figure BDA0002850812930000051
The steering matrix transverse vector is constructed in the form of α v (θ) and α h (θ, φ) in the channel steering vector α (θ, φ).
Figure BDA0002850812930000052
and portrait vector
Figure BDA0002850812930000053
The two do tensor product
Figure BDA0002850812930000054
Construct the channel steering vector
Figure BDA0002850812930000055

进一步的,所述步骤S4具体为:根据重构的信道导向矢量

Figure BDA0002850812930000056
以及信道复增益估计值
Figure BDA0002850812930000057
重构信道矩阵从而恢复下行信道信息。Further, the step S4 is specifically: according to the reconstructed channel steering vector
Figure BDA0002850812930000056
and the channel complex gain estimate
Figure BDA0002850812930000057
The channel matrix is reconstructed to restore downlink channel information.

与现有技术相比,本发明基于少量的、统一的训练序列的信道估计方法,多个用户可同时估计信道状态信息,有效减少导频开销,为实现更高的系统吞吐量打下基础。Compared with the prior art, the present invention is based on the channel estimation method based on a small number of uniform training sequences, and multiple users can estimate the channel state information at the same time, effectively reducing pilot frequency overhead, and laying a foundation for realizing higher system throughput.

附图说明Description of drawings

图1为实施例一提供的FDD下行多用户大规模MIMO系统中的信道估计方法流程图;1 is a flowchart of a channel estimation method in an FDD downlink multi-user massive MIMO system according to Embodiment 1;

图2为实施例一提供的FDD下行多用户大规模MIMO系统中的信道估计方法仿真图。FIG. 2 is a simulation diagram of a channel estimation method in an FDD downlink multi-user massive MIMO system according to the first embodiment.

具体实施方式Detailed ways

以下通过特定的具体实例说明本发明的实施方式,本领域技术人员可由本说明书所揭露的内容轻易地了解本发明的其他优点与功效。本发明还可以通过另外不同的具体实施方式加以实施或应用,本说明书中的各项细节也可以基于不同观点与应用,在没有背离本发明的精神下进行各种修饰或改变。需说明的是,在不冲突的情况下,以下实施例及实施例中的特征可以相互组合。The embodiments of the present invention are described below through specific specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and various details in this specification can also be modified or changed based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that the following embodiments and features in the embodiments may be combined with each other under the condition of no conflict.

实施例一Example 1

本实施例提供FDD下行多用户大规模MIMO系统中的信道估计方法,如图1所示,包括:This embodiment provides a channel estimation method in an FDD downlink multi-user massive MIMO system, as shown in FIG. 1 , including:

S1.针对多个用户设计训练序列以及处理接收信号;S1. Design training sequences for multiple users and process received signals;

S2.根据接收信号进行相应的角度估计;S2. Perform corresponding angle estimation according to the received signal;

S3.通过角度估计值构造信道导向矢量;S3. Construct a channel steering vector through the angle estimation value;

S4.根据构造的信道导向矢量重构下行信道状态信息。S4. Reconstruct downlink channel state information according to the constructed channel steering vector.

本实施例的FDD下行多用户大规模MIMO系统中的信道估计方法中,基于少量的、统一的训练序列的信道估计方法,用于下行多用户大规模MIMO系统,针对现有信道估计方法的不足进行了改进。In the channel estimation method in the FDD downlink multi-user massive MIMO system in this embodiment, the channel estimation method based on a small number of uniform training sequences is used in the downlink multi-user massive MIMO system, aiming at the shortcomings of the existing channel estimation methods. Improvements have been made.

具体应用案例如下:The specific application cases are as follows:

假设系统中有1个小区,1个基站,基站的天线数为128根,采样次数N取32,n取8,m取64。无线信号从基站到达每个用户的路径数Pk=6,每个用户接收天线为1个。平面阵列天线每行的天线数Mh=16,每列的天线数Mv=8,载波波长λ=3×108/1.5×109,阵元间距d=λ/2,信道矩阵中导向矢量中包含的仰角θ和方位角φ都是服从区间(-900,900]内的均匀分布。Assume that there is 1 cell and 1 base station in the system, the number of antennas of the base station is 128, the sampling times N is 32, n is 8, and m is 64. The number of paths of wireless signals from the base station to each user is P k =6, and each user has one receiving antenna. The number of antennas in each row of the planar array antenna M h = 16, the number of antennas in each column M v = 8, the carrier wavelength λ = 3 × 10 8 /1.5 × 10 9 , the array element spacing d = λ/2, the channel matrix is guided The elevation angle θ and azimuth angle φ contained in the vector are uniformly distributed in the interval (-900,900].

在步骤S1中,针对多个用户设计训练序列以及处理接收信号。In step S1, training sequences are designed for multiple users and received signals are processed.

基站处配备的128根天线采用平面阵(UPA Uniform-Planar-Array)排列方式,且M=16×8,其中每行有16根天线,每列有8根天线。第k个用户从基站到达接收端的信道具体表示为:The 128 antennas equipped at the base station are arranged in a planar array (UPA Uniform-Planar-Array), and M=16×8, in which there are 16 antennas in each row and 8 antennas in each column. The channel of the kth user from the base station to the receiver is specifically expressed as:

Figure BDA0002850812930000061
Figure BDA0002850812930000061

其中:Pk为基站与第k个用户的接收端之间的传播路径数,gk,l表示第k个用户的第l条路径的信道复增益,假设是服从零均值,单位方差的复高斯分布。α(θk,lk,l)表示第k个用户的第l条路径的信道的导向矢量,θk,l和φk,l分别表示第k个用户的第l条路径的仰角和方位角。Among them: P k is the number of propagation paths between the base station and the receiving end of the kth user, gk,l represents the channel complex gain of the lth path of the kth user, assuming that it obeys the zero mean, the unit variance complex Gaussian distribution. α(θ k,l , φ k,l ) represents the steering vector of the channel of the lth path of the kth user, θk ,l and φk ,l respectively represent the elevation angle of the lth path of the kth user and azimuth.

设计训练序列具体为:The design training sequence is as follows:

对所有的K个用户,在基站端选取训练序列的依据是:应当满足对第k个用户的第l条路径的信道导向矢量中的横向角度信息和纵向角度信息的采样次数满足平面阵列天线的行数和列数之比。For all K users, the basis for selecting the training sequence at the base station is: the sampling times of the horizontal angle information and the vertical angle information in the channel steering vector of the lth path of the kth user should satisfy the requirements of the planar array antenna. The ratio of the number of rows to the number of columns.

假设基站对所有的K个用户的发送信号均为SN,为128×32维的矩阵,具体表示为:Assuming that the transmitted signals of the base station to all K users are S N , which is a 128×32-dimensional matrix, which is specifically expressed as:

Figure BDA0002850812930000071
Figure BDA0002850812930000071

其中:Fn为8维的离散傅里叶(DFT)矩阵,0n为8维的零矩阵,0(M-7n)×n为72×8维的零矩阵,n满足4n=N,N为采样总数。Where: F n is an 8-dimensional discrete Fourier (DFT) matrix, 0 n is an 8-dimensional zero matrix, 0 (M-7n)×n is a 72×8-dimensional zero matrix, and n satisfies 4n=N, N is the total number of samples.

处理接收信号具体为:The processing of the received signal is as follows:

混入噪声后,第k个用户端的接收信号表示为:After mixing noise, the received signal of the kth user terminal is expressed as:

Figure BDA0002850812930000072
Figure BDA0002850812930000072

其中:ρ为信噪比,zk为第k个用户处混入的噪声,假设是服从零均值,单位方差的高斯噪声。Among them: ρ is the signal-to-noise ratio, z k is the noise mixed in at the kth user, assuming that it is Gaussian noise with zero mean and unit variance.

对发送信号

Figure BDA0002850812930000073
进行反DFT变换,则在第k个用户处的实际接收信号表示为:to send a signal
Figure BDA0002850812930000073
Perform inverse DFT transform, then the actual received signal at the kth user is expressed as:

Figure BDA0002850812930000074
Figure BDA0002850812930000074

其中,yk表示接收信号。where y k represents the received signal.

在步骤S2中,根据接收信号进行相应的角度估计。In step S2, corresponding angle estimation is performed according to the received signal.

预先创建一个关于角度的码本,其次遍历所有码本。A codebook of angles is created in advance, and then all codebooks are traversed.

S21.设定码本,将信道导向矢量α(θ,φ)中横向矢量αh(θ,φ)中的cosθsinφ和纵向矢量αv(θ)中的sinθ分别看成区间在(-1,1)中的两个整体。则在(-1,1)区间内分别创建两个64份的码本。S21. Set the codebook, and regard the cosθsinφ in the horizontal vector α h (θ, φ) in the channel steering vector α (θ, φ) and the sinθ in the longitudinal vector α v (θ) as the interval in (-1, 1) of the two wholes. Then, two 64 codebooks are created in the (-1,1) interval.

Figure BDA0002850812930000075
Figure BDA0002850812930000075

Figure BDA0002850812930000081
Figure BDA0002850812930000081

S22.第k个用户处的i(i=1,2,3,4,5,6)条路径的接收信号为yk(i),推出第k个用户的第l条路径的噪声表示如下:S22. The received signal of the i (i=1, 2, 3, 4, 5, 6) path at the kth user is y k (i), and the noise of the lth path of the kth user is deduced as follows :

Figure BDA0002850812930000082
Figure BDA0002850812930000082

其中:

Figure BDA0002850812930000083
Figure BDA0002850812930000084
表示每次遍历设定的角度码本的取值,
Figure BDA0002850812930000085
表示每次遍历设定的角度码本后的取值。in:
Figure BDA0002850812930000083
and
Figure BDA0002850812930000084
Indicates the value of the angle codebook set for each traversal,
Figure BDA0002850812930000085
Indicates the value after each traversal of the set angle codebook.

S23.结合预先设定的码本,将码字设定如下:S23. In combination with the preset codebook, the codeword is set as follows:

Figure BDA0002850812930000086
Figure BDA0002850812930000086

结合设定的码字

Figure BDA0002850812930000087
遍历预先设定的码本,从yk,r(i)中挑选出投影功率最大的码字,具体计算公式如下:Combining set codewords
Figure BDA0002850812930000087
Traverse the preset codebook, and select the codeword with the largest projection power from y k,r (i). The specific calculation formula is as follows:

Figure BDA0002850812930000088
Figure BDA0002850812930000088

则码字中的

Figure BDA0002850812930000089
Figure BDA00028508129300000810
就是第k个用户的第i条路径中的信道导向矢量中的角度信息的估计值,利用得到的角度信息的估计值,就可以重新弄构造信道的导向矢量
Figure BDA00028508129300000811
in the codeword
Figure BDA0002850812930000089
and
Figure BDA00028508129300000810
is the estimated value of the angle information in the channel steering vector in the i-th path of the kth user. Using the obtained estimated value of the angle information, the channel steering vector can be reconstructed
Figure BDA00028508129300000811

S24.结合预先设定的码本,遍历角度码本,挑选得到角度信息的估计值,从而进行对第k个用户的第i条路径增益的估计,具体计算公式如下:S24. Combine the preset codebook, traverse the angle codebook, select and obtain the estimated value of the angle information, so as to estimate the gain of the i-th path of the k-th user, and the specific calculation formula is as follows:

Figure BDA00028508129300000812
Figure BDA00028508129300000812

则信道复增益的估计值为:

Figure BDA00028508129300000813
Then the estimated value of the channel complex gain is:
Figure BDA00028508129300000813

在步骤S3中,通过角度估计值构造信道导向矢量。In step S3, a channel steering vector is constructed from the angle estimates.

结合角度信息估计值

Figure BDA0002850812930000091
通过信道导向矢量α(θ,φ)中的αv(θ)和αh(θ,φ)的形式构造出导向矩阵横向矢量
Figure BDA0002850812930000092
和纵向矢量
Figure BDA0002850812930000093
二者做张量积
Figure BDA0002850812930000094
可构造出信道导向矢量
Figure BDA0002850812930000095
Combined with angle information estimates
Figure BDA0002850812930000091
The steering matrix transverse vector is constructed in the form of α v (θ) and α h (θ, φ) in the channel steering vector α (θ, φ).
Figure BDA0002850812930000092
and portrait vector
Figure BDA0002850812930000093
The two do tensor product
Figure BDA0002850812930000094
The channel steering vector can be constructed
Figure BDA0002850812930000095

在步骤S4中,根据构造的信道导向矢量重构下行信道状态信息。In step S4, the downlink channel state information is reconstructed according to the constructed channel steering vector.

根据重构的信道导向矢量

Figure BDA0002850812930000096
以及信道复增益估计值
Figure BDA0002850812930000097
重构信道矩阵从而恢复下行信道信息。According to the reconstructed channel steering vector
Figure BDA0002850812930000096
and the channel complex gain estimate
Figure BDA0002850812930000097
The channel matrix is reconstructed to restore downlink channel information.

在本实施例中,为了分析所提出的FDD大规模天线系统中基于角度估计的信道估计方法的性能,定义了系统吞吐量为:C=log(1+ρ||hHf||2),f为预编码。In this embodiment, in order to analyze the performance of the proposed channel estimation method based on angle estimation in the FDD large-scale antenna system, the system throughput is defined as: C=log(1+ρ||h H f|| 2 ) , f is precoding.

如图2所示,是在上述举例的条件下,关于系统吞吐量的仿真图,其中“理想转态”是采用真实信道做预编码的系统吞吐量曲线,本实施例提出的方法采用重构信道做预编码的系统吞吐量曲线。从图2可以看出,本实施例提出的方法得到的信道估计值所做预编码,在同等条件下与真实信道所做预编码相比,系统的吞吐量相差都小于1个比特左右,且信噪比对本发明提出的方法影响不大,具有很好的系统性能。As shown in Figure 2, it is a simulation diagram of system throughput under the conditions of the above example, in which the "ideal transition" is the system throughput curve using the real channel for precoding, and the method proposed in this embodiment uses reconstruction System throughput curve for channel precoding. It can be seen from FIG. 2 that the precoding of the channel estimation value obtained by the method proposed in this embodiment is compared with the precoding of the real channel under the same conditions, and the difference in the throughput of the system is less than about 1 bit, and The signal-to-noise ratio has little influence on the method proposed by the present invention, and has good system performance.

与现有技术相比,本实施例为了更好的节省训练开销,在经过理论推导后,结合PRONY方程可以推出,对于大规模天线系统,针对阵列天线,其信道都是满足PRONY方程的形式。基于此,考虑到在多用户场景下,可以对不同的用户发送相同的训练序列同时进行训练,并结合角度估计方法获取完整的信道状态信息依此来减少导频开销和时间消耗。基于少量的、统一的训练序列的信道估计方法,多个用户可同时估计信道状态信息,有效减少导频开销,为实现更高的系统吞吐量打下基础。Compared with the prior art, in this embodiment, in order to better save training overhead, after theoretical derivation, it can be deduced in combination with the PRONY equation. For a large-scale antenna system, for an array antenna, the channels are all in a form that satisfies the PRONY equation. Based on this, considering that in a multi-user scenario, different users can be trained by sending the same training sequence at the same time, and combined with the angle estimation method to obtain complete channel state information to reduce pilot overhead and time consumption. The channel estimation method based on a small number of uniform training sequences allows multiple users to estimate channel state information simultaneously, effectively reducing pilot overhead and laying the foundation for achieving higher system throughput.

本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。The specific embodiments described herein are merely illustrative of the spirit of the invention. Those skilled in the art to which the present invention pertains can make various modifications or additions to the described specific embodiments or substitute in similar manners, but will not deviate from the spirit of the present invention or go beyond the definitions of the appended claims range.

Claims (9)

1.FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,包括步骤:1. The channel estimation method in the FDD downlink multi-user massive MIMO system, is characterized in that, comprises the steps: S1.针对多个用户设计训练序列以及处理接收信号;S1. Design training sequences for multiple users and process received signals; S2.根据接收信号进行相应的角度估计;S2. Perform corresponding angle estimation according to the received signal; S3.通过角度估计值构造信道导向矢量;S3. Construct a channel steering vector through the angle estimation value; S4.根据构造的信道导向矢量重构下行信道状态信息。S4. Reconstruct downlink channel state information according to the constructed channel steering vector. 2.根据权利要求1所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述步骤S1中多个用户设计训练序列以及处理接收信号中包括K个相互独立的单天线用户、一个基站、基站处配备M根天线。2. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 1, wherein, in the step S1, a plurality of users design training sequences and process the received signal, including K mutually independent single Antenna users, one base station, and M antennas are provided at the base station. 3.根据权利要求2所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述基站处配备的M根天线采用平面阵排列方式,且M=MhMv,其中每行有Mh根天线,每列有Mv根天线,则第k(k=1,2,3,…,K)个用户从基站到达接收端的信道具体表示为:3. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 2, wherein the M antennas equipped at the base station adopt a planar array arrangement, and M=M h M v , There are M h antennas in each row and M v antennas in each column, then the channel of the kth (k=1,2,3,...,K) user from the base station to the receiving end is specifically expressed as:
Figure FDA0002850812920000011
Figure FDA0002850812920000011
其中:Pk表示基站与第k个用户的接收端之间的传播路径数;gk,l表示第k个用户的第l条路径的信道复增益;α(θ,φ)表示信道的导向矢量;α(θk,lk,l)表示第k个用户的第l条路径的信道的导向矢量;θk,l和φk,l分别表示第k个用户的第l条路径的仰角和方位角。Where: P k represents the number of propagation paths between the base station and the receiver of the kth user; gk ,l represents the channel complex gain of the lth path of the kth user; α(θ,φ) represents the channel steering vector; α(θ k,lk,l ) represents the steering vector of the channel of the lth path of the kth user; θk ,l and φk ,l respectively represent the lth path of the kth user of elevation and azimuth.
4.根据权利要求3所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述信道的导向矢量α(θ,φ)表示为:4. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 3, wherein the steering vector α(θ, φ) of the channel is expressed as:
Figure FDA0002850812920000012
Figure FDA0002850812920000012
Figure FDA0002850812920000013
Figure FDA0002850812920000013
Figure FDA0002850812920000021
Figure FDA0002850812920000021
其中:d表示天线阵元之间的距离,假设d=λ/2,λ为载波波长;αv(θ)中包含平面阵列的纵向角度信息,αh(θ,φ)包含平面阵列的横向角度信息。Among them: d represents the distance between the antenna elements, assuming d=λ/2, λ is the carrier wavelength; α v (θ) contains the longitudinal angle information of the planar array, and α h (θ, φ) contains the horizontal direction of the planar array. angle information.
5.根据权利要求4所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述步骤S1中设计训练序列具体为:5. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 4, wherein the design training sequence in the step S1 is specifically: 假设基站对所有的K个用户的发送信号均为SN,为M×N维的矩阵,具体表示为:It is assumed that the transmitted signals of the base station to all K users are S N , which is an M×N-dimensional matrix, which is specifically expressed as:
Figure FDA0002850812920000022
Figure FDA0002850812920000022
其中:Fn表示n维的离散傅里叶DFT矩阵;0n为n维的零矩阵;0(M-7n)×n为(M-7n)×n维的零矩阵,n满足4n=N,N为采样总数。Where: F n represents an n-dimensional discrete Fourier DFT matrix; 0 n is an n-dimensional zero matrix; 0 (M-7n)×n is a (M-7n)×n-dimensional zero matrix, and n satisfies 4n=N , and N is the total number of samples.
6.根据权利要求5所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述步骤S1中处理接收信号具体为:6. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 5, wherein the processing of the received signal in the step S1 is specifically: 混入噪声后,第k(k=1,2,3,..,K)个用户端的接收信号表示为:After mixing with noise, the received signal of the kth (k=1, 2, 3, .., K) user terminal is expressed as:
Figure FDA0002850812920000023
Figure FDA0002850812920000023
其中:ρ为信噪比,zk为第k个用户处混入的噪声,假设是服从零均值,单位方差的高斯噪声;Among them: ρ is the signal-to-noise ratio, z k is the noise mixed in the kth user, assuming that it is Gaussian noise that obeys zero mean and unit variance; 对发送信号
Figure FDA0002850812920000024
进行反DFT变换,则在第k个用户处的实际接收信号表示为:
to send a signal
Figure FDA0002850812920000024
Perform inverse DFT transform, then the actual received signal at the kth user is expressed as:
Figure FDA0002850812920000031
Figure FDA0002850812920000031
其中,yk表示接收信号。where y k represents the received signal.
7.根据权利要求6所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述步骤S2具体包括:7. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 6, wherein the step S2 specifically comprises: S21.设定码本,将信道导向矢量α(θ,φ)中横向矢量αh(θ,φ)中的cosθsinφ和纵向矢量αv(θ)中的sinθ分别看成区间在(-1,1)中的两个整体;则在(-1,1)区间内分别创建两个m份的码本,且m=M/2;S21. Set the codebook, and regard the cosθsinφ in the horizontal vector α h (θ, φ) in the channel steering vector α (θ, φ) and the sinθ in the longitudinal vector α v (θ) as the interval in (-1, 1) two wholes; then create two m codebooks respectively in the (-1,1) interval, and m=M/2;
Figure FDA0002850812920000032
Figure FDA0002850812920000032
Figure FDA0002850812920000033
Figure FDA0002850812920000033
S22.第k个用户处的i(i=1,2,…,P)条路径的接收信号为yk(i),推出第k个用户的第l条路径的噪声,表示为:S22. The received signal of the i (i=1, 2, .
Figure FDA0002850812920000034
Figure FDA0002850812920000034
其中,
Figure FDA0002850812920000035
Figure FDA0002850812920000036
表示每次遍历设定的角度码本的取值;
Figure FDA0002850812920000037
表示每次遍历设定的角度码本后的取值;
in,
Figure FDA0002850812920000035
and
Figure FDA0002850812920000036
Indicates the value of the angle codebook set for each traversal;
Figure FDA0002850812920000037
Indicates the value after each traversal of the set angle codebook;
S23.结合预先设定的码本,将码字设定如下:S23. In combination with the preset codebook, the codeword is set as follows:
Figure FDA0002850812920000038
Figure FDA0002850812920000038
结合设定的码字
Figure FDA0002850812920000039
遍历预先设定的码本,从yk,r(i)中挑选出投影功率最大的码字,表示为:
Combining set codewords
Figure FDA0002850812920000039
Traverse the preset codebook, and select the codeword with the largest projection power from y k,r (i), which is expressed as:
Figure FDA00028508129200000310
Figure FDA00028508129200000310
码字中的
Figure FDA00028508129200000311
Figure FDA00028508129200000312
为第k个用户的第i条路径中的信道导向矢量中的角度信息的估计值,利用得到的角度信息的估计值,重新弄构造信道的导向矢量
Figure FDA0002850812920000041
in the codeword
Figure FDA00028508129200000311
and
Figure FDA00028508129200000312
is the estimated value of the angle information in the channel steering vector in the ith path of the kth user, and reconstructs the channel steering vector using the obtained estimated value of the angle information
Figure FDA0002850812920000041
S24.结合预先设定的码本,遍历角度码本,挑选得到角度信息的估计值,从而进行对第k个用户的第i条路径增益的估计,表示为:S24. Combine the preset codebook, traverse the angle codebook, select and obtain the estimated value of the angle information, so as to estimate the gain of the i-th path of the k-th user, which is expressed as:
Figure FDA0002850812920000042
Figure FDA0002850812920000042
则信道复增益的估计值为:
Figure FDA0002850812920000043
Then the estimated value of the channel complex gain is:
Figure FDA0002850812920000043
8.根据权利要求7所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述步骤S3具体为:8. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 7, wherein the step S3 is specifically: 结合角度信息估计值
Figure FDA0002850812920000044
通过信道导向矢量α(θ,φ)中的αv(θ)和αh(θ,φ)的形式构造出导向矩阵横向矢量
Figure FDA0002850812920000045
和纵向矢量
Figure FDA0002850812920000046
二者做张量积
Figure FDA0002850812920000047
构造出信道导向矢量
Figure FDA0002850812920000048
Combined with angle information estimates
Figure FDA0002850812920000044
The steering matrix transverse vector is constructed in the form of α v (θ) and α h (θ, φ) in the channel steering vector α (θ, φ).
Figure FDA0002850812920000045
and portrait vector
Figure FDA0002850812920000046
The two do tensor product
Figure FDA0002850812920000047
Construct the channel steering vector
Figure FDA0002850812920000048
9.根据权利要求8所述的FDD下行多用户大规模MIMO系统中的信道估计方法,其特征在于,所述步骤S4具体为:根据重构的信道导向矢量
Figure FDA0002850812920000049
以及信道复增益估计值
Figure FDA00028508129200000410
重构信道矩阵从而恢复下行信道信息。
9. The channel estimation method in the FDD downlink multi-user massive MIMO system according to claim 8, wherein the step S4 is specifically: according to the reconstructed channel steering vector
Figure FDA0002850812920000049
and the channel complex gain estimate
Figure FDA00028508129200000410
The channel matrix is reconstructed to restore downlink channel information.
CN202011526747.2A 2020-12-22 2020-12-22 Channel estimation method in FDD downlink multi-user large-scale MIMO system Expired - Fee Related CN112702093B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011526747.2A CN112702093B (en) 2020-12-22 2020-12-22 Channel estimation method in FDD downlink multi-user large-scale MIMO system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011526747.2A CN112702093B (en) 2020-12-22 2020-12-22 Channel estimation method in FDD downlink multi-user large-scale MIMO system

Publications (2)

Publication Number Publication Date
CN112702093A true CN112702093A (en) 2021-04-23
CN112702093B CN112702093B (en) 2021-11-30

Family

ID=75510202

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011526747.2A Expired - Fee Related CN112702093B (en) 2020-12-22 2020-12-22 Channel estimation method in FDD downlink multi-user large-scale MIMO system

Country Status (1)

Country Link
CN (1) CN112702093B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114244658A (en) * 2021-12-24 2022-03-25 杭州电子科技大学 A Channel Estimation Method Based on Double-Angle Estimation in Massive MIMO Systems

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160087699A1 (en) * 2004-08-11 2016-03-24 Interdigital Technology Corporation Channel sounding for improved system performance
CN106302274A (en) * 2016-08-26 2017-01-04 清华大学 A kind of extensive mimo system multiuser channel is estimated and tracking
CN106453163A (en) * 2016-10-11 2017-02-22 电子科技大学 Massive MIMO (Multiple Input Multiple Output) channel estimation method
CN108259072A (en) * 2018-01-10 2018-07-06 杭州电子科技大学 The method for reducing training sequence expense for the extensive MIMO downlink systems of FDD
CN109327918A (en) * 2018-11-19 2019-02-12 东南大学 A low-overhead FDD massive MIMO downlink channel reconstruction method
CN109802901A (en) * 2017-11-16 2019-05-24 上海交通大学 Based on the 3D mimo channel estimation method and system for reaching angular measurement
CN110086734A (en) * 2019-03-14 2019-08-02 西安电子科技大学 Channel tracking and restoration methods in a kind of extensive mimo system of time-varying

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160087699A1 (en) * 2004-08-11 2016-03-24 Interdigital Technology Corporation Channel sounding for improved system performance
CN106302274A (en) * 2016-08-26 2017-01-04 清华大学 A kind of extensive mimo system multiuser channel is estimated and tracking
CN106453163A (en) * 2016-10-11 2017-02-22 电子科技大学 Massive MIMO (Multiple Input Multiple Output) channel estimation method
CN109802901A (en) * 2017-11-16 2019-05-24 上海交通大学 Based on the 3D mimo channel estimation method and system for reaching angular measurement
CN108259072A (en) * 2018-01-10 2018-07-06 杭州电子科技大学 The method for reducing training sequence expense for the extensive MIMO downlink systems of FDD
CN109327918A (en) * 2018-11-19 2019-02-12 东南大学 A low-overhead FDD massive MIMO downlink channel reconstruction method
CN110086734A (en) * 2019-03-14 2019-08-02 西安电子科技大学 Channel tracking and restoration methods in a kind of extensive mimo system of time-varying

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
HAYDER ALMOSA; SUSANNA MOSLEH: "Downlink Channel Estimation with Limited Feedback for FDD Multi-User Massive MIMO with Spatial Channel Correlation", 《2018 IEEE INTERNATIONAL CONFERENCE ON COMMUNICATIONS (ICC)》 *
周蕾: "FDD大规模MIMO系统基于角度域稀疏性和互易性的信道估计与传输方法研究", 《中国优秀硕士学位论文全文数据库 (基础科学辑)》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114244658A (en) * 2021-12-24 2022-03-25 杭州电子科技大学 A Channel Estimation Method Based on Double-Angle Estimation in Massive MIMO Systems
CN114244658B (en) * 2021-12-24 2024-03-26 杭州电子科技大学 Channel estimation method based on multiple angle estimation in large-scale MIMO system

Also Published As

Publication number Publication date
CN112702093B (en) 2021-11-30

Similar Documents

Publication Publication Date Title
CN108933745B (en) Broadband channel estimation method based on super-resolution angle and time delay estimation
Love et al. Limited feedback precoding for spatial multiplexing systems
Xie et al. Geometric mean decomposition based hybrid precoding for millimeter-wave massive MIMO
EP2737640B1 (en) Apparatus and method for combining baseband processing and radio frequency beam steering in a wireless communication system
US8477872B2 (en) Multi-antenna system and method for transmitting and receiving information in multi-antenna system
Mao et al. Over-sampling codebook-based hybrid minimum sum-mean-square-error precoding for millimeter-wave 3D-MIMO
CN102546088B (en) A kind of block diagonalization method for precoding and device
CN107332596B (en) Zero forcing-based millimeter wave communication system hybrid precoding method
CN101232356A (en) Precoding method, system and apparatus in MIMO system
CN110943768B (en) Mixed precoding codebook joint design method of millimeter wave large-scale MIMO system
CN104113398B (en) MIMO blind Channel Estimation fuzziness minimizing technologies based on Orthogonal Space-Time Block Code
CN108063634B (en) An optimal regularized precoding method for low-precision quantized massive MIMO
CN109547082A (en) Mixing precoding optimization method based on the extensive antenna system of millimeter wave
Mahmood et al. 3-D antenna array structures for millimeter wave multi-user massive MIMO hybrid precoder design: A performance comparison
CN104243106A (en) Codebook construction method based of 3D (three dimensional) cross polarization channel based massive multiple input multiple output (MIMO) system
Yin et al. A partial channel reciprocity-based codebook for wideband FDD massive MIMO
CN107707284B (en) A Hybrid Precoding Method Based on Channel Statistics Codebook Quantization Feedback
CN112702092B (en) Channel estimation method in FDD downlink multi-user large-scale MIMO system
CN112702093B (en) Channel estimation method in FDD downlink multi-user large-scale MIMO system
CN102082635B (en) Codebook storage and information transceiving method and device using codebook
CN110943767A (en) Precoding design method based on channel part reciprocity in FDD large-scale MIMO system
Khorsandmanesh et al. Optimized precoding for MU-MIMO with fronthaul quantization
CN114338303B (en) Channel estimation method and system based on multidimensional Hankel matrix in large-scale MIMO system
CN105356922A (en) Nonuniform segmentation-based vertical-dimension codebook design method under 3D MIMO channel
CN114598574A (en) Millimeter wave channel estimation method based on deep learning

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211130

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