CN101815305A - 基于可信度的协同频谱感知方法 - Google Patents
基于可信度的协同频谱感知方法 Download PDFInfo
- Publication number
- CN101815305A CN101815305A CN201010109711A CN201010109711A CN101815305A CN 101815305 A CN101815305 A CN 101815305A CN 201010109711 A CN201010109711 A CN 201010109711A CN 201010109711 A CN201010109711 A CN 201010109711A CN 101815305 A CN101815305 A CN 101815305A
- Authority
- CN
- China
- Prior art keywords
- lambda
- gamma
- local
- sigma
- confidence level
- 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
Landscapes
- Mobile Radio Communication Systems (AREA)
Abstract
本发明的基于可信度的协同频谱感知方法,包括初始化、未授权用户本地感知、本地判决、感知信息量化、可信度融合、判决融合等步骤。本发明在先验信息未知的条件下给出了未授权用户本地感知的可信度,采用基站与未授权用户联合的混合门限判决策略,并通过量化本地感知结果来减小协同感知的网络开销,具有很好的感知性能。
Description
技术领域
本发明属于无线通信技术领域,涉及一种数字通信领域中的频谱感知实现方法,特别是在认知无线网络(cognitive radio network,CRN)中报告信道带宽受限情况下可采用的感知方法,具体地说是认知无线网络中未授权用户使用的基于可信度的协同频谱感知方法。
背景技术
无线电通信频谱是一种重要的国家资源,一般都是由国家机构控制并分配给授权用户。然而在这种传统的无线电通信系统中,频谱的分配是固定的。授权用户对其频谱具有独占性,其他用户不允许使用,这种分配方式有利于保证系统的服务质量(QoS)。由于频谱资源在不同位置不同时间段的利用有所不同,这种静态的分配方式浪费资源,降低了频谱的使用率。另一方面,无线服务如移动通信、公共安全、电视广播等的迅速发展对无线频谱资源的需求日益增大,无线频谱资源的缺乏问题日益严重。
认知无线电是一种用于提高无线电通信频谱利用率的新的智能技术。未授权用户可以感知周围的环境,利用已经分配给授权用户、但在某一特定的时刻和环境下并没有被占用的频带;能够根据输入激励的变化实时地调整其传输参数,在有限信号空间中以最优的方式有效传送信息,以实现无论何时何地都能保证通信的高可靠性和无线频谱的高利用率。认知无线电的首要任务是感知频谱环境,即频谱感知。
目前,单个未授权用户的频谱感知主要有能量检测、匹配滤波器检测、波形检测、循环谱检测和小波检测等几种方法。但由于实际信道中的多径和阴影效应,单个未授权用户频谱感知的性能并不乐观,因此需要靠不同用户间的协同频谱感知来对抗多径和阴影效应。目前协同频谱感知中对多个未授权用户的感知信息处理方式主要有两种:软信息融合和硬信息融合,软信息融合是指未授权用户的本地感知结果不做任何处理,直接发给基站,由基站进行融合处理;硬信息融合是指未授权用户对本地感知结果进行判决后发送判决结果给基站,基站再进行融合。软信息融合策略的感知性能较好,但是它比硬信息融合需要更大的感知开销。因此,在协同频谱感知中,如何有效利用未授权用户的本地感知信息,同时又能控制协同频谱感知所需的感知开销已成为一个重要的研究问题。现有技术无法实现报告信道带宽受限条件下的有效协同频谱感知。
发明内容
本发明的目的是针对报告信道带宽受限条件下无法实现频谱感知性能与感知开销很难折中的问题,提出一种基于可信度的协同频谱感知方法。通过有效利用多个未授权用户产生的本地感知信息,该方法能够实现频谱感知性能好,同时又满足感知开销要求,克服当前认知无线网络中协同频谱感知性能与感知开销之间难以同时满足要求的问题。
本发明的技术方案是:
一种基于可信度的协同频谱感知方法,它包括以下步骤:
步骤一:未授权用户在本地使用基于能量检测的频谱感知方法进行感知得到感知结果 其中y是未授权用户的接收信号,θ是未授权用户的本地感知结果,u是时间带宽积;
步骤二:
(a).设定本方法的本地虚警概率Qf C,选择软信息融合方法的虚警概率Qf S,根据 计算软信息融合方法的判决门限λ0,其中 N是系统中参与协同频谱感知的未授权用户数,然后检验Qf S是否满足以下两个条件;
其中Qf OR(λ′)是硬信息融合方法的虚警概率, 是硬信息融合方法的本地判决门限;
若Qf S满足检验条件,转步骤(b);否则,转步骤(a)重新选择选择软信息融合方法的虚警概率Qf S;
(b).根据设定的本地虚警概率Qf C和软信息融合方法的判决门限λ0计算未授权用户本地感知判决门限的均值为和基站可信度融合门限其中L是参与可信度融合的未授权用户数;得到未授权用户的两个本地判决门限分别为 Δλ为两个本地判决门限λ1,λ2与均值之差,它满足:
其中 σ2为未授权用户的本地噪声方差;
步骤三:未授权用户进行本地判决x=Φ(θ),判决规则如下:
其中H1代表授权用户存在,H0代表授权用户不存在, 是在本地感知结果θ落在本地判决门限λ1、λ2之间区域时对信号存在的可信度;
步骤四:未授权用户将本地感知和判决的处理结果通过报告信道传给基站,
(a).若本地判决为H1或H0,则直接将1比特判决结果即0或1发送给基站;
(b).若本地判决得到的是本地感知可信度αi,则将可信度进行M比特量化后再上传;其中,可信度M比特量化即 quantize(X,M)表示对X进行M比特量化;
步骤五:基站对收到的未授权用户本地判决结果xi按照长度进行分类并进行可信度融合,即基站将接收到的结果按照数据长度分为两类,即1比特数据b和M比特数据α,然后对所有M比特的可信度进行合并得到可信度融合结果ω,
步骤六:
本发明的有益效果:
本发明在报告信道带宽受限即感知开销受限的情况下,针对现有协同频谱感知中感知性能与感知开销难以调和的问题,在未授权用户进行本地判决时考虑了本地感知的可信度,未授权用户本地对感知信息进行预处理,采用基站与未授权用户联合的混合门限判决策略,使得协同频谱感知的性能和感知开销达到很好地折中。在授权用户先验信息未知的条件下,通过使用本发明中的选择本地判决门限和基站融合策略的方法,能够以远小于软信息融合策略的感知开销,达到与其接近的感知性能。
附图说明
图1为本发明所采用的系统模型。
图2为基于可信度的协同感知方法流程图。
图3为未量化情况下CCSS方法的仿真曲线。
图4为量化情况下CCSS方法的仿真曲线。
图5为CCSS方法的感知开销仿真曲线。
图6为本地判决门限的示意图。
具体实施方式
下面结合附图和实施例对本发明作进一步的说明。
如图1所示,一种认知无线网络中基于可信度的协同频谱感知方法,它包括以下步骤:
步骤一:未授权用户在本地使用基于能量检测的频谱感知方法进行感知得到感知结果 其中y是未授权用户的接收信号,θ是未授权用户的本地感知结果,u是时间带宽积;
该步骤中,所述的基于能量检测的频谱感知方法,具体是:未授权用户对本地信号经过一段时间的观察,累计得到能量值,未授权用户与设定的门限比较判决得到本地感知结果。
步骤二:
(a).设定本方法的本地虚警概率Qf C,选择软信息融合方法的虚警概率Qf S,根据 计算软信息融合方法的判决门限λ0,其中 N是系统中参与协同频谱感知的未授权用户数,然后检验Qf S是否满足以下两个条件;
其中Qf OR(λ′)是硬信息融合方法的虚警概率, 是硬信息融合方法的本地判决门限;
若Qf S满足检验条件,转步骤(b);否则,转步骤(a)重新选择选择软信息融合方法的虚警概率Qf S;
(b).根据设定的本地虚警概率Qf C和软信息融合方法的判决门限λ0计算未授权用户本地感知判决门限的均值为数值上等于硬信息融合方法的本地判决门限,则基站可信度融合门限为其中L是参与可信度融合的未授权用户数;得到未授权用户的两个本地判决门限分别为 Δλ为两个本地判决门限λ1,λ2与均值之差,如图6所示,它满足:
其中 σ2为未授权用户的本地噪声方差;
该步骤中,所述的计算本地判决门限的方法是指首先选择一个合适的软判决融合方法虚警概率值,利用软判决融合方法的虚警概率计算出两个本地判决门限的均值,然后根据基于可信度方法的虚警概率计算公式计算两个本地门限值λ1,λ2。
步骤三:未授权用户进行本地判决x=Φ(θ),判决规则如下:
其中H1代表授权用户存在,H0代表授权用户不存在, 是在本地感知结果θ落在本地判决门限λ1、λ2之间区域时对信号存在的可信度;
该步骤中,所述的未授权用户将本地感知的结果和本地判决门限λ1,λ2比较进行本地判决,是指若本地感知结果落在两个本地判决门限之间的区域,得到的是本地可信度;若本地感知结果落在两个判决门限之间的区域之外,则得到的是1比特本地判决结果;
未授权用户得到的本地感知可信度,具体是指在本地感知结果落在两个判决门限之间区域时对信号存在的可信程度α,定义为归一化的 值,其中θ是本地感知结果,λ1和λ2是本地判决门限。
步骤四:未授权用户将本地感知和判决的处理结果通过报告信道传给基站,
(a).若本地判决为H1或H0,则直接将1比特判决结果即0或1发送给基站;
(b).若本地判决得到的是本地感知可信度αi,则将可信度进行M比特量化后再上传;其中,可信度M比特量化即 quantize(X,M)表示对X进行M比特量化;
步骤五:基站对本地判决结果xi按照长度进行分类并进行可信度融合,即基站将接收到的结果按照数据长度分为两类,即1比特数据b和M比特数据α,然后对所有M比特的可信度进行合并得到可信度融合结果ω,
本步骤中,可信度融合,具体是指按照软信息融合的方法和可信度判决门限进行融合。
步骤六:
基站将可信度融合判决结果ω与其他所有1比特信息b使用OR准则进行判决融合得到融合结果x0,即x0=γ(ω,b)。
本步骤中,所述的基站将可信度融合结果和接收到的本地判决结果统一进行融合判决,具体是:若多个未授权用户的感知结果中有任何一个宣布授权用户存在,则融合结果判决为授权用户存在。
实施例一:
在本实施例中,如图1所示,协同频谱感知系统中有一个授权用户、N个未授权用户和一个基站,在授权用户周围,未授权用户随机分布在基站周围的一定区域内,且未授权用户之间的距离远小于未授权用户与授权用户之间的距离。
本实施例的流程如图2所示,基于可信度的协同频谱感知方法包括如下步骤:
(a).未授权用户使用能量感知方法进行本地感知,即经过一段时间的累计得到能量值 其中y(i)是未授权用户接收到的信号,u是时间带宽积;
(b).未授权用户进行本地判决x=Φ(θ),判决规则如下:
其中H1代表授权用户存在,H0代表授权用户不存在, 是在本地感知结果θ落在本地判决门限之间区域时对信号存在的可信度。
(c).未授权用户将本地感知和判决的处理结果通过报告信道传给基站,若本地判决为H1或H0,则直接将1比特判决结果发送给基站;.若本地判决得到的是本地感知可信度,则将可信度进行M比特量化后再上传。可信度M比特量化,即 quantize(X,M)表示对X进行M比特量化。
(d).基站对收到未授权用户本地判决结果按照长度进行分类并进行可信度融合,基站根据接收到的数据长度将接收到的结果按照比特数分为两类,即1比特数据b和M比特数据α。然后对所有M比特的可信度进行合并得到可信度融合结果ω,
(e).基站将可信度融合判决结果ω与其他所有1比特信息b使用OR准则进行判决融合得到融合结果x0,即x0=γ(ω,b),
根据方法的步骤,下面在给定基于可信度的协同频谱感知方法的常虚警概率(CFAR)的情况下,以AWGN信道为例,评估基于可信度的协同频谱感知方法的感知性能(ROC曲线)。其中硬信息融合方法分别使用AND和OR两种融合准则。
基于可信度的协同频谱感知方法的常虚警概率设定为 检测时间为100个符号周期,未授权用户数目为N=20,噪声方差为σ2=1,计算本地判决门限时选择软信息融合方法的虚警概率 经初始化计算,软信息融合的判决门限为λ0=2199.4,本地判决门限为λ1=34.3,λ2=185.6。
图3给出了本地判决未量化情况下的感知性能(检测概率Qd)随信噪比SNR变化情况的仿真曲线。由图3可以看出,多个未授权用户的协同频谱感知方法的性能优于单个用户频谱感知方法,软信息融合方法优于硬信息融合方法,基于可信度的协同频谱感知方法的性能明显优于硬信息融合方法,且与软信息融合方法的感知性能非常接近。
图4给出了本地判决量化情况下的感知性能(检测概率Qd)随信噪比SNR变化情况的仿真曲线。由图4可以看出,当软信息融合方法的量化间隔为2时,其频谱感知的性能比较接近未量化的情况;基于可信度的协同频谱感知方法的量化比特数为3时,其频谱感知的性能非常接近未量化的情况。
图5给出了当软信息融合方法的量化间隔为2和基于可信度的协同频谱感知方法的量化比特数为3时,本地判决量化情况下的感知开销W随信噪比SNR变化情况的仿真曲线。由图5可以看到,在信噪比较低时,基于可信度的协同频谱感知方法的感知开销只有软信息融合方法的感知开销的一半左右;而当信噪比升高时,基于可信度的协同频谱感知方法的感知开销是逐渐减少的,而软信息融合的感知开销是逐渐增加的。
本发明未涉及部分均与现有技术相同或可采用现有技术加以实现。
Claims (1)
1.一种基于可信度的协同频谱感知方法,其特征是它包括以下步骤:
步骤一:未授权用户在本地使用基于能量检测的频谱感知方法进行感知得到感知结果 其中y是未授权用户的接收信号,θ是未授权用户的本地感知结果,u是时间带宽积;
步骤二:
(a).设定本方法的本地虚警概率Qf C,选择软信息融合方法的虚警概率Qf S,根据 计算软信息融合方法的判决门限λ0,其中 N是系统中参与协同频谱感知的未授权用户数,然后检验Qf S是否满足以下两个条件;
其中Qf OR(λ′)是硬信息融合方法的虚警概率, 是硬信息融合方法的本地判决门限;
若Qf S满足检验条件,转步骤(b);否则,转步骤(a)重新选择选择软信息融合方法的虚警概率Qf S;
(b).根据设定的本地虚警概率Qf C和软信息融合方法的判决门限λ0计算未授权用户本地感知判决门限的均值为和基站可信度融合门限其中L是参与可信度融合的未授权用户数;得到未授权用户的两个本地判决门限分别为 Δλ为两个本地判决门限λ1,λ2与均值之差,它满足:
其中 σ2为未授权用户的本地噪声方差;
步骤三:未授权用户进行本地判决x=Φ(θ),判决规则如下:
其中H1代表授权用户存在,H0代表授权用户不存在, 是在本地感知结果θ落在本地判决门限λ1、λ2之间区域时对信号存在的可信度;步骤四:未授权用户将本地感知和判决的处理结果通过报告信道传给基站,
(a).若本地判决为H1或H0,则直接将1比特判决结果即0或1发送给基站;
(b).若本地判决得到的是本地感知可信度αi,则将可信度进行M比特量化后再上传;其中,可信度M比特量化即 quantize(X,M)表示对X进行M比特量化;
步骤五:基站对收到的未授权用户本地判决结果xi按照长度进行分类并进行可信度融合,即基站将接收到的结果按照数据长度分为两类,即1比特数据b和M比特数据α,然后对所有M比特的可信度进行合并得到可信度融合结果ω,
其中θ′=λ2-(λ2-λ1)*α,∑为噪声协方差矩阵,即 σ1 2,σ2 2,...,σN 2为未授权用户的本地噪声方差,是可信度融合门限;
步骤六:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101097114A CN101815305B (zh) | 2010-02-10 | 2010-02-10 | 基于可信度的协同频谱感知方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101097114A CN101815305B (zh) | 2010-02-10 | 2010-02-10 | 基于可信度的协同频谱感知方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101815305A true CN101815305A (zh) | 2010-08-25 |
CN101815305B CN101815305B (zh) | 2013-07-17 |
Family
ID=42622392
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101097114A Active CN101815305B (zh) | 2010-02-10 | 2010-02-10 | 基于可信度的协同频谱感知方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101815305B (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102123428A (zh) * | 2011-02-28 | 2011-07-13 | 华为技术有限公司 | 空闲频谱检测中合作用户的选择方法及装置 |
CN102223194A (zh) * | 2011-07-21 | 2011-10-19 | 电子科技大学 | 一种频谱感知方法 |
CN102291191A (zh) * | 2011-08-19 | 2011-12-21 | 电子科技大学 | 一种频谱感知方法 |
WO2012151732A1 (zh) * | 2011-05-06 | 2012-11-15 | 中国人民解放军理工大学 | 基于状态转移概率估计的感知时隙长度优化方法 |
CN103873171A (zh) * | 2014-03-27 | 2014-06-18 | 西安电子科技大学 | 基于多用户互协作的协作频谱感知决策融合方法 |
CN104254078A (zh) * | 2013-06-28 | 2014-12-31 | 清华大学 | 主用户多电平环境下的认知无线电系统联合频谱感知方法 |
CN105307185A (zh) * | 2014-07-11 | 2016-02-03 | 中国人民解放军理工大学 | 一种基于数据净化的群智协同频谱感知方法 |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103746752B (zh) * | 2013-12-18 | 2016-05-18 | 同济大学 | 一种基于分层狄利克雷过程的智能频谱感知方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1838572A (zh) * | 2006-04-26 | 2006-09-27 | 电子科技大学 | 一种基于可靠度的分布式频谱检测方法 |
CN101459445A (zh) * | 2008-12-29 | 2009-06-17 | 浙江大学 | 一种认知无线电系统中的合作频谱感知方法 |
CN101511092A (zh) * | 2009-01-04 | 2009-08-19 | 中国人民解放军理工大学 | 基于二重可信的协同频谱感知方法 |
CN101521896A (zh) * | 2009-04-17 | 2009-09-02 | 东南大学 | 认知无线电中的基于似然比的协同频谱感知方法 |
-
2010
- 2010-02-10 CN CN2010101097114A patent/CN101815305B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1838572A (zh) * | 2006-04-26 | 2006-09-27 | 电子科技大学 | 一种基于可靠度的分布式频谱检测方法 |
CN101459445A (zh) * | 2008-12-29 | 2009-06-17 | 浙江大学 | 一种认知无线电系统中的合作频谱感知方法 |
CN101511092A (zh) * | 2009-01-04 | 2009-08-19 | 中国人民解放军理工大学 | 基于二重可信的协同频谱感知方法 |
CN101521896A (zh) * | 2009-04-17 | 2009-09-02 | 东南大学 | 认知无线电中的基于似然比的协同频谱感知方法 |
Non-Patent Citations (1)
Title |
---|
《Personal, Indoor and Mobile Radio Communications, 2006 IEEE 17th International Symposium》 20060914 [Peng Qihang,Zeng Kun,Wang Jun,Li Shaoqian A Distributed Spectrum Sensing Scheme based on Credibility and Evidence Theory in CognitiveRadio Contex , * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102123428A (zh) * | 2011-02-28 | 2011-07-13 | 华为技术有限公司 | 空闲频谱检测中合作用户的选择方法及装置 |
CN102123428B (zh) * | 2011-02-28 | 2014-06-04 | 华为技术有限公司 | 空闲频谱检测中合作用户的选择方法及装置 |
WO2012151732A1 (zh) * | 2011-05-06 | 2012-11-15 | 中国人民解放军理工大学 | 基于状态转移概率估计的感知时隙长度优化方法 |
CN102223194A (zh) * | 2011-07-21 | 2011-10-19 | 电子科技大学 | 一种频谱感知方法 |
CN102223194B (zh) * | 2011-07-21 | 2014-04-09 | 电子科技大学 | 一种频谱感知方法 |
CN102291191A (zh) * | 2011-08-19 | 2011-12-21 | 电子科技大学 | 一种频谱感知方法 |
CN102291191B (zh) * | 2011-08-19 | 2013-08-14 | 电子科技大学 | 一种频谱感知方法 |
CN104254078A (zh) * | 2013-06-28 | 2014-12-31 | 清华大学 | 主用户多电平环境下的认知无线电系统联合频谱感知方法 |
CN104254078B (zh) * | 2013-06-28 | 2017-11-17 | 清华大学 | 主用户多电平环境下的认知无线电系统联合频谱感知方法 |
CN103873171A (zh) * | 2014-03-27 | 2014-06-18 | 西安电子科技大学 | 基于多用户互协作的协作频谱感知决策融合方法 |
CN105307185A (zh) * | 2014-07-11 | 2016-02-03 | 中国人民解放军理工大学 | 一种基于数据净化的群智协同频谱感知方法 |
CN105307185B (zh) * | 2014-07-11 | 2018-09-04 | 中国人民解放军理工大学 | 一种基于数据净化的群智协同频谱感知方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101815305B (zh) | 2013-07-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101815305B (zh) | 基于可信度的协同频谱感知方法 | |
CN103338082B (zh) | 一种基于“k秩”准则的双门限协作频谱感知方法 | |
CN103997743B (zh) | 一种认知无线电系统中基于有效容量的资源分配方法 | |
CN103338458B (zh) | 一种用于认知无线电系统的协作频谱感知方法 | |
CN106656374A (zh) | 一种基于双门限能量检测的协作宽带频谱感知方法 | |
CN103929255B (zh) | 一种基于多信道的认知用户能效优化方法 | |
CN104038296B (zh) | 一种认知无线电网络的协作频谱检测方法 | |
CN104105193B (zh) | 一种异构网络中基于Starckelberg博弈的功率分配方法 | |
CN109391972A (zh) | 支持载波聚合的终端到终端通信的资源选择方法及装置 | |
CN101753232B (zh) | 协作频谱检测方法和系统 | |
CN106059840A (zh) | 一种认知无线系统功率分配方法及装置 | |
CN101742539B (zh) | 一种无线资源配置方法和系统 | |
CN102291191B (zh) | 一种频谱感知方法 | |
CN103139860B (zh) | 一种异构认知无线电网络中基于速度自适应的联合切换方法 | |
CN102271022B (zh) | 一种基于最大广义特征值的频谱感知方法 | |
CN103338085B (zh) | 基于双门限值的两层协作频谱感知方法 | |
CN104796991B (zh) | 基于势博弈的ofdma系统的资源分配方法 | |
CN104270212A (zh) | 一种基于分组数据型序贯能量检测的信道频谱感知方法 | |
CN103780318A (zh) | 一种动态双门限协作频谱感知方法 | |
CN102547795B (zh) | 自适应频谱感知调度方法 | |
CN102291188B (zh) | 一种认知无线电系统中的合作频谱感知方法 | |
CN103024886B (zh) | 宽带认知无线电中感知时间和子载波功率联合分配的方法 | |
CN104022839B (zh) | 一种应用于合作频谱感知的信息融合判决方法 | |
CN105634624A (zh) | 一种面向认知无线网络的频域协作频谱感知方法 | |
CN103856946A (zh) | 基于差分能量检测的双门限协作频谱感知方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |