CN110753365B - Interference coordination method for heterogeneous cellular networks - Google Patents

Interference coordination method for heterogeneous cellular networks Download PDF

Info

Publication number
CN110753365B
CN110753365B CN201911042510.4A CN201911042510A CN110753365B CN 110753365 B CN110753365 B CN 110753365B CN 201911042510 A CN201911042510 A CN 201911042510A CN 110753365 B CN110753365 B CN 110753365B
Authority
CN
China
Prior art keywords
base station
user
channel
algorithm
genetic algorithm
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
CN201911042510.4A
Other languages
Chinese (zh)
Other versions
CN110753365A (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.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
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 University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN201911042510.4A priority Critical patent/CN110753365B/en
Publication of CN110753365A publication Critical patent/CN110753365A/en
Application granted granted Critical
Publication of CN110753365B publication Critical patent/CN110753365B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/02Arrangements for optimising operational condition
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W24/00Supervisory, monitoring or testing arrangements
    • H04W24/06Testing, supervising or monitoring using simulated traffic
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W28/00Network traffic management; Network resource management
    • H04W28/02Traffic management, e.g. flow control or congestion control
    • H04W28/0215Traffic management, e.g. flow control or congestion control based on user or device properties, e.g. MTC-capable devices
    • H04W28/0221Traffic management, e.g. flow control or congestion control based on user or device properties, e.g. MTC-capable devices power availability or consumption
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W28/00Network traffic management; Network resource management
    • H04W28/02Traffic management, e.g. flow control or congestion control
    • H04W28/0268Traffic management, e.g. flow control or congestion control using specific QoS parameters for wireless networks, e.g. QoS class identifier [QCI] or guaranteed bit rate [GBR]
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W52/00Power management, e.g. Transmission Power Control [TPC] or power classes
    • H04W52/04Transmission power control [TPC]
    • H04W52/30Transmission power control [TPC] using constraints in the total amount of available transmission power
    • H04W52/36Transmission power control [TPC] using constraints in the total amount of available transmission power with a discrete range or set of values, e.g. step size, ramping or offsets
    • H04W52/367Power values between minimum and maximum limits, e.g. dynamic range

Landscapes

  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Mobile Radio Communication Systems (AREA)

Abstract

本发明公开了一种异构蜂窝网络干扰协调方法,包括:在模拟的真实下行网络链路场景下,以最小化基站总发射功率消耗为目标,分别以用户QoS需求、宏基站最大发射功率、微基站最大发射功率、基站接入限制和信道复用限制为约束条件,确定优化问题;求解优化问题;对求解得到结果解码得到相应的用户接入、信道分配和功率控制策略。本发明的方法通过采用较高还原度和可靠度的网络场景,并基于该场景将用户接入策略、信道分配和功率控制策略协同考虑,建立了用户接入策略、信道分配和功率控制协同优化方案,并提出了具有较低复杂度的启发式算法求解该问题,可以有效保证用户差异化QoS需求以及基站最大发射功率约束的情况下,实现系统总功率消耗最小化。

Figure 201911042510

The invention discloses a heterogeneous cellular network interference coordination method, comprising: in a simulated real downlink network link scenario, aiming at minimizing the total transmission power consumption of the base station, taking the user QoS requirements, the maximum transmission power of the macro base station, The maximum transmit power of micro base station, base station access limit and channel multiplexing limit are used as constraints to determine the optimization problem; solve the optimization problem; decode the results obtained from the solution to obtain the corresponding user access, channel allocation and power control strategies. The method of the present invention establishes the collaborative optimization of user access strategy, channel allocation and power control by adopting a network scenario with high reduction degree and reliability, and considering the user access strategy, channel allocation and power control strategy cooperatively based on the scenario. A heuristic algorithm with lower complexity is proposed to solve the problem, which can effectively ensure the differentiated QoS requirements of users and the constraints of the maximum transmit power of the base station, and minimize the total power consumption of the system.

Figure 201911042510

Description

异构蜂窝网络干扰协调方法Interference coordination method for heterogeneous cellular networks

技术领域technical field

本发明属于移动通信领域,特别涉及一种异构蜂窝网络中的资源调度优化技术。The invention belongs to the field of mobile communication, and in particular relates to a resource scheduling optimization technology in a heterogeneous cellular network.

背景技术Background technique

随着移动用户需求的提升和应用的多样化,移动数据量以指数级增长的速度不断逼近现有蜂窝移动网络的容量上限,从而导致蜂窝网络过载而造成服务质量下降,因而需要进一步提升蜂窝移动网络的容量来满足用户不断增长的需求。一种解决方案是采用异构分层覆盖和增加基站部署,进而衍生出了密集异构网络,异构组网能够有效提升蜂窝小区的通信容量。With the increasing demands of mobile users and the diversification of applications, the amount of mobile data is approaching the upper limit of the capacity of the existing cellular mobile network at an exponential rate, which leads to the overload of the cellular network and the deterioration of service quality. Therefore, it is necessary to further improve the cellular mobile network The capacity of the network to meet the ever-increasing demands of users. One solution is to use heterogeneous layered coverage and increase base station deployment, which in turn derives a dense heterogeneous network. Heterogeneous networking can effectively improve the communication capacity of cells.

虽然密集异构组网一定程度上提高了移动通信网络的容量,但基站密集部署带来了严重的层内干扰,异构网络采用的频谱共享机制也会造成严重的跨层干扰,这导致异构网络频谱效率低下。此外,由于业务负载往往呈现局部热点化的特点,异构网络中多维度的资源分配也是一项巨大的挑战。若所有基站都以最大发射功率工作,不仅会造成能源浪费,能效降低,还会导致基站间干扰过大,信道质量降低,影响用户体验。通过优化用户与基站、信道的关联以及功率控制,在保证用户服务质量(Quality of Service,QoS)的前提下,减少基站间干扰、减小功耗,是本领域研究讨论的主要问题。Although dense heterogeneous networking improves the capacity of mobile communication networks to a certain extent, the dense deployment of base stations brings serious intra-layer interference, and the spectrum sharing mechanism used in heterogeneous networks also causes serious cross-layer interference, which leads to heterogeneous The spectral efficiency of the network is low. In addition, multi-dimensional resource allocation in heterogeneous networks is also a huge challenge due to the characteristics of local hotspots in service loads. If all base stations work at the maximum transmit power, it will not only cause energy waste and reduce energy efficiency, but also cause excessive interference between base stations, reduce channel quality, and affect user experience. By optimizing the association between users and base stations, channels, and power control, under the premise of ensuring user Quality of Service (QoS), reducing interference between base stations and reducing power consumption are the main issues discussed in this field.

目前关于异构蜂窝网络的研究范围十分广泛,包括蜂窝网络通信的相关技术、应用场景、实际应用中可能存在的问题(如接入策略、资源管理、功率控制)等。文献“李俭慧.密集异构网络的业务卸载技术研究[D].北京邮电大学,2017”提出了基于有效容量的延时业务卸载方案,在满足用户服务质量的限定条件下,最大化系统的能量效率;文献“王露.异构蜂窝网络的用户卸载策略研究[D].中国科学技术大学,2017”综合考虑用户归属关系、几乎空白子帧比例等参数构建了一个联合效用优化问题,并采用Gauss-Seidel方法求解。At present, the research scope of heterogeneous cellular networks is very wide, including related technologies of cellular network communication, application scenarios, and possible problems in practical applications (such as access strategy, resource management, power control), etc. The paper "Li Jianhui. Research on service offloading technology for dense heterogeneous networks [D]. Beijing University of Posts and Telecommunications, 2017" proposes a delayed service offloading scheme based on effective capacity, which maximizes the energy of the system under the limited conditions of user service quality. Efficiency; the literature "Wang Lu. Research on User Unloading Strategy for Heterogeneous Cellular Networks [D]. University of Science and Technology of China, 2017" comprehensively considers parameters such as user affiliation and the proportion of almost blank subframes to construct a joint utility optimization problem, and adopts Gauss-Seidel method to solve.

现有研究表明,在密集异构蜂窝网络中引入有效的资源调度策略能够提高系统的频谱利用率,减轻基站的功耗,合理的接入策略和信道分配、有效的干扰协调可以在一定程度上解决密集异构网络面临的一些实际问题;但实际通信场景中需要对不同用户的差异化QoS需求进行有效保障,且不同基站最大发射功率等性能存在差异。上述现有技术均未能很好地兼顾。Existing research shows that the introduction of an effective resource scheduling strategy in a dense heterogeneous cellular network can improve the spectrum utilization of the system and reduce the power consumption of the base station. It solves some practical problems faced by dense heterogeneous networks; however, in actual communication scenarios, it is necessary to effectively guarantee the differentiated QoS requirements of different users, and there are differences in performance such as the maximum transmit power of different base stations. None of the above-mentioned prior art can take into account well.

发明内容SUMMARY OF THE INVENTION

针对现有技术存在的上述问题,本发明提出了一种异构蜂窝网络干扰协调方法,具体包括如下步骤:Aiming at the above problems existing in the prior art, the present invention proposes a heterogeneous cellular network interference coordination method, which specifically includes the following steps:

步骤S1:在模拟的真实下行网络链路场景下,以最小化基站总发射功率消耗为目标,分别以用户QoS需求、宏基站最大发射功率、微基站最大发射功率、基站接入限制和信道复用限制为约束条件,确定优化问题;Step S1: In the simulated real downlink network scenario, with the goal of minimizing the total transmit power consumption of the base station, the user QoS requirements, the maximum transmit power of the macro base station, the maximum transmit power of the micro base station, the base station access limit and the channel complexity are respectively taken as the goal. Determine the optimization problem with constraints as constraints;

步骤S2:求解步骤S1确定优化问题;Step S2: Solve the optimization problem determined in step S1;

步骤S3:对步骤S2得到结果解码得到相应的用户接入、信道分配和功率控制策略。Step S3: Decode the result obtained in Step S2 to obtain a corresponding user access, channel allocation and power control strategy.

进一步地,所述步骤S1模拟的真实链路场景包括:M个移动用户、N个基站(其中L个宏基站和N-L个微基站),N个基站均匀部署构成小区,M个移动用户随机分布在小区内;所有基站同频部署,且频带被划分为I条正交信道,每个移动用户能且仅能连接一个基站,并从一条或多条信道获得服务,其他基站在相同基站发射的信号会对该用户产生干扰,不同信道间不存在干扰。Further, the real link scenario simulated in step S1 includes: M mobile users, N base stations (including L macro base stations and N-L micro base stations), N base stations are evenly deployed to form a cell, and M mobile users are randomly distributed. In a cell; all base stations are deployed on the same frequency, and the frequency band is divided into I orthogonal channels, each mobile user can connect to only one base station, and obtain services from one or more channels, and other base stations transmit at the same base station. The signal will interfere with the user, and there is no interference between different channels.

进一步地,所述步骤S1确定优化问题的为混合整数非线性优化问题。Further, the step S1 determines that the optimization problem is a mixed integer nonlinear optimization problem.

更进一步地,所述步骤S2具体利用基于遗传算法和粒子群算法的混合启发式搜索算法进行求解。Further, the step S2 is specifically solved by using a hybrid heuristic search algorithm based on a genetic algorithm and a particle swarm algorithm.

更进一步地,求解的具体步骤如下:Further, the specific steps of solving are as follows:

S21:使用相同格式编码遗传算法中基因序列和粒子群算法中位置信息,初始化遗传算法;S21: Use the same format to encode the gene sequence in the genetic algorithm and the position information in the particle swarm algorithm, and initialize the genetic algorithm;

S22:利用有限次遗传算法迭代求得一组可行解;S22: Obtain a set of feasible solutions by iteratively using a finite number of genetic algorithms;

S23:初始化粒子群算法,将遗传算法得到的结果输入粒子群优化算法求解;S23: Initialize the particle swarm optimization algorithm, and input the results obtained by the genetic algorithm into the particle swarm optimization algorithm to solve;

S24:将粒子群优化算法得到的解重新输入遗传算法,直至结果收敛或达到迭代次数上限。S24: Re-input the solution obtained by the particle swarm optimization algorithm into the genetic algorithm until the result converges or the upper limit of the number of iterations is reached.

更进一步地,步骤S21中所用遗传算法中染色体序列和粒子群优化算法中的位置信息为与用户接入策略和信道分配策略相对应的0-1整数型矩阵,以及和基站功率控制相关的实数矩阵,个体和粒子的适应度表达式即为问题的优化目标基站总功耗。Further, the chromosome sequence in the genetic algorithm used in step S21 and the position information in the particle swarm optimization algorithm are the 0-1 integer matrix corresponding to the user access strategy and the channel allocation strategy, and the real number related to the power control of the base station. The fitness expression of matrix, individual and particle is the total power consumption of the optimization target base station of the problem.

进一步地,在一次大循环中,对于遗传算法每次迭代依次对父代种群进行交叉操作、变异操作生成子代种群后,记录最优策略和最优适应度,然后执行选择操作进入下一次迭代,在有限次迭代后,强制停止遗传算法,输出当前种群。Further, in a large cycle, for each iteration of the genetic algorithm, the parent population is subjected to the crossover operation and mutation operation to generate the child population, and the optimal strategy and optimal fitness are recorded, and then the selection operation is performed to enter the next iteration. , after a finite number of iterations, the genetic algorithm is forced to stop and the current population is output.

进一步地,在每一次启动粒子群优化算法时,将所有粒子速度重置为0,并输入遗传算法得到当前种群的基因信息作为位置信息,迭代一次后,计算粒子速度并执行迭代。Further, every time the particle swarm optimization algorithm is started, the velocities of all particles are reset to 0, and the genetic information of the current population is obtained by inputting the genetic algorithm as the position information. After one iteration, the particle velocities are calculated and the iteration is performed.

本发明的有益效果:本发明的异构蜂窝网络干扰协调方法通过采用较高还原度和可靠度的网络场景,并基于该场景将用户接入策略、信道分配和功率控制策略协同考虑,建立了用户接入策略、信道分配和功率控制协同优化方案,并提出了具有较低复杂度的启发式算法求解该问题。本发明的方法针对异构蜂窝网络中的下行传输场景,可以在有效保证用户差异化QoS需求以及基站最大发射功率约束的情况下,实现系统总功率消耗最小化。Beneficial effects of the present invention: The heterogeneous cellular network interference coordination method of the present invention establishes a network scenario with a relatively high degree of restoration and reliability, and based on the scenario, the user access strategy, channel allocation and power control strategy are considered collaboratively. User access strategy, channel allocation and power control co-optimization scheme, and a heuristic algorithm with lower complexity is proposed to solve the problem. Aiming at the downlink transmission scenario in the heterogeneous cellular network, the method of the present invention can minimize the total power consumption of the system under the condition of effectively ensuring the differentiated QoS requirements of users and the constraint of the maximum transmit power of the base station.

附图说明Description of drawings

图1为本发明实施例的系统场景示意图。FIG. 1 is a schematic diagram of a system scenario according to an embodiment of the present invention.

图2为本发明实施例的方法的流程示意图。FIG. 2 is a schematic flowchart of a method according to an embodiment of the present invention.

具体实施方式Detailed ways

为便于本领域技术人员理解本发明的技术内容,下面结合附图对本发明内容进一步阐释。In order to facilitate those skilled in the art to understand the technical content of the present invention, the content of the present invention will be further explained below with reference to the accompanying drawings.

如图1所示为本发明的实施例所构建的网络场景,考虑下行链路,基站同频部署的场景。其中,MBS表示宏基站,SBS表示微基站,MU表示移动用户,AN×M和BM×I分别表示用户与基站的关联情况和用户占用信道的情况;虚线表示干扰信号,实线为有用信号。FIG. 1 shows a network scenario constructed by an embodiment of the present invention, considering the downlink and the scenario of co-frequency deployment of base stations. Among them, MBS represents the macro base station, SBS represents the micro base station, MU represents the mobile user, A N×M and BM×I respectively represent the association between the user and the base station and the channel occupied by the user; the dotted line represents the interference signal, and the solid line is useful. Signal.

本实施例中模拟真实下行链路通信场景,具体为:N个基站和M个用户组成的异构蜂窝网络下行链路系统,包含L个宏基站和N-L个微基站,另有M个移动用户随机分布在区域中。In this embodiment, a real downlink communication scenario is simulated, specifically: a heterogeneous cellular network downlink system composed of N base stations and M users, including L macro base stations and N-L micro base stations, and M mobile users randomly distributed in the area.

令集合

Figure BDA0002253248320000031
为基站集合,
Figure BDA0002253248320000032
为用户集合,
Figure BDA0002253248320000033
表示第n个基站,n<L时表示宏基站,反之则表示微基站,所有基站同频部署,存在同层干扰和跨层干扰;
Figure BDA0002253248320000034
表示第m个用户设备,假设每个用户只能关联一个基站,令an,m∈{0,1}为用户与基站的关联变量,an,m=1表示用户m接入了基站n,反之an,m=0。set of orders
Figure BDA0002253248320000031
is the set of base stations,
Figure BDA0002253248320000032
collection of users,
Figure BDA0002253248320000033
Indicates the nth base station. When n<L, it means a macro base station. Otherwise, it means a micro base station. All base stations are deployed on the same frequency, and there is co-layer interference and cross-layer interference;
Figure BDA0002253248320000034
Represents the mth user equipment, assuming that each user can only be associated with one base station, let an ,m ∈{0,1} be the association variable between the user and the base station, an n,m =1 indicates that user m accesses base station n , otherwise an ,m =0.

多个接入同一基站的用户会共享基站的系统资源,如以时分复用或频分复用获得服务,本实施例中将资源块作为基站系统资源的基本单位,对应于每个子信道,令

Figure BDA0002253248320000035
为基站资源块的集合,即每个基站共有I个子信道。bm,i表示用户m对第i个信道的占用情况,bm,i=1表示用户m占用了第i个子信道,并且会受到其他基站在这个信道上发射信号的干扰。令PS表示微基站最大发射功率,PM表示宏基站的最大发射功率。另外网络中存在集中控制器,可以实时感知基站和用户状态,并做出相应决策。Multiple users accessing the same base station will share the system resources of the base station, such as obtaining services by time division multiplexing or frequency division multiplexing.
Figure BDA0002253248320000035
is a set of base station resource blocks, that is, each base station has a total of I subchannels. b m,i indicates the occupancy of the i-th channel by user m, and b m,i =1 indicates that the user m occupies the i-th sub-channel and will be interfered by signals transmitted by other base stations on this channel. Let PS denote the maximum transmit power of the micro base station, and PM denote the maximum transmit power of the macro base station. In addition, there is a centralized controller in the network, which can sense the status of base stations and users in real time, and make corresponding decisions.

对于网络每一个调度周期,假设控制器根据当前用户位置信息,以及信道状态信息,实时计算出用户接入、信道分配和功率控制策略。相应策略服从的约束如下:For each scheduling period of the network, it is assumed that the controller calculates user access, channel allocation and power control strategies in real time according to the current user location information and channel state information. The corresponding policy obeys the following constraints:

用户接入策略:每个移动用户仅允许与一个基站建立连接,基站服务用户不设上限(但不能超过可分配子信道的数量)。User access strategy: each mobile user is only allowed to establish a connection with one base station, and the base station serves no upper limit (but cannot exceed the number of sub-channels that can be allocated).

信道分配策略:所有基站同频部署,频域被划分为多个正交子信道。基站可以为用户提供一个或多个子信道并保证每个用户至少占用一个子信道。Channel allocation strategy: All base stations are deployed on the same frequency, and the frequency domain is divided into multiple orthogonal sub-channels. The base station can provide users with one or more sub-channels and ensure that each user occupies at least one sub-channel.

功率控制策略:基站可将功率任意分配在每个子信道上,仅需满足总功率不超过该基站最大发射功率约束。Power control strategy: The base station can arbitrarily allocate power to each sub-channel, and only needs to satisfy the constraint that the total power does not exceed the maximum transmit power of the base station.

基于该网络场景,本发明的一种异构蜂窝网络干扰协调方法,具体包括以下步骤:Based on the network scenario, a heterogeneous cellular network interference coordination method of the present invention specifically includes the following steps:

S1、在模拟的真实下行网络场景下,以最小化基站总发射功率消耗为目标,分别以用户QoS需求、宏基站最大发射功率、微基站最大发射功率、基站接入限制和信道复用限制为约束条件,确定优化问题。具体为:以最小化系统总功率消耗为目标函数,分别以移动用户的QoS需求vm、宏基站最大发射功率pn1,i、微基站最大发射功率pn2,i、用户接入指示符an,m以及信道分配指示符bm,i的取值为约束条件,得到第一优化问题。S1. In the simulated real downlink network scenario, aiming at minimizing the total transmit power consumption of the base station, the user QoS requirements, the maximum transmit power of the macro base station, the maximum transmit power of the micro base station, the base station access limit and the channel reuse limit are respectively Constraints to determine the optimization problem. Specifically, the objective function is to minimize the total power consumption of the system, and the QoS requirement vm of the mobile user, the maximum transmit power p n1,i of the macro base station, the maximum transmit power p n2,i of the micro base station, and the user access indicator a are respectively taken as the objective function. The values of n,m and the channel assignment indicator b m,i are constraints, and the first optimization problem is obtained.

优化问题表达式如下:The optimization problem expression is as follows:

P1:

Figure BDA0002253248320000041
P1:
Figure BDA0002253248320000041

s.t.s.t.

C1:an,m∈{0,1}C1: a n,m ∈ {0,1}

C2:

Figure BDA0002253248320000042
C2:
Figure BDA0002253248320000042

C3:bm,i∈{0,1}C3: b m,i ∈ {0,1}

C4:

Figure BDA0002253248320000043
C4:
Figure BDA0002253248320000043

C5:

Figure BDA0002253248320000044
C5:
Figure BDA0002253248320000044

C6:

Figure BDA0002253248320000045
C6:
Figure BDA0002253248320000045

C7:

Figure BDA0002253248320000046
C7:
Figure BDA0002253248320000046

其中,P1表示优化问题的目标函数;约束条件C1中,an,m表示用户与基站的关联变量,an,m=1表示用户m接入了基站n,an,m=0则表示用户与建站没有建立连接;约束条件C2保证了用户至少且至多接入一个基站;约束条件C3中,bm,i表示用户m对第i个子信道的占用情况,bm,i=1表示用户m占用了第i个子信道,并且会受到其他基站在这个信道上上发射信号的干扰,bm,i=0表示用户与信道m与信道i无关联;约束条件C4保证了对于特定基站,它的每个子信道只能被一个移动用户占用(接入不同基站的用户可以复用同一信道);约束条件C5、C6中,

Figure BDA0002253248320000047
对应了宏基站n1和微基站n2在信道i上的发射功率,两约束分别保证了宏基站与微基站的总发射功率不会超过各自的发射功率上限PM/PS;约束条件C7保证了用户m得到的吞吐量Cm满足其QoS需求vm(最小数据速率)。最大发射功率PM=50w,PS=20w,用户差异化QoS需求vm随机分布于[105,3×105]bit/s。Among them, P1 represents the objective function of the optimization problem; in the constraint condition C1, an ,m represents the associated variable between the user and the base station, an n,m =1 means that the user m accesses the base station n, an n,m =0 means that The user does not establish a connection with the station establishment; Constraint C2 ensures that the user accesses at least one base station at most; In Constraint C3, b m,i represents the occupancy of the i-th subchannel by user m, and b m,i =1 represents the user m occupies the i-th subchannel, and will be interfered by signals transmitted by other base stations on this channel. b m,i = 0 means that the user is not associated with channel m and channel i; constraint C4 ensures that for a specific base station, it Each sub-channel can only be occupied by one mobile user (users accessing different base stations can reuse the same channel); in the constraints C5 and C6,
Figure BDA0002253248320000047
Corresponding to the transmit power of macro base station n 1 and micro base station n 2 on channel i, the two constraints respectively ensure that the total transmit power of macro base station and micro base station will not exceed their respective transmit power upper limit P M /PS ; Constraint C7 It is guaranteed that the throughput C m obtained by user m satisfies its QoS requirement vm (minimum data rate). The maximum transmit power P M =50w, P S =20w, and the user's differentiated QoS requirement vm is randomly distributed in [10 5 , 3×10 5 ]bit/s.

移动用户的QoS即用户从关联基站得到的数据速率,计算式为:The QoS of a mobile user is the data rate obtained by the user from the associated base station, and the calculation formula is:

Figure BDA0002253248320000051
Figure BDA0002253248320000051

其中,Cn,m,i为基站n在第i个子信道上为用户m提供的数据速率,表达式如下:Among them, C n,m,i is the data rate provided by base station n for user m on the ith subchannel, and the expression is as follows:

Figure BDA0002253248320000052
Figure BDA0002253248320000052

其中,B为每个正交子信道的带宽,N0为背景噪声的功率谱密度,Hn,m为信道增益模型:where B is the bandwidth of each orthogonal sub-channel, N 0 is the power spectral density of the background noise, and H n,m is the channel gain model:

Hn,m=PLn,m×HR×HS H n,m =PL n,m × H R ×HS

HR为由于多径效应产生的瑞利衰落,HS为随机的阴影衰落,具体地,PLn,m表达式为: HR is Rayleigh fading due to multipath effect, and H S is random shadow fading. Specifically, PL n,m is expressed as:

Figure BDA0002253248320000053
Figure BDA0002253248320000053

其中,dn,m表示用户同基站之间的距离。Among them, d n,m represents the distance between the user and the base station.

In,m,i为用户m在基站n的第i个子信道上受到的干扰,表达式如下:I n,m,i is the interference received by user m on the i-th subchannel of base station n, and the expression is as follows:

Figure BDA0002253248320000054
Figure BDA0002253248320000054

S2、利用遗传算法和粒子群算法联合求解上述混合整数非线性规划问题,S2. Use genetic algorithm and particle swarm algorithm to jointly solve the above mixed integer nonlinear programming problem,

步骤S1确定的问题是NP-hard问题,难以直接求解,故采用具有低复杂度的混合启发式搜索算法求解该问题。在算法的大循环内,反复运行遗传算法和粒子群算法,以达到良好的全局搜索性并加快收敛速度。具体流程如图2,具体步骤如下:The problem determined in step S1 is an NP-hard problem, which is difficult to solve directly, so a hybrid heuristic search algorithm with low complexity is used to solve the problem. In the large loop of the algorithm, the genetic algorithm and the particle swarm algorithm are repeatedly run to achieve a good global search and speed up the convergence. The specific process is shown in Figure 2, and the specific steps are as follows:

S21、使用相同格式编码遗传算法中基因序列和粒子群算法中位置信息以及适应度,初始化算法;具体实现如下S21. Use the same format to encode the gene sequence in the genetic algorithm and the position information and fitness in the particle swarm algorithm, and initialize the algorithm; the specific implementation is as follows

对于遗传算法,种群中每个个体代表一组策略,个体的染色体

Figure BDA0002253248320000055
Figure BDA0002253248320000056
为0-1整数型矩阵,
Figure BDA0002253248320000057
为实数矩阵,分别对应用户与基站的关联情况(用户接入第几个基站)、信道分配策略(用户占用接入基站的哪几个信道)、功率分配策略(基站在每个子信道上分配的功率)。For genetic algorithms, each individual in the population represents a set of strategies, and the chromosomes of the individual
Figure BDA0002253248320000055
Figure BDA0002253248320000056
is a 0-1 integer matrix,
Figure BDA0002253248320000057
It is a real number matrix, which corresponds to the association between the user and the base station (the number of base stations the user accesses), the channel allocation strategy (which channels the user occupies to access the base station), and the power allocation strategy (the number of channels allocated by the base station on each subchannel). power).

个体的适应度代表一组策略的性能好坏,在本问题中,以优化目标总功耗来衡量策略的性能表现:The fitness of an individual represents the performance of a set of strategies. In this problem, the total power consumption of the optimization target is used to measure the performance of the strategy:

Figure BDA0002253248320000061
Figure BDA0002253248320000061

粒子群优化算法中每个粒子代表一组策略,具有位置、速度和适应度三个属性。粒子的位置和适应度分别对应遗传算法个体的染色体和适应度,表示当前策略和性能指标;粒子的速度是由算法计算得出的预期获得更好适应度的策略变化方向。In the particle swarm optimization algorithm, each particle represents a set of strategies, and has three attributes: position, speed and fitness. The position and fitness of the particle correspond to the chromosome and fitness of the individual genetic algorithm, respectively, and represent the current strategy and performance index; the speed of the particle is calculated by the algorithm and is the direction of strategy change that is expected to obtain better fitness.

初始化种群数量和粒子群粒子数量NP、循环迭代次数T、遗传算法阶段迭代次数T1和粒子群优化阶段最大迭代次数T2,并设置算法参数:遗传算法中的子代种群数量NP1、交叉概率Pc、变异概率Pm,粒子群优化算法中的加速度ω、当前集体质心权值c1、历史最佳质心权值c2Initialize the number of populations and particle swarm particles NP, the number of loop iterations T, the number of iterations in the genetic algorithm stage T 1 and the maximum number of iterations in the particle swarm optimization stage T 2 , and set the algorithm parameters: the number of offspring populations in the genetic algorithm NP 1 , the crossover The probability P c , the mutation probability P m , the acceleration ω in the particle swarm optimization algorithm, the current collective centroid weight c 1 , and the historical best centroid weight c 2 .

S22、利用有限次遗传算法求得一组可行解;S22, obtain a set of feasible solutions by using the finite-order genetic algorithm;

每次遗传算法迭代的操作步骤为交叉操作、变异操作、选择操作。The operation steps of each genetic algorithm iteration are crossover operation, mutation operation and selection operation.

交叉操作:对于每一代种群,随机选择两个父代,以Pc的概率交换一部分染色体片段,并检测得到的新策略是否满足约束条件,若满足条件则得到一个子代个体。Crossover operation: For each generation of the population, randomly select two parents, exchange a part of the chromosome segments with the probability of P c , and check whether the obtained new strategy satisfies the constraints, and if the conditions are met, a child individual is obtained.

变异操作:随机选择一个个体,以概率Pm对其执行变异操作,抽取一个染色体片段在可行域内随机变化。Mutation operation: randomly select an individual, perform mutation operation on it with probability P m , and extract a chromosome segment to randomly change in the feasible region.

变异操作的规则(即变异算子)为:The rules of mutation operation (ie mutation operator) are:

Figure BDA0002253248320000062
Figure BDA0002253248320000062

Figure BDA0002253248320000063
Figure BDA0002253248320000063

其中,

Figure BDA0002253248320000064
代表抽取到的片段为整数型变量,
Figure BDA0002253248320000065
代表抽取到的片段为实数型变量,γ1、γ2均为[0,1]区间上的随机数(根据
Figure BDA0002253248320000066
规模大小,γ1、γ2可能为随机数列或矩阵),round计算符表示取与括号内数值最接近的整数,xmax、xmin、ymax、ymin分别表示抽取染色体片段每个分量可以取到的最大值或最小值。in,
Figure BDA0002253248320000064
Represents the extracted segment as an integer variable,
Figure BDA0002253248320000065
Represents that the extracted segments are real-number variables, and γ 1 and γ 2 are random numbers in the [0,1] interval (according to
Figure BDA0002253248320000066
Scale, γ 1 , γ 2 may be random number columns or matrices), the round operator means to take the integer closest to the value in the brackets, x max , x min , y max , y min respectively indicate that each component of the extracted chromosome segment can be The maximum or minimum value obtained.

循环以上步骤,共生成NP1个子代个体后,计算每个个体的适应度,并以适应度作为权值,使用轮盘赌算法选择NP个个体作为下一代种群。特别地,若当前个体中适应度最佳的个体被淘汰,则强制令该个体替换掉下一代种群中适应度最低的个体,保证算法朝更优方向进行。Repeat the above steps to generate NP 1 offspring individuals, calculate the fitness of each individual, and use the fitness as the weight, and use the roulette algorithm to select NP individuals as the next generation population. In particular, if the individual with the best fitness in the current individual is eliminated, the individual is forced to replace the individual with the lowest fitness in the next generation of the population to ensure that the algorithm proceeds in a better direction.

迭代T1次后,进入下一阶段。After iterating T 1 times, enter the next stage.

S23、初始化粒子群算法,将S22得到的结果输入粒子群优化算法求解;S23, initialize the particle swarm algorithm, and input the result obtained in S22 into the particle swarm optimization algorithm to solve;

粒子群算法中,每个粒子的位置代表一种方案,速度决定了方案迭代演进的速度,为保持算法一致性,位置的表达同遗传算法中染色体的表达,定义速度为每个粒子两次迭代位置的变化量,即:In particle swarm optimization, the position of each particle represents a plan, and the speed determines the speed of the iterative evolution of the plan. In order to maintain the consistency of the algorithm, the expression of the position is the same as the expression of the chromosome in the genetic algorithm, and the speed is defined as two iterations for each particle. The amount of change in position, that is:

Figure BDA0002253248320000071
Figure BDA0002253248320000071

Figure BDA0002253248320000072
Figure BDA0002253248320000072

pbestr,d为当前最佳个体的位置,gbestr,d为历史最优适应度出现的位置,

Figure BDA0002253248320000073
表示当前粒子的位置。pbest r,d is the position of the current best individual, gbest r,d is the position where the historical best fitness appears,
Figure BDA0002253248320000073
Indicates the current particle position.

上式中,

Figure BDA0002253248320000074
表示S23要求每次粒子群算法启动前将速度重置为0,输入遗传算法的解并进行一次迭代之后,重新加入速度分量进行优化。In the above formula,
Figure BDA0002253248320000074
Indicates that S23 requires that the velocity be reset to 0 before each particle swarm optimization algorithm is started, and after inputting the solution of the genetic algorithm and performing one iteration, the velocity component is re-added for optimization.

其中,

Figure BDA0002253248320000075
表示粒子的惯性,即算法会尝试保持上一次策略迭代的变化方向,以期获得更好的结果,
Figure BDA0002253248320000076
表示粒子会向当前集体质心靠拢,
Figure BDA0002253248320000077
表示粒子会向历史最佳的位置靠拢。in,
Figure BDA0002253248320000075
represents the inertia of the particle, i.e. the algorithm will try to maintain the direction of change from the last policy iteration in the hope of getting better results,
Figure BDA0002253248320000076
means that the particles will move closer to the current collective center of mass,
Figure BDA0002253248320000077
Indicates that the particle will move closer to the historically optimal position.

每次迭代更新所有粒子的位置和速度,其中整数型变量和实数型变量迭代表达式如下:The positions and velocities of all particles are updated in each iteration, and the iterative expressions of integer variables and real variables are as follows:

Figure BDA0002253248320000078
Figure BDA0002253248320000078

Figure BDA0002253248320000079
Figure BDA0002253248320000079

其中,x代表整数型变量,y代表连续变量。Among them, x represents an integer variable, and y represents a continuous variable.

计算所有粒子的适应度,并记当前适应度最高的粒子的位置为当前集体质心,如果该适应度优于历史最佳,则记当前位置为历史最佳质心,当前适应度为历史最佳适应度。Calculate the fitness of all particles, and record the position of the particle with the highest fitness as the current collective centroid. If the fitness is better than the best in history, record the current position as the best in history, and the current fitness as the best in history Spend.

粒子群优化迭代T2次后,进入下一阶段。After the particle swarm optimization iterates T 2 times, it enters the next stage.

S24、将S23得到的解重新输入遗传算法,直至结果收敛或达到迭代次数上限T;S24, re-input the solution obtained in S23 into the genetic algorithm, until the result converges or reaches the upper limit T of the number of iterations;

需要说明的是:由于遗传算法具有较好的全局搜索性,但收敛速度缓慢。粒子群优化算法收敛速度更快,但极易过早收敛到局部最优。使用遗传算法进行粗粒度的搜索,并将结果输入粒子群算法加速算法收敛。当算法外层循环达到次数上限或数代最佳个体适应度不再降低时,算法结束,输出最佳个体的基因信息和适应度。It should be noted that the genetic algorithm has good global search performance, but the convergence speed is slow. The particle swarm optimization algorithm converges faster, but it is easy to prematurely converge to a local optimum. Coarse-grained searches are performed using genetic algorithms, and the results are fed into particle swarm optimization to accelerate algorithm convergence. When the outer loop of the algorithm reaches the upper limit or the fitness of the best individual for several generations is no longer reduced, the algorithm ends, and the genetic information and fitness of the best individual are output.

S3、对最终得到历史最佳集体质心解码得到相应的用户接入、信道分配和功率控制策略,历史最佳适应度即为系统最小总功耗。S3. Decode the finally obtained historical best collective centroid to obtain a corresponding user access, channel allocation and power control strategy, and the historical best fitness is the minimum total power consumption of the system.

本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。Those of ordinary skill in the art will appreciate that the embodiments described herein are intended to assist readers in understanding the principles of the present invention, and it should be understood that the scope of protection of the present invention is not limited to such specific statements and embodiments. Various modifications and variations of the present invention are possible for those skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention shall be included within the scope of the claims of the present invention.

Claims (4)

1. A heterogeneous cellular network interference coordination method specifically comprises the following steps:
step S1: under a simulated real downlink network link scene, determining an optimization problem by respectively taking a user QoS requirement, a macro base station maximum transmitting power, a micro base station maximum transmitting power, a base station access limit and a channel multiplexing limit as constraint conditions with the aim of minimizing the total transmitting power consumption of a base station;
the real link scenario simulated in step S1 includes: the mobile communication system comprises M mobile users and N base stations, wherein the N base stations comprise L macro base stations and N-L micro base stations, the N base stations are uniformly deployed to form a cell, and the M mobile users are randomly distributed in the cell; all base stations are deployed at the same frequency, the frequency band is divided into I orthogonal channels, each mobile user can be connected with only one base station and obtains service from one or more channels, signals transmitted by other base stations in the same channel can generate interference on the user, and interference does not exist between different channels;
the step S1 determines that the optimization problem is specifically a mixed integer nonlinear optimization problem;
the mixed integer nonlinear optimization problem specifically comprises the following steps: minimizing the total power consumption of the system as an objective function, respectively based on the QoS requirement v of the mobile usermMacro base station maximum transmit power
Figure FDA0002904689030000011
Maximum transmitting power of micro base station
Figure FDA0002904689030000012
User access indicator an,mAnd a channel allocation indicator bm,iObtaining a first optimization problem by taking the value of (a) as a constraint condition;
the optimization problem expression is as follows:
Figure FDA0002904689030000013
s.t.
C1:an,m∈{0,1}
Figure FDA0002904689030000014
C3:bm,i∈{0,1}
Figure FDA0002904689030000015
Figure FDA0002904689030000016
Figure FDA0002904689030000017
Figure FDA0002904689030000018
wherein P1 represents the objective function of the optimization problem;
in constraint C1, an,mRepresenting the association variable of the user with the base station, an,m1 indicates that user m accesses base station n, an,mIf 0, it means that the user has not established connection with the base station;
constraint C2 ensures that a user has at least and at most one access to a base station;
in constraint C3, bm,iRepresenting the occupation of the ith sub-channel by user m, bm,i1 indicates that user m occupies the ith sub-channel and is interfered by signals transmitted by other base stations on the channel, bm,i0 means that user m is not associated with channel i;
constraint C4 ensures that each sub-channel of a particular base station can only be occupied by one mobile subscriber, and that subscribers accessing different base stations reuse the same channel;
in the constraints C5 and C6,
Figure FDA0002904689030000021
corresponds to a macro base station n1And a micro base station n2On the channel i, the two constraints respectively ensure that the total transmission power of the macro base station and the micro base station does not exceed the respective upper limit P of the transmission powerM/PS
Constraint C7 ensures throughput C for user mmMinimum data rate v to meet its QoS requirementsmMaximum transmission power PM=50w,PSUser differentiated QoS requirements v 20wmIs randomly distributed in [10 ]5,3×105]bit/s;
The QoS of a mobile user, i.e. the data rate the user gets from the associated base station, is calculated as:
Figure FDA0002904689030000022
wherein, Cn,m,iFor the data rate provided by base station n for user m on the ith subchannel, the expression is as follows:
Figure FDA0002904689030000023
where B is the bandwidth of each orthogonal sub-channel, N0Power spectral density, H, of background noisen,mFor the channel gain model:
Hn,m=PLn,m×HR×HS
HRis Rayleigh fading, H, due to multipath effectsSFor random shadow fading, in particular PLn,mThe expression is as follows:
Figure FDA0002904689030000024
wherein d isn,mRepresents the distance between the user and the base station;
In,m,ifor the interference experienced by user m on the ith subchannel of base station n, the expression is as follows:
Figure FDA0002904689030000025
step S2: solving the optimization problem determined in the step S1;
specifically, a hybrid heuristic search algorithm based on a genetic algorithm and a particle swarm algorithm is used for solving, and the solving specifically comprises the following steps:
s21: using the same format to encode gene sequences in the genetic algorithm and position information in the particle swarm algorithm, and initializing the genetic algorithm;
s22: iteratively solving a group of feasible solutions by utilizing a finite genetic algorithm;
s23: initializing a particle swarm algorithm, and inputting a result obtained by the genetic algorithm into the particle swarm optimization algorithm for solving;
s24: inputting the solution obtained by the particle swarm optimization algorithm into the genetic algorithm again until the result is converged or the upper limit of the iteration times is reached;
step S3: and decoding the result obtained in the step S2 to obtain the corresponding user access, channel allocation and power control strategy.
2. The method for coordinating interference in heterogeneous cellular networks according to claim 1, wherein the chromosome sequences in the genetic algorithm and the location information in the particle swarm optimization algorithm used in step S21 are 0-1 integer matrices corresponding to the user access policy and the channel allocation policy, and real matrices related to power control of the base station, and the fitness expression of the individuals and the particles is the total power consumption of the optimization target base station of the problem.
3. The method according to claim 1, wherein in a major loop, after performing crossover operation and mutation operation on a parent population to generate a child population for each iteration of the genetic algorithm, recording an optimal strategy and optimal fitness, performing selection operation to enter the next iteration, after a limited number of iterations, forcibly stopping the genetic algorithm, and outputting the current population.
4. The method according to claim 1, wherein every time the particle swarm optimization algorithm is started, all particle speeds are reset to 0, a genetic algorithm is input to obtain gene information of a current population as position information, and after one iteration, the particle speeds are calculated and the iteration is executed.
CN201911042510.4A 2019-10-30 2019-10-30 Interference coordination method for heterogeneous cellular networks Expired - Fee Related CN110753365B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911042510.4A CN110753365B (en) 2019-10-30 2019-10-30 Interference coordination method for heterogeneous cellular networks

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911042510.4A CN110753365B (en) 2019-10-30 2019-10-30 Interference coordination method for heterogeneous cellular networks

Publications (2)

Publication Number Publication Date
CN110753365A CN110753365A (en) 2020-02-04
CN110753365B true CN110753365B (en) 2021-05-04

Family

ID=69281141

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911042510.4A Expired - Fee Related CN110753365B (en) 2019-10-30 2019-10-30 Interference coordination method for heterogeneous cellular networks

Country Status (1)

Country Link
CN (1) CN110753365B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112422227B (en) * 2020-11-25 2022-05-20 湘潭大学 Method and device for scheduling satellite trusted navigation interference resources
CN113923693B (en) * 2021-11-03 2023-09-22 重庆理工大学 A method and device for optimizing energy efficiency of heterogeneous networks with simulated hybrid power supply
CN119233398A (en) * 2023-06-28 2024-12-31 中兴通讯股份有限公司 Base station power adjustment method, electronic device and storage medium
CN119031393B (en) * 2024-10-29 2025-02-14 华东交通大学 Safe computing unloading method for cache-assisted ultra-dense heterogeneous MEC network

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105392161A (en) * 2015-10-15 2016-03-09 北京航空航天大学 User accessing and power control combined optimization method in cooperative heterogeneous network
CN109219122A (en) * 2017-07-07 2019-01-15 大唐移动通信设备有限公司 A kind of setting method and device of transmission power
CN109831819A (en) * 2019-03-06 2019-05-31 重庆邮电大学 One kind being based on isomery cellular network sub-clustering SMDP base station dormancy method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105392161A (en) * 2015-10-15 2016-03-09 北京航空航天大学 User accessing and power control combined optimization method in cooperative heterogeneous network
CN109219122A (en) * 2017-07-07 2019-01-15 大唐移动通信设备有限公司 A kind of setting method and device of transmission power
CN109831819A (en) * 2019-03-06 2019-05-31 重庆邮电大学 One kind being based on isomery cellular network sub-clustering SMDP base station dormancy method

Also Published As

Publication number Publication date
CN110753365A (en) 2020-02-04

Similar Documents

Publication Publication Date Title
CN110753365B (en) Interference coordination method for heterogeneous cellular networks
CN109905918B (en) NOMA cellular Internet of vehicles dynamic resource scheduling method based on energy efficiency
CN112601284B (en) Downlink multi-cell OFDMA resource allocation method based on multi-agent deep reinforcement learning
Wang et al. Joint resource allocation and power control for D2D communication with deep reinforcement learning in MCC
CN110267338A (en) A joint resource allocation and power control method in D2D communication
CN104955077B (en) A kind of heterogeneous network cell cluster-dividing method and device based on user experience speed
CN106060876B (en) A kind of method of heterogeneous wireless network equally loaded
CN106792451B (en) A D2D Communication Resource Optimization Method Based on Multi-Swarm Genetic Algorithm
CN111787543B (en) 5G communication system resource allocation method based on improved wolf optimization algorithm
CN108965009B (en) Load known user association method based on potential game
Cheng et al. Joint user association and resource allocation in HetNets based on user mobility prediction
Chen et al. A deep learning based resource allocation scheme in vehicular communication systems
CN108366427A (en) Throughput of system and efficiency balance method based on power control in D2D communications
CN106231610A (en) Resource allocation methods based on sub-clustering in Femtocell double-layer network
CN105873214A (en) Resource allocation method of D2D communication system based on genetic algorithm
CN116347635A (en) A NB-IoT wireless resource allocation method based on NOMA and multi-agent reinforcement learning
Yu et al. Interference coordination strategy based on Nash bargaining for small‐cell networks
CN105578482B (en) A kind of cellular heterogeneous network resource allocation method
Zhang et al. Dynamic user-centric clustering for uplink cooperation in multi-cell wireless networks
Li et al. An energy-efficient cell planning strategy for heterogeneous network based on realistic traffic data
CN115866787A (en) Network resource allocation method integrating terminal direct transmission communication and multi-access edge calculation
CN103338453B (en) A kind of dynamic spectrum access method for hierarchical wireless network network and system
Jiang et al. Dueling double deep q-network based computation offloading and resource allocation scheme for internet of vehicles
Hou et al. Caching and resource allocation in small cell networks
CN103327541A (en) Service unloading method based on different QoS

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: 20210504

Termination date: 20211030

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