CN112929307A - 一种免授权通信下的分布式信道估计方法 - Google Patents
一种免授权通信下的分布式信道估计方法 Download PDFInfo
- Publication number
- CN112929307A CN112929307A CN202110217428.1A CN202110217428A CN112929307A CN 112929307 A CN112929307 A CN 112929307A CN 202110217428 A CN202110217428 A CN 202110217428A CN 112929307 A CN112929307 A CN 112929307A
- Authority
- CN
- China
- Prior art keywords
- channel estimation
- user
- channel
- parameters
- distributed
- 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
Links
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04L—TRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
- H04L25/00—Baseband systems
- H04L25/02—Details ; arrangements for supplying electrical power along data transmission lines
- H04L25/0202—Channel estimation
- H04L25/024—Channel estimation channel estimation algorithms
- H04L25/0242—Channel estimation channel estimation algorithms using matrix methods
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Power Engineering (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- Mobile Radio Communication Systems (AREA)
Abstract
本发明提供一种免授权通信下的分布式信道估计方法,包括以下步骤:S1:构建广义线性模型,并初始化信道估计参数和用户检测参数;S2:采用矢量形式的分布式广义期望传播算法求解广义线性模型,得到估计信道,并对信道估计参数进行更新,对用户检测与信道估计之间的信息进行转换,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测,得到活跃用户,并对用户检测参数进行更新;S3:判断是否达到预设的迭代结束条件;若是,则输出当前轮迭代得到的估计信道和活跃用户;若否,则返回步骤S2进行下一轮迭代。本发明提供一种免授权通信下的分布式信道估计方法,解决了目前信道估计得到的通信信道与真实信道差别较大的问题。
Description
技术领域
本发明涉及通信信道传输特性估计的技术领域,更具体的,涉及一种免授权通信下的分布式信道估计方法。
背景技术
通信信道通常指无线信号传输的媒介,比如生活中手机与运营商基站之间信号传输便存在通信信道,其受环境与距离位置等因素的影响。不同设备与基站之间真实的通信信道是不可知的,而通信信道是求解接收端信号的关键信息,因此信道估计通常作为建立通信的第一步,同时也是信号传输的关键一步,信道估计的准确性会影响后续信号传输的效率。
为了解决信道估计的问题,目前设备在与基站进行信号传输之前会通过发送特定信号序列来建立传输信道,理想情况下基站在接收到特定信号后通过简单计算便可得到通信信道的特性。但在实际应用中,设备发送序列内部通常具有强烈的相关性,比如zadoff-chu序列,由发送序列组成的矩阵在信道估计中作为一个重要的部分,现有算法要求该矩阵中的每一个元素即信号与信号之间无任何关系,如果不满足该要求会导致算法对信道估计的能力下降,导致计算得到的通信信道与真实信道差别较大,已不能满足信号传输的基本要求。
现有技术中,如2019-07-30公开的中国专利,一种自适应开销的活跃用户检测和信道估计方法,公开号为CN110071881A,利用估计到的活跃用户集合及其信道,通过计算接收信号误差评估估计质量,可以根据现实系统中活跃用户数量和信道环境,自适应调整帧结构中导频的时隙开销,以尽可能低的接入时延保证服务质量,实现超可靠的活跃用户检测和信道估计,但无法满足不同场景下的应用需求。
发明内容
本发明为克服目前信道估计得到的通信信道与真实信道差别较大的技术缺陷,提供一种免授权通信下的分布式信道估计方法。
为解决上述技术问题,本发明的技术方案如下:
一种免授权通信下的分布式信道估计方法,包括以下步骤:
S1:构建免授权通信中信道估计的广义线性模型(Generalized Linear Model),并初始化信道估计参数和用户检测参数;
S2:根据信道估计参数,对信道估计模型中矢量的消息采用矢量形式的分布式广义期望传播算法求解广义线性模型,得到估计信道,并对信道估计参数进行更新,
对用户检测与信道估计之间的信息进行转换,其中,消息从用户检测部分传到信道估计部分时使用标量消息到矢量消息的消息传递规则,消息从信道估计部分传到用户检测部分时使用矢量消息到标量消息的消息传递规则,
根据用户检测参数及用户检测模型中标量的消息,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测,得到活跃用户,并对用户检测参数进行更新;
S3:判断是否达到预设的迭代结束条件;
若是,则输出当前轮迭代得到的估计信道和活跃用户;
若否,则返回步骤S2进行下一轮迭代。
优选的,在步骤S1中,根据应用场景对信道估计准确度和计算时间的要求初始化信道估计参数和用户检测参数。
优选的,所述分布式广义期望传播算法的参数根据应用场景对信道估计准确度和计算时间的要求确定。
优选的,通过以下公式对信道估计部分与用户检测部分的信息进行转换:
μl←nl(xnl)表示从节点传输到节点xl的信息,μl→nl(xnl)表示从节点xl传输到节点的信息;μl←kl(xl)表示从节点δ(Zkl-Akxl)传输到节点xl的消息,μl→kl(xl)表示从节点xl传输到节点δ(Zkl-Akxl)的信息。
优选的,所述广义线性模型为:
其中,为接收信号,Qc(·)为复数的低精度量化器,导频矩阵 具有列内相关性,即来自于第n个用户的导频序列sn具有列内相关性,W是加性高斯白噪声,混合信道矩阵X=(x1,x2,…xN)T中的xn表示第n个用户与基站之间的传输损耗。
优选的,采用Ξn表示第n个用户是否活跃的随机性并用ξn表示其实现:
因此,混合信道矩阵X的每一个分量都相互独立且服从条件分布:
优选的,在步骤S2中,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测具体为:通过求解第n个用户ξn的后验概率并根据后验概率判断ξn是否为活跃用户。
优选的,所述后验概率为P(ξn|Y),表示在基站接收到信息的条件下用户活跃的概率;定义P(ξn=1|Y)表示第n个用户ξn为活跃用户的概率,P(ξn=0|Y)表示第n个用户ξn为不活跃用户的概率,其中,Y表示基站接收到的信息。
优选的,当后验概率为P(ξn=1|Y)>0.5时认为用户ξn为活跃用户,否则认为用户ξn为不活跃用户。
优选的,在步骤S3中,预设的迭代结束条件为:每轮迭代估计的信道准确度不再提高或提高小于预设的阈值。
与现有技术相比,本发明技术方案的有益效果是:
本发明提供了一种免授权通信下的分布式信道估计方法,通过结合标量形式的信念传播算法和矢量形式的分布式广义期望传播算法进行不断的迭代计算,实现准确的信道估计,能够满足发送序列具有相关性条件下的信道估计需求,且在该条件下可以正常工作并准确估计通信信道。
附图说明
图1为本发明的技术方案实施步骤流程图;
图2为本发明中信道估计部分与用户检测部分之间的信息交互示意图;
图3为本发明的通信模型示意图;
图4为本发明中信道估计的广义线性模型示意图;
图5为本发明中Hybrid DeGEC算法与现有的Hybrid GAMP算法的对比示意图;
图6为本发明中Hybrid DeGEC算法在不同分块数与现有的信道估计准确度的对比示意图;
图7为本发明中Hybrid DeGEC算法在不同分块数下的信道估计准确度与每次迭代所需时间的对比示意图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示,一种免授权通信下的分布式信道估计方法,包括以下步骤:
S1:构建免授权通信中信道估计的广义线性模型,并初始化信道估计参数和用户检测参数;
S2:根据信道估计参数,对信道估计模型中矢量的消息采用矢量形式的分布式广义期望传播算法求解广义线性模型,得到估计信道,并对信道估计参数进行更新,
对用户检测与信道估计之间的信息进行转换,其中,消息从用户检测部分传到信道估计部分时使用标量消息到矢量消息的消息传递规则,消息从信道估计部分传到用户检测部分时使用矢量消息到标量消息的消息传递规则,
根据用户检测参数及用户检测模型中标量的消息,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测,得到活跃用户,并对用户检测参数进行更新;在实际实施时,步骤S2还可以为:根据用户检测参数及用户检测模型中标量的消息,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测,得到活跃用户,并对用户检测参数进行更新,
对用户检测与信道估计之间的信息进行转换,其中,消息从用户检测部分传到信道估计部分时使用标量消息到矢量消息的消息传递规则,消息从信道估计部分传到用户检测部分时使用矢量消息到标量消息的消息传递规则,
根据信道估计参数,对信道估计模型中矢量的消息采用矢量形式的分布式广义期望传播算法求解广义线性模型,得到估计信道,并对信道估计参数进行更新;
S3:判断是否达到预设的迭代结束条件;
若是,则输出当前轮迭代得到的估计信道和活跃用户;
若否,则返回步骤S2进行下一轮迭代。
在具体实施过程中,通过将标量形式的信念传播算法BP和矢量形式的分布式广义期望传播算法DeGEC结合形成Hybrid DeGEC算法在免授权通信下进行分布式信道估计,能够满足发送序列具有相关性条件下的信道估计需求,且在该条件下可以正常工作并准确估计通信信道。
实施例2
更具体的,在步骤S1中,根据应用场景对信道估计准确度和计算时间的要求初始化信道估计参数和用户检测参数,从而灵活调整计算复杂度和估计准确度,以适应不同的应用环境条件。
更具体的,所述分布式广义期望传播算法的参数根据应用场景对信道估计准确度和计算时间的要求确定。
在具体实施过程中,通过调整分布式广义期望传播算法的分块数量参数,从而灵活调整信道估计的准确度和所需的计算时间,能够满足不同场景下的应用需求。
更具体的,通过以下公式对信道估计部分与用户检测部分的信息进行转换:
其中,信道估计部分的矢量消息通过公式转换为用户检测部分的标量消息,其方法是提取所有信道估计部分的矢量消息的相同分量μl←kl(xnl)并进行连乘获得新的标量消息μl←nl(xnl);用户检测部分的标量消息通过公式转换为信道估计部分的矢量消息,其方法是将所有用户检测部分的标量消息μl→nl(xnl)与除去相同方向信道估计部分μl←kl(xl)的矢量消息∏j≠kμl←jl(xl)进行连乘获得新的矢量消息μl→kl(xl);
μl←nl(xnl)表示从节点传输到节点xl的信息,μl→nl(xnl)表示从节点xl传输到节点的信息;μl←kl(xl)表示从节点δ(Zkl-Akxl)传输到节点xl的消息,μl→kl(xl)表示从节点xl传输到节点δ(Zkl-Akxl)的信息。如图2所示为信道估计部分与用户检测部分之间的信息交互。
在具体实施过程中,信道估计部分里传输的数据是矢量,而用户检测部分里传输的数据是标量,两个部分的数据维度并不相同,因此需要通过以上公式完成信道估计部分与用户检测部分之间的消息转换。
更具体的,在实际通信中,如图3所示,一个基站覆盖的范围内通常有多个设备,但同一时间内只有少量设备处于活跃状态,每个设备分配了特定的导频序列,活跃设备同时发送导频序列,基站通过接收得到的导频序列来恢复用户与基站之间的传输信道以及活跃用户的检测;如图4所示,所述广义线性模型为:
更具体的,采用Ξn表示第n个用户是否活跃的随机性并用ξn表示其实现:
因此,混合信道矩阵X的每一个分量都相互独立且服从条件分布:
更具体的,在步骤S2中,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测具体为:通过求解第n个用户ξn的后验概率并根据后验概率判断ξn是否为活跃用户。
更具体的,所述后验概率为P(ξn|Y),表示在基站接收到信息的条件下用户活跃的概率;定义P(ξn=1|Y)表示第n个用户ξn为活跃用户的概率,P(ξn=0|Y)表示第n个用户ξn为不活跃用户的概率,其中,Y表示基站接收到的信息。
更具体的,当后验概率为P(ξn=1|Y)>0.5时认为用户ξn为活跃用户,否则认为用户ξn为不活跃用户。
更具体的,在步骤S3中,预设的迭代结束条件为:每轮迭代估计的信道准确度不再提高或提高小于预设的阈值。
在具体实施过程中,还可以将迭代结束条件预设为达到最大的迭代轮数。
实施例3
Hybrid DeGEC算法的伪代码如下:
算法的执行首先需要根据使用环境对信道估计准确度和计算所需时间的要求选择不同的初始化参数K,对于信道估计准确度要求较高而计算时间要求不严格的场景K选择较小的值,对于要求计算时间减少但估计准确度不是很高的场景选择较大的k值。
信道估计算法的迭代过程首先从信道估计部分开始计算,即完成公式(4)到公式(17),随后在信道估计部分与用户检测部分的信息交互处完成信息从信道估计部分到用户检测部分的转换,即公式(18),之后在用户检测部分进行计算,并将得到的信息通过公式(2)和公式(3)传回信道估计部分,完成一轮迭代。
实施例4
图5为本实施例的Hybrid DeGEC算法与现有的Hybrid GAMP算法因子图对比。在本实施例中采用Hybrid DeGEC算法测试了具有30根天线、800个用户、发送信号序列长度为400并具有相关性、活跃用户比例为10%时的分布式MIMO通信情况,结果如图6所示,横坐标是迭代次数,纵坐标是估计信道与真实信道的差距,由图6可知,本实施例中的信道估计准确度比现有的Hybrid GAMP算法要高,并且通过调整本实施例中的设定参数分块数量K,可以获得不同信道估计准确度。除此之外,在图7中,横坐标为本实施例的Hybrid DeGEC算法中的参数分块数量K的不同取值,左边的纵坐标是每次循环迭代所需要的计算时间,右边的纵坐标是在第四轮迭代之后本实施例的Hybrid DeGEC算法的信道估计准确度;从图6可知参数K取值越小,信道估计的准确度越高,但所需要的计算时间也越长;参数K取值越大,信道估计的准确度下降,计算所需要的时间的呈指数下降。综合上述测试结果,Hybrid DeGEC算法可以通过调节参数K满足不同环境下对信道估计准确度与计算时间的要求。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (10)
1.一种免授权通信下的分布式信道估计方法,其特征在于,包括以下步骤:
S1:构建免授权通信中信道估计的广义线性模型,并初始化信道估计参数和用户检测参数;
S2:根据信道估计参数,对信道估计模型中矢量的消息采用矢量形式的分布式广义期望传播算法求解广义线性模型,得到估计信道,并对信道估计参数进行更新,
对用户检测与信道估计之间的信息进行转换,
根据用户检测参数及用户检测模型中标量的消息,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测,得到活跃用户,并对用户检测参数进行更新;
S3:判断是否达到预设的迭代结束条件;
若是,则输出当前轮迭代得到的估计信道和活跃用户;
若否,则返回步骤S2进行下一轮迭代。
2.根据权利要求1所述的一种免授权通信下的分布式信道估计方法,其特征在于,在步骤S1中,根据应用场景对信道估计准确度和计算时间的要求初始化信道估计参数和用户检测参数。
3.根据权利要求1所述的一种免授权通信下的分布式信道估计方法,其特征在于,所述分布式广义期望传播算法的参数根据应用场景对信道估计准确度和计算时间的要求确定。
7.根据权利要求1所述的一种免授权通信下的分布式信道估计方法,其特征在于,在步骤S2中,采用标量形式的信念传播算法对免授权通信下的用户进行活跃检测具体为:通过求解第n个用户ξn的后验概率并根据后验概率判断ξn是否为活跃用户。
8.根据权利要求7所述的一种免授权通信下的分布式信道估计方法,其特征在于,所述后验概率为P(ξn|Y),表示在基站接收到信息的条件下用户活跃的概率;定义P(ξn=1|Y)表示第n个用户ξn为活跃用户的概率,P(ξn=0|Y)表示第n个用户ξn为不活跃用户的概率,其中,Y表示基站接收到的信息。
9.根据权利要求8所述的一种免授权通信下的分布式信道估计方法,其特征在于,当后验概率为P(ξn=1|Y)>0.5时认为用户ξn为活跃用户,否则认为用户ξn为不活跃用户。
10.根据权利要求1所述的一种免授权通信下的分布式信道估计方法,其特征在于,在步骤S3中,预设的迭代结束条件为:每轮迭代估计的信道准确度不再提高或提高小于预设的阈值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110217428.1A CN112929307B (zh) | 2021-02-26 | 2021-02-26 | 一种免授权通信下的分布式信道估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110217428.1A CN112929307B (zh) | 2021-02-26 | 2021-02-26 | 一种免授权通信下的分布式信道估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112929307A true CN112929307A (zh) | 2021-06-08 |
CN112929307B CN112929307B (zh) | 2022-03-25 |
Family
ID=76172239
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110217428.1A Active CN112929307B (zh) | 2021-02-26 | 2021-02-26 | 一种免授权通信下的分布式信道估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112929307B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114024638A (zh) * | 2021-11-04 | 2022-02-08 | 哈尔滨工业大学 | 基于空间调制的上行免授权系统接收方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005074147A1 (en) * | 2004-01-30 | 2005-08-11 | Universite Laval | Multi-user adaptive array receiver and method |
CN107769824A (zh) * | 2017-10-27 | 2018-03-06 | 清华大学 | 联合多天线的稀疏贝叶斯学习的用户检测方法及系统 |
CN108540410A (zh) * | 2018-03-06 | 2018-09-14 | 东南大学 | 免调度传输系统联合活跃用户检测和信道估计方法 |
CN110011775A (zh) * | 2019-03-22 | 2019-07-12 | 安徽师范大学 | 联合实现活动用户检测及其信道估计方法及其系统 |
CN110071881A (zh) * | 2019-04-26 | 2019-07-30 | 北京理工大学 | 一种自适应开销的活跃用户检测和信道估计方法 |
CN111865845A (zh) * | 2020-06-04 | 2020-10-30 | 国网江西省电力有限公司供电服务管理中心 | 一种mimo用户检测与信道估计装置及方法 |
-
2021
- 2021-02-26 CN CN202110217428.1A patent/CN112929307B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005074147A1 (en) * | 2004-01-30 | 2005-08-11 | Universite Laval | Multi-user adaptive array receiver and method |
CN107769824A (zh) * | 2017-10-27 | 2018-03-06 | 清华大学 | 联合多天线的稀疏贝叶斯学习的用户检测方法及系统 |
CN108540410A (zh) * | 2018-03-06 | 2018-09-14 | 东南大学 | 免调度传输系统联合活跃用户检测和信道估计方法 |
CN110011775A (zh) * | 2019-03-22 | 2019-07-12 | 安徽师范大学 | 联合实现活动用户检测及其信道估计方法及其系统 |
CN110071881A (zh) * | 2019-04-26 | 2019-07-30 | 北京理工大学 | 一种自适应开销的活跃用户检测和信道估计方法 |
CN111865845A (zh) * | 2020-06-04 | 2020-10-30 | 国网江西省电力有限公司供电服务管理中心 | 一种mimo用户检测与信道估计装置及方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114024638A (zh) * | 2021-11-04 | 2022-02-08 | 哈尔滨工业大学 | 基于空间调制的上行免授权系统接收方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112929307B (zh) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
KR101681617B1 (ko) | 저밀도 확산 변조 검출을 위한 시스템 및 방법 | |
CN110233653B (zh) | 基于加权集成聚类算法的mimo系统的盲多径识别方法及系统 | |
CN109327850B (zh) | 基于梯度追踪和多步拟牛顿法技术的非正交多址接入系统多用户检测方法 | |
KR100975313B1 (ko) | 다중 입출력 무선통신 시스템에서 다변 다항식을 이용한신호검출 장치 및 방법 | |
Zhu et al. | A message-passing approach for joint channel estimation, interference mitigation, and decoding | |
Choi et al. | Robust learning-based ML detection for massive MIMO systems with one-bit quantized signals | |
Chataut et al. | Efficient and low complex uplink detection for 5G massive MIMO systems | |
CN112929307B (zh) | 一种免授权通信下的分布式信道估计方法 | |
Wen et al. | Joint channel-and-data estimation for large-MIMO systems with low-precision ADCs | |
Al-Jarrah et al. | Decision fusion in mobile wireless sensor networks using cooperative multiple symbol differential space time coding | |
Wei et al. | Message passing receiver design for uplink grant-free SCMA | |
Wang et al. | On massive IoT connectivity with temporally-correlated user activity | |
KR20100035369A (ko) | 신호 전송 제어 방법 및 장치 | |
CN108566227A (zh) | 一种多用户检测方法 | |
Idowu-Bismark et al. | A primer on MIMO detection algorithms for 5G communication network | |
CN107465472B (zh) | 一种基于路径合成的多径时延估计方法 | |
CN115967421A (zh) | 一种分布式天线系统下联合参数估计与信号重构的方法及系统 | |
Srinivas et al. | A semi-blind based channel estimator for pilot contaminated one-bit massive MIMO systems | |
Di Renzo et al. | Performance of cooperative multi-hop wireless systems over log-normal fading channels | |
Niazadeh et al. | An alternating minimization method for sparse channel estimation | |
Takeda et al. | Statistical mechanical analysis of the linear vector channel in digital communication | |
CN110086557B (zh) | 一种免调度多用户检测方法 | |
Al-Dulaimi et al. | Cooperative Spectrum Sensing for Three Secondary Users with Sequential Relaying for Cognitive Radio | |
CN116032317B (zh) | 一种免授权的联合活跃用户与数据检测方法 | |
David et al. | Feedback-controlled channel estimation with low-resolution ADCs in multiuser MIMO systems |
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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230918 Address after: 510000 rooms 704 and 705, building 3, No. 18, Shenzhou Road, Huangpu District, Guangzhou, Guangdong Province Patentee after: Guangzhou longjianda Electronics Co.,Ltd. Address before: 510090 Dongfeng East Road 729, Yuexiu District, Guangzhou City, Guangdong Province Patentee before: GUANGDONG University OF TECHNOLOGY |