受非对称死区输入影响的分数阶系统反步滑模控制方法与流程

1.本发明属于室内外无缝定位技术领域，特别是涉及受非对称死区输入影响的分数阶系统反步滑模控制方法。


背景技术：

2.随着分数阶微积分理论的发展，越来越多的学者开始在分数阶系统的控制中讨论非光滑非线性特性的影响，这些非线性特性的存在会使得系统性能变差，甚至会使系统变得不稳定。采用分数阶微积分对实际工程系统建模相比整数阶微积分理论有更精确的效果，由于分数阶算子的奇异性与复杂性，目前很多关于分数阶系统的控制研究中很难像整数阶系统构建统一形式的lyapunov函数，导致分数阶系统的稳定性分析一直是分数阶系统研究中的重点和难点，除此之外，很多分数阶系统的研究中，系统参数的事先假设已知，事实上，因为外界环境的干扰因素及建模的复杂性，系统的参数不可事先完全已知，因此，有关分数阶系统的自适应控制问题成为亟待解决的问题。反步法是一种设计控制器的迭代方法，该方法通过设计虚拟控制器及部分lyapunov函数迭代分析设计出实际控制器形式，所设计的控制器可以确保系统的全局稳定性，提升系统暂态性能。
3.基于以上背景，本发明基于反步法并结合滑模控制技术，考虑非对称死区输入对系统的影响，构建辅助分数阶系统产生虚拟信号抵消非线性输入的影响，同时设计虚拟控制器，结合分数阶频率分布模型，应用间接lyapunov稳定性分析理论分别验证趋近阶段和滑模阶段稳定性，仿真实例证明，本发明提出的一种受非对称死区输入影响的分数阶系统反步滑模控制方法可实现系统的自适应镇定控制，使得系统具有很好的抗干扰性能，鲁棒性强，系统的全部未知参数都得到有效辨识。


技术实现要素：

4.针对以上不足，为了实现受非对称死区输入影响的分数阶系统的自适应镇定控制，本发明基于反步法并结合滑模控制方案对分数阶系统的自适应镇定研究提出一种新的控制策略。该方案可以有效避免非线性输入对系统产生的影响，同时对系统未建模动态和所受到的外界干扰也能进行有效抑制。分数阶频率分布模型的应用可以将滑模控制中的趋近阶段和滑模阶段状态等效为无限维常微分方程，从而采用间接lyapunov稳定理论可以很好地验证系统全局稳定性。通过反步迭代设计虚拟控制器最终得到合适的实际控制器形式，实现系统全局镇定及全部未知参数的辨识。
5.本发明提供受非对称死区输入影响的分数阶系统反步滑模控制方法，具体包含如下步骤：
6.(1)确定一类具有非对称死区输入的严格反馈结构分数阶系统数学模型；
7.所述步骤(1)具有非对称死区输入的严格反馈结构分数阶系统数学模型为：
[0008][0009]
式(1)中，x＝(x1,x2,...,x
n
)
t
为系统状态变量；δ
i
为系统第i个方程参数矢量；g
i
(
·
)，f
i
(
·
)和f
i
(
·
)i＝1,2,...,n为已知的光滑非线性函数；
[0010]
所涉及如下形式的一类具有非对称死区输入的典型严格反馈分数阶系统：
[0011][0012][0013][0014][0015][0016]
式(2)典型结构中，α∈(0，1)，δf
i
(x)和d
i
(t)(i＝1,2,...,n)分别为未建模动态和外界干扰项，ψ(u(t))为非对称死区输入，描述如下：
[0017][0018]
其中，ρ

，ρ
‑
，u

，u
‑
为正实数，斜率参数ρ

和ρ
‑
有界，存在已知参数ρ1和ρ2满足：max{ρ

,ρ
‑
}＝ρ1，min{ρ

,ρ
‑
}＝ρ2，更进一步，非对称死区函数可改写为：
[0019]
ψ(u(t))＝ρ(t)u(t) δu(t)
ꢀꢀꢀ
(4)
[0020]
其中
[0021][0022][0023]
(2)为正确应用反步滑模控制方案，选择转换变量构建新系统；
[0024]
所述步骤(2)选择转换变量构建新系统，其选择如下：
[0025][0026]
式(7)中，i＝2,3,...,n，θ
j
，j＝1,2,...,n
‑
1为接下来设计的虚拟控制器，σ
j
，j＝1,2,...,n为下述辅助分数阶系统产生的虚拟信号；
[0027]
(3)为补偿非对称死区输入带来的不利影响，构建一个分数阶辅助系统；
[0028]
步骤(3)为抵消非对称死区输入带来的不利影响，构建如下形式的辅助分数阶系
统以产生上文提及的虚拟信号σ
j
：
[0029][0030]
其中c
i
＞0，c
n
＞0，所涉及的系统未建模动态和外界干扰有界且上界未知，即
[0031][0032]
式(9)中，β
i1
，β
i2
为未知正实数；
[0033]
(4)构造合适的分数阶形式的滑模面，并求导得到分数阶滑模态及虚拟控制器；
[0034]
所述步骤(4)，提及的分数阶滑模面具有以下形式：
[0035][0036]
其中，p＝1,2,...,n，通过对滑模变量s求一阶导数，当系统运行至滑模面时，得到如下形式的滑模态：
[0037]
d
α
ξ
p
＝
‑
ξ
p
‑
sgn(ξ
p
)
ꢀꢀꢀ
(11)
[0038]
根据期望滑模态动态描述方程，得到如下形式的虚拟控制器：
[0039][0040]
其中，分别为β
i1
，β
i2
估计值，i＝1,2,...,n；
[0041]
(5)基于滑模变量提出系统未知参数、未建模动态和外界干扰未知上界的自适应辨识律；
[0042]
步骤(5)，记为参数估计误差，提出如下形式未知参数辨识律：
[0043][0044]
(6)基于转换变量构建的新系统，建立分数阶频率分布模型，应用间接lyapunov稳定性分析理论，分析每个子系统稳定性，并得到实际控制器形式；
[0045]
步骤(6)，根据转换变量(7)和分数阶辅助系统(8)，得到如下形式的第一个新建子系统
[0046][0047]
将式(14)与参数辨识律(13)转换为分数阶频率分布模型
[0048][0049]
选择如下形式的lyapunov函数，分析第一个新建子系统稳定性；
[0050][0051]
对式(16)等号两侧分别对时间求一阶导数，并根据所设计的滑模面、不确定项上界、未知参数辨识律，推导得到
[0052][0053]
显然，经过间接lyapunov稳定理论验证，第一个新建子系统是渐近稳定的，同理，采用相同操作，前n
‑
1个子系统都是渐近稳定的，当所设计的实际控制器u(t)具有如下形式时，第n个子系统也为渐近稳定：
[0054][0055]
构建如下形式的lyapunov函数验证第n个子系统稳定性
[0056][0057]
对式(19)等号两侧分别对时间求一阶导数得
[0058][0059]
根据非对称死区函数的特点，存在一个分段正函数a(t)，使得ρ(t)/ρ2＝1 a(t)，根据所设计的实际控制器形式(18)，则有：
[0060][0061]
结合式(20)和式(21)，得
[0062][0063]
综上，所有转换变量构成的新子系统在实际控制器作用下都是渐近稳定的；
[0064]
(7)根据间接lyapunov稳定理论，构建合适稳定性分析函数，验证滑模态滑模阶段的稳定性，实现系统的自适应镇定控制。
[0065]
作为本发明进一步改进，步骤(7)，采用间接lyapunov稳定理论验证期望滑模态的稳定性，根据第一个期望滑模态并结合分数阶频率分布模型，得
[0066][0067]
对上述滑模态(23)，选择如下形式的lyapunov函数验证其稳定性
[0068][0069]
对式(24)等号两侧对时间求一阶导数，得
[0070][0071]
显然第一个期望滑模态渐近稳定，采用相同推导，得到第n个期望滑模态所对应lyapunov函数为
[0072][0073]
对式(26)等号两侧对时间求一阶导数，得
[0074]
[0075]
以上推导结果意味着n个期望滑模态均为渐近稳定，即在所设计的控制器u(t)的作用下，一类受非对称死区输入影响的分数阶系统可实现全局自适应镇定控制。
[0076]
本发明与现有技术相比，其显著优点为：
[0077]
(1)本发明首次将反步法和滑模控制技术进行结合，构建分数阶辅助系统，产生虚拟信号抵消非对称死区输入产生的不利影响，实现具有严格反馈结构的分数阶系统的全局镇定控制。
[0078]
(2)本发明提出的自适应镇定控制方案，可以很好的辨识系统全部未知参数、未建模动态和外界干扰项的上界，通过调整滑模面参数可以有效调节参数辨识速度，提高系统鲁棒性。
[0079]
(3)本发明基于分数阶频率分布模型，应用间接lyapunov稳定理论可以分别验证滑模控制趋近阶段和滑模阶段的渐近稳定性，可以很好避免分数阶算子难求整数阶导数问题，迭代得到的实际控制器实用性强，控制效果达到期望要求。
附图说明
[0080]
图1为本发明提出的一种受非对称死区输入影响的分数阶系统反步滑模控制方法的流程图；
[0081]
图2为非对称死区输入函数的结构图；
[0082]
图3为分数阶genesio
‑
tesi系统奇异吸引子；
[0083]
图4为分数阶genesio
‑
tesi系统未引入控制器前的状态轨迹曲线图；
[0084]
图5为激活控制器后转换变量构建的新系统的状态轨迹曲线图。
具体实施方式
[0085]
下面结合附图与具体实施方式对本发明作进一步详细描述：
[0086]
本发明基于反步法并结合滑模控制方案对分数阶系统的自适应镇定研究提出一种新的控制策略。该方案可以有效避免非线性输入对系统产生的影响，同时对系统未建模动态和所受到的外界干扰也能进行有效抑制。分数阶频率分布模型的应用可以将滑模控制中的趋近阶段和滑模阶段状态等效为无限维常微分方程，从而采用间接lyapunov稳定理论可以很好地验证系统全局稳定性。通过反步迭代设计虚拟控制器最终得到合适的实际控制器形式，实现系统全局镇定及全部未知参数的辨识。
[0087]
下面结合附图对本发明的应用原理作详细的说明。
[0088]
如图1所示,一种受非对称死区输入影响的分数阶系统反步滑模控制方法，主要包括以下步骤：
[0089]
针对具有严格反馈结构的分数阶系统数学描述方程，明确方程中各部分代表的含义，从系统实际运行环境出发，在模型中充分考虑未建模动态、外界干扰、非对称死区输入项，使得研究问题具有重要的实际意义；
[0090]
为补偿非对称死区输入给系统带来的不利影响，构建一个辅助分数阶系统产生虚拟信号来实现。选择合适的转换变量，构建新的子系统，以方便后续操作研究；
[0091]
根据所构建的新转换变量系统和辅助分数阶系统，设计分数阶形式的滑模面，基于分数阶分布频率模型，引入间接lyapunov稳定理论，逐步探讨转换变量构造的每个新子
系统的稳定性，从而得到虚拟控制器形式；
[0092]
为充分辨识系统未知参数、未建模动态和外界干扰项的未知上界，结合滑模变量和转换变量设计自适应辨识律，实现以上全部未知参数的辨识；
[0093]
通过对n个新子系统进行迭代验证其稳定性，得到综合的实际控制器形式，该控制器有很好的抑制外界干扰和系统不确定性的作用；
[0094]
实现对转换变量构造的新子系统的稳定性验证，也即完成滑模运动趋近阶段稳定性的验证，同样可将期望滑模态借助于分数阶频率分布模型转换为无限维常微分方程，采用间接lyapunov稳定理论验证滑模运动滑模阶段稳定性。实验结果表明本发明所提出控制方案是合理且有效的。
[0095]
针对一类具有严格反馈结构的分数阶系统，其数学模型描述为：
[0096][0097]
式(1)中，x＝(x1,x2,...,x
n
)
t
为系统状态变量；δ
i
为系统第i个方程参数矢量；g
i
(
·
)，f
i
(
·
)和f
i
(
·
)(i＝1,2,...,n)为已知的光滑非线性函数。本发明主要涉及如下形式的一类具有非对称死区输入的典型严格反馈分数阶系统：
[0098][0099]
式(2)典型结构中，α∈(0，1)，δf
i
(x)和d
i
(t)(i＝1,2,...,n)分别为未建模动态和外界干扰项，ψ(u(t))为非对称死区输入，描述如下：
[0100][0101]
其中，ρ

，ρ
‑
，u

，u
‑
为正实数，斜率参数ρ

和ρ
‑
有界，存在已知参数ρ1和ρ2满足：max{ρ

,ρ
‑
}＝ρ1，min{ρ

,ρ
‑
}＝ρ2，更进一步，非对称死区函数可改写为：
[0102]
ψ(u(t))＝ρ(t)u(t) δu(t)
ꢀꢀꢀ
(4)
[0103]
其中
[0104][0105]
[0106]
选择转换变量构建新系统，其选择如下：
[0107][0108]
式(7)中，i＝2,3,...,n，θ
j
(j＝1,2,...,n
‑
1)为接下来设计的虚拟控制器，σ
j
(j＝1,2,...,n)为下述辅助分数阶系统产生的虚拟信号。
[0109]
为抵消非对称死区输入带来的不利影响，构建如下形式的辅助分数阶系统以产生上文提及的虚拟信号σ
j
：
[0110][0111]
其中c
i
＞0，c
n
＞0。本发明中所涉及的系统未建模动态和外界干扰有界且上界未知，即
[0112][0113]
式(9)中，β
i1
，β
i2
为未知正实数。
[0114]
本发明中提及的分数阶滑模面具有以下形式：
[0115][0116]
其中，p＝1,2,...,n，通过对滑模变量s求一阶导数，当系统运行至滑模面时，得到如下形式的滑模态：
[0117]
d
α
ξ
p
＝
‑
ξ
p
‑
sgn(ξ
p
)
ꢀꢀꢀ
(11)
[0118]
根据期望滑模态动态描述方程，可得到如下形式的虚拟控制器：
[0119][0120]
其中，分别为β
i1
，β
i2
估计值，i＝1,2,...,n。
[0121]
记为参数估计误差，提出如下形式未知参数辨识律：
[0122][0123]
根据转换变量(7)和分数阶辅助系统(8)，得到如下形式的第一个新建子系统
[0124][0125]
将式(14)与参数辨识律(13)转换为分数阶频率分布模型
[0126][0127]
选择如下形式的lyapunov函数，分析第一个新建子系统稳定性
[0128]
对式(16)等号两侧分别对时间求一阶导数，并根据所设计的滑模面、不确定项上界、未知参数辨识律，可推导得到
[0129][0130]
显然，经过间接lyapunov稳定理论验证，第一个新建子系统是渐近稳定的，同理，采用相同操作，前n
‑
1个子系统都是渐近稳定的。当所设计的实际控制器u(t)具有如下形式时，第n个子系统也为渐近稳定：
[0131][0132]
构建如下形式的lyapunov函数验证第n个子系统稳定性
[0133]
[0134]
对式(19)等号两侧分别对时间求一阶导数可得
[0135][0136]
根据非对称死区函数的特点，存在一个分段正函数a(t)，使得ρ(t)/ρ2＝1 a(t)，根据所设计的实际控制器形式(18)，则有：
[0137][0138]
结合式(20)和式(21)，可得
[0139][0140]
综上，所有转换变量构成的新子系统在实际控制器作用下都是渐近稳定的。
[0141]
采用间接lyapunov稳定理论验证期望滑模态的稳定性，根据第一个期望滑模态形式并结合分数阶频率分布模型，可得
[0142][0143]
对上述滑模态(23)，选择如下形式的lyapunov函数验证其稳定性
[0144][0145]
对式(24)等号两侧对时间求一阶导数，可得
[0146][0147]
显然第一个期望滑模态渐近稳定，采用相同推导，可得到第n个期望滑模态所对应lyapunov函数为
[0148][0149]
对式(26)等号两侧对时间求一阶导数，可得
[0150][0151]
以上推导结果意味着n个期望滑模态均为渐近稳定，即在所设计的控制器u(t)的作用下，一类受非对称死区输入影响的分数阶系统可实现全局自适应镇定控制。
[0152]
图2所示为本发明中所提到的非对称死区非线性函数的结构示意图。
[0153]
图3
‑
5所示为本发明实例所涉及的吸引子、控制器激活前的原系统状态轨迹时间响应曲线、控制器激活后的新转换变量系统的时间响应曲线，在本发明的实施例中，选择分数阶genesio
‑
tesi系统为被控对象，数学描述为：
[0154][0155]
其中，a1＝1，a2＝1.1，a3＝
‑
0.232，a4＝1，δ3＝[a1,a2,a3,a4]
t
，δf(x)＝
‑
0.01cos(x3)，d(t)＝0.02sin(3t)，ψ(u(t))为非对称死区输入，描述如下：
[0156][0157]
控制方案参数m1＝m2＝m3＝5，c1＝c2＝c3＝2，η
31
＝5，η
32
＝2，系统初始条件为：x1(0)＝
‑
0.3，x2(0)＝0.1，x3(0)＝
‑
0.2，0.2，α＝0.8。分数阶genesio
‑
tesi系统混沌吸引子如图3所示，未激活前系统的状态轨迹如图4所示，显然，系统不稳定。当激活控制器u(t)，转换变量构成的新系统状态轨迹如图5所示，响应曲线显示本发明提出的控制方案可以有效实现一类受非对称死区输入影响的分数阶系统的全局镇定控制，控制方案有效可行。
[0158]
以上所述，仅是本发明的较佳实施例之一，并非是对本发明作任何其他形式的限制，而依据本发明的技术实质所作的任何修改或等同变化，仍属于本发明所要求保护的范围。