CN109120316B - 一种基于最大化最小相位差值的混合波束成形方法 - Google Patents

一种基于最大化最小相位差值的混合波束成形方法 Download PDF

Info

Publication number
CN109120316B
CN109120316B CN201811020499.7A CN201811020499A CN109120316B CN 109120316 B CN109120316 B CN 109120316B CN 201811020499 A CN201811020499 A CN 201811020499A CN 109120316 B CN109120316 B CN 109120316B
Authority
CN
China
Prior art keywords
vector
beam forming
user
phase difference
codebook
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
CN201811020499.7A
Other languages
English (en)
Other versions
CN109120316A (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.)
Dongying Zhengnuo Technology Service Co ltd
Original Assignee
Hangzhou Dianzi 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 Hangzhou Dianzi University filed Critical Hangzhou Dianzi University
Priority to CN201811020499.7A priority Critical patent/CN109120316B/zh
Publication of CN109120316A publication Critical patent/CN109120316A/zh
Application granted granted Critical
Publication of CN109120316B publication Critical patent/CN109120316B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/0413MIMO systems
    • 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/0408Diversity systems; Multi-antenna system, i.e. transmission or reception using multiple antennas using two or more spaced independent antennas using two or more beams, i.e. beam diversity
    • 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
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W16/00Network planning, e.g. coverage or traffic planning tools; Network deployment, e.g. resource partitioning or cells structures
    • H04W16/24Cell structures
    • H04W16/28Cell structures using beam steering

Abstract

本发明涉及一种基于最大化最小相位差值的混合波束成形方法。现有方法模拟波束成形部分需要大量码本集合用于恢复终端发送的信号,导致复杂度过高,难以实现。本发明首先利用信道增益门限和阵列响应向量的相关特性,从模拟波束成形向量码本中筛选出能够提供较高信号功率的向量集合。然后利用模拟波束成形向量与主要传播路径阵列响应向量的相位差衡量用户间干扰,提出采用最大化最小相位差值的准则抑制最强用户间干扰。本发明利用信道增益门限从模拟波束成形向量码本中筛选出能够提供较高信号功率的向量集合,保证了较大的接收功率。本发明方法能够减轻用户间的干扰,提高系统容量,同时保证较低复杂度。

Description

一种基于最大化最小相位差值的混合波束成形方法
技术领域
本发明属于无线通信技术领域,涉及一种对多用户系统降低干扰的方法,具体是毫米波大规模多输入多输出(multiple-input multiple-output,MIMO)多用户系统中的一种基于最大化最小相位差值的混合波束成形方法。
背景技术
随着无线通信的迅速发展,现有的低频段频谱不能满足人们日益增长的需求。使用未被充分利用的毫米波频段可满足人们高速通信的需求。毫米波波长较短,有利于布置大规模阵列,提高传输频带的利用率。在多用户毫米波大规模MIMO系统中,理论上可采用数字预编码获得最佳性能。但这需要每根天线连接一个单独的射频链路,在毫米波大规模MIMO系统中就需要大量射频链路,带来高功耗和高复杂度。为解决这个困难,通常在毫米波大规模MIMO系统中采用混合波束成形,其所需的射频链路数量远远少于天线数量。但是,多用户毫米波大规模MIMO系统中,用户间存在干扰。当干扰较强的时候,用户的通信质量将得不到保证。同时,降低干扰可能会导致接收功率较小。因此,需要在功率和干扰之间取得合理折中,从而提高系统容量。
现有方法在多用户毫米波大规模MIMO系统面临着问题。基于波束控制的混合波束成形侧重接收功率最大化,不能有效抑制用户间干扰。基于施密特正交化的混合波束成形在有限精度的相移条件下不能有效抑制用户间干扰。
发明内容
本发明的目的是针对现有混合波束成形方法不能兼顾接收功率和降低干扰的问题,提供一种基于最大化最小相位差值的混合波束成形方法。
本发明方法基于单个小区内,包含一个基站和多个用户的上行毫米波通信系统。基站天线阵列为均匀线性阵,且包含几十或几百根天线。用户为单天线。基站采用混合波束成形。基站端模拟波束成形码本有限。对任何用户而言,码本设置为有限个天线阵列响应向量的集合
Figure GDA0002409779270000011
其中,c为实数,M为基站端天线数,
Figure GDA0002409779270000012
为可实现相移器相位,即模拟波束成形角,限定模拟波束成形角的范围为[0,π)。
Figure GDA0002409779270000013
为天线阵列响应向量,
Figure GDA0002409779270000014
λ为载波波长,d为基站端相邻天线间距,j为复数,满足j2=-1。根据模拟波束成形矩阵WRF和信道特性矩阵H采用最小均方误差准则计算数字处理矩阵
Figure GDA0002409779270000021
表示WRF的共轭转置。
本发明方法的具体步骤是:
步骤1、功率保证阶段:
计算每个用户满足功率要求的码本集合:基站模拟波束成形码本为有限,首先通过路径增益门限的设置,选出高于路径增益门限的路径,接着根据这些路径对应的波达角从码本中选择与其相位差值最小的向量形成每个用户的候选模拟波束成形向量集合。具体方法是:
1.1.计算路径增益门限:
用户k的信道向量为hk
Figure GDA0002409779270000022
Lk为路径数,αkl为相应的路径增益,
Figure GDA0002409779270000023
为波达角;计算用户各路径增益模的均值并将其设置为路径增益门限εk
Figure GDA0002409779270000024
1.2.计算满足功率要求的模拟波束成形向量集合:
对于用户k,首先计算选出路径增益模高于εk的路径,即满足|αkl|>εk,然后计算对应路径波达角与码本中,向量相位差的绝对值小于
Figure GDA0002409779270000025
的向量,即满足
Figure GDA0002409779270000026
满足上述两个条件的集合即为满足功率要求的模拟波束成形向量集合Fk
Figure GDA0002409779270000027
其中θk为满足功率要求的模拟波束成形向量角,
Figure GDA0002409779270000028
为θk的集合;
Figure GDA0002409779270000029
步骤2、干扰降低阶段:
2.1.信道特性向量的近似:
将hk近似为
Figure GDA00024097792700000210
Figure GDA00024097792700000211
为最大路径增益,
Figure GDA00024097792700000212
为其对应的相位;
2.2.最大化最小相位差:
对于用户k,计算Pk中元素与其他k-1个用户最大路径增益对应波达角相位差的最小值Dkl
Figure GDA00024097792700000213
其中,
Figure GDA00024097792700000214
Tk={m|m≠k,m=1,···,K};Tk为其他k-1个用户的索引集合,m为用户索引;
Figure GDA0002409779270000031
Pk表示集合
Figure GDA0002409779270000032
去除已经选出的最优波束成形角度
Figure GDA0002409779270000033
后剩下的集合,t为实数,表示已选出的最优波束成形角的索引;
计算出所有Dkl,找出最大值所对应的
Figure GDA0002409779270000034
中元素作为最优波束成形角度
Figure GDA0002409779270000035
Figure GDA0002409779270000036
得到模拟波束成形向量,
Figure GDA0002409779270000037
[WRF]k表示模拟波束成形矩阵WRF的第k列;按照k=1,···,K的顺序重复步骤1和2得到完整的WRF矩阵。
本发明方法的创新关键是:利用路径增益门限及路径波达角与码本最小相位差,保证了较高的接收信号功率;最大化最小相位差避免了最强干扰。
本发明利用信道增益门限从模拟波束成形向量码本中筛选出能够提供较高信号功率的向量集合,保证了较大的接收功率。利用模拟波束成形向量与主要传播路径阵列响应向量的相位差衡量用户间干扰,采用最大化最小相位差值的准则抑制最强用户间干扰,进一步提高了传输效率。同时,最大化最小相位差算法,优化范围小,计算复杂度低,具有很好的可实现性。本发明方法能够减轻用户间的干扰,提高系统容量,同时保证较低复杂度。
附图说明
图1是本发明实施例中系统和速率与终端数量关系的仿真图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面结合附图对本发明作进一步的详细描述。
毫米波大规模多输入多输出多用户系统中的基于最大化最小相位差值的混合波束成形方法,应用场景如下:
单个小区内,K个单天线用户同时向基站发送数据,基站配有M根天线和K个射频链路用于接收信号,阵列为均匀线形阵。
发送信号向量为s∈CK×1,基站接收到的信号经过混合波束成形矩阵处理,即WRF∈CM×K,及WRF∈CK×K来恢复向量s,恢复后的向量表示为
Figure GDA0002409779270000038
ρ表示发送信号功率,n∈CM×1为接收噪声向量。对任何用户而言,码本设置为有限个天线阵列响应向量的集合
Figure GDA0002409779270000039
Figure GDA00024097792700000310
为天线阵列响应向量,
Figure GDA00024097792700000311
λ为载波波长,d为基站端相邻天线间距,j为复数,满足j2=-1。根据模拟波束成形矩阵WRF和信道特性矩阵H采用最小均方误差准则计算数字处理矩阵
Figure GDA0002409779270000041
表示WRF的共轭转置。
该方法的具体步骤是:
步骤1、功率保证阶段:
计算每个用户满足功率要求的码本集合:基站模拟波束成形码本为有限,首先通过路径增益门限的设置,选出高于路径增益门限的路径,接着根据这些路径对应的波达角从码本中选择与其相位差值最小的向量形成每个用户的候选模拟波束成形向量集合。具体方法是:
1.1.计算路径增益门限:
用户k的信道向量为hk
Figure GDA0002409779270000042
Lk为路径数,αkl为相应的路径增益,
Figure GDA0002409779270000043
为波达角;计算用户各路径增益模的均值并将其设置为路径增益门限εk
Figure GDA0002409779270000044
1.2.计算满足功率要求的模拟波束成形向量集合:
对于用户k,首先计算选出路径增益模高于εk的路径,即满足|αkl|>εk,然后计算对应路径波达角与码本中,向量相位差的绝对值小于
Figure GDA0002409779270000045
的向量,即满足
Figure GDA0002409779270000046
满足上述两个条件的集合即为满足功率要求的模拟波束成形向量集合Fk
Figure GDA0002409779270000047
其中θk为满足功率要求的模拟波束成形向量角,
Figure GDA0002409779270000048
为θk的集合;
Figure GDA0002409779270000049
步骤2、干扰降低阶段:
2.1.信道特性向量的近似:
由于hk为多个阵列响应向量与路径增益乘积累加的形式,不能进行相位比较;将hk近似为
Figure GDA00024097792700000410
为最大路径增益,
Figure GDA00024097792700000411
为其对应的相位;
2.2.最大化最小相位差:
对于用户k,计算Pk中元素与其他k-1个用户最大路径增益对应波达角相位差的最小值Dkl
Figure GDA00024097792700000412
其中,
Figure GDA00024097792700000413
Tk={m|m≠k,m=1,···,K};Tk为其他k-1个用户的索引集合,m为用户索引;
Figure GDA0002409779270000051
Pk表示集合
Figure GDA0002409779270000052
去除已经选出的最优波束成形角度
Figure GDA0002409779270000053
后剩下的集合,t为实数,表示已选出的最优波束成形角的索引;
计算出所有Dkl,找出最大值所对应的
Figure GDA0002409779270000054
中元素作为最优波束成形角度
Figure GDA0002409779270000055
Figure GDA0002409779270000056
得到模拟波束成形向量,
Figure GDA0002409779270000057
[WRF]k表示模拟波束成形矩阵WRF的第k列;按照k=1,···,K的顺序重复步骤1和2得到完整的WRF矩阵。
仿真试验的结果如图1所示,对终端数量的影响进行了仿真。为了直观地体现出本发明方法的优越性,将本方法的仿真结果与现有的波束控制、施密特正交化的修正方法进行了对比。
由图1可以看出,系统和速率随着用户数的增多递增,所提出算法性能一直优于施密特正交化修正算法和波束控制算法。并且,施密特正交化修正算法和波束控制算法的算法性能随着用户的增多渐趋平稳,说明对干扰的消除不够好。反观本发明所提出的算法,其性能随着用户数的增长稳步上升,这也反映出所提出算法在降低用户间干扰方面性能的优越性。本发明方法明显地提高了系统的和速率,也就是说该方法能够显著降低干扰。

Claims (1)

1.一种基于最大化最小相位差值的混合波束成形方法,该方法基于单个小区内,包含一个基站和多个用户的上行毫米波通信系统;基站天线阵列为均匀线性阵,用户为单天线,基站采用混合波束成形;对任何用户而言,码本设置为有限个天线阵列响应向量的集合
Figure FDA0002409779260000011
其中c为实数,M为基站端天线数,
Figure FDA0002409779260000012
为可实现相移器相位,即模拟波束成形角,限定模拟波束成形角的范围为[0,π),
Figure FDA0002409779260000013
为天线阵列响应向量,
Figure FDA0002409779260000014
λ为载波波长,d为基站端相邻天线间距,j为复数,满足j2=-1,根据模拟波束成形矩阵WRF和信道特性矩阵H采用最小均方误差准则计算数字处理矩阵
Figure FDA0002409779260000015
Figure FDA0002409779260000016
表示WRF的共轭转置;其特征在于该方法的具体步骤是:
步骤1、功率保证阶段:
计算每个用户满足功率要求的码本集合:基站模拟波束成形码本为有限,首先通过路径增益门限的设置,选出高于路径增益门限的路径,接着根据这些路径对应的波达角从码本中选择与其相位差值最小的向量形成每个用户的候选模拟波束成形向量集合;具体方法是:
1.1.计算路径增益门限:
用户k的信道向量为hk
Figure FDA0002409779260000017
Lk为路径数,αkl为相应的路径增益,
Figure FDA0002409779260000018
为波达角;计算用户各路径增益模的均值并将其设置为路径增益门限εk
Figure FDA0002409779260000019
1.2.计算满足功率要求的模拟波束成形向量集合:
对于用户k,首先计算选出路径增益模高于εk的路径,即满足|αkl|>εk,然后计算对应路径波达角与码本中,向量相位差的绝对值小于
Figure FDA00024097792600000110
的向量,即满足
Figure FDA00024097792600000111
满足上述两个条件的集合即为满足功率要求的模拟波束成形向量集合Fk
Figure FDA00024097792600000112
其中θk为满足功率要求的模拟波束成形向量角,
Figure FDA00024097792600000113
为θk的集合;
Figure FDA0002409779260000021
步骤2、干扰降低阶段:
2.1.信道特性向量的近似:
将hk近似为
Figure FDA0002409779260000022
Figure FDA0002409779260000023
为最大路径增益,
Figure FDA0002409779260000024
为其对应的相位;
2.2.最大化最小相位差:
对于用户k,计算Pk中元素与其他k-1个用户最大路径增益对应波达角相位差的最小值Dkl
Figure FDA0002409779260000025
其中,
Figure FDA0002409779260000026
Tk={m|m≠k,m=1,···,K};Tk为其他k-1个用户的索引集合,m为用户索引;
Figure FDA0002409779260000027
Pk表示集合
Figure FDA0002409779260000028
去除已经选出的最优波束成形角度
Figure FDA0002409779260000029
后剩下的集合,t为实数,表示已选出的最优波束成形角的索引;
计算出所有Dkl,找出最大值所对应的
Figure FDA00024097792600000213
中元素作为最优波束成形角度
Figure FDA00024097792600000210
Figure FDA00024097792600000211
得到模拟波束成形向量,
Figure FDA00024097792600000212
[WRF]k表示模拟波束成形矩阵WRF的第k列;按照k=1,···,K的顺序重复步骤1和2得到完整的WRF矩阵。
CN201811020499.7A 2018-09-03 2018-09-03 一种基于最大化最小相位差值的混合波束成形方法 Active CN109120316B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811020499.7A CN109120316B (zh) 2018-09-03 2018-09-03 一种基于最大化最小相位差值的混合波束成形方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811020499.7A CN109120316B (zh) 2018-09-03 2018-09-03 一种基于最大化最小相位差值的混合波束成形方法

Publications (2)

Publication Number Publication Date
CN109120316A CN109120316A (zh) 2019-01-01
CN109120316B true CN109120316B (zh) 2020-05-22

Family

ID=64861816

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811020499.7A Active CN109120316B (zh) 2018-09-03 2018-09-03 一种基于最大化最小相位差值的混合波束成形方法

Country Status (1)

Country Link
CN (1) CN109120316B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109787828B (zh) * 2019-01-18 2021-08-31 杭州电子科技大学 认知swipt最佳节点选择与波束成型联合设计方法
CN111698005B (zh) * 2019-03-15 2022-07-29 华为技术有限公司 用于获取波束权值的方法和装置
CN110492912B (zh) * 2019-07-22 2020-09-08 杭州电子科技大学 一种基于分组优化的混合波束成形方法
CN110365388B (zh) * 2019-07-31 2022-07-29 东南大学 一种低复杂度毫米波多播波束成形方法
US20230170959A1 (en) * 2021-12-01 2023-06-01 Mediatek Inc. Method and apparatus for hybrid beamforming with autonomous beamformers in mobile communications

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107566303A (zh) * 2017-07-27 2018-01-09 东华大学 一种基于贝叶斯压缩感知的毫米波信道估计方法
CN108306667A (zh) * 2017-12-25 2018-07-20 南京理工大学 大规模mimo中基于线性搜索的混合到达角估计

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9444534B2 (en) * 2012-02-06 2016-09-13 Samsung Electronics Co., Ltd. Apparatus and method for low complexity spatial division multiple access in a millimeter wave mobile communication system
US9148325B2 (en) * 2013-05-23 2015-09-29 Ridha HAMILA System and methods for compensation of I/Q imbalance in beamforming OFDM systems
EP3123626B1 (en) * 2014-03-24 2019-05-01 LG Electronics Inc. Method of performing a hybrid beamforming in a wireless communication system and apparatus therefor
CN107453795B (zh) * 2017-08-21 2020-06-30 东南大学 多用户毫米波通信系统的波束分配方法及其装置和系统
CN108270471A (zh) * 2017-12-27 2018-07-10 东南大学 大规模天线阵列系统中低反馈混合预编码与接收合并方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107566303A (zh) * 2017-07-27 2018-01-09 东华大学 一种基于贝叶斯压缩感知的毫米波信道估计方法
CN108306667A (zh) * 2017-12-25 2018-07-20 南京理工大学 大规模mimo中基于线性搜索的混合到达角估计

Also Published As

Publication number Publication date
CN109120316A (zh) 2019-01-01

Similar Documents

Publication Publication Date Title
CN109120316B (zh) 一种基于最大化最小相位差值的混合波束成形方法
Li et al. Hybrid precoding and combining design for millimeter-wave multi-user MIMO based on SVD
CN107046434B (zh) 大规模mimo系统模数混合预编码方法
CN107453795B (zh) 多用户毫米波通信系统的波束分配方法及其装置和系统
CN100392998C (zh) 结合空时分组编码的智能天线下行波束形成方法
CN110048752B (zh) 毫米波大规模mimo中的一种混合预编码方法
CN112039565B (zh) 基于分布式部分连接的大规模mimo混合预编码方法
CN108736943B (zh) 一种适用于大规模mimo系统的混合预编码方法
CN111049557B (zh) 基于统计信道信息的毫米波mimo系统混合预编码方法
CN110138425B (zh) 低复杂度阵列天线多输入多输出系统混合预编码算法
CN107809275B (zh) 一种基于毫米波mimo系统的有限反馈混合预编码方法
CN110650103B (zh) 利用冗余字典加强稀疏性的透镜天线阵列信道估计方法
CN109560845B (zh) 一种低复杂度通用混合预编码方法
CN111953393A (zh) 大规模mimo混合预编码器及匹配关系
Singh et al. A survey on hybrid beamforming in MmWave massive MIMO system
CN112636794A (zh) 一种基于贪心算法的波束赋形方法
CN110233649B (zh) 一种基于毫米波mimo系统的动态子阵列设计方法
CN109361434B (zh) 基站协作传输的毫米波mimo混合预编码方法
Wang et al. Hybrid beamforming with time delay compensation for millimeter wave MIMO frequency selective channels
CN108566236B (zh) 用户终端、基站以及混合波束成形传输方法和系统
CN113708811A (zh) 一种毫米波大规模mimo系统中的混合预编码设计方法
CN110492912B (zh) 一种基于分组优化的混合波束成形方法
CN108683441B (zh) 混合预编码中的多用户波束赋形方法
CN112398513A (zh) 一种massive MIMO系统的波束赋形方法
Yan et al. Multiple RF chains assisted parallel beam search for mmWave hybrid beamforming 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

Effective date of registration: 20210511

Address after: 710000 Shaanxi Xi'an economic and Technological Development Zone, Fengcheng four road, Ming Guang road, southeast corner of the new business building 1808

Patentee after: SHAANXI ZHUANYI INTELLECTUAL PROPERTY OPERATION Co.,Ltd.

Address before: 310018 No. 2 street, Xiasha Higher Education Zone, Hangzhou, Zhejiang

Patentee before: HANGZHOU DIANZI University

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20230109

Address after: Room 403-01, Block D, Software Park, No. 228, Nanyi Road, Dongying Development Zone, Shandong 257091

Patentee after: Dongying Zhengnuo Technology Service Co.,Ltd.

Address before: 710000 Shaanxi Xi'an economic and Technological Development Zone, Fengcheng four road, Ming Guang road, southeast corner of the new business building 1808

Patentee before: SHAANXI ZHUANYI INTELLECTUAL PROPERTY OPERATION CO.,LTD.

TR01 Transfer of patent right