1.飞轮储能用磁悬浮轴承时滞系统变转速切换扰动抑制方法，其特征在于：包括，将时滞因素τ纳入考虑构建磁悬浮轴承闭环控制系统微分方程组，进而转化为状态方程形式的时滞模型；

磁悬浮轴承转子系统在不同转速范围下的受力分析描述成马尔可夫链的统计特性，基于所述时滞模型采用Markov随机跳变过程来对磁悬浮轴承系统进行建模，将全转速划分为低速段、中速段和高速段，依次分别对应磁悬浮轴承Markov系统描述的三个模态，得到磁悬浮轴承系统的时滞随机跳变模型：式中，x(t)为磁悬浮轴承系统的状态量位移及位移变化率，u(t)为磁悬浮轴承系统的输入量电磁铁控制电流，v(t)为磁悬浮轴承系统的质量不平衡扰动，x(t‑τ)为磁悬浮轴承系统各环节引起的时滞状态量，z(t)为系统性能输出，A(rt)、B(rt)、E(rt)、K(rt)和C(rt)都是系统矩阵；

基于上述时滞随机跳变模型利用H∞控制理论推导三种模态下相应线性矩阵不等式的满足条件，在该满足条件下磁悬浮轴承系统在全转速变化范围内随机稳定，通过统计获取上述三个模态相互切换的基本转移概率矩阵，将磁悬浮轴承系统的自身参数和相应的基本转移概率矩阵引入线性矩阵不等式的满足条件，求取对应三个模态的H∞状态反馈控制器的参数——比例时间系数和微分系数，将上述H∞状态反馈控制器的参数带入所述时滞随机跳变模型得到模态依赖H∞状态反馈控制器；

所述时滞因素τ＝τc+τa+τm+τs，其中τc为控制器滞后时间系数，τa为功率放大器滞后时间系数，τm为执行器滞后时间系数，τs为传感器滞后时间系数，磁悬浮轴承系统属于高速动态系统，控制器时滞系数τc和执行器滞后时间系数τm都是时间变量，因此，将磁悬浮轴承系统时滞参数τ(t)定义为有界时变函数：τ(t)＝τ0+η(t)，|η(t)|＜σ，式中τ0、σ和μ都为常数，可通过对反馈回路中时滞进行测算来确定，则有τ(t)∈[h1，h2]，其中h1为系统时滞的下界值，h2为系统时滞的上界值；

rt＝1即模态1，rt＝2即模态2，rt＝3即模态3，模态1对应于磁悬浮系统高速段，模态2对应于磁悬浮系统中速段，模态3对应于磁悬浮系统低速段，π11表示转子转速保持在高速段的概率，π12表示转子转速由高速段跳变到中速段的概率，π21表示转子转速由中速段提高到高速段的概率，π22表示转子转速保持在中速段的概率，π23表示转子转速从中速段跳变到低速段的概率，π33表示转子转速保持在低速段的概率，π32表示转子转速由低速段提高到中速段的概率，ω1、ω2和ω3分别表示转子低速阈值、中速阈值和高速阈值，磁悬浮轴承系统的时滞随机跳变模型的基本转移概率矩阵可表示为：由于飞轮转子的转速的变化是一个连续的过程，其不可能从低速直接跳变到高速状态，也不可能从高速直接跳变到低速状态，因此π13＝π31＝0；

线性矩阵不等式的满足条件包括：给定标量0≤h1≤h2，γ＞0，τ1i，τ2i，τ3i，对于转移概率已知的磁悬浮轴承系统的时滞随机跳变模型若具有对称矩阵Pi＞0，Q1i＞0，Q1＞0，Q2＞0，Q3＞0，R1＞0，R2＞0和矩阵Hi，Mi，Ni，Wi，Si满足下式线性矩阵不等式：当d(t)＝h1，a＝2，b＝1时

当d(t)＝h2，a＝1，b＝2时

式中

ψ112＝‑τ1iBiWi+HiBi，

ψ24＝‑s11+s13‑2bR2+s12‑s14，

ψ26＝‑2s12+2s14+6bR2，

ψ210＝‑s11+s13，

ψ211＝‑s12+s14，

ψ33＝‑Q1‑R1‑4aR2，

ψ34＝s11+s13‑s12‑s14，

ψ35＝2s12+2s14，

ψ310＝s11+s13，

ψ311＝s12+s14，

ψ44＝‑Q2‑4bR2，

ψ712＝‑τ2iBiWi+MiBi，

ψ1212＝He(‑τ3iWi)，

‑1

则磁悬浮轴承系统的在全转速变化范围内随机稳定，且满足H∞性能指标，Ki＝Wi Ni。

2.根据权利要求1所述的飞轮储能用磁悬浮轴承时滞系统变转速切换扰动抑制方法，其特征在于：当rt＝1即模态1时，磁悬浮轴承系统低速运行，运动学模型中既不考虑不平衡振动力也不考虑陀螺效应，此时磁悬浮轴承系统的时滞随机跳变模型为：u1(t)＝K1x(t‑τ)

z1(t)＝C1x(t)

当rt＝2即模态2时，磁悬浮轴承系统中速运行，这时受力分析仅考虑不平衡振动力，此时磁悬浮轴承系统的时滞随机跳变模型为：u2(t)＝K2x(t‑τ)

z2(t)＝C2x(t)

当rt＝3即模态3时，磁悬浮轴承系统高速运行，转子系统受力分析同时考虑不平衡振动力和陀螺效应，此时磁悬浮轴承系统的时滞随机跳变模型为：