Method of identifying the parameters of a permanent magnet synchronous motor
Patent Information
- Authority / Receiving Office
- EP · EP
- Patent Type
- Applications
- Current Assignee / Owner
- Filing Date
- 2024-04-30
- Publication Date
- 2026-03-11
AI Technical Summary
Existing methods for identifying parameters of permanent magnet synchronous motors (PMSMs) face challenges such as rank deficiency, inverter nonlinearity, and actuation delays, leading to inaccurate parameter estimation and increased complexity and cost, especially in industrial settings where motors operate in parallel.
A method that identifies steady states and uses data from two operating conditions to minimize estimation errors, compensates for inverter nonlinearity by identifying distorted voltages, and adjusts for actuation delays, ensuring accurate estimation of stator winding resistances and rotor magnet flux without invasive signal injection or new drive control schemes.
This approach provides a non-invasive, cost-effective method for accurate parameter identification, reducing estimation errors and maintaining motor performance, applicable in industrial environments without requiring experimental testing or new control schemes.
Smart Images

Figure IB2024054188_07112024_PF_FP_ABST
Abstract
Description
[0001] “METHOD OF IDENTIFYING THE PARAMETERS OF A PERMANENT MAGNET
[0002] SYNCHRONO US MOTOR ”
[0003] Description
[0004] Technical Field
[0005] The present invention relates to electric motors and control systems therefor, and more particularly it concerns a method of identifying the parameters of a permanent magnet synchronous motor (hereinafter also referred to with the acronym PMSM, from the English denomination 'Permanent Magnet Synchronous Motor'), especially an isotropic motor with surface-mounted magnets.
[0006] Preferably, though not exclusively, the invention finds application in identifying parameters of motors used in capping devices for applying caps or closures to containers such as bottles and the like.
[0007] Background Art
[0008] Knowledge of the parameters of a permanent magnet synchronous motor, such as stator winding inductances and resistances and flux generated by the rotor magnets, is required at several stages of the motor life, e.g. design validation, control system tuning, monitoring of operating conditions and fault diagnostics. Considering that these parameters can vary considerably during the operation of the machine in which the motor is inserted, the ability to monitor their evolution over time represents a very interesting opportunity in the field of preventive and predictive maintenance and fault prevention.
[0009] Numerous methods for estimating the parameters of a PMSM are known in literature, and these can be classified into two main families: off-line methods and on-line methods.
[0010] Off-line methods are based on the a posteriori processing of data acquired through various motor tests, which are generally carried out with the motor disconnected from its machine.
[0011] In contrast, on-line methods are based on algorithms implemented in real time by the control unit of the electrical drive of the machine using the motor, and exploit the measurement data commonly available in such drives (phase currents, rotor position and rotational speed, supply voltage references fed to the inverter), which are collected and processed during ordinary operation of the motor, without disconnecting it from its machine.
[0012] While ensuring high accuracy, off-line methods have the problem of requiring the use of experienced personnel to disconnect the motor from the application in which it is used and the use of specific laboratory instrumentation, which makes them particularly costly in an industrial production environment where several motors operate in parallel on different machines. In addition, they require continuous downtime to disconnect and reconnect the motors to be tested, which is detrimental to the productivity of the plant in which the motors are used.
[0013] Known on-line methods, while exempt from the above-mentioned problems, are themselves affected by two distinct problems.
[0014] The first problem is that they can hardly be implemented from scratch on commercial drives not specifically designed for this purpose. In fact, the software of commercial drive control units is generally not modifiable to implement customized identification algorithms. Furthermore, the computational load required by on-line identification algorithms may be too much for the control unit.
[0015] The second problem is that the estimation accuracy of on-line methods is marred by the rank deficiency of the system being identified and the deficiency of measurement of the voltages actually applied to the motor.
[0016] Rank deficiency is due to the fact that the number of parameters to be identified is greater than the number of motor equations available at each step of sampling of the measured quantities (current, position, speed and voltage), which leads, among other things, to non-unambiguous parameter identification. Two different approaches have been proposed to overcome this problem.
[0017] The first approach is based on reducing the parameters to be identified, said reduction being achieved by associating part of these parameters with nameplate values or with values acquired from other estimates and measurements. However, this first approach has poor accuracy, as the variations that these parameters undergo as the operating conditions of the motor (current, speed and temperature) change are not taken into account.
[0018] The second approach, used also in the invention, is instead based on increasing the available equations. For this purpose, it was proposed to use additional data obtained by perturbing the operating point of the motor with the inj ection of current, voltage or position signals. An example is described in CN111030534A. Apart from the fact that signal injection in commercial drive systems is not always possible, the critical issue with this solution is its invasiveness. Indeed, signal injection can reduce the performance of the motor drive in terms of energy efficiency and control accuracy.
[0019] Regarding the problem of the lack of measurements of the voltages actually applied to the motor, on-line parametric identification algorithms generally use, as estimates of such lacking measurements, the reference supply voltages sent by the control unit to the inverter. However, because of the non-linearity of the inverter and implementation delays there is a discrepancy between actual voltages and reference voltages, and thus this procedure can give rise to non- negligible estimation errors. In principle, it has been shown that the non-linearity of the inverter can be compensated for by introducing, into the motor model equations, an additional parameter to be identified. The discrepancy caused by actuation delays increases as the ratio of the electrical speed of the motor to the frequency of the inverter pulse-width modulation increases. This is generally overlooked in the literature concerning on-line identification methods, but it can cause significant estimation errors in medium-low to high speed applications.
[0020] CN113141140A discloses a method for the on-line identification of the parameters of a PMSM with surface-mounted magnets which takes into account the problems of rank deficiency and actuation delay. The problem of rank deficiency is solved by using data acquired during a single operating condition of the motor. This known solution is inherently inaccurate because it does not include compensation for the non-linearity of the inverter. Furthermore, the use of data acquired during a single operating condition of the motor does not allow for the best possible minimisation of estimation errors.
[0021] CN108183648A discloses a method for the on-line identification of the parameters of a PMSM with compensation of the non-linearity of the inverter and actuation delays. The known method identifies the parameters by performing experimental DC tests on the motor and adopts an ad hoc electrical drive control scheme to compensate for the non-linearity of the inverter. The method is thus invasive and is rendered complex and costly by the need to design a new drive control scheme, which, moreover, must generally be adapted on a case-by-case basis.
[0022] BRESCIAELIAET AL: “Parameter Identification of PMSMs Considering VSI Nonlinearity with Coupled Adaline NNs”, 2023 IEEE 3RD International Conference on Power Electronics and Computer Applications (ICPECA), IEEE, 29 January 2023, pages 259-265 XP034319234 discloses a method of identifying the parameters of a permanent magnet synchronous motor. Other methods are disclosed in BRESCIAELIAET AL: “Automated Parameter Identification of SPMSMs Based on Two Steady States Using Cloud Computing Resources” 2021 International Conference on Electrical, Computer and Energy Technologies (ICECET), IEEE, 9 December 2021, pages 1-6, XP034082182, and LI FUMIN ET AL. : “Optimal FOPID Error Voltage Control Dead-Time Compensation Based on FOPI current control for PMSM Servo System”, 2023 International Conference on Fractional Differentiation and its Applications (ICFDA), IEEE, 14 March 2023, pages 1-6, XP034362840.
[0023] Summary of Invention
[0024] The object of the present invention is to provide a method of identifying the parameters of a permanent magnet synchronous motor that solves the problems of rank deficiency, inverter nonlinearity and actuation delays and overcomes the drawbacks of prior art.
[0025] This object is achieved with a method wherein the steady states of the motor are identified and the operating conditions of said steady states, namely current, speed and temperature, of said steady states are extracted, and wherein, for identifying the at least the stator winding resistances and the flux generated by the rotor magnets, the data of two operating conditions of different steady states are used, which operating conditions are selected so as to minimise estimation errors for the parameter being identified.
[0026] Preferably, the steady state identification is performed by using the R-statistic algorithm, and the operating condition extraction includes slicing the data of each steady state into data subsets relating to time intervals in which the motor temperature is such that temperature-dependent parameters of the motor are substantially constant, the parameter identification using, for each operating condition, the data of each subset.
[0027] Advantageously, the compensation of the non-linearity of the inverter includes identifying a voltage, hereinafter referred to as distorted voltage, representative of the effects of such nonlinearity, such an identification being carried out by minimising, for all samples of each operating condition of each steady state, the difference between a voltage actually applied to the motor and a reference voltage.
[0028] According to a preferred feature of the invention, the minimization of the estimation errors for the stator winding resistances and the flux generated by the rotor magnets includes minimising an error comprising a first error component due to the difference in the resistance and flux values in the two operating conditions, and a second error component due to an error in the compensation of the actuation delays by which the reference voltage is affected because of the inverter non-linearity.
[0029] Advantageously, in this case, the first error component is determined from a rough estimation of the stator winding resistances and the flux generated by the rotor magnets obtained by using predetermined values of motor temperature and rotor speed coefficients.
[0030] According to another feature of the invention, the minimisation of the estimation errors for the stator winding resistances and the flux generated by the rotor magnets includes imposing the constraints that the majorant of the estimation error is lower than a desired fraction of the respective motor parameter, and that the value of a parameter allowing evaluation of the linear independence of the equations describing a mathematical method of the motor in each of the two operating conditions selected is not close to a value denoting the lack of independence of said equations.
[0031] A further object of the invention is a capping device actuated by a permanent magnet synchronous motor and controlled by a control system, wherein said control system comprises at least one processing unit programmed for implementing the method according to the invention.
[0032] A still further object of the invention is a computer program product loaded into the memory of at least one processing unit, in particular a processing unit belonging to a control system of a device for capping containers such as bottles and the like actuated by a permanent magnet synchronous motor, the product including portions of software code for implementing the aforesaid method when the product is executed on the ate least one processing unit.
[0033] Brief Description of Drawings
[0034] These and other features and advantages of the present invention will become evident from the following description of a preferred embodiment given by way of non-limiting example with reference to the annexed figure, which shows a general flowchart of the method according to the invention.
[0035] Description of Embodiments
[0036] A method of identifying the parameters of an isotropic PMSM with surface-mounted magnets is described below. The method is based on a static mathematical model of the motor, in particular as far as the coordinates d, q are concerned. This model requires using only measurement data obtained during the steady state of the motor, i.e. in an operating condition in which both the absorbed currents and the speed are constant over time.
[0037] That being stated, the first step 101 of the method according to invention consists of collecting data from the control unit of the electrical drive of the motor. As mentioned, the data are those commonly available in PMSM control units, i.e. phase currents, rotor angular position and rotational speed, motor temperature and supply voltage references supplied to the inverter.
[0038] The collected data are used in the subsequent step 102 for identifying the steady states of the motor. To discriminate the data obtained during the steady state from the data obtained during the transient state, an algorithm known as R-statistic, applied to two different measured variables of the motor, i.e. q-axis current and mechanical speed ct>r of the rotor. For each sample k, the algorithm assesses independently whether these two variables are in the transient state or in the steady state, by also using the former samples k-N, where N is selected by the user. The formula of the R-statistic for identifying the steady states of mechanical speed coralone is shown below: where cor,n indicates the measured speed to which some noise has been added to improve detection performance. This algorithm is in fact designed to analyse variables affected by noise. The motor speed is considered as being in the steady state at step k if RUr(k) < Rcrt, where Rcrt is an index close to 1 (typical values are 1.2 -1.4) arbitrarily selected by the user with a trial-and- error method.
[0039] To facilitate understanding of the subsequent steps of the method, the equations describing the discrete time model of the motor within the reference system d, q are given below: where: k indicates the kthsample; u*d and u*q are the reference voltages of the d-axis and q-axis, respectively; iq is the current of the q-axis; co is the electrical speed of the rotor (mechanical speed multiplied by the number of polar pairs); Lq, R and \| / mare the stator inductance of the q- axis, the stator resistance and the concatenated magnetic flux of the rotor, respectively; Vdead, also known as distorted voltage, is an unknown voltage that takes into account the non-linearity phenomena introduced by the inverter; and the coefficients Dj and Dq, which according to literature are obtained by studying pulse width modulation inverters, are functions of the phase currents and rotor position and can therefore be calculated from the available measurements.
[0040] According to the invention, the voltage Vdead is also identified to compensate for the nonlinearity phenomena.
[0041] It should be noted that the equations do not contain the d-axis current id, because in isotropic motors this current is usually kept equal to 0 with an appropriate control action: indeed, in these motors, as is known to the person skilled in the art, the current id does not contribute to developing torque and therefore to doing useful work, but it only increases the power loss, thereby reducing the efficiency of the motor. In addition, there is only one inductance component, because in isotropic motors, under conditions of magnetic non-saturation, the inductances Ld and Lq of the two axes are the same.
[0042] After identifying the steady states, at step 103 there is effected a compensation of the actuation delay, i.e. of the fact that the control voltages of the inverter calculated by the control unit at instant k-Ts, where Tsis the sampling time of the control unit, are actually applied to the motor at instant (k+l,5)Ts. This delay is due to two different contributions. The first contribution is related to the digital implementation of the control algorithm and is due to the fact that the voltages calculated at a certain sampling instant are stored in memory to be then sent to the inverter at the instant of the subsequent sampling: this contribution is therefore a pure delay equal to Ts. The second contribution is related to the fact that the inverter receiving the voltages at step k+1 can in turn apply them to the motor after 0.5 Tson average. Because of this actuation delay, in the mathematical model of the motor it is not correct to associate the voltages detected at step k with the other variables measured at the same step k, because, as mentioned, the voltages at step k are applied only after 1.5 Tsand therefore are not actually equal to those that will be applied after 1.5 TS. In fact, the voltages at step k are calculated by taking as reference the electric position 9(k) of the rotor at said step k, but, when they are applied after 1.5 TS, the rotor will have performed a rotation of about:
[0043] A9 = 1.5 Tsro(k) = 1 ,5(9(k)-9(k-l)), where co(k) is the electrical speed measure at step k. The rotation A9 is known because both 9(k) and 0(k- 1) are known.
[0044] The compensation of the actuation delay then consists of replacing the reference voltages u*d,u*q with the voltages actually applied ud, uq, which are given by the projections of the reference voltages on a system of orthogonal axes rotated precisely by A9. The assessment of the actually applied voltages ud, uq, can be made by using the following rotation matrix: where udq and u*dq are column vectors in which the first element is the d-axis voltage, and the other element is the q-axis voltage.
[0045] At the subsequent step 104, the operating conditions, i.e. speed, current and temperature of the motor, are extracted from the steady states. It should be noted that different operating conditions can share the same values of speed and current (i.e., they can belong to the same steady state), but will then have a certain difference in temperature A9. However, for the reasons set out below, the operating conditions to be actually used for identifying the stator resistance and the flux of the magnets must belong to different steady states. By using two steady states, two equation systems (la), (lb) can be written, so that there will be four equations with four unknowns and the rank deficiency problem is solved without resorting to signal injection.
[0046] To extract the operating conditions, thermal transients must be taken into account. Indeed, the R-statistic algorithm can easily detect only quite fast transients, such as electrical or mechanical transients, and actually acts onto the current and speed of the motor to discard measurements obtained during transients of these quantities. However, the motor is also affected by thermal transients, which are much longer than the electrical and mechanical ones. In each of the steady states of the motor detected by the algorithm there can therefore be temperature variations that have considerable influence on two of the parameters to be identified (resistance and flux). This means that these two parameters vary during the same steady state if temperature changes. It is therefore appropriate to slice the data of each steady state into subsets relating to more or less uniform temperature ranges, in which the two parameters can be considered constant. The different operating conditions, uniquely identified by temperature, current and speed, are then obtained from said subsets. The data of these operating conditions are then used for identifying parameters at subsequent steps.
[0047] To take into account the slicing of the data of each steady state into data subsets, in order to carry out these subsequent steps it will be necessary to calculate the average values u, co, iq, 0, Dqof the various variables in the different subsets, where 0 is the temperature of the motor.
[0048] Once the operating conditions have been extracted, for each of them the identification of the stator inductance (step 105) is immediate, based on the relation (la) which, not considering Vdead which is a periodic function of the period with zero mean when id = 0, gives where j indicates the operating condition.
[0049] At this point, Vdead can be identified and therefore the non-linearity of the inverter can be compensated (step 106) by solving, for each operating condition, the minimization problem represented by the relation: where i, h and j indicate the 1thsteady state, the h*11subset and the jthoperating condition, respectively, and ksj ]-, and l<ej ]-, indicate the first sample and last sample of the subset h of the steady state I, respectively. Clearly, given relation (la), relation (4) represents the application of the method of least squares to the difference between the actually applied voltage ud and the reference voltage u*d- It should be noted that the effectiveness of research depends on the availability of samples within a period of DdVdead, which decreases as the engine speed increases because the sampling time Tsis fixed. At the subsequent step 107, identification of R and \| / m is carried out. As mentioned, the data of two operating conditions are used for this purpose, and therefore there will be two equations (lb).
[0050] The system of the two equations can be written in matrix form as: where the symbolsAand - above a quantity indicate the estimated value of such quantity and the mean value obtained as a result of splicing the data of each steady state into data subsets, and a and 0 indicate the two operating conditions that are used, which will be selected according to a criterion illustrated below. To solve the system, it is necessary that the two equations are linearly independent of each other, which involves that the determinant of the matrix of currents and speeds in the second member of the relation (5) is 0 and therefore the following condition is met
[0051] Meeting the condition (6) clearly requires that the operating conditions used belong to different steady states.
[0052] Under these conditions, by introducing (6) into the equation (lb) relating to each of the two operating conditions, R and \| / m (or, more precisely, the estimated values thereof R and ipm) can be obtained according to the relations:
[0053] The use of the above relations leads to estimation errors because the parameter values generally vary as the operating condition varies, and thus a different pair of resistance and flux values can be associated with each operating condition. Therefore, when parameters are to be identified in a certain operating condition a (main condition), the second operating condition 0 (auxiliary condition) is selected in such a way that the estimation error is minimal.
[0054] To take into account the dependency of the parameters on the operating conditions and the presence of errors in compensating the actuation delay and the non-linearity of the inverter, for each operating condition the equation (lb) can be re-written by introducing a term eUq representative of voltage errors, thereby obtaining
[0055] For simplicity, only one equation has been written, omitting indication of the operating condition.
[0056] As far as resistance is concerned, by taking into account (8) and breaking down the error term into two components ER and sR^j, which are due to the resistance value difference in the two operating conditions and the error of compensation of the q-axis voltage, respectively, the equations (7a) can be re-written as: i.e. the values of R are given by the actual value of R in the two operating conditions and two error terms. These in turn are given by:
[0057] Similarly, the equations for the values of i|rmcan be rewritten, which values will also be given by the actual value of \| / m in the two operating conditions and the two error terms due to the difference in the concatenated flux value in the two operating conditions and the error of compensation of the q-axis voltage.
[0058] From the equations (10) and the analogous equations for rpm it can be seen that the estimation errors depend on speed co and current iq. As mentioned, the two operating conditions are selected in such a way as to minimize the estimation error, and more precisely the error majorants. Referring to the estimation of R in the operating condition a, the error majorants are given by with where £R, ERMarethe values of the aforementioned error components, which values are obtained by using approximate estimates Ri, ipmi (i = a, P) of R and \| / m. These are calculated by using only temperature and speed data according to the relations where Ro and i|io indicate the approximate estimates of R and \| / m at speed 0 and temperature 20°C, ao is the temperature coefficient of copper resistivity, Po is the assumed quadratic frequency coefficient and SpMO is the assumed temperature coefficient of the permanent magnets. The use of precautionary values for £Uq, Po and «PM0 ensures that ERtot> ER + ^RM-
[0059] To obtain Ro, the operating conditions a, p are selected as follows: with the constraints: where lj jm an(j q jmare positive adjustment parameters. The relation (14) originates from the fact that Ro is calculated by correcting the estimate RQ of the resistance in the condition a. Therefore, it is convenient to find the condition a in such a way that the values Ra and RQ are close to each other, i.e. in such a way that temperature and speed in the equations (13) are close to 20°C and 0, respectively. The temperature correction can be obtained easily, because aO does not depend on the motor used, but this does not apply to the frequency correction, because po depends on the motor. To avoid adding precautionary factors to the relation (14), one only makes 65a close to 0 (first constraint for the relation (14)). In addition, errors £Raand fRjyj have to be minimised. The first of these two errors cannot be explicitly minimised, because the estimates of the parameters are unknown, but it can be made small by selecting the two operating conditions in such a way that |1 - r| in the first of the equations (12) is large (second constraint for the relation (14)). Instead, ERMis minimized by considering it as a target function. Once the conditions a, 0 are selected according to the relation (14), Ra is calculated according to the relation (7a). Finally, a temperature and frequency correction gives RQ:
[0060] With similar considerations, to obtain ijro the operating conditions a, 0 are selected as follows: with the constraints: where ®limis an adjustment parameter. Nocis the total number of operating conditions. Once the conditions a, 0 are selected according to the relation (16), r|rm0 is calculated according to the relation (7b). Finally, a temperature and frequency correction gives vpO
[0061] The majorant of the error of estimation of the concatenated flux of the rotor in the operating condition 0 can in turn be calculated with the procedure described for the stator resistance.
[0062] Once the approximate estimates have been obtained, error majorants can be calculated and explicitly minimised. Regarding the resistance R, given a main operating condition a for which this has to be determined, the auxiliary condition 0 will be selected in such a way as to minimize £Ratob namely with the constraints r < srl or r > Sr2 and £Rjtot <xR^a, where:
[0063] - sriand sr2 are constants that are arbitrarily selected on a case-by-case basis, the first constant being less than 1 and the second constant being more than 1, which avoid selecting r too close to 1 and thus effectively guarantee the linear independence of the two equations describing the model under the two operating conditions;
[0064] - xR instead is a parameter that is also arbitrarily selected on a case-by-case basis and has a value between 0 and 1, which requires that the error majorant be less than a given fraction of the approximate estimate of Ra.
[0065] The foregoing constraints, too, are met only in the case of operating conditions belonging to different steady states.
[0066] An estimate in the main operating condition may be rejected if there is no such auxiliary operating condition that the error majorant £Rtot issmaller than the given fraction of Rot.
[0067] Once the auxiliary condition p has been found, actual identification is obtained with the equation (7a).
[0068] A similar procedure will lead to the identification of \| / m-
[0069] The invention effectively overcomes the drawbacks of prior art. In fact, it offers a totally non-invasive method that requires neither experimental testing nor designing a new drive control scheme, but is based solely on the analysis of data generated by a standard drive during normal operation of the motor. Moreover, using the data of two operating conditions for identification allows finding parameter values that actually ensure minimisation of the estimation error.
[0070] Furthermore, it is clear from the foregoing description that the invention does not use the nameplate values of the parameters. First of all, this ensures greater applicability of the method. In fact, in industrial plants where motors of different sizes and different manufacturers are involved, it is not trivial to have data sheets available for all the types of motors present. This also makes the method more robust. In fact, the nameplate values are referred to measurements effected at the beginning of the lifecycle of the motor and in predetermined environmental and operating conditions, which will generally differ from the actual conditions of the application in which the motor operates. In addition, because of the wear of various components of the PMSM, the parameters will vary over time.
[0071] It is apparent that the above description has been given solely as a non-limiting example and that variations and modifications are possible without departing from the scope of protection of the invention as defined in the appended claims.
Claims
Claims1. Method of identifying the parameters of a permanent magnet synchronous motor, more particularly stator winding inductances and resistances and flux generated by the rotor magnets, the method providing for compensating the non-linearity of the inverter supplying the motor with power and of the actuation delays, characterised in that the steady states of the motor are identified (102) and the operating conditions, namely current, speed and temperature, of said steady states are extracted (104), and in that, for identifying the parameters, the data of two operating conditions of different steady states are used, which operating conditions are selected so as to minimise estimation errors for the parameter being identified.
2. Method according to claim 1, wherein the steady state identification (102) is performed by using the R-statistic algorithm, and the operating condition extraction (104) includes slicing the data of each steady state into data subsets relating to time intervals in which the motor temperature is such that temperature-dependent parameters of the motor are substantially constant, the parameter identification using, for each operating condition, the data of each of said subsets.
3. Method according to claim 2, wherein the compensation of the non-linearity of the inverter (106) includes identifying a voltage representative of the effects of such non-linearity, such an identification being carried out by minimising, for all samples of each operating condition of each steady state, the difference between a voltage actually applied to the motor and a reference voltage.
4. Method according to claim 1, wherein the minimisation of the estimation errors for the stator winding resistances and the flux generated by the rotor magnets includes minimising an error comprising a first error component due to the difference in the resistance and flux values in the two operating conditions, and a second error component due to an error in the compensation of the actuation delays by which the reference voltage is affected because of the inverter non-linearity.
5. Method according to claim 4, wherein the first error component is determined from a rough estimation of the stator winding resistances and the flux generated by the rotor magnets obtained by using predetermined values of motor temperature and rotor speed coefficients.
6. Method according to any preceding claim, wherein the minimisation of the estimation errors for the stator winding resistances and the flux generated by the rotor magnets includes imposing the constraints that the majorant of the estimation error is lower than a desired fraction of the respective motor parameter, and that the value of a parameter which allowsevaluation of the linear independence of the two equations describing a mathematical method of the motor in the two operating conditions selected and is linked to current and speed in the two conditions is not close to a unitary value denoting the lack of independence of said equations.
7. Method according to any preceding claim, for identifying the parameters of a permanent magnet synchronous motor employed in a device for capping containers.
8. Capping device actuated by a permanent magnet synchronous motor and controlled by a control system, characterised in that said control system comprises at least one processing unit programmed for implementing the method according to any preceding claim.
9. Computer program product loaded into the memory of a processing unit, in particular a processing unit belonging to a control system of a device for capping containers such as bottles and the like actuated by a permanent magnet synchronous motor, the product including portions of software code for implementing the method according to any of claims 1 to 7 when the product is executed on the processing unit.