现代控制理论理论教案
绪论
【教学目的】 了解现代控制理论的基本原理及方法,以便进行系统分析与设计,同时为进一步学习现代控制理论打下较扎实的基础。
【教学重点】 了解控制理论发展的三个阶段并掌握各阶段的主要任务。
【教学方法及手段】 课堂教学
【课外作业】 阅读教材
【学时分配】 2学时
【教学
内容
财务内部控制制度的内容财务内部控制制度的内容人员招聘与配置的内容项目成本控制的内容消防安全演练内容
】
本教材绪论部分主要讲述了以下几个问题:
一、控制理论发展简况
1)古典控制理论:研究对象以单输入、单输出线性定常系统为主,以传递函数为系统的基本描述,以频率法和根轨迹法为主要分析与设计手段。
2)现代控制理论以状态状态空间模型为基础,可研究多输入、多输出、时变、非线性等各种对象;研究系统内部结构的关系提出了能控性、能观测性等重要概念,提出了不少设计方法。
3)大系统与智能控制阶段。
二、现代控制理论的基本内容
(1)线性多变量系统理论。这是现代控制理论中最基础、最成熟的部分。它揭示系统的内在想律,从能控性、能观测性两个基本概念出发,研究系统的极点配置、状态观测器设计和抗干扰问题的一般理论。
(2)最优控制理论。在被控对象数学模型已知的情况下,寻求一个最优控制规律(或最优控制函数),使系统从某一个初始状态到达最终状态并使控制系统的性能在某种意义下是最优的。
(3)最优估计理论。在对象数学模型已知的情况下,最优估计理论研究的问题是如何从被噪声污染的观测数据中,确定系统的状态,并使这种估计在某种意义下是最优的。由于噪声是随机的,而且是非乎稳随机过程(随机序列),这种憎况下的状态估计是卡尔曼提出和解决的,故又称卡尔曼滤波。这种滤波方法是保证状态估计为线性无偏最小估计误差方差的估计。
(4)系统辨识与参数估计。这是基于对象的输入、输出数据、在希望的估计准则下,建立与对象等价的动态系统(即建立对象的数学模型),由于效学模型一船地说,是由阶致和参数决定的。因此,要决定系统的阶数和参数(即参数估计)。
三、本课程的基本任务
该课程是工业自动化专业的一门重要的专业基础课程。通过这门课的学习了解现代控制理论的基本原理及方法,以便进行系统分析与设计,同时为进一步学习现代控制理论打下较扎实的基础。所谓系统分析,就是指在规定的条件下,对数学模型已知的性能进行分析。系统分析包括定量分析和定性分析。定量分析是通过系统对某一个输人信号的实际响应来进行的;定性分析则研究系统能控性、能观测性、稳定性和关联性等一般特性。各种设计方法往往来源于系统分析。因此,系统分析是十分重要的。所谓系统设计,就是构造一个能完成给定任务的系统,这个系统具有所希望的瞬态,稳态性能以及抗干扰性能。一般地说,设计过程不是一个简单的一次能完成的过程,而是一个逐步完善的过程。在这个过程中,有可能引入补偿器或调整某些参数。
第一章 控制系统的数学模型
第一节 状态空间表达式
【教学目的】 了解状态空间描述的基本概念,掌握根据物理机理来建立状态空间表达式。掌握状态空间表达式的建立方法。
【教学重点】 基本概念的剖析与掌握。
【教学难点】 掌握状态变量是确定系统状态的最小一组变量。
【教学方法及手段】 课堂教学
【课外作业】1.1
【学时分配】2学时
【教学内容】
一、状态、状态变量和状态空间
通过RLC电路讲清楚状态、状态变量、状态空间的基本概念。
二、状态空间表达式
通过RLC电路的状态方程的建立将其分析结果推广到一般情况,可得到以下各种情况:
1)多输入、多输出(MIMO)线性定常系统:
=
(1-1)
其中A为系统矩阵,B为输入矩阵,C为输出矩阵,D为直接传输矩阵或称关联矩阵。
2)单输入、单输出(SISO)系统:
EMBED Equation.3
(1-2)
3)多输入、多输出(MIMO)线性时变系统:
EMBED Equation.3
(1-3)
4)非线性时变系统:
EMBED Equation.3
Y=
(1-4)
5)非线性定常系统:
=
(1-5)
三、状态变量的选取
1)同一系统可以取不同的状态变量;
2)状态变量的选取是非唯一的;
3)系统状态变量的数目是唯一的。
四、状态空间表达式建立的举例
通过质量、弹簧、阻尼器系统和直流他励电动机的状态空间表达式的建立以了解实际系统的建模步骤及思想。
第二节 由微分方程求状态空间表达式
【教学目的】 掌握根据系统微分方程建立状态空间表达式的方法。
【教学重点】 状态方程的建立。
【教学难点】 不同形式状态方程的建立。
【教学方法及手段】 课堂教学
【课外作业】 1.5
【学时分配】 2学时
【教学内容】
微分方程中不含有输入信号导数项
一般情况下,系统的输入和输出关系由n阶微分方程描述:
(1-6)
(1-7)
二、微分方程含有输入信号的导数项
(1-8)
第三节 传递函数矩阵
【教学目的】 掌握系统传递函数矩阵也是线性定常系统的一种描述。
【教学重点】 系统传递函数矩阵的求解
【教学难点】 由状态空间表达式求系统传递函数阵
【教学方法及手段】 课堂教学
【课外作业】 1.7
【学时分配】 1学时
【教学内容】
1) 单输入、单输出系统的传递函数:
(1-39)讲例1-4
2) 传递函数矩阵:
;(1-43)讲例1-5,例1-5的特点为两输入、两输出系统,这有别于单位输入、单输出系统。
3) 闭环系统传递函数矩阵
(1-9)
4) 传递函数描述和状态空间描述的比较。见P19
第四节 离散系统的数学描述
【教学目的】 了解离散系统空间表达式的建立方法。
【教学重点】 差分方程、 脉冲传递函数化为离散系统状态空间表达式。
【教学难点】 离散系统空间表达式与连续系统表达式的区别。
【教学方法及手段】 课堂教学
【课外作业】 1-8,1-9
【学时分配】 1学时
【教学内容】
1)差分方程中不含有输入量差分项的状态空间表达式的建立;
以三维为例,
,
2)差分方程中含有输入量差分项的状态空间表达式的建立;
G、H、C同上
讲清例1-6并要求画出状态图
3)脉冲传递函数(矩阵)
, (1-10)
通过例1-7搞清离散系统的传递函数矩阵的求法。
第五节 线性变换
【教学目的】 通过研究线性变换关系得到便于应用且简单的状态空间表达式
【教学重点】 各种
标准
excel标准偏差excel标准偏差函数exl标准差函数国标检验抽样标准表免费下载红头文件格式标准下载
的状态空间表达式,如能控、能观、对角、约旦型。
【教学难点】 非奇异变换阵的选取
【教学方法及手段】 课堂教学
【课外作业】 1-11,1-12
【学时分配】 2学时
【教学内容】
1)等价系统方程
线性定常系统的方程为
EMBED Equation.3
通过线性变换
EMBED Equation.3 ,
,
,
,
于是转换后的系统方程为:
2)线性变换的基本特性
a、 线性变换不改变系统特征值;b、线性变换不改变系统的传递函数矩阵。
3)化系统矩阵A为标准形
a、 化A为对角阵;讲例1-8,1-9
b、 化A为约当阵
例:考虑由下式确定的系统:
试求其状态空间表达式之能控标准形、能观测标准形和对角线标准形。
解:能控标准形为:
能观测标准形为:
对角线标准形为:
讲例1-10化A为约旦型。
小 结
本章介绍了状态空间描述和传递函数短阵描述。介绍了从状态变量的定义、状态变量的选取到建立状态空间表达式的整个过程,对于线性定常系统,在初始松弛情况下,也可以来用传递函数矩阵描述。这两种描述在系统分析和设计中都有应用。至于采用何种描述,应视所研究的问题以及时这两种描述的熟悉程度而定。
一个系统,状态变量的数目是唯一的,而状态变量的选取是非唯一的。选取不同助状态变量,建立的状态空间表达式亦异。它们之间可以通过线性变换进行转换。本章介绍了线性变换定义、基本持性以及应用变换的方法获得几种标准形。线性变换的方法相当重要,本门课程很多章节中均要应用。
传递函数矩阵的描述与状态变量选择无关,即系统状态变量的不同选择,传递函数(短阵)是不改变的。
第二章 线性控制系统的运动分析
第一节 线性定常系统齐次状态方程的解
状态转移矩阵(由定义求
,由拉普拉斯变换求
)
【教学目的】 了解状态转移矩阵的基本概念及求法
【教学重点】 状态转移矩阵的两种求法
【教学难点】 由拉普拉斯变换求
【教学方法及手段】 课堂教学
【课外作业】 2.1
【学时分配】 2学时
【教学内容】
1)齐次方程
EMBED Equation.3 的解为
;
2)状态转移矩阵的基本性质.P41
[例2.1] 试求如下线性定常系统
的状态转移矩阵Ф(t)和状态转移矩阵的逆Ф-1(t)。
[解] 对于该系统,
其状态转移矩阵由下式确定
由于
其逆矩阵为
因此
由于Ф-1(t)=Ф(-t),故可求得状态转移矩阵的逆为
第二节 状态转移矩的求法(凯莱-哈密尔顿(Caley-Hamilton)
定理, 对角线标准形与Jordan标准形法)
【教学目的】 了解状态转移矩阵的另外两种求法
【教学重点】 对角线标准形
【教学难点】 凯莱-哈密尔顿(Caley-Hamilton)定理
【教学方法及手段】 课堂教学
【课外作业】2.2,2.4
【学时分配】2学时
【教学内容】1)化eAt为A的有限项法(Caley-Hamilton定理法);
2)对角线标准形与Jordan标准形法
例:
[解] 该矩阵的特征方程为
因此,矩阵A有三个相重特征值λ=1。可以证明,矩阵A也将具有三重特征向量(即有两个广义特征向量)。易知,将矩阵A变换为Jordan标准形的变换矩阵为
矩阵P的逆为
EMBED Equation.3
注意到
可得
=
例:
试用前面介绍的两种方法计算
[解] 方法一 由于A的特征值为0和-2(λ1=0,λ2= -2),故可求得所需的变换矩阵P为
P =
因此,由式(2.10)可得
方法二 由于
可得
因此
第三节 线性定常系统非齐次状态方程的解
第四节 线性离散系统的运动分析
【教学目的】掌握线性定常系统非齐次状态方程的解及线性离散系统的运动分析
【教学重点】 线性定常系统非齐次状态方程的解
【教学难点】 线性离散系统的运动分析
【教学方法及手段】 课堂教学
【课外作业】2.6,2.20
【学时分配】2学时
【教学内容】
给定线性定常系统非齐次状态方程为
Σ:
(2-28)
其中,
,且初始条件为
。
将方程(2.28)写为
在上式两边左乘e-At,可得
将上式由O积分到t,得
故可求出其解为
(2-31)
(2.32)
式中
为系统的状态转移矩阵。
[例2.2] 求下列系统的时间响应:
式中,u(t)为t = 0时作用于系统的单位阶跃函数,即u(t)=1(t)。
[解] 对该系统
状态转移矩阵
已在前例中求得,即
因此,系统对单位阶跃输入的响应为:
或
如果初始状态为零,即X(0)=0,可将X(t)简化为
小 结
本章对系统运动的分析是通过求系统方程的解来进行的。状态方程是矩阵微分(差分)方程,输出方程是矩阵代数方程。因此,求系统方程的解的关键在于求状态方程的解。
线性系统方程曲解是借助状态转移矩阵来表示的。本章介绍了状态转移矩阵的定义、基本性质和求解方法。重点介绍了线性定常系统
状态转移矩阵的四种计算方法。有了状态转移矩阵,就可以求出系统在初始状
态激励下的自由运动(齐次状态方程的解)以及在输入向量作用下的强迫运动(非齐次状态方程的解)。应当指出,系统自由运动轨线的形态是由状态转移短阵决定的,也就是由A唯一决定的。然而对一个系统来说,A是一定的,因此只有靠人为地采取
措施
《全国民用建筑工程设计技术措施》规划•建筑•景观全国民用建筑工程设计技术措施》规划•建筑•景观软件质量保证措施下载工地伤害及预防措施下载关于贯彻落实的具体措施
(如第五章的状态反馈和输出反馈)来改造自由运动的形态。
状态x(t)求出后,即可求出系统的输出y(t)。不同的输入向量,响应y(t)不同。但是只要有了y(t)就可以按经典控制理论中介绍的时域分析法来定量地分桥系统的性能。由于这个响应y(t)是针对某个控制u(t)而言的,这就为用u(t)来达到希望的y(t)形态提供了可能(见第五章第六节)。
实验一 状态空间控制模型系统仿真及状态方程求解
【实验目的】 借助Matlab工具在计算机上实现前二章所讲的
重要内容。
【实验重点】 各种不同状态方程之间的转换及状态方程的求解、
【实验难点】 对角型、约旦型、模态型的转换
【教学方法及手段】 上机实验。
【课外作业】认真写实验报告,复习巩固实验内容
【学时分配】2学时
第三章 控制系统的能控性和能观测性
【教学目的】 掌握系统的能控性的概念及其判据。
【教学重点】 线性连续定常系统的能控性的判断。
【教学难点】 能控性判据概念的理解。
【教学方法及手段】 课堂教学
【课外作业】3-1
【学时分配】2学时
【教学内容】
在经典控制理论中,着眼点在于研究对系统输出的控制。对于一个单输入—单输出系统来说,系统的输出量既是被控量,又是观测量。因此,输出量明显地受输人信号控制,同时,也能观测,即系统不存在能控、不能控和能观测、不能观测的问题。 现代控制理论着眼点在于研究系统状态的控制和观测。这时就遇到系统的能控性和能观测性问题了。
能控性(controllability)和能观测性(observability)深刻地揭示了系统的内部结构关系,由R.E.Kalman于60年代初首先提出并研究的这两个重要概念,在现代控制理论的研究与实践中,具有极其重要的意义,事实上,能控性与能观测性通常决定了最优控制问题解的存在性。例如,在极点配置问题中,状态反馈的的存在性将由系统的能控性决定;在观测器设计和最优估计中,将涉及到系统的能观测性条件。
通过例3-1、3-2可知,研究系统的状态变量与输人信号之间的关系时,存在能控与不能控的问题。
系统能观测问题是研究测量输出变量y去确定系统状态变量的问题。通过例3-3可知,状态x存在能观测和不能观测的问题。
至此,我们可以知道,在基于状态空间描述的现代控制理论中,存在状态能控性和能观测性问题。这是两个反映系统构造特性的基本概念。在本章中,我们的讨论将限于线性系统。将首先给出能控性与能观测性的定义,然后推导出判别系统能控和能观测性的若干判据。
第一节 能控性及其判据
一、线性定常系统的能控性及其判据
(一)能控性定义
线性定常系统状态方程为
=
(3-1)
其中x、u分别为n、r维向量,A、B为满足矩阵运算的常值矩阵。若给定系统的一个初始状态可为x(t0)(t0可为0),如果在的有限时间区间[t0,t1]内,存在容许控制u(t)使x(t1)=0,则称系统状态在t0时刻是能控的;如果系统对任意一个初始状态都能控,则称系统是状态完全能控的,简称系统是状态能控的或系统是能控的。
由这个定义可知:
1)系统能控性定义中的初始状态x(t0)是状态空间中任意的非零有限点,控制的目标是状态空间坐标原点(有的文献称为达原点的能控性)。
2)如果在时间区间[t0,t1]内存在容许控制u(t),使系统从状态空间坐标原点推向预先指定的状态x(t1),则称为状态能达性。由于连续系统状态转移矩阵是非奇异的.因此可以证明系统能控性与能达性是等价的。 ·
3)在能控性研究中,我们考察的并不是x(t0)推向x(t1)=0的时变形式,而是考察能控状态在状态空间中的分布。很显然,只有整个状态空间中所有的有限点都是能控的,系统才是能控的。
(二)能控性判据
判据一:若式(3-1)系统能控,则
能控性矩阵
满秩。即
判据一的证明从略,结合具体例子介绍其方法。
[例] 考虑由下式确定的系统:
由于
即Q为奇异,所以该系统是状态不能控的。
[例] 考虑由下式确定的系统:
(3-2)
对于该情况,
即Q为非奇异,因此系统是状态能控的。
判据二:由于状态能控的条件是A的特征向量互异,关于定常系统能控性的判据很多。除了上述的代数判据外,本小节将给出一种相当直观的方法,这就是从标准形的角度给出的判据。
考虑(3-2)的线性系统。 如果A的特征向量互不相同,则可找到一个非奇异线性变换矩阵P可将A阵转换为对角阵,当且仅当转换后的输入矩阵
没有一行的所有元素均为零时,系统才是状态能控的。如果式(3-2)中的矩阵A不具有互异的特征向量,则不能将其化为对角线形式。在这种情况下,可将A化为Jordan标准形。例如,若A的特征值分别λ1,λ1,λ1,λ4,λ4,λ6,…,λn,并且有n - 3个互异的特征向量,那么A的Jordan标准形为
(3-3)
与每个Jordan块最后一行相对应的
的任一行元素不其中,在主对角线上的3×3和2×2子矩阵称为Jordan块。
假设能找到一个变换矩阵S,使得
如果利用
x = S z
(3-4)
定义一个新的状态向量z,将式(3-4)代入式(3-2)中,可得到
(3-5)
从而式(3-5)确定的系统的状态能控性条件可表述为,当且仅当:(1)式(3-3)中的矩阵J中没有两个Jordan块与同一特征值有关;(2)全为零;(3)对应于不同特征值的
的每一行的元素不全为零时,则系统是状态能控的。
[例] 下列系统是状态能控的:
下列系统是状态不能控的:
判据三: 用传递函数矩阵表达的状态能控性条件
状态能控的条件也可用传递函数或传递矩阵描述。
状态能控性的充要条件是在传递函数或传递函数矩阵中不出现相约现象。如果发生相约,那么在被约去的模态中,系统不能控。
[例] 考虑下列传递函数:
显然,在此传递函数的分子和分母中存在可约的因子(s+2.5)(因此少了一阶)。由于有相约因子,所以该系统状态不能控。
当然,将该传递函数写为状态方程,可得到同样的结论。状态方程为
由于
即能控性矩阵
的秩为1,所以可得到状态不能控的同样结论。
第二节 能观测性及其判据
【教学目的】 掌握系统的能观性的概念及其判据。
【教学重点】 线性连续定常系统的能观性的判断。
【教学难点】 能观性判据概念的理解。
【教学方法及手段】 课堂教学
【课外作业】3-2
【学时分配】2学时
【教学内容】
现在讨论线性系统的能观测性。考虑零输入时的状态空间表达式
(3-6)
式中,
。
如果每一个状态x(to)都可通过在有限时间间隔to≤t≤t1内,由y(t)观测值确定,则称系统为(完全)能观测的。本节仅讨论线性定常系统。不失一般性,设to=0。
能观测性的概念非常重要,这是由于在实际问题中,状态反馈控制遇到的困难是一些状态变量不易直接量测。因而在构造控制器时,必须首先估计出不可量测的状态变量。在“系统综合”部分我们将指出,当且仅当系统是能观测时,才能对系统状态变量进行观测或估计。
在下面讨论能观测性条件时,我们将只考虑由式(3-6)给定的零输入系统。这是因为,若采用如下状态空间表达式
则
从而
由于矩阵A、B、C和D均为已知,u(t)也已知,所以上式右端的最后两项为已知,因而它们可以从被量测值y(t)中消去。因此,为研究能观测性的充要条件,只考虑式(3-6)所描述的零输入系统就可以了。
判据一、线性定常系统状态能观测性及其判据
考虑由式(3-6)所描述的线性定常系统。将其重写为
易知,其输出向量为
将
写为A的有限项的形式,即
因而
或
(3-7)
显然,如果系统是能观测的,那么在0≤t≤t1时间间隔内,给定输出y(t),就可由式(3-7)唯一地确定出x(0)。可以证明,这就要求nm×n维能观测性矩阵
的秩为n。
由上述分析,我们可将能观测的充要条件表述为:由式(3-6)所描述的线性定常系统,当且仅当n×nm维能观测性矩阵
的秩为n,即
时,该系统才是能观测的。此为判据一。
[例] 试判断由式
EMBED Equation.3
所描述的系统是否为能控和能观测的。
[解] 由于能控性矩阵
的秩为2,即
,故该系统是状态能控的。
对于输出能控性,可由系统输出能控性矩阵的秩确定。由于
的秩为1,即
,故该系统是输出能控的。
为了检验能观测性条件,我们来验算能观测性矩阵的秩。由于
的秩为2,
,故此系统是能观测的。
判据二、状态能观测性条件的标准形判据
考虑由式(3.13)和(3.14)所描述的线性定常系统,将其重写为
(3-8)
设非奇异线性变换矩阵P可将A化为对角线矩阵,如果m×n维矩阵
的任一列中都不含全为零的元素,那么系统是能观测的。
如果不能将式(3-8)变换为对角线标准形,则可利用一个合适的线性变换矩阵P,将其中的系统矩阵A变换为Jordan标准形。系统能观测的充要条件为:(1)J中没有两个Jordan块与同一特征值有关;(2)与每个Jordan块的第一行相对应的矩阵
列中,没有一列元素全为零;(3)与相异特征值对应的矩阵
列中,没有一列包含的元素全为零。为了说明条件(2),在下例中,对应于每个Jordan块的第一行的
列之元素用下划线表示。
[例] 下列系统是能观测的:
显然,下列系统是不能观测的:
判据三、用传递函数矩阵表达的能观测性条件
类似地,能观测性条件也可用传递函数或传递函数矩阵表达。此时能观测性的充要条件是:在传递函数或传递函数矩阵中不发生相约现象。如果存在相约,则约去的模态其输出就不能观测了。
[例3.6] 证明下列系统是不能观测的。
式中
[解] 由于能观测性矩阵
注意到
即
,故该系统是不能观测的。
事实上,在该系统的传递函数中存在相约因子。由于X1(s)和U (s)之间的传递函数为
又Y (s)和X1(s)之间的传递函数为
故Y(s)与U(s)之间的传递函数为
显然,分子、分母多项式中的因子(s+1)可以约去。这意味着,该系统是不能观测的,或者说一些不为零的初始状态x(0)不能由y(t)的量测值确定。
第四节 离散系统的能控性和能观测性
【教学目的】 掌握离散系统的能控性、能观性的概念及其判据。
【教学重点】 线性定常离散系统的能控性、能观性的判断。
【教学难点】 能控性、能观性判据概念的理解。
【教学方法及手段】 课堂教学
【课外作业】3-9,3-10
【学时分配】2学时
【教学内容】
一、能控性定义
关于离散系统的能控性和能观测性问题,几乎与连续系统完全类似地有一套相应的理论和方法。因此本节只作扼要地介绍。线性定常离散系统方程为
(3-9)
其中x(k)、u(k)、y(k)分别为n、r、m维向量,G、H、C为满足矩阵运算的矩阵。
对系统(3.33)的任一初始状态x(0),存在k>0,在有限时间区间[0,k]内,存在容许控制序列u(k),使得x(k)=0,则称系统是状态完全能控的,简称系统是能控的。由于在k时刻,有x(k)=0,有的书称为第k步能控。如果在有限时间区间[0,k]内,存在容许控制序列u(k),将系统从状态空间坐标原点x(k)=0推向预先指定的状态x(k),则称为能达性。在连续系统中,系统的能达性与能控性是等价的,而离散系统的能达性与的控性之间关系如何呢?离散系统与连续系统略有差别。在离散系统中,如果系数矩阵G是非奇异的.则能达性与能控性等价。也就是说,离散系统中的能达性和能控性等价是有条件的。
二、能控性判据
线性定常离散系统能控的充要条件为
满秩。
例3-11 线性定常离散系统状态方程为
试判别系统的能控性。
解:
故系统能控。
三、能观性定义
对于系统(3—9),根据有限个采样周期y(k),可以唯一地确定系统的任一初始状态x(0),则称系统是状态完全能观测的,简称系统是能观测的,有的书称为第k步能观测。同样也可以讨论系统的能检测性,而且离散系统的能检测性、能观测性之间的关系与连续系统略有差别。在离散系统中,只有系数短阵G是非奇异时,能检测性与能观测性才是等价的,也就是说,离散系统的能检测性和能观测性是有条件的等价。
四、能观测性判据
系统(3-9)能观测的充分必要条件是
型能观测性矩阵
的秩为
。即
例3-12 线性定常离散系统方程为
试判别系统的能观测性。
故系统能观测。
第五节 对偶原理
下面讨论能控性和能观测性之间的关系。为了阐明能控性和能观测性之间明显的相似性,这里将介绍由R.E.Kalman提出的对偶原理。
考虑由下述状态空间表达式描述的系统S1:
(3-10)
式中,
。
系统能控性是研究输入u(t)与状态x(t)之间的关系,而能观测性是研究输出y(t)与状态x(t)之间的关系。通过上面的讨论可以看到,能控性与能观测性,无论在概念上还是判据的形式上,都很相似。它给人们一个启示,即能控性与能观测性之间存在某种内在的联系。这个联系就是卡尔曼提出的对偶性。
现在.我们来构造一个系统
(3-11)
其状态图如下所示。
对偶系统有两个基本特征:
1.对偶的两个系统传递函数阵互为转置
2.对偶的两个系统特征值相同
3. 线性定常系统(3—10)和系统(3—11)为对偶系统,系统(3—10)的能控性等价于系统(3—11)的能观测性;而系统(3—10)的能观测性与系统(3—11)的能控性等价。这就是对偶原理。
例3—14 线性定常系统方程为
试判别系统能观测性。
解 该题可以宜接检查能观测性短阵的秩来判grj系统能观测性。但是为了熟悉对偶原理的应用,下面用检查其对偶系统能控性来判别系统能观测性。该题的对偶系统为
能控性矩阵
对偶系统能控。根据对偶原理知,原系统能观测。
实际上
,系统能观测,与按对偶原理判别结果一致。
小 结
能控性和能观测性是系统定性分析的重要内容之一。本章介绍能控性和能观测性的定义,导出了线性系统能控性、能观测性的定理。其中定理3—6和定理3—14是本章两个基本结果。因为导出它们所用的假定最少(只需假定连续性),因此可以最广泛地应用。若引入附加假定(连续可微性),则得到定理3—8和定理3—15,它们虽仅给出充分条件,但易于应用。对于线性定常系统,可以得到系统能控和能观穗的充分必要条件。如果将能控性、能观测性的定理一一对应列出,持会发现其间的对偶性。对偶原理搭起了控制问题和估计问题的桥梁,在理论和实际两方面具有很大意义。
第四章 控制系统的稳定性
稳定性是系统定性分析的又一个重要内容。
工程
路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理
实际中,可以应用的系统必须是稳定的。不稳定的系统是不能付诸实用的。在系统分析和设计中,不可避免地会遇到稳定性问题。
【教学目的】 掌握系统稳定性的概念及判据
【教学重点】 李亚甫诺夫意义下稳定性的定义
【教学难点】 有关稳定性几个重要概念的理解
【教学方法及手段】课堂教学
【课外作业】 4-1
【学时分配】 2学时
【教学内容】
随着科学技术的发展以及航空、航天工业发展的需要,控制问题由线性、定赏、单输入单输出系统问题刚E线性、时变、多输入多输出系统向题延伸。使得稳定性问题分析的复杂程度急剧地增加。那些在经典控制理论中行之有效的穆定性分析方法在此无能为力。必须寻求其它方法。
李亚甫诺夫在1892年发表了《运动稳定性一般问题》
论文
政研论文下载论文大学下载论文大学下载关于长拳的论文浙大论文封面下载
,建立了运动稳定性的一般理论和方法。他把分析常微分方程组稳定性的所有方法归纳为两种。第一种方法是求出常微分方程的解,分析系统的稳定性,这是一种间接方法;第二种方法是不需要求解激微分方程的解而能提供稳定性的信息。这是一种直接方法。由于求解非线性时变微分方程组的解是很困难的,甚至是不可能的。因此,李亚甫诺夫第二法就显得特别的重要。该方法研究系统稳定性是建立在这样一个事实之上的,即系统的一个平衡状态若为渐近稳定时,在外界作用下,系统能量要发生变化。但是,系统贮存的能量必将随着时间的增长而衰减,直至趋于平稳状态而使能量趋于最小值。
书上通过例-1以一个机械平移系统为例来说明这个问题。
李亚甫诺夫构造了所谓广义能量函数,称之为李亚甫诺夫函数,记成V(x,t)。当李亚甫诺夫函数不显含时间t,就记成V(x)。通过研究V(x,t)或V(x)及其沿系统状态软线运动随时间的变化率的定号性就可以给出系统稳定性的信息。
图4-2以二维系统为例,说明了能量函数E(x1,x2)随时间t的的增加而而连续减小。
李亚甫诺夫第二法是研究系统乎衡状态稳定性的。什么是系统平衡状态呢? 在例4—1中,x1=0,x2=0称为平衡状态。一般地说,系统的状态方程为
(4—1)
其初始状态为x(0)。系统的状态轨线x(t)是随着时间而变化的。当且仅当
有,
则称
为系统的平衡状态。由此可见,当状态轨线x(t)达到平衡状态时,如果系统不加输入,则状态就永远停留在平衡状态。
第一节 李亚甫诺夫意义下稳定性的定义
一、 稳定
二、渐近稳定
李亚甫诺夫意义下渐近稳定就是经典控制理论中所说的稳定。工程中的系统都要求是李亚甫诺夫意义下渐近稳定。
三、大范围渐近稳定
渐近稳定性是系统的一个局部稳定性概念。如果对于状态空间中,初始状态
是整个状态空间中的任何点,而从
出发的状态轨线
有
(4-2)
则称
=0为李亚甫诺夫意义下大范围渐近稳定或李亚甫诺夫意义下全局渐近稳定。当大范围渐近稳定与初始时刻
选择无关时,则称一致大范围渐近稳定。很显然,对于大范围渐近稳定的系统,其必要条件是整个状态空间中,只存在一个平衡状态。对于线性系统,只要系统
=0是渐近稳定的,则一定
是大范围渐近稳定的。
四、不稳定
第三节 李亚甫诺夫第二法
【教学目的】 掌握李亚甫诺夫第二法判断系统稳定性的几个定理。
【教学重点】 一致大范例渐近稳定的概念
【教学难点】 二次型定号性与稳定性的关系
【教学方法及手段】课堂教学
【课外作业】复习所讲内容
【学时分配】2学时
【教学内容】
在本章第一节中已简单地介绍了李亚甫诺夫第二法研究系统稳定性的基本方法,即构造一个与系统状态x有关的标量函数V(x,t)来表征系统的广义能量。V(x,t)称为李亚甫诺夫函数。研究V(x,t)及其沿状态轨线随时间的变化率的定号性,就可以得到有关系统的稳定性信息。换句话说,对一个系统来说,如果我们能够构造V(x,t),就能判断该系统的运动
核定性。
本节介绍李亚甫诺夫第二法判断系统稳定性的几个定理。
定理4-1 设系统状态方程为
(4-3)
在平衡状态
的某领域内,标量函数
,具有连续一阶导数,且满足
1、
2、
则
是一致渐近稳定的。如果
则
是一致大范例渐近稳定的。
定理4-2 设系统状态方程为
(4-4)
在平衡状态
的某领域内,标量函数
,具有连续一阶导数,且满足
1、
2、
3、除
平衡状态外,还有
的点,但不会在整条状态轨迹上有
,则
是一致渐近稳定的。如果
则
是一致大范例渐近稳定的。
定理4-3 设系统状态方程为
(4-5)
在平衡状态
的某领域内,标量函数
,具有连续一阶导数,且满足
1、
2、
则
是一致稳定的。
定理4-4 设系统状态方程为
(4-6)
在平衡状态
的某领域内,标量函数
,具有连续一阶导数,且满足
1、
2、
则
是不稳定的。
第四节 线性连续系统的稳定性
【教学目的】 掌握李亚甫诺夫第二法判断线性连续及离散系统的稳定性。
【教学重点】 连续系统李亚甫诺夫方程
【教学难点】 离散系统李亚甫诺夫方程
【教学方法及手段】 课堂教学
【课外作业】 4-4、4-5、4-6
【学时分配】 2学时
【教学内容】
设系统状态方程为
假定A是非奇异矩阵,这时系统存在唯一的平衡状态
。李亚甫诺夫函数
为状态变量
的二次型形式,即
(4-7)
其中P为
型正定的对称常值矩阵。显然有
,当要求
为渐近稳定时,
应为负定的。令
(4-8)
式中的Q阵为正定对称矩阵且满足
,式4-8称为李亚甫诺夫方程。且Q是正定矩阵。式(4-8)称为李亚甫诺夫方程。因为Q是正定短阵,则
<0,这就意味着沿
的任意轨线x(t),V(x)随时间单调减小,当
时,V(x)最终将趋于零。根据李亚甫诺夫稳定性定理4.1可知,
是一致渐近稳定的。因为是线性系统,故
是一致大范围渐近稳定的。
在运用李亚甫诺夫方程(4—8)判别系统
的稳定性时,先指定正定阵Q,再按式(4-7)求出P阵,然后检查P阵的正定性。由于Q阵的形式可以任意给定,并且最终的判断结果与正定阵Q的不同选择无关。故最方便也是最简单的选择是选取Q=I(单位阵)。这时李亚南诺夫方程就成为
(4—20)
根据式(4-7)求出P阵,用赛尔维斯特判据来检验其正定性,当P阵是正定阵时,则
为一致渐近稳定的,并且是一致大范围渐近稳定的。详细讲解例4-6。
线性定常离散系统的稳定性
线性定常离散系统的状态方程为
(4-9)
是系统的平衡状态。下面用李亚甫诺夫第二法来研究系统
的渐近稳定性问题。
对于式(4—21)描述的线性定常离散系统。假设G为
型奇异常阵,
是唯一的平衡状态。选取李亚甫诺夫函数为
式中,P为
正定的对称常值短阵。
显然有 V[
]>0,
而
的差分为
式中Q为正定对称常阵。而
(4-23)
称为李亚甫诺夫方程。与线性定常连续系统类似,判别系统
的渐近稳定性时,通常是给出一个正定对称常P,然后用式(4.23)求出P阵,并验证其正定性。如果P阵是正定的,则
为一致渐近稳定的,且是一致大范围斯近稳定的。讲解书上例4-7。
小 结
稳定性与能控性、能观测性一样都是系统的重要特性。本章介绍了李亚甫诺夫意义下稳定性的定义和李亚甫诺夫第二法分析系统平衡状态稳定性的定理,同时介绍了线性定常系统零状态响应的稳定性及其与平衡状态稳定性的关系,即平衡状态渐近稳定包含了BIBO稳定,而BIBO稳定的系统未必是平衡状态惭近稳定,只有当系统能控又能观测时,BIBO稳定的系统才是平衡状态渐近稳定。
线性定常系统的稳定性可以由传递函数的极点或由A的特征值来分析。也可以用李亚甫诺夫第二法来分析。应该指出,到目前为止,还没有构造李亚甫诺夫函数的一般方法,而靠经验与技巧。由于李亚甫诺夫第二法给出的结果是非线性系统稳定性的充分条件,所以,对某个系统而盲,构造不出李亚甫诺夫函数,我们不能说,该系统不稳定,只能说,无法提供有关系统稳定性的信息。
稳定性分析方法同样可应用到离散系统中去。只是线性离散系统的李亚甫诺夫方程形式和线性连续系统略有不同。
实验二 能控能观判据及稳定性判据
【实验目的】 借助Matlab工具在计算机上实现能控性及能观性判据
【实验重点】 李亚甫诺夫稳定性判据
【实验难点】 能控性、能观性判据
【教学方法及手段】 上机实验。
【课外作业】认真写实验报告,复习巩固实验内容
【学时分配】2学时
第五章 线性定常系统的综合
第一节 引 言 第二节 状态反馈和输出反馈
【教学目的】掌握状态反馈及输出反馈的概念
【教学重点】 状态反馈的意义
【教学难点】 状态反馈的的作用
【教学方法及手段】课堂教学
【课外作业】复习所讲内容
【学时分配】2学时
【教学内容】
在第二章,研究的是在己知系统的结构和参数情况下系统的运动,从而了解系统的运动形态。第三章介绍了系统的能控性和能观测性。第四章是系统稳定性问题。如果将上述研究的内容概括起来说,就是在已知系统的结构和参数情况下,研究系统的性能或特性,即所谓系统分析问题。 本章将研究线性定常系统的综合。这是一个与系统分析相反的命题,是在给定被控对象的情况下,通过设计控制器的结构和参数,使系统满足预先规定的性能指标要求。采用的方法是先涵量系统的状态,再用状态来确定被控对象上所加的控制输人,从而构成状态反馈系统。对状态反馈系统来说,能控性和能观测性同祥具有很重要的意义。采用状态反馈,对系统能控性和能观测性有无影响呢?这是本章讨论的重要内容之一。同时研究一个能控的系统,引入状态反蚀可以任意配置状态反馈系统的极点,保证系统具有所希望的瞬态性能和稳态性能;对于系统的状态变量无法测量但又要用它来实现反馈的情况,将介绍在系统朗观测条件下,通过状态重构方法。设计状态观测器。用重构状态实现状态反馈。本章还将研究用状态反馈进行系统解耦。
在经典控制理论中,利用系统的输出进行反馈,构成输出负反馈系统,可以得到较满意的系统性能;减小于扰对系统的影响;减小被控对象参数变化对系统性能的影响。因此,输出反馈控制得到了广泛的应用。在现代控制理论中,为了达到希望的控制要求,也采用反馈控制方法来构成反馈系统。这里采用的反馈控制有状态反馈和输出反馈两种。
一、状态反馈
线性定常系统方程为
(5-1)
其中状态x、输入u和输出y分别为n、r、m维向量。A、B、C、D为满足短阵运算的短阵。
假定有可能设置n个传感器,使全部状态变量均可用于反馈。其反馈控制律为
u=V-Kx (5-2)
其中 K为r
n型反馈增益矩阵;V为r维输入向量。构成的状态反馈系统如图5—1所示。
状态反馈系统方程为
=Ax十B(V—Kx)=(A一BK)x+BV
y=(C—DK)x+DV (5-3)
由方程(5—3)可知:
1)状态反馈不增加新的状态变
量o
2)状态反馈对输入矩阵B和直
接传畅矩阵D无影响o
3)系统的系数矩阵由A变成
(A-BK)。
4)输出矩阵由C变成(C—DK)。
系统的瞬态性能主要由系数矩
阵决定。A、B阵是已知的,不能改
变。K阵可以在一个很宽范围内选
择。因此,通过适当的方法选择反馈
阵K,就可以使系统达到希望的控制
目的。
二、输出反馈
在工程实践中,输出反馈也是常
用的。对方程(5—1)所描述的线性定
常系统,采用输出反馈控制律为
u=V-Hy (5-4)
其中 H为r
m型常值短阵。输出反馈系统的结构图如图5—2所示。
第3节 状态反馈系统的能控性和能观测性
【教学目的】掌握状态反馈对原系统的影响及具体设计方法。
【教学重点】 状态反馈对原系统能观测性及能控性的影响。
【教学难点】 状态反馈的设计
【教学方法及手段】课堂教学
【课外作业】5-1,5-2
【学时分配】2学时
【教学内容】
线性定常系统方程为
=Ax+Bu
y=Cx (5-6)
其中 x、u、y维数同前。如果引人状态反馈
u=V—Kx (5-7)
其中 v、K意义同前,则状态反馈系统方程为
=(A-BK)x+Bv
y=Cx (5-8)
对状态反馈系统来说,能控性和能观测性同样具有很重要的意义。那么,引入状态反馈的系统能控性、能观测性与未引入状态反馈的情况下的系统能控性、能观测性有什么关系呢?换句话说,状态反馈对系统能控性、能观测性有无影响呢?这个问题的结论是状态反馈不改变系统的能控性,但可能改变系统能观测性。
定理5.1 对于任何常值反馈阵K,状态反馈系统能控的充分必要条件是式(5-6)的系统能控。
例5-1 系统方程为
y=
可以验证,该系统能控、能观测。
现在引入状态反馈
u=V-
状态反馈系统方程为
rankQc=2, 该系统能控.
rankQ0=rank
即状态反馈系统不能观测,即引入状态反馈改变了这个系统的能观测性。
这个例子说明原系统是能观测的,引入状态反馈后,改变了原系统的能观测性。也可能出现相反的情况,即原系统不能观测,引入状态反馈后,状态反馈系统变成能观测了。一般地说,当用状态反馈配置的系统极点与原系统零点相同,即出现零、极点相消时,状态反馈就改变了系统的能现测性。
第四节 极点配置
状态反馈系统购稳定性和瞬态性能主要是由系统极点决定的。如果引人状态反馈将系统
的极点配置在S左半平面的希望位置上,则可以得到满意的系统特性。一个系统引入状态反
馈可以任意配置极点的条件是原系统能控。现在介绍单输入—单输出系统的极点配置。
系统方程为
=Ax+Bu
y=Cx
其中x为n维状态向量,u、y为标量。A、B、C分别为n
n、n
l、1
n型矩阵。
状态反馈律为
u=V—Kx
其中 K为1
n型的常值矩阵。
状态反馈系统方程为
=(A-bK)x+bv
y=Cx
状态图如图5—3所示。
若式(5—10)的系统能控,则引人状态反
馈可以任意配置式(5—12)的状态反馈系统
的极点。由式(5—12)可知,A、b一定,配置
系统的极点,就是确定K阵。通过计算合
适的K阵,将系统极点配置在S平面上所
希望的位置。确定K阵有许多方法。先
介绍一种常用的方法。
例5—2线性定常系统状态方程为
=
x+
引入状态反馈配置系统的极点为S1,2=-1
j;s 3=—2,试确定反馈阵K。
解:设反馈阵K=
A- bK=
_
=
希望的状态反馈系统特征多项式为
=
EMBED Equation.3 和
的S同次幂系数应相等,得到
EMBED Equation.3 =4
=4,
=1即K=[4 4 1]. 状态反蚀系统的状态图如图5-4所示
例5—3 设系统如图5—5所示。
采用输出反馈,如图5-6a所示
开环传递函数为g(s)=
令kF=10F,画出系统的根执迹团如图5-6b所示。可见当kF值增大到一定值后,闭环系统就变成不稳定的了。
现在采用状态反馈。结构图如图5—7a所示。
利用结构变换,图5—7a可变成如图5—7b形式。可见状态反馈系统等效为具有输出微分反馈的系统。当状态反馈系统的极点选择为s1,2=一0.708
时,求得k0=0.702,k1=-3.393,k2=4.77状态反馈系统开环传递函数为
gk(s)=
令kF=9.54F,画出状态反馈系统的根轨迹如图5—7c所示。比较图5—7c和图5—6b可知,当增益kF>0的所有情况,状态反馈系统皆稳定。由于g(s)的极点很接近gk(s)的零点,所以即使增益变化,状态反馈系统的极点位置变动很少。这表明状态反馈系统有好的相对稳定性以及系统响应受F变化的影响小。
第4节 状态重构和状态观测器
【教学目的】掌握状态重构和状态观测器的设计原理及方法。
【教学重点】 状态观测器的设计。
【教学难点】 状态重构概念的理解
【教学方法及手段】课堂教学
【课外作业】5-4
【学时分配】2学时
【教学内容】
引入状态反馈可以得到较好的系统性能。而实现状态反馈的前提是状态变量必须能用传感器测量得到。但是由于种种原因,状态变量并不是都可测量得到。例如,系统中的某些状态基于系统的结构特性或者是状态变量本身无物理意义,而无法测得;有些状态变量虽然可以测量得到,但应用的传感器价格很贵;有些状态信号很脆弱,在测量点易混进噪声,使得这些状态实际上难以应用。上述情况表明,我们得不到实际能应用的系统状态变量。运用状态反馈又必须有可应用的状态变量,怎么办呢?能否通过系统的输入量和输出量来构造系统的状态呢?回答是肯定的。可以根据系统的输入量、输出量和系统的结构、参数来实现系统的状态重构。实现状态重构的系统称为状态观测器。
系统方程为
=Ax+Bu
y=Cx (5-9)
x(t0)=x(0)
其中 x、u、y分别为n、r、m维向量,A,B,C为满足短阵运算的矩阵。设系统的状态不能得到。解决系统状态重构的一个直观想法是构造一个系统,该系统的输入为式(5—33)的系统的输入u(t);系统的结构、参数与式(5—33)的系统相同。于是得到如下的系统方程
=A
+Bu
y=C
小 结
本章是的控性与能观测性的应用。当线性定常系统能按时,通过引入状态反馈就能任意配置状态反馈系统的特征值,从而得到希望的系统性能;输出比例反馈系统是不能实现系统特征值任意配置的。采用输出动态反馈也可以实现系统特征值的任意配置。这一点在本章第四节中,通过对状态反馈系统进行结构变换.得到等价的输出动态反浊系统,已经证实。当系统状态不能够取得时,若系统能观测,可以设计状态观测器,用重构的状态向量实现反馈,构成具有状态观澜器的状态反馈系统。
实验三 状态反馈及状态观测器设计
【实验目的】在计算机上实现状态反馈和状态观测器的设计。
【实验重点】 状态反馈的设计。
【实验难点】 状态观测器的设计
【教学方法及手段】上机实验
【课外作业】认真写实验报告,复习巩固实验内容
【学时分配】2学时
15
_1176788856.unknown
_1176874289.unknown
_1182407254.unknown
_1182408517.unknown
_1182408829.unknown
_1182408954.unknown
_1182409025.unknown
_1182409295.unknown
_1182410232.unknown
_1916272007.unknown
_1182409838.unknown
_1182409291.unknown
_1182408966.unknown
_1182408969.unknown
_1182408973.unknown
_1182408957.unknown
_1182408943.unknown
_1182408949.unknown
_1182408930.unknown
_1182408685.unknown
_1182408687.unknown
_1182408812.unknown
_1182408686.unknown
_1182408617.unknown
_1182408683.unknown
_1182408684.unknown
_1182408618.unknown
_1182408615.unknown
_1182408616.unknown
_1182408518.unknown
_1182407496.unknown
_1182408411.unknown
_1182408515.unknown
_1182408516.unknown
_1182408513.unknown
_1182408514.unknown
_1182408512.unknown
_1182408267.unknown
_1182408268.unknown
_1182408266.unknown
_1182407359.unknown
_1182407385.unknown
_1182407389.unknown
_1182407374.unknown
_1182407266.unknown
_1182407355.unknown
_1182407351.unknown
_1182407258.unknown
_1176876524.unknown
_1176905825.unknown
_1176907039.unknown
_1176907210.unknown
_1176907708.unknown
_1176908745.unknown
_1182404397.unknown
_1182405648.unknown
_1182406090.unknown
_1182404614.unknown
_1176908795.unknown
_1176907952.unknown
_1176907502.unknown
_1176907174.unknown
_1176907186.unknown
_1176907056.unknown
_1176907155.unknown
_1176906268.unknown
_1176906639.unknown
_1176906837.unknown
_1176906893.unknown
_1176906667.unknown
_1176906350.unknown
_1176906002.unknown
_1176906157.unknown
_1176903301.unknown
_1176904504.unknown
_1176904962.unknown
_1176904997.unknown
_1176905617.unknown
_1176904522.unknown
_1176904527.unknown
_1176903489.unknown
_1176903515.unknown
_1176903392.unknown
_1176901127.unknown
_1176903143.unknown
_1176903225.unknown
_1176901985.unknown
_1176900324.unknown
_1176900574.unknown
_1176899716.unknown
_1176875614.unknown
_1176876245.unknown
_1176876325.unknown
_1176876357.unknown
_1176876287.unknown
_1176876163.unknown
_1176876198.unknown
_1176876059.unknown
_1176874641.unknown
_1176875085.unknown
_1176875359.unknown
_1176874800.unknown
_1176874452.unknown
_1176869406.unknown
_1176869444.unknown
_1176869519.unknown
_1176870891.unknown
_1176870949.unknown
_1176871009.unknown
_1176869593.unknown
_1176870663.unknown
_1176870757.unknown
_1176869516.unknown
_1176869517.unknown
_1176869515.unknown
_1176869514.unknown
_1176869408.unknown
_1176869409.unknown
_1176869411.unknown
_1176869407.unknown
_1176869057.unknown
_1176869211.unknown
_1176869320.unknown
_1176869132.unknown
_1176868666.unknown
_1176868878.unknown
_1176789740.unknown
_1176868344.unknown
_1176868461.unknown
_1176868524.unknown
_1176868407.unknown
_1176790210.unknown
_1176789664.unknown
_1176789691.unknown
_1176789074.unknown
_1176789129.unknown
_1176788960.unknown
_1176469705.unknown
_1176613519.unknown
_1176666225.unknown
_1176697170.unknown
_1176697896.unknown
_1176698062.unknown
_1176697865.unknown
_1176666867.unknown
_1176696326.unknown
_1176666569.unknown
_1176665578.unknown
_1176666153.unknown
_1176666188.unknown
_1176666084.unknown
_1176665091.unknown
_1176665549.unknown
_1176664777.unknown
_1176482861.unknown
_1176554112.unknown
_1176611126.unknown
_1176612332.unknown
_1176554328.unknown
_1176483316.unknown
_1176554041.unknown
_1176483242.unknown
_1176471954.unknown
_1176482022.unknown
_1176482090.unknown
_1176479806.unknown
_1176471649.unknown
_1176471884.unknown
_1176471618.unknown
_1176466869.unknown
_1176468144.unknown
_1176469286.unknown
_1176469360.unknown
_1176469601.unknown
_1176469327.unknown
_1176469214.unknown
_1176469247.unknown
_1176469179.unknown
_1176467547.unknown
_1176467739.unknown
_1176467967.unknown
_1176467613.unknown
_1176467478.unknown
_1176441444.unknown
_1176454745.unknown
_1176456552.unknown
_1176466077.unknown
_1176466244.unknown
_1176456973.unknown
_1176454869.unknown
_1176456081.unknown
_1176454833.unknown
_1176441803.unknown
_1176441853.unknown
_1176454637.unknown
_1176441586.unknown
_1176441587.unknown
_1168084795.unknown
_1168085074.unknown
_1168085114.unknown
_1176441260.unknown
_1176441382.unknown
_1176441383.unknown
_1168085903.unknown
_1168085957.unknown
_1168085964.unknown
_1168085117.unknown
_1168085104.unknown
_1168085110.unknown
_1168085078.unknown
_1168084817.unknown
_1168084825.unknown
_1168084835.unknown
_1168084857.unknown
_1168084867.unknown
_1168084830.unknown
_1168084821.unknown
_1168084806.unknown
_1168084813.unknown
_1168084798.unknown
_1168084090.unknown
_1168084114.unknown
_1168084784.unknown
_1168084790.unknown
_1168084123.unknown
_1168084208.unknown
_1168084780.unknown
_1168084194.unknown
_1168084118.unknown
_1168084104.unknown
_1168084109.unknown
_1168084099.unknown
_1168080889.unknown
_1168083727.unknown
_1168083925.unknown
_1168083717.unknown
_1168080881.unknown
_1168080885.unknown
_1033251839.unknown
_1168080873.unknown
_1033331331.unknown
_1033251827.unknown