CN114157392B - 一种分布式irs辅助通信系统安全传输的优化方法 - Google Patents

一种分布式irs辅助通信系统安全传输的优化方法 Download PDF

Info

Publication number
CN114157392B
CN114157392B CN202111436627.8A CN202111436627A CN114157392B CN 114157392 B CN114157392 B CN 114157392B CN 202111436627 A CN202111436627 A CN 202111436627A CN 114157392 B CN114157392 B CN 114157392B
Authority
CN
China
Prior art keywords
irs
base station
optimization
phase shift
beam forming
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.)
Active
Application number
CN202111436627.8A
Other languages
English (en)
Other versions
CN114157392A (zh
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.)
Anhui Normal University
Original Assignee
Anhui Normal 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 Anhui Normal University filed Critical Anhui Normal University
Priority to CN202111436627.8A priority Critical patent/CN114157392B/zh
Publication of CN114157392A publication Critical patent/CN114157392A/zh
Application granted granted Critical
Publication of CN114157392B publication Critical patent/CN114157392B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/04013Intelligent reflective surfaces
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04KSECRET COMMUNICATION; JAMMING OF COMMUNICATION
    • H04K1/00Secret communication
    • H04K1/02Secret communication by adding a second signal to make the desired signal unintelligible
    • 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/06Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station
    • H04B7/0613Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission
    • H04B7/0615Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal
    • H04B7/0617Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas at the transmitting station using simultaneous transmission of weighted versions of same signal for beam forming
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04KSECRET COMMUNICATION; JAMMING OF COMMUNICATION
    • H04K3/00Jamming of communication; Counter-measures
    • H04K3/40Jamming having variable characteristics
    • H04K3/44Jamming having variable characteristics characterized by the control of the jamming waveform or modulation type
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04KSECRET COMMUNICATION; JAMMING OF COMMUNICATION
    • H04K3/00Jamming of communication; Counter-measures
    • H04K3/80Jamming or countermeasure characterized by its function
    • H04K3/82Jamming or countermeasure characterized by its function related to preventing surveillance, interception or detection
    • H04K3/825Jamming or countermeasure characterized by its function related to preventing surveillance, interception or detection by jamming

Landscapes

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

Abstract

本发明公开一种分布式IRS辅助通信系统安全传输的优化方法,主要步骤包括:S1:构建一种由分布式IRS辅助多输入单输出的无线通信系统;S2:通信场景考虑为直通链路,基站通过分布式IRS将信息传输给系统中的合法用户;S3:基站在发送有用信号同时发送人工干扰噪声;在基站总的发射能量受限场景下,构建基于波束成形、人工噪声、以及IRS相移参数为优化参数的最大化系统安全速率和的优化方程。基于交替优化算法对基站的主动波束形成、人工噪声和IRS相移参数进行解耦合优化求解;S4:采用连续凸近似和半定限松弛来优化基站的主动波束形成、人工噪声;S5:采用流形优化或最小‑最大化算法进行求解。本发明可以显著提高通信系统的安全性。

Description

一种分布式IRS辅助通信系统安全传输的优化方法
技术领域
本发明属于无线通信技术领域,涉及一种分布式IRS辅助通信系统安全传输的优化方法。
背景技术
随着第五代(5G)无线网络的广泛部署,网络的总体架构变得越来越庞大,越来越复杂。对网络数据的速率和能效、可靠性和低延迟等的要求也越来越高。然而,目前现有的技术可能无法完全满足这些需求。鉴于此,急需开发新的技术来降低网络成本、复杂性和能源消耗。作为最具有代表性的新技术之一,智能反射面(IRS,Intelligent ReflectingSurface)通过智能地改变传播环境,在系统设计中提供了额外的自由度。IRS由大量的被动反射元件组成,可以通过可编程控制器巧妙地调整相位移,并被动地反映电磁波(EM)的相位。IRS反射元件通常由小而低成本的组件组成,它们可以有效地反映接收到的信号,而不需要专用的射频处理或再传输。此外,IRS还通常具有价廉、重量轻、易于部署的优点,这使得它们很容易在墙壁、天花板、建筑立面、广告面板等上安装或移除。
另一方面,随着无线网络环境越来越复杂,传统无线通信系统中的信息安全技术难以保证系统的安全性。但随着通信新技术的引入,利用IRS结合物理层技术解决无线接入的安全性具有很强的理论和现实意义。为此,在分布式IRS辅助无线网络通信系统中,多个IRS协同工作将基站的信息安全有效的传输到合法用户端。在基站的发射信号中增加了人工噪声(AN,Artificial Noise),用来损害窃听者的信道。这不仅提高了合法用户的接收功率,而且抑制了窃听者的接收功率,从而增强了系统物理层的安全性。
发明内容
本发明的目的在于提供一种分布式IRS辅助通信系统安全传输的优化方法。本发明将多IRS引入到无线系统中,通过联合优化基站的波束形成、人工噪声和IRS相移参数,减弱反射到窃听者的信号并增强反射到合法用户的信号,从而最大化系统总安全率(SSR,System Secrecy Rate)。
为达到上述目的,本发明提供如下一种分布式IRS辅助通信系统安全传输的优化方法,该方法包括以下步骤:
S1:构建一种由分布式IRS辅助多输入单输出的无线通信系统,系统由一个多天线基站,K个单天线用户,一个单天线的窃听者以及L个IRS构成,其中每个IRS中含有M个反射单元;
S2:通信场景考虑为直通链路(基站到用户)由于障碍物阻挡或相距较远,基站通过分布式IRS将信息传输给系统中的合法用户;
S3:为了提高系统的安全性,基站在发送有用信号同时发送人工干扰噪声(AN);在基站总的发射能量受限场景下,构建基于波束成形、人工噪声、以及IRS相移参数为优化参数的最大化系统安全速率和的优化方程。基于交替优化算法对基站的主动波束形成、人工噪声和IRS相移参数进行解耦合优化求解;
S4:给定IRS相移参数,优化主动波束形成、人工噪声优化计算:采用连续凸近似(SCA,Successive Convex Approximation)和半定限松弛(SDR,Semi-definiteRelaxation)来优化基站的主动波束形成、人工噪声;
S5:给定基站的主动波束形成、人工噪声,优化针对IRS相移参数:采用流形优化(MO,Manifold Optimization)或最小-最大化(MM,Minorize-Maximization)算法进行求解。不断交替迭代优化求解收敛。
在所述S1中,包括以下内容:
第k个用户和窃听者的接收信号分别为:
其中,Gr,l表示基站和第l个IRS之间的等效基带通道向量,分别表示IRS与用户k之间以及IRS与窃听者之间的信道向量,/>表示第l个IRS的相移矩阵,其中/>表示第l个IRS的第m个相移。nk和ne分别是第k个用户和窃听者接收端的加性高斯白噪声。基站处发送的信号表示为/>其中wk和sk分别为第k个用户的波束形成向量以及用户k的期望信号。z是人工干扰向量用来干扰窃听者。
根据香农公式,第k个合法用户和窃听者的信息速率分别为:
其中,采用单位带宽,/>分别是第k个合法用户和窃听者的噪声功率。
在所述S2中,为提高系统总安全率(SSR,System Secrecy Rate),提出优化问题:
其中,Z=zzH,Pmax是基站的最大发射功率。C1表示基站的最大发射功率约束;C2表示IRS恒模限制;C3是人工干扰信息的约束。
在所述S3中,采用交替优化算法对基站的主动波束形成、人工噪声和IRS相移参数进行解耦合,即分为两个子问题:
(1)子问题1::主动波束形成、人工噪声的优化:
s,t.C1,C3,C4. (6)
其中,
(2)子问题2::IRS相移参数的优化:
s.t.C2. (7)
其中,
在所述S4中,针对基站的主动波束形成、人工噪声的优化,固定IRS相移参数,采用连续凸近似(SCA)和半定松弛(SDR)来优化基站的主动波束形成、人工噪声:
采用连续凸近似将问题(6)转化如下:
s.t.C1,C3
C4:Rank(W)≤1. (8)
其中,
Wi和Zi为W和Z的任意可行解。
注意到问题(8)的剩余非凸性源于一级约束C4。为了解决这个问题,我们通过半定松弛(SDR)技术利用优化工具CVX来消除约束C4。从而得到最优的W和Z值。
在所述S5中,针对IRS相移参数,固定基站的主动波束形成、人工噪声,)采用流形优化(MO)或最小-最大化(MM)算法进行求解。大幅提高传统的优化算法收敛速度。不断交替迭代直到系统总保密率收敛:
将问题(7)转化为:
s.t.C2. (9)
利用交替优化迭代得到优化解,进一步通过高斯分解得到优化方程的最终解:主动波束形成、人工噪声值以及IRS相移参数。
与现有技术相比,本发明具有以下的有益技术效果:
1,本发明提供了一种分布式IRS辅助通信系统模型,解决了受高大建筑物阻挡或偏远山区的小区通信问题。
2,提供一种分布式IRS辅助通信系统安全传输方案,解决了非视距场景下安全性通信问题。
3,提供一种分布式IRS辅助通信系统安全传输的优化方法。通过算法优化,获得了较高的系统总保密率。
附图说明
图1为本发明的系统模型图;
图2为本发明的算法流程图;
图3为本发明算法与基础算法的安全性收敛比较图;
图4为本发明算法与MM算法的安全性收敛比较图。
具体实施方式
现结合附图对本发明的具体实施做进一步的详细说明:
如图1所示,分布式IRS辅助通信系统包括:一个多天线基站,K个单天线用户,一个单天线窃听者和L个IRS,其中基站含有Nt根天线,每个IRS中含有M个反射单元。采用下行传输方式。设基站与IRS,IRS与用户、窃听者之间的信道系数分别为: 所有信道均服从小尺度莱斯衰落模型。设系统每个信道带宽为单位带宽。第l个IRS的相移矩阵/>其中/>表示第l个IRS的第m个相移。为提高系统总安全率,提出优化问题:
C3:Z≥0.
其中,分别代表第k个合法用户和窃听者的信息速率,/>w为基站的主动波束形成向量,Z=zzH,其中z是人工干扰向量,由基站生成发送用来阻碍窃听者监听。假设系统采用单位带宽,/>分别是第k个用户和窃听者的噪声功率。Pmax是基站的最大发射功率。C1表示基站的最大发射功率约束;C2表示IRS恒模限制;C3是人工干扰信息的约束。
由于P是一个非凸表达式,且C2为单位模量约束。这很难直接使用现有的算法来获得该优化问题的全局最优解。因此采用交替迭代优化该问题,具体思路如下:1)先固定IRS相移参数,采用连续凸近似(SCA)和半有限松弛(SDR)来优化基站的波束形成、人工噪声;2)再固定基站的波束形成、人工噪声,采用流形优化(MO)或最小-最大化(MM)算法进行求解。不断交替迭代1,2直到系统总保密率收敛。
1)先固定IRS相移参数,采用连续凸近似(SCA)和半有限松弛(SDR)来优化基站的波束形成、人工噪声:
在给定相移矩阵Φl情况下,原问题可重新描述为:
s.t.C1,C3
C4:Rank(W)≤1..
其中,W=wwH
采用连续凸近似方法,构造G1,G2的全局解:
即对于任何可行点Wi,可微凸函数G1(W)满足以下不等式:
对于任何可行点Wi、Zi,可微凸函数G2(W,Z)满足以下不等式:
其中,为G1的全局解,/>为G2的全局解。则原问题转化为:
s.t.C1,C3,C4.
注意到该问题的剩余非凸性源于一级约束C4。为了解决这个问题,我们通过半定松弛(SDR)技术利用优化工具CVX来消除约束C4。从而得到最优的W和Z值。
2)再固定基站的主动波束形成、人工噪声,采用流形优化(MO)或最小-最大化(MM)算法进行求解。
其中,
为解决该问题,引用定理:E为任意正定矩阵,有以下函数:
其中,δ(Y)=log|Y|-Tr(YE)+M,M为无关量。Y的最优解可以表示为Yopt=E-1
由此在问题中先固定Φl值,引入变量令/> 将问题转换如下:
其中,在计算出/>值情况下,忽略常数项,简化P2
在给定值的情况下,进一步将原问题转化为:
s.t.C2.
采用流形优化(MO)或最小-最大化(MM)算法求解IRS相移参数。利用交替优化迭代得到优化解,进一步通过高斯分解得到优化方程的最终解:主动波束形成、人工噪声值以及IRS相移参数。具体的算法流程如图2。
下面结合仿真对本发明的实际效果做详细的描述。
1)仿真条件
如图1所示,设基站含有4根天线,两个IRS,每个IRS有16个智能反射单元,基站到IRS、IRS到用户、窃听者的路径损耗均为2.2。假设基站坐标为(0,0,0),IRS1坐标为(300,0,10)、IRS2坐标为(300,20,10),用户1坐标为(280,0,1.5)、用户2坐标为(300,0,1.5),窃听者坐标为(320,0,1.5)。
2)仿真结果
在本实施例中,图3,图4分别给出了本发明中采用MO算法布置IRS相移与分布式随机布置IRS相移的收敛性差异、本发明中采用MO算法布置IRS相移与采用MM算法布置IRS相移的收敛性差异,由图3可观测出:当分布式IRS的反射单元数M=16、64时本发明提出的算法所实现的系统总安全率均优于随机布置IRS相移的系统总安全率。由图4可观测出:当分布式IRS的反射单元数M=16、64时本发明提出的采用MO算法布置IRS相移与采用MM算法布置IRS相移所实现的系统总安全率的差别不大,且收敛所需迭代次数相近。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所在的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种分布式IRS辅助通信系统安全传输的优化方法,其特征在于:该方法包括以下步骤:
S1:构建一种由分布式IRS辅助多输入单输出的无线通信系统,系统由一个多天线基站,K个单天线用户,一个单天线的窃听者已经L个IRS构成,其中每个IRS中含有M个反射单元;
S2:通信场景考虑为从基站到用户的直通链路,由于障碍物阻挡或相距较远,基站通过分布式IRS将信息传输给系统中的合法用户;
S3:为了提高系统的安全性,基站在发送有用信号同时发送人工干扰噪声;在基站总的发射能量受限场景下,构建基于波束成形、人工噪声、以及IRS相移参数为优化参数的最大化系统安全速率和的优化方程,基于交替优化算法对基站的主动波束形成、人工噪声和IRS相移参数进行解耦合优化求解,为提高系统总安全率,提出优化问题:
C3:Z≥0.
其中,Z=zzH,Pmax是基站的最大发射功率,C1表示基站的最大发射功率约束;C2表示IRS恒模限制;C3是人工干扰信息的约束;
采用交替优化算法对基站的主动波束形成、人工噪声和IRS相移参数进行解耦合,即分为两个子问题:
(1)子问题1::主动波束形成、人工噪声的优化:
s.t.C1,C3. (6)
其中,
(2)子问题2::IRS相移参数的优化:
s.t.C2. (7)
其中,
S4:给定IRS相移参数,优化主动波束形成、人工噪声优化计算:采用连续凸近似和半定限松弛来优化基站的主动波束形成、人工噪声;
S5:给定基站的主动波束形成、人工噪声,优化针对IRS相移参数:采用流形优化或最小-最大化算法进行求解,不断交替迭代优化求解收敛。
2.根据权利要求1所述的一种分布式IRS辅助通信系统安全传输的优化方法,其特征在于:所述S1中,包括以下内容:
第k个合法用户和窃听者的接收信号分别为:
其中,Gr,l表示基站和第l个IRS之间的等效基带通道向量,分别表示IRS与用户k之间以及IRS与窃听者之间的信道向量,/>表示第l个IRS的相移矩阵,其中/>表示第l个IRS的第m个相移,nk和ne分别是第k个用户和窃听者接收端的加性高斯白噪声,基站处发送的信号表示为/>其中wk和sk分别为第k个用户的波束形成向量以及用户k的期望信号,z是人工干扰向量用来干扰窃听者,
根据香农公式,第k个合法用户和窃听者的信息速率分别为:
其中,采用单位带宽,/>分别是第k个合法用户和窃听者的噪声功率。
3.根据权利要求1所述的一种分布式IRS辅助通信系统安全传输的优化方法,其特征在于:在所述S4中,针对基站的主动波束形成、人工噪声的优化,固定IRS相移参数,采用连续凸近似(SCA)和半定松弛(SDR)来优化基站的主动波束形成、人工噪声:
采用连续凸近似将问题(6)转化如下:
s.t.C1,C3
C4:Rank(W)≤1. (8)
其中,
和Zi为W和Z的任意可行解。
4.根据权利要求1所述的一种分布式IRS辅助通信系统安全传输的优化方法,其特征在于:在所述S5中,针对IRS相移参数,固定基站的主动波束形成、人工噪声,采用流形优化进行求解,大幅提高传统的优化算法收敛速度,不断交替迭代直到系统总保密率收敛:
将问题(7)转化为:
s.t.C2.(9)
利用交替优化迭代得到优化解,进一步通过高斯分解得到优化方程的最终解:主动波束形成、人工噪声值以及IRS相移参数。
CN202111436627.8A 2021-11-30 2021-11-30 一种分布式irs辅助通信系统安全传输的优化方法 Active CN114157392B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111436627.8A CN114157392B (zh) 2021-11-30 2021-11-30 一种分布式irs辅助通信系统安全传输的优化方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111436627.8A CN114157392B (zh) 2021-11-30 2021-11-30 一种分布式irs辅助通信系统安全传输的优化方法

Publications (2)

Publication Number Publication Date
CN114157392A CN114157392A (zh) 2022-03-08
CN114157392B true CN114157392B (zh) 2023-11-28

Family

ID=80784379

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111436627.8A Active CN114157392B (zh) 2021-11-30 2021-11-30 一种分布式irs辅助通信系统安全传输的优化方法

Country Status (1)

Country Link
CN (1) CN114157392B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115314094B (zh) * 2022-05-30 2024-05-24 安徽师范大学 一种基于级联IRSs与合法用户协同安全传输策略分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112073102A (zh) * 2020-09-08 2020-12-11 郑州大学 一种基于irs的安全波束成形方法与装置
CN113037346A (zh) * 2021-03-12 2021-06-25 重庆邮电大学 Irs与人工噪声辅助的mimo系统物理层安全设计方法
CN113315547A (zh) * 2021-05-28 2021-08-27 北京邮电大学 智能反射面辅助多小区的鲁棒联合传输波束优化方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9537556B2 (en) * 2014-07-11 2017-01-03 Huawei Technologies Canada Co., Ltd. Systems and methods for optimized beamforming and compression for uplink MIMO cloud radio access networks
CN111355520B (zh) * 2020-03-10 2022-03-08 电子科技大学 一种智能反射表面辅助的太赫兹安全通信系统设计方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112073102A (zh) * 2020-09-08 2020-12-11 郑州大学 一种基于irs的安全波束成形方法与装置
CN113037346A (zh) * 2021-03-12 2021-06-25 重庆邮电大学 Irs与人工噪声辅助的mimo系统物理层安全设计方法
CN113315547A (zh) * 2021-05-28 2021-08-27 北京邮电大学 智能反射面辅助多小区的鲁棒联合传输波束优化方法

Also Published As

Publication number Publication date
CN114157392A (zh) 2022-03-08

Similar Documents

Publication Publication Date Title
CN112243252B (zh) 一种面向无人机中继网络的安全传输增强方法
CN113507304B (zh) 一种智能反射面辅助的无人机安全通信方法
CN113949427B (zh) 一种多用户无线网络安全能效优化设计方法和系统
CN113098652A (zh) 一种智能反射面辅助的功率域noma安全通信系统设计方法
CN114222289A (zh) 一种智能反射面辅助全双工无线携能网络保密通信方法
CN116015374B (zh) 智能反射面辅助的非正交多址接入携能网络安全波束赋形方法
CN114665937B (zh) 一种多输入多输出收发机设计方法及装置
CN115021792B (zh) 一种无线通信系统安全传输方法
CN113938891A (zh) 一种反射面辅助的用户节点不可信noma网络安全通信方法
CN113490223A (zh) 一种基于irs的无线供能反向散射通信方法及系统
CN114006645A (zh) 一种认知无人机中继辅助的安全传输方法及系统
CN115484607A (zh) 一种ris辅助的swipt无线系统安全通信方法
CN114157392B (zh) 一种分布式irs辅助通信系统安全传输的优化方法
CN114626229A (zh) 智能反射面辅助的无线通信与能量收集系统性能分析方法
CN114268351A (zh) 基于无人机智能反射面辅助的安全通信方法
CN117081636A (zh) 可重构智能表面辅助主动干扰的发射功率优化方法及装置
CN114337902B (zh) 一种irs辅助的毫米波多小区间干扰的抑制方法
CN116095699A (zh) 利用双边缘计算的高安全卸载方法、系统、终端及其介质
CN115189730A (zh) 双智能可重构表面辅助双向方向调制的功率分配设计方法
Peng et al. Enhancing underwater sensor network security with coordinated communications
CN113891286A (zh) 无人机辅助通信系统的保密通信性能优化装置及方法
Wang Optimal joint dual transmitter receiver diversity and power control for wireless networks
Li et al. Active RIS‐Assisted Transmission Design for Wireless Secrecy Network with Energy Harvesting
CN115208491B (zh) 一种基于分布式智能反射面的功率域noma通信系统设计方法
Zheng et al. Secure Offloading in NOMA-Enabled Multi-Access Edge Computing Networks

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