CN107196721B - 针对时间异步且存在频偏情况下的ofdm信号频谱感知方法 - Google Patents

针对时间异步且存在频偏情况下的ofdm信号频谱感知方法 Download PDF

Info

Publication number
CN107196721B
CN107196721B CN201710472617.7A CN201710472617A CN107196721B CN 107196721 B CN107196721 B CN 107196721B CN 201710472617 A CN201710472617 A CN 201710472617A CN 107196721 B CN107196721 B CN 107196721B
Authority
CN
China
Prior art keywords
signal
sampling
ofdm
received signal
frequency offset
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
CN201710472617.7A
Other languages
English (en)
Other versions
CN107196721A (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.)
Ningbo University Science Park Development Co ltd
Original Assignee
Ningbo 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 Ningbo University filed Critical Ningbo University
Priority to CN201710472617.7A priority Critical patent/CN107196721B/zh
Publication of CN107196721A publication Critical patent/CN107196721A/zh
Application granted granted Critical
Publication of CN107196721B publication Critical patent/CN107196721B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04BTRANSMISSION
    • H04B17/00Monitoring; Testing
    • H04B17/30Monitoring; Testing of propagation channels
    • H04B17/382Monitoring; Testing of propagation channels for resource allocation, admission control or handover

Landscapes

  • Physics & Mathematics (AREA)
  • Electromagnetism (AREA)
  • Engineering & Computer Science (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Signal Processing (AREA)
  • Monitoring And Testing Of Transmission In General (AREA)

Abstract

本发明公开了一种针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其对来自监测信道的接收信号进行采样,得到采样信号;然后,根据采样信号中的所有采样点的采样值,并在时间异步且存在频偏情况下,利用接收信号中的OFDM信号的循环前缀的自相关性,估计得到噪声功率;再根据采样信号中的采样点的采样值,计算采样信号的协方差矩阵;接着,根据噪声功率和采样信号的协方差矩阵,计算检验统计量;最后,通过比较检验统计量与获取的判决门限的大小,判定在监测信道内有无授权用户信号,以实现频谱感知;优点是能够有效地提高OFDM信号的频谱感知性能,且计算复杂度低。

Description

针对时间异步且存在频偏情况下的OFDM信号频谱感知方法
技术领域
本发明涉及一种认知无线电***中的频谱感知方法,尤其是涉及一种基于噪声功率估计的针对时间异步且存在频偏情况下的OFDM信号频谱感知方法。
背景技术
随着各种无线通信业务以及新兴通信技术(比如ZigBee技术、快速接入模块、LTE技术、超宽带网络等)的快速发展,频谱资源的重要性和稀缺性日益凸显。大量的研究表明,当前的频谱资源短缺现状并不是由于频率资源物理上的不足,而主要是由于固定频谱分配管理机制导致频谱利用率低造成的。认知无线电正是试图通过提高频谱利用率,从本质上解决无线通信资源越来越紧缺的问题而提出的。认知无线电***能够实时感知周围的通信环境,识别出可利用的空闲信道,然后根据频谱感知结果自适应地调整认知无线网络的***参数,使得认知无线电***具有智能化辨识与改变频谱使用机会的能力。为了防止对已有通信***产生干扰,认知无线电***必须能够有效可靠地识别出空闲信道,因此频谱感知是认知无线电中的关键技术之一。
正交频分复用(Orthogonal frequency division multiplexing,OFDM)技术具有频谱利用率高的特点,并且是被当前和将来无线通信标准广泛采用的技术之一,因此对OFDM信号的频谱感知(即判断信道中是否存在OFDM信号)研究具有非常重要的意义。现有的针对OFDM信号的频谱感知方法可以分为频域检测和时域检测两类。其中,频域检测需要计算采样信号的频谱,因此具有较大的计算量;时域检测主要利用OFDM信号中循环前缀的自相关特性实现频谱感知,计算量较低。Zeng等人于2009年在《Spectrum sensingalgorithms for cognitive radio based on statistical covariance(认知无线电中基于统计协方差的频谱感知方法)》中提出的基于协方差矩阵频谱感知方法,其能够利用OFDM信号中的相关性来实现频谱感知,具有计算量低的特点,但是该方法没有考虑时间异步且存在频偏的具体情况,所以无法有效利用OFDM信号的循环前缀自相关特性。
发明内容
本发明所要解决的技术问题是提供一种基于噪声功率估计的针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其能够有效地提高OFDM信号的频谱感知性能,且计算复杂度低。
本发明解决上述技术问题所采用的技术方案为:一种针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其特征在于它的处理过程为:首先,对来自监测信道的接收信号进行采样,得到采样信号;然后,根据采样信号中的所有采样点的采样值,并在时间异步且存在频偏情况下,利用接收信号中的OFDM信号的循环前缀的自相关性,估计得到噪声功率;再根据采样信号中的采样点的采样值,计算采样信号的协方差矩阵;接着,根据噪声功率和采样信号的协方差矩阵,计算检验统计量;最后,通过比较检验统计量与获取的判决门限的大小,判定在监测信道内有无授权用户信号,以实现频谱感知。
该OFDM信号频谱感知方法具体包括以下步骤:
步骤一:利用认知无线电***中的采样模块对来自监测信道的接收信号进行M次采样,得到由M个采样点的采样值构成的采样信号,其中,M=K×N,K表示接收信号中的OFDM信号中包含的OFDM符号的总个数,K≥1,N表示接收信号中的OFDM信号中的任意一个OFDM符号内的子载波的总个数,N=Nc+Nd,Nc表示接收信号中的OFDM信号的循环前缀的长度,Nd表示接收信号中的OFDM信号中的任意一个OFDM符号内的有效数据的子载波的总个数,Nd≥Nc>1,且设Nc和Nd都为偶数;
步骤二:根据采样信号中的所有采样点的采样值,并在时间异步且存在频偏情况下,利用接收信号中的OFDM信号的循环前缀的自相关性,估计得到噪声功率,记为
Figure BDA0001327452390000021
Figure BDA0001327452390000022
其中,ρη表示接收信号的相关系数的估计值,
Figure BDA0001327452390000023
表示接收信号的功率的估计值;
步骤三:根据采样信号中的采样点的采样值,计算采样信号的协方差矩阵,记为Rx,Rx为一个L×L维的矩阵,将Rx中第i行第q列的元素记为Rx(i,q),
Figure BDA0001327452390000031
其中,L∈[2,M-1],1≤i≤L,1≤q≤L,x(p)表示采样信号中的第p个采样点的采样值,x*(p-q+i)为x(p-q+i)的共扼,x(p-q+i)表示采样信号中的第p-q+i个采样点的采样值;
步骤四:根据
Figure BDA0001327452390000032
和Rx,计算检验统计量,记为T,
Figure BDA0001327452390000033
步骤五:获取判决门限,记为λ;然后比较检验统计量T与判决门限λ的大小,如果T>λ,则判定在监测信道内有授权用户信号;如果T≤λ,则判定在监测信道内无授权用户信号,监测信道处于空闲状态,从而实现频谱感知。
所述的步骤二中,
Figure BDA0001327452390000034
其中,Re()为取实部操作函数,
Figure BDA0001327452390000035
Figure BDA0001327452390000036
x(2n-1)表示采样信号中的第2n-1个采样点的采样值,x*(2n-1+Nd)为x(2n-1+Nd)的共轭,x(2n-1+Nd)表示采样信号中的第2n-1+Nd个采样点的采样值,x(2n)表示采样信号中的第2n个采样点的采样值,x*(2n+Nd)为x(2n+Nd)的共轭,x(2n+Nd)表示采样信号中的第2n+Nd个采样点的采样值,j为虚数表示形式,符号“∠”为角度表示符号。
所述的步骤二中,
Figure BDA0001327452390000037
其中,1≤m≤M,x(m)表示采样信号中的第m个采样点的采样值,符号“| |”为取绝对值符号。
与现有技术相比,本发明的优点在于:
1)本发明方法不需要计算协方差矩阵的特征值和采样信号的频谱,因此其具有计算复杂度低,操作简单的特点。
2)本发明方法在存在频偏的情况下,将采样信号分成奇偶两部分来计算接收信号的相关系数,从而消除频偏对估计的噪声功率精确度的影响,提高了OFDM信号的频谱感知性能。
3)与传统的噪声功率估计器仅在噪声样本中有效相比,本发明方法估计噪声功率的过程充分利用了OFDM信号的循环前缀的自相关性,从而能够准确地估计出当前感知时隙内的噪声功率。
附图说明
图1为本发明方法的流程框图;
图2为当Nd=90、Nc=10且在接收信噪比为-10分贝时,分别利用本发明方法与Zeng等人提出的协方差法得到的ROC曲线的比较示意图;
图3为当Nd=80、Nc=20且在接收信噪比为-10分贝时,分别利用本发明方法与Zeng等人提出的协方差法得到的ROC曲线的比较示意图。
具体实施方式
以下结合附图实施例对本发明作进一步详细描述。
本发明提出的一种针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其流程框图如图1所示,其处理过程为:首先,对来自监测信道的接收信号进行采样,得到采样信号;然后,根据采样信号中的所有采样点的采样值,并在时间异步且存在频偏情况下,利用接收信号中的OFDM信号的循环前缀的自相关性,估计得到噪声功率;再根据采样信号中的采样点的采样值,计算采样信号的协方差矩阵;接着,根据噪声功率和采样信号的协方差矩阵,计算检验统计量;最后,通过比较检验统计量与获取的判决门限的大小,判定在监测信道内有无授权用户信号(即OFDM信号),以实现频谱感知。
本发明的针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其具体包括以下步骤:
步骤一:利用认知无线电***中的采样模块对来自监测信道的接收信号进行M次采样,得到由M个采样点的采样值构成的采样信号,其中,M=K×N,K表示接收信号中的OFDM信号中包含的OFDM符号的总个数,K≥1,N表示接收信号中的OFDM信号中的任意一个OFDM符号内的子载波的总个数,N=Nc+Nd,Nc表示接收信号中的OFDM信号的循环前缀的长度,Nd表示接收信号中的OFDM信号中的任意一个OFDM符号内的有效数据的子载波的总个数,Nd≥Nc>1,且设Nc和Nd都为偶数。
步骤二:根据采样信号中的所有采样点的采样值,并在时间异步且存在频偏情况下,利用接收信号中的OFDM信号的循环前缀的自相关性,估计得到噪声功率,记为
Figure BDA0001327452390000051
Figure BDA0001327452390000052
其中,ρη表示接收信号的相关系数的估计值,
Figure BDA0001327452390000053
表示接收信号的功率的估计值。
在此具体实施例中,步骤二中,
Figure BDA0001327452390000054
其中,Re()为取实部操作函数,
Figure BDA0001327452390000055
Figure BDA0001327452390000056
x(2n-1)表示采样信号中的第2n-1个采样点的采样值,x*(2n-1+Nd)为x(2n-1+Nd)的共轭,x(2n-1+Nd)表示采样信号中的第2n-1+Nd个采样点的采样值,x(2n)表示采样信号中的第2n个采样点的采样值,x*(2n+Nd)为x(2n+Nd)的共轭,x(2n+Nd)表示采样信号中的第2n+Nd个采样点的采样值,j为虚数表示形式,符号“∠”为角度表示符号。
在此具体实施例中,步骤二中,
Figure BDA0001327452390000057
其中,1≤m≤M,x(m)表示采样信号中的第m个采样点的采样值,符号“| |”为取绝对值符号。
步骤三:根据采样信号中的采样点的采样值,计算采样信号的协方差矩阵,记为Rx,Rx为一个L×L维的矩阵,将Rx中第i行第q列的元素记为Rx(i,q),
Figure BDA0001327452390000058
其中,L∈[2,M-1],如取
Figure BDA0001327452390000059
1≤i≤L,1≤q≤L,x(p)表示采样信号中的第p个采样点的采样值,x*(p-q+i)为x(p-q+i)的共扼,x(p-q+i)表示采样信号中的第p-q+i个采样点的采样值。
步骤四:根据
Figure BDA00013274523900000510
和Rx,计算检验统计量,记为T,
Figure BDA00013274523900000511
步骤五:利用现有技术获取判决门限,记为λ;然后比较检验统计量T与判决门限λ的大小,如果T>λ,则判定在监测信道内有授权用户信号;如果T≤λ,则判定在监测信道内无授权用户信号,监测信道处于空闲状态,从而实现频谱感知。
通过以下仿真以进一步说明本发明方法的可行性和有效性。
假设来自监测信道的接收信号中的OFDM信号时间异步且存在频偏,接收信号中的OFDM信号中的任意一个OFDM符号的有效数据内的子载波的总个数Nd=90,接收信号中的OFDM信号的循环前缀的长度Nc=10,接收信号中的OFDM信号中的任意一个OFDM符号内的子载波的总个数N=Nc+Nd=10+90=100,取K=14,则采样信号由M=K×N=14×100=1400个采样点的采样值构成。设取L=5。
图2给出了当Nd=90、Nc=10且在接收信噪比为-10分贝时,分别利用本发明方法与Zeng等人提出的协方差法得到的ROC曲线。从图2中可以看出,利用本发明方法得到的ROC曲线远高于利用Zeng等人提出的协方差法得到的ROC曲线,充分表明了本发明方法的频谱感知性能优于Zeng等人提出的协方差法。
假设来自监测信道的接收信号中的OFDM信号时间异步且存在频偏,接收信号中的OFDM信号中的任意一个OFDM符号的有效数据内的子载波的总个数Nd=80,接收信号中的OFDM信号的循环前缀的长度Nc=20,接收信号中的OFDM信号中的任意一个OFDM符号内的子载波的总个数N=Nc+Nd=20+80=100,取K=14,则采样信号由M=K×N=14×100=1400个采样点的采样值构成。设取L=5。
图3给出了当Nd=80、Nc=20且在接收信噪比为-10分贝时,分别利用本发明方法与Zeng等人提出的协方差法得到的ROC曲线。从图3中可以看出,利用本发明方法得到的ROC曲线远高于利用Zeng等人提出的协方差法得到的ROC曲线,充分表明了本发明方法的频谱感知性能优于Zeng等人提出的协方差法。
另外,图2和图3所示的ROC曲线称为感受性曲线,横坐标为虚警概率,纵坐标为检测概率;曲线下方的面积越大,则表示该曲线所对应的频谱感知方法的检测性能越好。

Claims (3)

1.一种针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其特征在于它的处理过程为:首先,对来自监测信道的接收信号进行采样,得到采样信号;然后,根据采样信号中的所有采样点的采样值,并在时间异步且存在频偏情况下,利用接收信号中的OFDM信号的循环前缀的自相关性,估计得到噪声功率;再根据采样信号中的采样点的采样值,计算采样信号的协方差矩阵;接着,根据噪声功率和采样信号的协方差矩阵,计算检验统计量;最后,通过比较检验统计量与获取的判决门限的大小,判定在监测信道内有无授权用户信号,以实现频谱感知;
该针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其特征在于它具体包括以下步骤:
步骤一:利用认知无线电***中的采样模块对来自监测信道的接收信号进行M次采样,得到由M个采样点的采样值构成的采样信号,其中,M=K×N,K表示接收信号中的OFDM信号中包含的OFDM符号的总个数,K≥1,N表示接收信号中的OFDM信号中的任意一个OFDM符号内的子载波的总个数,N=Nc+Nd,Nc表示接收信号中的OFDM信号的循环前缀的长度,Nd表示接收信号中的OFDM信号中的任意一个OFDM符号内的有效数据的子载波的总个数,Nd≥Nc>1,且设Nc和Nd都为偶数;
步骤二:根据采样信号中的所有采样点的采样值,并在时间异步且存在频偏情况下,利用接收信号中的OFDM信号的循环前缀的自相关性,估计得到噪声功率,记为
Figure FDA0002357748400000011
Figure FDA0002357748400000012
其中,ρη表示接收信号的相关系数的估计值,
Figure FDA0002357748400000013
表示接收信号的功率的估计值;
步骤三:根据采样信号中的采样点的采样值,计算采样信号的协方差矩阵,记为Rx,Rx为一个L×L维的矩阵,将Rx中第i行第q列的元素记为Rx(i,q),
Figure FDA0002357748400000014
其中,L∈[2,M-1],1≤i≤L,1≤q≤L,x(p)表示采样信号中的第p个采样点的采样值,x*(p-q+i)为x(p-q+i)的共扼,x(p-q+i)表示采样信号中的第p-q+i个采样点的采样值;
步骤四:根据
Figure FDA0002357748400000021
和Rx,计算检验统计量,记为T,
Figure FDA0002357748400000022
步骤五:获取判决门限,记为λ;然后比较检验统计量T与判决门限λ的大小,如果T>λ,则判定在监测信道内有授权用户信号;如果T≤λ,则判定在监测信道内无授权用户信号,监测信道处于空闲状态,从而实现频谱感知。
2.根据权利要求1所述的针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其特征在于所述的步骤二中,
Figure FDA0002357748400000023
其中,Re()为取实部操作函数,
Figure FDA0002357748400000024
Figure FDA0002357748400000025
x(2n-1)表示采样信号中的第2n-1个采样点的采样值,x*(2n-1+Nd)为x(2n-1+Nd)的共轭,x(2n-1+Nd)表示采样信号中的第2n-1+Nd个采样点的采样值,x(2n)表示采样信号中的第2n个采样点的采样值,x*(2n+Nd)为x(2n+Nd)的共轭,x(2n+Nd)表示采样信号中的第2n+Nd个采样点的采样值,j为虚数表示形式,符号“∠”为角度表示符号。
3.根据权利要求1或2所述的针对时间异步且存在频偏情况下的OFDM信号频谱感知方法,其特征在于所述的步骤二中,
Figure FDA0002357748400000026
其中,1≤m≤M,x(m)表示采样信号中的第m个采样点的采样值,符号“| |”为取绝对值符号。
CN201710472617.7A 2017-06-21 2017-06-21 针对时间异步且存在频偏情况下的ofdm信号频谱感知方法 Active CN107196721B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710472617.7A CN107196721B (zh) 2017-06-21 2017-06-21 针对时间异步且存在频偏情况下的ofdm信号频谱感知方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710472617.7A CN107196721B (zh) 2017-06-21 2017-06-21 针对时间异步且存在频偏情况下的ofdm信号频谱感知方法

Publications (2)

Publication Number Publication Date
CN107196721A CN107196721A (zh) 2017-09-22
CN107196721B true CN107196721B (zh) 2020-07-28

Family

ID=59879302

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710472617.7A Active CN107196721B (zh) 2017-06-21 2017-06-21 针对时间异步且存在频偏情况下的ofdm信号频谱感知方法

Country Status (1)

Country Link
CN (1) CN107196721B (zh)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101245460B1 (ko) * 2011-05-04 2013-03-19 성균관대학교산학협력단 Ofdm신호 검파를 위한 사이클로스테이셔너리티 기반 스펙트럼 센싱 기술의 성능 분석 방법
CN103248443A (zh) * 2013-05-06 2013-08-14 宁波大学 时间异步且循环前缀长度已知情况下的ofdm频谱感知方法
CN105281854A (zh) * 2015-11-05 2016-01-27 西安电子科技大学 一种基于非圆信号的局部最大功效不变检验频谱感知方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR101245460B1 (ko) * 2011-05-04 2013-03-19 성균관대학교산학협력단 Ofdm신호 검파를 위한 사이클로스테이셔너리티 기반 스펙트럼 센싱 기술의 성능 분석 방법
CN103248443A (zh) * 2013-05-06 2013-08-14 宁波大学 时间异步且循环前缀长度已知情况下的ofdm频谱感知方法
CN105281854A (zh) * 2015-11-05 2016-01-27 西安电子科技大学 一种基于非圆信号的局部最大功效不变检验频谱感知方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
《Development and performance evaluation of blind eigenvalue to noise power ratio detector for co-operative spectrum sensing》;Abhijeet A. Masal等;《2014 International Conference on Electronics and Communication Systems (ICECS)》;20140908;第1-5页 *
《Spectrum Sensing Using Weighted Covariance Matrix in Rayleigh Fading Channels》;Ming Jin等;《IEEE Transactions on Vehicular Technology》;20151130;第64卷(第11期);第5137-5148页 *
《基于OFDM信号导频周期性的频谱感知方法及USRP实现》;姚俊腾等;《无线通信技术》;20170516;第26卷(第1期);第1-4页 *

Also Published As

Publication number Publication date
CN107196721A (zh) 2017-09-22

Similar Documents

Publication Publication Date Title
CN106130942B (zh) 一种基于循环谱的无线通信信号调制识别及参数估计方法
WO2006096728A2 (en) System and method for ranging
CN108322277B (zh) 一种基于协方差矩阵反特征值的频谱感知方法
CN101557378A (zh) Ofdm***中导频发送、信道估计和噪声功率估计方法
CN104270234B (zh) 一种基于欠奈奎斯特采样的宽带信号检测识别方法
CN109698724A (zh) 入侵检测方法、装置、设备及存储介质
CN108989259B (zh) 无线综测仪窄带物理上行共享信道的时偏估计方法及***
CN106357574A (zh) 基于顺序统计量的bpsk/qpsk信号调制盲识别方法
CN101635695B (zh) 一种确定最大多径时延的方法及装置
CN109738868B (zh) 一种基于信道识别的外辐射源雷达非平稳杂波抑制方法
CN107196720B (zh) 针对时间同步且不存在频偏情况下的ofdm信号频谱感知方法
CN107196721B (zh) 针对时间异步且存在频偏情况下的ofdm信号频谱感知方法
WO2010080995A1 (en) Methods and systems for time tracking in ofdm systems
CN103346849A (zh) 一种抵抗模仿主用户信号的恶意攻击的频谱感知方法
CN116132237B (zh) 一种基于周期卷积与分形的ofdm信号检测识别方法
CN102238123A (zh) 测距码检测方法和装置
CN107465473B (zh) 针对时间同步且存在频偏情况下的ofdm信号频谱感知方法
CN107222870B (zh) 针对时间异步且不存在频偏情况下的ofdm信号频谱感知方法
CN103248442B (zh) 一种时间同步情况下的ofdm信号频谱感知方法
Algriree et al. A CR-5G network based on multi-user for various waveforms detection
US20110007791A1 (en) Method for the detection and generation of a useful signal and associated devices and communications system
CN103248443B (zh) 时间异步且循环前缀长度已知情况下的ofdm频谱感知方法
Wan et al. The modified iterative detector/estimator algorithm for sparse channel estimation
CN110691055A (zh) Oqam/ofdm中的时频偏联合估计方法
CN104936190B (zh) 基于稀疏系数信息估计的自适应宽带频谱压缩感知方法

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

Effective date of registration: 20221223

Address after: No. 777, Zhongguan West Road, Zhuangshi Street, Ningbo City, Zhejiang Province 315000

Patentee after: Ningbo University Science Park Development Co.,Ltd.

Address before: 315211, Fenghua Road, Jiangbei District, Zhejiang, Ningbo 818

Patentee before: Ningbo University

TR01 Transfer of patent right