CN102707295A - 一种检测北斗卫星b1频段民用信号发生互相关的方法 - Google Patents

一种检测北斗卫星b1频段民用信号发生互相关的方法 Download PDF

Info

Publication number
CN102707295A
CN102707295A CN2012102129003A CN201210212900A CN102707295A CN 102707295 A CN102707295 A CN 102707295A CN 2012102129003 A CN2012102129003 A CN 2012102129003A CN 201210212900 A CN201210212900 A CN 201210212900A CN 102707295 A CN102707295 A CN 102707295A
Authority
CN
China
Prior art keywords
satellite
simple crosscorrelation
cross
correlation
detection
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
Application number
CN2012102129003A
Other languages
English (en)
Other versions
CN102707295B (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.)
Shanghai Xingsi Semiconductor Co ltd
Original Assignee
XIAN HUAXUN MICROELECTRONIC CO Ltd
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 XIAN HUAXUN MICROELECTRONIC CO Ltd filed Critical XIAN HUAXUN MICROELECTRONIC CO Ltd
Priority to CN 201210212900 priority Critical patent/CN102707295B/zh
Publication of CN102707295A publication Critical patent/CN102707295A/zh
Application granted granted Critical
Publication of CN102707295B publication Critical patent/CN102707295B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Position Fixing By Use Of Radio Waves (AREA)
  • Radio Relay Systems (AREA)

Abstract

本发明公开了一种检测北斗卫星B1频段民用信号发生互相关的方法，包括接收机在成功捕获到一颗新卫星，并通过校验和位同步之后，开始做互相关的检测，若该卫星的测量值满足预先设定好的互相关检测标准，则认为已检测到了互相关，否则认为互相关不存在。这样就可以尽早检测出互相关的存在，报告该卫星捕获失败，而不必再转入跟踪阶段，甚至造成假锁,从而提高接收机效率；另一方面，可告知基带信号处理模块开启互相关模式，做抑制互相关的处理。

Description

一种检测北斗卫星B1频段民用信号发生互相关的方法
技术领域
本发明涉及一种检测北斗导航系统1-37号卫星B1频点民用信号发生互相关的方法。
背景技术
北斗B1频点民码信号为短伪随机码序列设计的Gold码，以完成快速捕获。通常情况下此类短随机码具有良好的自相关和互相关特性，可满足大多数应用场合的要求,但在强弱信号共存的情况下,强弱卫星信号间的互相关效应将导致民用接收机无法正常捕获和跟踪弱卫星信号。
假设接收信号中仅包含1颗强信号卫星和1颗弱信号卫星，接收信号模型为：
Figure BDA00001809261500011
其中A1是弱信号幅度，A2是强信号幅度，当捕获弱信号的时候，经下变频、解扩和相干积分，得到I、Q两个支路信号分别为：
Figure BDA00001809261500012
式中第一项代表自相关，第二项代表互相关，最后一项是白噪声。从中可看出互相关效应的大小主要与A1、A2（强弱卫星信号强度）、τ2强信号的伪码延迟和ωd2强信号的多普勒频率有关。北斗B1频点普通测距码是码长为2046chip周期为1ms的Gold码，其互相关函数的旁瓣比自相关峰低33dB左右，33dB是互相关函数的标准差，也就是说互相关函数的旁瓣大多数都比自相关峰低33dB左右，也有某些旁瓣仅比自相关峰低29.8dB，最坏的情况是有多普勒频率时，旁瓣比自相关峰只低24.6dB，这时再加上强弱信号接收功率的影响，将会导致互相关。
发明内容
本发明要解决的问题是提供一种检测北斗B1频点民用信号发生互相关的方法，尽早准确的检测出互相关的存在，报告该卫星捕获失败，而不必再转入跟踪阶段，甚至造成假锁，从而提高接收机效率；另一方面，可告知基带信号处理模块开启互相关模式，做抑制互相关的处理。
为达到以上目的，本发明是采取如下技术方案予以实现的。
一种检测北斗卫星B1频段民用信号发生互相关的方法，包括下述步骤：
步骤1：接收机捕获到一颗新卫星后，做强信号或中弱信号校验和位同步，校验和位同步成功之后，开始互相关的检测，否则重新捕获卫星；
步骤2：把已跟踪到的各颗卫星的测量值放入互相关检测数据库，以便与待检测卫星的测量值作比较；
步骤3：待检测卫星的载噪比C/N0与预先设定好的阈值作比较，若大于此阈值，则认为互相关不存在，待检测卫星捕获成功，转入跟踪阶段；若小于或等于此阈值，则认为互相关可能存在，继续以下步骤的检测；
步骤4：待检测卫星的测量值与互相关检测数据库中各卫星的测量值作比较，若检测到互相关存在，则停止互相关的检测，报告该卫星捕获失败，并开启互相关模式，抑制互相关；
步骤5：如果对互相关数据库中的卫星轮循一次之后，仍未检测到互相关，那么继续检测该卫星是否与GEO卫星发生了互相关；
这时只需要判断待检测卫星的载噪比C/N0和多普勒频率Doppler Offset是否满足相应的GEO卫星互相关检测标准即可，若符合判断条件，则认为检测到互相关，报告该卫星捕获失败，开启互相关模式，抑制互相关；否则认为互相关不存在，该卫星捕获成功，转入跟踪阶段。
进一步的，所述已跟踪到的卫星至少有一颗。
进一步的，所述步骤3中将预先设定的阈值对完成捕获、检验和位同步的卫星进行互相关检测，预先设定的阈值为35dB-Hz。
进一步的，所述步骤4中互相关检测数据库中各卫星的测量值包括载噪比C/N0、多普勒频率Doppler Offset和导航比特之后的第几毫秒值Bit Offset；
其中，多普勒频率Doppler Offset数值换算为FreqDiff值进行判断，FreqDiff值按照下式进行换算：
FreqDiff=(Doppler1-Doppler2)%1000
式中，Doppler1为待检测卫星的多普勒频率，Doppler2为互相关数据库中某颗卫星的多普勒频率。
进一步的，当待检测卫星和互相关数据库中某颗卫星的载噪比C/N0的差值和导航比特之后的第几毫秒Bit Offset差值，以及FreqDiff值，符合互相关检测标准中某一组判断条件时，就认为检测到互相关；
所述判断条件如下：
Table1互相关检测判断标准
  C/N0DiffMin   FreqDiffMin   FreqDiffMax   Bit OffsetMin
  22   0   25   0
  28   25   40   0
  32   40   70   0
所述步骤5中GEO卫星的互相关检测判断标准：
Table 1GEO卫星互相关检测判断标准
  C/N0DiffMin   FreqDiffMin   FreqDiffMax   Bit OffsetMin
  30   0   220   0
进一步的，待检测卫星的多普勒频率与互相关数据库中某颗卫星的多普勒频率二者的差，再对1000取模，所得的值应属于某一特定范围；所述FreqDiffMin即为该范围的下限值，所述FreqDiffMax为该范围的上限值。
进一步的，所述互相关数据库中某颗卫星的第几毫秒值bit offset除GEO卫星的第几毫秒值Bit Offset为0或1以外，其他卫星的第几毫秒值bit offset均为0~19的整数。
本发明的优点是，该方法能够准确有效的检测出北斗导航系统1-37号卫星发射的B1频点民用信号之间是否发生了互相关，检测标准的确定综合考虑了各种因素，有非常严格的理论依据；此外，根据GEO卫星的特性，使用特定的检测标准对它进行单独检测，这样能更有效的检测出互相关。
附图说明
图1是本发明检测卫星信号是否发生互相关的流程。
具体实施方式
以下结合具体实施例对本发明作进一步的详细描述。
一种检测卫星信号发生互相关的方法，参见图1所示，包括下述步骤：
步骤1：接收机捕获一颗新卫星信号，做强信号或中信号的校验及位同步，若校验及位同步成功，则进入下一步；若校验及位同步失败，则返回，重新捕获新卫星信号；
步骤2：把已跟踪到的各颗卫星的测量值放入互相关检测数据库，以便与待检测卫星的测量值作比较，此处已跟踪到的卫星至少有一颗；
步骤3：待检测卫星的载噪比C/N0与预先设定好的阈值作比较，此处预先设定的阈值为35dB-Hz；若大于此阈值，则认为互相关不存在，待检测卫星捕获成功，转入跟踪阶段；若小于或等于此阈值，则认为互相关可能存在，继续以下步骤的检测；
步骤4：待检测卫星的测量值与互相关检测数据库中各卫星的测量值作比较，若检测到互相关存在，则停止互相关的检测，报告该卫星捕获失败，并开启互相关模式，抑制互相关；
互相关检测数据库中各卫星的测量值包括载噪比C/N0、多普勒频率DopplerOffset和导航比特之后的第几毫秒值Bit Offset；
其中，多普勒频率Doppler Offset数值换算为FreqDiff值进行判断，FreqDiff值按照下式进行换算：
FreqDiff=(Doppler1-Doppler2)%1000
式中，Doppler1为待检测卫星的多普勒频率，Doppler2为互相关数据库中某颗卫星的多普勒频率。
当待检测卫星和互相关数据库中某颗卫星的载噪比C/N0的差值和导航比特之后的第几毫秒Bit Offset差值，以及FreqDiff值，符合互相关检测标准中某一组判断条件时，就认为检测到互相关；
所述判断条件如下：
Table1互相关检测判断标准
  C/N0DiffMin   FreqDiffMin   FreqDiffMax   Bit OffsetMin
  22   0   25   0
  28   25   40   0
  32   40   70   0
所述步骤5中GEO卫星的互相关检测判断标准：
Table 1GEO卫星互相关检测判断标准
  C/N0DiffMin   FreqDiffMin   FreqDiffMax   Bit OffsetMin
  30   0   220   0
上述FreqDiffMin,FreqDiffMin分别表示待检测卫星的多普勒频率与互相关数据库中某颗卫星的多普勒频率二者的差，再对1000取模，所得的值；所述FreqDiffMin即为该范围的下限值，所述FreqDiffMax为该范围的上限值。
互相关数据库中某颗卫星的第几毫秒值bit offset除GEO卫星的第几毫秒值Bit Offset为0或1以外，其他卫星的第几毫秒值bit offset均为0~19的整数。
步骤5：如果对互相关数据库中的卫星轮循一次之后，仍未检测到互相关，那么继续检测该卫星是否与GEO卫星发生了互相关；
这时只需要判断待检测卫星的载噪比C/N0和多普勒频率Doppler Offset是否满足相应的GEO卫星互相关检测标准即可，若符合判断条件，则认为检测到互相关，报告该卫星捕获失败，开启互相关模式，抑制互相关；否则认为互相关不存在，该卫星捕获成功，转入跟踪阶段。
该检测卫星信号发生互相关的方法和流程如上述步骤1-5所述，这里主要阐述确定互相关检测标准的方法。
假设接收信号中仅包含1颗强信号卫星和1颗弱信号卫星，接收信号模型为：
Figure BDA00001809261500061
当捕获弱信号的时候，经下变频、解扩和相干积分，得到I、Q两个支路信号分别为：
Figure BDA00001809261500062
上述公式中各符号的意义在背景技术部分已有说明。注意观察上两式的第二项为互相关部分，且可以合写为：
Figure BDA00001809261500063
其中，T代表北斗B1频点普通测距码的周期，N为自然数，两个Gold码相乘结果是另一个Gold码，但有不同的延时，是一个新的Gold码，周期同样为1ms，
因此，可把上式表示为傅里叶级数的形式：
Figure BDA00001809261500071
其中，ak是1个周期的Gold码基带波形的傅里叶系数，
Figure BDA00001809261500072
从式中可以看出，互相关的影响由离散谱ak和
Figure BDA00001809261500073
组成，而ak由两个伪码的延时差决定，
Figure BDA00001809261500074
则由多普率频率决定，经分析有以下结论：
强信号的多普勒频率等于1KHz的整数倍时，互相关的影响最为严重，程度由第k根离散谱ak决定；当多普勒频率不是1KHz的整数倍时，互相关的影响与积分时间NT有关，当积分时间为1ms的整数倍时，积分结果中最多有2根谱线发生作用，该谱线使得
Figure BDA00001809261500075
当待检测卫星与互相关数据库中某颗卫星的FreqDiff在[0,25Hz]范围内时，说明二者的多普勒频率的差非常接近1KHz的整数倍，很有可能发生了较为严重的互相关；接下来再比较二者的C/N0，如果待检测卫星和相关通道数据库某颗已跟踪到的卫星的C/N0相差20dB-Hz以内时，再比较一下测量值Bit Offset，如果二者的Bit Offset值也相等，那么就认为检测到了互相关。
当检测发生互相关的卫星是否为GEO卫星时，确定检测标准的方法和思路与上述的方法基本一致，但还需考虑到GEO卫星相对于用户来说，多普勒频率一般不会超过220Hz，判断时直接看待检测卫星的多普率频率与1KHz整数倍的接近程度；此外，GEO卫星的相关积分时间较短，在确定C/N0的判断条件时应考虑到这一因素。
以上所述仅为本发明较佳实施例，然其并非用以限定本发明的范围，任何书须本项技术的人员，在不脱离本发明的精神和范围内，可在此基础上做进一步的改进和变化，因为本发明的保护范围当以本申请的权利要求书所界定的范围为准。

Claims (7)

1.一种检测北斗卫星B1频段民用信号发生互相关的方法，其特征在于，包括下述步骤：
步骤1：接收机捕获到一颗新卫星后，做强信号或中弱信号校验和位同步，校验和位同步成功之后，开始互相关的检测，否则重新捕获卫星；
步骤2：把已跟踪到的各颗卫星的测量值放入互相关检测数据库，以便与待检测卫星的测量值作比较；
步骤3：待检测卫星的载噪比C/N0与预先设定好的阈值作比较，若大于此阈值，则认为互相关不存在，待检测卫星捕获成功，转入跟踪阶段；若小于或等于此阈值，则认为互相关可能存在，继续以下步骤的检测；
步骤4：待检测卫星的测量值与互相关检测数据库中各卫星的测量值作比较，若检测到互相关存在，则停止互相关的检测，报告该卫星捕获失败，并开启互相关模式，抑制互相关；
步骤5：如果对互相关数据库中的卫星轮循一次之后，仍未检测到互相关，那么继续检测该卫星是否与GEO卫星发生了互相关；
这时只需要判断待检测卫星的载噪比C/N0和多普勒频率Doppler Offset是否满足相应的GEO卫星互相关检测标准即可，若符合判断条件，则认为检测到互相关，报告该卫星捕获失败，开启互相关模式，抑制互相关；否则认为互相关不存在，该卫星捕获成功，转入跟踪阶段。
2.根据权利要求1所述的检测北斗卫星B1频段民用信号发生互相关的方法，其特征在于，所述已跟踪到的卫星至少有一颗。
3.根据权利要求1所述的检测北斗卫星B1频段民用信号发生互相关的方法，其特征在于，所述步骤3中将预先设定的阈值对完成捕获、检验和位同步的卫星进行互相关检测，预先设定的阈值为35dB-Hz。
4.根据权利要求1所述的检测北斗卫星B1频段民用信号发生互相关的方法，其特征在于，所述步骤4中互相关检测数据库中各卫星的测量值包括载噪比C/N0、多普勒频率Doppler Offset和导航比特之后的第几毫秒值Bit Offset；
其中，多普勒频率Doppler Offset数值换算为FreqDiff值进行判断，FreqDiff值按照下式进行换算：
FreqDiff=(Doppler1-Doppler2)%1000
式中，Doppler1为待检测卫星的多普勒频率，Doppler2为互相关数据库中某颗卫星的多普勒频率。
5.根据权利要求4所述的检测北斗卫星B1频段民用信号发生互相关的方法，其特征在于，当待检测卫星和互相关数据库中某颗卫星的载噪比C/N0的差值和导航比特之后的第几毫秒Bit Offset差值，以及FreqDiff值，符合互相关检测标准中某一组判断条件时，就认为检测到互相关；
所述判断条件如下：
Table1互相关检测判断标准
  C/N0DiffMin   FreqDiffMin   FreqDifMax   Bit OffsetMin   22   0   25   0   28   25   40   0   32   40   70   0
所述步骤5中GEO卫星的互相关检测判断标准：
Table 1GEO卫星互相关检测判断标准
  C/N0DiffMin   FreqDiffMin   FreqDiffMax   Bit OffsetMin   30   0   220   0
6.根据权利要求5所述的检测北斗卫星B1频段民用信号发生互相关的方法，其特征在于，待检测卫星的多普勒频率与互相关数据库中某颗卫星的多普勒频率二者的差，再对1000取模，所得的值应属于某一特定范围；所述FreqDiffMin即为该范围的下限值，所述FreqDiffMax为该范围的上限值。
7.根据权利要求5所述的检测北斗卫星B1频段民用信号发生互相关的方法，其特征在于，所述互相关数据库中某颗卫星的第几毫秒值bit offset除GEO卫星的第几毫秒值Bit Offset为0或1以外，其他卫星的第几毫秒值bit offset均为0~19的整数。
CN 201210212900 2012-06-26 2012-06-26 一种检测北斗卫星b1频段民用信号发生互相关的方法 Active CN102707295B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210212900 CN102707295B (zh) 2012-06-26 2012-06-26 一种检测北斗卫星b1频段民用信号发生互相关的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210212900 CN102707295B (zh) 2012-06-26 2012-06-26 一种检测北斗卫星b1频段民用信号发生互相关的方法

Publications (2)

Publication Number Publication Date
CN102707295A true CN102707295A (zh) 2012-10-03
CN102707295B CN102707295B (zh) 2013-12-25

Family

ID=46900224

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210212900 Active CN102707295B (zh) 2012-06-26 2012-06-26 一种检测北斗卫星b1频段民用信号发生互相关的方法

Country Status (1)

Country Link
CN (1) CN102707295B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10281585B2 (en) 2015-12-07 2019-05-07 Stmicroelectronics S.R.L. Method for the acquisition of a GEO (geostationary earth orbit) satellite and corresponding receiver
CN116626716A (zh) * 2023-07-21 2023-08-22 齐鲁空天信息研究院 一种北斗信号的跟踪监测方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6633255B2 (en) * 1995-10-09 2003-10-14 Qualcomm Inc. Method for open loop tracking GPS signals
CN101206255A (zh) * 2007-12-13 2008-06-25 西安华迅微电子有限公司 一种捕获gps弱信号的方法
CN101852849A (zh) * 2009-04-02 2010-10-06 深圳艾科创新微电子有限公司 一种gps弱信号捕获的方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6633255B2 (en) * 1995-10-09 2003-10-14 Qualcomm Inc. Method for open loop tracking GPS signals
CN101206255A (zh) * 2007-12-13 2008-06-25 西安华迅微电子有限公司 一种捕获gps弱信号的方法
CN101852849A (zh) * 2009-04-02 2010-10-06 深圳艾科创新微电子有限公司 一种gps弱信号捕获的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
张炳琪 等: "北斗一号卫星信号高灵敏度捕获方法", 《北京理工大学学报》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10281585B2 (en) 2015-12-07 2019-05-07 Stmicroelectronics S.R.L. Method for the acquisition of a GEO (geostationary earth orbit) satellite and corresponding receiver
CN116626716A (zh) * 2023-07-21 2023-08-22 齐鲁空天信息研究院 一种北斗信号的跟踪监测方法
CN116626716B (zh) * 2023-07-21 2023-09-29 齐鲁空天信息研究院 一种北斗信号的跟踪监测方法

Also Published As

Publication number Publication date
CN102707295B (zh) 2013-12-25

Similar Documents

Publication Publication Date Title
CN101295999B (zh) 一种随机接入前导的检测方法
CN104143997B (zh) 一种基于扩频体制的多址干扰抗错锁判决方法
CN108008419A (zh) 基于fpga的抗转发式欺骗干扰方法及其检测系统
CN107786240B (zh) 多重载波频率跳变的数字通信电台装置
CN104793221A (zh) 卫星信号接收机中实现快速捕获失锁信号的方法
CN104270167A (zh) 一种基于多维特征神经网络的信号检测和估计方法
CN109088838A (zh) 一种高动态下直扩dpsk信号的伪码—多普勒快速捕获方法
CN102546499A (zh) 一种实线性调频信号的分数阶信道化接收方法
CN102707295B (zh) 一种检测北斗卫星b1频段民用信号发生互相关的方法
CN103427871B (zh) 一种抑制Smart AGC抗干扰能力的转发式干扰方法
CN102307055B (zh) Dsss频域干扰检测方法
CN102655419B (zh) 一种捕获扩频信号的判决函数计算方法
CN103259561B (zh) 一种大型建筑物变形遥测技术中多址干扰抵消方法
CN102185628A (zh) 基于自适应功率积累的扩频码相位捕获装置及其捕获方法
O'Mahony et al. Constrained optimization of coincidence detector parameters for GPS acquisition
CN105871530A (zh) 一种抗高动态频偏变化的定时跟踪方法
CN103024795B (zh) 一种td-scdma中同步定时异常的检测方法
CN107202996B (zh) 基于多点相关法的卫星导航抗欺骗干扰设计方法
CN101562463B (zh) 一种扩频捕获方法及其实现的接收端
CN102984719A (zh) 基于预处理的主同步信号检测方法、装置及终端
CN105072075A (zh) 一种多门限判决的ofdm同步方法
CN105119631B (zh) 一种基于频谱检测的扩频接收机多址干扰抑制方法
CN104849735B (zh) 一种组合导航接收机的频率假锁检测装置和方法
CN105656511A (zh) 一种适应于有频偏和低信噪比环境下的差分相关捕获方法
CN103501187A (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
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20121003

Assignee: Jiangsu Ling bit Microelectronic Technology Co.,Ltd.

Assignor: Xi'an Huaxun Microelectronics Co.,Ltd.

Contract record no.: 2014320000224

Denomination of invention: Method for detecting cross-correlation of civil signals in BeiDou satellite B1 frequency band

Granted publication date: 20131225

License type: Exclusive License

Record date: 20140319

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model
C56 Change in the name or address of the patentee

Owner name: XI'AN AEROSPACE HUAXUN TECHNOLOGY CO., LTD.

Free format text: FORMER NAME: XIAN HUAXUN MICROELECTRONIC CO., LTD.

CP03 Change of name, title or address

Address after: 710075 hi tech Zone, Shaanxi, Xi'an province 58, No. three, Hui International three floor

Patentee after: XI'AN AEROSPACE HUAXUN TECHNOLOGY Co.,Ltd.

Address before: 710075, Xi'an hi tech Zone, Shaanxi Province three science and Technology Road 58, Hao Hui International third floor

Patentee before: Xi'an Huaxun Microelectronics Co.,Ltd.

TR01 Transfer of patent right

Effective date of registration: 20230309

Address after: Room 401 and Room 402, Building 5, No. 690, Bibo Road, Pilot Free Trade Zone, Pudong New Area, Shanghai, March 2012

Patentee after: Shanghai Xingsi Semiconductor Co.,Ltd.

Address before: 710075 Floor 3, Huihao International, No. 58, Keji 3rd Road, High-tech Zone, Xi'an, Shaanxi Province

Patentee before: XI'AN AEROSPACE HUAXUN TECHNOLOGY Co.,Ltd.

TR01 Transfer of patent right