CN111193448A - Surface-mounted permanent magnet synchronous motor load torque observation method based on extended Kalman filter - Google Patents

Surface-mounted permanent magnet synchronous motor load torque observation method based on extended Kalman filter Download PDF

Info

Publication number
CN111193448A
CN111193448A CN202010062836.XA CN202010062836A CN111193448A CN 111193448 A CN111193448 A CN 111193448A CN 202010062836 A CN202010062836 A CN 202010062836A CN 111193448 A CN111193448 A CN 111193448A
Authority
CN
China
Prior art keywords
permanent magnet
load torque
magnet synchronous
synchronous motor
current
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
CN202010062836.XA
Other languages
Chinese (zh)
Other versions
CN111193448B (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.)
Jiangsu Simand Electric Co Ltd
Original Assignee
Jiangsu Simand Electric 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 Jiangsu Simand Electric Co Ltd filed Critical Jiangsu Simand Electric Co Ltd
Priority to CN202010062836.XA priority Critical patent/CN111193448B/en
Publication of CN111193448A publication Critical patent/CN111193448A/en
Application granted granted Critical
Publication of CN111193448B publication Critical patent/CN111193448B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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/13Observer control, e.g. using Luenberger observers or Kalman filters
    • 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/20Estimation of torque
    • 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
    • H02P2207/055Surface mounted magnet motors

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)
  • Control Of Electric Motors In General (AREA)

Abstract

The invention provides a surface-mounted permanent magnet synchronous motor load torque observation method based on an extended Kalman filter, which is characterized in that a mathematical model of a permanent magnet synchronous motor for load torque observation is obtained according to a mathematical model under a two-phase rotating coordinate system of the permanent magnet synchronous motor on the basis of vector control, and is substituted into an extended Kalman filter algorithm to realize the observation of the load torque of the permanent magnet synchronous motor. The surface-mounted permanent magnet synchronous motor load torque observation method based on the extended Kalman filter realizes real-time monitoring of the surface-mounted permanent magnet synchronous motor load torque by using the extended Kalman filter algorithm, can realize sensorless monitoring of the permanent magnet synchronous motor load torque, and has the advantages of simple structure, low cost and stronger anti-interference capability.

Description

Surface-mounted permanent magnet synchronous motor load torque observation method based on extended Kalman filter
Technical Field
The invention relates to the technical field of permanent magnet synchronous motors, in particular to a surface-mounted permanent magnet synchronous motor load torque observation method based on an extended Kalman filter.
Background
Permanent Magnet Synchronous Motors (PMSM) have the advantages of high power, high energy density and simple structure, and are increasingly widely applied in numerous industrial fields. In controlling the operation of the motor, the motor is controlled to provide accurate torque control over a full speed range and a full load range. In addition, by comparing the given command torque with the actual torque of the drive motor system, the operation state of the drive motor system can be indirectly reflected, and whether the fault exists or not can be deduced. Therefore, the real-time monitoring of the torque is not only a requirement for ensuring the control performance, but also an important safety index.
The traditional torque measuring method measures the torque through a torque measuring instrument, and has the defects of high measuring cost, large influence of instrument precision, limitation of installation conditions in some occasions, inconvenience in maintenance, incapability of realizing real-time online monitoring of the load torque and the like.
Disclosure of Invention
In order to solve the problems in the prior art, the invention provides the surface-mounted permanent magnet synchronous motor load torque observation method based on the extended Kalman filter, overcomes the defects in the prior art, and has the advantages of predictability, self-adaptive capacity, anti-interference performance, simple structure, easiness in software implementation, low cost and the like.
In order to achieve the purpose, the invention adopts the following technical scheme:
the surface-mounted permanent magnet synchronous motor load torque observation method based on the extended Kalman filter comprises the following steps:
step 1: obtaining actual rotation degree of a surface-mounted permanent magnet synchronous motor and three-phase current i of the motorabcThree-phase current iabcObtaining the current i under a static two-phase coordinate system through Clarke transformationαAnd iβObtaining the current i under a rotating two-phase coordinate system through Park transformationdAnd iq
Step 2: the difference value of the given rotating speed of the motor and the actual rotating speed of the motor is adjusted by a rotating speed Adjuster (ASR) to determine the given current i of the q axisq *D-axis current is given current id *
And step 3: setting q axis to current iq *And a feedback current iqAnd d-axis given current id *And a feedback current idThe difference value of the d-axis voltage and the q-axis voltage u are obtained after being adjusted by A Current Regulator (ACR)d、uq
And 4, step 4: d and q axis real time voltage ud、uqα and β axis real-time voltage u is obtained through inverse PARK conversionα、uβ
Step 5, real-time voltage u of α and β axesα、uβSVPWM modulation is carried out to obtain a pulse width modulation waveform, the pulse width modulation waveform is sent to an inverter to control a permanent magnet synchronous motor, and real-time three-phase current i is obtainedabcAnd three phase voltage uabc
Step 6: will measure the real-time current id、iqReal time voltage ud、uqAnd real-time electrical angular velocity weAnd sending the load torque to a load torque observer based on an Extended Kalman Filter (EKF) to realize the real-time observation of the load torque of the permanent magnet synchronous motor.
Furthermore, in the surface-mounted permanent magnet synchronous motor, the magnetic path structure of the rotor is symmetrical, and the magnetic permeability of the permanent magnet material is close to that of the air gap, so that the L-shaped magnetic path structure is L-shapedd=LqUnder a two-phase rotating coordinate system d-q, the voltage equation of the permanent magnet synchronous motor is as follows:
Figure BDA0002375047300000031
in the formula ud、uq、id、iqVoltages and currents of a d axis and a q axis of the two-phase rotating coordinate system are respectively; rs、L、ψfStator resistance, quadrature-direct axis inductance and rotor permanent magnet flux linkage; w is aeIs the electrical angular velocity of the motor;
transforming equation (1) into a current equation yields:
Figure BDA0002375047300000032
and the electromagnetic torque is:
Figure BDA0002375047300000033
wherein p isnThe number of pole pairs of the permanent magnet synchronous motor is;
the mechanical equation of motion of the motor is:
Figure BDA0002375047300000034
wherein wmIs the mechanical angular velocity of the motor, and J is the moment of inertia; b is the damping coefficient, TLIs a load torque, and an electrical angle we=pnwm
So the mechanical motion equation of the transformed motor is:
Figure BDA0002375047300000035
in practical situations, the motor load torque fluctuation time is far longer than the motor control system dynamic process time, so that the load torque is calculated by regarding as a steady-state value, and derivation of the load torque can obtain the following equation:
Figure BDA0002375047300000036
the following can be written by the formulas (2), (5) and (6):
Figure BDA0002375047300000041
equation (7) can be written as the following equation of state:
Figure BDA0002375047300000042
wherein x ═ idiqweTL]TIs a state variable, an input variable u ═ uduq]TAnd input variable y ═ idiq]T
The following formulas (7) and (8) can be obtained:
Figure BDA0002375047300000043
Figure BDA0002375047300000044
Figure BDA0002375047300000045
in order to construct the digital system of the extended Kalman filter state observer, discretization processing can obtain
Figure BDA0002375047300000051
Finishing to obtain:
Figure BDA0002375047300000052
where T is the sampling period, the discretized state equation (14) is a deterministic equation, but in practical systems, model parameters are uncertain and variable, measurement noise inevitably exists in stator voltages and currents, and these uncertainties are incorporated into the system noise vector V and measurement noise W, so equation (14) can be changed to:
Figure BDA0002375047300000053
where V (k) is system noise and W (k) is measurement noise, both of which are zero-mean white noise.
Further, the specific process of step 6 is as follows:
the discretized mathematical model of the EKF algorithm is as follows:
Figure BDA0002375047300000054
wherein x (k +1) represents a state estimation value at the k +1 moment, x (k) represents a state estimation value at the k moment, V (k) is a system noise vector, and W (k) is measurement noise;
the state estimation of the extended kalman filter observer is mainly divided into two states: prediction and correction
(1) And (3) state prediction:
Figure BDA0002375047300000055
in the formula (I), the compound is shown in the specification,
Figure BDA0002375047300000061
and
Figure BDA0002375047300000062
indicates the (k +1) -th predicted value,
Figure BDA0002375047300000063
represents the k-th estimated value;
the above equation is mainly based on the input u (k) and last state estimation
Figure BDA0002375047300000064
Predicting a vector at the (k +1) time;
(2) calculating a covariance matrix:
Figure BDA0002375047300000065
in the formula (I), the compound is shown in the specification,
Figure BDA0002375047300000066
is a covariance matrix, mainly to find the gain matrix K (K +1),
Figure BDA0002375047300000067
an error covariance matrix at the moment k is obtained, and Q is a covariance matrix of system noise V; and the Jacobian matrix F (k) is:
Figure BDA0002375047300000068
the calculation result is as follows:
Figure BDA0002375047300000069
(3) calculating a gain matrix:
Figure BDA00023750473000000610
wherein, K (K +1) is a gain matrix which is mainly used for completing the correction of state vector estimation, and R is a covariance matrix of a measurement noise vector W;
(4) and (3) state vector estimation:
Figure BDA00023750473000000613
in the formula (I), the compound is shown in the specification,
Figure BDA00023750473000000611
is the state vector estimation value of the (k +1) th time, y (k +1) is the measurement state vector,
Figure BDA00023750473000000612
for the predicted output state vector, the above equation completes the state vector estimation of the (k +1) th;
(5) calculating an estimation error covariance matrix:
Figure BDA0002375047300000071
in the formula (I), the compound is shown in the specification,
Figure BDA0002375047300000072
the covariance matrix reflects the magnitude of the state estimation error at this time, and is called in the next state estimation, so that iterative operation is performed, and the load torque observed value at each moment is obtained.
Compared with the prior art, the invention has the beneficial technical effects that: according to the observation method for the load torque of the surface-mounted permanent magnet synchronous motor based on the extended Kalman filter, the observer in the form of an algorithm replaces the traditional parameter measurement method in the physical form, the load torque of the surface-mounted permanent magnet synchronous motor is monitored in real time by using the algorithm of the extended Kalman filter, and the observation method has the advantages of low cost, simple structure, easiness in software implementation, strong self-adaption capability and anti-jamming capability and the like.
Drawings
FIG. 1 is a schematic diagram of a vector control system according to an embodiment of the present invention;
FIG. 2 is a schematic diagram of a load torque observation process based on EKF algorithm according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the embodiments of the present invention will be clearly and completely described below with reference to the drawings of the embodiments of the present invention. It is to be understood that the embodiments described are only a few embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the described embodiments of the invention, are within the scope of the invention.
The technical scheme of the invention is further explained in detail by combining the drawings in the specification.
As shown in fig. 1, the method for observing the load torque of the surface-mounted permanent magnet synchronous motor based on the extended kalman filter comprises the following steps:
step 1: obtaining actual rotation degree of a surface-mounted permanent magnet synchronous motor and three-phase current i of the motorabcThree-phase current iabcObtaining the current i under a static two-phase coordinate system through Clarke transformationαAnd iβObtaining the current i under a rotating two-phase coordinate system through Park transformationdAnd iq
Step 2: the difference value of the given rotating speed of the motor and the actual rotating speed of the motor is adjusted by a rotating speed Adjuster (ASR) to determine the given current i of the q axisq *D-axis current is given current id *
And step 3: setting q axis to current iq *And contrary toCurrent feed iqAnd d-axis given current id *And a feedback current idThe difference value of the d-axis voltage and the q-axis voltage u are obtained after being adjusted by A Current Regulator (ACR)d、uq
And 4, step 4: d and q axis real time voltage ud、uqα and β axis real-time voltage u is obtained through inverse PARK conversionα、uβ
Step 5, real-time voltage u of α and β axesα、uβSVPWM modulation is carried out to obtain a pulse width modulation waveform, the pulse width modulation waveform is sent to an inverter to control a permanent magnet synchronous motor, and real-time three-phase current i is obtainedabcAnd three phase voltage uabc
Step 6: will measure the real-time current id、iqReal time voltage ud、uqAnd real-time electrical angular velocity weAnd sending the load torque to a load torque observer based on an Extended Kalman Filter (EKF) to realize the real-time observation of the load torque of the permanent magnet synchronous motor.
In the surface-mounted permanent magnet synchronous motor, the rotor magnetic circuit structure is symmetrical, and the magnetic conductivity of the permanent magnet material is close to that of the air gap, so that the L-shaped magnetic circuit structure is L-shapedd=LqUnder a two-phase rotating coordinate system d-q, the voltage equation of the permanent magnet synchronous motor is as follows:
Figure BDA0002375047300000081
in the formula ud、uq、id、iqVoltages and currents of a d axis and a q axis of the two-phase rotating coordinate system are respectively; rs、L、ψfStator resistance, quadrature-direct axis inductance and rotor permanent magnet flux linkage; w is aeIs the electrical angular velocity of the motor;
transforming equation (1) into a current equation yields:
Figure BDA0002375047300000091
and the electromagnetic torque is:
Figure BDA0002375047300000092
wherein p isnThe number of pole pairs of the permanent magnet synchronous motor is;
the mechanical equation of motion of the motor is:
Figure BDA0002375047300000093
wherein wmIs the mechanical angular velocity of the motor, and J is the moment of inertia; b is the damping coefficient, TLIs a load torque, and an electrical angle we=pnwm
So the mechanical motion equation of the transformed motor is:
Figure BDA0002375047300000094
in practical situations, the motor load torque fluctuation time is far longer than the motor control system dynamic process time, so that the load torque is calculated by regarding as a steady-state value, and derivation of the load torque can obtain the following equation:
Figure BDA0002375047300000095
the following can be written by the formulas (2), (5) and (6):
Figure BDA0002375047300000101
equation (7) can be written as the following equation of state:
Figure BDA0002375047300000102
wherein x ═ idiqweTL]TIs a state variable, an input variable u ═ uduq]TAnd input variable y ═ idiq]T
The following formulas (7) and (8) can be obtained:
Figure BDA0002375047300000103
Figure BDA0002375047300000104
Figure BDA0002375047300000105
in order to construct the digital system of the extended Kalman filter state observer, discretization processing can obtain
Figure BDA0002375047300000111
Finishing to obtain:
Figure BDA0002375047300000112
where T is the sampling period, the discretized state equation (14) is a deterministic equation, but in practical systems, model parameters are uncertain and variable, measurement noise inevitably exists in stator voltages and currents, and these uncertainties are incorporated into the system noise vector V and measurement noise W, so equation (14) can be changed to:
Figure BDA0002375047300000113
where V (k) is system noise and W (k) is measurement noise, both of which are zero-mean white noise.
The specific process of step 6 is as follows:
the discretized mathematical model of the EKF algorithm is as follows:
Figure BDA0002375047300000114
wherein x (k +1) represents a state estimation value at the k +1 moment, x (k) represents a state estimation value at the k moment, V (k) is a system noise vector, and W (k) is measurement noise;
the state estimation of the extended kalman filter observer is mainly divided into two states: prediction and correction
(1) And (3) state prediction:
Figure BDA0002375047300000115
in the formula (I), the compound is shown in the specification,
Figure BDA0002375047300000121
and
Figure BDA0002375047300000122
indicates the (k +1) -th predicted value,
Figure BDA0002375047300000123
represents the k-th estimated value;
the above equation is mainly based on the input u (k) and last state estimation
Figure BDA00023750473000001213
Predicting a vector at the (k +1) time;
(2) calculating a covariance matrix:
Figure BDA0002375047300000124
in the formula (I), the compound is shown in the specification,
Figure BDA0002375047300000125
is a covariance matrix, mainly to find the gain matrix K (K +1),
Figure BDA0002375047300000126
an error covariance matrix at the moment k is obtained, and Q is a covariance matrix of system noise V; and the Jacobian matrix F (k) is:
Figure BDA0002375047300000127
the calculation result is as follows:
Figure BDA0002375047300000128
(3) calculating a gain matrix:
Figure BDA0002375047300000129
wherein, K (K +1) is a gain matrix which is mainly used for completing the correction of state vector estimation, and R is a covariance matrix of a measurement noise vector W;
(4) and (3) state vector estimation:
Figure BDA00023750473000001210
in the formula (I), the compound is shown in the specification,
Figure BDA00023750473000001211
is the state vector estimation value of the (k +1) th time, y (k +1) is the measurement state vector,
Figure BDA00023750473000001212
for the predicted output state vector, the above equation completes the state vector estimation of the (k +1) th;
(5) calculating an estimation error covariance matrix:
Figure BDA0002375047300000131
in the formula (I), the compound is shown in the specification,
Figure BDA0002375047300000132
the covariance matrix reflects the magnitude of the state estimation error at this time, and is called in the next state estimation, so that iterative operation is performed, and the load torque observed value at each moment is obtained.
According to the observation method for the load torque of the surface-mounted permanent magnet synchronous motor based on the extended Kalman filter, the observer in the form of an algorithm replaces the traditional parameter measurement method in the physical form, the load torque of the surface-mounted permanent magnet synchronous motor is monitored in real time by using the algorithm of the extended Kalman filter, and the observation method has the advantages of low cost, simple structure, easiness in software implementation, strong self-adaption capability and anti-jamming capability and the like.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (3)

1. The surface-mounted permanent magnet synchronous motor load torque observation method based on the extended Kalman filter is characterized by comprising the following steps:
step 1: obtaining actual rotation degree of a surface-mounted permanent magnet synchronous motor and three-phase current i of the motorabcThree-phase current iabcObtaining the current i under a static two-phase coordinate system through Clarke transformationαAnd iβObtaining the current i under a rotating two-phase coordinate system through Park transformationdAnd iq
Step 2: the difference value of the given rotating speed of the motor and the actual rotating speed of the motor is adjusted by a rotating speed Adjuster (ASR) to determine the given current i of the q axisq *D-axis current is given current id *
And step 3: setting q axis to current iq *And a feedback current iqAnd d-axis given current id *And a feedback current idThe difference value of the d-axis voltage and the q-axis voltage u are obtained after being adjusted by A Current Regulator (ACR)d、uq
And 4, step 4: d and q axis real time voltage ud、uqα and β axis real-time voltage u is obtained through inverse PARK conversionα、uβ
Step 5, real-time voltage u of α and β axesα、uβSVPWM modulation is carried out to obtain a pulse width modulation waveform, and the pulse width modulation waveform is sent to the inverseThe transformer controls the permanent magnet synchronous motor to obtain real-time three-phase current iabcAnd three phase voltage uabc
Step 6: will measure the real-time current id、iqReal time voltage ud、uqAnd real-time electrical angular velocity weAnd sending the load torque to a load torque observer based on an Extended Kalman Filter (EKF) to realize the real-time observation of the load torque of the permanent magnet synchronous motor.
2. The extended Kalman filter-based surface-mounted permanent magnet synchronous motor load torque observation method according to claim 1, characterized in that:
in the surface-mounted permanent magnet synchronous motor, the magnetic path structure of the rotor is symmetrical, and the magnetic conductivity of the permanent magnet material is close to that of the air gap, so that the L-shaped magnetic flux path structure is L-shapedd=LqUnder a two-phase rotating coordinate system d-q, the voltage equation of the permanent magnet synchronous motor is as follows:
Figure FDA0002375047290000021
in the formula ud、uq、id、iqVoltages and currents of a d axis and a q axis of the two-phase rotating coordinate system are respectively; rs、L、ψfStator resistance, quadrature-direct axis inductance and rotor permanent magnet flux linkage; w is aeIs the electrical angular velocity of the motor;
transforming equation (1) into a current equation yields:
Figure FDA0002375047290000022
and the electromagnetic torque is:
Figure FDA0002375047290000023
wherein p isnThe number of pole pairs of the permanent magnet synchronous motor is;
the mechanical equation of motion of the motor is:
Figure FDA0002375047290000024
wherein wmIs the mechanical angular velocity of the motor, and J is the moment of inertia; b is the damping coefficient, TLIs a load torque, and an electrical angular velocity we=pnwm
So the mechanical motion equation of the transformed motor is:
Figure FDA0002375047290000025
in practical situations, the motor load torque fluctuation time is far longer than the motor control system dynamic process time, so that the load torque is calculated by regarding as a steady-state value, and derivation of the load torque can obtain the following equation:
Figure FDA0002375047290000026
the following can be written by the formulas (2), (5) and (6):
Figure FDA0002375047290000031
equation (7) can be written as the following equation of state:
Figure FDA0002375047290000032
wherein x ═ idiqweTL]TIs a state variable, an input variable u ═ uduq]TAnd input variable y ═ idiq]T
The following formulas (7) and (8) can be obtained:
Figure FDA0002375047290000033
Figure FDA0002375047290000034
Figure FDA0002375047290000035
in order to construct the digital system of the extended Kalman filter state observer, discretization processing can obtain
Figure FDA0002375047290000041
Finishing to obtain:
Figure FDA0002375047290000042
where T is the sampling period, the discretized state equation (14) is a deterministic equation, but in practical systems, model parameters are uncertain and variable, measurement noise inevitably exists in stator voltages and currents, and these uncertainties are incorporated into the system noise vector V and measurement noise W, so equation (14) can be changed to:
Figure FDA0002375047290000043
where V (k) is system noise and W (k) is measurement noise, both of which are zero-mean white noise.
3. The extended Kalman filter-based surface-mounted permanent magnet synchronous motor load torque observation method according to claim 1, characterized in that: the specific process of step 6 is as follows:
the discretized mathematical model of the EKF algorithm is as follows:
Figure FDA0002375047290000044
wherein x (k +1) represents a state estimation value at the k +1 moment, x (k) represents a state estimation value at the k moment, V (k) is a system noise vector, and W (k) is measurement noise;
the state estimation of the extended kalman filter observer is mainly divided into two states: prediction and correction
(1) And (3) state prediction:
Figure FDA0002375047290000051
in the formula (I), the compound is shown in the specification,
Figure FDA0002375047290000052
and
Figure FDA0002375047290000053
indicates the (k +1) -th predicted value,
Figure FDA0002375047290000054
represents the k-th estimated value;
the above equation is mainly based on the input u (k) and last state estimation
Figure FDA00023750472900000512
Predicting a vector at the (k +1) time;
(2) calculating a covariance matrix:
Figure FDA0002375047290000055
in the formula (I), the compound is shown in the specification,
Figure FDA0002375047290000056
is a covariance matrix, mainly to find the gain matrix K (K +1),
Figure FDA0002375047290000057
an error covariance matrix at the moment k is obtained, and Q is a covariance matrix of system noise V; and the Jacobian matrix F (k) is:
Figure FDA0002375047290000058
the calculation result is as follows:
Figure FDA0002375047290000059
(3) calculating a gain matrix:
Figure FDA00023750472900000510
wherein, K (K +1) is a gain matrix which is mainly used for completing the correction of state vector estimation, and R is a covariance matrix of a measurement noise vector W;
(4) and (3) state vector estimation:
Figure FDA00023750472900000511
in the formula (I), the compound is shown in the specification,
Figure FDA0002375047290000061
is the state vector estimation value of the (k +1) th time, y (k +1) is the measurement state vector,
Figure FDA0002375047290000062
for the predicted output state vector, the above equation completes the state vector estimation of the (k +1) th;
(5) calculating an estimation error covariance matrix:
Figure FDA0002375047290000063
in the formula (I), the compound is shown in the specification,
Figure FDA0002375047290000064
the covariance matrix reflects the magnitude of the state estimation error, and is called in the next state estimation, so that iterative operation is performed to obtain each timeMoment load torque observations.
CN202010062836.XA 2020-01-20 2020-01-20 Permanent magnet synchronous motor load torque observation method based on extended Kalman filter Active CN111193448B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010062836.XA CN111193448B (en) 2020-01-20 2020-01-20 Permanent magnet synchronous motor load torque observation method based on extended Kalman filter

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010062836.XA CN111193448B (en) 2020-01-20 2020-01-20 Permanent magnet synchronous motor load torque observation method based on extended Kalman filter

Publications (2)

Publication Number Publication Date
CN111193448A true CN111193448A (en) 2020-05-22
CN111193448B CN111193448B (en) 2022-11-04

Family

ID=70708339

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010062836.XA Active CN111193448B (en) 2020-01-20 2020-01-20 Permanent magnet synchronous motor load torque observation method based on extended Kalman filter

Country Status (1)

Country Link
CN (1) CN111193448B (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112865642A (en) * 2021-01-18 2021-05-28 南京航空航天大学 Harmonic suppression system and method for line-controlled steering permanent magnet synchronous motor
CN112910328A (en) * 2021-01-22 2021-06-04 绍兴敏动科技有限公司 Permanent magnet synchronous motor acceleration arrangement method based on torque observation compensation
CN113315439A (en) * 2021-06-02 2021-08-27 江苏科技大学 Improved control method of ship propulsion motor
CN113708684A (en) * 2021-08-31 2021-11-26 哈尔滨理工大学 Permanent magnet synchronous motor control method and device based on extended potential observer
CN116101364A (en) * 2022-12-23 2023-05-12 吉林大学 Control method of power-assisted steering motor

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109412488A (en) * 2018-10-29 2019-03-01 天水电气传动研究所有限责任公司 A kind of permanent magnet synchronous motor dynamic matrix control method
CN110138301A (en) * 2019-04-29 2019-08-16 江苏大学 A kind of permanent magnet synchronous motor position-sensorless control method based on spreading kalman filter
CN113131817A (en) * 2021-04-09 2021-07-16 陕西科技大学 Online parameter identification system and method for permanent magnet synchronous motor

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109412488A (en) * 2018-10-29 2019-03-01 天水电气传动研究所有限责任公司 A kind of permanent magnet synchronous motor dynamic matrix control method
CN110138301A (en) * 2019-04-29 2019-08-16 江苏大学 A kind of permanent magnet synchronous motor position-sensorless control method based on spreading kalman filter
CN113131817A (en) * 2021-04-09 2021-07-16 陕西科技大学 Online parameter identification system and method for permanent magnet synchronous motor

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
谷涛等: ""基于EKF的永磁同步电机无位置传感器控制"", 《电源技术》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112865642A (en) * 2021-01-18 2021-05-28 南京航空航天大学 Harmonic suppression system and method for line-controlled steering permanent magnet synchronous motor
CN112910328A (en) * 2021-01-22 2021-06-04 绍兴敏动科技有限公司 Permanent magnet synchronous motor acceleration arrangement method based on torque observation compensation
CN113315439A (en) * 2021-06-02 2021-08-27 江苏科技大学 Improved control method of ship propulsion motor
CN113708684A (en) * 2021-08-31 2021-11-26 哈尔滨理工大学 Permanent magnet synchronous motor control method and device based on extended potential observer
CN113708684B (en) * 2021-08-31 2022-09-30 哈尔滨理工大学 Permanent magnet synchronous motor control method and device based on extended potential observer
CN116101364A (en) * 2022-12-23 2023-05-12 吉林大学 Control method of power-assisted steering motor

Also Published As

Publication number Publication date
CN111193448B (en) 2022-11-04

Similar Documents

Publication Publication Date Title
CN111193448B (en) Permanent magnet synchronous motor load torque observation method based on extended Kalman filter
Tsuji et al. A sensorless vector control system for induction motors using q-axis flux with stator resistance identification
CN102420561B (en) Speed sensorless vector control method on basis of cascaded high voltage inverter
KR101087581B1 (en) Sensorless control method of permanent magnet synchronous motor
CN109194219B (en) Method and system for controlling permanent magnet synchronous motor based on model-free nonsingular terminal sliding mode
CN107070341A (en) Permagnetic synchronous motor method for suppressing torque ripple based on Robust Iterative Learning Control
CN111181458A (en) Surface-mounted permanent magnet synchronous motor rotor flux linkage observation method based on extended Kalman filter
CN103872951A (en) Permanent magnet synchronous motor torque control method based on sliding mode flux linkage observer
CN108092567A (en) A kind of Speed control of permanent magnet synchronous motor system and method
CN108183651A (en) A kind of control method of permanent magnet linear synchronous motor position-sensor-free
CN110995102A (en) Direct torque control method and system for permanent magnet synchronous motor
CN108306570A (en) Direct torque control method for permanent magnetic synchronous electric machine and system
Barambones Caramazana et al. Sliding mode position control for real-time control of induction motors
Chen et al. Adaptive second-order active-flux observer for sensorless control of pmsms with mras-based vsi nonlinearity compensation
Badini et al. MRAS-based speed and parameter estimation for a vector-controlled PMSM drive
Stănică et al. A brief review of sensorless motors position control
Wu et al. Sensorless model reference adaptive control of permanent magnet synchronous motor for industrial robots
CN114465543B (en) Permanent magnet synchronous motor sensorless control method
Mahawar et al. Parmanent magnet synchronous motor control based on model reference adaptive approach
CN111130409B (en) Ultra-low speed high-precision positioning control method of frameless torque motor
Zheng et al. Sensorless speed control for brushless DC motors system using sliding-mode controller and observers
Brock et al. Reducing Energy Losses for Fan Applications with V/f control of PMSMs
Saberi et al. A full range permanent magnet synchronous motor position and speed estimation using adaptive non-singular fast terminal sliding mode observer
Liao et al. Speed control of dual sensorless interior permanent-magnet synchronous motors based on the sliding mode observer principle
Zhang et al. Current Prediction Control Strategy of Permanent Magnet Synchronous Motor Based on Extended Kalman Filter

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