Lyapunov函数的一些构造方法
-Lyapunov指数的计算方法
【总结】Lyapunov指数的计算方法非线性理论近期为了把计算LE的一些问题弄清楚，看了有7～9本书！下面以吕金虎《混沌时间序列分析及其应用》、马军海《复杂非线性系统的重构技术》为主线，把目前已有的LE计算方法做一个汇总！1. 关于连续系统Lyapunov指数的计算方法连续系统LE的计算方法主要有定义方法、Jacobian方法、QR分解方法、奇异值分解方法，或者通过求解系统的微分方程，得到微分方程解的时间序列，然后利用时间序列（即离散系统）的LE求解方法来计算得到。
关于连续系统LE的计算，主要以定义方法、Jacobian方法做主要介绍内容。
（1）定义法定义法求解Lyapunov指数.JPG关于定义法求解的程序，和matlab板块的“连续系统LE求解程序”差不多。
以Rossler系统为例Rossler系统微分方程定义程序function dX = Rossler_ly(t,X)% Rossler吸引子，用来计算Lyapunov指数% a=0.15,b=0.20,c=10.0% dx/dt = -y-z,% dy/dt = x+ay,% dz/dt = b+z(x-c),a = 0.15;b = 0.20;c = 10.0;x=X(1); y=X(2); z=X(3);% Y的三个列向量为相互正交的单位向量Y = [X(4), X(7), X(10);X(5), X(8), X(11);X(6), X(9), X(12)];% 输出向量的初始化，必不可少dX = zeros(12,1);% Rossler吸引子dX(1) = -y-z;dX(2) = x+a*y;dX(3) = b+z*(x-c);% Rossler吸引子的Jacobi矩阵Jaco = [0 -1 -1;1 a 0;z 0 x-c];dX(4:12) = Jaco*Y;求解LE代码：% 计算Rossler吸引子的Lyapunov指数clear;yinit = [1,1,1];orthmatrix = [1 0 0;0 1 0;0 0 1];a = 0.15;b = 0.20;c = 10.0;y = zeros(12,1);% 初始化输入y(1:3) = yinit;y(4:12) = orthmatrix;tstart = 0; % 时间初始值tstep = 1e-3; % 时间步长wholetimes = 1e5; % 总的循环次数steps = 10; % 每次演化的步数iteratetimes = wholetimes/steps; % 演化的次数mod = zeros(3,1);lp = zeros(3,1);% 初始化三个Lyapunov指数Lyapunov1 = zeros(iteratetimes,1); Lyapunov2 = zeros(iteratetimes,1); Lyapunov3 = zeros(iteratetimes,1);for i=1:iteratetimestspan = tstart:tstep:(tstart + tstep*steps); [T,Y] = ode45('Rossler_ly', tspan, y);% 取积分得到的最后一个时刻的值y = Y(size(Y,1),:);% 重新定义起始时刻tstart = tstart + tstep*steps;y0 = [y(4) y(7) y(10);y(5) y(8) y(11);y(6) y(9) y(12)];%正交化y0 = ThreeGS(y0);% 取三个向量的模mod(1) = sqrt(y0(:,1)'*y0(:,1));mod(2) = sqrt(y0(:,2)'*y0(:,2));mod(3) = sqrt(y0(:,3)'*y0(:,3));y0(:,1) = y0(:,1)/mod(1);y0(:,2) = y0(:,2)/mod(2);y0(:,3) = y0(:,3)/mod(3);lp = lp+log(abs(mod));%三个Lyapunov指数Lyapunov1(i) = lp(1)/(tstart);Lyapunov2(i) = lp(2)/(tstart);Lyapunov3(i) = lp(3)/(tstart);y(4:12) = y0';end% 作Lyapunov指数谱图i = 1:iteratetimes;plot(i,Lyapunov1,i,Lyapunov2,i,Lyapunov3)程序中用到的ThreeGS程序如下：%G-S正交化function A = ThreeGS(V) % V 为3*3向量v1 = V(:,1);v2 = V(:,2);v3 = V(:,3);a1 = zeros(3,1);a2 = zeros(3,1);a3 = zeros(3,1);a1 = v1;a2 = v2-((a1'*v2)/(a1'*a1))*a1;a3 = v3-((a1'*v3)/(a1'*a1))*a1-((a2'*v3)/(a2'*a2))*a2;A = [a1,a2,a3];计算得到的Rossler系统的LE为———— 0.063231 0.092635 -9.8924Wolf文章中计算得到的Rossler系统的LE为————0.09 0 -9.77需要注意的是——定义法求解的精度有限，对有些系统的计算往往出现计果和理论值有偏差的现象。
lyapunov函数
lyapunov函数Lyapunov函数是一种用于研究系统稳定性的重要工具，它可以用来检验系统耐受外界干扰的能力，以及系统发生振荡现象的可能性。
它是由俄罗斯数学家安德烈利亚普诺夫（Andrey Lyapunov）最早提出的，是一种重要的动态系统的稳定性理论。
Lyapunov函数也称为质能函数（或者拉普拉斯函数），它是一个定义在特定空间中的实值函数，它能够浓缩动态系统状态并同时反映系统的稳定性。
简而言之，Lyapunov函数有助于发现动态系统中某个状态点是否是稳定的。
Lyapunov函数的定义为：对于一个系统，它可以根据任意一个参数来评估系统状态点的稳定性。
当系统从状态点A移动到状态点B 时，系统的Lyapunov函数能够帮助我们了解更多的系统的行为：第一，该函数能够衡量不同状态点之间的影响，即该函数能够测量系统从A点变为B点的过程中，A点对B点的影响有多大。
第二，它能够表明系统能否在未来某个时刻稳定地维持本身的状态点。
由于Lyapunov函数可以衡量系统状态点的稳定性，因此可以用来实施控制策略，来防止系统振荡。
Lyapunov函数可以作为一个智能控制系统中的一个重要分量，它能够有效地检测和对应外部的环境因素，进而把外部的环境因素转换为控制指令，以便调节系统的状态。
Lyapunov函数在动态系统建模和分析方面，具有无可比拟的效用。
除了上述提到的用于检验系统稳定性以外，它还可以用来检查系统控制功能、建立系统模型和探索各种可能性，以此探究系统的行为特征，它对数学原理研究和应用研究都是非常重要的。
Lyapunov函数的研究仍在不断发展，越来越多的研究者将Lyapunov函数与其他技术融合在一起，以便更好地理解系统的行为，提高控制策略的可靠性，更好地探索异常情况的发生可能性，以及更多的分析细节。
总之，Lyapunov函数是一种重要的动态系统理论，它能够帮助我们检验系统的稳定性，为智能控制系统提供重要参考，在模型建立和控制策略制定等方面具有重要意义，而且它也在不断发展，以适应不断变化的环境和情况。
常微分方程中李雅普诺夫函数构造方法
武夷学院学报JOURNAL OF WUYI UNIVERSITY第40卷第3期2021年3月Vol.40 No.3Mar. 2021常微分方程中李雅普诺夫函数构造方法陆求赐】，张宋传2,王学彬2(1.武夷学院人文与教师教育学院，福建武夷山354300； 2.武夷学院数计学院，福建武夷山354300)扌商 要：李雅普诺夫直接法,就是在不求方程组解的情况下，构造一个李雅普诺夫函数,通过微分方程组所计算岀来的全导数的符号性质，来判断微分方程组零解的稳定性，但至今仍没有一种统一的方法来构造李雅普诺夫函数。
通过 对一个含参数例子的分析,介绍几种常见且适用的李雅普诺夫函数的构造方法:首次积分法、能量函数法、分离变量法、待定系数法和二次型矩阵法等。
关键词：常微分方程;李雅普诺夫函数;零解的稳定性;构造方法中图分类号：O175.26 文献标识码:A 考虑如下自治的微分方程组牯(x )(1)这里X =(x 1,x 2，…x ”)表示n ("逸2)维向量,i =1,2，…,n ,假设f (0)=0(0表示向量),且f (x )在某邻域G ：||x ||臆A (A 为正常数)内有连续的偏导数，从而方程组(1)的由其初值条件x (t °)=x 。
所决定的解在邻域G 内存在且唯一°李雅普诺夫第二方法:通过方程组(1)来构造一个特殊的函数V (x ),并假设V (x )关于所有变元的偏导数存在且连续，将V (x )对变量t 求导，并将方程组(1)代入得：d V § dV dx , § 坠卩 f(2)dt =移 dx, dt=移 dx ,)这个式(2)就是函数V (x )沿方程组(1)求得的全导数，现通过式(2)的符号来判定方程组(1)零解的稳定 性,这就是李雅普诺夫第二方法的思路,这个函数收稿日期：2020-07-21基金项目：武夷学院校科研基金项目(XL201408)；福建省教育厅科技项目(JA15512、JAT160519)；福建省自然科学基金(2016J01682)；高级引进人才 科研启动基金(YJ201802)资助°作者简介:陆求赐(1975-)，男,汉族,副教授，主要从事基础数学教学和微分方程的研究°文章编号：1674-2109(2021)03-0016-05V (x )就称为李雅普诺夫函数，简称V 函数［1-2］°李雅普诺夫函数作为常微分方程的重要内容之一,学好它可以为进入常微分方程领域的研究打好基础。
Lyapunov函数的一些构造方法
dL 并且 1 dt
y v 0 x x* , * * , w w* . y v M
OK
Lyapunov 函数的构造
三、Lyapunov函数形式的推广
* y u y u x* x*v* v u v* 其中 F * du * du du. * * x y v u u ky u x
推广 形式：
u w* w* u du,
w
w w*
u w* u w du, * du, w f u
*
w( w* ) du, f u
w
w*
f (u ) f ( w* ) du. u
Lyapunov 函数的构造
四、一类SI模型的全局稳定性
Lyapunov 函数的一些构造方法
李建全
空军工程大学(西安)
E-mail: jianq_li@
Lyapunov 函数的构造
一、引言 二、带有免疫的病毒动力学模型 三、Lyapunov函数形式的推广 四、一类SI模型的全局稳定性
Lyapunov 函数的构造
一、引言
全局渐近稳定性 (GAS)是动力系统研究的重要 问题。 对于高维自治微分系统 ，证明GAS最常用的方法是 Lyapunov - LaSalle定理， 即若自治系统 * x' f x, x Rn
函数 2、 情形f w 1. 定义Lyapunov
*
dL 则 1 dt
* x x p 0 * * x 2 w w * x x M c
2
x* xy*v y v* x v 3 x x* yv* y *v 0,
Liapunov函数的构造
Liapunov 函数的构造摘要：Liapunov 函数是一种判定微分方程零解稳定性的重要方法，所以本文首先介绍了Liapunov 函数以及判断微分方程的稳定性定理，然后着重介绍了Liapunov 函数的几种构造方法，包括常系数线性系统的巴尔巴欣公式、线性类比法.通过这两种构造方法，我们将初步了解Liapunov 函数的构造在判断微分方程零解稳定性中的重要作用. 关键词：Liapunov 函数；零解；稳定性.引言在常微分方程中，稳定性理论研究是很重要的一部分，即研究当时间趋于无穷时，其解的形态将会怎样变化，他在自然科学、工程力学、环境生态、社会经济等方面有着重要的应用。
在本章第一节中介绍了稳定性的相关定义，也介绍了对于可以求得微分方程的解析解时，如何利用定义判断其零解的稳定性。
但是在实际问题中提出的微分方程往往是很复杂的，无法求得其解析解，这就需要从方程本身来判断零解的稳定性Liapunov 直接方法就是求解这一问题的有效途径。
本文先引入Liapunov 函数，即V 函数的定义，以及Liapunov 稳定性的定理，然后介绍几种构造Liapunov 函数方法。
1 Liapunov 稳定性的定理1.1 V 函数设函数(x)V 在nR 中原点的某邻域U 中有定义,(x)V 在U 中连续可微，且满足(0)0V =定义 1.1若除原点外对所有x U ∈均有(x)0((x)0)V V ><，则称(x)V 为正定函数（负定函数）；若除原点外对所有x U ∈均有(x)0((x)0)V V ≥≤，则称(x)V 为半正定函数（半负定函数）；若在U 中原点的任一邻域内(x)V 既可以取正，也可以取负，则称(x)V 为变号函数.例如，222123(x)x x x V =++是3R 中的正定函数，但在4R 中确实半正定函数，而222123(x)x x x V =+-是3R 中的变号函数.一般(x)V 函数的符号判断十分困难，通常把(x)V 在原点展开为Taylor 级数12(x)V (x)V (x)V (x),m m m V ++=+++其中12V (x),V (x),V (x)m m m ++分C 是x 的m 次，1m +，2m +齐次函数，根据(x)V 展开式中的最低次项的系数，通常就可以判断(x)V 在原点邻域内的符号.因为再原点附近其他项都可以视为第一项的高阶无穷小. 1.2 Liapunov 稳定性定理设n 维自治微分方程(t)(x),(0)0dx f f dt==（1.1） 的解为12(t)(x (t),x (t),x (t))Tn x = ，为了研究方程（1.1）零解的稳定性，考察随时间变化时(x(t))V 的变化情况，将(x(t))V 视为t 的复合函数，关于t 求导可得11(x(t))(x)(x)n n k k k k k kdx dV V Vf V dt x dt x ==∂∂==∂∂∑∑ （1.2） 式（1.2）称为函数(x)V 沿着方程（1.1）轨线的全导数.介绍了(x)V 以及其全导数后，接下来简单介绍下Liapunov 稳定性理论的几个定理. 定理1.1若有原点的邻域U 和一个正定（负定）函数(x)V ，使得其全导数(x)V是半负定（半正定），则称系统（1.1）的零解是稳定的；特别地，当(x)V是负定（正定）时，系统（1.1）的零解是渐进稳定的.定理1.2设在原点的邻域U 内有函数(x)V ，它沿着方程（1.1）的轨线的全导数(x)V是正定（负定）的，而(x)V 本身不是半负定（半正定）的，则方程（1.1）的零解是不稳定的. 这两个定理是直接通过构造Liapunov 函数，来判断方程的零解是否稳定的，定理在本章第2节已经详细的证明过，这里不再做证明.2 Liapunov 函数的构造第一节中所讨论的两个定理都是一个函数稳定的充分条件，即存在一个()V x ，和它的全导数满足定理1.1时，系统的零解是稳定的. 满足定理1.2的条件时，系统的零解是不稳定的. 在使用Liapunov 函数判定稳定性时应当注意，当找不到满足稳定性定理的条件的函数()V x 时，并无法断言此系统的零解是不稳定的,并且构造的Liapunov 函数不同时，判断零解是否渐进稳定以及吸引域的大小也会有些差异.再利用Liapunov 方法判断系统零解稳定性时，需要明确满足一定条件的Liapunov 函数是否存在，即当系统的零解有某种稳定性时，满足这个稳定性定理的()V x 是否存在，这就是上述定理1.1和定理1.2的逆命题.是否成立. 2.1 Liapunov 函数的存在性考虑微分方程组d (t,x),(t,0)0xf f dt==（2.1）记(){}0,|,G t x t t x h =≥≤，设(t,)f x 在G 连续，关于x 满足Liapunov 条件.令00(t)(t,t ,)x x φ=是方程组（2.1）满足00(t )x φ=的解.定理2.1 若方程组2.1的零解是稳定的，则有正定函数(t,)V x ，使得其全导数(t,)V x是半负定的.证 首先根据方程组（2.1）的零解构造出正定函数(t,)V x ，在验证(t,)V x是半负定.取02(t t )0(t,)(1)(t ,t,)V x e x φ--=+020,()min(t ,t,)t t x y hW x y φ≥≤≤=其中0(t ,t,)x φ表示方程组（2.1）在t 时刻过x 的解在0t 时刻的位置坐标. 显然有02(t t )0(t,0)(1)(t ,t,0)0V e φ--=+=因为方程组（2.1）的零解是稳定的，所以对任意的0ε>，有0δ>，使得当x δ<时，有00(t,t ,),x t t φε<≥（2.2）于是当x h ε≤≤时，对所有的0t t ≥有0(t ,t,)0x φδ≥>.否则就有0x ，0x h ε≤≤，以及10t t ≥，使得010(t ,t ,)x φδ<.由式（2.2）得1001010(t ,t ,(t ,t ,x )),t t φφε<≥又因为1001001000(t ,t ,(t ,t ,x ))(t ,t ,x )x ,x φφφε==≥与式（2.2）矛盾. 由此得(x)W 是正定函数，显然有2010(t,x)(t ,t ,x )(x)V W φ≥≥即(t,x)V 是正定函数，另一方面002(t t )000002(t t )(t,(t,t ,x ))(1)(t ,t,(t,t ,x ))(1)V e ex φφφ----=+=+所以有02(t t )000(t,(t,t ,x ))0dV e x dtφ--=-<即(t,x)V是半负定的.例1研究下述微分方程组零解的稳定性.1221dx x dtdx xdt⎧=-⎪⎪⎨⎪=⎪⎩ 满足初值问题10102020(t ),(t )x x x x ==.解容易解得满足上述初值问题的解为101020100200201020100200(t,t ,,)cos(t t )sin(t t ),(t,t ,,)sin(t t )cos(t t ).x x x x x x x x φφ=---⎧⎨=-+-⎩ 所以有定理2.1，可以取函数0(t t )221210122012(t,,)(1e )((t ,t,,)(t ,t,,))V x x x x x x φφ--=++0(t t )2212(1e )()x x --=++显然0(t t )2222121212(t,,)(1e )()V x x x x x x --=++≥+即通过这种方法构造的12(t,,)V x x 是正定函数.并求得其全导数为：0(t t )221212(t,,)e ()V x x x x --=-+即12(t,,)V x x 正定，12(t,,)V x x半负定，所以上述方程组的零解是稳定的. 2.2 常系数线性系统的巴尔巴欣公式对于常系数线性微分方程组d ,(),,n ij n n xAx A a x R dt⨯==∈（2.3） 这里可以假设该线性自治系统的Liapunov 函数为一个二次型，不妨设其为()TV x x Bx =，则他沿系统（2.3）的全导数为()()(A B BA).TT T d V x x Bx x x dt==+显然其全导数依然是一个二次型.这样，可以从()V x的二次型出发，利用A B BA C T+=进而来确定二次型()V x ，再根据()V x 和()V x的符号来判断方程组（2.3）零解的稳定性.下面就利用这一思想介绍二阶方程组的巴尔巴欣公式.对于二维微分方程组11111222211222dx a x a x dtdx a x a xdt⎧=+⎪⎪⎨⎪=+⎪⎩（2.4）可以给出一个二次型221111212222()2W x w x w x x w x =++进而来求这样的二次型221111212222()2V x v x v x x v x =++使得()V x 沿着方程组（2.4）解曲线的全导数满足(2.4)()2()dV x W x dt =这里的()W x 是选取的，因此111222,,w w w 可视为已知量，利用()W x ，()V x 的定义和上述全导数的等式可以求得未知量111222,,v v v 满足下述方程组：111121121112111122122122121212222222()a v a v w a v a a v a v w a v a v w +=⎧⎪+++=⎨⎪+=⎩当112211221221()()0a a a a a a ∆=+-≠时，可以求解出111222,,v v v ，即可得出()V x2211221111211212112221221222201()20x x x x w a a V x w a a a a w a a =-∆+通过这个过程，可知构造出了一个Liapunov 函数()V x ，并求出其沿着自治系统（2.4）解曲线的全导数()2()V x W x =，从而可以根据()V x 和()W x 的符号来判断系统（2.4）零解的稳定性.例2 讨论微分方程组1121124426dx x x dtdx x xdt⎧=-+⎪⎪⎨⎪=-⎪⎩ 解取221212(,)()W x x x x =-+，利用上述的巴尔巴欣公式可得2211220214201()16004102146x x x x V x --=---()222221122121211(786)(32)2020x x x x x x x x =++==+++ 显然()V x 是正定的，()2()V x W x =是负定的，故上述方程组的零解是渐进稳定的. 2.3 线性类比法线性类比法是将一些非线性系统当作线性系统，用类比的方法构造出需要的Liapunov 函数.例3设1()f x 连续可导，(0)0f =，讨论微分方程组111222211222()dx f x a x dtdx a x a xdt⎧=+⎪⎪⎨⎪=+⎪⎩（2.5） 零解的稳定性.解当1111()f x a x =时，上述非线性系统（2.5）就是线性系统11111222211222dx a x a x dtdx a x a xdt⎧=+⎪⎪⎨⎪=+⎪⎩（2.6） 线性系统（2.6）的特征方程为2112211221221()0a a a a a a λλ-++-=容易看出，当1122112212210,0a a a a a a +<->时，特征方程有两个根都有负实部，线性系统（2.6）的零解是渐进稳定的.非线性系统（2.5）的零解稳定性无法用特征根的方法来判断，但可以用类似于线性系统的Liapunov 函数去判断其稳定性.事实上，对线性系统（2.6）取2121122112212211(,)()()V x x a a a a a a x =+-则12(,)V x x是半负定的函数，利用巴尔巴欣公式得221211221221122112211(,)()()22V x x a a a a x a x a x =-+-12(,)V x x 是正定函数，所以线性系统（2.6）的零解是稳定的.对比非线性系统（2.5）和线性系统（2.6），线性系统（2.6）中的11a 相当于非线性系统（2.5）中的11()f x x .所以很自然地猜想，当112222122111()()0,0f x f x a a a a x x +<->（2.7）成立时，非线性系统（2.5）的零解很可能是稳定的，注意到12112212211112212211101()()2x a a a a x a a a a x dx -=-⎰由此类比构造与线性系统类似的V 函数121122212211122112201()1(,)()()2x f x V x x a a a x dx a x a x x =-+-⎰ 当式（2.7）的条件成立时，12(,)V x x 正定，计算全导数得12221122111122(2.5)(,)(())(())V x x a f x a a x f x a x =-++22112222112212211222()((())())a x a x a f x a x a a x a x -+-+21122122122111()()()()0f x f x a a a a x x x =-+≤ 所以，非线性系统（2.5）的零解是稳定的.3 总结本文中我们仅仅给出了两种构造Liapunov 函数的方法，也很容易看出即便是这两种发放也都有其局限性. 仅仅对于常系数线性系统才可以使用巴尔巴欣公式，并且当微分方程组的阶数太大时，利用巴尔巴欣公式将会产生巨大的计算量.线性类比法需要我们先找到非线性系统对应的线性系统的Liapunov 函数，然后类比构造出非线性系统的Liapunov 函数，但是在实际问题中，往往很难进行类比，更多的要靠个人的经验总结.在教材上还介绍了能量函数法、分离变量法、变梯度法，微分方程组的形式结构不同，就要用不同的方法，遗憾的是至今仍未形成构造Liapunov 函数的通用方法，找到通用方法最大的困难是依赖于一个未知的V 函数的存在，然而构造这种V 函数又没有一般的规律可循，这是一种技巧性问题，更大的可能只能靠个人的经验.但是，我相信随着学科的发展，数学家们会继续研究构造Liapunov 函数的一般方法，Liapunov 函数的构造定会取得重大突破.4 参考文献[1]马知恩，周义仓，李承治. 常微分方程定性与稳定性方法[M]. 北京:科学出版社,2015.6. [2]张志芬，丁同仁，黄文灶，董镇喜. 微分方程定性理论[M]. 北京: 科学出版社,1985. [3]王高雄，周之铭，朱思铭等. 常微分方程[M]. 北京：高等教育出版社,2006.7.。
离散周期lyapunov方程和离散周期riccati方程的迭代算法
摘要作为线性时变系统的最简单形式，线性周期系统由于其广泛的应用，一直是学者们研究的热点。
线性周期系统，是一类系数矩阵带有周期性的线性系统，在各个领域中都有着广泛的应用。
为了研究离散周期系统的稳定性问题，离散周期Lyapunov方程的求解就显得至关重要。
同样，在进行离散周期系统的线性二次最优状态反馈控制器的设计时，需要用到离散周期Riccati方程的解。
基于这样的研究背景，本文针对离散周期系统下的Lyapunov方程和Riccati方程，给出了其求解的迭代算法。
针对离散周期Lyapunov方程，推导出了相应的迭代算法，分别对零初始条件和任意初始条件的情况给出了严谨的收敛性证明，并通过数值仿真验证了算法的有效性。
并且将最新估计信息的思想引入了迭代算法，得到了新的基于最新估计信息的迭代算法，同样对给出了算法在零初始条件下和非零初始条件下，迭代算法的严谨的收敛性证明，利用数值仿真例子证明了算法是有效并且收敛的。
并且通过对两种算法的数值仿真对比发现，基于最新估计信息的迭代算法的收敛速度要快于原始的迭代算法，从而验证了加入最新估计信息的迭代算法的优越性。
针对推导出的离散周期Riccati方程的迭代算法，给出了其在零初始条件下的收敛性证明，并通过数值仿真验证了算法的有效性，同样，为了改进算法，加入了最新估计信息，得到了新的基于最新估计信息的迭代算法。
同样对该算法的收敛性进行了严谨的证明与数值仿真验证，说明了该算法是有效可用的。
针对两种方程的迭代算法，为了研究最新估计信息对迭代算法的影响程度，引入了加权的思想，得到了带权重因子的新的迭代算法，并进行了收敛性证明。
通过数值仿真，给出了不同权重因子下的收敛性曲线，通过对比可以看出当全部使用最新估计信息时，算法的收敛速度最快，由此可见，加入最新估计信息能有效提高迭代算法的收敛速度。
关键词：离散周期系统；Lyapunov方程；Riccati方程；迭代算法AbstractAs the simplest form of time-varying linear systems, periodic linear systems have been attracting much attention during the past several decades. This is partially because this type of systems has very wide application. To investigate the stabilization problem of the periodic linear systems, it is important to achieve the solution of the periodic Lyapunov matrix equation. Similarly, the design of linear quadratic optimal state feedback controller based on the robust control is related to the stabilizing positive definite solution of Riccati equation. Based on this research background, we propose iterative algorithms for solving discrete-time periodic Lyapunov matrix equation and discrete-time periodic Riccati matrix equation.Iterative algorithms for discrete periodic Lyapunov equations are derived, respectively to the zero initial conditions and arbitrary initial conditions. And the proof of convergence is given. The effectiveness of the algorithm is verified by numerical simulation. And the latest information estimation theory is into the iterative algorithm, the proof of the convergence is also given. The validity of the algorithm is verified by numerical simulations. Finally, the simulation analysis of the two algorithms find that the convergence rate of the iterative algorithm based on the estimation of the latest information is faster than the original algorithm. It proves the superiority of the iterative algorithm adding the latest information of the estimation.Iterative algorithm for discrete periodic Riccati equations is derived, given the zero initial condition of convergence, and the effectiveness of the algorithm is verified through numerical simulation. In order to improve the algorithm with the latest estimate information, a new iterative algorithm based on the information of the latest estimation is given. The convergence of the new algorithm is proved and the validity of the algorithm is verified by numerical simulation. Through numerical simulation, the convergence curves of different weighting factors are given. It found that using the latest estimate information, the convergence speed is the fastest. Therefore, adding the latest estimation information can effectively improve the convergence speed of iterative algorithm.Key words：discrete-time linear periodic system，periodic Lyapunov equations，periodic Riccati equations，iterative algorithms目录摘要 (I)ABSTRACT ..................................................................................................................... I I 第1章绪论 . (1)1.1课题的来源及研究的背景意义 (1)1.2国内外在该方向上的研究现状及分析 (2)1.3本文的主要研究内容 (6)第2章离散周期系统Lyapunov方程快速迭代算法 (8)2.1相关的概念与性质 (8)2.2原始迭代算法 (9)2.2.1显式迭代算法 (9)2.2.2数值仿真 (12)2.3基于最新估计信息的迭代算法 (16)2.3.1显示迭代算法 (16)2.3.2数值仿真 (19)2.4本章小结 (24)第3章离散周期Riccati方程的迭代算法 (25)3.1相关的概念与性质 (25)3.2问题的描述 (25)3.3原始迭代算法 (25)3.3.1显示迭代算法 (26)3.3.2数值仿真 (28)3.4基于最新估计信息的迭代算法 (29)3.4.1显示迭代算法 (30)3.4.2数值仿真 (32)3.5本章小结 (34)第4章离散周期Riccati方程的加权最新估计迭代算法 (35)4.1 带加权因子的快速迭代算法 (35)4.2数值仿真 (37)4.3本章小结 (39)结论 (40)参考文献 (41) (45)致谢 (46)第1章绪论1.1课题的来源及研究的背景意义随着对控制系统的研究越来越深入，人们发现，许多生活中的系统是线性周期系统。
lyapunov-schmidt方法
Lyapunov-Schmidt方法是一种用于非线性方程组的求解的数值方法。
它是由俄罗斯数学家Aleksandr Lyapunov和德国数学家Ernst Schmidt分别在19世纪和20世纪提出的。
这种方法在处理非线性问题时非常有效，并且在应用数学和工程领域得到了广泛的应用。
Lyapunov-Schmidt方法的核心思想是将原始的非线性方程组转化成一系列线性方程组，从而简化求解过程。
这种方法的优势在于可以通过有限步骤来逼近非线性方程组的解，从而大大提高了求解效率。
下面我们将详细介绍Lyapunov-Schmidt方法的原理和应用。
1. Lyapunov-Schmidt方法的原理Lyapunov-Schmidt方法的原理是通过引入一组正交归一的特征函数，将原始的非线性方程组转化为一系列正交归一的线性方程组。
这样一来，原始的非线性方程组就被分解成了一系列互相独立的线性方程组，从而使得求解过程变得更加简单和高效。
2. Lyapunov-Schmidt方法的应用Lyapunov-Schmidt方法在科学和工程领域有着广泛的应用。
比如在物理学中，通过Lyapunov-Schmidt方法可以求解复杂的非线性波动方程，从而对物质的运动和变形进行研究。
在工程领域，Lyapunov-Schmidt方法可以用于求解具有非线性特性的结构力学问题，如弹性体的变形和弹性波的传播等。
3. 使用案例我们以一个简单的非线性方程组为例来说明Lyapunov-Schmidt方法的求解过程。
假设我们有一个非线性方程组：f(x, y) = 0g(x, y) = 0我们希望求解这个方程组的解。
我们可以通过Lyapunov-Schmidt方法将原始的非线性方程组转化为一系列正交归一的线性方程组：Φ1(x, y) = 0λ1(Φ1x + Φ1y) = 0Φ2(x, y) = 0λ2(Φ2x + Φ2y) = 0...我们可以通过求解这一系列线性方程组来逼近原始的非线性方程组的解。
lyapunov指数计算
lyapunov指数计算
Lyapunov指数（Lyapunov exponent）是一种用于描述动态系统混沌性质的指标。
在数学上，Lyapunov指数是描述线性化系统的稳定性的指标，它可以判断非线性系统是否具有混沌性质。
计算Lyapunov指数的基本过程如下：
1.首先，选择一个合适的初始状态，并计算该状态在系统中的轨迹。
2.然后，选取一个邻域，计算在该邻域内的状态与初始状态的差异随时间的变化情况。
3.对于每个时间步长，计算邻域中的点向初始状态点移动的距离与时间的比值。
4. 重复上述步骤，直到获得足够的数据，然后计算Lyapunov指数。
5. Lyapunov指数表示在该系统中相邻轨迹的指数级别分离速度。
具体计算Lyapunov指数的过程比较复杂，一般需要借助计算机进行模拟和计算。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
推广 形式：
u w* w* u du,
w
w w*
u w* u w du, * du, w f u
*
w
w
w*
f (u ) f ( w* ) du, f u
w
w*
f (u ) f ( w* ) du. u
Lyapunov 函数的构造
四、一类SI模型的全局稳定性
x*
u
y*
u
ky*
v*
u
c
w*
?
u
dL 则 1 dt
* p 0 x x * * x 2 w w 0 * c ww* x x M
2
x* xy*v y v* x v 3 x x* yv* y *v 0,
Lyapunov 函数的一些构造方法
李建全
空军工程大学(西安)
E-mail: jianq_li@
Lyapunov 函数的构造
一、引言 二、带有免疫的病毒动力学模型 三、Lyapunov函数形式的推广 四、一类SI模型的全局稳定性
Lyapunov 函数的构造
一、引言
全局渐近稳定性 (GAS)是动力系统研究的重要 问题。 对于高维自治微分系统 ，证明GAS最常用的方法是 Lyapunov - LaSalle定理， 即若自治系统 * x' f x, x Rn
Lyapunov 函数的构造
二、带有免疫的病毒动力学模型
x 健康细胞； y 被感染细胞； v 游离病毒； w 免疫细胞.
其中f w可分别取w,1, w /1 w.
M有惟一正平衡点 假设系统 E* x* , y* , v* , w* ,
下面通过构造 Lyapunov 函数来证明其GAS.
* *
dL 并且 1 dt
y v 0 x x* , * * , w w* . y v M
OK
Lyapunov 函数的构造
函数 3、 情形f w w / 1 w. 定义Lyapunov * * x ux y u y p w 1 u u w* x*v* v u v* L du du du du
结论：
Lyapunov 函数的构造
由于
则(IN)变为
Lyapunov 函数的构造
(IN) 定义Lyapunov函数 则有
OK
Lyapunov 函数的构造
补充：
取
条件：
Lyapunov 函数的构造
上面介绍的构造Lyapunov函数方法也可 用于如下模型：
和
A class of Lyapunov functions and the global stability of some epidemic models with nonlinear incidence, Journal of Applied Analysis and Computation , 6(2016)38-46
Lyapunov 函数的构造
谢 谢 ！
函数 2、 情形f w 1. 定义Lyapunov
*
dL 则 1 dt
* x x p 0 * * x 2 w w * x x M c
2ຫໍສະໝຸດ x* xy*v y v* x v 3 x x* yv* y *v 0,
x * y * * * v *
* * *
函数 1、 情形f w w. 定义Lyapunov
Nonlinear Analysis—Real World Applications: 12(2011)2163-2173; 13(2012)2006-2016.
Lyapunov 函数的构造
* y u y p w u w* 2 p u x* x*v* v u v* w w du ,, L* du * du du * * * x y v u u ky u 2 ccw u x
Lyapunov 函数的构造
* w p u w ux u y x v u v du, L du du du w* * x y v c u u u ky u dL1 x x* p0 w w* * 2 * 则 x 2 * c w w dt M x x * * * x xy v y v * * x v 3 x x* yv* y *v 0, v dL1 * y * 并且 0 x x , * * , w w . OK y v dt M
在区域内有惟一平衡点 E *， 并且能找到Lyapunov 函数L， 使得dL/dt * 在内关于E *负定或半负定， 同时系统(*)在集 合 dL/dt * 0 上的最大不变集为E * ， 则 E * 在内GAS。
证明的关键在于构造恰当的Lyapunov函数， 并证明其沿系统解的全导数是负定或半负定的。
* *
dL 并且 1 dt
y v 0 x x* , * * , w w* . y v M
OK
Lyapunov 函数的构造
三、Lyapunov函数形式的推广
* y u y u x* x*v* v u v* 其中 F * du * du du. * * x y v u u ky u x