CN104333423B - 一种基于基尼相关的信号检测电路及方法 - Google Patents

一种基于基尼相关的信号检测电路及方法 Download PDF

Info

Publication number
CN104333423B
CN104333423B CN201410455266.5A CN201410455266A CN104333423B CN 104333423 B CN104333423 B CN 104333423B CN 201410455266 A CN201410455266 A CN 201410455266A CN 104333423 B CN104333423 B CN 104333423B
Authority
CN
China
Prior art keywords
signal
input
sgn
parallel
noise
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
CN201410455266.5A
Other languages
English (en)
Other versions
CN104333423A (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.)
Zhongshan Yimeijie Intelligent Technology Co ltd
Original Assignee
Guangdong University of Technology
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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201410455266.5A priority Critical patent/CN104333423B/zh
Publication of CN104333423A publication Critical patent/CN104333423A/zh
Application granted granted Critical
Publication of CN104333423B publication Critical patent/CN104333423B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明涉及一种基于基尼相关的信号检测电路及方法,该电路包括两个比较器阵列、一个减法器阵列、两个乘法器阵列,两个加法器树,一个除法器以及一个寄存器;比较器、减法器阵列分别为一个n × n的矩阵(n为信号长度),矩阵上的每一个元素分别作比较运算和减法运算;乘法器阵列为n2个2输入的乘法器;加法器树为一个n2个输入的加法器;除法器完成2输入除法运算;寄存器寄存相关运算的结果。当环境噪声中含有脉冲噪声成分时,匹配滤波器和积矩相关系数基本失效,而基尼相关在含有脉冲成分的噪声环境中表现出其极佳的健壮性,包括十分接近真实值的数学期望以及较小的标准差。

Description

一种基于基尼相关的信号检测电路及方法
技术领域
本发明涉及通信技术领域,特别涉及一种用于雷达和声呐测距及无线电通信的同步捕获中基于基尼相关的信号检测电路及方法。
背景技术
捕获雷达和声呐的回波信号进行延时检测,以及在无线电通信中信号同步都需要对接收信号进行信号检测,判断特定的信号是否已经被接收机所接收。因此,对接收信号的信号检测是一个十分重要的问题。目前,信号延时检测的手段以匹配滤波器及以其为基础的改进型滤波器为主。匹配滤波器是一种在输出信号信噪比最大化意义下的最优线性滤波器。由于其形式简单,易于实现,且理论证明完备,得到了广泛的应用。但是匹配滤波器是在假设环境噪声为独立的加性高斯白噪声的前提下进行论证和设计的。而在真实环境中,噪声中经常会出现加性高斯白噪声以外的成分。其中,脉冲噪声是最为典型的环境噪声成分之一,可以由云层放电,冰层破裂,生物的活动等原因产生。由于接收到的信号中混入了脉冲噪声,对应的一段时间内信噪比会迅速下降。在这一情况之下,匹配滤波器的性能也会相应地下降,甚至失效。检测接收到的信号中是否存在感兴趣的信号片段,本质上是能否在较高的概率下成功地判断接收到的信号中是否包含着与已知的信号相似度高的信号片段的问题。因此,我们可以使用相关分析中的相关系数来处理这一类问题。文献中用于相关分析的经典方法是使用积矩相关系数(Pearson’s Product Moment CorrelationCoefficient,PPMCC)。而相关文献已经证明,PPMCC对于信号中的脉冲噪声极其敏感无法在脉冲噪声环境中很好地工作。
为减少脉冲噪声对信号的干扰而产生的估计误差,有效地对信号进行检测,本发明首次采用基尼相关(Gini Correlation)来对雷达和声呐信号以及无线电通信信号进行处理。令表示n对原始数据,并由此得到一组新的数据对序列其中X(1)<...<X(n)为X的序统计量,相应的Y[i]则为X(i)的伴随。假设Xj位于序列中的第K个位置,则定义数字K为Xj的秩次,记为Pj。类似地可以定义Yj的秩次并记为Qj。令分别代表Xi和Yi的算术平均。则基尼相关的定义如下:
其中Y表示未受噪声干扰的信道数据,而X则为受噪声污染的信道数据。由上式可见,基尼相关只利用数据对中其中一路信号的具体数据值,而另一路信号则使用其秩次信息。因为无论是雷达和声呐测距或无线电通信的同步捕获,都是使用一路已知的信号作为模板与一路接收到的待检测信号做相关分析并做出判决。所以已知的信号模板可以认为是完全“干净”的数据,应该得到充分的利用;而接收到的数据因为包含有噪声,所以要尽可能地减少噪声的影响。基尼相关的定义很好地适应了这一要求,既充分利用了有效的数据值,也利用了数据的秩次信息减少了噪声的干扰。在本发明中,我们利用了基尼相关对数值与秩次组合的数据信息的有效处理能力在尽可能减小噪声干扰的同时有效地对信号延时进行检测。
发明内容
本发明的主要目的在于,针对上述现有技术中的不足,提供一种基于基尼相关的信号检测电路及方法。
本发明解决现有技术问题所采用的技术方案是:一种基于基尼相关的信 号检测电路,包括两个比较器阵列、一个减法器阵列、两个乘法器阵列,两个加法器树,一个除法器以及一个寄存器;
信号x1……xn、y1……yn同时分别送入相应的比较器、减法器阵列,其中信号x1……xn并行完成n2次比较操作,即sgn(xi-xj),其中(i,j=1……n),而信号y1……yn并行完成n2次比较操作和减法操作,即sgn(yi-yj)和yi-yj,其中(i,j=1……n),以上过程同步进行,n为信号长度;
得到的比较,减法结果同时送入相应的乘法器阵列,分别同步完成n2次并行2输入乘法操作,即sgn(xi-xj)(yi-yj)和sgn(yi-yj)(yi-yj),其中(i,j=1……n);
得到的乘法结果同时送入相应的加法器树,分别同步完成n2个输入的并行加法操作,即其中(i,j=1……n);
得到的加法结果送入2输入除法器中,完成其中(i,j=1……n);
得到的相关运算结果送入寄存器中。
进一步的,比较器、减法器阵列为一个n×n的矩阵(n为信号长度),矩阵上的每一个元素分别作比较运算和减法运算;乘法器阵列为n2个2输入的乘法器;加法器树为一个n2个输入的加法器;除法器完成2输入除法运算;寄存器寄存相关运算的结果。
一种基于基尼相关的信号检测方法,包括步骤:
S1、信号x1……xn、y1……yn同时分别送入相应的比较器、减法器阵列,其中信号x1……xn并行完成n2次比较操作,即sgn(xi-xj),其中(i,j=1……n),而信号y1……yn并行完成n2次比较操作和减法操作,即sgn(yi-yj)和yi-yj,其中(i,j=1……n),以上过程同步进行,n为信号长度;
S2、步骤S1得到的比较,减法结果同时送入相应的乘法器阵列,分别同步完成n2次并行2输入乘法操作,即sgn(xi-xj)(yi-yj)和sgn(yi-yj)(yi-yj),其中, (i,j=1……n);
S3、步骤S2得到的乘法结果同时送入相应的加法器树,分别同步完成n2个输入的并行加法操作,即其中(i,j=1……n);
S4、步骤S3得到的加法结果送入2输入除法器中,完成其中(i,j=1……n);
S5、步骤S4得到的相关运算结果送入寄存器中。
本发明的有益效果是:当环境噪声中含有脉冲噪声成分时,匹配滤波器和积矩相关系数基本失效,而基尼相关在含有脉冲成分的噪声环境中表现出其极佳的健壮性,包括十分接近真实值的数学期望以及较小的标准差。这说明了基尼相关是一种在含有脉冲噪声干扰的环境中进行信号检测的有效工具,而且与目前普遍使用的方法相比具有极佳的性能优势。
附图说明
图1为本发明检测电路结构图;
图2为基尼相关进行信号检测的示意图;
本发明目的的实现、功能特点及优点将结合实施例,参照附图做进一步说明。
具体实施方式
以下将结合附图及具体实施例详细说明本发明的技术方案,以便更清楚、直观地理解本发明的发明实质。
信号检测是信号处理的一项基本内容,回波信号延时检测和通信信号同步捕获又是雷达、声呐和现代通信技术所必不可少的环节。本发明是针对现有常用信号延时检测技术无法很好地克服信号由于在传播中被加入的脉冲噪 声带来的干扰影响信号检测的判决而提出的。匹配滤波器可以使通过的信号得到最大的信噪比输出SNR=2Es/N0(其中Es为输入有用信号能量,N0/2为加性高斯白噪声的双边功率谱密度),并在此标准下作为最优的线性滤波器而在信号延时检测中得到最为广泛的应用。但对于其性能的理论论证过程是假设环境噪声是单纯的加性高斯白噪声,而实际的环境噪声中经常会包含有脉冲噪声分量。此时,一定时间内信噪比会大幅下降,进而严重影响匹配滤波器的性能。从相关分析的角度来看,经典的相关系数也存在对脉冲噪声敏感和因不能充分利用有效信息而带来的估计误差的问题。基尼相关可以通过利用被噪声污染的信号的秩次信息来减少噪声的干扰,同时充分利用干净信号的数据值信息避免有效信息的丢失减少估计误差。因此把基尼相关引入到信号处理领域作为一种信号检测的手段可以有效地解决脉冲噪声背景下的信号延时检测问题。并且通过推导得到基尼相关的一种等价表达形式,以此表达式为基础提出一种FPGA计算电路的结构设计。
具体的,定义函数H(t),当t>0时,H(t)=1;当t≤0时,H(t)=0。则对(1)式中的P和Q分别有:
符号函数Sgn(t),当t>0,Sgn(t)=1;当t=0,Sgn(t)=0;当t<0,Sgn(t)=-1。则:
其中Ti≠Tj′
表示(1)式中rG的分子,把(2)式和(4)式的关系代入得:
因为当i=j时Sgn(Xi-Xj)=0,所以:
交换下标i和j得:
由(5)式及(6)式得:
同理,令表示(1)式中rG的分母,把(2)式和(4)式的关系代入并作类似于(6)式的变换可得:
则有:
根据式(9)式可以设计出相应的本发明的电路。
参照图1所示,本发明提供了一种基于基尼相关的信号检测电路,包括两个比较器阵列、一个减法器阵列、两个乘法器阵列,两个加法器树,一个除法器以及一个寄存器;比较器、减法器阵列分别为一个n×n的矩阵(n为信号长度),矩阵上的每一个元素分别作比较运算和减法运算;乘法器阵列为n2个2输入的乘法器;加法器树为一个n2个输入的加法器;除法器完成2输入除法运算;寄存器寄存相关运算的结果。
具体的,信号x1……xn、y1……yn同时分别送入相应的比较器、减法器阵列,其中信号x1……xn并行完成n2次比较操作,即sgn(xi-xj),其中(i,j=1……n),而信号y1……yn并行完成n2次比较操作和减法操作,即sgn(yi-yj)和yi-yj,其中(i,j=1……n),以上过程同步进行,n为信号长度;
得到的比较,减法结果同时送入相应的乘法器阵列,分别同步完成n2次并行2输入乘法操作,即sgn(xi-xj)(yi-yj)和sgn(yi-yj)(yi-yj),其中(i,j=1……n);
得到的乘法结果同时送入相应的加法器树,分别同步完成n2个2输入的并行加法操作,即其中(i,j=1……n);
得到的加法结果送入2输入除法器中,完成其中(i,j=1……n);
得到的相关运算结果送入寄存器中。
一种基于基尼相关的信号检测方法,包括步骤:
S1、信号x1……xn、y1……yn同时分别送入相应的比较器、减法器阵列,其中信号x1……xn并行完成n2次比较操作,即sgn(xi-xj),其中(i,j=1……n),而信号y1……yn并行完成n2次比较操作和减法操作,即sgn(yi-Yj)和yi-yj,其中 (i,j=1……n),以上过程同步进行,n为信号长度;
S2、步骤S1得到的比较,减法结果同时送入相应的乘法器阵列,分别同步完成n2次并行2输入乘法操作,即sgn(xi-xj)(yi-yj)和sgn(yi-yj)(yi-yj),其中(i,j=1……n);
S3、步骤S2得到的乘法结果同时送入相应的加法器树,分别同步完成n2个输入的并行加法操作,即其中(i,j=1……n);
S4、步骤S3得到的加法结果送入2输入除法器中,完成其中(i,j=1……n);
S5、步骤S4得到的相关运算结果送入寄存器中。
本发明首次把基尼相关引入到信号处理领域,用以解决目前常用信号检测技术无法很好地处理环境噪声中的脉冲噪声分量对双信道相关分析的影响的问题。现有信号检测技术中,以匹配滤波器最为常用。但是其理论基础基于环境噪声为单纯的加性高斯白噪声的假设与现实应中的系统工作环境不完全相符。其中脉冲噪声分量是环境噪声中常见的成分。由于脉冲噪声的瞬时能量极大,在一定时间段里,大大降低了信号的信噪比,所以即使理想的匹配滤波器可以输出系统所能达到的最大信噪比也无法避免有用信号湮没在噪声中。因此在考虑脉冲噪声的情况下,以信噪比作为标准已然不合适。
相关系数作为相关分析中的一个强有力的工具也可以作为信号检测中评价两路信号相似程度的一个度量。文献中经典的相关系数是积矩相关系数。但是当这一相关系数被应用在一路信号是干净的信号模板而另一路是被噪声污染的信号时存在不足之处。积矩相关系数因为对脉冲噪声十分敏感而不能作为一种有效的度量。为了减少噪声,尤其是脉冲噪声的干扰,文献中指出可以使用相应数据的秩次信息而不直接使用数据值来应对。而已知的有用信息应该全部都得到应用,避免信息丢失引起误差。由此,基尼相关的定义本身就体现着这样的特点。
通过蒙特卡洛实验对比匹配滤波器、积矩相关系数、和基尼相关可以说明基尼相关在处理混杂着脉冲噪声的噪声背景下的工作性能。其示意图如图2所示,图2中(a)表示信号发射器发射的信号,也就是已知的信号模板。发出的信号进过传输由接收器接收,在理想状态之下,接收到的信号应是与发射信号波形完全一致但存在一定延时K0的信号如图2中(b)所示;但是由于环境中存在着各种噪声,主要是高斯白噪声和经常出现的脉冲噪声,接收器接收到的实际信号为噪声与信号的叠加如图2中(c)所示。为了检测出接收信号的延时,使用已知的信号模板作为扫描窗口如图2中(d)所示,按顺序从接收到的信号中每次选取一段等长度的信号与信号模板做相关分析得到相关值,最后判定所有相关值中的最大值所对应的时间点K0为想要检测的信号所在的位置如图2中(e)所示。
实验设计如下:发射信号(信号模板)Y设为:
其均值和方差分别为E(Y)=0和Var(Yr)=1。噪声部分由混合高斯模型产生:
(1-ε)N(0,σ2)+εN(0,σ′2),
其中ε=0.02表示脉冲噪声在整个噪声背景中发生的概率,σ′=104>>σ表示脉冲噪声。接收信号为有用信号与噪声相加之和。在此设定值下接收信号的信噪比定义为
蒙特卡洛实验次数设为10000次。设定回波信号位于接收信号的第1000个时间点即有k0=1000。实验目的为验证各种方法能否在不同的信噪比下 检测出隐藏在背景噪声中的有用信号。实验的示意图如图1所示。相关实验数据如表1所示。
SNR 匹配滤波器 积矩相关系数 基尼相关
0 748.20±461.72 753.62±455.13 999.98±0.89
-1 743.86±457.64 741.69±454.23 999.99±0.96
-2 757.11±457.24 749.62±452.87 999.99±1.06
-3 750.58±461.76 7S3.23±455.54 999.99±1.17
-4 756.80±456.89 749.51±452.75 999.97±1.63
-5 756.97±457.10 749.47±453.04 999.96±2.01
-6 750.91±461.86 753.33±455.27 1000.03±4.19
-7 756.89±456.90 749.66±452.73 999.87±6.35
-8 756.93±456.86 749.70±452.73 999.90±10.76
-9 751.90±462.03 753.08±455.56 999.99±15.72
-10 756.84±456.81 749.48±452.80 999.98±22.34
表1对比试验数据
表1中的数据显示方式是多次蒙特卡洛试验的“均值±标准差”。实验结果表明,当环境噪声中含有脉冲噪声成分时,匹配滤波器和积矩相关系数基本失效,而基尼相关在含有脉冲成分的噪声环境中表现出其极佳的健壮性,包括十分接近真实值的均值以及较小的标准差。这说明了基尼相关是一种在含有脉冲噪声干扰的环境中进行信号检测的有效工具,而且与目前普遍使用的方法相比具有极佳的性能优势。
以上所述仅为本发明的优选实施例,并非因此限制其专利范围,凡是利用本发明说明书及附图内容所作的等效结构或等效流程变换,直接或间接运用在其他相关的技术领域,均同理包括在本发明的专利保护范围内。

Claims (3)

1.一种基于基尼相关的信号检测电路,其特征在于,包括两个比较器阵列、一个减法器阵列、两个乘法器阵列,两个加法器树,一个除法器以及一个寄存器;
信号x1……xn、y1……yn同时分别送入相应的比较器、减法器阵列,其中信号x1……xn并行完成n2次比较操作,即sgn(xi−xj) ,而信号y1……yn并行完成n2次比较操作和减法操作,即sgn(yi−yj)和yi−yj,以上过程同步进行,n为信号长度;
得到的比较,减法结果同时送入相应的乘法器阵列,分别同步完成n2次并行2输入乘法操作,即sgn(xi−xj)(yi−yj)和sgn(yi−yj)(yi−yj),其中(i,j=1……n);
得到的乘法结果同时送入相应的加法器树,分别同步完成n2个输入的并行加法操作,即,其中(i,j=1……n);
得到的加法结果送入2输入除法器中,完成
得到的相关运算结果送入寄存器中。
2.根据权利要求1所述的基于基尼相关的信号检测电路,其特征在于:比较器、减法器阵列为分别为一个n × n的矩阵(n为信号长度),矩阵上的每一个元素分别作比较运算和减法运算;乘法器阵列为n2个2输入的乘法器;加法器树为一个n2个输入的加法器;除法器完成2输入除法运算;寄存器寄存相关运算的结果。
3.一种基于基尼相关的信号检测方法,其特征在于,包括步骤:
S1、信号x1……xn、y1……yn同时分别送入相应的比较器、减法器阵列,其中信号x1……xn并行完成n2次比较操作,即sgn(xi−xj),而信号y1……yn并行完成n2次比较操作和减法操作,即sgn(yi−yj)和yi−yj,以上过程同步进行,n为信号长度;
S2、 步骤S1得到的比较,减法结果同时送入相应的乘法器阵列,分别同步完成n2次并行2输入乘法操作,即sgn(xi−xj)(yi−yj)和sgn(yi−yj)(yi−yj) ;
S3、步骤S2得到的乘法结果同时送入相应的加法器树,分别同步完成n2个输入的并行加法操作,即
S4、步骤S3得到的加法结果送入2输入除法器中,完成
S5、步骤S4得到的相关运算结果送入寄存器中。
CN201410455266.5A 2014-09-09 2014-09-09 一种基于基尼相关的信号检测电路及方法 Expired - Fee Related CN104333423B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410455266.5A CN104333423B (zh) 2014-09-09 2014-09-09 一种基于基尼相关的信号检测电路及方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410455266.5A CN104333423B (zh) 2014-09-09 2014-09-09 一种基于基尼相关的信号检测电路及方法

Publications (2)

Publication Number Publication Date
CN104333423A CN104333423A (zh) 2015-02-04
CN104333423B true CN104333423B (zh) 2017-05-10

Family

ID=52408085

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410455266.5A Expired - Fee Related CN104333423B (zh) 2014-09-09 2014-09-09 一种基于基尼相关的信号检测电路及方法

Country Status (1)

Country Link
CN (1) CN104333423B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106648901A (zh) * 2016-12-29 2017-05-10 广东工业大学 一种多通道信号相关性分析方法及系统
CN110083865B (zh) * 2019-03-25 2022-11-22 广东工业大学 一种抗脉冲噪声的已知信号检测方法及检测电路
CN110389745B (zh) * 2019-05-15 2023-10-13 广东工业大学 一种基于基尼回归的信号检测与估计电路及方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1254033C (zh) * 2002-12-20 2006-04-26 华为技术有限公司 一种多径搜索方法和装置
US7471717B2 (en) * 2004-05-17 2008-12-30 Yi Ping Huang Apparatus and method for acquiring spread-spectrum signals

Also Published As

Publication number Publication date
CN104333423A (zh) 2015-02-04

Similar Documents

Publication Publication Date Title
CN107220606B (zh) 基于一维卷积神经网络的雷达辐射源信号的识别方法
CN104392086B (zh) 一种基于皮尔逊秩次变量相关系数的信号检测电路及方法
CN107356921B (zh) 一种频率分集阵列雷达基于一次频偏目标定位的方法
CN110619296A (zh) 一种基于奇异分解的信号降噪方法
CN107884751A (zh) 一种利用单通道接收信号作信源数目估计的方法
CN104333423B (zh) 一种基于基尼相关的信号检测电路及方法
CN103235294A (zh) 一种基于外辐射源定位的微弱信号分离估计方法
CN104502905B (zh) 基于多通道顺序统计的距离扩展目标检测方法
CN112162244A (zh) 一种相关噪声和随机丢包环境下的事件触发目标跟踪方法
CN102833016B (zh) 一种在多径环境下针对同源信号时差的测定方法
CN107390196A (zh) 基于多基地雷达的快起伏目标双门限恒虚警检测方法
CN107121662A (zh) 基于空域稀疏表示的单站无源定位方法
CN108872961A (zh) 基于低门限的雷达微弱目标检测方法
CN115616503A (zh) 一种基于卷积神经网络模型的雷达干扰信号类型识别方法
CN107390167A (zh) 一种多通道接收机信号模拟方法
CN110954862B (zh) 一种稀疏贝叶斯框架下基于全局窄带模型的辐射源直接定位方法
CN110083865B (zh) 一种抗脉冲噪声的已知信号检测方法及检测电路
CN104200093A (zh) 一种积分延拓抑制局域均值分解端点效应的方法
CN102087313B (zh) 一种卫星搜救信号的频率估计方法
CN104408027A (zh) 一种基于广义协方差和张量分解的欠定盲辨识方法
CN106353743A (zh) 匹配于等效形状参数的近最优雷达目标检测方法
CN104363078B (zh) 基于鲁棒竞争聚类的欠定系统实正交空时分组码盲识别方法
CN105244040B (zh) 一种音频信号一致性对比方法
CN107607913A (zh) 基于对数累积量的海杂波Pareto分布参数估计方法
CN105487049B (zh) 一种检测识别非直达超宽带信号的方法及系统

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
CB03 Change of inventor or designer information

Inventor after: Xu Weichao

Inventor after: Ma Rubao

Inventor after: Lin Baohong

Inventor after: Zhang Yun

Inventor before: Ma Rubao

Inventor before: Xu Weichao

Inventor before: Lin Baohong

Inventor before: Zhang Yun

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20180726

Address after: 528400 room 1310, digital exhibition hall, 16 Convention and Exhibition Road, Zhongshan Torch Development Zone, Guangdong

Patentee after: Zhongshan Yimeijie Intelligent Technology Co.,Ltd.

Address before: 510006 100 Panyu West Road, Panyu University, Guangzhou, Guangdong

Patentee before: Guangdong University of Technology

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

Granted publication date: 20170510

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