CN114325164B - Multi-fault diagnosis method for single-phase three-level rectifier - Google Patents

Multi-fault diagnosis method for single-phase three-level rectifier Download PDF

Info

Publication number
CN114325164B
CN114325164B CN202111417465.3A CN202111417465A CN114325164B CN 114325164 B CN114325164 B CN 114325164B CN 202111417465 A CN202111417465 A CN 202111417465A CN 114325164 B CN114325164 B CN 114325164B
Authority
CN
China
Prior art keywords
matrix
fault
state
state variable
voltage
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.)
Active
Application number
CN202111417465.3A
Other languages
Chinese (zh)
Other versions
CN114325164A (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.)
Hefei University of Technology
Original Assignee
Hefei University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Hefei University of Technology filed Critical Hefei University of Technology
Priority to CN202111417465.3A priority Critical patent/CN114325164B/en
Publication of CN114325164A publication Critical patent/CN114325164A/en
Application granted granted Critical
Publication of CN114325164B publication Critical patent/CN114325164B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Testing Of Short-Circuits, Discontinuities, Leakage, Or Incorrect Line Connections (AREA)

Abstract

The invention provides a multi-fault diagnosis method for a single-phase three-level rectifier, and belongs to the field of power fault diagnosis. The method comprises the following steps: establishing a hybrid logic dynamic model, establishing a state space expression comprising an actuator fault and two voltage sensing fault systems, carrying out state augmentation transformation to construct a new system, carrying out primary coordinate transformation, carrying out matrix transformation, carrying out secondary coordinate transformation, calculating the state quantity of a reduced-order system, designing a self-adaptive reduced-order sliding-mode observer, and giving a self-adaptive diagnosis threshold value to carry out multi-fault diagnosis. The reduced-order sliding mode observer provided by the invention does not need a system to provide a very accurate dynamic model, only needs to reasonably design a sliding mode surface by utilizing a tracking error of a track, and can adaptively switch an approach rate, thereby accelerating the approach speed, reducing buffeting of sliding mode motion, improving the diagnosis accuracy, and simultaneously diagnosing 3 different faults.

Description

Multi-fault diagnosis method for single-phase three-level rectifier
Technical Field
The invention relates to the field of power fault diagnosis, in particular to a multi-fault diagnosis method for a single-phase three-level rectifier.
Background
In an electric traction drive system, a single-phase three-level rectifier plays an increasingly important role, wherein the single-phase three-level rectifier is one of the core components of the traction drive system of a high-speed train. Aiming at a single-phase three-level rectifier, on one hand, the output end of the rectifier is provided with two voltage sensors, and the aging and damage of the voltage sensors can be caused by the long-time running of a high-speed train and the interference of the external environment, so that the measurement feedback data of the sensors is abnormal, and serious casualties and property loss are caused; on the other hand, the single-phase three-level rectifier can also have actuator faults at the same time, the actuator faults generated in the actual transportation process can cause abnormal functions of a traction transmission system of the high-speed train, serious traffic accidents are caused, and meanwhile, the difficulty of fault diagnosis can be greatly increased.
The diagnosis methods for the voltage sensor fault and the actuator fault of the single-phase three-level rectifier mainly comprise the following two diagnosis methods:
1. the method is based on a data driving method, and the corresponding fault diagnosis is realized by acquiring fault data information and analyzing and processing the data information; corresponding papers and patents such as "a data-drive a vehicle to an activator and sensor fault detection, isolation and estimation in discrete-time linear systems" and "a data-driven transmission sensor fault diagnosis method" (application publication No. CN 202011201104.0) do not need to establish an accurate mathematical model, but need a large amount of data bases, and need to perform complex data processing on the bases, and the algorithms are complex, and have large workload and high difficulty.
2. The method is based on an analytic model, compares the information quantity obtained from the mathematical model of the actual system to be diagnosed with the actual measurement by knowing and establishing the mathematical model of the actual system, and carries out fault diagnosis by analyzing the residual error; meanwhile, if multi-fault diagnosis is performed, fault decoupling is performed, a corresponding mathematical model or observer is established to obtain an estimated quantity, and a residual error analysis is performed after the estimated quantity is compared with an actual measured quantity, wherein the corresponding papers are as follows:
a Simultaneous project activator and sensor fault evaluation for non-linear Lipschitz systems aims at performing complex fault decoupling on faults of an actuator and a sensor to form two subsystems and then respectively designing an observer for the two subsystems to perform fault diagnosis;
Sensor-Fault-Estimation-Based Tolerant Control for Single-Phase Two-Level PWM Rectifier in Electric conduction System, which aims at multi-Fault diagnosis of faults of an actuator and a Sensor, although the Fault System is subjected to step reduction processing, the diagnosis object is a Single-Phase Two-Level Rectifier instead of a mainstream Single-Phase three-Level Rectifier;
when a Sliding Mode Observer is established to diagnose the Fault of a Sensor, the adopted Approach law also generates buffeting when the Approach rate is slow in the observation process, so that the observed value is inaccurate;
the paper only expands the Fault amount of one voltage sensor, namely, only carries out Fault diagnosis of one voltage sensor when carrying out multi-Fault diagnosis of faults of an actuator and a sensor.
In summary, the single-phase three-level rectifier plays an extremely important role in the electric traction drive system, and the existing diagnosis methods for the sensor fault and the actuator fault of the single-phase three-level rectifier have disadvantages, so that the technical problems to be solved in the whole research field are solved by solving the disadvantages of the prior art.
Disclosure of Invention
The invention aims to provide a multi-fault diagnosis method for a single-phase three-level rectifier, aiming at the problems in the background art. Specifically, the problem that an accurate system model is difficult to establish is solved by designing a reduced-order sliding-mode observer and using proper sliding-mode control, but the method can generate jitter characteristics in the control process, and further designs a self-adaptive approach rate, and the jitter characteristics in the control process are reduced by the self-adaptive switching approach rate, so that the accuracy of system fault diagnosis is further improved; meanwhile, a reduced-order sliding mode observer is used, the harsh precondition established by the observer is overcome, and the fault diagnosis of two voltage sensors and one actuator is carried out simultaneously; finally, during fault diagnosis, the bounded voltage disturbance quantity and the bounded external disturbance which are possibly generated during measurement of the voltage sensor are also considered, so that the accuracy of the diagnosis method is improved, and the diagnosis precision is improved.
In order to achieve the aim, the invention provides a multi-fault diagnosis method for a single-phase three-level rectifier, and a circuit topology related to the method comprises a network side voltage source U s Network side equivalent inductor L s And net side equivalent resistance R s Rectifier bridge, two identical support capacitors C d1 ,C d2 A DC side load and two identical voltage sensors; support capacitor C d1 And a support capacitor C d2 After being connected in series, the capacitor C is connected in parallel between a direct current positive bus P and a direct current negative bus Q1 of a direct current side load d1 And a support capacitor C d1 The contact point of (2) is marked as a direct current bus midpoint O; two identical voltage sensors are respectively marked as a voltage sensor 1 and a voltage sensor 2, and the voltage sensor 1 is connected with a supporting capacitor C d1 At both ends of which the voltage sensor 2 is connected to the supporting capacitor C d2 Both ends of (a);
the rectifier bridge is divided into two-phase bridge arms, and the two-phase bridge arms are connected with the direct-current side load in parallel; recording two-phase bridge arms as bridge arms k, wherein k is a bridge sequence, and k = a, b; in the two-phase bridge arms, each phase of bridge arm comprises 4 switching tubes with reverse connection diodes and two clamping diodes, namely a rectifier bridge comprises 8 switching tubes with reverse connection diodes and 4 clamping diodes, and the 8 switching tubes form an actuator of a single-phase three-level rectifier; marking 8 switch tubes as switch tubes V γ denotes a serial number of the switching tube, γ =1,2,3,4, and 4 clamping diodes are denoted as clamping diodes D ckρ ρ is the serial number of the clamp diode, ρ =1,2; in each of the two-phase arms, a switching tube V k1 Switch tube V k2 Switch tube V k3 And a switch tube V k4 Are sequentially connected in series, wherein, the switch tube V k2 And a switching tube V k3 Is marked as the input point tau of the rectifier bridge k K = a, b; in each of the two-phase arms, a clamping diode D ck1 The cathode of the switch tube is connected with the switch tube V k1 And a switching tube V k2 Between, the clamping diode D ck1 Anode of (2) is connected to a clamping diode D ck2 Cathode of (2), clamping diode D ck2 Anode of the switch tube is connected with the switch tube V k3 And a switching tube V k4 And clamping diode D ck1 And a clamping diode D ck2 The connecting point of the direct current bus is connected with the midpoint O of the direct current bus;
the network side equivalent inductance L s One end of is connected with an input point tau of the rectifier bridge a The other end is connected with the equivalent resistance R of the network side in sequence s Grid side voltage source U s The other end of the series network side voltage source is connected with an input point tau of the rectifier bridge b
The multi-fault diagnosis method comprises the steps of fault diagnosis of an actuator of a single-phase three-level rectifier and weak fault diagnosis of a multi-voltage sensor, and specifically comprises the following steps:
step 1, establishing a hybrid logic dynamic model of a single-phase three-level rectifier, and calculating a phase voltage U of an input end of a rectifier bridge ab Is estimated value of
Figure BDA0003371171010000059
Sampling the current of the network side and recording the current of the network side as the current i of the network side s Sampling support capacitor C d1 And a support capacitor C d2 And is denoted as DC voltage u 1 ,u 2 Sampling the DC side voltage U dc (ii) a Establishing a mixed logic dynamic model of the single-phase three-level rectifier, and calculating the phase voltage U of the input end of the rectifier bridge ab Is estimated value of
Figure BDA0003371171010000051
The input end phase voltage U of the rectifier bridge ab For input point tau of rectifier bridge a And the input point tau of the rectifier bridge b A voltage in between;
the expression of the hybrid logic dynamic model of the single-phase three-level rectifier is as follows:
Figure BDA0003371171010000052
Figure BDA0003371171010000053
wherein the content of the first and second substances,
Figure BDA0003371171010000054
is an estimate of the voltage of the a-phase,
Figure BDA0003371171010000055
is an estimate of the b-phase voltage, T h Is a DC voltage u 1 Mixed logic dynamic function of, T l Is a DC voltage u 2 A hybrid logical dynamic function of (a);
the input end phase voltage U of the rectifier bridge ab Is estimated value of
Figure BDA0003371171010000056
The expression of (a) is:
Figure BDA0003371171010000057
step 2, establishing a multi-fault-containing single-phase three-level rectifier system state space expression, and recording the expression as a multi-fault system 1, wherein the expression is as follows:
Figure BDA0003371171010000058
wherein x (t) is a state variable of the multi-fault system 1 and is marked as a primary state variable x (t),
Figure BDA0003371171010000061
x (t) is n 1 Dimensional state variables, noted
Figure BDA0003371171010000062
t is a variable of time and is,
Figure BDA0003371171010000063
the first derivative of the primary state variable x (t),
Figure BDA0003371171010000064
wherein
Figure BDA0003371171010000065
Is a network side current i s The derivative of (a) is determined,
Figure BDA0003371171010000066
are respectively DC voltage u 1 ,u 2 A derivative of (a); u (t) is the input of the multi-fault system 1 and is recorded as primary system input u (t), u (t) = u s U (t) is n 2 Dimensional state variables, noted
Figure BDA0003371171010000067
f a (t) actuator failure for multiple failure System 1, denoted as actuator failure f a (t),f a (t) is n 3 Dimension state variable, as
Figure BDA0003371171010000068
Eta (t) is the harmonic disturbance quantity of the multi-fault system 1 and is recorded as harmonic disturbance eta (t), and eta (t) is n 4 Dimension state variable, as
Figure BDA0003371171010000069
y (t) is the output of the multi-fault system 1 and is recorded as the system output y (t), and y (t) is n 5 Dimension state variable, as
Figure BDA00033711710100000610
f s (t) Voltage sensor Fault for multiple Fault System 1, noted System Voltage sensor Fault f s (t),f s (t) is n 6 Dimension state variable, as
Figure BDA00033711710100000611
Wherein
Figure BDA00033711710100000612
For a fault in voltage sensor 1, it is noted as a fault f in voltage sensor 1 u1 (t),f u1 (t) is n 7 Dimensional state variables, noted
Figure BDA00033711710100000613
f u2 (t) is a voltage sensor 2 failure, denoted as voltage sensor 2 failure f u2 (t),f u2 (t) is n 8 Dimensional state variables, noted
Figure BDA00033711710100000614
A 1 The state matrix for the multiple fault system 1, denoted as primary state matrix A 1
Figure BDA00033711710100000615
Wherein L is equivalent inductance L at network side s R is the equivalent resistance R of the network side s Resistance value of C 1 ,C 2 Are respectively a support capacitor C d1 And a support capacitor C d2 The capacitance value of (a); b is 1 The input matrix for the multiple fault system 1 is marked as a primary input matrix B 1
Figure BDA0003371171010000071
G 1 For actuator failure f a (t) the coefficient matrix is recorded as a primary actuator fault coefficient matrix,
Figure BDA0003371171010000072
N 1 the coefficient matrix of harmonic disturbance eta (t) is recorded as a primary disturbance coefficient matrix N 1
Figure BDA0003371171010000073
C y1 For multiple fault systems1 output matrix, denoted primary output matrix
Figure BDA0003371171010000074
Figure BDA0003371171010000075
F 1 For system voltage sensor failure f s (t) coefficient matrix, which is recorded as primary voltage sensor fault coefficient matrix F 1
Figure BDA0003371171010000076
System voltage sensor fault f s (t) actuator failure f a (t) and the harmonic disturbance η (t) are both continuous and bounded functions satisfying the time variable t, denoted | | f s (t)||≤r s ,||f a (t)||≤r a ,||η(t)||≤r η Wherein | | | f s (t) | | is f s Norm of (t, | f a (t) | | is f a The norm of (t), where | | | η (t) | | is the norm of η (t), r s For system voltage sensor fault f s Upper limit of (t), r a For actuator failure f a Upper limit of (t), r η Is the supremum of the harmonic disturbance eta (t), and r s 、r a And r η Are all normal numbers, and are recorded as r s >0,r a >0,r η >0;
Step 3, introducing system state variable augmentation transformation to the multi-fault system 1, and defining the following augmentation matrixes and vectors, wherein the expression is as follows:
Figure BDA0003371171010000077
Figure BDA0003371171010000081
Figure BDA0003371171010000082
a new system after the multi-fault system 1 is subjected to system state variable amplification transformation can be obtained, and the new system is marked as a multi-fault system 2, wherein the expression of the state space of the new system is as follows:
Figure BDA0003371171010000083
wherein z is 1 (t) is the state variable of the multiple fault system 2, noted as the secondary state variable z 1 (t),z 1 (t) is n 1 +n 3 +n 6 Dimensional space vector, denoted as
Figure BDA0003371171010000084
Is a quadratic state variable z 1 (ii) the derivative of (t),
Figure BDA0003371171010000085
wherein
Figure BDA0003371171010000086
Is a failure f of the voltage sensor 1 u1 (ii) the derivative of (t),
Figure BDA0003371171010000087
for failure of the voltage sensor 2
Figure BDA0003371171010000088
A derivative of (a); f. of 1 (t) is a fault variable of the multi-fault system 2, and is recorded as a system fault quantity f 1 (t);
E 1 As the second state variable derivative
Figure BDA00033711710100000812
Is recorded as a first derivative coefficient matrix E 1 In a
Figure BDA0003371171010000089
In the expression (2) of (a),
Figure BDA00033711710100000810
represents n 1 The dimension-unit matrix is a matrix of the dimension units,
Figure BDA00033711710100000811
represents n 3 A dimension unit matrix;
A 2 the state matrix for the multiple fault system 2, denoted as the quadratic state matrix A 2
B 2 The input matrix for the multiple fault system 2, denoted as the quadratic input matrix B 2
F 2 Is the amount of system failure f 1 (t) coefficient matrix, denoted as quadratic failure coefficient matrix F 2 In a
Figure BDA0003371171010000091
In the expression of (a) in (b),
Figure BDA0003371171010000092
represents n 7 The dimension-unit matrix is a matrix of dimension units,
Figure BDA0003371171010000093
represents n 8 A dimension unit matrix;
N 2 the coefficient matrix of harmonic disturbance eta (t) in the multi-fault system 2 is recorded as a secondary disturbance coefficient matrix N 2
Figure BDA0003371171010000094
Output matrix for multiple fault system 2
Figure BDA0003371171010000095
Is recorded as a quadratic output matrix
Figure BDA0003371171010000096
Step 4, carrying out one-time coordinate transformation on the multi-fault system 2
Step 4.1, outputting the secondary output matrix
Figure BDA0003371171010000097
Transposing, recording the transposed output matrix as
Figure BDA0003371171010000098
Handle
Figure BDA0003371171010000099
The matrix decomposition into the product of an orthogonal matrix and an upper triangular matrix is recorded as
Figure BDA00033711710100000910
U is an orthogonal matrix and Q is an upper triangular matrix; let P = U again T ,J=Q T ,U T Is the transpose of the orthogonal matrix U, Q T Is the transpose matrix of the upper triangular matrix Q, then the obtained
Figure BDA00033711710100000911
Wherein P satisfies P T =P -1 ,P T Is the transpose of the matrix P, P -1 Is the inverse matrix of the matrix P, and takes P as a primary coordinate transformation matrix;
step 4.2, introduce coordinate transformation z 2 (t)=Pz 1 (t), the multi-fault system 2 can be converted into a new system through coordinate transformation, the new system is recorded as a multi-fault system 3, and the state space expression of the new system is as follows:
Figure BDA00033711710100000912
wherein z is 2 (t) is the state variable of the multiple fault system 3, noted as the third state variable z 2 (t);
Figure BDA00033711710100000913
Is a cubic state variable z 2 (t) a derivative of;
E 2 is a third state variable derivative
Figure BDA00033711710100000914
Coefficient matrix of (2), noted as second derivativeCoefficient matrix E 2 ,E 2 =E 1 P T ;A 3 The state matrix of the multi-fault system 3 is marked as a cubic state matrix A 3 ,A 3 =A 2 P T ;B 3 The input matrix of the multi-fault system 3 is recorded as a cubic input matrix B 3 ;F 3 For faults f of the multiple fault system 3 1 (t) coefficient matrix, denoted as cubic failure coefficient matrix F 3 ;N 3 The coefficient matrix of the harmonic disturbance eta (t) of the multi-fault system 3 is recorded as a third-order disturbance coefficient matrix N 3 (ii) a J is the output matrix of the multiple fault system 3, denoted as the cubic output matrix J, and is denoted as the block matrix, i.e., J = (J) 1 0), wherein J 1 Left blocking submatrix of J, J 1 Belongs to n 5 Dimensional space, is denoted as
Figure BDA0003371171010000101
And J 1 Is a non-singular matrix;
step 5, matrix transformation is carried out on the multi-fault system 3
Step 5.1, the second derivative coefficient matrix E 2 Conversion into a block matrix, i.e. E 2 =(E 21 E 22 ) In which E 21 Is a coefficient matrix E of the second derivative 2 Left blocking submatrix of (E) 21 Is (n) 1 +n 3 +n 6 )×n 5 Dimensional space, is denoted as
Figure BDA0003371171010000102
E 22 Is a coefficient matrix E of the second derivative 2 Right blocking submatrix of (E) 22 Is (n) 1 +n 3 +n 6 )×(n 1 +n 3 +n 6 -n 5 ) Dimensional space, is denoted as
Figure BDA0003371171010000103
Then designing a transformation matrix xi, recording as the matrix transformation matrix xi,
Figure BDA0003371171010000104
wherein xi 11 Being the upper part of the submatrix xi of the transform matrix xi 11 Is n 5 ×(n 1 +n 3 +n 6 ) Dimensional space, is denoted as
Figure BDA0003371171010000105
Ξ 21 A lower blocking submatrix of the transform matrix xi, xi 21 Belong to (n) 1 +n 3 +n 6 -n 5 )×(n 1 +n 3 +n 6 ) Dimensional space, is denoted as
Figure BDA0003371171010000106
And xi 11 Satisfies E 22 T Ξ 11 T =0, solution formula E 22 T Ξ 11 T =0, then may be xi 11 Wherein the matrix E 22 T Is a coefficient matrix E of the second derivative 2 Right blocking submatrix E 22 Transpose of (2), matrix xi 11 T An upper blocking submatrix xi being a transform matrix xi 11 Transposing; xi 21 =E 22 T (E 22 E 22 T ) -1 Matrix (E) 22 E 22 T ) -1 Is a matrix E 22 E 22 T The inverse of (c);
step 5.2, the left and right state space expressions of the multi-fault system 3 are multiplied by the transformation matrix xi to obtain a new transformed system, the new system is marked as a multi-fault system 4, and the state space expressions are as follows:
Figure BDA0003371171010000111
wherein z is 3 (t) is the state variable of the multiple fault system 4, noted as the quartile state variable z 3 (t);
Figure BDA0003371171010000112
Is a quartic state variable z 3 (t) a derivative of;
E 3 is the fourth derivative of the state variable
Figure BDA0003371171010000113
Is recorded as a coefficient matrix of the third derivative E 3
Figure BDA0003371171010000114
Wherein E 31 Is a coefficient matrix E of the third derivative 3 Upper left blocking sub-matrix of (E) 32 Is a coefficient matrix E of the third derivative 3 The lower left block sub-matrix of (a),
Figure BDA0003371171010000115
is n 1 +n 3 +n 6 -n 5 A dimension unit matrix; a. The 4 The state matrix for the multiple fault system 4, denoted as the fourth state matrix A 4
Figure BDA0003371171010000116
Wherein A is 411 Is a four-times state matrix A 4 Upper left blocking sub-matrix of A 412 Is a four-times state matrix A 4 Upper right blocking sub-matrix of, A 421 Is a four-times state matrix A 4 Left lower blocking submatrix of, A 422 Is a four-times state matrix A 4 The lower right blocking submatrix; b is 4 The input matrix for the multiple fault system 4 is denoted as the fourth input matrix B 4
Figure BDA0003371171010000117
Wherein B is 41 Is a four-input matrix B 4 Upper blocking matrix of, B 42 Is a four-input matrix B 4 A lower block matrix of (a); f 4 For faults f of the multiple fault system 4 1 (t) coefficient matrix, denoted as Quaternary failure coefficient matrix F 4
Figure BDA0003371171010000118
Wherein F 41 Is a four-fault coefficient matrix F 41 Upper block matrix of F 42 Is a four-fault coefficient matrix F 4 A lower blocking matrix of (a); n is a radical of 4 For sensor failureThe coefficient matrix of harmonic disturbance eta (t) of the system 4 is recorded as a fourth disturbance coefficient matrix N 4
Figure BDA0003371171010000119
Wherein N is 41 Is a four-times disturbance coefficient matrix N 41 Upper blocking submatrix of (2), N 42 Is a four-times disturbance coefficient matrix N 4 A lower block submatrix of (a);
step 6, carrying out secondary coordinate transformation on the multi-fault system 4
Step 6.1, a coordinate transformation matrix T is designed and recorded as a quadratic coordinate transformation matrix T, and T is expressed as a block matrix, namely
Figure BDA0003371171010000121
Wherein therein
Figure BDA0003371171010000122
Is n 5 The dimension-unit matrix is a matrix of dimension units,
Figure BDA0003371171010000123
is n 1 +n 3 +n 6 -n 5 Dimension unit matrix, L is a free matrix, denoted as free matrix L, which belongs to (n) 1 +n 3 +n 6 -n 5 )×n 5 Dimensional space, as
Figure BDA0003371171010000124
Step 6.2 introduction of coordinate transformation z 4 (t)=Tz 3 (t) a transformed new system is obtained, and the new system is marked as a multi-fault system 5, and the state space expression of the new system is as follows:
Figure BDA0003371171010000125
wherein z is 4 (t) is the state variable of the multiple fault system 5, noted as the quintic state variable z 4 (t) for the fifth state variable z 4 (t) blocking, i.e.
Figure BDA0003371171010000126
z 41 (t) is a five-state variable z 4 (t) upper block subvector z 41 (t),z 41 (t) is n 5 Dimension vector, z 42 (t) is a five-state variable z 4 (t) lower block subvector z 42 (t),z 42 (t) is n 1 +n 3 +n 6 -n 5 A dimension vector;
Figure BDA0003371171010000127
is a quintic state variable z 4 Derivative of (t), i.e.
Figure BDA0003371171010000128
Is z 4 (t) an upper block vector z 41 (ii) the derivative of (t),
Figure BDA0003371171010000129
is z of the upper block vector 42 (t) derivative of; is transformed by the above coordinate z 4 (t)=Tz 3 (t) it can be seen that,
Figure BDA00033711710100001210
E 4 fifth order state variable derivative
Figure BDA00033711710100001211
A coefficient matrix of (a), and
Figure BDA00033711710100001212
A 5 the state matrix for the multi-fault system 5 is denoted as the five-time state matrix A 5
Figure BDA00033711710100001213
Wherein A is 511 Is a quintic state matrix A 5 Upper left blocking sub-matrix of A 512 Is a quintic state matrix A 5 Upper right blocking sub-matrix of (A) 521 Is a quintic state matrix A 5 Left lower blocking submatrix of (A) 522 Is a quintic state matrix A 5 The lower right blocking submatrix; t is -1 An inverse matrix which is a quadratic coordinate transformation matrix T; b is 5 The input matrix for the multiple fault system 5 is denoted as the quintic input matrix B 5
Figure BDA0003371171010000131
Wherein B is 51 For five inputs of matrix B 5 Upper blocking submatrix of (2), B 52 For five inputs of matrix B 5 A lower block submatrix of (a); f 5 The coefficient matrix for the fault F (t) of the multiple fault system 5 is noted as the quintic fault coefficient matrix F 5
Figure BDA0003371171010000132
Wherein F 51 Is a quintic fault coefficient matrix F 5 Upper blocking submatrix of F 52 Is a quintic fault coefficient matrix F 5 A lower block submatrix of (a); n is a radical of 5 The coefficient matrix of harmonic disturbance eta (t) of the multi-fault system 5 is recorded as a quintic disturbance coefficient matrix N 5
Figure BDA0003371171010000133
Wherein N is 51 For a quintic disturbance coefficient matrix N 5 Upper blocking submatrix of (2), N 52 For a quintic disturbance coefficient matrix N 5 Lower blocking submatrix of (J) * Is the output matrix of the multiple fault system 5 and knows the following relationship:
Figure BDA0003371171010000134
step 7, designing a self-adaptive reduced-order sliding mode observer
Step 7.1, calculating the state quantity theta of the order-reduced system according to the multi-fault system 5, and recording the state quantity theta as the state quantity theta of the order-reduced system, wherein the calculation formula is as follows:
Θ=(LE 31 +E 32 -L)z 41 (t)+z 42 (t)
and solving a state space expression of the reduced-order system according to the reduced-order system state quantity theta, wherein the state space expression is as follows:
Figure BDA0003371171010000135
wherein the content of the first and second substances,
Figure BDA0003371171010000136
is the derivative of the state quantity theta of the reduced-order system, and delta is the state expression z of the reduced-order system 41 The coefficient matrix of (t) is marked as a reduced system coefficient matrix delta, and the expression is as follows:
Δ=-(LA 512 +A 522 )(LE 51 +E 52 -L)+L(A 511 -A 512 L)+A 511 -A 522 L
step 7.2, designing a self-adaptive reduced order sliding mode observer aiming at a state space expression of a reduced order system, wherein the expression is as follows:
Figure BDA0003371171010000141
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003371171010000142
is the estimated value of the state quantity theta of the reduced order system and is recorded as the estimated value of the state quantity theta of the reduced order system
Figure BDA0003371171010000143
Is an estimate of the state quantity of the reduced order system
Figure BDA0003371171010000144
Is a sliding mode gain matrix, Γ = (LF) 51 +F 52 ) And v is the approach law,
Figure BDA0003371171010000145
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003371171010000146
in order to be able to vary the parameter 1,
Figure BDA0003371171010000147
tanh () is a hyperbolic tangent function,
Figure BDA0003371171010000148
is variable parameter 2, and
Figure BDA0003371171010000149
σ is a variable parameter 3, and σ is an element (0, 1), θ is a variable parameter 4, and θ > 1, ρ is a variable parameter 5, β > 0, Ψ is a positive definite symmetric matrix, χ is a diagonal matrix,
Figure BDA00033711710100001410
epsilon is a constant number term, epsilon > 1 Θ (t) is the error of the estimation,
Figure BDA00033711710100001411
let S = e Θ (t), wherein S is a sliding mode surface of the designed adaptive reduced-order sliding mode observer;
step 7.3, obtaining a free matrix L by solving a Lyapunov equation, wherein the expression of the Lyapunov equation is as follows:
(LA 412 +A 422 ) T P+P(LA 412 +A 422 )=-I
wherein (LA) 412 +A 422 ) T Is LA 412 +A 422 P is a positive definite symmetric matrix, and I is a unit matrix;
step 8, fault diagnosis of the actuator and the sensor is carried out
Step 8.1, sampling the value of the system output y (t) of the multiple fault system 1 in step 1 and substituting the sampled value into z in step 6.2 41 (t)=z 31 (t)=J 1 -1 y (t), the five state variables z can be obtained 4 (t) upper block subvector z 41 (t) value, again according to the known five state variable z 4 (t) upper block vector z 41 Value of (t), Θ = (LE) in step 7.1 31 +E 32 -L)z 41 (t)+z 42 (t), step 7.2
Figure BDA00033711710100001412
And the error e is estimated in step 7.2 Θ (t) =0, and z can be obtained 42 (t), the expression of which is:
Figure BDA0003371171010000151
then the five state variables z are obtained 4 (t) upper block subvector z 41 (t) value and quintic state variable z 4 (t) lower block subvector z 42 (t) value substituted into step 6.2
Figure BDA0003371171010000152
The state variable z can be calculated five times 4 (t) a value;
step 8.2, starting from the multiple fault system 2, according to a coordinate transformation z 2 (t)=Pz 1 (t) and quadratic coordinate transformation z 4 (t)=Tz 3 (t) and the matrix transformation does not change the cubic state variable z in the sensor failure system 3 in step 5 2 (t) obtaining z 3 (t)=z 2 (t) thereby obtaining z 4 (t)=TPz 1 (t) calculating a secondary state variable z by inverse operation of the matrix 1 (t)=P T T -1 z 4 (t), and the five state variables z calculated in step 8.1 are further added 4 (t) substitution of formula z 1 (t)=P T T -1 z 4 (t) obtaining a secondary state variable z 1 (t);
Step 8.3, mixing
Figure BDA0003371171010000153
Is recorded as a primary state variable estimation value,
Figure BDA0003371171010000154
is recorded as the estimated value of the actuator fault
Figure BDA0003371171010000155
Is recorded as the fault estimation value of the voltage sensor 1
Figure BDA0003371171010000156
Is recorded as the fault estimation value of the voltage sensor 2
Figure BDA0003371171010000157
Calculating a state variable estimate
Figure BDA0003371171010000158
Actuator fault estimation
Figure BDA0003371171010000159
Voltage sensor fault 1 estimation
Figure BDA00033711710100001510
And voltage sensor 2 fault estimation
Figure BDA00033711710100001511
The specific calculation formula is as follows:
Figure BDA00033711710100001512
Figure BDA00033711710100001513
Figure BDA00033711710100001514
Figure BDA00033711710100001515
step 8.4, diagnosing faults of the actuator, the voltage sensor 1 and the voltage sensor 2, and giving a self-adaptive diagnosis threshold value T th
Defining fault location quantity Z 1
Figure BDA00033711710100001516
And positioned as follows:
if Z is 1 =0, no actuator failure has occurred in the multiple fault system 1;
if Z is 1 =1, the multi-fault system 1 has an actuator fault;
defining fault location Z 2
Figure BDA0003371171010000161
And positioned as follows:
if Z is 2 =0, no voltage sensor 1 fault occurred in the multiple fault system 1;
if Z is 2 =1, multiple fault system 1 has voltage sensor 1 fault;
defining fault location Z 3
Figure BDA0003371171010000162
And positioned as follows:
if Z is 3 =0, no voltage sensor 2 failure occurred in the multiple fault system 1;
if Z is 3 =1, the multiple fault system 1 has a voltage sensor 2 fault.
Preferably, the DC voltage u in step 1 1 Hybrid logic dynamic function T of h And a DC voltage u 2 Hybrid logic dynamic function T of l The calculation process of (c) is as follows:
recording the switching function of the k-phase bridge arm as S k K = a, b, then:
Figure BDA0003371171010000163
Figure BDA0003371171010000164
recording the pulse control signal of the switching tube as v K = a, b, γ =1,2,3,4, then the k-phase bridge arm switching function S k =v k1 v k2 -v k3 v k4
Preferably, the diagnostic adaptive threshold T of step 8.4 th The expressions are respectively as follows:
T th =E 1112 )
wherein E is 1 Is a constant 1, Γ 1 Is a constant number 2, and Γ 1 ∈(1,2),ζ 1 ζ is a bounded external perturbation 2 Is a bounded voltage perturbation.
Compared with the prior art, the invention has the beneficial effects that:
1. the fault state and the fault value of the system can be accurately estimated by designing the reduced-order sliding-mode observer, and the defect of poor robustness for weak fault diagnosis in the existing analytical model-based method is overcome;
2. by designing a self-adaptive sliding mode approach rate and a self-adaptive switching approach rate, buffeting influence in observation of the reduced-order sliding mode observer is reduced, the approach rate is accelerated, and the accuracy of system fault diagnosis is further improved;
3. through the increase of system state variables, the fault diagnosis of two voltage sensors and one actuator is carried out simultaneously, and the types of fault diagnosis are increased;
4. by designing the self-adaptive diagnosis threshold, bounded voltage disturbance possibly existing in the fault value of the voltage sensor is considered, and the accuracy of fault diagnosis is improved.
Drawings
FIG. 1 is a topology of a single phase three level rectifier in an example of the invention;
FIG. 2 is a schematic diagram of a single phase three level rectifier fault diagnostic method of the present invention;
FIG. 3 is a flow chart of a single phase three level rectifier fault diagnostic method of the present invention;
FIG. 4 shows an actuator failure f in this embodiment a (t) actuator estimation value
Figure BDA0003371171010000171
And an adaptive threshold T th A simulation graph of (a);
FIG. 5 shows the failure of the voltage sensor 1 in this example
Figure BDA0003371171010000172
Voltage sensor
1 fault estimation
Figure BDA0003371171010000173
And an adaptive threshold T th A simulated oscillogram of (c);
FIG. 6 shows a failure f of the voltage sensor 2 in this example u2 (t), voltage sensor 2 Fault estimation value
Figure BDA0003371171010000174
And an adaptive threshold T th A simulated waveform diagram of (c).
Detailed Description
Fig. 1 is a topology diagram of a single phase three level rectifier in an embodiment of the present invention. It can be seen from the figure that the circuit topology according to the invention comprises a network-side voltage source U s Network side equivalent inductance s and network side equivalent resistance R s Rectifier bridge, two identical support capacitors C d1 ,C d2 A direct current side load and two identical voltage sensors; support capacitor C d1 And a support capacitor C d2 After being connected in series, the capacitor C is connected in parallel between a direct current positive bus P and a direct current negative bus Q1 of a direct current side load d1 And a support capacitor C d2 The contact point of (2) is marked as a direct current bus midpoint O; two identical voltage sensors are respectively marked as a voltage sensor 1 and a voltage sensor 2, and the voltage sensor 1 is connected with a supporting capacitor C d1 At both ends of which the voltage sensor 2 is connected to the supporting capacitor C d2 At both ends of the same. In FIG. 1, S V1 Is a voltage sensor 1, S V2 For the voltage sensor 2, two sensors are used for the measurement of the voltage value.
The rectifier bridge is divided into two-phase bridge arms, and the two-phase bridge arms are connected with the direct-current side load in parallel; recording two-phase bridge arms as bridge arms k, wherein k is a bridge sequence, and k = a, b; in two-phase bridge arms, each phase of bridge arm comprises 4 switching tubes with reverse-connection diodes, two clamping diodes, namely a rectifier bridge comprises 8 switching tubes with reverse-connection diodes and 4 clamping diodesThe 8 switching tubes form an actuator of the single-phase three-level rectifier; marking 8 switch tubes as switch tubes V γ denotes a serial number of the switching tube, γ =1,2,3,4, and 4 clamping diodes are denoted as clamping diodes D ckρ ρ is the serial number of the clamp diode, ρ =1,2; in each of the two-phase arms, a switching tube V k1 And a switch tube V k1 Switch tube V k3 Switch tube V k4 Are sequentially connected in series, wherein, the switch tube V k2 And a switching tube V k3 Is recorded as the rectifier bridge input point tau k K = a, b; in each of the two-phase arms, a clamping diode D ck1 The cathode of the switch tube is connected with the switch tube V k1 And a switching tube V k2 Between, the clamping diode D ck1 Anode of (2) is connected to a clamping diode D ck2 Cathode of (2), clamping diode D ck2 Anode of the switch tube is connected with the switch tube V k3 And a switching tube V k4 And clamping diode D ck1 And a clamping diode D ck2 The connecting point of the direct current bus is connected with the midpoint O of the direct current bus.
The network side equivalent inductor L s One end of is connected with an input point tau of the rectifier bridge a The other end is connected with the equivalent resistance R of the network side in sequence s Grid side voltage source U s The other end of the series network side voltage source is connected with an input point tau of the rectifier bridge b
Fig. 2 is a schematic diagram of a multi-fault diagnosis method for a single-phase three-level rectifier according to the present invention, fig. 3 is a flowchart of the multi-fault diagnosis method for the single-phase three-level rectifier according to the present invention, and as can be seen from fig. 2 to fig. 3, the multi-fault diagnosis method includes an actuator fault diagnosis for the single-phase three-level rectifier and a multi-voltage sensor weak fault diagnosis, and specifically includes the following steps:
step 1, establishing a hybrid logic dynamic model of a single-phase three-level rectifier, and calculating a phase voltage U of an input end of a rectifier bridge ab Is estimated value of
Figure BDA0003371171010000191
Sampling the current of the network side and recording the current of the network side as the current i of the network side s Sampling ofSupport capacitor C d1 And a support capacitor C d2 And is denoted as DC voltage u 1 ,u 2 Sampling the DC side voltage U dc (ii) a Establishing a mixed logic dynamic model of the single-phase three-level rectifier, and calculating the phase voltage U of the input end of the rectifier bridge ab Is estimated value of
Figure BDA0003371171010000192
The input end phase voltage U of the rectifier bridge ab For input point tau of rectifier bridge a And the input point tau of the rectifier bridge b The voltage in between.
The expression of the hybrid logic dynamic model of the single-phase three-level rectifier is as follows:
Figure BDA0003371171010000193
Figure BDA0003371171010000194
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003371171010000195
is an estimate of the voltage of the a-phase,
Figure BDA0003371171010000196
is an estimate of the b-phase voltage, T h Is a direct current voltage u 1 Mixed logical dynamic function of (2), T l Is a DC voltage u 2 Mixed logical dynamic functions of (1).
The input end phase voltage U of the rectifier bridge ab Is estimated value of
Figure BDA0003371171010000197
The expression of (a) is:
Figure BDA0003371171010000198
in this embodiment, the straight lineCurrent voltage u 1 Hybrid logic dynamic function T of h And a DC voltage u 2 Mixed logic dynamic function T of l The calculation process of (2) is as follows:
recording the switching function of the k-phase bridge arm as S k K = a, b, then:
Figure BDA0003371171010000201
Figure BDA0003371171010000202
recording the pulse control signal of the switching tube as v K = a, b, γ =1,2,3,4, then the k-phase bridge arm switching function S k =v k1 v k2 -v k3 v k4
Step 2, establishing a multi-fault-containing single-phase three-level rectifier system state space expression, and recording the expression as a multi-fault system 1, wherein the expression is as follows:
Figure BDA0003371171010000203
wherein x (t) is a state variable of the multi-fault system 1 and is marked as a primary state variable x (t),
Figure BDA0003371171010000204
x (t) is n 1 Dimension state variable, as
Figure BDA0003371171010000205
t is a variable of time and is,
Figure BDA0003371171010000206
the first derivative of the primary state variable x (t),
Figure BDA0003371171010000207
wherein
Figure BDA0003371171010000208
Is a net side current i s The derivative of (a) of (b),
Figure BDA0003371171010000209
are respectively DC voltage u 1 ,u 2 A derivative of (a); u (t) is the input of the multiple fault system 1 and is recorded as one system input u (t), u (t) = u s U (t) is n 2 Dimensional state variables, noted
Figure BDA00033711710100002010
f a (t) actuator failure for multiple failure System 1, denoted as actuator failure f a (t),f a (t) is n 3 Dimensional state variables, noted
Figure BDA00033711710100002011
Eta (t) is the harmonic disturbance quantity of the multi-fault system 1 and is recorded as harmonic disturbance eta (t), and eta (t) is n 4 Dimension state variable, as
Figure BDA00033711710100002012
y (t) is the output of the multi-fault system 1 and is denoted as the system output y (t), and y (t) is n 5 Dimensional state variables, noted
Figure BDA00033711710100002013
f s (t) Voltage sensor Fault for multiple Fault System 1, noted System Voltage sensor Fault f s (t),f s (t) is n 6 Dimensional state variables, noted
Figure BDA0003371171010000211
Wherein
Figure BDA0003371171010000212
For a failure of voltage sensor 1, note as a failure f of voltage sensor 1 u1 (t),f u1 (t) is n 7 Dimension state variable, as
Figure BDA0003371171010000213
f u2 (t) is a voltage sensor 2 failure, denoted as voltage sensor 2 failure f u2 (t),f u2 (t) is n 8 Dimension state variable, as
Figure BDA0003371171010000214
A 1 The state matrix for the multiple fault system 1, denoted as primary state matrix A 1
Figure BDA0003371171010000215
Wherein L is equivalent inductance L at network side s R is the equivalent resistance R of the network side s Resistance value of C 1 ,C 2 Are respectively a support capacitor C d1 And a support capacitor C d2 The capacitance value of (a); b is 1 The input matrix for the multiple fault system 1 is denoted as primary input matrix B 1
Figure BDA0003371171010000216
G 1 For actuator failure f a (t) the coefficient matrix is recorded as a primary actuator fault coefficient matrix,
Figure BDA0003371171010000217
N 1 the coefficient matrix of harmonic disturbance eta (t) is recorded as a primary disturbance coefficient matrix N 1
Figure BDA0003371171010000218
C y1 The output matrix of the multi-fault system 1 is recorded as a primary output matrix
Figure BDA0003371171010000219
Figure BDA00033711710100002110
F 1 For system voltage sensor fault f s (t) coefficient matrix, denoted as primary voltage sensor failure coefficient matrix F 1
Figure BDA00033711710100002111
System voltage sensor fault f s (t) actuator failure f a (t) and harmonic disturbances η (t) are both continuous and bounded functions satisfying a time variable t, denoted as | | f s (t)||≤r s ,||f a (t)||≤r a ,||η(t)||≤r η Wherein | | | f s (t) | | is f s Norm of (t, | f a (t) | | is f a The norm of (t), where | | | η (t) | | is the norm of η (t), r s For system voltage sensor fault f s Upper limit of (t), r a For actuator failure f a Upper limit of (t), r η Is a supremum of harmonic disturbance η (t), and r s 、r a And r η Are all normal numbers and are recorded as r s >0,r a >0,r η >0。
In the present embodiment, R =0.34 Ω, L =2.2 × 10 -3 H,C 1 =16×10 -3 F,C 2 =16×10 -3 F,
Figure BDA0003371171010000221
n 1 =3,n 2 =1,n 3 =1,n 4 =1,n 5 =5,n 6 =2,n 7 =1,n 8 =1,
Figure BDA0003371171010000222
f u1 (t)=0,t<10s,
Figure BDA0003371171010000223
t≥10s,
f u2 (t)=0,t<10s,
Figure BDA0003371171010000224
t≥10s。
System voltage sensor fault f s (t) actuator failure f a (t) and harmonic disturbance η (t) are satisfied with respect toContinuous and bounded function of the amount of time t, denoted | | f s (t)||≤r s ,||f a (t)||≤r a ,||η(t)||≤r η Wherein | | | f s (t) | | is f s Norm of (t, | f a (t) | | is f a The norm of (t), where | | | η (t) | | is the norm of η (t), r s For system voltage sensor fault f s Upper limit of (t), r a For actuator failure f a Upper limit of (t), r η Is the supremum of the harmonic disturbance eta (t), and r s 、r a And r η Are all normal numbers, and are recorded as r s >0,r a >0,r η >0。
Step 3, introducing system state variable augmentation transformation to the multi-fault system 1, and defining the following augmentation matrixes and vectors, wherein the expression is as follows:
Figure BDA0003371171010000225
Figure BDA0003371171010000231
Figure BDA0003371171010000232
a new system of the multi-fault system 1 after system state variable augmentation transformation can be obtained, and the new system is marked as a multi-fault system 2, and the state space expression of the new system is as follows:
Figure BDA0003371171010000233
wherein z is 1 (t) is the state variable of the multiple fault system 2, noted as the secondary state variable z 1 (t),z 1 (t) is n 1 +n 3 +n 6 Dimensional space vector, is
Figure BDA0003371171010000234
Is a quadratic state variable z 1 (ii) the derivative of (t),
Figure BDA0003371171010000235
wherein
Figure BDA0003371171010000236
Is a fault f of the voltage sensor 1 u1 (ii) the derivative of (t),
Figure BDA0003371171010000237
for failure of the voltage sensor 2
Figure BDA0003371171010000238
A derivative of (a); f. of 1 (t) is a fault variable of the multi-fault system 2, and is recorded as a system fault quantity f 1 (t);
E 1 As the second state variable derivative
Figure BDA0003371171010000239
Is recorded as a first derivative coefficient matrix E 1 In a
Figure BDA00033711710100002310
In the expression of (a) in (b),
Figure BDA00033711710100002311
represents n 1 The dimension-unit matrix is a matrix of dimension units,
Figure BDA00033711710100002312
represents n 3 A dimension unit matrix;
A 2 the state matrix for the multi-fault system 2 is denoted as a quadratic state matrix A 2
B 2 The input matrix for the multiple fault system 2, denoted as the quadratic input matrix B 2
F 2 Is the amount of system failure f 1 (t) coefficient matrix, denoted as quadratic failure coefficient matrix F 2 In a
Figure BDA0003371171010000241
In the expression of (a) in (b),
Figure BDA0003371171010000242
represents n 7 The dimension-unit matrix is a matrix of dimension units,
Figure BDA0003371171010000243
represents n 8 A dimension unit matrix;
N 2 the coefficient matrix of harmonic disturbance eta (t) in the multi-fault system 2 is marked as a secondary disturbance coefficient matrix N 2
Figure BDA0003371171010000244
Output matrix for multiple fault system 2
Figure BDA0003371171010000245
Is recorded as a quadratic output matrix
Figure BDA0003371171010000246
The parameters in this example are as follows:
Figure BDA0003371171010000247
Figure BDA0003371171010000248
Figure BDA0003371171010000249
Figure BDA00033711710100002410
step 4, carrying out one-time coordinate transformation on the multi-fault system 2
Step 4.1, outputting the secondary output matrix
Figure BDA00033711710100002411
Transposing, recording the transposed output matrix as
Figure BDA0003371171010000251
Handle
Figure BDA0003371171010000252
The matrix is decomposed into the product of an orthogonal matrix and an upper triangular matrix, and the product is recorded as
Figure BDA0003371171010000253
U is an orthogonal matrix and Q is an upper triangular matrix; let P = U again T ,J=Q T ,U T Is the transpose of the orthogonal matrix U, Q T Is the transpose of the upper triangular matrix Q, then the result is
Figure BDA0003371171010000254
Wherein P satisfies P T =P -1 ,P T Is the transpose of the matrix P, P -1 Is the inverse of matrix P, which is denoted as the primary coordinate transformation matrix.
Step 4.2, introducing coordinate transformation z 2 (t)=Pz 1 (t), the multi-fault system 2 can be transformed into a new system through coordinate transformation, the new system is recorded as a multi-fault system 3, and the state space expression is as follows:
Figure BDA0003371171010000255
wherein z is 2 (t) is the state variable of the multiple fault system 3, noted as the third state variable z 2 (t);
Figure BDA0003371171010000256
Is a cubic state variable z 2 (t) derivative of;
E 2 is the third order state variable derivative
Figure BDA0003371171010000257
Is denoted as the second derivative coefficient matrix E 2 ,E 2 =E 1 P T ;A 3 The state matrix of the multi-fault system 3 is marked as a cubic state matrix A 3 ,A 3 =A 2 P T ;B 3 The input matrix of the multi-fault system 3 is recorded as a cubic input matrix B 3 ;F 3 For faults f of multiple-fault systems 3 1 (t) coefficient matrix, denoted as cubic failure coefficient matrix F 3 ;N 3 The coefficient matrix of the harmonic disturbance eta (t) of the multi-fault system 3 is recorded as a third-order disturbance coefficient matrix N 3 (ii) a J is the output matrix of the multiple fault system 3, denoted as the cubic output matrix J, and is denoted as the block matrix, i.e., J = (J) 1 0), wherein J 1 Left blocking submatrix of J, J 1 Belongs to n 5 Dimensional space, is denoted as
Figure BDA0003371171010000258
And J 1 Is a non-singular matrix.
In the present embodiment, it is preferred that,
Figure BDA0003371171010000261
Figure BDA0003371171010000262
Figure BDA0003371171010000263
Figure BDA0003371171010000264
Figure BDA0003371171010000265
Figure BDA0003371171010000266
step 5, matrix transformation is carried out on the multi-fault system 3
Step 5.1, coefficient matrix E of second derivative 2 Conversion into a block matrix, i.e. E 2 =(E 21 E 22 ) In which E 21 Is a coefficient matrix E of the second derivative 2 Left blocking submatrix of (E) 21 Is (n) 1 +n 3 +n 6 )×n 5 Dimensional space, is denoted as
Figure BDA0003371171010000271
E 22 Is a coefficient matrix E of the second derivative 2 Right blocking submatrix, E 22 Is (n) 1 +n 3 +n 6 )×(n 1 +n 3 +n 6 -n 5 ) Dimensional space, is denoted as
Figure BDA0003371171010000272
Then designing a transformation matrix xi, which is recorded as a matrix transformation matrix xi,
Figure BDA0003371171010000273
wherein xi 11 Being an upper blocking submatrix of the transform matrix xi, xi 11 Is n 5 ×(n 1 +n 3 +n 6 ) Dimensional space, is denoted as
Figure BDA0003371171010000274
Ξ 21 Being a lower blocking submatrix of the transform matrix xi, xi 21 Belong to (n) 1 +n 3 +n 6 -n 5 )×(n 1 +n 3 +n 6 ) Dimensional space, as
Figure BDA0003371171010000275
And xi 11 Satisfies E 22 T Ξ 11 T =0, solution of formula E 22 T Ξ 11 T =0, then may be xi 11 Wherein the matrix E 22 T Is a coefficient matrix E of the second derivative 2 Right blocking submatrix E 22 Transpose of (2), matrix xi 11 T An upper blocking submatrix xi being a transform matrix xi 11 Transposing; xi 21 =E 22 T (E 22 E 22 T ) -1 Matrix (E) 22 E 22 T ) -1 Is a matrix E 22 E 22 T The inverse of (c).
Step 5.2, the left and right state space expressions of the multi-fault system 3 are multiplied by the transformation matrix xi to obtain a new transformed system, the new system is marked as a multi-fault system 4, and the state space expressions are as follows:
Figure BDA0003371171010000276
wherein z is 3 (t) is the state variable of the multiple fault system 4, noted as the quartile state variable z 3 (t);
Figure BDA0003371171010000277
Is a quartic state variable z 3 (t) derivative of;
E 3 is the fourth derivative of the state variable
Figure BDA0003371171010000278
Is noted as a coefficient matrix of the third derivative E 3
Figure BDA0003371171010000279
Wherein E 31 Is a coefficient matrix E of the third derivative 3 Upper left blocking sub-matrix of (E) 32 Is a coefficient matrix E of the third derivative 3 The lower left block sub-matrix of (a),
Figure BDA00033711710100002710
is n 1 +n 3 +n 6 -n 5 A dimension unit matrix; a. The 4 Is a state matrix of the multi-fault system 4,is denoted as a quartic state matrix A 4
Figure BDA00033711710100002711
Wherein A is 411 Is a four-times state matrix A 4 Upper left blocking sub-matrix of, A 411 Is a four-times state matrix A 4 Upper right blocking sub-matrix of (A) 421 Is a four-times state matrix A 4 Left lower blocking submatrix of, A 422 Is a four-times state matrix A 4 The lower right blocking submatrix; b is 4 The input matrix for the multiple fault system 4 is denoted as the fourth input matrix B 4
Figure BDA0003371171010000281
Wherein B is 41 Is a four-input matrix B 4 Upper block matrix of, B 42 Is a four-input matrix B 4 A lower block matrix of (a); f 4 For faults f of the multiple fault system 4 1 (t) coefficient matrix, denoted as Quaternary failure coefficient matrix F 4
Figure BDA0003371171010000282
Wherein F 41 Is a four-fault coefficient matrix F 41 Upper block matrix of F 42 Is a four-fault coefficient matrix F 4 A lower block matrix of (a); n is a radical of 4 The coefficient matrix for the harmonic disturbance η (t) of the sensor fault system 4 is noted as the fourth-order disturbance coefficient matrix N 4
Figure BDA0003371171010000283
Wherein N is 41 Is a four-times disturbance coefficient matrix N 41 Upper blocking submatrix of (2), N 42 Is a four-times disturbance coefficient matrix N 4 The lower block submatrix of (1).
In the present embodiment, it is preferred that,
Figure BDA0003371171010000284
Figure BDA0003371171010000285
Figure BDA0003371171010000286
Figure BDA0003371171010000291
step 6, carrying out secondary coordinate transformation on the multi-fault system 4
Step 6.1 design coordinate transformation matrix T, record as quadratic coordinate transformation matrix T, and express T as block matrix, that is
Figure BDA0003371171010000292
Wherein therein
Figure BDA00033711710100002912
Is n 5 The dimension-unit matrix is a matrix of the dimension units,
Figure BDA0003371171010000293
is n 1 +n 3 +n 6 -n 5 Dimension unit matrix, L is a free matrix, denoted as free matrix L, which belongs to (n) 1 +n 3 +n 6 -n 5 )×n 5 Dimensional space, is denoted as
Figure BDA0003371171010000294
Step 6.2 introduction of coordinate transformation z 4 (t)=Tz 3 (t) a transformed new system is obtained, and the new system is marked as a multi-fault system 5, and the state space expression of the new system is as follows:
Figure BDA0003371171010000295
wherein z is 4 (t) is the state variable of the multiple fault system 5, noted as the quintic state variable z 4 (t) for the quintic state variable z 4 (t) Is divided into blocks, i.e.
Figure BDA0003371171010000296
z 41 (t) is a quintic state variable z 4 (t) upper block subvector z 41 (t),z 41 (t) is n 5 Dimension vector, z 42 (t) is a five-state variable z 4 (t) lower block subvector z 42 (t),z 42 (t) is n 1 +n 3 +n 6 -n 5 A dimension vector;
Figure BDA0003371171010000297
is a five-fold state variable z 4 Derivative of (t), i.e.
Figure BDA0003371171010000298
Is z 4 (t) upper block vector z 41 (ii) the derivative of (t),
Figure BDA0003371171010000299
is z of the upper block vector 42 (t) derivative of; is transformed by the above coordinate z 4 (t)=Tz 3 (t) it can be seen that,
Figure BDA00033711710100002910
E 4 is the fifth derivative of the state variable
Figure BDA00033711710100002911
A coefficient matrix of (a), and
Figure BDA0003371171010000301
A 5 the state matrix for the multi-fault system 5 is denoted as the five-time state matrix A 5
Figure BDA0003371171010000302
Wherein A is 511 Is a quintic state matrix A 5 Upper left blocking sub-matrix of, A 512 Is a quintic state matrix A 5 Upper right blocking sub-matrix of, A 521 Is a quintic state matrix A 5 Is divided into left lower blocksSub-matrix, A 522 Is a quintic state matrix A 5 The lower right blocking submatrix; t is -1 An inverse matrix of the quadratic coordinate transformation matrix T; b 5 The input matrix for the multiple fault system 5 is marked as five-times input matrix B 5
Figure BDA0003371171010000303
Wherein B is 51 For five inputs of matrix B 5 Upper blocking submatrix of (2), B 52 For five inputs of matrix B 5 A lower block submatrix of (a); f 5 The coefficient matrix for the fault F (t) of the multiple fault system 5 is noted as the quintic fault coefficient matrix F 5
Figure BDA0003371171010000304
Wherein F 51 Is a quintic fault coefficient matrix F 5 Upper blocking submatrix of F 52 Is a quintic fault coefficient matrix F 5 A lower blocking submatrix of (a); n is a radical of hydrogen 5 The coefficient matrix of harmonic disturbance eta (t) of the multi-fault system 5 is recorded as a quintic disturbance coefficient matrix N 5
Figure BDA0003371171010000305
Wherein N is 51 For a quintic disturbance coefficient matrix N 5 Upper blocking submatrix of (2), N 52 For a quintic disturbance coefficient matrix N 5 Lower blocking submatrix, J * Is the output matrix of the multiple fault system 5, and the following relationship can be known:
Figure BDA0003371171010000306
in the present embodiment, it is preferred that,
Figure BDA0003371171010000307
Figure BDA0003371171010000311
Figure BDA0003371171010000312
step 7, designing a self-adaptive reduced-order sliding mode observer
Step 7.1, calculating the state quantity theta of the order-reduced system according to the multi-fault system 5, and recording the state quantity theta as the state quantity theta of the order-reduced system, wherein the calculation formula is as follows:
Θ=(LE 31 +E 32 -L)z 41 (t)+z 42 (t)
and solving a state space expression of the reduced-order system according to the reduced-order system state quantity theta, wherein the state space expression is as follows:
Figure BDA0003371171010000313
wherein, the first and the second end of the pipe are connected with each other,
Figure BDA0003371171010000314
is the derivative of the state quantity theta of the reduced-order system, and delta is the state expression z of the reduced-order system 41 The coefficient matrix of (t) is marked as a reduced system coefficient matrix delta, and the expression is as follows:
Δ=-(LA 512 +A 522 )(LE 51 +E 52 -L)+L(A 511 -A 512 L)+A 511 -A 522 L
step 7.2, designing a self-adaptive reduced order sliding mode observer aiming at a state space expression of a reduced order system, wherein the expression is as follows:
Figure BDA0003371171010000315
wherein the content of the first and second substances,
Figure BDA0003371171010000321
is the estimated value of the state quantity theta of the reduced order system and is recorded as the estimated value of the state quantity theta of the reduced order system
Figure BDA0003371171010000322
Is an estimate of the state quantity of the reduced order system
Figure BDA0003371171010000323
Is a sliding mode gain matrix, Γ = (LF) 51 +F 52 ) V is the approach law,
Figure BDA0003371171010000324
wherein the content of the first and second substances,
Figure BDA0003371171010000325
in order to be able to vary the parameter 1,
Figure BDA0003371171010000326
tanh () is a hyperbolic tangent function,
Figure BDA0003371171010000327
is variable parameter 2, and
Figure BDA0003371171010000328
σ is a variable parameter 3, and σ is an element (0, 1), θ is a variable parameter 4, and θ > 1, ρ is a variable parameter 5, β > 0, Ψ is a positive definite symmetric matrix, χ is a diagonal matrix,
Figure BDA0003371171010000329
epsilon is a constant number term, epsilon > 1 Θ (t) is the error of the estimation,
Figure BDA00033711710100003210
let S = e Θ (t), wherein S is a sliding mode surface of the designed adaptive reduced-order sliding mode observer;
step 7.3, obtaining a free matrix L by solving a Lyapunov equation, wherein the expression of the Lyapunov equation is as follows:
(LA 412 +A 422 ) T P+P(LA 412 +A 422 )=-I
wherein (LA) 412 +A 422 ) T Is LA 412 +A 422 P is a positive definite symmetric matrix, and I is a unit matrix.
In the present embodiment, it is preferred that,
Figure BDA00033711710100003211
σ=0.1,ρ=2,θ=20,
Figure BDA00033711710100003212
step 8, fault diagnosis of the actuator and the sensor is carried out
Step 8.1, sampling the value of the system output y (t) of the multiple fault system 1 in step 1 and substituting the sampled value into z in step 6.2 41 (t)=z 31 (t)=J 1 -1 y (t), the state variable z can be obtained five times 4 (t) upper block subvector z 41 (t) value, again according to the known five state variable z 4 (t) upper block vector z 41 Value of (t), Θ = (LE) in step 7.1 31 +E 32 -L)z 41 (t)+z 42 (t), step 7.2
Figure BDA0003371171010000331
And the error e is estimated in step 7.2 Θ (t) =0, and z can be obtained 42 (t) the expression is:
Figure BDA0003371171010000332
then the five state variables z are obtained 4 (t) upper block subvector z 41 (t) value and quintic state variable z 4 (t) lower block subvector z 42 (t) value substituted into step 6.2
Figure BDA0003371171010000333
The five state variable z can be calculated 4 (t) a value;
step 8.2, starting from the multiple fault system 2, according to a coordinate transformation z 2 (t)=Pz 1 (t) and quadratic coordinate transformation z 4 (t)=Tz 3 (t) and the matrix transformation does not change the cubic state variable z in the sensor fault system 3 in step 5 2 (t) obtaining z 3 (t)=z 2 (t) thereby obtaining z 4 (t)=TPz 1 (t) calculating a secondary state variable z by inverse operation of the matrix 1 (t)=P T T -1 z 4 (t), and the five state variables z calculated in step 8.1 are further added 4 (t) substitution of formula z 1 (t)=P T T -1 z 4 (t) obtaining a secondary state variable z 1 (t);
Step 8.3, mixing
Figure BDA0003371171010000334
Is recorded as a primary state variable estimation value,
Figure BDA0003371171010000335
recording as the estimated value of the actuator fault
Figure BDA0003371171010000336
Is recorded as the fault estimation value of the voltage sensor 1
Figure BDA0003371171010000337
Is recorded as the fault estimation value of the voltage sensor 2
Figure BDA0003371171010000338
Calculating a primary state variable estimate
Figure BDA0003371171010000339
Actuator fault estimation
Figure BDA00033711710100003310
Voltage sensor fault 1 estimation
Figure BDA00033711710100003311
And voltage sensor 2 fault estimation
Figure BDA00033711710100003312
The specific calculation formula is as follows:
Figure BDA00033711710100003313
Figure BDA00033711710100003314
Figure BDA00033711710100003315
Figure BDA00033711710100003316
step 8.4, diagnosing faults of the actuator, the voltage sensor 1 and the voltage sensor 2, and giving a self-adaptive diagnosis threshold value T th
In this embodiment, the adaptive threshold T is diagnosed th The expressions are respectively as follows:
T th =E 1112 )
wherein E is 1 Is a constant 1, Γ 1 Is a constant number 2, and Γ 1 ∈(1,2),ζ 1 For a bounded external disturbance, ζ 2 Is a bounded voltage perturbation. Specifically, in the present embodiment, E 1 =0.01,Γ 1 =1.01,ζ 1 =0.02sin(10t),
Figure BDA0003371171010000347
Defining fault location quantity Z 1
Figure BDA0003371171010000341
And positioned as follows:
if Z is 1 =0, no actuator failure occurred in the multiple fault system 1;
if Z is 1 =1, the multi-fault system 1 has an actuator fault;
defining fault location quantity Z 2
Figure BDA0003371171010000342
And positioned as follows:
if Z is 2 =0, no voltage sensor 1 failure has occurred in the multiple fault system 1;
if Z is 2 =1, multiple fault system 1 has voltage sensor 1 fault;
defining fault location quantity Z 3
Figure BDA0003371171010000343
And positioned as follows:
if Z is 3 =0, no voltage sensor 2 failure occurred in the multiple fault system 1;
if Z is 3 =1, the multiple fault system 1 has a voltage sensor 2 fault.
And the multi-fault diagnosis of the single-phase three-level rectifier is finished.
In order to prove the technical effect of the invention, the invention is simulated.
FIG. 4 shows an actuator failure f in this example a (t) actuator failure estimation value
Figure BDA0003371171010000344
And an adaptive threshold T th A simulated waveform diagram of (c). As can be seen from this graph, no actuator failure occurred before 8s, and the estimated value of actuator failure after 8s
Figure BDA0003371171010000345
The actuator fault f can be well estimated a (t), and actuator failure estimation
Figure BDA0003371171010000346
Has exceeded adaptive diagnostic threshold T th I.e. an actuator failure has occurred.
FIG. 5 shows the failure of the voltage sensor 1 in this example
Figure BDA0003371171010000357
Voltage sensor
1 fault estimation
Figure BDA0003371171010000351
Subject to adaptive threshold T th A simulated waveform diagram of (c). As can be seen from this graph, no failure occurred in the voltage sensor 1 before 8s, and the estimated value of the failure in the voltage sensor 1 after 8s
Figure BDA0003371171010000352
The failure f of the voltage sensor 1 can be well estimated u1 (t) and voltage sensor 1 failure estimation value
Figure BDA0003371171010000353
Has exceeded the adaptive diagnostic threshold T th Namely, the voltage sensor 1 malfunction occurs.
FIG. 6 shows the failure f of the voltage sensor 2 in this example u2 (t), voltage sensor 2 Fault estimation value
Figure BDA0003371171010000354
And an adaptive threshold T th The simulated waveform of (2). As can be seen from this graph, no failure occurred in the voltage sensor 2 before 8s, and after 8s, a failure occurred in the voltage sensor 2, and an estimated value of the failure in the voltage sensor 2
Figure BDA0003371171010000355
The voltage sensor fault 2f can be well estimated u2 (t) and voltage sensor 2 failure estimation
Figure BDA0003371171010000356
Has exceeded the adaptive diagnostic threshold T th I.e. a voltage sensor 2 failure has occurred.

Claims (3)

1. A multi-fault diagnosis method for a single-phase three-level rectifier relates to a circuit topology structure comprising a network side voltage source U s Network side equivalent inductor L s And net side equivalent resistance R s Rectifier bridge, two identical support capacitors C d1 ,C d2 A DC side load and two identical voltage sensors; support capacitor C d1 And a support capacitor C d2 A DC positive bus P and a DC negative bus Q connected in parallel with a DC side load after being connected in series 1 Between, support the capacitor C d1 And a support capacitor C d2 The contact point of (2) is marked as a direct current bus midpoint O; two identical voltage sensors are respectively marked as a voltage sensor 1 and a voltage sensor 2, and the voltage sensor 1 is connected with a supporting capacitor C d1 At both ends of which the voltage sensor 2 is connected to the supporting capacitor C d2 Two ends of (a);
the rectifier bridge is divided into two-phase bridge arms, and the two-phase bridge arms are connected with the direct-current side load in parallel; recording two-phase bridge arms as bridge arms k, wherein k is a bridge sequence, and k = a, b; in the two-phase bridge arms, each phase of bridge arm comprises 4 switching tubes with reverse connection diodes and two clamping diodes, namely a rectifier bridge comprises 8 switching tubes with reverse connection diodes and 4 clamping diodes, and the 8 switching tubes form an actuator of a single-phase three-level rectifier; marking 8 switch tubes as switch tubes V γ denotes the serial number of the switching tube, γ =1,2,3,4, and 4 clamping diodes are denoted as clamping diodes D ckρ ρ is the number of the clamping diode, ρ =1,2; in each of the two-phase arms, a switching tube V k1 Switch tube V k2 Switch tube V k3 Switch tube V k4 Are sequentially connected in series, wherein, the switch tube V k2 And a switching tube V k3 Is recorded as the rectifier bridge input point tau k K = a, b; in each of the two-phase arms, a clamping diode D ck1 The cathode of the switch tube is connected with the switch tube V k1 And a switching tube V k2 Between, the clamping diode D ck1 Anode of (2) is connected to a clamping diode D ck2 Cathode of (2), clamping diode D ck2 Anode of the switch tube is connected with the switch tube V k3 And a switching tube V k4 And clamping diode D ck1 And a clamping diode D ck2 The connecting point of the direct current bus is connected with the midpoint O of the direct current bus;
the network side equivalent inductance L s One end of is connected with an input point tau of the rectifier bridge a The other end is connected with the equivalent resistance R of the network side in sequence s Grid side voltage source U s The other end of the series network side voltage source is connected with an input point tau of the rectifier bridge b
The multi-fault diagnosis method is characterized by comprising the following steps of performing fault diagnosis on an actuator of a single-phase three-level rectifier and performing weak fault diagnosis on a multi-voltage sensor:
step 1, establishing a hybrid logic dynamic model of a single-phase three-level rectifier, and calculating a phase voltage U of an input end of a rectifier bridge ab Is estimated by
Figure FDA0003911108320000021
The network side current is sampled and recorded as the network side current i s Sampling support capacitor C d1 And a support capacitor C d2 And is denoted as DC voltage u 1 ,u 2 Sampling the DC side voltage U dc (ii) a Establishing a mixed logic dynamic model of the single-phase three-level rectifier, and calculating the phase voltage U of the input end of the rectifier bridge ab Is estimated value of
Figure FDA0003911108320000022
The input end phase voltage U of the rectifier bridge ab For the input point tau of the rectifier bridge a And the input point tau of the rectifier bridge b A voltage in between;
the expression of the hybrid logic dynamic model of the single-phase three-level rectifier is as follows:
Figure FDA0003911108320000023
Figure FDA0003911108320000024
wherein, the first and the second end of the pipe are connected with each other,
Figure FDA0003911108320000025
as an estimate of the voltage of the a-phase,
Figure FDA0003911108320000026
Is an estimate of the b-phase voltage, T h Is a DC voltage u 1 Mixed logic dynamic function of, T l Is a direct current voltage u 2 A hybrid logical dynamic function of (a);
the input end phase voltage U of the rectifier bridge ab Is estimated value of
Figure FDA0003911108320000027
The expression of (a) is:
Figure FDA0003911108320000028
step 2, establishing a multi-fault-containing single-phase three-level rectifier system state space expression, and recording the expression as a multi-fault system 1, wherein the expression is as follows:
Figure FDA0003911108320000031
wherein x (t) is a state variable of the multi-fault system 1 and is marked as a primary state variable x (t),
Figure FDA0003911108320000032
x (t) is n 1 Dimension state variable, as
Figure FDA0003911108320000033
t is a variable of time and t is,
Figure FDA0003911108320000034
the first derivative of the primary state variable x (t),
Figure FDA0003911108320000035
wherein
Figure FDA0003911108320000036
Is a net side current i s The derivative of (a) of (b),
Figure FDA0003911108320000037
are respectively DC voltage u 1 ,u 2 A derivative of (a); u (t) is the input of the multi-fault system 1 and is recorded as primary system input u (t), u (t) = u s U (t) is n 2 Dimension state variable, as
Figure FDA0003911108320000038
f a (t) actuator failure for multiple failure System 1, denoted as actuator failure f a (t),f a (t) is n 3 Dimensional state variables, noted
Figure FDA0003911108320000039
Eta (t) is the harmonic disturbance quantity of the multi-fault system 1 and is recorded as harmonic disturbance eta (t), and eta (t) is n 4 Dimension state variable, as
Figure FDA00039111083200000310
y (t) is the output of the multi-fault system 1 and is denoted as the system output y (t), and y (t) is n 5 Dimension state variable, as
Figure FDA00039111083200000311
f s (t) Voltage sensor Fault for multiple Fault System 1, noted System Voltage sensor Fault f s (t),f s (t) is n 6 Dimension state variable, as
Figure FDA00039111083200000312
Wherein
Figure FDA00039111083200000313
For a fault in voltage sensor 1, it is noted as a fault f in voltage sensor 1 u1 (t),f u1 (t) is n 7 Dimension state variable, as
Figure FDA00039111083200000314
f u2 (t) is a voltage sensor 2 failure, denoted as voltage sensor 2 failure f u2 (t),f u2 (t) is n 8 Dimension state variable, as
Figure FDA00039111083200000315
A 1 The state matrix for the multi-fault system 1 is marked as a primary state matrix A 1
Figure FDA00039111083200000316
Wherein L is equivalent inductance L at network side s R is the equivalent resistance R of the network side s Resistance value of C 1 ,C 2 Are respectively a support capacitor C d1 And a support capacitor C d2 The capacitance value of (a); b 1 The input matrix for the multiple fault system 1 is marked as a primary input matrix B 1
Figure FDA0003911108320000041
G 1 For actuator failure f a (t) the coefficient matrix is recorded as a primary actuator fault coefficient matrix,
Figure FDA0003911108320000042
N 1 the coefficient matrix of harmonic disturbance eta (t) is recorded as a primary disturbance coefficient matrix N 1
Figure FDA0003911108320000043
C y1 The output matrix of the multi-fault system 1 is recorded as a primary output matrix
Figure FDA0003911108320000044
Figure FDA0003911108320000045
F 1 For system voltage sensor fault f s (t) coefficient matrix, memoryIs a primary voltage sensor fault coefficient matrix F 1
Figure FDA0003911108320000046
System voltage sensor fault f s (t) actuator failure f a (t) and the harmonic disturbance η (t) are both continuous and bounded functions satisfying the time variable t, denoted | | f s (t)||≤r s ,||f a (t)||≤r a ,||η(t)||≤r η Wherein | | | f s (t) | | is f s Norm of (t, | f a (t) | | is f a The norm of (t), where | | | η (t) | | is the norm of η (t), r s For system voltage sensor fault f s Upper limit of (t), r a For actuator failure f a Upper limit of (t), r η Is a supremum of harmonic disturbance η (t), and r s 、r a And r η Are all normal numbers and are recorded as r s >0,r a >0,r η >0;
Step 3, introducing system state variable augmentation transformation to the multi-fault system 1, and defining the following augmentation matrixes and vectors, wherein the expression is as follows:
Figure FDA0003911108320000047
Figure FDA0003911108320000051
Figure FDA0003911108320000052
a new system of the multi-fault system 1 after system state variable augmentation transformation can be obtained, and the new system is marked as a multi-fault system 2, and the state space expression of the new system is as follows:
Figure FDA0003911108320000053
wherein z is 1 (t) is the state variable of the multiple fault system 2, noted as the secondary state variable z 1 (t),z 1 (t) is n 1 +n 3 +n 6 Dimensional space vector, is
Figure FDA0003911108320000054
Figure FDA0003911108320000055
Is a quadratic state variable z 1 (ii) the derivative of (t),
Figure FDA0003911108320000056
wherein
Figure FDA0003911108320000057
Is a failure f of the voltage sensor 1 u1 (ii) the derivative of (t),
Figure FDA0003911108320000058
for failure of the voltage sensor 2
Figure FDA0003911108320000059
A derivative of (a); f. of 1 (t) is a fault variable of the multi-fault system 2, and is recorded as a system fault quantity f 1 (t);
E 1 As derivatives of state variables of second order
Figure FDA00039111083200000510
Is denoted as the first derivative coefficient matrix E 1 In a
Figure FDA00039111083200000511
In the expression of (a) in (b),
Figure FDA00039111083200000512
represents n 1 Dimension unitThe matrix is a matrix of a plurality of pixels,
Figure FDA00039111083200000513
represents n 3 A dimension unit matrix;
A 2 the state matrix for the multi-fault system 2 is denoted as a quadratic state matrix A 2
B 2 The input matrix of the multi-fault system 2 is marked as a secondary input matrix B 2
F 2 Is the amount of system failure f 1 (t) coefficient matrix, denoted as quadratic failure coefficient matrix F 2 In a
Figure FDA0003911108320000061
In the expression of (a) in (b),
Figure FDA0003911108320000062
represents n 7 The dimension-unit matrix is a matrix of the dimension units,
Figure FDA0003911108320000063
represents n 8 A dimension unit matrix;
N 2 the coefficient matrix of harmonic disturbance eta (t) in the multi-fault system 2 is recorded as a secondary disturbance coefficient matrix N 2
Figure FDA0003911108320000064
Output matrix for multiple fault system 2
Figure FDA0003911108320000065
Is recorded as a quadratic output matrix
Figure FDA0003911108320000066
Step 4, carrying out one-time coordinate transformation on the multi-fault system 2
Step 4.1, outputting the secondary output matrix
Figure FDA0003911108320000067
Transposing, recording the transposed output matrix as
Figure FDA0003911108320000068
Handle
Figure FDA0003911108320000069
The matrix is decomposed into the product of an orthogonal matrix and an upper triangular matrix, and the product is recorded as
Figure FDA00039111083200000610
U is an orthogonal matrix and Q is an upper triangular matrix; let P = U again T ,J=Q T ,U T Is the transpose of the orthogonal matrix U, Q T Is the transpose of the upper triangular matrix Q, then the result is
Figure FDA00039111083200000611
Wherein P satisfies P T =P -1 ,P T Is the transpose of the matrix P, P -1 Is the inverse matrix of the matrix P, and takes P as a primary coordinate transformation matrix;
step 4.2, introduce coordinate transformation z 2 (t)=Pz 1 (t), the multi-fault system 2 can be converted into a new system through coordinate transformation, the new system is recorded as a multi-fault system 3, and the state space expression of the new system is as follows:
Figure FDA00039111083200000612
wherein z is 2 (t) is the state variable of the multi-fault system 3, and is recorded as a cubic state variable z 2 (t);
Figure FDA00039111083200000613
Is a cubic state variable z 2 (t) derivative of;
E 2 is a third state variable derivative
Figure FDA00039111083200000614
Is denoted as the second derivative coefficient matrix E 2 ,E 2 =E 1 P T ;A 3 The state matrix for the multi-fault system 3, denoted as the cubic state matrix A 3 ,A 3 =A 2 P T ;B 3 The input matrix of the multi-fault system 3 is recorded as a cubic input matrix B 3 ;F 3 For faults f of the multiple fault system 3 1 (t) coefficient matrix, denoted as cubic failure coefficient matrix F 3 ;N 3 The coefficient matrix of the harmonic disturbance eta (t) of the multi-fault system 3 is recorded as a third-order disturbance coefficient matrix N 3 (ii) a J is the output matrix of the multiple fault system 3, denoted as the cubic output matrix J, and is denoted as the block matrix, i.e., J = (J) 1 0), wherein J 1 Left blocking submatrix of J, J 1 Belong to n 5 Dimensional space, is denoted as
Figure FDA0003911108320000071
And J 1 Is a non-singular matrix;
step 5, matrix transformation is carried out on the multi-fault system 3
Step 5.1, the second derivative coefficient matrix E 2 Conversion into a block matrix, i.e. E 2 =(E 21 E 22 ) In which E 21 Is a coefficient matrix E of the second derivative 2 Left blocking submatrix of (E) 21 Is (n) 1 +n 3 +n 6 )×n 5 Dimensional space, is denoted as
Figure FDA0003911108320000072
E 22 Is a coefficient matrix E of the second derivative 2 Right blocking submatrix, E 22 Is (n) 1 +n 3 +n 6 )×(n 1 +n 3 +n 6 -n 5 ) Dimensional space, is denoted as
Figure FDA0003911108320000073
Then designing a transformation matrix xi, which is recorded as a matrix transformation matrix xi,
Figure FDA0003911108320000074
wherein xi 11 Being an upper blocking submatrix of the transform matrix xi, xi 11 Is n 5 ×(n 1 +n 3 +n 6 ) Dimensional space, as
Figure FDA0003911108320000075
Ξ 21 Being a lower blocking submatrix of the transform matrix xi, xi 21 Belong to (n) 1 +n 3 +n 6 -n 5 )×(n 1 +n 3 +n 6 ) Dimensional space, is denoted as
Figure FDA0003911108320000076
And xi 11 Satisfies E 22 T Ξ 11 T =0, solution of formula E 22 T Ξ 11 T =0, then may be xi 11 Wherein the matrix E 22 T Is a coefficient matrix E of the second derivative 2 Right blocking submatrix E 22 Transpose of (2), matrix xi 11 T An upper blocking submatrix xi being a transform matrix xi 11 Transposing; xi 21 =E 22 T (E 22 E 22 T ) -1 Matrix (E) 22 E 22 T ) -1 Is a matrix E 22 E 22 T The inverse of (1);
step 5.2, the left and right state space expressions of the multi-fault system 3 are multiplied by the transformation matrix xi to obtain a new transformed system, the new system is marked as a multi-fault system 4, and the state space expressions are as follows:
Figure FDA0003911108320000081
wherein z is 3 (t) is the state variable of the multiple fault system 4, noted as the quartic state variable z 3 (t);
Figure FDA0003911108320000082
Is a quartic state variable z 3 (t) derivative of;
E 3 is the fourth derivative of the state variable
Figure FDA0003911108320000083
Is recorded as a coefficient matrix of the third derivative E 3
Figure FDA0003911108320000084
Wherein E 31 Is a coefficient matrix E of the third derivative 3 Upper left blocking sub-matrix of (E) 32 Is a coefficient matrix E of the third derivative 3 The lower left block sub-matrix of (a),
Figure FDA0003911108320000085
is n 1 +n 3 +n 6 -n 5 A dimension unit matrix; a. The 4 The state matrix of the multi-fault system 4 is marked as a fourth state matrix A 4
Figure FDA0003911108320000086
Wherein A is 411 Is a four-times state matrix A 4 Upper left blocking sub-matrix of, A 412 Is a four-times state matrix A 4 Upper right blocking sub-matrix of, A 421 Is a four-times state matrix A 4 Left lower blocking submatrix of, A 422 Is a four-times state matrix A 4 The lower right blocking submatrix; b is 4 An input matrix for the multiple fault system 4, denoted as a quadruple input matrix B 4
Figure FDA0003911108320000087
Wherein B is 41 Is a four-input matrix B 4 Upper block matrix of, B 42 Is a four-input matrix B 4 A lower block matrix of (a); f 4 For faults f of the multiple fault system 4 1 (t) coefficient matrix, denoted as Quaternary failure coefficient matrix F 4
Figure FDA0003911108320000088
Wherein F 41 Is a four-fault coefficient matrix F 41 Upper block matrix of F 42 Is a four-fault coefficient matrix F 4 A lower block matrix of (a); n is a radical of 4 The coefficient matrix for the harmonic disturbance η (t) of the sensor fault system 4 is noted as the fourth-order disturbance coefficient matrix N 4
Figure FDA0003911108320000089
Wherein N is 41 Is a four-times disturbance coefficient matrix N 41 Upper blocking submatrix of (2), N 42 Is a four-times disturbance coefficient matrix N 4 A lower blocking submatrix of (a);
step 6, carrying out secondary coordinate transformation on the multi-fault system 4
Step 6.1 design coordinate transformation matrix T, record as quadratic coordinate transformation matrix T, and express T as block matrix, that is
Figure FDA0003911108320000091
Wherein
Figure FDA0003911108320000092
Is n 5 The dimension-unit matrix is a matrix of dimension units,
Figure FDA0003911108320000093
is n 1 +n 3 +n 6 -n 5 Dimension unit matrix, L is a free matrix, denoted as free matrix L, which belongs to (n) 1 +n 3 +n 6 -n 5 )×n 5 Dimensional space, is denoted as
Figure FDA0003911108320000094
Step 6.2 introduction of coordinate transformation z 4 (t)=Tz 3 (t) obtaining a transformed new system, and marking the new system as a multi-fault system 5, wherein the expression of the state space of the new system is as follows:
Figure FDA0003911108320000095
wherein z is 4 (t) is the state variable of the multiple fault system 5, noted as the quintic state variable z 4 (t) for the fifth state variable z 4 (t) blocking, i.e.
Figure FDA0003911108320000096
z 41 (t) is a quintic state variable z 4 (t) upper block subvector z 41 (t),z 41 (t) is n 5 Dimension vector, z 42 (t) is a five-state variable z 4 (t) lower block subvector z 42 (t),z 42 (t) is n 1 +n 3 +n 6 -n 5 A dimension vector;
Figure FDA0003911108320000097
is a five-fold state variable z 4 Derivative of (t), i.e.
Figure FDA0003911108320000098
Figure FDA0003911108320000099
Is z 4 (t) an upper block vector z 41 (ii) the derivative of (t),
Figure FDA00039111083200000910
is z of the upper block vector 42 (t) a derivative of; is transformed by the above coordinate z 4 (t)=Tz 3 (t) it can be known that,
Figure FDA00039111083200000911
E 4 fifth order state variable derivative
Figure FDA00039111083200000912
A coefficient matrix of (a), and
Figure FDA00039111083200000913
A 5 the state matrix for the multi-fault system 5 is denoted as the five-time state matrix A 5
Figure FDA00039111083200000914
Wherein A is 511 Is a quintic state matrix A 5 Upper left blocking sub-matrix of, A 512 Is a quintic state matrix A 5 Upper right blocking sub-matrix of, A 521 Is a quintic state matrix A 5 Left lower blocking submatrix of, A 522 Is a quintic state matrix A 5 The lower right blocking submatrix; t is a unit of -1 An inverse matrix of the quadratic coordinate transformation matrix T; b is 5 The input matrix for the multiple fault system 5 is marked as five-times input matrix B 5 ,
Figure FDA0003911108320000101
Wherein B is 51 For five inputs of matrix B 5 Upper blocking submatrix of (2), B 52 For five inputs of matrix B 5 A lower block submatrix of (a); f 5 The coefficient matrix for the fault F (t) of the multiple fault system 5 is noted as the quintic fault coefficient matrix F 5
Figure FDA0003911108320000102
Wherein F 51 Is a quintic fault coefficient matrix F 5 Upper blocking submatrix of F 52 Is a quintic fault coefficient matrix F 5 A lower block submatrix of (a); n is a radical of hydrogen 5 The coefficient matrix of the harmonic disturbance eta (t) of the multi-fault system 5 is recorded as a quintic disturbance coefficient matrix N 5
Figure FDA0003911108320000103
Wherein N is 51 For a quintic disturbance coefficient matrix N 5 Upper block submatrix of, N 52 For a quintic disturbance coefficient matrix N 5 Lower blocking submatrix, J * Is the output matrix of the multiple fault system 5 and knows the following relationship:
Figure FDA0003911108320000104
step 7, designing a self-adaptive reduced-order sliding mode observer
Step 7.1, calculating the state quantity theta of the order-reduced system according to the multi-fault system 5, and recording the state quantity theta as the state quantity theta of the order-reduced system, wherein the calculation formula is as follows:
Θ=(LE 31 +E 32 -L)z 41 (t)+z 42 (t)
and solving a state space expression of the reduced-order system according to the state quantity theta of the reduced-order system, wherein the state space expression is as follows:
Figure FDA0003911108320000105
wherein the content of the first and second substances,
Figure FDA0003911108320000106
is the derivative of the state quantity theta of the reduced-order system, and delta is the state expression z of the reduced-order system 41 The coefficient matrix of (t) is marked as a reduced system coefficient matrix delta, and the expression is as follows:
Δ=-(LA 512 +A 522 )(LE 51 +E 52 -L)+L(A 511 -A 512 L)+A 511 -A 522 L
step 7.2, designing a self-adaptive reduced order sliding mode observer aiming at a state space expression of a reduced order system, wherein the expression is as follows:
Figure FDA0003911108320000111
wherein the content of the first and second substances,
Figure FDA0003911108320000112
is the estimated value of the state quantity theta of the reduced order system and is recorded as the estimated value of the state quantity of the reduced order system
Figure FDA0003911108320000113
Figure FDA0003911108320000114
Is an estimate of the state quantity of the reduced order system
Figure FDA0003911108320000115
Is a sliding mode gain matrix, Γ = (LF) 51 +F 52 ) And v is the approach law,
Figure FDA0003911108320000116
wherein the content of the first and second substances,
Figure FDA0003911108320000117
in order to be able to vary the parameter 1,
Figure FDA0003911108320000118
tanh () is a hyperbolic tangent function,
Figure FDA0003911108320000119
is variable parameter 2, and
Figure FDA00039111083200001110
sigma is a variable parameter 3, sigma belongs to (0, 1), theta is a variable parameter 4, theta is more than 1, rho is a variable parameter 5, beta is more than 0, psi is a positive definite symmetrical matrix, and chi is a diagonal matrix,
Figure FDA00039111083200001111
epsilon is a constant number term, epsilon > 1 Θ (t) is the error of the estimation,
Figure FDA00039111083200001112
let S = e Θ (t), wherein S is a sliding mode surface of the designed adaptive reduced-order sliding mode observer;
step 7.3, obtaining a free matrix L by solving a Lyapunov equation, wherein the expression of the Lyapunov equation is as follows:
(LA 412 +A 422 ) T P+P(LA 412 +A 422 )=-I
wherein (LA) 412 +A 422 ) T Is LA 412 +A 422 P is a positive definite symmetric matrix, and I is a unit matrix;
step 8, fault diagnosis of the actuator and the sensor is carried out
Step 8.1, sampling the value of the system output y (t) of the multiple fault system 1 in step 1 and substituting the sampled value into z in step 6.2 41 (t)=z 31 (t)=J 1 -1 y (t), the state variable z can be obtained five times 4 (t) upper block subvector z 41 (t) value, again according to the known five state variable z 4 (t) an upper block vector z 41 Value of (t), Θ = (LE) in step 7.1 31 +E 32 -L)z 41 (t)+z 42 (t), step 7.2
Figure FDA00039111083200001113
And the error e is estimated in step 7.2 Θ (t) =0, and z can be obtained 42 (t) the expression is:
Figure FDA0003911108320000121
then the five state variables z are obtained 4 (t) an upper block subvector z 41 (t) value and quintic state variable z 4 (t) lower block subvector z 42 (t) value is substituted into step 6.2
Figure FDA0003911108320000122
The state variable z can be calculated five times 4 (t) a value;
step 8.2, starting from the multiple fault system 2, according to a coordinate transformation z 2 (t)=Pz 1 (t) and quadratic coordinate transformation z 4 (t)=Tz 3 (t) and the matrix transformation does not change the cubic state variable z in the sensor failure system 3 in step 5 2 (t) obtaining z 3 (t)=z 2 (t) to obtain z 4 (t)=TPz 1 (t) calculating a secondary state variable z by inverse operation of the matrix 1 (t)=P T T - 1 z 4 (t), and the five state variables z calculated in step 8.1 are further added 4 (t) substitution of formula z 1 (t)=P T T -1 z 4 (t) obtaining a secondary state variable z 1 (t);
Step 8.3, mixing
Figure FDA0003911108320000123
Is recorded as a primary state variable estimation value,
Figure FDA0003911108320000124
recording as the estimated value of the actuator fault
Figure FDA0003911108320000125
Figure FDA0003911108320000126
Is recorded as the fault estimation value of the voltage sensor 1
Figure FDA0003911108320000127
Figure FDA0003911108320000128
Is recorded as the fault estimation value of the voltage sensor 2
Figure FDA0003911108320000129
Calculating a state variable estimate
Figure FDA00039111083200001210
Actuator fault estimation
Figure FDA00039111083200001211
Voltage sensor fault 1 estimation
Figure FDA00039111083200001212
And voltage sensor 2 fault estimation
Figure FDA00039111083200001213
The specific calculation formula is as follows:
Figure FDA00039111083200001214
Figure FDA00039111083200001215
Figure FDA00039111083200001216
Figure FDA00039111083200001217
step 8.4, diagnosing faults of the actuator, the voltage sensor 1 and the voltage sensor 2, and giving a self-adaptive diagnosis threshold value T th
Defining fault location quantity Z 1
Figure FDA00039111083200001218
And positioned as follows:
if Z is 1 =0, no actuator failure occurred in the multiple fault system 1;
if Z is 1 =1, the multi-fault system 1 has an actuator fault;
defining fault location quantity Z 2
Figure FDA0003911108320000131
And positioned as follows:
if Z is 2 =0, no voltage sensor 1 fault occurred in the multiple fault system 1;
if Z is 2 =1, multiple fault system 1 has voltage sensor 1 fault;
defining fault location quantity Z 3
Figure FDA0003911108320000132
And positioned as follows:
if Z is 3 =0, no voltage sensor 2 failure occurred in the multiple fault system 1;
if Z is 3 =1, the multiple fault system 1 has a voltage sensor 2 fault.
2. The method according to claim 1, wherein the DC voltage u of step 1 is the DC voltage u 1 Hybrid logic dynamic function T of h And a DC voltage u 2 Mixed logic dynamic function T of l The calculation process of (2) is as follows:
recording the switching function of the k-phase bridge arm as S k K = a, b, then:
Figure FDA0003911108320000133
Figure FDA0003911108320000134
recording the pulse control signal of the switching tube as v K = a, b, γ =1,2,3,4, then the k-phase bridge arm switching function S k =v k1 v k2 -v k3 v k4
3. Method for multiple fault diagnosis of a single-phase three-level rectifier according to claim 1, characterized in that the diagnosis adaptive threshold T of step 8.4 is defined th The expressions are respectively as follows:
T th =E 1112 )
wherein E is 1 Is a constant 1, Γ 1 Is a constant number 2, and Γ 1 ∈(1,2),ζ 1 For a bounded external disturbance, ζ 2 Is a bounded voltage perturbation.
CN202111417465.3A 2021-11-24 2021-11-24 Multi-fault diagnosis method for single-phase three-level rectifier Active CN114325164B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111417465.3A CN114325164B (en) 2021-11-24 2021-11-24 Multi-fault diagnosis method for single-phase three-level rectifier

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111417465.3A CN114325164B (en) 2021-11-24 2021-11-24 Multi-fault diagnosis method for single-phase three-level rectifier

Publications (2)

Publication Number Publication Date
CN114325164A CN114325164A (en) 2022-04-12
CN114325164B true CN114325164B (en) 2023-03-10

Family

ID=81045740

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111417465.3A Active CN114325164B (en) 2021-11-24 2021-11-24 Multi-fault diagnosis method for single-phase three-level rectifier

Country Status (1)

Country Link
CN (1) CN114325164B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114825281B (en) * 2022-04-22 2024-03-26 合肥工业大学 Multi-fault estimation method of staggered parallel Boost PFC system
CN114914888B (en) * 2022-04-25 2024-03-26 合肥工业大学 Multi-fault estimation method for single-phase three-level rectifier

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0053229A1 (en) * 1980-12-01 1982-06-09 VEB Starkstrom-Anlagenbau Leipzig-Halle Method and circuit for detecting and breaking disturbing electric arcs
KR20020016182A (en) * 2000-08-24 2002-03-04 권영한 Apparatus and method for diagnosing error of phase control rectifier
DE10112569A1 (en) * 2001-03-15 2002-10-02 Bosch Gmbh Robert Fault diagnosis method for three phase stationary or motor vehicle generators for early detection of broken connections between generator and rectifier to prevent consequential damage
JP2005114440A (en) * 2003-10-06 2005-04-28 Japan Aviation Electronics Industry Ltd Acceleration sensor of capacitance detection type capable of diagnosing malfunction
CN203933382U (en) * 2014-05-28 2014-11-05 合肥华耀电子工业有限公司 A kind of APFC of tape jam measuring ability
WO2019091914A1 (en) * 2017-11-10 2019-05-16 Renault S.A.S Method for controlling a three-phase vienna rectifier when a power switch is faulty
CN110609194A (en) * 2019-06-29 2019-12-24 南京理工大学 Three-phase rectifier open-circuit fault diagnosis method based on voltage space vector
CN111751760A (en) * 2020-06-12 2020-10-09 武汉大学 Three-phase rectifier power tube fault diagnosis method and device based on current signals
CN111983508A (en) * 2020-07-09 2020-11-24 华中科技大学 T-type three-phase four-wire rectifier fault real-time detection and positioning method and system
CN113031570A (en) * 2021-03-18 2021-06-25 哈尔滨工业大学 Rapid fault estimation method and device based on self-adaptive unknown input observer
CN113281680A (en) * 2021-05-21 2021-08-20 合肥工业大学 Open-circuit fault diagnosis method for single-phase three-level rectifier of high-speed rail traction system

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107832687A (en) * 2017-10-27 2018-03-23 武汉大学 Fault diagnostic method for transformer winding based on wireless identification sensing
CN108416103A (en) * 2018-02-05 2018-08-17 武汉大学 A kind of method for diagnosing faults of electric automobile of series hybrid powder AC/DC convertor
US10746803B2 (en) * 2018-12-06 2020-08-18 Hamilton Sunstrand Corporation Fault detection and isolation in generator modules

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0053229A1 (en) * 1980-12-01 1982-06-09 VEB Starkstrom-Anlagenbau Leipzig-Halle Method and circuit for detecting and breaking disturbing electric arcs
KR20020016182A (en) * 2000-08-24 2002-03-04 권영한 Apparatus and method for diagnosing error of phase control rectifier
DE10112569A1 (en) * 2001-03-15 2002-10-02 Bosch Gmbh Robert Fault diagnosis method for three phase stationary or motor vehicle generators for early detection of broken connections between generator and rectifier to prevent consequential damage
JP2005114440A (en) * 2003-10-06 2005-04-28 Japan Aviation Electronics Industry Ltd Acceleration sensor of capacitance detection type capable of diagnosing malfunction
CN203933382U (en) * 2014-05-28 2014-11-05 合肥华耀电子工业有限公司 A kind of APFC of tape jam measuring ability
WO2019091914A1 (en) * 2017-11-10 2019-05-16 Renault S.A.S Method for controlling a three-phase vienna rectifier when a power switch is faulty
CN110609194A (en) * 2019-06-29 2019-12-24 南京理工大学 Three-phase rectifier open-circuit fault diagnosis method based on voltage space vector
CN111751760A (en) * 2020-06-12 2020-10-09 武汉大学 Three-phase rectifier power tube fault diagnosis method and device based on current signals
CN111983508A (en) * 2020-07-09 2020-11-24 华中科技大学 T-type three-phase four-wire rectifier fault real-time detection and positioning method and system
CN113031570A (en) * 2021-03-18 2021-06-25 哈尔滨工业大学 Rapid fault estimation method and device based on self-adaptive unknown input observer
CN113281680A (en) * 2021-05-21 2021-08-20 合肥工业大学 Open-circuit fault diagnosis method for single-phase three-level rectifier of high-speed rail traction system

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
An experimental research on comparison of two kinds of voltage sag generators;Wu Jian 等;《Proceedings of The 7th International Power Electronics and Motion Control Conference》;20120806;全文 *
基于极端树与堆栈式稀疏自编码算法的电力电子电路故障诊断;张亚茹 等;《电子测量技术》;20191123;第42卷(第22期);全文 *

Also Published As

Publication number Publication date
CN114325164A (en) 2022-04-12

Similar Documents

Publication Publication Date Title
CN114325164B (en) Multi-fault diagnosis method for single-phase three-level rectifier
Xiao et al. Event-triggered H∞ filtering of discrete-time switched linear systems
CN110346736B (en) NPC three-level inverter fault diagnosis method based on improved treelet transformation
CN113281680B (en) Open-circuit fault diagnosis method for single-phase three-level rectifier of high-speed rail traction system
CN112101457B (en) PMSM demagnetizing fault diagnosis method based on torque signal fuzzy intelligent learning
CN110441643A (en) Inverter power pipe open circuit fault diagnostic method in control system for permanent-magnet synchronous motor
CN113534000B (en) New energy automobile driving system inverter and current sensor fault diagnosis method
CN113285481B (en) Grid-connected converter inductance parameter online estimation method, prediction control method and system
CN110456217B (en) MMC fault positioning method based on WPD-FOA-LSSVM dual model
CN113534035B (en) Tiny fault diagnosis method for current sensor of new energy electric vehicle driving system
CN112310991B (en) Energy storage converter improved finite set model prediction control method based on disturbance observation
Guo et al. Fault diagnosis of modular multilevel converter based on adaptive chirp mode decomposition and temporal convolutional network
CN111983414B (en) Open-circuit fault diagnosis method and system for rail train traction converter
CN114280526B (en) Digital differential traceability system and method for electronic transformer calibrator
CN112731193B (en) NPC inverter multi-class fault diagnosis method based on sliding-mode observer
CN111413646A (en) Real-time detection method for open-circuit fault of bridge arm of three-phase voltage source inverter
Shi et al. Digital twin approach for IGBT parameters identification of a three-phase DC-AC inverter
CN116702060A (en) Multi-level inverter power device fault diagnosis method
CN112948755B (en) Telemetry sine parameter interpretation method
CN114914888B (en) Multi-fault estimation method for single-phase three-level rectifier
Torki et al. Sparse recovery diagnosis method applied to hybrid dynamical system: The case of three-phase dc-ac inverter for wind turbine
CN112198458B (en) Method and system for detecting open-circuit fault of three-phase voltage source inverter in real time
Suskis et al. Converter state-space model estimation using dynamic mode decomposition
CN110297426B (en) Parallel type half-cycle repetitive control method adopting periodic interference differential suppression strategy
CN114325380B (en) Fault diagnosis method for permanent magnet synchronous motor driving system

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant