系统稳定性判定分析(一)---- 常系数线性系统内部稳定性

2024-05-01 15:20

本文主要是介绍系统稳定性判定分析(一)---- 常系数线性系统内部稳定性,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

文章目录

  • 线性系统稳定性分析
    • 线性定常系统(线性时不变/自治系统)稳定性分析
    • 线性定常系统内部稳定性示例
  • 参考文献

  从上一节 系统稳定性的介绍中可以得知,分析系统内部稳定(Lyapunov意义下稳定)时可不考虑系统的输出情况,可直接通过系统的状态方程分析系统的稳定性。系统的状态方程根据构建形式的不同,可分为线性系统与非线性系统。如下基于由简入繁的原则,首先整理线性系统的稳定性分析过程。

本文内容主要基于Antsaklis, P. J., & Michel, A. N. (1997). Linear systems (Vol. 8). New York: McGraw-Hill. 一书。

线性系统稳定性分析

   ~~~~     线性系统根据状态方程中是否显含时间变量 t t t(状态变量本身为时间 t t t的函数,因此状态方程一定与时间 t t t有关。若状态方程中显示的包含时间 t t t的表达式,则表示状态方程显含时间变量 t t t。否则,状态方程不显含时间变量 t t t),可将线性系统划分为线性自治系统(或线性时不变系统) x ˙ = f ( x ( t ) ) , (1) \dot{\pmb{x}}=f(\pmb{x}(t)),\tag{1} x˙=f(x(t)),(1)与线性非自治系统(线性时变系统) x ˙ = f ( t , x ( t ) ) , (2) \dot{\pmb{x}}=f(t,\pmb{x}(t)),\tag{2} x˙=f(t,x(t)),(2)通常为表述简单,系统 (1) 与 (2) 可直接省略括号内的时间变量 t t t,表示为 x ˙ = f ( x ) \dot{\pmb{x}}=f(\pmb{x}) x˙=f(x) x ˙ = f ( t , x ) \dot{\pmb{x}}=f(t,\pmb{x}) x˙=f(t,x)
   ~~~~     如下首先分析线性自治系统的稳定性。

线性定常系统(线性时不变/自治系统)稳定性分析

   ~~~~     考虑如下线性时不变系统 x ˙ = A x . (3) \dot{\pmb{x}}=A\pmb{x}.\tag{3} x˙=Ax.(3)定义 x e \pmb{x}_e xe 为系统(3)平衡点,则系统(3)李雅普诺夫意义下的稳定性可通过如下定理进行判定:

定理1. 系统(3)的平衡点 x e \pmb{x}_e xe稳定的,当且仅当其系数矩阵 A A A (或者系统(3)的雅可比矩阵) 的所有特征值的实部小于等于0,且每个实部为零的特征值都有一个相关的 1 阶 Jordan 块。系统(3)的平衡点 x e \pmb{x}_e xe渐近稳定的,当且仅当其系数矩阵 A A A (或者系统(3)的雅可比矩阵) 的所有特征值的实部小于0
相反地,系统(3)的平衡点 x e \pmb{x}_e xe不稳定的,当且仅当其系数矩阵 A A A (或者其雅可比矩阵) 至少存在一个实部大于0的特征值,或者其实部为零的特征值所对应的Jordan块的阶数大于1阶。

补充: 有关特征值,Jordan块的介绍可参见矩阵分析:特征值,相似度对角化,Jordan标准形。
注: 通常我们所说的一个系统的收敛是与系统平衡点的渐近稳定有关,关于两者有如下定义:

定义1. 系统(3)的平衡点 x e \pmb{x}_e xe 是渐近稳定的当且仅当其(1) 是稳定的,(2) 当 t → ∞ t\to \infty t 时,系统(3)的解趋近于平衡点 x e \pmb{x}_e xe

线性定常系统内部稳定性示例

【示例一】分析如下系统在李雅普诺夫意义下的稳定性:
[ x ˙ 1 x ˙ 2 ] = [ 0 1 − 1 0 ] [ x 1 x 2 ] . (4) \begin{bmatrix} \dot x_1\\ \dot x_2 \end{bmatrix}=\begin{bmatrix} 0&1\\ -1&0 \end{bmatrix}\begin{bmatrix} x_1\\ x_2\end{bmatrix}.\tag{4} [x˙1x˙2]=[0110][x1x2].(4)通过计算可以得知系统(4)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T,其特征值为 λ 1 , 2 = ± i \lambda_{1,2} = \pm i λ1,2=±i,其Jordan块为 [ i 0 0 − i ] \begin{bmatrix} i&0\\ 0&-i \end{bmatrix} [i00i]。根据定理1,系统(4)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T是稳定的。系统(4)在初值为(0,0)条件下的运行轨迹如下所示:
在这里插入图片描述
系统(4)在初值为(0,1)条件下的运行轨迹如下所示:
在这里插入图片描述
结合以上两图可知,系统(4)的平衡点是稳定的,整体运行轨迹在平衡点的一定范围内波动。

【示例二】分析如下系统在李雅普诺夫意义下的稳定性:
[ x ˙ 1 x ˙ 2 ] = [ 0 1 0 0 ] [ x 1 x 2 ] . (5) \begin{bmatrix} \dot x_1\\ \dot x_2 \end{bmatrix}=\begin{bmatrix} 0&1\\ 0&0 \end{bmatrix}\begin{bmatrix} x_1\\ x_2\end{bmatrix}.\tag{5} [x˙1x˙2]=[0010][x1x2].(5)通过计算可以得知系统(5)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T,其特征值为 λ 1 , 2 = 0 \lambda_{1,2} = 0 λ1,2=0,其Jordan块为 [ 0 1 0 0 ] \begin{bmatrix} 0&1\\ 0&0 \end{bmatrix} [0010],特征值0所对应的Jordan块的阶数为2。根据定理1,系统(5)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T是不稳定的。系统(5)在初值为(0,1)条件下的运行轨迹如下所示:
在这里插入图片描述
【示例三】分析如下系统在李雅普诺夫意义下的稳定性:
[ x ˙ 1 x ˙ 2 ] = [ 2.8 9.6 9.6 − 2.8 ] [ x 1 x 2 ] . (6) \begin{bmatrix} \dot x_1\\ \dot x_2 \end{bmatrix}=\begin{bmatrix} 2.8&9.6\\ 9.6&-2.8 \end{bmatrix}\begin{bmatrix} x_1\\ x_2\end{bmatrix}.\tag{6} [x˙1x˙2]=[2.89.69.62.8][x1x2].(6)通过计算可以得知系统(6)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T,其特征值为 λ 1 , 2 = ± 10 \lambda_{1,2} = \pm10 λ1,2=±10,其Jordan块为 [ 10 0 0 − 10 ] \begin{bmatrix} 10&0\\ 0&-10 \end{bmatrix} [100010],存在一个实部大于0的特征值。根据定理1,系统(6)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T是不稳定的。系统(6)在初值为(0,1)条件下的运行轨迹如下所示:在这里插入图片描述
【示例四】分析如下系统在李雅普诺夫意义下的稳定性:
[ x ˙ 1 x ˙ 2 ] = [ − 1 0 − 1 − 2 ] [ x 1 x 2 ] . (7) \begin{bmatrix} \dot x_1\\ \dot x_2 \end{bmatrix}=\begin{bmatrix} -1&0\\ -1&-2 \end{bmatrix}\begin{bmatrix} x_1\\ x_2\end{bmatrix}.\tag{7} [x˙1x˙2]=[1102][x1x2].(7)通过计算可以得知系统(7)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T,其特征值为 λ 1 = − 1 \lambda_1= -1 λ1=1 λ 2 = − 2 \lambda_2= -2 λ2=2其Jordan块为 [ − 1 0 0 − 2 ] \begin{bmatrix} -1&0\\ 0&-2 \end{bmatrix} [1002],存在一个实部大于0的特征值。根据定理1,系统(7)的平衡点 x e = [ x 1 e , x 2 e ] T = [ 0 , 0 ] T \pmb{x}_e=[x_{1e}, x_{2e}]^T=[0, 0]^T xe=[x1e,x2e]T=[0,0]T是渐近稳定的。系统(7)在初值为(0,1)条件下的运行轨迹如下所示:
在这里插入图片描述

注: 上述四个示例的代码如下,需适当修改函数中的状态方程表达式,以及主函数中的时间变量的区间范围。

function dxdt = vdp1(t,x)%dxdt = [x(2);0];
%  dxdt = [x(2);-x(1)];
% dxdt = [-x(1);-x(1)-2*x(2)]; %dxdt = [x(2);-2*sin(x(1))];dxdt = [2.8*x(1)+9.6*x(2);9.6*x(1)-2.8*x(2)];
end
%[t,y] = ode45(odefun,tspan,y0)
[t,x] = ode45(@vdp1,[0 20],[0,1]);
plot(t,x)
xlabel('Time t');
ylabel('Solution x');
legend('x_1','x_2')

参考文献

[1] Antsaklis, P. J., & Michel, A. N. (1997). Linear systems (Vol. 8). New York: McGraw-Hill.
[2] 线性系统稳定性一般定理、齐次线性系统稳定性
[3] 矩阵分析:特征值,相似度对角化,Jordan标准形

这篇关于系统稳定性判定分析(一)---- 常系数线性系统内部稳定性的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/952011

相关文章

找不到Anaconda prompt终端的原因分析及解决方案

《找不到Anacondaprompt终端的原因分析及解决方案》因为anaconda还没有初始化,在安装anaconda的过程中,有一行是否要添加anaconda到菜单目录中,由于没有勾选,导致没有菜... 目录问题原因问http://www.chinasem.cn题解决安装了 Anaconda 却找不到 An

Spring定时任务只执行一次的原因分析与解决方案

《Spring定时任务只执行一次的原因分析与解决方案》在使用Spring的@Scheduled定时任务时,你是否遇到过任务只执行一次,后续不再触发的情况?这种情况可能由多种原因导致,如未启用调度、线程... 目录1. 问题背景2. Spring定时任务的基本用法3. 为什么定时任务只执行一次?3.1 未启用

Linux系统之主机网络配置方式

《Linux系统之主机网络配置方式》:本文主要介绍Linux系统之主机网络配置方式,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、查看主机的网络参数1、查看主机名2、查看IP地址3、查看网关4、查看DNS二、配置网卡1、修改网卡配置文件2、nmcli工具【通用

Linux系统之dns域名解析全过程

《Linux系统之dns域名解析全过程》:本文主要介绍Linux系统之dns域名解析全过程,具有很好的参考价值,希望对大家有所帮助,如有错误或未考虑完全的地方,望不吝赐教... 目录一、dns域名解析介绍1、DNS核心概念1.1 区域 zone1.2 记录 record二、DNS服务的配置1、正向解析的配置

C++ 各种map特点对比分析

《C++各种map特点对比分析》文章比较了C++中不同类型的map(如std::map,std::unordered_map,std::multimap,std::unordered_multima... 目录特点比较C++ 示例代码 ​​​​​​代码解释特点比较1. std::map底层实现:基于红黑

Linux系统中配置静态IP地址的详细步骤

《Linux系统中配置静态IP地址的详细步骤》本文详细介绍了在Linux系统中配置静态IP地址的五个步骤,包括打开终端、编辑网络配置文件、配置IP地址、保存并重启网络服务,这对于系统管理员和新手都极具... 目录步骤一:打开终端步骤二:编辑网络配置文件步骤三:配置静态IP地址步骤四:保存并关闭文件步骤五:重

Spring、Spring Boot、Spring Cloud 的区别与联系分析

《Spring、SpringBoot、SpringCloud的区别与联系分析》Spring、SpringBoot和SpringCloud是Java开发中常用的框架,分别针对企业级应用开发、快速开... 目录1. Spring 框架2. Spring Boot3. Spring Cloud总结1. Sprin

Spring 中 BeanFactoryPostProcessor 的作用和示例源码分析

《Spring中BeanFactoryPostProcessor的作用和示例源码分析》Spring的BeanFactoryPostProcessor是容器初始化的扩展接口,允许在Bean实例化前... 目录一、概览1. 核心定位2. 核心功能详解3. 关键特性二、Spring 内置的 BeanFactory

Windows系统下如何查找JDK的安装路径

《Windows系统下如何查找JDK的安装路径》:本文主要介绍Windows系统下如何查找JDK的安装路径,文中介绍了三种方法,分别是通过命令行检查、使用verbose选项查找jre目录、以及查看... 目录一、确认是否安装了JDK二、查找路径三、另外一种方式如果很久之前安装了JDK,或者在别人的电脑上,想

MyBatis-Plus中Service接口的lambdaUpdate用法及实例分析

《MyBatis-Plus中Service接口的lambdaUpdate用法及实例分析》本文将详细讲解MyBatis-Plus中的lambdaUpdate用法,并提供丰富的案例来帮助读者更好地理解和应... 目录深入探索MyBATis-Plus中Service接口的lambdaUpdate用法及示例案例背景