CN106527144B - 磁悬浮系统的抗干扰轨迹跟踪降阶控制方法 - Google Patents
磁悬浮系统的抗干扰轨迹跟踪降阶控制方法 Download PDFInfo
- Publication number
- CN106527144B CN106527144B CN201611145515.6A CN201611145515A CN106527144B CN 106527144 B CN106527144 B CN 106527144B CN 201611145515 A CN201611145515 A CN 201611145515A CN 106527144 B CN106527144 B CN 106527144B
- Authority
- CN
- China
- Prior art keywords
- matrix
- magnetic suspension
- suspension system
- controlled device
- follows
- 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
Links
- 239000000725 suspension Substances 0.000 title claims abstract description 82
- 238000000034 method Methods 0.000 title claims abstract description 50
- 239000011159 matrix material Substances 0.000 claims abstract description 119
- 229910000831 Steel Inorganic materials 0.000 claims description 18
- 239000010959 steel Substances 0.000 claims description 18
- 238000006073 displacement reaction Methods 0.000 claims description 14
- 238000010276 construction Methods 0.000 claims description 10
- 239000000203 mixture Substances 0.000 claims description 10
- 238000000354 decomposition reaction Methods 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 230000017105 transposition Effects 0.000 claims description 6
- KUEUWHJGRZKESU-UHFFFAOYSA-N Niceritrol Chemical compound C=1C=CN=CC=1C(=O)OCC(COC(=O)C=1C=NC=CC=1)(COC(=O)C=1C=NC=CC=1)COC(=O)C1=CC=CN=C1 KUEUWHJGRZKESU-UHFFFAOYSA-N 0.000 claims description 4
- 230000005611 electricity Effects 0.000 claims description 4
- 230000005484 gravity Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical group [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 claims description 3
- 230000001133 acceleration Effects 0.000 claims description 3
- 230000003252 repetitive effect Effects 0.000 claims description 3
- 230000009131 signaling function Effects 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000002401 inhibitory effect Effects 0.000 abstract 1
- 230000008859 change Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000007812 deficiency Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000007667 floating Methods 0.000 description 2
- 238000005339 levitation Methods 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000005265 energy consumption Methods 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 230000002035 prolonged effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Feedback Control In General (AREA)
Abstract
磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，包括以下步骤：获取磁悬浮系统中被控对象的运动微分方程组，并通过非线性滤波器对其进行线性化处理，得到被控对象状态空间方程；根据被控对象与给定跟踪信号组成的扩阶系统建立线性矩阵方程组，并通过其解来构造一个基于输出反馈的结构化控制器，实现对给定信号的跟踪；在未知扰动和噪声的影响下，通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷输出反馈控制器改善系统抑制未知扰动和噪声的能力。
Description
技术领域
本发明涉及磁悬浮技术领域，具体涉及一种磁悬浮系统的抗干扰轨迹跟踪降阶控制方法。
背景技术
磁悬浮是通过电磁力作用于物体使其克服自身重力保持悬浮的一种新型技术。由于磁悬浮技术具有无接触、功耗低、输出大、污染少等特点，该技术能有效的延长机械设备使用寿命，降低能量损耗，还能应用于真空，高温等特殊环境中。因此与常规技术相比，磁悬浮技术在实际工程实践中得到了更广泛的关注和应用，如磁悬浮列车、高速磁悬浮电机、磁悬浮轴承等。
在对磁悬浮系统实施控制时，由于其系统的不稳定和非线性增加了对其控制的难度。通常在传统工程中采用PID等经典控制算法，该类方法参数调节简单并且易于实现，但磁悬浮系统具有参数不确定性，该类算法很难在高精度和高可靠性的应用工程中达到实际要求。为弥补此类经典控制算法的不足，另一类复杂算法，如智能控制和非线性控制等被尝试应用于磁悬浮领域。此类算法可以很好的解决磁悬浮系统的非线性和模型误差问题，但由于此类算法自身的复杂性，因此很难在实际的磁悬浮系统控制中得到应用。此外，在磁悬浮系统的实际控制中，通常存在来自工作环境和系统内部未知扰动和噪声，而鲁棒控制算法可以有效地抑制扰动和噪声对系统的影响，但该类算法设计的控制器通常阶数较高，所以构造比较复杂且计算量较大。
发明内容
本发明所要解决上述现有技术的不足，提供一种磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，该方法可以通过一个低阶控制器改善磁悬浮系的统对未知扰动和噪声的抑制能力，易于实现，降低了成本，使被控对象的轨迹按给定信号变化，实现对信号的跟踪。
本发明为解决上述技术问题提供了如下解决方案：本发明设计了一种磁悬浮系统的抗干扰轨迹跟踪降阶控制方法。其中，磁悬浮系统包括控制器部分和被控对象；被控对象包括激光位移传感器、功率放大器、电磁铁及小钢球。控制器部分包括降阶H无穷控制器和非线性滤波器，其工作原理如下：先通过求解线性矩阵方程组构造一个结构化控制器，实现对给定信号的跟踪；再通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；同时利用非线性滤波器消除电磁力和模型误差对被控对象的影响。
磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，具体包括：
1)、在磁悬浮系统中获取所述被控对象的运动微分方程组，对其进行线性化处理，并通过非线性滤波器得到磁悬浮系统在消除电磁力和模型误差后的被控对象状态空间方程；
2)、在磁悬浮系统中得到被控对象与给定跟踪信号所组成的扩阶系统，并通过求解线性矩阵方程组构造一个结构化控制器，实现对给定信号的跟踪；
3)、在能量有界的未知扰动和噪声的影响下，通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；
进一步，步骤1)中，获取被控对象的状态空间方程的方法包括以下步骤：
(1.1)所述被控对象中磁悬浮系统的运动微分方程组如下：
其中，s为钢球的位移，i为电磁铁的控制电流，s0为钢球处于平衡状态时的位移，A为电磁铁中铁芯的导磁截面积，N为电磁铁的线圈匝数，R为电磁铁的线圈电阻，U为电磁铁的电压，i0为钢球处于平衡状态时电磁铁的控制电流，m为钢球的质量，g为重力加速度，μ0为空气磁导率，F(i,s)为非线性电磁力，L为电磁铁的静态电感；
(1.2)由于磁悬浮系统有一定的可控范围，所以可在磁悬浮系统的平衡点s0附近对其进行线性化处理；通过式(1)，可得被控对象中磁悬浮系统的运动微分方程如下：
其中，
接着，通过式(2)，可以得到被控对象的状态空间方程如下：
其中，即小钢球的位移和速度，u为控制增益，f为外部扰动和模型误差，y为观测输出，Ag为动态矩阵，Bu为输入矩阵，Bf为外部扰动和模型误差的常数矩阵，Cg为观测输出矩阵；
(1.3)使用非线性滤波器消除电磁力和模型误差对被控对象的影响；具体步骤如下：
(1.3.1)磁悬浮系统中的非线性滤波器如下：
其中，z1为钢球位移的估计量，z2为钢球速度的估计量，z3为电磁力和模型误差f的估计量，e为输出误差，ω为系统控制带宽，bu为已知系统参数，y为系统输出，u0为控制增益，Kd为干扰补偿增益；
(1.3.2)通过非线性滤波器对被控对象电磁力和模型误差f的估计量z3，可由式(4)，将磁悬浮系统中电磁力和模型误差Bff消除，同时系统在能量有界的未知扰动和噪声d的影响下的状态空间方程如下所示：
其中，z为外部扰动的输出，Cz为输出矩阵，Dzu和Dzd为常数矩阵；
进一步，步骤2)中，在磁悬浮系统中得到被控对象与给定跟踪信号所组成的扩阶系统，并通过求解线性矩阵方程组构造一个结构化控制器，实现对给定信号的跟踪；具体包括：
(2.1)在磁悬浮系统中，获取被控对象与给定跟踪信号组成的闭环状态空间方程；具体包括：
(2.1.1)引入给定的跟踪信号，其状态空间方程如下：
其中，xω为给定信号的状态量，ξ为给定信号的输入，ωi为给定信号的输出，Ai为动态矩阵，Bi为输入矩阵，Ci为输出矩阵，Di为常数矩阵。
(2.1.2)将跟踪信号作用于被控对象，由于d为能量有界的外部扰动，不会影响对系统对跟定信号的跟踪，为方便叙述，在这部分可不做考虑，所以由式(5)、(6)，使被控对象的状态空间方程变为如下：
其中，e0为系统的控制输出，即跟踪误差，Ce为跟踪误差的输出矩阵，Dω和Du为跟踪误差的常数矩阵，Dgω为观测输出的常数矩阵；
(2.1.3)通过式(6)、(7)，可以得到被控对象与跟踪信号的闭环系统如下：
其中，A11＝Ag，B12＝Bu，C11＝Ce，C12＝DωCi，C21＝Cg，C22＝DgωCi，D11＝DωDi，D12＝Du，D21＝DgωDi；
(2.2)通过求解线性矩阵方程组得到一个结构化控制器，获得被控对象与给定跟踪信号组成的扩阶系统；具体包括：
(2.2.1)在磁悬浮系统中，构造一个结构化控制器：
其中，为动态矩阵，为输入矩阵，为输出矩阵，Ak，Bk，均为常数矩阵，且Πi和Xi为线性矩阵方程组的解；
(2.2.2)通过上述步骤可得磁悬浮系统中被控对象与给定跟踪信号组成的扩阶系统如下：
其中，为此扩阶系统的状态量，eG为此扩阶系统的控制输出，yG为此扩阶系统的观测输出；
更进一步，步骤3)中，在能量有界的未知扰动和噪声的影响下，通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；具体包括：
(3.1)在能量有界的未知扰动和噪声d的影响下，通过式(10)，可得到d与磁悬浮系统所组成的闭环系统状态空间方程如下：
其中，z为外部扰动的输出，为观测输出，为动态矩阵，为输入矩阵，为观测输出矩阵，为输入矩阵，为输出矩阵，和为常数矩阵；
(3.2)根据磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，通过摄动法修正标准型线性矩阵不等式组并对其解值进行迭代；具体包括：
(3.2.1)定义：
其中，T表示矩阵的转置，⊥表示矩阵的核空间，I为单位矩阵，γ为给定的性能指标常数，P和Q为待定未知矩阵；
常数εB＞0，εr＞0，满足：且||N1+N2||表示矩阵的欧式范数，通过式(12)～(14)，可得到线性矩阵不等式组如下：
(3.2.2)对线性矩阵不等式组的解进行迭代，具体步骤如下：
(3.2.2.a)定义矩阵：M＝-N2+N1，对M进行奇异值分解，可得：
M＝ZΣV* (16)
其中，*表示矩阵的共轭转置，Σ为奇异值矩阵，Z和V为酉矩阵；
(3.2.2.b)进一步，定义：
其中，P0＝P，+表示矩阵的Moore-Penrose逆，W11，W12和W22为常数矩阵；
(3.2.2.c)将P1-Q-1进行特征值分解，可得到如下：
P1-Q-1＝Θdiag(λ1,...λl,λl+1,...λn)ΘT (19)
其中，Θ为酉矩阵，diag()表示对角矩阵，λ1，λl，λl+1，λn为矩阵P1-Q-1的特征值并且为降序排列；
(3.2.2.d)更近一步，将满足εr≥λl+1≥...≥λn的特征值置为0，可通过式(19)，得到并得到如下条件：
其中，rank()表示矩阵的秩，若所求的能使式(20)的条件成立，则满足方法要求；否则，令P0＝P1，重复式(17)～(19)，直到结果满足式(20)的条件，求得和Q；
(3.3)通过上述所求和Q，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；具体步骤如下：
(3.3.1)定义：
其中， 和为常数矩阵且满足
通过式(21)～(22)，可以得到一个线性矩阵不等式如下：
其中，为控制器中的参数矩阵，且
所以，可通过上述步骤获得控制器参数，并最后实现了磁悬浮系统的抗干扰轨迹跟踪降阶控制方法。
本发明设计的一种磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，该方法可以通过一个低阶控制器改善磁悬浮系的统对未知扰动和噪声的抑制能力，易于实现，降低了成本，使被控对象的轨迹按给定信号变化，实现对信号的跟踪。
本发明的优点是：针对磁悬浮系统自身的开环不稳定性与系统的非线性，通过一个低阶控制器改善磁悬浮系的统对未知扰动和噪声的抑制能力，同时使磁悬浮系统具备了一定的信号跟踪能力，具有一定的现实意义和理论价值，拓宽了磁悬浮系统的应用领域。
附图说明
图1是本发明方法流程图。
图2是本发明方法的实时实验效果图。
具体实施方式
为使本发明的目的、技术方案和具体效果更加清晰，下面结合附图和实际实验数据对本发明的技术方案作进一步描述。
如图1所示，本发明设计了一种磁悬浮系统的抗干扰轨迹跟踪降阶控制方法。其中，磁悬浮系统包括控制器部分和被控对象；被控对象包括激光位移传感器、功率放大器、电磁铁及小钢球。控制器部分包括降阶H无穷控制器和非线性滤波器，其工作原理如下：先通过求解线性矩阵方程组构造一个结构化控制器，实现对给定信号的跟踪；再通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；同时利用非线性滤波器消除电磁力和模型误差对被控对象的影响。
磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，具体包括：
1)、在磁悬浮系统中获取所述被控对象的运动微分方程组，对其进行线性化处理，并通过非线性滤波器得到磁悬浮系统在消除电磁力和模型误差后的被控对象状态空间方程；
2)、在磁悬浮系统中得到被控对象与给定跟踪信号所组成的扩阶系统，并通过求解线性矩阵方程组构造一个结构化控制器，实现对给定信号的跟踪；
3)、在能量有界的未知扰动和噪声的影响下，通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；
进一步，步骤1)中，获取被控对象的状态空间方程的方法包括以下步骤：
(1.1)所述被控对象中磁悬浮系统的运动微分方程组如下：
其中，s为钢球的位移，i为电磁铁的控制电流，s0为钢球处于平衡状态时的位移，电磁铁中铁芯的导磁截面积A＝0.00159m2，电磁铁的线圈匝数N＝2450匝，电磁铁的线圈电阻R＝13.8Ω，U为电磁铁的电压，钢球处于平衡状态时电磁铁的控制电流i0＝0.3943A，钢球的质量m＝94g，重力加速度g＝9.8N/kg，空气磁导率μ0＝4π×10-7H/m，F(i,x)为非线性电磁力，电磁铁的静态电感L＝135mH；
(1.2)由于磁悬浮系统有一定的可控范围，所以可在磁悬浮系统的平衡点s0附近对其进行线性化处理；通过式(1)，可得被控对象中磁悬浮系统的运动微分方程如下：
其中，
接着，通过式(2)，可以得到被控对象的状态空间方程如下：
其中，即小钢球的位移和速度，u为控制增益，f为外部扰动和模型误差，y为观测输出，动态矩阵为输入矩阵为外部扰动和模型误差的常数矩阵为Bf＝0.124，输出矩阵为
(1.3)使用非线性滤波器消除电磁力和模型误差对被控对象的影响；具体步骤如下：
(1.3.1)磁悬浮系统中的非线性滤波器如下：
其中，z1为钢球位移的估计量，z2为钢球速度的估计量，z3为电磁力和模型误差f的估计量，e为输出误差，y为系统输出，u0为控制增益，系统控制带宽ω＝410Hz，已知系统参数bu＝7.6367，干扰补偿增益为Kd＝-0.0162；
(1.3.2)通过非线性滤波器对被控对象电磁力和模型误差f的估计量z3，可由式(4)，将磁悬浮系统中电磁力和模型误差Bff消除，同时系统在能量有界的未知扰动和噪声d的影响下的状态空间方程如下所示：
其中，z为外部扰动的输出，Cz为输出矩阵，Dzu和Dzd为常数矩阵；
进一步，步骤2)中，在磁悬浮系统中得到被控对象与给定跟踪信号所组成的扩阶系统，并通过求解线性矩阵方程组构造一个结构化控制器，实现对给定信号的跟踪；具体包括：
(2.1)在磁悬浮系统中，获取被控对象与给定跟踪信号组成的闭环状态空间方程；具体包括：
(2.1.1)引入给定的跟踪信号，其状态空间方程如下：
其中，xω为给定信号的状态量，ξ为给定信号的输入，ωi为给定信号的输出，动态矩阵为输入矩阵为输出矩阵为常数矩阵为
(2.1.2)将跟踪信号作用于被控对象，由于d为能量有界的外部扰动，不会影响对系统对跟定信号的跟踪，为方便叙述，在这部分可不做考虑，所以由式(5)、(6)，使被控对象的状态空间方程变为如下：
其中，e0为系统的控制输出，即跟踪误差，跟踪误差的输出矩阵为Ce＝[1 0]，跟踪误差的常数矩阵为Dω＝[-1 0]和Du＝0，观测输出的常数矩阵为
(2.1.3)通过式(6)、(7)，可以得到被控对象与跟踪信号的闭环系统如下：
其中，A11＝Ag，B12＝Bu，C11＝Ce，C12＝DωCi，C21＝Cg，C22＝DgωCi，D11＝DωDi，D12＝Du，D21＝DgωDi；
(2.2)通过求解线性矩阵方程组得到一个结构化控制器，获得被控对象与给定跟踪信号组成的扩阶系统；具体包括：
(2.2.1)在磁悬浮系统中，构造一个结构化控制器：
其中，为动态矩阵，为输入矩阵，为输出矩阵，Ak，Bk，均为常数矩阵，且Πi＝[-256.7863 0]和为线性矩阵方程组的解；
(2.2.2)通过上述步骤可得磁悬浮系统中被控对象与给定跟踪信号组成的扩阶系统如下：
其中，为此扩阶系统的状态量，eG为此扩阶系统的控制输出，yG为此扩阶系统的观测输出；
更进一步，步骤3)中，在能量有界的未知扰动和噪声的影响下，通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；具体包括：
(3.1)在能量有界的未知扰动和噪声d的影响下，通过式(10)，可得到d与磁悬浮系统所组成的闭环系统状态空间方程如下：
其中，z为外部扰动的输出，为观测输出，动态矩阵为输入矩阵为观测输出矩阵为外部扰动的输入矩阵为外部扰动的输出矩阵为外部扰动输出的常数矩阵为和观测输出的常数矩阵为
(3.2)根据磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，通过摄动法修正标准型线性矩阵不等式组并对其解值进行迭代；具体包括：
(3.2.1)定义：
其中，T表示矩阵的转置，⊥表示矩阵的核空间，I为单位矩阵，给定的性能指标常数为γ＝0.8688，P和Q为待定未知矩阵；
常数εB＝0.1，εr＝0.01，满足：且||N1+N2||表示矩阵的欧式范数，通过式(12)～(14)，可得到线性矩阵不等式组如下：
(3.2.2)对线性矩阵不等式组的解进行迭代，具体步骤如下：
(3.2.2.a)定义矩阵：
对M进行奇异值分解，可得：
M＝ZΣV* (16)
其中，*表示矩阵的共轭转置，奇异值矩阵为酉矩阵为
(3.2.2.b)进一步，定义：
其中，P0＝P，+表示矩阵的Moore-Penrose逆，常数矩阵：
(3.2.2.c)将P1-Q-1进行特征值分解，可得到如下：
P1-Q-1＝Θdiag(λ1,...λl,λl+1,...λn)ΘT (19)
其中，为酉矩阵，diag()表示对角矩阵，λ1＝7.9179，λ2＝0.9113，λ3＝1.658×10-15，λ4＝1.8729×10-15为矩阵P1-Q-1的特征值并且为降序排列；
(3.2.2.d)更近一步，将满足εr≥λl+1≥...≥λn的特征值置为0，可通过式(19)，得到并得到如下条件：
其中，rank()表示矩阵的秩，若所求的能使式(20)的条件成立，则满足方法要求；否则，令P0＝P1，重复式(17)～(19)，直到结果满足式(20)的条件且求得
(3.3)通过上述所求和Q，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；具体步骤如下：
(3.3.1)定义：
其中， 常数矩阵为： 且满足
通过式(21)～(22)，可以得到一个线性矩阵不等式如下：
其中，为控制器中的参数矩阵，且
所以，如图2所示，可通过上述步骤获得控制器参数，并最后实现了磁悬浮系统的抗干扰轨迹跟踪降阶控制方法。
本发明设计的一种磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，该方法可以通过一个低阶控制器改善磁悬浮系的统对未知扰动和噪声的抑制能力，易于实现，降低了成本，使被控对象的轨迹按给定信号变化，实现对信号的跟踪。
本发明的优点是：针对磁悬浮系统自身的开环不稳定性与系统的非线性，通过一个低阶控制器改善磁悬浮系的统对未知扰动和噪声的抑制能力，同时使磁悬浮系统具备了一定的信号跟踪能力，具有一定的现实意义和理论价值，拓宽了磁悬浮系统的应用领域。
以上结合附图详细说明和陈述了本发明的实施方式，但并不局限于上述方式。在本领域的技术人员所具备的知识范围内，只要以本发明的构思为基础，还可以做出多种变化和改进。
Claims (1)
1.磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，具体包括以下步骤：
1)、在磁悬浮系统中获取被控对象的运动微分方程组，对其进行线性化处理，并通过非线性滤波器得到磁悬浮系统在消除电磁力和模型误差后的被控对象状态空间方程；具体包括：
(1.1)所述被控对象中磁悬浮系统的运动微分方程组如下：
其中，s为钢球的位移，i为电磁铁的控制电流，s0为钢球处于平衡状态时的位移，A为电磁铁中铁芯的导磁截面积，N为电磁铁的线圈匝数，R为电磁铁的线圈电阻，U为电磁铁的电压，i0为钢球处于平衡状态时电磁铁的控制电流，m为钢球的质量，g为重力加速度，μ0为空气磁导率，F(i,s)为非线性电磁力，L为电磁铁的静态电感；
(1.2)由于磁悬浮系统有一定的可控范围，所以可在磁悬浮系统的平衡点s0附近对其进行线性化处理；通过式(1)，可得被控对象中磁悬浮系统的运动微分方程如下：
其中，
接着，通过式(2)，可以得到被控对象的状态空间方程如下：
其中，即小钢球的位移和速度，u为控制增益，f为外部扰动和模型误差，y为观测输出，Ag为动态矩阵，Bu为输入矩阵，Bf为外部扰动和模型误差的常数矩阵，Cg为观测输出矩阵；
(1.3)使用非线性滤波器消除电磁力和模型误差对被控对象的影响；具体步骤如下：
(1.3.1)磁悬浮系统中的非线性滤波器如下：
其中，z1为钢球位移的估计量，z2为钢球速度的估计量，z3为电磁力和模型误差f的估计量，e为输出误差，ω为系统控制带宽，bu为已知系统参数，y为系统输出，u0为控制增益，Kd为干扰补偿增益；
(1.3.2)通过非线性滤波器对被控对象电磁力和模型误差f的估计量z3，可由式(4)，将磁悬浮系统中电磁力和模型误差Bff消除，同时系统在能量有界的未知扰动和噪声d的影响下的状态空间方程如下所示：
其中，z为外部扰动的输出，Cz为输出矩阵，Dzu和Dzd为常数矩阵；
2)、在磁悬浮系统中得到被控对象与给定跟踪信号所组成的扩阶系统，并通过求解线性矩阵方程组构造一个结构化控制器，实现对给定信号的跟踪；具体包括：
(2.1)在磁悬浮系统中，获取被控对象与给定跟踪信号组成的闭环状态空间方程；具体包括：
(2.1.1)引入给定的跟踪信号，其状态空间方程如下：
其中，xω为给定信号的状态量，ξ为给定信号的输入，ωi为给定信号的输出，Ai为动态矩阵，Bi为输入矩阵，Ci为输出矩阵，Di为常数矩阵；
(2.1.2)将跟踪信号作用于被控对象，由于d为能量有界的外部扰动，不会影响对系统对跟定信号的跟踪，为方便叙述，在这部分可不做考虑，所以由式(5)、(6)，使被控对象的状态空间方程变为如下：
其中，e0为系统的控制输出，即跟踪误差，Ce为跟踪误差的输出矩阵，Dω和Du为跟踪误差的常数矩阵，Dgω为观测输出的常数矩阵；
(2.1.3)通过式(6)、(7)，可以得到被控对象与跟踪信号的闭环系统如下：
其中，A11＝Ag，B12＝Bu，C11＝Ce，C12＝DωCi，C21＝Cg，C22＝DgωCi，D11＝DωDi，D12＝Du，D21＝DgωDi；
(2.2)通过求解线性矩阵方程组得到一个结构化控制器，获得被控对象与给定跟踪信号组成的扩阶系统；具体包括：
(2.2.1)在磁悬浮系统中，构造一个结构化控制器：
其中，为动态矩阵，为输入矩阵，为输出矩阵，Ak，Bk，均为常数矩阵，且Πi和Xi为线性矩阵方程组的解；
(2.2.2)通过上述步骤可得磁悬浮系统中被控对象与给定跟踪信号组成的扩阶系统如下：
其中，为此扩阶系统的状态量，eG为此扩阶系统的控制输出，yG为此扩阶系统的观测输出；
3)、在能量有界的未知扰动和噪声的影响下，通过摄动法修正标准型线性矩阵不等式组，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；具体包括：
(3.1)在能量有界的未知扰动和噪声d的影响下，通过式(10)，可得到d与磁悬浮系统所组成的闭环系统状态空间方程如下：
其中，z为外部扰动的输出，yG为观测输出，为动态矩阵，为输入矩阵，为观测输出矩阵，为输入矩阵，为输出矩阵，和为常数矩阵；
(3.2)根据磁悬浮系统的抗干扰轨迹跟踪降阶控制方法，通过摄动法修正标准型线性矩阵不等式组并对其解值进行迭代；具体包括：
(3.2.1)定义：
其中，T表示矩阵的转置，⊥表示矩阵的核空间，I为单位矩阵，γ为给定的性能指标常数，P和Q为待定未知矩阵；
常数εB＞0，εr＞0，满足：且||N1+N2||表示矩阵的欧式范数，通过式(12)～(14)，可得到线性矩阵不等式组如下：
(3.2.2)对线性矩阵不等式组的解进行迭代，具体步骤如下：
(3.2.2.a)定义矩阵：M＝-N2+N1，对M进行奇异值分解，可得：
M＝ZΣV* (16)
其中，*表示矩阵的共轭转置，Σ为奇异值矩阵，Z和V为酉矩阵；
(3.2.2.b)进一步，定义：
其中，P0＝P，+表示矩阵的Moore-Penrose逆，W11，W12和W22为常数矩阵；
(3.2.2.c)将P1-Q-1进行特征值分解，可得到如下：
P1-Q-1＝Θdiag(λ1,...λl,λl+1,...λn)ΘT (19)
其中，Θ为酉矩阵，diag()表示对角矩阵，λ1，λl，λl+1，λn为矩阵P1-Q-1的特征值并且为降序排列；
(3.2.2.d)更近一步，将满足εr≥λl+1≥...≥λn的特征值置为0，可通过式(19)，得到并得到如下条件：
其中，rank()表示矩阵的秩，若所求的能使式(20)的条件成立，则满足方法要求；否则，令P0＝P1，重复式(17)～(19)，直到结果满足式(20)的条件，求得和Q；
(3.3)通过上述所求和Q，构造降阶H无穷控制器改善系统抑制未知扰动和噪声的能力；具体步骤如下：
(3.3.1)定义：
其中， M和N为常数矩阵且满足
通过式(21)～(22)，可以得到一个线性矩阵不等式如下：
其中，为控制器中的参数矩阵，且
所以，可通过上述步骤获得控制器参数，并最后实现了磁悬浮系统的抗干扰轨迹跟踪降阶控制方法。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611145515.6A CN106527144B (zh) | 2016-12-13 | 2016-12-13 | 磁悬浮系统的抗干扰轨迹跟踪降阶控制方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201611145515.6A CN106527144B (zh) | 2016-12-13 | 2016-12-13 | 磁悬浮系统的抗干扰轨迹跟踪降阶控制方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106527144A CN106527144A (zh) | 2017-03-22 |
| CN106527144B true CN106527144B (zh) | 2019-02-01 |
Family
ID=58342421
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201611145515.6A Active CN106527144B (zh) | 2016-12-13 | 2016-12-13 | 磁悬浮系统的抗干扰轨迹跟踪降阶控制方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106527144B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110794680B (zh) * | 2019-11-14 | 2022-04-05 | 浙江工业大学 | 一种基于扩张状态观测器的磁悬浮球系统预测跟踪控制方法 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP4395556B2 (ja) * | 2004-11-22 | 2010-01-13 | 独立行政法人 宇宙航空研究開発機構 | 磁力支持装置の制御定数自動調整方法 |
| CN102707626B (zh) * | 2012-06-25 | 2015-03-18 | 北京航空航天大学 | 一种基于位移刚度力超前补偿的自动平衡磁悬浮转子系统 |
| CN104503237B (zh) * | 2014-12-15 | 2017-02-08 | 北京航空航天大学 | 一种基于傅里叶变换的磁悬浮转子谐波振动控制方法 |
| CN106200380B (zh) * | 2016-07-15 | 2019-01-08 | 浙江工业大学 | 基于非线性观测的磁悬浮系统跟踪控制方法 |
-
2016
- 2016-12-13 CN CN201611145515.6A patent/CN106527144B/zh active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN106527144A (zh) | 2017-03-22 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Fei et al. | Adaptive sliding mode control of dynamic systems using double loop recurrent neural network structure | |
| CN108762096B (zh) | 一种基于离散型非线性级联扩张状态观测器的控制力矩陀螺框架系统扰动抑制方法 | |
| CN105182984B (zh) | 飞行器俯仰姿态的线性自抗扰控制器设计与参数整定方法 | |
| CN110119087B (zh) | 一种有向通信下二阶多智能体系统固定时一致性跟踪方法 | |
| CN107491081B (zh) | 一种抗干扰四旋翼无人机姿态控制方法 | |
| CN109333529A (zh) | 含预定义性能的多单臂机械手输出一致控制器及设计方法 | |
| CN106200380B (zh) | 基于非线性观测的磁悬浮系统跟踪控制方法 | |
| CN115773569B (zh) | 基于自抗扰解耦的海洋平台通风系统风量控制方法 | |
| CN106527144B (zh) | 磁悬浮系统的抗干扰轨迹跟踪降阶控制方法 | |
| CN108775370A (zh) | 一种磁流变阻尼器的磁滞补偿控制方法 | |
| Sun et al. | A design for two-wheeled self-balancing robot based on Kalman filter and LQR | |
| Ouyang et al. | Adaptive LADRC parameter optimization in magnetic levitation | |
| Chen et al. | Model-guided data-driven decentralized control for magnetic levitation systems | |
| CN112436166A (zh) | 质子交换膜燃料电池空气供给系统的滑模控制方法 | |
| CN111692210A (zh) | 一种基于查表法的主动磁轴承控制器的构造方法 | |
| CN107272416A (zh) | 一类线性参数变化系统动态量化h∞控制方法 | |
| CN111025904A (zh) | 一种基于扩张状态观测器的磁悬浮球连续滑模控制方法 | |
| CN110095985A (zh) | 一种观测器设计方法和抗干扰控制系统 | |
| Al-Bakri et al. | Adaptive Model Predictive Control for a Magnetic Suspension System under Initial Position Dispersions and Voltage Disturbances. | |
| Wang et al. | Application of H∞ control based on mixed sensitivity in the electro-hydraulic load simulator | |
| CN115202209A (zh) | 径向混合磁轴承的线性自抗扰控制器消除扰动方法 | |
| Zhang et al. | Measurement and control of magnetic flux signal in a maglev system | |
| Chen et al. | Design of magnetic levitation ball control based on co-simulation of SIMULINK and ADAMS | |
| Huang et al. | Brain emotional learning based intelligent controller for nonlinear system | |
| Shen et al. | An improved chattering-free sliding mode control with finite time convergence for reentry vehicle |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | 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 |