CN114061632A - 一种用于补偿指定次谐波的高精度磁编码器解码方法 - Google Patents

一种用于补偿指定次谐波的高精度磁编码器解码方法 Download PDF

Info

Publication number
CN114061632A
CN114061632A CN202111228651.2A CN202111228651A CN114061632A CN 114061632 A CN114061632 A CN 114061632A CN 202111228651 A CN202111228651 A CN 202111228651A CN 114061632 A CN114061632 A CN 114061632A
Authority
CN
China
Prior art keywords
signals
harmonic
decoupling
sine
signal
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
CN202111228651.2A
Other languages
English (en)
Other versions
CN114061632B (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.)
Shangda Electric Technology Jiaxing Co ltd
Original Assignee
Shangda Electric Technology Jiaxing 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 Shangda Electric Technology Jiaxing Co ltd filed Critical Shangda Electric Technology Jiaxing Co ltd
Priority to CN202111228651.2A priority Critical patent/CN114061632B/zh
Publication of CN114061632A publication Critical patent/CN114061632A/zh
Application granted granted Critical
Publication of CN114061632B publication Critical patent/CN114061632B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01DMEASURING NOT SPECIALLY ADAPTED FOR A SPECIFIC VARIABLE; ARRANGEMENTS FOR MEASURING TWO OR MORE VARIABLES NOT COVERED IN A SINGLE OTHER SUBCLASS; TARIFF METERING APPARATUS; MEASURING OR TESTING NOT OTHERWISE PROVIDED FOR
    • G01D5/00Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable
    • G01D5/12Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable using electric or magnetic means
    • G01D5/244Mechanical means for transferring the output of a sensing member; Means for converting the output of a sensing member to another variable where the form or nature of the sensing member does not constrain the means for converting; Transducers not specially adapted for a specific variable using electric or magnetic means influencing characteristics of pulses or pulse trains; generating pulses or pulse trains
    • G01D5/24471Error correction
    • G01D5/24476Signal processing
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/0003Control strategies in general, e.g. linear type, e.g. P, PI, PID, using robust control
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P21/00Arrangements or methods for the control of electric machines by vector control, e.g. by control of field orientation
    • H02P21/14Estimation or adaptation of machine parameters, e.g. flux, current or voltage
    • H02P21/18Estimation of position or speed
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P25/00Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details
    • H02P25/02Arrangements or methods for the control of AC motors characterised by the kind of AC motor or by structural details characterised by the kind of motor
    • H02P25/022Synchronous motors
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02PCONTROL OR REGULATION OF ELECTRIC MOTORS, ELECTRIC GENERATORS OR DYNAMO-ELECTRIC CONVERTERS; CONTROLLING TRANSFORMERS, REACTORS OR CHOKE COILS
    • H02P2207/00Indexing scheme relating to controlling arrangements characterised by the type of motor
    • H02P2207/05Synchronous machines, e.g. with permanent magnets or DC excitation
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E40/00Technologies for an efficient electrical power generation, transmission or distribution
    • Y02E40/40Arrangements for reducing harmonics

Abstract

本发明公开了一种用于补偿指定次谐波的高精度磁编码器解码方法,包括,利用双同步旋转坐标系将磁编码器的正余弦信号vC和vS进行双同步坐标变换,得到四个正、负序分量;通过正反向解耦网络对四个正、负序分量进行解耦,获得解耦分量,并通过低通滤波器获得各解耦分量的模值;利用各解耦分量的模值计算基波信号的幅值和相位,并结合锁相环输出角度重构基波信号vS_1和vC_1;根据两路正余弦信号和基波信号获得包含各次谐波的高频信号vS_har和vC_har,并先后对其进行双同步坐标变换、解耦和滤波操作,获得指定次数的谐波信号vS_n和vC_n;将指定次数的谐波信号反馈到两路正余弦信号中并减去,形成闭环;本发明实现了对角度的解算,同时提高了角度结算精度,能够保证系统的动静态性能。

Description

一种用于补偿指定次谐波的高精度磁编码器解码方法
技术领域
本发明涉及磁编码器解码的技术领域,尤其涉及一种用于补偿指定次谐波的高精度磁编码器解码方法。
背景技术
永磁同步电机(Permanent-Magnet Synchronous Motor,PMSM)具有功率密度高、结构简单、噪声小、效率高等特点,广泛应用于航空航天、国防、工农业生产和日常生活的各个领域。
在永磁同步电机控制系统中,需要获得精确的转子位置信息来实现高动态性能控制。获得转子位置信息的一般方法是通过机械式位置传感器直接检测,其中包括霍尔传感器、光电编码器、旋转变压器、磁编码器等。位置传感器作为伺服系统中的反馈驱动模块,其性能好坏和灵敏度直接决定着整个伺服驱动系统的性能、精度和分辨率。这几种传感器中磁编码器具有响应速度快,抗震能力强,环境适应能力强,成本低,具有较高的分辨率等优点,但是精度相对较低。而在物联网发展的推进下,未来对于磁编码器发展方向必然是高精尖。在高精度、高动态性能要求的伺服系统中,必须实时和精确地测量转子的位置和转速,因此提高传感器的性能显得尤为重要。
磁编码器的信号解码方法可分为硬件解码和软件解码两种。硬件解码包括使用分立器件搭建解码电路和使用专用解码芯片完成解码。磁编码器常见的软件解码算法包括反正切法,CORDIC算法,查表法和锁相环跟踪算法。为了提高结算精度,一般在解码前先对信号进行处理,考虑各种误差补偿。现有技术大多仅考虑幅值偏差、相位偏差和直流偏差,对于谐波误差的考虑较少,难以对幅值偏差、相位偏差和谐波误差进行消除。
发明内容
本部分的目的在于概述本发明的实施例的一些方面以及简要介绍一些较佳实施例。在本部分以及本申请的说明书摘要和发明名称中可能会做些简化或省略以避免使本部分、说明书摘要和发明名称的目的模糊,而这种简化或省略不能用于限制本发明的范围。
鉴于上述现有存在的问题,提出了本发明。
因此,本发明提供了一种用于补偿指定次谐波的高精度磁编码器解码方法,能够消除指定次数谐波的交叉反馈网络,实现对幅值偏差、相位偏差和谐波误差进行消除,提高了解算位置信息精度。
为解决上述技术问题,本发明提供如下技术方案:包括,利用双同步旋转坐标系将磁编码器的正余弦信号vC和vS进行双同步坐标变换,得到d、q轴上的四个正、负序分量;通过正反向解耦网络对所述四个正、负序分量进行解耦,获得解耦分量
Figure BDA0003315155970000021
Figure BDA0003315155970000022
并通过低通滤波器获得各解耦分量的模值
Figure BDA0003315155970000023
Figure BDA0003315155970000024
利用所述各解耦分量的模值
Figure BDA0003315155970000025
Figure BDA0003315155970000026
计算基波信号的幅值和相位,并结合锁相环输出角度重构基波信号vS_1和vC_1;根据两路正余弦信号和所述基波信号获得包含各次谐波的高频信号vS_har和vC_har,并先后对其进行双同步坐标变换、解耦和滤波操作,获得指定次数的谐波信号vS_n和vC_n;将所述指定次数的谐波信号vS_n和vC_n反馈到两路正余弦信号中并减去,形成闭环。
作为本发明所述的用于补偿指定次谐波的高精度磁编码器解码方法的一种优选方案,其中:所述双同步坐标变换包括,所述磁编码器的正余弦信号vC和vS为:
Figure BDA0003315155970000027
将所述磁编码器的正余弦信号vC和vS合成为一个幅值变化的旋转矢量vA,利用双同步旋转坐标系将所述旋转矢量vA分解到正反向旋转的两个坐标系上,得到d、q轴上的四个正、负序分量:
Figure BDA0003315155970000028
其中,V1c和V1s为正余弦信号的基波幅值,σ1c和σ1s为正余弦信号的初始相位,θ表示理想角度,即锁相环输出角度;Vnc和Vns为正余弦信号中各次谐波的幅值,σnc和σns为正余弦信号中各次谐波的初始相位,n为谐波次数;
Figure BDA0003315155970000029
Figure BDA00033151559700000210
表示分解矢量
Figure BDA00033151559700000211
Figure BDA00033151559700000212
分解到正负向旋转坐标系下的d、q轴分量,
Figure BDA00033151559700000213
Figure BDA00033151559700000214
表示
Figure BDA00033151559700000215
Figure BDA00033151559700000216
两个矢量的幅值。
作为本发明所述的用于补偿指定次谐波的高精度磁编码器解码方法的一种优选方案,其中:所述解耦包括,
Figure BDA0003315155970000031
作为本发明所述的用于补偿指定次谐波的高精度磁编码器解码方法的一种优选方案,其中:计算基波信号的幅值和相位包括,对所述各解耦分量的模值
Figure BDA0003315155970000032
Figure BDA0003315155970000033
求平方和,再通过开方运算进而获得基波信号的幅值
Figure BDA0003315155970000034
将解耦分量的模值
Figure BDA0003315155970000035
除以解耦分量的模值
Figure BDA0003315155970000036
再通过反正切函数获得基波信号的相位
Figure BDA0003315155970000037
Figure BDA0003315155970000038
作为本发明所述的用于补偿指定次谐波的高精度磁编码器解码方法的一种优选方案,其中:所述重构基波信号vS_1和vC_1包括,
Figure BDA0003315155970000039
作为本发明所述的用于补偿指定次谐波的高精度磁编码器解码方法的一种优选方案,其中:所述高频信号vS_har和vC_har包括,将所述两路正余弦信号和所述基波信号进行相减,获得包含各次谐波的高频信号vS_har和vC_har
Figure BDA00033151559700000310
作为本发明所述的用于补偿指定次谐波的高精度磁编码器解码方法的一种优选方案,其中:所述指定次数的谐波信号vS_n和vC_n包括,先后对高频信号vS_har和vC_har进行双同步坐标变换、解耦和滤波操作,所需角度θn通过将锁相环输出角度θ乘上所要补偿谐波的次数n得到,获得指定次数的谐波信号vS_n和vC_n
Figure BDA0003315155970000041
其中,
Figure BDA0003315155970000042
表示指定次数的谐波正负序信号的幅值;
Figure BDA0003315155970000043
表示指定次数的谐波正负序信号的相位。
作为本发明所述的用于补偿指定次谐波的高精度磁编码器解码方法的一种优选方案,其中:还包括,将所述指定次数的谐波信号vS_n和vC_n反馈到两路正余弦信号中并减去,获得指定次数的谐波补偿信号v'S和v'C,而后对所述指定次数的谐波补偿信号v'S和v'C进行双同步坐标变换、解耦和滤波操作,形成闭环;通过不断反馈,逐步衰减指定次数的谐波信号;其中,指定次数的谐波补偿信号v'S和v'C为:
Figure BDA0003315155970000044
本发明的有益效果:本发明实现了对角度的解算,同时考虑了基波、谐波的幅值和相位,并对其进行了准确的信号重构,对于指定次谐波能够较好的衰减,进而提高了角度结算精度,同时还能保证系统的动静态性能。
附图说明
为了更清楚地说明本发明实施例的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其它的附图。其中:
图1为本发明第一个实施例所述的用于补偿指定次谐波的高精度磁编码器解码方法的整体结构框图;
图2为本发明第一个实施例所述的用于补偿指定次谐波的高精度磁编码器解码方法的双同步旋转坐标系;
图3为本发明第一个实施例所述的用于补偿指定次谐波的高精度磁编码器解码方法的正负向解耦网络的结构框图;
图4为本发明第一个实施例所述的用于补偿指定次谐波的高精度磁编码器解码方法的信号重构的结构框图;
图5为本发明第二个实施例所述的用于补偿指定次谐波的高精度磁编码器解码方法的转子位置同实际转子位置误差波形与传统方法误差波形的对比示意图。
具体实施方式
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合说明书附图对本发明的具体实施方式做详细的说明,显然所描述的实施例是本发明的一部分实施例,而不是全部实施例。基于本发明中的实施例,本领域普通人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明的保护的范围。
在下面的描述中阐述了很多具体细节以便于充分理解本发明,但是本发明还可以采用其他不同于在此描述的其它方式来实施,本领域技术人员可以在不违背本发明内涵的情况下做类似推广,因此本发明不受下面公开的具体实施例的限制。
其次,此处所称的“一个实施例”或“实施例”是指可包含于本发明至少一个实现方式中的特定特征、结构或特性。在本说明书中不同地方出现的“在一个实施例中”并非均指同一个实施例,也不是单独的或选择性的与其他实施例互相排斥的实施例。
本发明结合示意图进行详细描述,在详述本发明实施例时,为便于说明,表示器件结构的剖面图会不依一般比例作局部放大,而且所述示意图只是示例,其在此不应限制本发明保护的范围。此外,在实际制作中应包含长度、宽度及深度的三维空间尺寸。
同时在本发明的描述中,需要说明的是,术语中的“上、下、内和外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本发明和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本发明的限制。此外,术语“第一、第二或第三”仅用于描述目的,而不能理解为指示或暗示相对重要性。
本发明中除非另有明确的规定和限定,术语“安装、相连、连接”应做广义理解,例如:可以是固定连接、可拆卸连接或一体式连接;同样可以是机械连接、电连接或直接连接,也可以通过中间媒介间接相连,也可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1
参照图1~4,为本发明的第一个实施例,该实施例提供了一种用于补偿指定次谐波的高精度磁编码器解码方法,包括:
S1:利用双同步旋转坐标系将磁编码器的正余弦信号vC和vS进行双同步坐标变换,得到d、q轴上的四个正、负序分量。
磁编码器的正余弦信号vC和vS为:
Figure BDA0003315155970000061
其中,V1c和V1s为正余弦信号的基波幅值,σ1c和σ1s为正余弦信号的初始相位,θ表示理想角度,即锁相环输出角度;Vnc和Vns为正余弦信号中各次谐波的幅值,σnc和σns为正余弦信号中各次谐波的初始相位,n为谐波次数。
将磁编码器的正余弦信号vC和vS合成为一个幅值变化的旋转矢量vA,利用双同步旋转坐标系将旋转矢量vA分解到正反向旋转的两个坐标系上,得到两个幅值恒定但旋转方向相反的矢量
Figure BDA0003315155970000062
Figure BDA0003315155970000063
参照图2,它包含两个旋转方向相反的坐标系,其中,dq+1表示以角速度
Figure BDA0003315155970000064
正向旋转的坐标系,dq-1表示以角速度
Figure BDA0003315155970000065
反向旋转的坐标系,vA为不对称分量,θ'为估计角度,θ为理想角度;根据图2,可以得到d、q轴上的四个正、负序分量:
Figure BDA0003315155970000066
Figure BDA0003315155970000067
其中,
Figure BDA0003315155970000068
Figure BDA0003315155970000069
表示分解矢量
Figure BDA00033151559700000610
Figure BDA00033151559700000611
分解到正负向旋转坐标系下的d、q轴分量,
Figure BDA00033151559700000612
Figure BDA00033151559700000613
表示
Figure BDA00033151559700000614
Figure BDA00033151559700000615
两个矢量的幅值;上标T表示转置。
从上式可以看出,旋转变压器输出矢量vA在正向旋转坐标轴上的投影包含由反向旋转矢量耦合产生的波动分量,在反向坐标轴上的投影与之相反;为了消除正负坐标旋转耦合产生的波动分量,需要使用正反向解耦网络进行解耦。S2:通过正反向解耦网络对四个正、负序分量进行解耦,获得解耦分量
Figure BDA0003315155970000071
Figure BDA0003315155970000072
Figure BDA0003315155970000073
并通过低通滤波器获得各解耦分量的模值
Figure BDA0003315155970000074
Figure BDA0003315155970000075
参照图3,图3为正负向解耦网络的结构框图,正负向解耦网络的作用是消除正负坐标旋转耦合产生的波动分量;其中,
Figure BDA0003315155970000076
Figure BDA0003315155970000077
表示指定次数谐波分解到正负向旋转坐标系下的d、q轴分量,
Figure BDA0003315155970000078
Figure BDA0003315155970000079
表示指定信号正负序分量的模值,
Figure BDA00033151559700000710
Figure BDA00033151559700000711
为通过正反向解耦网络后的解耦正负序分量,即解耦分量,nθ'表示估计角度乘上指定次数。
具体的,解耦的具体过程为:将正向旋转矢量的q轴解耦分量
Figure BDA00033151559700000712
通过锁相环求解角度,即通过一个比例积分环节,得到估计转速
Figure BDA00033151559700000713
再通过一个积分环节即可得到估计角度θ',当锁相环达到稳态,即
Figure BDA00033151559700000714
为零时,此时输出角度θ'与实际转子位置θ达到同步状态,即θ'=θ,则:
Figure BDA00033151559700000715
在双同步坐标系中,将上式展开为:
Figure BDA00033151559700000716
通过解耦将反向旋转矢量耦合产生的频率为2ω的波动分量消除,解耦思路即将波动分量重构出来并减去;其中,
Figure BDA00033151559700000717
Figure BDA00033151559700000718
四个恒定分量可以将
Figure BDA00033151559700000719
Figure BDA00033151559700000720
通过低通滤波器得到,角度由锁相环求得;而四个恒定分量又可以反馈到反向解耦的环节中,形成一个解耦网络,如图3所示,消除正负坐标旋转耦合产生的波动分量得到获得解耦分量
Figure BDA0003315155970000081
Figure BDA0003315155970000082
Figure BDA0003315155970000083
较佳的是,将解耦后的
Figure BDA0003315155970000084
送入锁相环,而锁相环输出角度有反馈回双同步旋转变换和解耦环节中形成闭环,实现双同步旋转坐标系锁相环求解角度,且在消除了频率为2ω的波动分量后,锁相环输出角度更加准确。
S3:利用各解耦分量的模值
Figure BDA0003315155970000085
Figure BDA0003315155970000086
计算基波信号的幅值和相位,并结合锁相环输出角度重构基波信号vS_1和vC_1
对各解耦分量的模值
Figure BDA0003315155970000087
Figure BDA0003315155970000088
求平方和,再通过开方运算进而获得基波信号的幅值
Figure BDA0003315155970000089
将解耦分量的模值
Figure BDA00033151559700000810
除以解耦分量的模值
Figure BDA00033151559700000811
再通过反正切函数获得基波信号的相位
Figure BDA00033151559700000812
Figure BDA00033151559700000813
将解算出的角度θ'导入,这样即可近似重构出基波信号vS_1和vC_1
Figure BDA00033151559700000814
S4:根据两路正余弦信号和基波信号获得包含各次谐波的高频信号vS_har和vC_har,并先后对其进行双同步坐标变换、解耦和滤波操作,获得指定次数的谐波信号vS_n和vC_n
将两路正余弦信号和基波信号进行相减,获得包含各次谐波的高频信号vS_har和vC_har
Figure BDA0003315155970000091
进一步的,先后对高频信号vS_har和vC_har进行双同步坐标变换、解耦和滤波操作,所需角度θn通过将锁相环输出角度θ乘上所要补偿谐波的次数n得到,获得指定次数的谐波信号vS_n和vC_n
具体的,(1)对高频信号vS_har和vC_har进行过双同步坐标变换,对于指定次数谐波含量同样可以合成为一个幅值变化的旋转矢量vA n,在双旋转坐标系下分解得到两个幅值恒定但旋转方向相反的矢量
Figure BDA0003315155970000092
Figure BDA0003315155970000093
得到:
Figure BDA0003315155970000094
其中,
Figure BDA0003315155970000095
Figure BDA0003315155970000096
表示
Figure BDA0003315155970000097
Figure BDA0003315155970000098
分解到正负向旋转坐标系下的d、q轴分量,
Figure BDA0003315155970000099
Figure BDA00033151559700000910
表示
Figure BDA00033151559700000911
Figure BDA00033151559700000912
两个矢量的幅值,对于非指定次数谐波上式中不包含直流分量,在后续滤波环节不会有输出,能够实现针对某一次谐波进行处理。
(2)将指定次数谐波分解的四个正负序分量
Figure BDA00033151559700000913
Figure BDA00033151559700000914
通过正负向解耦网络,并同样通过低通滤波器可以得到四个的恒定分量
Figure BDA00033151559700000915
Figure BDA00033151559700000916
即可得到指定次谐波的正负序信号的幅值
Figure BDA00033151559700000917
和相位
Figure BDA00033151559700000918
Figure BDA00033151559700000919
进而可重构出指定次数的谐波信号vS_n和vC_n
Figure BDA00033151559700000920
S5:将指定次数的谐波信号vS_n和vC_n反馈到两路正余弦信号中并减去,形成闭环。
将指定次数的谐波信号vS_n和vC_n反馈到两路正余弦信号中并减去,获得指定次数的谐波补偿信号v'S和v'C,而后对指定次数的谐波补偿信号v'S和v'C进行双同步坐标变换、解耦和滤波操作,形成闭环;
通过不断反馈,逐步衰减指定次数的谐波信号;
其中,指定次数的谐波补偿信号v'S和v'C为:
Figure BDA0003315155970000101
实施例2
为了对本方法中采用的技术效果加以验证说明,本实施例选择传统无谐波补偿的双同步旋转坐标系锁相环角度跟踪算法与本方法进行仿真对比,以验证本方法所具有的真实效果。
根据下面仿真实例对本发明所设计的方法进行验证。
本实施例以一组包含2、3、4、5次谐波的正余弦信号为仿真对象,在MATLAB中进行仿真实验,实验系统中设定仿真步长Ts=1e-5s,基波信号频率为50hz,正弦信号幅值为1.011,相位偏移为1°,余弦信号幅值0.991,相位偏移-0.5°;二次谐波正弦分量幅值为0.011,相位为偏移2°,余弦分量幅值0.009,相位偏移1°;三次谐波正弦分量幅值为0.012,相位偏移为1°,余弦分量幅值0.008,相位偏移-1°;四次谐波正弦分量幅值为0.011,相位偏移为1°,余弦分量幅值0.009,相位偏移2°;五次谐波正弦分量幅值为0.009,相位偏移为1°,余弦分量幅值0.011,相位偏移-2°。
提取恒定直流分量的巴特沃斯低通滤波器阶数为2,截止频率设为20Hz;首先用传统的双同步旋转坐标系锁相环角度跟踪算法进行结算,在0.3s时将谐波补偿环节加入,在0.6s时将频率突变为70hz;其中,补偿环节对二次和三次谐波进行补偿,由MATLAB的反正切环节提供理想角度作为参考,实验得到改进前后两种方法的角度误差结果如图5所示。
从图5可见,角度误差稳定在-0.015-0.023rad之间,在0.3s加入反馈环节,反馈环节对系统的动态性能无影响,而且能够降低角度误差,在-0.009-0.008rad之间;在0.6s突然加速,频率上升为70hz,可以看出转速突变对估算角度影响不大,仅在突加速瞬间误差稍增大,角度误差在0.09rad左右,估计角度能快速收敛于实际角度,恢复到稳定状态,说明本方法具有良好的抗转速扰动能力,能够减小转子位置误差波动,位置观测更加准确,同时也能提高系统的动静态性能。
应当认识到,本发明的实施例可以由计算机硬件、硬件和软件的组合、或者通过存储在非暂时性计算机可读存储器中的计算机指令来实现或实施。所述方法可以使用标准编程技术-包括配置有计算机程序的非暂时性计算机可读存储介质在计算机程序中实现,其中如此配置的存储介质使得计算机以特定和预定义的方式操作——根据在具体实施例中描述的方法和附图。每个程序可以以高级过程或面向对象的编程语言来实现以与计算机系统通信。然而,若需要,该程序可以以汇编或机器语言实现。在任何情况下,该语言可以是编译或解释的语言。此外,为此目的该程序能够在编程的专用集成电路上运行。
此外,可按任何合适的顺序来执行本文描述的过程的操作,除非本文另外指示或以其他方式明显地与上下文矛盾。本文描述的过程(或变型和/或其组合)可在配置有可执行指令的一个或多个计算机系统的控制下执行,并且可作为共同地在一个或多个处理器上执行的代码(例如,可执行指令、一个或多个计算机程序或一个或多个应用)、由硬件或其组合来实现。所述计算机程序包括可由一个或多个处理器执行的多个指令。
进一步,所述方法可以在可操作地连接至合适的任何类型的计算平台中实现,包括但不限于个人电脑、迷你计算机、主框架、工作站、网络或分布式计算环境、单独的或集成的计算机平台、或者与带电粒子工具或其它成像装置通信等等。本发明的各方面可以以存储在非暂时性存储介质或设备上的机器可读代码来实现,无论是可移动的还是集成至计算平台,如硬盘、光学读取和/或写入存储介质、RAM、ROM等,使得其可由可编程计算机读取,当存储介质或设备由计算机读取时可用于配置和操作计算机以执行在此所描述的过程。此外,机器可读代码,或其部分可以通过有线或无线网络传输。当此类媒体包括结合微处理器或其他数据处理器实现上文所述步骤的指令或程序时,本文所述的发明包括这些和其他不同类型的非暂时性计算机可读存储介质。当根据本发明所述的方法和技术编程时,本发明还包括计算机本身。计算机程序能够应用于输入数据以执行本文所述的功能,从而转换输入数据以生成存储至非易失性存储器的输出数据。输出信息还可以应用于一个或多个输出设备如显示器。在本发明优选的实施例中,转换的数据表示物理和有形的对象,包括显示器上产生的物理和有形对象的特定视觉描绘。
如在本申请所使用的,术语“组件”、“模块”、“系统”等等旨在指代计算机相关实体,该计算机相关实体可以是硬件、固件、硬件和软件的结合、软件或者运行中的软件。例如,组件可以是,但不限于是:在处理器上运行的处理、处理器、对象、可执行文件、执行中的线程、程序和/或计算机。作为示例,在计算设备上运行的应用和该计算设备都可以是组件。一个或多个组件可以存在于执行中的过程和/或线程中,并且组件可以位于一个计算机中以及/或者分布在两个或更多个计算机之间。此外,这些组件能够从在其上具有各种数据结构的各种计算机可读介质中执行。这些组件可以通过诸如根据具有一个或多个数据分组(例如,来自一个组件的数据,该组件与本地系统、分布式系统中的另一个组件进行交互和/或以信号的方式通过诸如互联网之类的网络与其它系统进行交互)的信号,以本地和/或远程过程的方式进行通信。
应说明的是,以上实施例仅用以说明本发明的技术方案而非限制,尽管参照较佳实施例对本发明进行了详细说明,本领域的普通技术人员应当理解,可以对本发明的技术方案进行修改或者等同替换,而不脱离本发明技术方案的精神和范围,其均应涵盖在本发明的权利要求范围当中。

Claims (8)

1.一种用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:包括,
利用双同步旋转坐标系将磁编码器的正余弦信号vC和vS进行双同步坐标变换,得到d、q轴上的四个正、负序分量;
通过正反向解耦网络对所述四个正、负序分量进行解耦,获得解耦分量
Figure FDA0003315155960000011
Figure FDA0003315155960000012
Figure FDA0003315155960000013
并通过低通滤波器获得各解耦分量的模值
Figure FDA0003315155960000014
Figure FDA0003315155960000015
利用所述各解耦分量的模值
Figure FDA0003315155960000016
Figure FDA0003315155960000017
计算基波信号的幅值和相位,并结合锁相环输出角度重构基波信号vS_1和vC_1
根据两路正余弦信号和所述基波信号获得包含各次谐波的高频信号vS_har和vC_har,并先后对其进行双同步坐标变换、解耦和滤波操作,获得指定次数的谐波信号vS_n和vC_n
将所述指定次数的谐波信号vS_n和vC_n反馈到两路正余弦信号中并减去,形成闭环。
2.如权利要求1所述的用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:所述双同步坐标变换包括,
所述磁编码器的正余弦信号vC和vS为:
Figure FDA0003315155960000018
将所述磁编码器的正余弦信号vC和vS合成为一个幅值变化的旋转矢量vA,利用双同步旋转坐标系将所述旋转矢量vA分解到正反向旋转的两个坐标系上,得到d、q轴上的四个正、负序分量:
Figure FDA0003315155960000019
其中,V1c和V1s为正余弦信号的基波幅值,σ1c和σ1s为正余弦信号的初始相位,θ表示理想角度,即锁相环输出角度;Vnc和Vns为正余弦信号中各次谐波的幅值,σnc和σns为正余弦信号中各次谐波的初始相位,n为谐波次数;
Figure FDA0003315155960000021
Figure FDA0003315155960000022
表示分解矢量
Figure FDA0003315155960000023
Figure FDA0003315155960000024
分解到正负向旋转坐标系下的d、q轴分量,
Figure FDA0003315155960000025
Figure FDA0003315155960000026
表示
Figure FDA0003315155960000027
Figure FDA0003315155960000028
两个矢量的幅值。
3.如权利要求2所述的用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:所述解耦包括,
Figure FDA0003315155960000029
4.如权利要求2或3所述的用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:计算基波信号的幅值和相位包括,
对所述各解耦分量的模值
Figure FDA00033151559600000210
Figure FDA00033151559600000211
求平方和,再通过开方运算进而获得基波信号的幅值
Figure FDA00033151559600000212
将解耦分量的模值
Figure FDA00033151559600000213
除以解耦分量的模值
Figure FDA00033151559600000214
再通过反正切函数获得基波信号的相位
Figure FDA00033151559600000215
Figure FDA00033151559600000216
5.如权利要求4所述的用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:所述重构基波信号vS_1和vC_1包括,
Figure FDA00033151559600000217
6.如权利要求5所述的用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:所述高频信号vS_har和vC_har包括,
将所述两路正余弦信号和所述基波信号进行相减,获得包含各次谐波的高频信号vS_har和vC_har
Figure FDA00033151559600000218
7.如权利要求6所述的用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:所述指定次数的谐波信号vS_n和vC_n包括,
先后对高频信号vS_har和vC_har进行双同步坐标变换、解耦和滤波操作,所需角度θn通过将锁相环输出角度θ乘上所要补偿谐波的次数n得到,获得指定次数的谐波信号vS_n和vC_n
Figure FDA0003315155960000031
其中,
Figure FDA0003315155960000032
表示指定次数的谐波正负序信号的幅值;
Figure FDA0003315155960000033
表示指定次数的谐波正负序信号的相位。
8.如权利要求7所述的用于补偿指定次谐波的高精度磁编码器解码方法,其特征在于:还包括,
将所述指定次数的谐波信号vS_n和vC_n反馈到两路正余弦信号中并减去,获得指定次数的谐波补偿信号v'S和v'C,而后对所述指定次数的谐波补偿信号v'S和v'C进行双同步坐标变换、解耦和滤波操作,形成闭环;
通过不断反馈,逐步衰减指定次数的谐波信号;
其中,指定次数的谐波补偿信号v'S和v'C为:
Figure FDA0003315155960000034
CN202111228651.2A 2021-10-21 2021-10-21 一种用于补偿指定次谐波的高精度磁编码器解码方法 Active CN114061632B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111228651.2A CN114061632B (zh) 2021-10-21 2021-10-21 一种用于补偿指定次谐波的高精度磁编码器解码方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111228651.2A CN114061632B (zh) 2021-10-21 2021-10-21 一种用于补偿指定次谐波的高精度磁编码器解码方法

Publications (2)

Publication Number Publication Date
CN114061632A true CN114061632A (zh) 2022-02-18
CN114061632B CN114061632B (zh) 2024-03-19

Family

ID=80235213

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111228651.2A Active CN114061632B (zh) 2021-10-21 2021-10-21 一种用于补偿指定次谐波的高精度磁编码器解码方法

Country Status (1)

Country Link
CN (1) CN114061632B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117516596A (zh) * 2024-01-05 2024-02-06 华中科技大学 一种磁电编码器的高精度在线补偿方法

Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11103527A (ja) * 1997-09-29 1999-04-13 Tokyo Electric Power Co Inc:The 高調波補償方式
US20060076480A1 (en) * 2004-10-13 2006-04-13 Mitutoyo Corporation Encoder output signal correction apparatus and method
US20060077083A1 (en) * 2004-10-13 2006-04-13 Mitutoyo Corporation Encoder output signal correction apparatus and method
US20090105979A1 (en) * 2007-10-23 2009-04-23 Gerardo Escobar Valderrama Fixed reference frame phase-locked loop (FRF-PLL) for unbalanced line voltage conditions
CN102095915A (zh) * 2010-11-24 2011-06-15 重庆大学 一种采用多同步参考坐标系变换的电压信号检测装置
CN103257271A (zh) * 2013-05-16 2013-08-21 南京工程学院 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法
CN103776471A (zh) * 2014-01-28 2014-05-07 华中科技大学 一种基于双同步旋转坐标系的磁编码器
CN105162137A (zh) * 2015-09-12 2015-12-16 南昌航空大学 复杂电网下无功及谐波电流检测算法
CN106053918A (zh) * 2016-05-19 2016-10-26 重庆大学 多通道的电网电压同步信号提取装置及方法
US20170284826A1 (en) * 2016-03-31 2017-10-05 Delta Tau Data Systems, Inc. Method and Apparatus for Encoder Interpolation
CN108599261A (zh) * 2018-04-24 2018-09-28 上海电力学院 基于非线性pi和解耦双同步坐标系锁相环的锁相方法
CN108631656A (zh) * 2018-05-30 2018-10-09 阳光电源股份有限公司 一种旋转变压器解码方法及系统
CN109579880A (zh) * 2018-12-26 2019-04-05 上海英威腾工业技术有限公司 带有自适应补偿功能的磁编码器
CN110061615A (zh) * 2019-06-03 2019-07-26 上海理工大学 逆变器非线性特性的定子电流谐波补偿方法
US20190245458A1 (en) * 2018-02-03 2019-08-08 S&C Electric Company Three phase inverter grid voltage feed forward method under unbalanced operating conditions
US20190379205A1 (en) * 2017-02-27 2019-12-12 Magnetic Pumping Solutions, Llc Sogi-based pll for grid connection and motor control
CN110768666A (zh) * 2019-10-28 2020-02-07 南京工程学院 基于卡尔曼滤波器的双同步坐标系解耦的锁相环系统及方法
CN111245304A (zh) * 2018-11-29 2020-06-05 安徽美芝精密制造有限公司 补偿方法、补偿装置、电机和存储介质
CN112698095A (zh) * 2021-01-22 2021-04-23 中南大学 一种基波和任意次谐波的高精度检测方法
CN112737450A (zh) * 2020-12-24 2021-04-30 上海大学 一种用于spmsm转子位置估计的高频注入补偿方法
CN112829602A (zh) * 2021-01-04 2021-05-25 宝能(西安)汽车研究院有限公司 车辆扭矩控制方法、装置以及车辆
CN113131820A (zh) * 2020-01-13 2021-07-16 中国科学院电工研究所 一种永磁同步电机周期性误差的分析方法及装置

Patent Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH11103527A (ja) * 1997-09-29 1999-04-13 Tokyo Electric Power Co Inc:The 高調波補償方式
US20060076480A1 (en) * 2004-10-13 2006-04-13 Mitutoyo Corporation Encoder output signal correction apparatus and method
US20060077083A1 (en) * 2004-10-13 2006-04-13 Mitutoyo Corporation Encoder output signal correction apparatus and method
US20090105979A1 (en) * 2007-10-23 2009-04-23 Gerardo Escobar Valderrama Fixed reference frame phase-locked loop (FRF-PLL) for unbalanced line voltage conditions
CN102095915A (zh) * 2010-11-24 2011-06-15 重庆大学 一种采用多同步参考坐标系变换的电压信号检测装置
CN103257271A (zh) * 2013-05-16 2013-08-21 南京工程学院 一种基于stm32f107vct6的微电网谐波与间谐波检测装置及检测方法
CN103776471A (zh) * 2014-01-28 2014-05-07 华中科技大学 一种基于双同步旋转坐标系的磁编码器
CN105162137A (zh) * 2015-09-12 2015-12-16 南昌航空大学 复杂电网下无功及谐波电流检测算法
US20170284826A1 (en) * 2016-03-31 2017-10-05 Delta Tau Data Systems, Inc. Method and Apparatus for Encoder Interpolation
CN106053918A (zh) * 2016-05-19 2016-10-26 重庆大学 多通道的电网电压同步信号提取装置及方法
US20190379205A1 (en) * 2017-02-27 2019-12-12 Magnetic Pumping Solutions, Llc Sogi-based pll for grid connection and motor control
US20190245458A1 (en) * 2018-02-03 2019-08-08 S&C Electric Company Three phase inverter grid voltage feed forward method under unbalanced operating conditions
CN108599261A (zh) * 2018-04-24 2018-09-28 上海电力学院 基于非线性pi和解耦双同步坐标系锁相环的锁相方法
CN108631656A (zh) * 2018-05-30 2018-10-09 阳光电源股份有限公司 一种旋转变压器解码方法及系统
CN111245304A (zh) * 2018-11-29 2020-06-05 安徽美芝精密制造有限公司 补偿方法、补偿装置、电机和存储介质
CN109579880A (zh) * 2018-12-26 2019-04-05 上海英威腾工业技术有限公司 带有自适应补偿功能的磁编码器
CN110061615A (zh) * 2019-06-03 2019-07-26 上海理工大学 逆变器非线性特性的定子电流谐波补偿方法
CN110768666A (zh) * 2019-10-28 2020-02-07 南京工程学院 基于卡尔曼滤波器的双同步坐标系解耦的锁相环系统及方法
CN113131820A (zh) * 2020-01-13 2021-07-16 中国科学院电工研究所 一种永磁同步电机周期性误差的分析方法及装置
CN112737450A (zh) * 2020-12-24 2021-04-30 上海大学 一种用于spmsm转子位置估计的高频注入补偿方法
CN112829602A (zh) * 2021-01-04 2021-05-25 宝能(西安)汽车研究院有限公司 车辆扭矩控制方法、装置以及车辆
CN112698095A (zh) * 2021-01-22 2021-04-23 中南大学 一种基波和任意次谐波的高精度检测方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
CHEN DONGDONG 等: "On-Line Compensation of Resolver Periodic Error for PMSM Drives", IEEE TRANSACTIONS ON INDUSTRY APPLICATIONS, vol. 55, no. 6, pages 5990 - 6000, XP011755146, DOI: 10.1109/TIA.2019.2934089 *
THUONG NGOC-CONG TRAN 等: "Improving the Accuracy of an Absolute Magnetic Encoder by Using Harmonic Rejection and a Dual-Phase-Locked Loop", IEEE TRANS. INDUSTRIAL ELECTRONICS, vol. 66, no. 7, pages 5476 - 5486, XP011712577, DOI: 10.1109/TIE.2018.2869366 *
文武松;张颖超;王璐;詹天文;龙江涛;: "解耦双同步坐标系下单相锁相环技术", 电力系统自动化, no. 20 *
李伟伟 等: "考虑信号幅值和正交误差的旋转变压器解码算法", 电机与控制应用, vol. 42, no. 8, pages 17 - 23 *
李岩;于爽;夏加宽;姜昊驰;: "一种新型磁阻式正余弦旋转变压器解码方法", 电气工程学报, no. 09 *
罗皓文等: "一种基于同步旋转坐标系的指定次谐波补偿控制", 湖北工业大学学报, vol. 28, no. 1, pages 51 - 55 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117516596A (zh) * 2024-01-05 2024-02-06 华中科技大学 一种磁电编码器的高精度在线补偿方法
CN117516596B (zh) * 2024-01-05 2024-04-26 华中科技大学 一种磁电编码器的高精度在线补偿方法

Also Published As

Publication number Publication date
CN114061632B (zh) 2024-03-19

Similar Documents

Publication Publication Date Title
Zheng et al. A capacitive rotary encoder based on quadrature modulation and demodulation
JP4294558B2 (ja) 角度検出信号処理装置
CN113670340B (zh) 基于x/y信号相位差辨识的半球谐振陀螺驻波方位角测量方法及系统
US20070242043A1 (en) Contactless electron joystick of universal joint structure using single hall sensor
CN106610586B (zh) 一种基于并行forc和相位滞后-超前补偿的磁悬浮转子谐波电流抑制方法
CN113595453B (zh) 一种高响应旋转变压器解码方法及系统
Lara et al. Position error compensation in quadrature analog magnetic encoders through an iterative optimization algorithm
Zhang et al. Composite state observer for resolver-to-digital conversion
CN112859942B (zh) 一种陀螺框架伺服系统高精度位置辨识方法
CN114061632A (zh) 一种用于补偿指定次谐波的高精度磁编码器解码方法
Zhang et al. Absolute position acquisition for linear synchronous motor with passive mover
Qin et al. Angle Tracking observer with improved accuracy for resolver-to-digital conversion
Zhang et al. Automatic calibration of resolver signals via state observers
CN111245304B (zh) 补偿方法、补偿装置、电机和存储介质
CN117110643A (zh) 速率积分半球谐振陀螺角速率输出方法
CN116929322A (zh) 一种全角模式半球谐振陀螺精确控制方法及系统
CN116412837A (zh) 半球谐振陀螺力平衡模式校准电信号的生成方法
Rahman et al. Rotor Speed and Position Estimation Analysis of Interior PMSM Machines in Low and Medium-High Speed Regions Adopting an Improved Flux Observer for Electric Vehicle Applications
Wang et al. Study on high precision magnetic encoder based on PMSM sensorless control
Huang et al. An Improved Full-Speed Domain Sensorless Control Scheme for Permanent Magnet Synchronous Motor Based on Hybrid Position Observer and Disturbance Rejection Optimization
Yu et al. An Optimized Time Sequence for Sensorless Control of IPMSM Drives via High-Frequency Square-Wave Signal Injection Scheme
Karabeyli et al. Enhancing the accuracy for the open-loop resolver to digital converters
Ali LMI-based state feedback control structure for resolving grid connectivity issues in DFIG-based WT systems
JP4708070B2 (ja) 角度検出信号処理装置
Li et al. Model analysis and resonance suppression of wide-bandwidth inertial reference system

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