CN108241297B - 一类奇异时滞系统的渐近稳定性判别方法、维持方法 - Google Patents

一类奇异时滞系统的渐近稳定性判别方法、维持方法 Download PDF

Info

Publication number
CN108241297B
CN108241297B CN201810006540.9A CN201810006540A CN108241297B CN 108241297 B CN108241297 B CN 108241297B CN 201810006540 A CN201810006540 A CN 201810006540A CN 108241297 B CN108241297 B CN 108241297B
Authority
CN
China
Prior art keywords
time
lag
stability
singular
condition
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
CN201810006540.9A
Other languages
English (en)
Other versions
CN108241297A (zh
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.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
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 University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201810006540.9A priority Critical patent/CN108241297B/zh
Publication of CN108241297A publication Critical patent/CN108241297A/zh
Application granted granted Critical
Publication of CN108241297B publication Critical patent/CN108241297B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • 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)
  • Supply And Distribution Of Alternating Current (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一类奇异时滞系统的渐近稳定性判别方法、维持方法，涉及自动控制技术领域。解决了现有技术存在无法全面的描述该类奇异时滞依赖系统的稳定性的技术问题。该奇异时滞系统的渐近稳定性判别方法，包括：利用数学分解方法对奇异时滞系统的时滞稳定性进行分析，以得出维持一类奇异时滞依赖系统稳定性的条件；依据维持该类奇异时滞依赖系统稳定性的条件，判别该类奇异时滞依赖系统的稳定性。该奇异时滞系统的渐近稳定性维持方法，包括：利用维持本发明提供的一类奇异时滞依赖系统稳定性的条件，以提高该类奇异时滞系统的渐近稳定性。本发明用于全面的描述并提高一类奇异时滞依赖系统的稳定性。

Description

一类奇异时滞系统的渐近稳定性判别方法、维持方法
技术领域
本发明涉及自动控制技术领域，尤其涉及一类奇异时滞系统的渐近稳定性判别方法以及该类奇异时滞系统的渐近稳定性维持方法。
背景技术
在各类物理、生物和工业系统中时滞现象是普遍存在的，时滞的存在使得系统分析和综合变得更加复杂,同时时滞的存在也是系统不稳定和系统性能变差的重要原因。奇异系统是一类比正常系统更具广泛形式的动力系统，比正常系统能更好地描述物理系统,奇异系统模型广泛存在于社会生产的各个领域中，如：电力系统、经济系统、机器人系统、宇航工程和生物系统。因此，奇异时滞系统在关于控制问题的数学模型中有重要作用。
根据依赖于时滞的大小可以把奇异时滞系统的稳定条件分成两类，不包含时滞信息的系统称为时滞独立系统，包含时滞信息的系统称为时滞依赖系统。
现有技术中奇异时滞系统的时滞独立稳定性已被广泛研究，例如：使用特征函数方法，对时滞独立稳定性的条件进行了讨论。
本申请人发现：现有技术至少存在以下技术问题：
现有技术中奇异时滞依赖系统的稳定性无法完全描述，由此无法准确地判断奇异时滞依赖系统的稳定性，也无法有效地提高奇异时滞依赖系统的稳定性。
发明内容
本发明的至少一个目的是提出一类奇异时滞系统的渐近稳定性判别方法以及该类奇异时滞系统的渐近稳定性维持方法，解决了现有技术存在无法全面的描述该类奇异时滞依赖系统的稳定性的技术问题。
为实现上述目的，本发明提供了以下技术方案：
本发明实施例提供的一类奇异时滞系统的渐近稳定性判别方法，包括以下步骤:
步骤A、利用数学分解方法对一类奇异时滞系统的时滞稳定性进行分析，以得出维持该类奇异时滞依赖系统稳定性的条件；
步骤B、依据所述维持一类奇异时滞依赖系统稳定性的条件，判别所述奇异时滞依赖系统的稳定性。
优选地，所述数学分解方法为Schur分解方法。
优选地，所述条件为充分必要条件。
优选地，所述一类奇异时滞系统为线性时滞系统。
优选地，所述步骤A，包括以下步骤：
步骤A1、得到非奇异时滞系统的稳定的条件；
步骤A2、将非奇异时滞系统的稳定的条件与Schur分解方法相结合来分析一类奇异时滞系统的时滞稳定性，以得出维持该类奇异时滞依赖系统稳定性的条件。
优选地，所述步骤A2，包括以下步骤:
步骤A21：假设条件det(σE-A)≠0成立，矩阵(δE-A)-1存在，系统
Figure BDA0001538799860000023
其中E,A∈Rn×n，δ∈C，E为奇异矩阵；两边都左乘(δE-A)-1，于是得到
Figure BDA0001538799860000024
式子；
步骤A22：因为矩阵(δE-A)-1存在，令
Figure BDA0001538799860000021
证明
Figure BDA0001538799860000022
故有
Figure BDA0001538799860000031
因此将式子
Figure BDA0001538799860000032
改写成
Figure BDA0001538799860000033
步骤A23：因为
Figure BDA0001538799860000034
为方阵，存在酉矩阵U使得
Figure BDA0001538799860000035
成立；其中T是一个上三角矩阵，它的对角元素是
Figure BDA0001538799860000036
的特征值；选取适当的U阵使得T前m个对角元素都是非零的特征值α1,α2,…,αm，而剩余的对角线元素是零特征值αm+1,αm+2,…,αn；
步骤A24：令x(t)＝Uz(t)则
Figure BDA0001538799860000037
写成
Figure BDA0001538799860000038
的形式；
步骤A25：令z(t)＝[z1(t),....,zm(t),zm+1(t),....zn(t)]T系统
Figure BDA0001538799860000039
的最后一个标量方程仅包含一个变量zn(t)，第n个方程是
Figure BDA00015387998600000310
并得一致性条件zn(t)≡0,t≥-τ；
步骤A26：令zn(t)≡0,t≥-τ，在第n-1个方程中有
Figure BDA00015387998600000311
显然zn-1(t)≡0,t≥-τ是成立的；现在已证得标量元素zn(t)≡zn-1(t)≡0,t≥-τ，同理也证得zk(t)≡0,t≥-τ,k＝m+1,m+2,...,n.
步骤A27：在系统
Figure BDA00015387998600000317
中，对于k＝m+1,m+2,....n；令zk(t)＝0得到一个m维的系统
Figure BDA00015387998600000312
其中R的对角元素α1,α2,…,αm都不为零，并且有
Figure BDA00015387998600000313
亦有[zm+1(t),....,zn(t)]T＝0；
步骤A28：因为矩阵R是非奇异的，系统
Figure BDA00015387998600000314
化为一个时滞系统，当且仅当它的特征方程
Figure BDA00015387998600000315
的所有特征根都在左半平面上时，该系统是渐近稳定的；
当且仅当m个标量方程
Figure BDA00015387998600000316
是渐近稳定的，系统
Figure BDA0001538799860000041
的所有特征根位于左半平面；
步骤A29：有式子x(t)＝Uz(t)＝U[z1(t),...,zm(t),0...,0]T并且矩阵U是非奇异的，系统
Figure BDA0001538799860000042
是渐近稳定的等价于系统
Figure BDA0001538799860000043
是渐近稳定的等价于
Figure BDA0001538799860000044
式中m个标量方程渐近稳定
优选地，所述步骤B，具体为：
若条件det(σE-A)≠0成立，其中：E,A∈Rn×n，δ∈C；令ak≠0，k＝1,2,...,m，ak＝0，k＝m+1,m+2,...,n，并且
Figure BDA0001538799860000045
其中uk,vk∈R；定义αk＝λk(Eδ)；当且仅当对于所有的ak≠0，k＝1,2,...,m，有条件
Figure BDA0001538799860000046
成立时，判定系统
Figure BDA0001538799860000047
是渐近稳定的。
优选地，所述奇异时滞系统为奇异时滞依赖系统。
本发明实施例提供的一类奇异时滞系统的渐近稳定性维持方法，包括以下步骤：
根据本发明任一技术方案提供的一类奇异时滞系统的渐近稳定性判别方法判别一类奇异时滞系统的渐近稳定性；
利用所述维持一类奇异时滞依赖系统稳定性的条件，以提高奇异时滞系统的渐近稳定性。
基于上述技术方案，本发明实施例至少可以产生如下技术效果：
本发明可以利用数学分解方法对一类奇异时滞系统的时滞稳定性进行分析，以得出维持该类奇异时滞依赖系统稳定性的条件，由此全面地描述了该类奇异时滞依赖系统的稳定性(本发明优选方案中使用Schur分解方法对线性时滞系统的时滞稳定性进行分析，给出了一类奇异时滞依赖系统的充分必要条件，结果表明该判别方法能完全描述该类奇异时滞依赖系统稳定性)，然后依据维持该类奇异时滞依赖系统稳定性的条件，判别该类奇异时滞依赖系统的稳定性，解决了现有技术存在无法全面的描述该类奇异时滞依赖系统的稳定性的技术问题。
附图说明
此处所说明的附图用来提供对本发明的进一步理解，构成本申请的一部分，本发明的示意性实施例及其说明用于解释本发明，并不构成对本发明的不当限定。在附图中：
图1为本发明实施例所提供的一类奇异时滞系统的渐近稳定性判别方法的主要流程的示意图；
图2为应用本发明实施例所提供的一类奇异时滞系统的渐近稳定性判别方法的该类奇异时滞系统的渐近稳定性维持方法的主要流程的示意图。
具体实施方式
下面可以参照附图图1～图2以及文字内容理解本发明的内容以及本发明与现有技术之间的区别点。下文通过附图以及列举本发明的一些可选实施例的方式，对本发明的技术方案(包括优选技术方案)做进一步的详细描述。需要说明的是：本实施例中的任何技术特征、任何技术方案均是多种可选的技术特征或可选的技术方案中的一种或几种，为了描述简洁的需要本文件中无法穷举本发明的所有可替代的技术特征以及可替代的技术方案，也不便于每个技术特征的实施方式均强调其为可选的多种实施方式之一，所以本领域技术人员应该知晓：可以将本发明提供的任一技术手段进行替换或将本发明提供的任意两个或更多个技术手段或技术特征互相进行组合而得到新的技术方案。本实施例内的任何技术特征以及任何技术方案均不限制本发明的保护范围，本发明的保护范围应该包括本领域技术人员不付出创造性劳动所能想到的任何替代技术方案以及本领域技术人员将本发明提供的任意两个或更多个技术手段或技术特征互相进行组合而得到的新的技术方案。
本发明实施例提供了一种可以简单、全面地描述一类奇异时滞依赖系统的稳定性，进而可以有效地提高该类奇异时滞依赖系统的稳定性的奇异时滞系统的渐近稳定性判别方法以及该类奇异时滞系统的渐近稳定性维持方法。本发明基于线性时滞系统的稳定性判别方法，使用Schur分解的方法给出了一类奇异时滞系统的渐近稳定性的判别方法，该类奇异时滞系统的稳定性可以被完全描述。
下面结合图1～图2对本发明提供的技术方案进行更为详细的阐述。
如图1所示，本发明实施例所提供的奇异时滞系统的渐近稳定性判别方法，包括以下步骤:
步骤A、利用数学分解方法对奇异时滞系统的时滞稳定性进行分析，以得出维持一类奇异时滞依赖系统稳定性的条件；
步骤B、依据维持一类奇异时滞依赖系统稳定性的条件，判别该类奇异时滞依赖系统的稳定性。
本发明可以利用数学分解方法对一类奇异时滞系统的时滞稳定性进行分析，以得出维持该类奇异时滞依赖系统稳定性的条件，由此全面地描述了该类奇异时滞依赖系统的稳定性(本发明优选方案中使用Schur分解方法对线性时滞系统的时滞稳定性进行分析，给出了一类奇异时滞依赖系统的充分必要条件，结果表明该判别方法能完全描述该类奇异时滞依赖系统稳定性)，然后依据维持该类奇异时滞依赖系统稳定性的条件，判别该类奇异时滞依赖系统的稳定性。
作为可选地实施方式，数学分解方法为Schur分解方法。Schur分解是应用Schur定理(或称：舒尔定理)进行分解的数学方法。
作为可选地实施方式，条件为充分必要条件。必要条件越全面、越充分，则对一类奇异时滞依赖系统的稳定性的描述越全面准确。
作为可选地实施方式，一类奇异时滞系统为线性时滞系统。线性时滞系统更便于采用数学分解方法进行分析。
作为可选地实施方式，步骤A，包括以下步骤：
步骤A1、得到非奇异时滞系统的稳定的条件(结论)；
步骤A2、将非奇异时滞系统的稳定的条件与Schur分解方法相结合来分析一类奇异时滞系统的时滞稳定性，以得出维持该类奇异时滞依赖系统稳定性的条件。
作为可选地实施方式，奇异时滞系统为奇异时滞依赖系统。奇异时滞依赖系统包含时滞信息，其稳定性更具有研究价值。
作为可选地实施方式，步骤A2，包括以下步骤:
步骤A21：假设条件det(σE-A)≠0成立，矩阵(δE-A)-1存在，系统
Figure BDA0001538799860000071
其中E,A∈Rn×n，δ∈C，E为奇异矩阵；两边都左乘(δE-A)-1，于是得到
Figure BDA0001538799860000072
式子；
步骤A22：因为矩阵(δE-A)-1存在，令
Figure BDA0001538799860000073
得出(可以证明)
Figure BDA0001538799860000074
并得出
Figure BDA0001538799860000075
因此将式子
Figure BDA0001538799860000076
改写成
Figure BDA0001538799860000077
步骤A23：因为
Figure BDA0001538799860000078
为方阵，存在酉矩阵U使得
Figure BDA0001538799860000079
成立；其中T是一个上三角矩阵，它的对角元素是
Figure BDA00015387998600000710
的特征值；选取适当的U阵使得T前m个对角元素都是非零的特征值α1,α2,…,αm，而剩余的对角线元素是零特征值αm+1,αm+2,…,αn；
步骤A24：令x(t)＝Uz(t)，则
Figure BDA00015387998600000711
写成
Figure BDA00015387998600000712
的形式；
步骤A25：令z(t)＝[z1(t),....,zm(t),zm+1(t),....zn(t)]T系统
Figure BDA00015387998600000713
的最后一个标量方程仅包含一个变量zn(t)，第n个方程是
Figure BDA0001538799860000081
并得一致性条件zn(t)≡0,t≥-τ；
步骤A26：令zn(t)≡0,t≥-τ，在第n-1个方程中有
Figure BDA0001538799860000082
显然zn-1(t)≡0,t≥-τ是成立的；现在已证得标量元素zn(t)≡zn-1(t)≡0,t≥-τ，同理也证得zk(t)≡0,t≥-τ,k＝m+1,m+2,...,n.
步骤A27：在系统
Figure BDA0001538799860000083
中，对于k＝m+1,m+2,....n；令zk(t)＝0得到一个m维的系统
Figure BDA0001538799860000084
其中R的对角元素α1,α2,…,αm都不为零，并且有
Figure BDA0001538799860000085
亦有[zm+1(t),....,zn(t)]T＝0；
步骤A28：因为矩阵R是非奇异的，系统
Figure BDA0001538799860000086
化为一个时滞系统，当且仅当它的特征方程
Figure BDA0001538799860000087
的所有特征根都在左半平面上时，该系统是渐近稳定的；
当且仅当m个标量方程
Figure BDA0001538799860000088
是渐近稳定的，系统
Figure BDA0001538799860000089
的所有特征根位于左半平面；
步骤A29：有式子x(t)＝Uz(t)＝U[z1(t),...,zm(t),0...,0]T并且矩阵U是非奇异的，系统
Figure BDA00015387998600000810
是渐近稳定的等价于系统
Figure BDA00015387998600000811
是渐近稳定的等价于
Figure BDA00015387998600000812
式中m个标量方程渐近稳定
下面更为详细地说明本发明步骤A2：
步骤A21：如果对于任意的ε＞0，存在δ＝δ(ε)＞0即所有连续函数
Figure BDA00015387998600000813
满足一致性条件并且有
Figure BDA00015387998600000814
对于所有的t≥0，系统
Figure BDA00015387998600000815
的初值问题的解
Figure BDA00015387998600000816
满足
Figure BDA00015387998600000817
则奇异时滞系统的平凡解是稳定的。如果系统
Figure BDA0001538799860000091
是稳定的并且满足
Figure BDA0001538799860000092
则它的平凡解是渐近稳定的。
步骤A22：E,A∈Rn×n，δ∈C。假设条件det(σE-A)≠0成立，故矩阵(δE-A)-1存在。系统
Figure BDA00015387998600000915
其中矩阵E是奇异的，两边都左乘(δE-A)-1，于是得到
Figure BDA00015387998600000916
式子。
因为矩阵(δE-A)-1存在，令
Figure BDA0001538799860000093
可以证明
Figure BDA0001538799860000094
故有
Figure BDA0001538799860000095
因此式子
Figure BDA0001538799860000096
可以改写成
Figure BDA0001538799860000097
步骤A23：如果F∈Cn×n，那么存在酉矩阵U∈Cn×n,有UHFU＝T。其中T∈Cn×n是上三角矩阵，可以选择适当的U使得F的特征值出现在对角元素上的任意位置。
Figure BDA0001538799860000098
为方阵，所以存在酉矩阵U使得
Figure BDA0001538799860000099
成立。其中T是一个上三角矩阵，它的对角元素是
Figure BDA00015387998600000910
的特征值。选取适当的U阵使得T前面的m个对角元素都是非零的特征值α1,α2,…,αm，而剩余的对角线元素是零特征值αm+1,αm+2,…,αn。
步骤A24：令x(t)＝Uz(t),其中U为酉矩阵，有UHU＝I，则
Figure BDA00015387998600000911
可以写成
Figure BDA00015387998600000912
的形式；
步骤A25：令z(t)＝[z1(t),....,zm(t),zm+1(t),....zn(t)]T系统
Figure BDA00015387998600000913
的最后一个标量方程仅包含一个变量zn(t)，第n个方程是
Figure BDA00015387998600000914
并可得一致性条件zn(t)≡0,t≥-τ。
步骤A26：令zn(t)≡0,t≥-τ，在第n-1个方程中有
Figure BDA0001538799860000101
显然zn-1(t)≡0,t≥-τ是成立的。现证得标量元素zn(t)≡zn-1(t)≡0,t≥-τ，同理也可证得zk(t)≡0,t≥-τ,k＝m+1,m+2,...,n.
步骤A27：在系统
Figure BDA0001538799860000102
中，对于k＝m+1,m+2,....n。令zk(t)＝0可以得到一个m维的系统
Figure BDA0001538799860000103
其中R的对角元素α1,α2,…,αm都不为零，并且有
Figure BDA0001538799860000104
亦有[zm+1(t),....,zn(t)]T＝0。
步骤A28：因为矩阵R是非奇异的，系统
Figure BDA0001538799860000105
可化为一个时滞系统，当且仅当它的特征方程
Figure BDA0001538799860000106
的所有特征根都在左半平面上时，该系统是渐近稳定的。
当且仅当m个标量方程
Figure BDA0001538799860000107
是渐近稳定的，系统
Figure BDA0001538799860000108
的所有特征根位于左半平面。
步骤A29：有x(t)＝Uz(t)
，
Figure BDA0001538799860000109
和[zm+1(t),....,zn(t)]T＝0，故可得x(t)＝Uz(t)＝U[z1(t),...,zm(t),0...,0]T。因为矩阵U是非奇异的，故系统
Figure BDA00015387998600001010
是渐近稳定的等价于系统
Figure BDA00015387998600001011
是渐近稳定的同样等价于
Figure BDA00015387998600001012
式中m个标量方程渐近稳定。
作为可选地实施方式，步骤B，具体为：
若条件det(σE-A)≠0成立，其中：E,A∈Rn×n，δ∈C；令ak≠0，k＝1,2,...,m，ak＝0，k＝m+1,m+2,...,n，并且
Figure BDA00015387998600001013
其中uk,vk∈R；定义αk＝λk(Eδ)；当且仅当对于所有的ak≠0，k＝1,2,...,m，有条件
Figure BDA00015387998600001014
成立时，判定系统
Figure BDA00015387998600001015
是渐近稳定的。
下面更为详细地说明本发明步骤B即奇异时滞系统稳定性的判别的简单方法：
步骤B1：系统
Figure BDA0001538799860000111
的稳定性由m个标量微分方程
Figure BDA0001538799860000112
决定，其中αk＝λk(Eδ),βk＝δλk(Eδ)-1。
因为ak≠0，k＝1,2,...,m，系统
Figure BDA0001538799860000113
可以改写成
Figure BDA0001538799860000114
即为
Figure BDA0001538799860000115
步骤B2：矩阵A的特征值可以表示成λk(A)＝ak+ibk，其中ak,bk∈R,对于k＝1,2,...,n，当且仅当对于所有特征值λk(A),条件
Figure BDA0001538799860000116
成立时，时滞系统
Figure BDA0001538799860000117
是渐近稳定的。
步骤B3：
Figure BDA0001538799860000118
是一个时滞微分方程。E,A∈Rn×n。假设条件det(σE-A)≠0成立，令ak≠0，k＝1,2,...,m，ak＝0，k＝m+1,m+2,...,n，并且
Figure BDA0001538799860000119
其中uk,vk∈R。
当且仅当对于所有的ak≠0，k＝1,2,...,m，条件
Figure BDA00015387998600001110
成立时，系统
Figure BDA00015387998600001111
是渐近稳定的。
F的第j个特征值由λj(F)表示，FH表示F的共轭转置，符合supt∈[-τ,0]表示f(t)在t∈[-τ,0]上的上确界。
从证明结论可以看出，本发明给出了一类奇异时滞系统的简单的稳定性的判别方法，因为给出的判别方法是充分必要的，因此系统
Figure BDA00015387998600001112
的稳定性可以被完全描述。
如图2所示，本发明实施例提供的一类奇异时滞系统的渐近稳定性维持方法，包括以下步骤：
根据本发明任一技术方案提供的一类奇异时滞系统的渐近稳定性判别方法判别该类奇异时滞系统的渐近稳定性；
利用维持该类奇异时滞依赖系统稳定性的条件，以提高该类奇异时滞系统的渐近稳定性。
上述方法应用了本发明提供的一类奇异时滞系统的渐近稳定性判别方法提高了该类奇异时滞系统的渐近稳定性。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制；尽管参照较佳实施例对本发明进行了详细的说明，所属领域的普通技术人员应当理解：依然可以对本发明的具体实施方式进行修改或者对部分技术特征进行等同替换；而不脱离本发明技术方案的精神，其均应涵盖在本发明请求保护的技术方案范围当中。

Claims (8)

1.一类奇异时滞系统的渐近稳定性判别方法，其特征在于，包括以下步骤:
步骤A、利用数学分解方法对一类奇异时滞系统的时滞稳定性进行分析，以得出维持该类奇异时滞依赖系统稳定性的条件；
步骤B、依据所述维持一类奇异时滞依赖系统稳定性的条件，判别所述该类奇异时滞依赖系统的稳定性；
其中，所述步骤A，包括以下步骤：
步骤A1、得到非奇异时滞系统的稳定的条件；
步骤A2、将非奇异时滞系统的稳定的条件与Schur分解方法相结合来分析一类奇异时滞系统的时滞稳定性，以得出该类维持奇异时滞依赖系统稳定性的条件；
所述步骤A2，包括以下步骤:
步骤A21：假设条件det(σE-A)≠0成立，矩阵(δE-A)-1存在，系统
Figure FDA00025859688500000110
其中E,A∈Rn×n，δ∈C，E为奇异矩阵；两边都左乘(δE-A)-1，于是得到
Figure FDA0002585968850000011
式子；
步骤A22：因为矩阵(δE-A)-1存在，令
Figure FDA0002585968850000012
得出
Figure FDA0002585968850000013
并得出
Figure FDA0002585968850000014
因此将式子
Figure FDA0002585968850000015
改写成
Figure FDA0002585968850000016
步骤A23：因为
Figure FDA0002585968850000017
为方阵，存在酉矩阵U使得
Figure FDA0002585968850000018
成立；其中T是一个上三角矩阵，它的对角元素是
Figure FDA0002585968850000019
的特征值；选取适当的U阵使得T前m个对角元素都是非零的特征值α1,α2,…,αm，而剩余的对角线元素是零特征值αm+1,αm+2,…,αn；
步骤A24：令x(t)＝Uz(t)，则
Figure FDA0002585968850000021
写成
Figure FDA0002585968850000022
的形式；
步骤A25：令z(t)＝[z1(t),....,zm(t),zm+1(t),....zn(t)]T，系统
Figure FDA0002585968850000023
的最后一个标量方程仅包含一个变量zn(t)，第n个方程是
Figure FDA0002585968850000024
并得一致性条件zn(t)≡0,t≥-τ；
步骤A26：令zn(t)≡0,t≥-τ，在第n-1个方程中有
Figure FDA0002585968850000025
显然zn-1(t)≡0,t≥-τ是成立的；现在已证得标量元素zn(t)≡zn-1(t)≡0,t≥-τ，同理也证得zk(t)≡0,t≥-τ,k＝m+1,m+2,...,n.
步骤A27：在系统
Figure FDA0002585968850000026
中，对于k＝m+1,m+2,....n；令zk(t)＝0得到一个m维的系统
Figure FDA0002585968850000027
其中R的对角元素α1,α2,…,αm都不为零，并且有
Figure FDA0002585968850000028
亦有[zm+1(t),....,zn(t)]T＝0；
步骤A28：因为矩阵R是非奇异的，系统
Figure FDA0002585968850000029
化为一个时滞系统，当且仅当它的特征方程
Figure FDA00025859688500000210
的所有特征根都在左半平面上时，该系统是渐近稳定的；
当且仅当m个标量方程
Figure FDA00025859688500000211
是渐近稳定的，系统
Figure FDA00025859688500000212
的所有特征根位于左半平面；
步骤A29：有式子x(t)＝Uz(t)＝U[z1(t),...,zm(t),0...,0]T并且矩阵U是非奇异的，系统
Figure FDA00025859688500000213
是渐近稳定的等价于系统
Figure FDA00025859688500000214
是渐近稳定的等价于
Figure FDA00025859688500000215
式中m个标量方程渐近稳定。
2.根据权利要求1所述的一类奇异时滞系统的渐近稳定性判别方法，其特征在于，所述数学分解方法为Schur分解方法。
3.根据权利要求1所述的一类奇异时滞系统的渐近稳定性判别方法，所述条件为充分必要条件。
4.根据权利要求1所述的一类奇异时滞系统的渐近稳定性判别方法，其特征在于，所述奇异时滞系统为线性时滞系统。
5.根据权利要求1所述的一类奇异时滞系统的渐近稳定性判别方法，其特征在于，所述步骤B，具体为：
若条件det(σE-A)≠0成立，其中：E,A∈Rn×n，δ∈C；令ak≠0，k＝1,2,...,m，ak＝0，k＝m+1,m+2,...,n，并且
Figure FDA0002585968850000031
其中uk,vk∈R；定义αk＝λk(Eδ)；当且仅当对于所有的ak≠0，k＝1,2,...,m，有条件
Figure FDA0002585968850000032
成立时，判定系统
Figure FDA0002585968850000033
是渐近稳定的。
6.根据权利要求5所述的一类奇异时滞系统的渐近稳定性判别方法，其特征在于，所述步骤B，具体包括：
步骤B1：系统
Figure FDA0002585968850000034
的稳定性由m个标量微分方程
Figure FDA0002585968850000035
决定，其中αk＝λk(Eδ),βk＝δλk(Eδ)-1；
因为ak≠0，k＝1,2,...,m，系统
Figure FDA0002585968850000036
改写成
Figure FDA0002585968850000037
即为
Figure FDA0002585968850000038
步骤B2：矩阵A的特征值表示成λk(A)＝ak+ibk，其中ak,bk∈R,对于k＝1,2,...,n，当且仅当对于所有特征值λk(A),条件
Figure FDA0002585968850000039
时滞系统
Figure FDA00025859688500000310
是渐近稳定的；
步骤B3：
Figure FDA00025859688500000311
是一个时滞微分方程；E,A∈Rn×n；假设条件det(σE-A)≠0成立，令ak≠0，k＝1,2,...,m，ak＝0，k＝m+1,m+2,...,n，并且
Figure FDA0002585968850000041
其中uk,vk∈R；
当且仅当对于所有的ak≠0，k＝1,2,...,m，条件
Figure FDA0002585968850000042
成立时，系统
Figure FDA0002585968850000043
是渐近稳定的。
7.根据权利要求1所述的一类奇异时滞系统的渐近稳定性判别方法，其特征在于，所述奇异时滞系统为奇异时滞依赖系统。
8.一类奇异时滞系统的渐近稳定性维持方法，其特征在于，包括以下步骤：
根据权利要求1-7任一所述的一类奇异时滞系统的渐近稳定性判别方法判别该类奇异时滞系统的渐近稳定性；
利用所述维持一类奇异时滞依赖系统稳定性的条件，以提高该类奇异时滞系统的渐近稳定性。
CN201810006540.9A 2018-01-04 2018-01-04 一类奇异时滞系统的渐近稳定性判别方法、维持方法 Active CN108241297B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810006540.9A CN108241297B (zh) 2018-01-04 2018-01-04 一类奇异时滞系统的渐近稳定性判别方法、维持方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810006540.9A CN108241297B (zh) 2018-01-04 2018-01-04 一类奇异时滞系统的渐近稳定性判别方法、维持方法

Publications (2)

Publication Number Publication Date
CN108241297A CN108241297A (zh) 2018-07-03
CN108241297B true CN108241297B (zh) 2020-11-13

Family

ID=62699336

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810006540.9A Active CN108241297B (zh) 2018-01-04 2018-01-04 一类奇异时滞系统的渐近稳定性判别方法、维持方法

Country Status (1)

Country Link
CN (1) CN108241297B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101645601A (zh) * 2009-08-27 2010-02-10 天津大学 电力系统时滞依赖型鲁棒稳定的判别方法
KR101090341B1 (ko) * 2010-08-03 2011-12-07 인하대학교 산학협력단 ＩＥＥＥ．８０２．１５．４ａ 첩 스프레드 스펙트럼 기반 저 복잡도 거리측정 알고리즘
CN102801158A (zh) * 2012-07-31 2012-11-28 山东电力集团公司电力经济技术研究院 基于Pade近似的时滞电力系统特征值计算与稳定性判别方法
CN103576556A (zh) * 2013-11-18 2014-02-12 浙江理工大学 一种稳定奇异时滞电路系统的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106099921B (zh) * 2016-07-21 2018-09-28 天津大学 一种电力系统时滞稳定裕度快速求解方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101645601A (zh) * 2009-08-27 2010-02-10 天津大学 电力系统时滞依赖型鲁棒稳定的判别方法
KR101090341B1 (ko) * 2010-08-03 2011-12-07 인하대학교 산학협력단 ＩＥＥＥ．８０２．１５．４ａ 첩 스프레드 스펙트럼 기반 저 복잡도 거리측정 알고리즘
CN102801158A (zh) * 2012-07-31 2012-11-28 山东电力集团公司电力经济技术研究院 基于Pade近似的时滞电力系统特征值计算与稳定性判别方法
CN103576556A (zh) * 2013-11-18 2014-02-12 浙江理工大学 一种稳定奇异时滞电路系统的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
时滞奇异系统的稳定性分析;叶菲;《中国优秀硕士学位论文全文数据库-信息科技辑》;20140915(第09期);第6、12、17、18、21、52页 *

Also Published As

Publication number Publication date
CN108241297A (zh) 2018-07-03

Similar Documents

Publication Publication Date Title
Deng et al. Output-feedback stabilization of stochastic nonlinear systems driven by noise of unknown covariance
Zhang et al. A vehicle ABS adaptive sliding-mode control algorithm based on the vehicle velocity estimation and tyre/road friction coefficient estimations
Yu et al. Recursive identification for Hammerstein–Wiener systems with dead-zone input nonlinearity
Grima Linear-noise approximation and the chemical master equation agree up to second-order moments for a class of chemical systems
Wu et al. Incremental modeling of a new high-order polynomial surrogate model
CN110286595B (zh) 一类受饱和非线性输入影响的分数阶系统自适应控制方法
CN111159636B (zh) 一种柔性多体系统动力学半解析灵敏度分析方法
Li et al. Observer-based reliable exponential stabilization and H∞ control for switched systems with faulty actuators: an average dwell time approach
Liu et al. Output feedback control for stochastic nonholonomic systems with growth rate restriction
CN108196446A (zh) 模型未知的双电机负载的动态规划最优控制方法
Li et al. Design of exponential state estimator for neural networks with distributed delays
CN108241297B (zh) 一类奇异时滞系统的渐近稳定性判别方法、维持方法
CN115773569A (zh) 基于自抗扰解耦的海洋平台通风系统风量控制方法
Tang et al. The error bounds of combined forecasting
CN109995278B (zh) 一种考虑输入受限的电机伺服系统自调节控制方法
CN104730925A (zh) 一种输入限幅pi控制方法
CN111781830A (zh) 基于估计状态扰动观测的混合阶机电系统协同控制方法
Hu et al. Stability of Runge-Kutta methods for delay differential systems with multiple delays
Bruggeman et al. Time-dependent hierarchical regulation analysis: deciphering cellular adaptation
Shafai et al. Robust stability and stabilization of uncertain delay systems
EP3168578B1 (en) Estimating system parameters from sensor measurements
Sun et al. Large time behaviors of solutions to the unipolar hydrodynamic model of semiconductors with physical boundary effect
CN105094114B (zh) 一种确定lpv控制器的稳定性的方法
Teng Permanence and stability in non-autonomous logistic systems with infinite delay
CN106444383A (zh) 前件不匹配的tsfmb时滞系统控制方法

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