CN107607037A - A kind of magnetism encoder scaling method based on transregional tabulating method of arc tangent - Google Patents

A kind of magnetism encoder scaling method based on transregional tabulating method of arc tangent Download PDF

Info

Publication number
CN107607037A
CN107607037A CN201710838855.5A CN201710838855A CN107607037A CN 107607037 A CN107607037 A CN 107607037A CN 201710838855 A CN201710838855 A CN 201710838855A CN 107607037 A CN107607037 A CN 107607037A
Authority
CN
China
Prior art keywords
angle
measured
value
section
encoder
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
CN201710838855.5A
Other languages
Chinese (zh)
Other versions
CN107607037B (en
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 Frequency Control Technology Co.,Ltd.
Original Assignee
Harbin University of Science and Technology
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 Harbin University of Science and Technology filed Critical Harbin University of Science and Technology
Priority to CN201710838855.5A priority Critical patent/CN107607037B/en
Publication of CN107607037A publication Critical patent/CN107607037A/en
Application granted granted Critical
Publication of CN107607037B publication Critical patent/CN107607037B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Transmission And Conversion Of Sensor Element Output (AREA)

Abstract

The invention discloses a kind of magnetism encoder scaling method based on transregional tabulating method of arc tangent,It is related to magnetism encoder manufacturing field,The collection of signal is carried out by magnetism encoder first,By stages arc tangent inquiring arithmetic is recycled to calculate the measurement angle of magnetism encoder,By way of 360 ° of whole circumference is divided into 8 45 ° of subintervals,Avoid larger be also easy to produce of difficulty in computation when arc-tangent value has positive and negative infinite and calculate the problem of wrong,Finally benchmark is used as using 16 precision photoelectric encoder measured values,Demarcation amendment is carried out to the measurement angle that cutting method anyway is tried to achieve based on transregional average tabulating method,Error is eliminated by way of across intervals linear interpolation,This method can overcome that tradition has positive and negative infinite value in blanking method anyway and magnet steel position has the problem of machine error causes measurement accuracy to decline,Angle calibration is carried out to magnetism encoder.

Description

A kind of magnetism encoder scaling method based on transregional tabulating method of arc tangent
Technical field:
The present invention relates to a kind of magnetism encoder scaling method based on transregional tabulating method of arc tangent, belongs to magnetoelectricity coding Device manufacturing technology field.
Background technology:
Encoder is used for measurement angle, is the core parts for realizing motor control, is widely used in mechanical engineering, machine The high-tech sectors such as people, aviation, precision optical instrument, vital effect is played in modern industry.Wherein magnetoelectricity encodes The strong interference immunity of device, to the adaptable of the adverse circumstances such as dust and misty rain, code disk can be repeatedly recorded and wipe Remove, and the cost manufactured is relatively low relative to photoelectric encoder, is adapted to small lot production, industrial environment encodes to magnetoelectricity at present The demand of device is growing.
The precision of magnetism encoder is less than photoelectric encoder, and often it is demarcated using photoelectric encoder, passes through compensation The deviation of synchronism output is to improve the measurement accuracy of magnetism encoder between the two.Arc tangent algorithm is that magnetism encoder is main Angle computation method, this method sampled using Hall element to magnetic field signal, then carries out follow-up signal calculating processing, And obtain current angle value.
But traditional arc tangent algorithm is that difficulty in computation is larger when positive infinite sum bears infinite and easily draws in arc-tangent value Rise and calculate mistake, and mechanical erection deviation inevitably be present in the installation of Hall element, and magnet steel magnetization process are present partially Difference, the magnetic field signal as caused by magnet steel are not ideal sinusoidal waveforms, and the magnetic sensitive characteristic of each Hall element is different, can cause There is deviation in the magnetoelectricity signal of collection, it is digital quantity that magnetoelectricity signal, which is needed by analog-to-digital conversion, in addition, can be introduced in this process Noise causes angle calculation error, and above-mentioned error can cause the angle-measurement accuracy of magnetism encoder to decline.
The content of the invention:
In view of the above-mentioned problems, the technical problem to be solved in the present invention is to provide a kind of based on transregional tabulating method of arc tangent Magnetism encoder scaling method, magnetism encoder angle is calculated by by stages arc tangent inquiring arithmetic, and based on transregional Average tabulating method carries out angle calibration amendment, to improve the stability and measurement accuracy of magnetism encoder calculating.
Above-mentioned purpose is mainly realized by following scheme:
A kind of magnetism encoder scaling method based on transregional tabulating method of arc tangent of the present invention, it is characterised in that:Institute The specific implementation process for stating method is:
Step 1:The collection of magnetoelectricity signal
Magnetism encoder is pasted onto electricity using generating device of the single pair pole magnet steel of surface-mount type as signal, single pair pole magnet steel On the end face of machine rotating shaft, during axis of rotation, single pair pole magnet steel is with axis of rotation and produces sinusoidal magnetic field, four Hall element (A+、 A-、B+, B-) interval 90 ° be evenly distributed in Hall Plate, Hall Plate be arranged on electric motor end cap on, the uniform axis of Hall element Conllinear with shaft axis, Hall element is arranged towards magnet steel, magnetoelectricity analog signal is acquired using Hall component, 4 tunnels Hall signal passes through difference, and analog-to-digital conversion calculates d-q axle magnetoelectricity signals;
In view of the non-linear of magnetic field, Hall element A+The magnetoelectricity signal of collectionContaining DC bias and higher hamonic wave, Formula (1) can be used to represent:
Wherein θ is measurement angle,Σ is DC bias, △ c1For higher hamonic wave, on A- Hall elements The magnetoelectricity signal of sensingFormula (2) can be used to represent:
Similarly, B+, B- collection signal be respectively
By A+、A-、B+, the magnetoelectricity signal diameter that collects of B- Halls subtract each other to obtain d-q axle magnetoelectricity signals, such as formula (3) institute Show:
Step 2:Magnetism encoder angle is calculated using by stages arc tangent inquiring arithmetic
Magnetism encoder angle is calculated such as shown in (4) using arc tangent formula:
Whole circumference is divided into 8 sections, the angular range of each minizone is between 0 ° to 45 °, and arc-tangent value is 0 To between 1, can avoid in calculating process due to arc-tangent value occur it is just infinite and negative infinite and caused by result of calculation mistake Situation, by comparing positive and negative and its absolute value the size of d-q axle hall signals, judge residing angular interval and with accordingly Arc tangent formula calculates relative angle θ of the angle to be measured in residing section0', work as θ0' not in 0 ° to 45 ° of interval range when Need by adding and subtracting periodic transformation in 0 ° to 45 ° interval range, as shown in formula (5), i.e.,
θ00' ± nT, wherein T=45 °, n is integer (5)
When arc tangent formula used is arctan (dh/qh) when, relative angle adds the angle lower limit θ in residing sectionlI.e. It can obtain the actual angle θ of angle to be measuredT, as shown in formula (6);
θT0l (6)
When arc tangent formula used is arctan (qh/dh) when, the angle upper limit θ in residing sectionhSubtract relative angle i.e. It can obtain the actual angle θ of angle to be measuredT, as shown in formula (7);
θTh0 (7)
Such as when d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 0~45 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl=0 °, The then actual angle θ of angle to be measuredT0+0°;
When d axis signals is just, q axis signals are unsatisfactory for for just | dh|<|qh| when, angle to be measured is between 45~90 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 90 °, then the actual angle θ of angle to be measuredT=90 ° of-θ0
When d axis signals is just, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured be in 90~135 ° it Between, relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 135 °, then the actual angle θ of angle to be measuredT=135 ° of-θ0
When d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 135~180 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl= 135 °, then the actual angle θ of angle to be measuredT0+135°;
When d axis signals are negative, q axis signals are negative, and are met | dh|<|qh| when, angle to be measured is between 180~225 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl= 180 °, then the actual angle θ of angle to be measuredT0+180°;
When d axis signals are negative, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured be in 225~270 ° it Between, relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 270 °, then the actual angle θ of angle to be measuredT=270 ° of-θ0
When d axis signals is just, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured be in 270~315 ° it Between, relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 315 °, then the actual angle θ of angle to be measuredT=315 ° of-θ0
When d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 315~360 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl= 180 °, then the actual angle θ of angle to be measuredT0+315°;
Step 3:Angle calibration amendment based on transregional average tabulating method
There is angular deviation in the angle obtained by arc tangent algorithm, the angle value being calculated is carried out into angle calibration, Using the photoelectric encoder measured value of 16 precision (0-65535 represents 0-360 °) as benchmark, the magnetoelectricity that corrected Calculation comes out is compiled The measured value of code device;
Magnetism encoder, motor, photoelectric encoder are coaxially mounted to demarcate on turntable, motor rotation, by magnetoelectricity The angle value synchronism output that encoder and photoelectric encoder measure, obtains the angular deviation of magnetism encoder and photoelectric encoder Value, is handled the deviation of output by way of across intervals linear interpolation, such as to collection point MjCorresponding deviation is entered During row processing, first to MjThe adjacent four collection point M in left sidej-1, Mj-2, Mj-3, Mj-4Corresponding deviation asks for deviation average θp1, then To MjThe adjacent four collection point M in right sidej+1, Mj+2, Mj+3, Mj+4Corresponding deviation asks for error mean θp2, utilize θp1, θp2Line Property interpolation goes out the angle error value θ of target pointp0, as shown in formula (8), compiled all interpolation results as magnetism encoder and photoelectricity Angle compensation form is made in offset between code device, takes the offset in angle compensation form to be used as angle compensation amount by looking into The output result of magnetism encoder is modified, to improve the measurement accuracy of magnetism encoder.
The magnetism encoder angle obtained using arctangent computation corresponds to the integer type data of 16, due to high 10 data Be present fluctuation in stable and low 6 data, therefore only included in angle compensation form with the high ten digit value k (0-1023) of angle value For the deviation corresponding to collection point, low 6 bit value b (0-63) is obtained by local mesh subdivision, i.e., using the high 10 of angle value Numerical value tables look-up to obtain error correction values C (k) to form, and using low 6 bit value b (0-63) to the mistake that is checked under same high 10 Poor correction value is averaging processing, and corresponding low 6 bit value is distributed to after obtained average result is added up successively, is now passed through The final correction value C of magnetism encoder angular error after tabling look-up and segmentingF(k) formula (9) can be written as:
The then measured value θ of final magnetism encoder angleFBy being calculated using by stages arc tangent inquiring arithmetic The actual angle θ of angle to be measuredTWith obtaining the final correction value C of angular error based on transregional average tabulating methodF(k) plus and, such as Shown in formula (10):
θFT+CF(k) (10)
Beneficial effects of the present invention are:
1st, by stages arc tangent inquiring arithmetic is proposed, is overcome in traditional arc tangent algorithm, due to arc-tangent value Difficulty in computation is larger when bearing infinite for positive infinite sum and easily causes the situation for calculating mistake, improves the stabilization of arctangent computation Property.
2nd, transregional average tabulating method is proposed, angle compensation value is optimized for a smooth curve, avoids angle The introducing of system noise in compensation process, magnetism encoder and photoelectric encoder noise are eliminated to being compensated in angle compensation forms The harmful effect of value, the angle jitter of magnetism encoder is avoided, improve the measurement accuracy of magnetism encoder.
3rd, angle compensation is carried out by the way of angle compensation form using looking into take, reduce the calculating time of control system, carry The high response speed of encoder.
Brief description of the drawings:
For ease of explanation, the present invention is described in detail by following specific implementations and accompanying drawing.
Fig. 1 is the demarcation schematic flow sheet that the inventive method obtains;
Fig. 2 is that the computation interval of by stages arc tangent inquiring arithmetic in the present invention divides schematic diagram;
Fig. 3 is the magnetoelectricity angular encoder and light for carrying out angle calibration amendment in the present invention without transregional average tabulating method The angle difference figure of photoelectric coder;
Fig. 4 is the magnetoelectricity angular encoder and light for carrying out angle calibration amendment in the present invention without transregional average tabulating method The once differentiation value figure of the angle difference of photoelectric coder;
Fig. 5 is the magnetoelectricity angular encoder and light for carrying out angle calibration amendment in the present invention by transregional average tabulating method The angle difference figure of photoelectric coder;
Fig. 6 is the magnetoelectricity angular encoder and light for carrying out angle calibration amendment in the present invention by transregional average tabulating method The once differentiation value figure of the angle difference of photoelectric coder.
Embodiment:
It is specific below by what is shown in accompanying drawing to make the object, technical solutions and advantages of the present invention of greater clarity Embodiment describes the present invention.However, it should be understood that these descriptions are merely illustrative, and it is not intended to limit the model of the present invention Enclose.In addition, in the following description, the description to known features and technology is eliminated, to avoid unnecessarily obscuring the present invention's Concept.
As shown in Figure 1, Figure 2, shown in Fig. 3, Fig. 4, Fig. 5, Fig. 6, present embodiment uses following technical scheme:One kind is based on The magnetism encoder scaling method of transregional tabulating method of arc tangent, it is characterised in that:The specific implementation process of methods described is:
Step 1:The collection of magnetoelectricity signal
Magnetism encoder is pasted onto electricity using generating device of the single pair pole magnet steel of surface-mount type as signal, single pair pole magnet steel On the end face of machine rotating shaft, during axis of rotation, single pair pole magnet steel is with axis of rotation and produces sinusoidal magnetic field, four Hall element (A+、 A-、B+, B-) interval 90 ° be evenly distributed in Hall Plate, Hall Plate be arranged on electric motor end cap on, the uniform axis of Hall element Conllinear with shaft axis, Hall element is arranged towards magnet steel, magnetoelectricity analog signal is acquired using Hall component, 4 tunnels Hall signal passes through difference, and analog-to-digital conversion calculates d-q axle magnetoelectricity signals;
In view of the non-linear of magnetic field, Hall element A+The magnetoelectricity signal of collectionContaining DC bias and higher hamonic wave, Formula (1) can be used to represent:
Wherein θ is measurement angle,Σ is DC bias, △ c1For higher hamonic wave, on A- Hall elements The magnetoelectricity signal of sensingFormula (2) can be used to represent:
Similarly, B+, B- collection signal be respectively
By A+、A-、B+, the magnetoelectricity signal diameter that collects of B- Halls subtract each other to obtain d-q axle magnetoelectricity signals, such as formula (3) institute Show:
Step 2:Magnetism encoder angle is calculated using by stages arc tangent inquiring arithmetic
Magnetism encoder angle is calculated such as shown in (4) using arc tangent formula:
Whole circumference is divided into 8 sections, the angular range of each minizone is between 0 ° to 45 °, and arc-tangent value is 0 To between 1, can avoid in calculating process due to arc-tangent value occur it is just infinite and negative infinite and caused by result of calculation mistake Situation, by comparing positive and negative and its absolute value the size of d-q axle hall signals, judge residing angular interval and with accordingly Arc tangent formula calculates relative angle θ of the angle to be measured in residing section0', work as θ0' not in 0 ° to 45 ° of interval range when Need by adding and subtracting periodic transformation in 0 ° to 45 ° interval range, as shown in formula (5), i.e.,
θ00' ± nT, wherein T=45 °, n is integer (5)
When arc tangent formula used is arctan (dh/qh) when, relative angle adds the angle lower limit θ in residing sectionlI.e. It can obtain the actual angle θ of angle to be measuredT, as shown in formula (6);
θT0l (6)
When arc tangent formula used is arctan (qh/dh) when, the angle upper limit θ in residing sectionhSubtract relative angle i.e. It can obtain the actual angle θ of angle to be measuredT, as shown in formula (7);
θTh0 (7)
Such as when d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 0~45 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl=0 °, The then actual angle θ of angle to be measuredT0+0°;
When d axis signals is just, q axis signals are unsatisfactory for for just | dh|<|qh| when, angle to be measured is between 45~90 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 90 °, then the actual angle θ of angle to be measuredT=90 ° of-θ0
When d axis signals is just, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured be in 90~135 ° it Between, relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 135 °, then the actual angle θ of angle to be measuredT=135 ° of-θ0
When d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 135~180 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl= 135 °, then the actual angle θ of angle to be measuredT0+135°;
When d axis signals are negative, q axis signals are negative, and are met | dh|<|qh| when, angle to be measured is between 180~225 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl= 180 °, then the actual angle θ of angle to be measuredT0+180°;
When d axis signals are negative, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured be in 225~270 ° it Between, relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 270 °, then the actual angle θ of angle to be measuredT=270 ° of-θ0
When d axis signals is just, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured be in 270~315 ° it Between, relative angle θ of the angle to be measured in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh= 315 °, then the actual angle θ of angle to be measuredT=315 ° of-θ0
When d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 315~360 °, Relative angle θ of the angle to be measured in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl= 180 °, then the actual angle θ of angle to be measuredT0+315°;
Step 3:Angle calibration amendment based on transregional average tabulating method
There is angular deviation in the angle obtained by arc tangent algorithm, the angle value being calculated is carried out into angle calibration, Using the photoelectric encoder measured value of 16 precision (0-65535 represents 0-360 °) as benchmark, the magnetoelectricity that corrected Calculation comes out is compiled The measured value of code device;
Magnetism encoder, motor, photoelectric encoder are coaxially mounted to demarcate on turntable, motor rotation, by magnetoelectricity The angle value synchronism output that encoder and photoelectric encoder measure, obtains the angular deviation of magnetism encoder and photoelectric encoder Value, is handled the deviation of output by way of across intervals linear interpolation, such as to collection point MjCorresponding deviation is entered During row processing, first to MjThe adjacent four collection point M in left sidej-1, Mj-2, Mj-3, Mj-4Corresponding deviation asks for deviation average θp1, then To MjThe adjacent four collection point M in right sidej+1, Mj+2, Mj+3, Mj+4Corresponding deviation asks for error mean θp2, utilize θp1, θp2Line Property interpolation goes out the angle error value θ of target pointp0, as shown in formula (8), compiled all interpolation results as magnetism encoder and photoelectricity Angle compensation form is made in offset between code device, takes the offset in angle compensation form to be used as angle compensation amount by looking into The output result of magnetism encoder is modified, to improve the measurement accuracy of magnetism encoder.
The magnetism encoder angle obtained using arctangent computation corresponds to the integer type data of 16, due to high 10 data Be present fluctuation in stable and low 6 data, therefore only included in angle compensation form with the high ten digit value k (0-1023) of angle value For the deviation corresponding to collection point, low 6 bit value b (0-63) is obtained by local mesh subdivision, i.e., using the high 10 of angle value Numerical value tables look-up to obtain error correction values C (k) to form, and using low 6 bit value b (0-63) to the mistake that is checked under same high 10 Poor correction value is averaging processing, and corresponding low 6 bit value is distributed to after obtained average result is added up successively, is now passed through The final correction value C of magnetism encoder angular error after tabling look-up and segmentingF(k) formula (9) can be written as:
The then measured value θ of final magnetism encoder angleFBy being calculated using by stages arc tangent inquiring arithmetic The actual angle θ of angle to be measuredTWith obtaining the final correction value C of angular error based on transregional average tabulating methodF(k) plus and, such as Shown in formula (10):
θFT+CF(k) (10)
The general principle and principal character and advantages of the present invention of the present invention has been shown and described above.The technology of the industry Personnel are it should be appreciated that the present invention is not limited to the above embodiments, and the simply explanation described in above-described embodiment and specification is originally The principle of invention, without departing from the spirit and scope of the present invention, various changes and modifications of the present invention are possible, these changes Change and improvement all fall within the protetion scope of the claimed invention.The claimed scope of the invention by appended claims and its Equivalent thereof.

Claims (1)

  1. A kind of 1. magnetism encoder scaling method based on transregional tabulating method of arc tangent, it is characterised in that:The tool of methods described Body implementation process is:
    Step 1:The collection of magnetoelectricity signal
    Using generating device of the single pair pole magnet steel of surface-mount type as signal, single pair pole magnet steel is pasted onto motor and turned magnetism encoder On the end face of axle, during axis of rotation, single pair pole magnet steel is with axis of rotation and produces sinusoidal magnetic field, four Hall element (A+、A-、 B+, B-) interval 90 ° be evenly distributed in Hall Plate, Hall Plate be arranged on electric motor end cap on, the uniform axis of Hall element with turn Axle axis collinear, Hall element is arranged towards magnet steel, magnetoelectricity analog signal is acquired using Hall component, 4 road Halls Signal passes through difference, and analog-to-digital conversion calculates d-q axle magnetoelectricity signals;
    In view of the non-linear of magnetic field, Hall element A+The magnetoelectricity signal of collectionContaining DC bias and higher hamonic wave, formula can be used (1) represent:
    Wherein θ is measurement angle,Σ is DC bias, △ c1For higher hamonic wave, sense on A- Hall elements Magnetoelectricity signalFormula (2) can be used to represent:
    Similarly, B+, B- collection signal be respectively
    By A+、A-、B+, the magnetoelectricity signal diameter that collects of B- Halls subtract each other to obtain d-q axle magnetoelectricity signals, as shown in formula (3):
    Step 2:Magnetism encoder angle is calculated using by stages arc tangent inquiring arithmetic
    Magnetism encoder angle is calculated such as shown in (4) using arc tangent formula:
    Whole circumference is divided into 8 sections, the angular range of each minizone is between 0 ° to 45 °, and arc-tangent value is 0 to 1 Between, can avoid in calculating process due to arc-tangent value occur it is just infinite and negative infinite and caused by result of calculation mistake feelings Condition, by comparing positive and negative and its absolute value the size of d-q axle hall signals, judge residing angular interval and with anti-accordingly Tangent formula calculates relative angle θ of the angle to be measured in residing section0', work as θ0' not in 0 ° to 45 ° of interval range when need Will be by adding and subtracting periodic transformation in 0 ° to 45 ° interval range, as shown in formula (5), i.e.,
    θ00' ± nT, wherein T=45 °, n is integer (5)
    When arc tangent formula used is arctan (dh/qh) when, relative angle adds the angle lower limit θ in residing sectionlIt can obtain To the actual angle θ of angle to be measuredT, as shown in formula (6);
    θT0l (6)
    When arc tangent formula used is arctan (qh/dh) when, the angle upper limit θ in residing sectionhSubtracting relative angle can obtain To the actual angle θ of angle to be measuredT, as shown in formula (7);
    θTh0 (7)
    Such as when d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 0~45 °, to be measured Relative angle θ of the angle in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl=0 °, then treat The actual angle θ of measuring angleT0+0°;
    When d axis signals is just, q axis signals are unsatisfactory for for just | dh|<|qh| when, angle to be measured is between 45~90 °, to be measured Relative angle θ of the angle in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh=90 °, then The actual angle θ of angle to be measuredT=90 ° of-θ0
    When d axis signals is just, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured is between 90~135 °, is treated Relative angle θ of the measuring angle in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh=135 °, The then actual angle θ of angle to be measuredT=135 ° of-θ0
    When d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 135~180 °, to be measured Relative angle θ of the angle in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl=135 °, then The actual angle θ of angle to be measuredT0+135°;
    When d axis signals are negative, q axis signals are negative, and are met | dh|<|qh| when, angle to be measured is between 180~225 °, to be measured Relative angle θ of the angle in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl=180 °, then The actual angle θ of angle to be measuredT0+180°;
    When d axis signals are negative, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured is between 225~270 °, is treated Relative angle θ of the measuring angle in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh=270 °, The then actual angle θ of angle to be measuredT=270 ° of-θ0
    When d axis signals is just, q axis signals are negative, and are unsatisfactory for | dh|<|qh| when, angle to be measured is between 270~315 °, is treated Relative angle θ of the measuring angle in the section0' can be by arctan (qh/dh) try to achieve, the angle upper limit θ in residing sectionh=315 °, The then actual angle θ of angle to be measuredT=315 ° of-θ0
    When d axis signals is just, q axis signals be negative, and satisfaction | dh|<|qh| when, angle to be measured is between 315~360 °, to be measured Relative angle θ of the angle in the section0' can be by arctan (dh/qh) try to achieve, the angle lower limit θ in residing sectionl=180 °, then The actual angle θ of angle to be measuredT0+315°;
    Step 3:Angle calibration amendment based on transregional average tabulating method
    There is angular deviation in the angle obtained by arc tangent algorithm, the angle value being calculated is carried out into angle calibration, by 16 The magnetism encoder that the photoelectric encoder measured value of position precision (0-65535 represents 0-360 °) comes out as benchmark, corrected Calculation Measured value;
    Magnetism encoder, motor, photoelectric encoder are coaxially mounted to demarcate on turntable, motor rotation, magnetoelectricity encoded The angle value synchronism output that device and photoelectric encoder measure, obtains the angu-lar deviation of magnetism encoder and photoelectric encoder, leads to The mode crossed across intervals linear interpolation is handled the deviation of output, and such as count collection MjAt corresponding deviation During reason, first to MjThe adjacent four collection point M in left sidej-1, Mj-2, Mj-3, Mj-4Corresponding deviation asks for deviation average θp1, then to Mj The adjacent four collection point M in right sidej+1, Mj+2, Mj+3, Mj+4Corresponding deviation asks for error mean θp2, utilize θp1, θp2It is linear to insert Mend out the angle error value θ of target pointp0, as shown in formula (8), using all interpolation results as magnetism encoder and photoelectric encoder Between offset angle compensation form is made, take the offset in angle compensation form to be used as angle compensation amount to magnetic by looking into The output result of photoelectric coder is modified, to improve the measurement accuracy of magnetism encoder.
    The magnetism encoder angle obtained using arctangent computation corresponds to the integer type data of 16, due to high 10 data stabilizations It is to adopt only to be included and low 6 data have fluctuation, therefore in angle compensation form with the high ten digit value k (0-1023) of angle value The corresponding deviation of collection point, low 6 bit value b (0-63) is obtained by local mesh subdivision, i.e., using high 10 bit value of angle value Table look-up to obtain error correction values C (k) to form, and the error checked under same high 10 is repaiied using low 6 bit value b (0-63) On the occasion of being averaging processing, corresponding low 6 bit value is distributed to after obtained average result is added up successively, now by tabling look-up And the final correction value C of magnetism encoder angular error after subdivisionF(k) formula (9) can be written as:
    The then measured value θ of final magnetism encoder angleFFor by the angle measurement being calculated using by stages arc tangent inquiring arithmetic The actual angle θ of degreeTWith obtaining the final correction value C of angular error based on transregional average tabulating methodF(k) plus and, such as formula (10) It is shown:
    θFT+CF(k) (10)。
CN201710838855.5A 2017-09-18 2017-09-18 Magnetoelectric encoder calibration method based on arc tangent trans-interval tabulation method Active CN107607037B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710838855.5A CN107607037B (en) 2017-09-18 2017-09-18 Magnetoelectric encoder calibration method based on arc tangent trans-interval tabulation method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710838855.5A CN107607037B (en) 2017-09-18 2017-09-18 Magnetoelectric encoder calibration method based on arc tangent trans-interval tabulation method

Publications (2)

Publication Number Publication Date
CN107607037A true CN107607037A (en) 2018-01-19
CN107607037B CN107607037B (en) 2020-05-19

Family

ID=61060726

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710838855.5A Active CN107607037B (en) 2017-09-18 2017-09-18 Magnetoelectric encoder calibration method based on arc tangent trans-interval tabulation method

Country Status (1)

Country Link
CN (1) CN107607037B (en)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108180933A (en) * 2018-02-02 2018-06-19 哈尔滨理工大学 It is a kind of based on magnetism encoder automatic correcting method of the permanent magnet synchronous motor without sensor speed control
CN108871181A (en) * 2018-05-07 2018-11-23 哈尔滨理工大学 A kind of multipair pole magnetoelectric encoder dynamic multiwindow interval prediction angle sorting method
CN108896074A (en) * 2018-06-04 2018-11-27 哈尔滨理工大学 Eight Halls are distributed magnetism encoder high-resolution angle value calculation method and device
CN109059968A (en) * 2018-07-02 2018-12-21 哈尔滨理工大学 Finite angle angle value divided method and device based on angle value growth
CN109099938A (en) * 2018-07-04 2018-12-28 哈尔滨理工大学 The angular interval scan angle angle value jump point suppressing method and device tabled look-up based on number of poles
CN110793491A (en) * 2019-11-25 2020-02-14 中国航空工业集团公司北京航空精密机械研究所 High-precision absolute type external synchronization angle capturing device and method
CN110940371A (en) * 2019-12-13 2020-03-31 浙江禾川科技股份有限公司 Calibration method, device and equipment of rotary magnetoelectric encoder
CN111006696A (en) * 2019-12-16 2020-04-14 横店集团英洛华电气有限公司 Magnetic encoder and angle calculation method thereof
CN111044982A (en) * 2019-12-23 2020-04-21 珠海纳睿达科技有限公司 Radar azimuth positioning method
CN111043951A (en) * 2018-10-15 2020-04-21 大银微***股份有限公司 Torsion sensing mechanism of rotating shaft
CN111146986A (en) * 2019-12-30 2020-05-12 深圳市越疆科技有限公司 Position locating method and device of magnetic encoder, electronic equipment and computer readable storage medium
CN111504360A (en) * 2020-05-19 2020-08-07 哈尔滨理工大学 Time coordinate-based magnetoelectric encoder angle precision division method and device
CN111578970A (en) * 2020-05-21 2020-08-25 哈尔滨理工大学 Angle precision division method and device for multi-pair-electrode magnetoelectric encoder based on nucleic acid sequence comparison idea
CN111693073A (en) * 2020-06-28 2020-09-22 哈尔滨理工大学 Bidirectional redundant magnetoelectric encoder and redundancy detection method thereof
CN111721329A (en) * 2020-07-07 2020-09-29 哈尔滨理工大学 three-Hall magnetoelectric encoder and arc-tangent-free calculation angle calculation method
CN111750903A (en) * 2020-07-07 2020-10-09 哈尔滨理工大学 Winding integrated magnetoelectric encoder and independent calibration method thereof
CN111811546A (en) * 2020-07-26 2020-10-23 哈尔滨理工大学 Electromagnetic encoder for generating magnetism by electrified coil and resolving method thereof
CN111854807A (en) * 2020-05-15 2020-10-30 哈尔滨理工大学 Domain angle magnetoelectric encoder and signal logic extension fine-dividing method thereof
CN112067023A (en) * 2020-07-26 2020-12-11 哈尔滨理工大学 two-Hall magnetoelectric encoder and angle value temperature drift compensation method thereof
CN114440940A (en) * 2022-03-17 2022-05-06 哈尔滨理工大学 double-Hall magnetoelectric encoder angle estimation method and device based on special-shaped gear
CN114866669A (en) * 2022-04-26 2022-08-05 杭州海康威视数字技术股份有限公司 Video camera and voice coil motor driving method for video camera

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101876559A (en) * 2009-04-30 2010-11-03 浙江关西电机有限公司 Position detecting device and signal processing device and method thereof
CN102686979A (en) * 2009-12-28 2012-09-19 法国欧陆汽车公司 Method for determining the position of a magnetic element using linear hall-effect sensors and associated device
US20140015457A1 (en) * 2012-07-12 2014-01-16 Kia Motors Corporation System and method for calibrating offset of motor resolver
CN103837169A (en) * 2014-02-28 2014-06-04 哈尔滨工业大学 Self-correcting device and method for magneto-electricity encoder and magneto-electricity encoder

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101876559A (en) * 2009-04-30 2010-11-03 浙江关西电机有限公司 Position detecting device and signal processing device and method thereof
CN102686979A (en) * 2009-12-28 2012-09-19 法国欧陆汽车公司 Method for determining the position of a magnetic element using linear hall-effect sensors and associated device
US20140015457A1 (en) * 2012-07-12 2014-01-16 Kia Motors Corporation System and method for calibrating offset of motor resolver
CN103837169A (en) * 2014-02-28 2014-06-04 哈尔滨工业大学 Self-correcting device and method for magneto-electricity encoder and magneto-electricity encoder

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郝双晖等: "绝对式多极磁电轴角编码器的设计", 《光学精密工程》 *

Cited By (33)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108180933A (en) * 2018-02-02 2018-06-19 哈尔滨理工大学 It is a kind of based on magnetism encoder automatic correcting method of the permanent magnet synchronous motor without sensor speed control
CN108180933B (en) * 2018-02-02 2023-07-25 哈尔滨理工大学 Magneto-electric encoder self-correction method based on sensorless speed control of permanent magnet synchronous motor
CN108871181A (en) * 2018-05-07 2018-11-23 哈尔滨理工大学 A kind of multipair pole magnetoelectric encoder dynamic multiwindow interval prediction angle sorting method
CN108896074A (en) * 2018-06-04 2018-11-27 哈尔滨理工大学 Eight Halls are distributed magnetism encoder high-resolution angle value calculation method and device
CN108896074B (en) * 2018-06-04 2020-11-03 哈尔滨理工大学 Method and device for resolving high-resolution angle value of eight-Hall distributed magnetoelectric encoder
CN109059968A (en) * 2018-07-02 2018-12-21 哈尔滨理工大学 Finite angle angle value divided method and device based on angle value growth
CN109059968B (en) * 2018-07-02 2020-09-25 哈尔滨理工大学 Limited angle value subdivision method and device based on angle value growth
CN109099938A (en) * 2018-07-04 2018-12-28 哈尔滨理工大学 The angular interval scan angle angle value jump point suppressing method and device tabled look-up based on number of poles
CN109099938B (en) * 2018-07-04 2020-11-24 哈尔滨理工大学 Method and device for suppressing jumping points of angle value during angle interval scanning based on pole number table look-up
CN111043951A (en) * 2018-10-15 2020-04-21 大银微***股份有限公司 Torsion sensing mechanism of rotating shaft
CN110793491A (en) * 2019-11-25 2020-02-14 中国航空工业集团公司北京航空精密机械研究所 High-precision absolute type external synchronization angle capturing device and method
CN110940371A (en) * 2019-12-13 2020-03-31 浙江禾川科技股份有限公司 Calibration method, device and equipment of rotary magnetoelectric encoder
CN111006696A (en) * 2019-12-16 2020-04-14 横店集团英洛华电气有限公司 Magnetic encoder and angle calculation method thereof
CN111044982B (en) * 2019-12-23 2021-09-28 广东纳睿雷达科技股份有限公司 Radar azimuth positioning method
CN111044982A (en) * 2019-12-23 2020-04-21 珠海纳睿达科技有限公司 Radar azimuth positioning method
CN111146986A (en) * 2019-12-30 2020-05-12 深圳市越疆科技有限公司 Position locating method and device of magnetic encoder, electronic equipment and computer readable storage medium
CN111146986B (en) * 2019-12-30 2022-08-12 深圳市越疆科技有限公司 Position locating method and device of magnetic encoder, electronic equipment and computer readable storage medium
CN111854807A (en) * 2020-05-15 2020-10-30 哈尔滨理工大学 Domain angle magnetoelectric encoder and signal logic extension fine-dividing method thereof
CN111504360B (en) * 2020-05-19 2022-07-05 哈尔滨理工大学 Time coordinate-based magnetoelectric encoder angle precision division method and device
CN111504360A (en) * 2020-05-19 2020-08-07 哈尔滨理工大学 Time coordinate-based magnetoelectric encoder angle precision division method and device
CN111578970A (en) * 2020-05-21 2020-08-25 哈尔滨理工大学 Angle precision division method and device for multi-pair-electrode magnetoelectric encoder based on nucleic acid sequence comparison idea
CN111693073B (en) * 2020-06-28 2022-01-25 哈尔滨理工大学 Bidirectional redundant magnetoelectric encoder and redundancy detection method thereof
CN111693073A (en) * 2020-06-28 2020-09-22 哈尔滨理工大学 Bidirectional redundant magnetoelectric encoder and redundancy detection method thereof
CN111721329B (en) * 2020-07-07 2021-11-23 哈尔滨理工大学 three-Hall magnetoelectric encoder and arc-tangent-free calculation angle calculation method
CN111750903A (en) * 2020-07-07 2020-10-09 哈尔滨理工大学 Winding integrated magnetoelectric encoder and independent calibration method thereof
CN111721329A (en) * 2020-07-07 2020-09-29 哈尔滨理工大学 three-Hall magnetoelectric encoder and arc-tangent-free calculation angle calculation method
CN112067023A (en) * 2020-07-26 2020-12-11 哈尔滨理工大学 two-Hall magnetoelectric encoder and angle value temperature drift compensation method thereof
CN111811546A (en) * 2020-07-26 2020-10-23 哈尔滨理工大学 Electromagnetic encoder for generating magnetism by electrified coil and resolving method thereof
CN111811546B (en) * 2020-07-26 2022-05-24 哈尔滨理工大学 Electromagnetic encoder for generating magnetism by electrified coil and resolving method thereof
CN114440940A (en) * 2022-03-17 2022-05-06 哈尔滨理工大学 double-Hall magnetoelectric encoder angle estimation method and device based on special-shaped gear
CN114440940B (en) * 2022-03-17 2022-11-01 哈尔滨理工大学 double-Hall magnetoelectric encoder angle estimation method and device based on special-shaped gear
CN114866669A (en) * 2022-04-26 2022-08-05 杭州海康威视数字技术股份有限公司 Video camera and voice coil motor driving method for video camera
CN114866669B (en) * 2022-04-26 2024-06-07 杭州海康威视数字技术股份有限公司 Video camera and voice coil motor driving method for video camera

Also Published As

Publication number Publication date
CN107607037B (en) 2020-05-19

Similar Documents

Publication Publication Date Title
CN107607037A (en) A kind of magnetism encoder scaling method based on transregional tabulating method of arc tangent
CN102111158B (en) Device for subdividing sine signal and cosine signal of position sensor and coding data, and implementation method thereof
CN104568041B (en) A kind of direct-reading gage for eliminating adjacent runner magnetic interference
CN114858192B (en) double-Hall magnetoelectric encoder based on double-wheel structure and angle resolving method thereof
CN107588793B (en) Magnetic angle sensor calibration method based on discrete sine-cosine transform
CN108180933B (en) Magneto-electric encoder self-correction method based on sensorless speed control of permanent magnet synchronous motor
CN106767956B (en) Magnetic induction absolute value encoder of high-speed high-precision machine tool spindle and measuring gear thereof
JP6331177B1 (en) Angle sensor system
US20110025312A1 (en) Magnetic encoder and method of detecting absolute rotational position
CN203881354U (en) Combined type encoder
CN103983290A (en) Composite type absolute value encoder
CN109781150B (en) Control method of magnetic encoder, magnetic encoder and knitting machine
CN104613991B (en) Encoder grating disc and photoelectric encoder
CN201748928U (en) Tunnel magnetoresistance effect magnetic encoder
CN201858990U (en) Device for subdividing sine-cosine signals and encoding data of position sensor
CN103983291B (en) Coding method of photomagnetic coder
CN210071933U (en) Signal sampling circuit
CN112304210B (en) Hall angle sensor angle calibration method
CN101982736A (en) Code disk of electro-optic rotary encoder
CN103604447A (en) Method for realizing high-resolution incremental bus-based optical-electricity encoder
CN104655156B (en) A kind of error compensating method of the absolute type code-disc of matrix form coded system
CN205333002U (en) Measurement turns to angle sensor detection device
CN205986831U (en) Decoder and because electrical resolver of this decoder for electrical resolver
CN114001768B (en) Self-calibration device of magnetoelectric encoder
JP2010078340A (en) Method for compensating error of encoder

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
TR01 Transfer of patent right

Effective date of registration: 20211014

Address after: 201108 No. 535, Shennan Road, Xinzhuang Industrial Zone, Minhang District, Shanghai

Patentee after: Shanghai Frequency Control Technology Co.,Ltd.

Address before: 150080 No. 52, Xuefu Road, Nangang District, Heilongjiang, Harbin

Patentee before: HARBIN University OF SCIENCE AND TECHNOLOGY