CN113238492A - Parameter identification method for ARMAX model of vane motor - Google Patents
Parameter identification method for ARMAX model of vane motor Download PDFInfo
- Publication number
- CN113238492A CN113238492A CN202110398512.8A CN202110398512A CN113238492A CN 113238492 A CN113238492 A CN 113238492A CN 202110398512 A CN202110398512 A CN 202110398512A CN 113238492 A CN113238492 A CN 113238492A
- Authority
- CN
- China
- Prior art keywords
- model
- identification
- armax
- continuous rotation
- hydraulic servo
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 37
- 230000003190 augmentative effect Effects 0.000 claims abstract description 14
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012795 verification Methods 0.000 claims description 5
- 101000802640 Homo sapiens Lactosylceramide 4-alpha-galactosyltransferase Proteins 0.000 claims description 3
- 102100035838 Lactosylceramide 4-alpha-galactosyltransferase Human genes 0.000 claims description 3
- 238000004088 simulation Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 238000010408 sweeping Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B17/00—Systems involving the use of models or simulators of said systems
- G05B17/02—Systems involving the use of models or simulators of said systems electric
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
The invention discloses a parameter identification method for an ARMAX model of a blade continuous rotation electro-hydraulic servo motor, belonging to the field of parameter identification, and the identification method comprises the following steps: 1. establishing an ARMAX identification model of a blade continuous rotation electro-hydraulic servo motor system; 2. determining the amplitude, frequency and sweep frequency time of an input signal; 3. determining the system order and the noise order of the ARMAX discrete model by adopting a final prediction error order-fixing method; 4. optimizing an ARMAX identification model of the blade continuous rotation motor by using an augmented least square method; 5. comparing the output data of the identification model with the actually acquired output data; 6. and verifying the accuracy of the identification model.
Description
Technical Field
The invention relates to the field of system parameter identification, in particular to a method for identifying parameters of an ARMAX model of a blade continuous rotation electro-hydraulic servo motor.
Background
The flight attitude simulation turntable is semi-physical simulation equipment with important national defense strategic significance and economic value, and the quality of the performance of the flight attitude simulation turntable is directly related to the fidelity of an aircraft simulation result. The blade continuous rotation electro-hydraulic servo motor has the advantages of large transmission torque, high response speed, high control precision, wide speed regulation range, high dynamic position rigidity and steady-state speed rigidity, strong anti-interference capability and the like.
The research on the continuous rotation electrohydraulic servo motor of the blade needs to identify relevant parameters, and at present, the identification of the parameters of the continuous rotation electrohydraulic servo motor of the blade does not always have an accurate mathematical model due to more parameter variables and the irregularity of parameter change, so that the parameter identification of a continuous rotation motor system by adopting a discrete identification system model is easier to obtain an accurate system model. The augmented least square method is a more common identification algorithm in system identification, and for an electro-hydraulic servo system, because the state space variable of the electro-hydraulic servo system has more parameters to be identified and the estimation precision required by the system is higher, the calculated amount is large; in addition, the electro-hydraulic servo system is a nonlinear system essentially, and the influence of nonlinear factors needs to be considered when a complete mathematical model of the electro-hydraulic servo system is constructed; the traditional blade continuous rotation electro-hydraulic servo motor identification usually ignores the complexity influence on identification caused by the uncertainty and nonlinearity of a system, so that the problems of low parameter fitting degree and low identification accuracy exist.
Therefore, the identification method which can simplify the mathematical model, optimize the calculation logic and save the calculation time on the premise of ensuring the accuracy is considered.
Disclosure of Invention
The invention provides a parameter identification method for an ARMAX model of a blade continuous rotation electro-hydraulic servo motor, which is used for identifying the process of the parameters of the blade continuous rotation electro-hydraulic servo motor and carrying out online real-time identification and optimization on the parameters by using an augmented least square method. The method effectively overcomes the influence of system nonlinearity and uncertain factors, thereby effectively improving the control precision of the electro-hydraulic servo system.
The technical scheme of the invention is as follows:
a parameter identification method for an ARMAX model of a blade continuous rotation electro-hydraulic servo motor is disclosed, and the basic idea is as follows: firstly, selecting and preprocessing input signals, secondly, determining the order of the ARMAX discrete model by using a final prediction error order method, thirdly, optimizing the ARMAX identification model of the blade continuous rotation motor by using an augmented least square method, and finally, verifying the identification result. The method specifically comprises the following steps:
the method comprises the following steps: establishing an ARMAX model expression:
A(z-1)z(k)=z-dB(z-1)θ(k)+D(z-1)e(k) (1)
wherein D (z)-1) e (k) ═ epsilon (k); ε (k) is colored noise, D (z)-1) To generate output, θ (k) is the excitation signal, d is the delay times, e (k) is white noise, and { θ (k) }, { z (k) } are input and output for the identification model, respectively.
Step two: optimizing an ARMAX identification model of the continuous rotary motor by adopting an augmented least square method:
in the order of model na、nbAnd ndUnder the determined conditions, the parameter vector of the model is:
step three: unbiased estimation of β with augmented least squares estimation:
substituting the undetectable noise e (k-1), e (k-2), e (k-3) in the data vector h (k) to obtain:
step four: the augmented least squares recursion algorithm is as follows:
P(k)=[I-K(k)hT(k)]P(k-1) (8)
where Λ (k) is a positive definite weighting matrix, k (k) is a gain matrix, and:
step five: establishing a state space model of a blade continuous rotation electro-hydraulic servo motor system:
after the identification data is processed, the identification data of the first 80s is subjected to system identification, and the data of the second 20s is used as model verification data. According to the obtained model order, obtaining a three-order discrete model as follows:
the transfer function from the discrete transfer function is:
step six: and identifying system parameters and verifying identification results.
The invention is mainly characterized in that:
1. according to the parameter identification method for the blade continuous rotation electro-hydraulic servo motor, the ARMAX parameter of the continuous rotation electro-hydraulic servo motor can be updated on line by using the current sampling input and output;
2. and accurate ARMAX model identification parameters can be provided in time when the model parameters change.
Drawings
FIG. 1 is a diagram of an input signal of an identification system.
FIG. 2 is a block diagram of an ARMAX recognition model.
FIG. 3 is a comparison graph of the recognition result and the actual result.
Fig. 4 is a graph of actual output versus swept output.
Fig. 5 is a comparison graph of the verification segment.
Detailed Description
The control method of the present invention will be described in detail with reference to the accompanying drawings.
The method comprises the following steps: as shown in FIG. 1, the initial frequency of the input signal of the recognition system, i.e., the sweep frequency signal, is 0.2Hz, with a peak value of 2 °.
Step two: as shown in FIG. 2, a block diagram of an ARMAX identification model with an excitation signal of θ (k) and colored noise ε (k) added to the output, A (z)-1)z(k)=z-dB(z-1)θ(k)+D(z-1) e (k). Wherein D (z)-1) e (k) ═ epsilon (k), epsilon (k) is a linear system D (z) driven by white noise-1) And the output, D (z), produced-1) Does not affect the convergence of the transfer function. From fig. 2, the ARMAX model of the system dynamics is:
A(z-1)z(k)=z-dB(z-1)θ(k)+D(z-1)e(k) (1)
in the formula: { θ (k) }, { z (k) } are respectively input and output of the identification model, d is the delay times, and e (k) is white noise; ε (k) is colored noise.
Step three: determining the order of the model by using a final prediction error order-determining method, listing FPE values of different orders, and taking na=3,nb=1,nd=3。
Step four: optimizing an ARMAX identification model of the continuous rotary motor by adopting an augmented least square method:
in the order of model na、nbAnd ndHas been confirmedUnder certain conditions, the parameter vector of the model is:
step five: write (2) as augmented least squares format:
z(k)=hT(k)β(k)+e(k) (4)
step six: carrying out unbiased estimation on beta by using an augmented least square estimation method, and substituting the undetectable noise e (k-1), e (k-2) and e (k-3) in the data vector h (k) to obtain:
step seven: the augmented least squares recursion algorithm is as follows:
P(k)=[I-K(k)hT(k)]P(k-1) (9)
where Λ (k) is a positive definite weighting matrix, k (k) is a gain matrix, and:
step eight: establishing a state space model of the continuous rotation electro-hydraulic servo motor system:
after the identification data is processed, the identification data of the first 80s is subjected to system identification, and the data of the second 20s is used as model verification data. According to the obtained model order, obtaining a three-order discrete model as follows:
the transfer function from the discrete transfer function is:
step nine: and identifying system parameters and verifying identification results.
After the system parameters of the blade continuous rotation motor are identified according to the identification data, the comparison of the identification results shows that the fitting degree of the identification results and the system identification output results is 87.76%, as shown in fig. 3. The same sweep frequency signal is adopted for sweeping the frequency of the system after identification, the result after sweeping the frequency is compared with the actual identification result, and the comparison of the obtained sweep frequency result and the identification result shows that the fitting degree of the sweep frequency output and the identification output is 71.35 percent, as shown in fig. 4. Fig. 5 is a comparison graph of the verification segment.
Claims (7)
1. A parameter identification method for an ARMAX model of a blade continuous rotation electro-hydraulic servo motor is characterized by comprising the following steps:
the method comprises the following steps: establishing an ARMAX identification model of the blade continuous rotation electro-hydraulic servo motor system according to the identification data of the continuous rotation motor system;
step two: the selection of the input signal has great influence on the identification result, and an input signal capable of exciting a system must be selected;
step three: determining the system order and the noise order of the ARMAX discrete model by adopting a final prediction error order-fixing method;
step four: optimizing an ARMAX identification model of the continuous rotary motor by adopting an augmented least square method;
step five: and comparing the output data of the identification model with the actually acquired output data to verify the accuracy of the identification model.
2. The parameter identification method of the ARMAX model for the vane continuous rotation electro-hydraulic servo motor as claimed in claim 1, wherein the ARMAX identification model is established as follows:
A(z-1)z(k)=z-dB(z-1)θ(k)+D(z-1) e (k), wherein D (z)-1)e(k)=ε(k) (1)
In the formula: { θ (k) }, { z (k) } are respectively input and output of the identification model, d is the delay times, and e (k) is white noise; ε (k) is colored noise.
3. The parameter identification method for the ARMAX model of the blade continuous rotation electro-hydraulic servo motor as claimed in claim 1, wherein the augmented least squares model is established as follows:
in the order of model na、nbAnd ndUnder the determined conditions, the parameter vector of the model is:
4. the method for identifying the parameters of the ARMAX model of the blade continuous rotation electro-hydraulic servo motor as claimed in claim 3, wherein the β is estimated unbiased by using an augmented least squares estimation method:
substituting the undetectable noise e (k-1), e (k-2), e (k-3) in the data vector h (k) to obtain:
5. the method for identifying the parameters of the ARMAX model of the blade continuous rotation electro-hydraulic servo motor as claimed in claim 4, wherein the method comprises the following steps:
P(k)=[I-K(k)hT(k)]P(k-1) (8)
where Λ (k) is a positive definite weighting matrix, k (k) is a gain matrix, and:
6. the parameter identification method for the ARMAX model of the blade continuous rotation electro-hydraulic servo motor is characterized in that the establishment of a state space model of a blade continuous rotation electro-hydraulic servo motor system is carried out;
after the identification data is processed, the identification data of the first 80s is subjected to system identification, and the data of the second 20s is used as model verification data. According to the obtained model order, obtaining a three-order discrete model as follows:
the transfer function from the discrete transfer function is:
7. the method for identifying parameters of the ARMAX model of the vane continuous rotation electro-hydraulic servo motor as claimed in claim 1, wherein the analysis identifies a degree of fit of the result to the output result.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110398512.8A CN113238492A (en) | 2021-04-12 | 2021-04-12 | Parameter identification method for ARMAX model of vane motor |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110398512.8A CN113238492A (en) | 2021-04-12 | 2021-04-12 | Parameter identification method for ARMAX model of vane motor |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113238492A true CN113238492A (en) | 2021-08-10 |
Family
ID=77128169
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110398512.8A Pending CN113238492A (en) | 2021-04-12 | 2021-04-12 | Parameter identification method for ARMAX model of vane motor |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113238492A (en) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005242580A (en) * | 2004-02-25 | 2005-09-08 | Osaka Prefecture | Parameter estimation method, data prediction method, parameter estimation device, data prediction device, and computer program |
CN102540891A (en) * | 2012-01-17 | 2012-07-04 | 中冶南方工程技术有限公司 | Recursive extended least squares algorithm-based crystallizer ARMAX (Auto Regressive Moving Average Exogenous) model identification method |
CN103261987A (en) * | 2010-10-22 | 2013-08-21 | 斯奈克玛 | Method and device for monitoring a feedback loop of a variable-eometry actuator system of a jet engine |
JP2016151811A (en) * | 2015-02-16 | 2016-08-22 | 国立大学法人大阪大学 | Parameter identification method of system |
CN108879786A (en) * | 2018-08-15 | 2018-11-23 | 浙江运达风电股份有限公司 | The discrimination method and device of wind power generating set main component frequency and damping ratio |
CN111293686A (en) * | 2020-02-29 | 2020-06-16 | 上海电力大学 | ARMAX system identification-based real-time evaluation method for inertia of power system |
CN112627953A (en) * | 2020-12-22 | 2021-04-09 | 上海海事大学 | Ship SCR system ammonia injection amount control method based on ARMAX and MMPC |
-
2021
- 2021-04-12 CN CN202110398512.8A patent/CN113238492A/en active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2005242580A (en) * | 2004-02-25 | 2005-09-08 | Osaka Prefecture | Parameter estimation method, data prediction method, parameter estimation device, data prediction device, and computer program |
CN103261987A (en) * | 2010-10-22 | 2013-08-21 | 斯奈克玛 | Method and device for monitoring a feedback loop of a variable-eometry actuator system of a jet engine |
CN102540891A (en) * | 2012-01-17 | 2012-07-04 | 中冶南方工程技术有限公司 | Recursive extended least squares algorithm-based crystallizer ARMAX (Auto Regressive Moving Average Exogenous) model identification method |
JP2016151811A (en) * | 2015-02-16 | 2016-08-22 | 国立大学法人大阪大学 | Parameter identification method of system |
CN108879786A (en) * | 2018-08-15 | 2018-11-23 | 浙江运达风电股份有限公司 | The discrimination method and device of wind power generating set main component frequency and damping ratio |
CN111293686A (en) * | 2020-02-29 | 2020-06-16 | 上海电力大学 | ARMAX system identification-based real-time evaluation method for inertia of power system |
CN112627953A (en) * | 2020-12-22 | 2021-04-09 | 上海海事大学 | Ship SCR system ammonia injection amount control method based on ARMAX and MMPC |
Non-Patent Citations (3)
Title |
---|
孔祥臻 等: "伺服阀控液压马达系统的参数模型辨识", 机床与液压, vol. 36, no. 12, pages 105 - 106 * |
王晓晶 等: "Self-correcting wavelet neural network control of continuous rotary electro-hydraulic servo motor", HIGH TECHNOLOGY LETTERS, vol. 27, no. 1, pages 26 - 37 * |
王晓晶 等: "连续回转马达电液伺服系统辨识及控制", 哈尔滨工程大学学报, vol. 32, no. 8, pages 1045 - 1051 * |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106487297B (en) | A kind of PMSM parameter identification method based on Unscented kalman filtering algorithm | |
US10331810B2 (en) | Method for determining a model of an output quantity of a technical system | |
CN110826021B (en) | Robust identification and output estimation method for nonlinear industrial process | |
CN110276139B (en) | Induction motor key parameter obtaining method based on depth confidence strategy gradient algorithm | |
CN109194225B (en) | Online identification method for parameters of doubly-fed motor | |
CN110705041B (en) | EASI-based linear structure working modal parameter identification method | |
CN115733675A (en) | Distributed filtering method based on induction motor system | |
CN111506996A (en) | Self-adaptive identification method of turntable servo system based on identification error limitation | |
CN110968961A (en) | Parameter identification method for continuous rotation electro-hydraulic servo motor | |
CN109687792B (en) | Online identification optimization method for traction motor rotor parameters facing vector control system | |
CN117040339B (en) | Vector decoupling method, device and equipment for permanent magnet synchronous motor and storage medium | |
CN113238492A (en) | Parameter identification method for ARMAX model of vane motor | |
CN110084324B (en) | Kalman filtering parameter self-adaptive updating method based on extreme learning machine | |
CN116085245A (en) | Online compressor performance prediction method and system based on OS-ELM | |
CN114337427B (en) | Rotary inertia identification method of recursive least square method with forgetting factor | |
CN116976209A (en) | Data-driven unmanned aerial vehicle system frequency domain identification method and system based on joint decision | |
CN113688474B (en) | Electric steering engine simulation modeling method based on parameter identification | |
CN114136538B (en) | Temperature modeling method for pressure sensor calibration device based on random variation decibel leaf learning | |
CN113486523B (en) | Global identification method for linear variable parameter vibration system | |
CN115102444A (en) | Self-adaptive integral sliding mode prediction control method for permanent magnet synchronous motor | |
CN118381406B (en) | Vector decoupling control method for permanent magnet synchronous motor and related equipment | |
CN111682532A (en) | Excitation system uncompensated phase-frequency characteristic online modeling method and storage medium | |
CN114744941B (en) | Permanent magnet synchronous motor permanent magnet demagnetization online monitoring method and system | |
CN113325709B (en) | Adaptive control method under variable load condition of rotary table | |
CN118611518A (en) | Raccoon optimization algorithm-based permanent magnet synchronous motor parameter identification method |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20210810 |