CN108871181B - 一种多对极磁电编码器动态多窗口区间预测角度细分方法 - Google Patents
一种多对极磁电编码器动态多窗口区间预测角度细分方法 Download PDFInfo
- Publication number
- CN108871181B CN108871181B CN201810424798.0A CN201810424798A CN108871181B CN 108871181 B CN108871181 B CN 108871181B CN 201810424798 A CN201810424798 A CN 201810424798A CN 108871181 B CN108871181 B CN 108871181B
- Authority
- CN
- China
- Prior art keywords
- window
- angle value
- theta
- pair
- pole
- 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
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B7/00—Measuring arrangements characterised by the use of electric or magnetic techniques
- G01B7/30—Measuring arrangements characterised by the use of electric or magnetic techniques for measuring angles or tapers; for testing the alignment of axes
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Transmission And Conversion Of Sensor Element Output (AREA)
Abstract
一种多对极磁电编码器动态多窗口区间预测角度细分方法,用于准确判断磁电编码器多对极角度对应的极数;本发明的方法包括:(1)解算单对极及多对极磁电编码器角度值;(2)依据多对极磁电编码器角度值建立动态多窗口;(3)依据动态多窗口解算编码器的多对极角度值对应的极数;(4)依据动态多窗口、单对极角度值、多对极角度值进行多对极磁电编码器角度值细分。
Description
技术领域
本发明总体上涉及磁电编码器,具体地说,涉及一种用于多对极磁电编码器动态预测窗口角度细分方法。
背景技术
随着工业化的加速发展,对于角位移传感器的分辨率,精度等指标有着较高的要求。目前磁电编码器和光电编码器作为角位移传感器被广泛使用。光电编码器的精度高,但是体积较大,易碎。相对而言,磁电编码器可以在粉尘、油污的环境中工作,但是分辨率难以提高。为了提高磁电编码器的精度,学者做了大量的研究。常用的磁电编码器结构通常包括定子、转子和传感器,转子上固定有永久磁铁,形成磁路系统。永久磁铁随着转子旋转,在旋转过程中形成旋转的磁场信号,磁电编码器信号检测板相对于电机定子处于静止状态,通过信号检测板上的霍尔器件对变化的磁场信号进行采集,通过信号处理解算得到当前转子角度位置,为了提高磁电编码器角度值分辨率,采用单对极磁钢与多对极磁钢组合的方式提高角度值分辨率,多对级磁钢旋转一周产生多周期信号磁场,单对极磁钢旋转一周产生单周期磁场。通过单周期磁场确定当前角度值的绝对位置,多对极产生的角度值信号对单对极信号解算得到的角度值进行角度细分,从而提高角度值分辨率。
然而,磁电式角位移传感器角度值的解算是通过模拟信号采集,然后进行数字转化解算得到的,这种解算过程会引入模数转换噪声,以及硬件电路引入的电源噪声。在电机转子机械位置固定的情况下,解算出的磁电编码器角度值仍然会在一定范围内往复跳动。这种现象会造成磁电编码器过零点位置难以准确判断。对于采用单对极与多对极磁钢组合的磁电编码器而言,转子旋转一周,多对极角度值会多次过零点,过零点时角度值会在0和最大角度值(65536LSB)间往复跳动。
发明内容
针对上述问题,本发明提出一种多对极磁电编码器动态多窗口区间预测角度细分方法,通过动态多窗口角度区间预测方法实现了磁电编码器多对极角度极数的准确判断,实现了磁电编码器角度值的精确细分,为实现高分辨率磁电编码器提供了技术保证。
本发明解决其技术问题的方案是:
一种多对极磁电编码器动态多窗口区间预测角度细分方法,该方法包括以下步骤:
(1)解算单对极及多对极磁电编码器角度值;
(2)依据多对极磁电编码器角度值建立动态多窗口;
(3)依据动态多窗口解算编码器的多对极角度值对应的极数;
(4)依据动态多窗口、单对极角度值、多对极角度值进行多对极磁电编码器角度值细分。
本发明是这样构思的:
步骤一:解算单对极及多对极磁电编码器角度值
以16对极磁电编码器为例进行说明,其中16对极磁钢是32个N、S轴向充磁磁片形成的环形磁钢;为了解算单对极与多对极磁电编码器角度值,通过霍尔信号采集板上的单对极霍尔与多对极霍尔获得单对极d轴霍尔值d1、单对极q轴霍尔值q1、多对极d轴霍尔值d2、多对极q轴霍尔值q2,利用单对极d轴霍尔值d1与单对极q轴霍尔值q1计算单对极角度θ单,利用多对极d轴霍尔值d2与多对极q轴霍尔值q2计算多对极角度θ多,角度计算采用反正切计算方法,其数学模型如公式(1)所示:
步骤二:依据多对极磁电编码器角度值建立动态多窗口
在得到单对极角度值θ单及多对极角度值θ多后,为了确定单对极角度值θ单对应的多对极角度值θ多的极数N,采用动态窗口进行多对极角度值θ多的极数N的判断,动态窗口θ1的数学模型如式(2)所示,其中多对极对应的当前窗口角度值为θ2∈(0~65536LSB):
依据动态窗口θ1的数学模型建立多对极角度值对应的当前窗口角度值θ2、上窗口角度值θ3与下窗口角度值θ4的数学模型如式(3)所示:
步骤三:依据动态多窗口解算编码器的多对极角度值对应的极数
以单对极磁电编码器角度值θ单(0-65536LSB)为横坐标,多对极编码器的当前窗口角度值θ2(0-65536LSB)、上窗口角度值θ3(0-65536LSB)以及下窗口角度值θ4(0-65536LSB)对应的窗口极数映射项为纵坐标,以此确定出每个单对极磁电编码器角度值θ单对应的多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2以及下窗口极数N3的数值,在获得映射关系后,依据多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2以及下窗口极数N3对当前单对极角度值θ单对应的当前窗口角度值θ2的最终的极数N进行判断;根据多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3的大小关系以及当前窗口角度值θ2所处的象限(0-16384LSB,16384-32768LSB,32768-49152LSB,49152-65536LSB)对当前窗口角度值θ2的最终极数N进行计算,其具体过程为:
首先将单对极角度值θ单投影在512(29)个区间内,根据式(2)建立动态多窗口;然后建立多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3与512个单对极区间角度值间的映射关系并制成表格;在实际工作状态下,通过单对极角度值θ单的高9位查找该表格,得到多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3;在获得当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3后执行逻辑判断,如果上窗口极数N2大于当前窗口极数N1,并且当前多对极角度值θ2位于第一象限(0-16384LSB)内,则最终的多对极编码器的当前窗口角度值θ2的极数N为N1+1;如果当前窗口极数N1大于下窗口极数N3,并且多对极编码器的当前窗口角度值θ2位于第四象限(49152-65536LSB)内,则多对极编码器的当前窗口角度值θ2的最终的极数N为N1-1;
步骤四:依据动态多窗口、单对极角度值、多对极角度值进行多对极磁电编码器角度值细分
通过上述方法对多对极编码器的当前窗口角度值θ2所处的最终的极数N进行判断计算,实现多对极编码器角度值θ多在整圈圆周上的细分,细分后的角度值θ5的数学模型如公式(4)所示:
θ5=65536×N+θ2 (4)
本示例为16对极磁电编码器,因此最后细分后的磁电编码器的最大角度分辨率为65535*16=1048576(LSB)。
本发明的有益效果是:
1.采用多对极磁电编码器动态多窗口区间预测角度细分方法,为提高磁电编码器分辨率提供技术基础。
2.动态窗口的建立消除了磁电编码器角度值噪声引发的区间判断错误问题,有效避免了多对极区间误判的问题,提高了磁电编码器的精度。
3.采用多对极磁电编码器动态多窗口区间预测角度细分方法的磁电编码器可以提高角度值分辨率,使得控制系统的增益变大,并且在大增益的情况下不会出现振动,提高控制系统的响应速度。
附图说明
附图1:磁电编码器的结构示意图;
附图2:编码器的d轴与q轴示意图;
附图3:动态窗口示意图;
附图4:编码器的单对极角度值与多对极的当前窗口角度值、上窗口角度值以及下窗口角度值的对应关系示意图;
附图5:查表原理示意图;
附图6:多对极编码器的当前窗口角度值的最终极数逻辑判断流程图;
附图7:细分后的多对极角度值示意图;
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了,下面通过附图中示出的具体实施例来描述本发明。但是应该理解,这些描述只是示例性的,而并非要限制本发明的范围。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本发明的概念。
如图1、图2、图3、图4、图5、图6、图7所示,本具体实施方式采用以下技术方案:一种多对极磁电编码器动态多窗口区间预测角度细分方法,其特征在于:所述方法的具体实现过程为:
1.一种多对极磁电编码器动态多窗口区间预测角度细分方法,其特征在于:所述方法的具体实现过程为:
步骤一:解算单对极及多对极磁电编码器角度值
以16对极磁电编码器为例进行说明,其中16对极磁钢是32个N、S轴向充磁磁片形成的环形磁钢,磁电编码器的结构如图1所示;为了解算单对极与多对极磁电编码器角度值,通过霍尔信号采集板上的单对极霍尔与多对极霍尔获得单对极d轴霍尔值d1、单对极q轴霍尔值q1、多对极d轴霍尔值d2、多对极q轴霍尔值q2,利用单对极d轴霍尔值d1与单对极q轴霍尔值q1计算单对极角度θ单,利用多对极d轴霍尔值d2与多对极q轴霍尔值q2计算多对极角度θ多,其中编码器的d轴与q轴如图2所示;角度计算采用反正切计算方法,其数学模型如公式(1)所示:
步骤二:依据多对极磁电编码器角度值建立动态多窗口
在得到单对极角度值θ单及多对极角度值θ多后,为了确定单对极角度值θ单对应的多对极角度值θ多的极数N,采用动态窗口进行多对极角度值θ多的极数N的判断,动态窗口θ1的数学模型如式(2)所示,其中多对极对应的当前窗口角度值为θ2∈(0~65536LSB):
依据上述公式(1)、(2)建立的动态窗口,如图3所示;
依据动态窗口θ1的数学模型建立多对极角度值对应的当前窗口角度值θ2、上窗口角度值θ3与下窗口角度值θ4的数学模型如式(3)所示:
步骤三:依据动态多窗口解算编码器的多对极角度值对应的极数
如图4所示,以单对极磁电编码器角度值θ单(0-65536LSB)为横坐标,多对极编码器的当前窗口角度值θ2(0-65536LSB)、上窗口角度值θ3(0-65536LSB)以及下窗口角度值θ4(0-65536LSB)对应的窗口极数映射项为纵坐标,以此确定出每个单对极磁电编码器角度值θ单对应的多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2以及下窗口极数N3的数值,在获得映射关系后,依据多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2以及下窗口极数N3对当前单对极角度值θ单对应的当前窗口角度值θ2的最终的极数N进行判断;根据多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3的大小关系以及当前窗口角度值θ2所处的象限(0-16384LSB,16384-32768LSB,32768-49152LSB,49152-65536LSB)对当前窗口角度值θ2的最终极数N进行计算,其具体过程为:
首先将单对极角度值θ单投影在512(29)个区间内,根据式(2)建立动态多窗口;然后建立多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3与512个单对极区间角度值间的映射关系并制成表格;在实际工作状态下,通过单对极角度值θ单的高9位查找该表格,得到多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3,其中查表原理如图5所示;在获得当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3后执行逻辑判断,如果上窗口极数N2大于当前窗口极数N1,并且当前多对极角度值θ2位于第一象限(0-16384LSB)内,则最终的多对极编码器的当前窗口角度值θ2的极数N为N1+1;如果当前窗口极数N1大于下窗口极数N3,并且多对极编码器的当前窗口角度值θ2位于第四象限(49152-65536LSB)内,则最终的多对极编码器的当前窗口角度值θ2的极数N为N1-1,其中多对极编码器的当前窗口角度值θ2的最终极数N逻辑判断流程图如图6所示;
步骤四:依据动态多窗口、单对极角度值、多对极角度值进行多对极磁电编码器角度值细分
通过上述方法对多对极编码器的当前窗口角度值θ2所处的最终极数N进行判断计算,实现多对极编码器角度值θ多在整圈圆周上的细分,细分后的角度值θ5的数学模型如公式(4)所示:
θ5=65536×N+θ2 (4)
本示例为16对极磁电编码器,因此最后细分后的磁电编码器的最大角度分辨率为65535*16=1048576(LSB),此时依据该方法获得的细分后的多对极角度值θ5如图7所示。
以上显示和描述了本发明的基本原理和主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下,本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等效物界定。
Claims (1)
1.一种多对极磁电编码器动态多窗口区间预测角度细分方法,其特征在于:所述方法的具体实现过程为:
步骤一:解算单对极及多对极磁电编码器角度值
16对极磁电编码器,其中16对极磁钢是32个N、S轴向充磁磁片形成的环形磁钢;为了解算单对极与多对极磁电编码器角度值,通过霍尔信号采集板上的单对极霍尔与多对极霍尔获得单对极d轴霍尔值d1、单对极q轴霍尔值q1、多对极d轴霍尔值d2、多对极q轴霍尔值q2,利用单对极d轴霍尔值d1与单对极q轴霍尔值q1计算单对极角度θ单,利用多对极d轴霍尔值d2与多对极q轴霍尔值q2计算多对极角度θ多,角度计算采用反正切计算方法,其数学模型如公式(1)所示:
步骤二:依据多对极磁电编码器角度值建立动态多窗口
在得到单对极角度值θ单及多对极角度值θ多后,为了确定单对极角度值θ单对应的多对极角度值θ多的极数N,采用动态窗口进行多对极角度值θ多的极数N的判断,动态窗口θ1的数学模型如式(2)所示,其中多对极对应的当前窗口角度值为θ2∈(0~65536LSB):
依据动态窗口θ1的数学模型建立多对极角度值对应的当前窗口角度值θ2、上窗口角度值θ3与下窗口角度值θ4的数学模型如式(3)所示:
步骤三:依据动态多窗口解算编码器的多对极角度值对应的极数
以单对极磁电编码器角度值θ单为横坐标,θ单的坐标范围为0-65536LSB,多对极编码器的当前窗口角度值θ2、上窗口角度值θ3以及下窗口角度值θ4对应的窗口极数映射项为纵坐标,θ2的坐标范围为0-65536LSB,θ3的坐标范围为0-65536LSB,θ4的坐标范围为0-65536LSB,以此确定出每个单对极磁电编码器角度值θ单对应的多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2以及下窗口极数N3的数值,在获得映射关系后,依据多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2以及下窗口极数N3对当前单对极角度值θ单对应的当前窗口角度值θ2的最终的极数N进行判断;根据多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3的大小关系以及当前窗口角度值θ2所处的象限0-16384LSB,16384-32768LSB,32768-49152LSB,49152-65536LSB对当前窗口角度值θ2的最终极数N进行计算,其具体过程为:
首先将单对极角度值θ单投影在512个区间内,根据式(2)建立动态多窗口;然后建立多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3与512个单对极区间角度值间的映射关系并制成表格;在实际工作状态下,通过单对极角度值θ单的高9位查找该表格,得到多对极编码器的当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3;在获得当前窗口角度值θ2的当前窗口极数N1、上窗口极数N2、下窗口极数N3后执行逻辑判断,如果上窗口极数N2大于当前窗口极数N1,并且当前多对极角度值θ2位于第一象限0-16384LSB内,则最终的多对极编码器的当前窗口角度值θ2的极数N为N1+1;如果当前窗口极数N1大于下窗口极数N3,并且多对极编码器的当前窗口角度值θ2位于第四象限49152-65536LSB内,则多对极编码器的当前窗口角度值θ2的最终的极数N为N1-1;
步骤四:依据动态多窗口、单对极角度值、多对极角度值进行多对极磁电编码器角度值细分
通过上述方法对多对极编码器的当前窗口角度值θ2所处的最终的极数N进行判断计算,实现多对极编码器角度值θ多在整圈圆周上的细分,细分后的角度值θ5的数学模型如公式(4)所示:
θ5=65536×N+θ2 (4)
16对极磁电编码器,因此最后细分后的磁电编码器的最大角度分辨率为65535*16=1048576LSB。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810424798.0A CN108871181B (zh) | 2018-05-07 | 2018-05-07 | 一种多对极磁电编码器动态多窗口区间预测角度细分方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810424798.0A CN108871181B (zh) | 2018-05-07 | 2018-05-07 | 一种多对极磁电编码器动态多窗口区间预测角度细分方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108871181A CN108871181A (zh) | 2018-11-23 |
CN108871181B true CN108871181B (zh) | 2020-11-10 |
Family
ID=64327310
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810424798.0A Active CN108871181B (zh) | 2018-05-07 | 2018-05-07 | 一种多对极磁电编码器动态多窗口区间预测角度细分方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108871181B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111854807A (zh) * | 2020-05-15 | 2020-10-30 | 哈尔滨理工大学 | 一种关于有域角磁电编码器及其信号逻辑延伸精分方法 |
CN111521839B (zh) * | 2020-06-11 | 2022-10-28 | 哈尔滨理工大学 | 一种基于多对极磁钢及辅助定子绕组测速磁电编码器 |
CN111721329B (zh) * | 2020-07-07 | 2021-11-23 | 哈尔滨理工大学 | 一种三霍尔磁电编码器及免反正切计算角度解算方法 |
CN112067023A (zh) * | 2020-07-26 | 2020-12-11 | 哈尔滨理工大学 | 一种两霍尔磁电编码器及其角度值温漂补偿方法 |
CN112665531B (zh) * | 2020-11-30 | 2023-07-14 | 天津津航技术物理研究所 | 一种多对级旋变坐标变换解角方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE10320990A1 (de) * | 2003-05-09 | 2004-11-25 | Dr. Johannes Heidenhain Gmbh | Induktiver Drehwinkelsensor und damit ausgestatteter Drehgeber |
CN203396406U (zh) * | 2013-08-20 | 2014-01-15 | 杭州日鼎控制技术有限公司 | 编码器的细分装置 |
CN106225813B (zh) * | 2016-07-07 | 2018-08-17 | 航天鑫创自控装备发展有限公司 | 编码器磁钢结构、编码器以及窗口余数区间判断矫正算法 |
CN107607037B (zh) * | 2017-09-18 | 2020-05-19 | 哈尔滨理工大学 | 一种基于反正切跨区间制表法的磁电编码器标定方法 |
-
2018
- 2018-05-07 CN CN201810424798.0A patent/CN108871181B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108871181A (zh) | 2018-11-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108871181B (zh) | 一种多对极磁电编码器动态多窗口区间预测角度细分方法 | |
CN107607037B (zh) | 一种基于反正切跨区间制表法的磁电编码器标定方法 | |
CN108896074B (zh) | 八霍尔分布磁电编码器高分辨率角度值解算方法及装置 | |
CN111551194B (zh) | 基于整数除法判断的多对极磁电编码器角度精分方法及装置 | |
CN109099938B (zh) | 基于极数查表的角度区间扫描角度值跳点抑制方法及装置 | |
CN105007016B (zh) | 一种基于旋转变压器的永磁同步电机测速方法 | |
CN112362089B (zh) | 一种多对极磁电编码器及其高分辨率高可靠角度解算方法 | |
CN103222168B (zh) | 一种伺服电机和伺服控制系统 | |
CN202364092U (zh) | 一种伺服电机和伺服控制系统 | |
CN109059968B (zh) | 基于角度值生长的有限角度值细分方法及装置 | |
CN112117079B (zh) | 编码器磁体结构、编码器、电机、电气设备及车辆 | |
CN103023399A (zh) | 旋转变压器转子零位角度标定系统 | |
CN112113585A (zh) | 编码器及编码器绝对角度的检测方法 | |
CN106655959A (zh) | 一种无刷直流电机转速测量方法 | |
CN111721329A (zh) | 一种三霍尔磁电编码器及免反正切计算角度解算方法 | |
CN106208539A (zh) | 一种磁电式编码器 | |
CN107529384B (zh) | 基于线性霍尔元件的微飞轮测速方法及装置 | |
CN113916265A (zh) | 用于永磁同步直线电机的霍尔位置传感信号的处理方法 | |
CN202940767U (zh) | 旋转变压器转子零位角度标定系统 | |
CN110212818B (zh) | 一种永磁电机角度检测方法 | |
CN105720876B (zh) | 一种永磁同步电机的内功率因数角检测方法 | |
CN106612083A (zh) | 电机转子初始位置和电机转子位置检测方法及装置 | |
CN115441605A (zh) | 一种六相开关磁阻电机及无位置传感器转子位置估算方法和系统 | |
CN111578970B (zh) | 基于核酸序列比对思想的多对极磁电编码器角度精分方法及装置 | |
CN112729100A (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: 20230902 Address after: No. 83, Qunxing 1st Road, Suzhou Industrial Park, Suzhou City, Jiangsu Province, 215000 Patentee after: Suzhou Greenspai Power System Co.,Ltd. Address before: 150080 No. 52, Xuefu Road, Nangang District, Heilongjiang, Harbin Patentee before: HARBIN University OF SCIENCE AND TECHNOLOGY |
|
TR01 | Transfer of patent right |