随机控制系统仿真.ppt

《随机控制系统仿真.ppt》由会员分享，可在线阅读，更多相关《随机控制系统仿真.ppt（65页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、随机控制系统仿真现在学习的是第1页，共65页按照原则，控制系统可分为确定性系统和随机系统两大类。确定性系统：系统参数、指令、干扰均为已知量随机系统：系统参数、指令、干扰都是随机的随机系统在一定条件一定条件下可以简化为确定性系统。如果系统的一次行为不能不能代表多次运行性能，则该系统不能简化为确定性系统，此时，必须按随机系统进行仿真必须按随机系统进行仿真。多数武器系统为随机系统。引言随机系统仿真的目的是获取状态变量、输出等的统计特性。现在学习的是第2页，共65页7.1概率和随机过程l一、随机控制系统系统模型：若a、b为零阵，x0为已知确定值，则该系统为确定性系统。若a、b、u（t）、x0均为不确定

2、的，则该系统为随机系统。现在学习的是第3页，共65页l二、基本概念1、随机事件：可能发生，可能不发生的事件。2、随机过程：随机系统状态变量和输出的变化过程。这是一组曲线。既含有一次实验中x随t变化的过程，也包含多次实验在相同时刻x的值。3、样本函数：随机过程任一条状态变量或输出变量曲线。xi(t)Li=1是随机过程。xk(t)、x1(t)是样本函数4、随机变量：某时刻的状态变量或输出变量称为随机变量。现在学习的是第4页，共65页随机过程样本函数随机变量现在学习的是第5页，共65页l三、连续随机变量的统计特性1、概率密度p(y)：每个y值发生的可能性。2、概率分布函数：变量值小于或等于y的随机变

3、量的概率。3、均值（一阶矩/数学期望）：E（y）E(y1+y2+yn)=E(y1)+E(y2)+E(yn)现在学习的是第6页，共65页4、均方值E（y2）与均方根sqrt（E(y2)）5、方差2y与方差根y：表示随机变量偏离均值的程度。方差又称二阶中心矩。2y=E(y-E(y)2=E(y2)-E2(y)2(y1+y2+yn)=2(y1)+2(y2)+2(yn)随机系统仿真的目的是获取状态变量、输出等的统计特性。现在学习的是第7页，共65页 用仿真方法获取随机变量的统计特性时，随机变量只有有限个数值，此时，各统计特性的计算公式为：均值：均方值：方差：现在学习的是第8页，共65页l四、典型概率分布

4、1、均匀分布 概率密度：p(y)=1/(b-a)，a=y=b =0 ，其他 概率分布函数：P(y)=0 ，y=a =(y-a)/(b-a)，ayb 均值：E（y）=(a+b)/2，方差：2y=(b-a)2/12。现在学习的是第9页，共65页2、正态分布：自然界中最常见的一种分布。又称高斯分布，钟形分布。2y 为方差，m为均值。概率密度：现在学习的是第10页，共65页中心极限定理：正态概率分布随机变量可以用中心极限定理：正态概率分布随机变量可以用无数个无数个任意分布任意分布的的相互相互独立独立的随机变量合成。的随机变量合成。在在工程工程上，可用上，可用1012个互相独立的个互相独立的均匀分布均匀

5、分布合成一个正态分布。合成一个正态分布。概率分布函数：现在学习的是第11页，共65页l五、随机过程的统计特性平稳随机过程：统计特性不随时间变化。各态历经平稳随机过程：一个样本函数的时间统计特性等于随机过程 统计特性。常见的随机过程都是各态历经的平稳随机过程，故可用样本函数的统计特性来表征随机过程的统计特性。现在学习的是第12页，共65页1、均值：yk(t)为样本，Y(t)为随机过程。2、方差：现在学习的是第13页，共65页3、自相关函数：衡量随机过程功率强弱的尺度。反映Y（t1）与Y（t2）的相关性。RY(t1，t2)大，则Y(t)变化平缓，可预见性大，RY(t1，t2)小，则Y(t)变化剧烈

6、，可预见性小。现在学习的是第14页，共65页4、功率谱密度（能量谱密度）是RY(t1，t2)的傅氏变换，反映Y（t）中不同频率的能量。5、白色随机过程：最简单最简单的随机过程 最理想最理想的随机过程 均值：E（yk）=0，方差：2，自相关函数：Ry()=2()功率谱密度：Sy()=2，各频率分量能量相等，带宽无限。白色随机过程的自相关函数及功率谱密度现在学习的是第15页，共65页l六、随机过程作用下线性系统的响应1、输入输出的时域关系系统模型为脉冲响应函数h（t）y(t)与x(t)为卷积关系均方值：表明输出的均方值与输入的自相关函数有关对白色随机过程：现在学习的是第16页，共65页2、输入输出

7、的频域关系系统模型为传递函数Sy()=|H(j )|2 Sx()若Sx()=2，|H(j )|2不为常值，则y(t)为有色噪声。例：H(j )=1/（j T+1)，Sx()=2 则，y(t)为有色噪声。现在学习的是第17页，共65页白色随机过程作用下线性系统的频域响应现在学习的是第18页，共65页7.2随机控制系统仿真的专门问题l随机数产生l随机数概率分布的转变l随机数概率功率谱密度的转变l系统随机参数误差的产生l系统随机初值的产生l系统随机干扰的产生现在学习的是第19页，共65页一、随机数产生产生伪随机数序列=改变概率分布、功率谱密度=加 入随机系统=进行仿真1、公式：计算机产生伪随机序列，

8、乘同余法公式：x(i+1)=(*x(i+1)=(*x x(i)(mod)(i)(mod)其中：=2k/2，且=8N+/-3，K为字长，N为整数。=2n，且0 2k，n为整数2、特点（1）伪随机序列，周期为/4（2）均匀分布，值域：0 2n-1现在学习的是第20页，共65页3、统计特性（1）概率密度：p（x）=1/，0=x T =2（1-/T）/12 ，|=T现在学习的是第21页，共65页当 较小时，Sx（）近似为常值（5）功率谱密度：在T时，Rx()=0其自相关函数的最大值为：Rx（0）=2T/12=9.6077e+016 0.01时，Rx()=0现在学习的是第25页，共65页乘同余法产生的伪

9、随机序列的功率谱密度功率谱密度：Sx（0）=2*T/12=9.607679205057059e+014低频时，Sx（）近似为常值。T越小，为常值的频带越大。现在学习的是第26页，共65页l二、随机数概率分布的改变多数情况下需要将多数情况下需要将均匀分布改为正态分布。均匀分布改为正态分布。中心极限定理：正态概率分布随机变量可以用中心极限定理：正态概率分布随机变量可以用无数个无数个任意分任意分布布的的相互独立相互独立的随机变量合成。的随机变量合成。在在工程工程上，可用上，可用1012个互相独立的个互相独立的均匀分布均匀分布合成一个正态合成一个正态分布。分布。现在学习的是第27页，共65页1、产生值

10、域为（0，2k-1）的均匀分布随机数序列x1(i)。2、将该序列标称化，值域变为（0，1）方法：x1(i)*1/(2k-1)3、均值移位，值域变为（-0.5，0.5）方法：x1(i)-0.5此时，E（x）=0，2x=1/12 重复1-3步，产生n个不相关的均匀分布随机数序列，记为 x2(i)、x3(i)xn(i)。4、用上述n个不相关的均匀分布随机数序列构成一个正态分布随机数序列。现在学习的是第28页，共65页序列1：x11,x21,xm1 序列2：x12,x22,xm2.序列n：x1n,x2n,xmn一般，取n=12，则合成的随机数序列 E（x）=0，2x=1合成：x1,x2,xn现在学习的

11、是第29页，共65页l三、随机数功率谱密度的改变利用线性系统做成形滤波器 Sy()=|H(j )|2*2T/12，其中，Sy()为需要的功率谱密度，2T/12 为伪随机序列的功率谱密度。则根据Sy()可求出线性系统的频率特性H(j )。现在学习的是第30页，共65页l四、仿真中一些随机参数的产生1、参数误差及随机初值的产生 二者均为正态分布，所给参数为最大允许误差A。任务：根据最大允许误差产生随机参数误差或随机初值序列。现在学习的是第31页，共65页分析分析：由正态分布的性质，一般认为最大允许误差 A=4，为随机序列的方差根。即：随机序列的方差为（取等号情况）：2=（A/4）2 （1）若产生1

12、2个值域为（-b，b）的独立的均匀分布随机序列，其方差均为 2=b2/3则合成后所得正态分布序列的方差为 2=4b2 （2)考虑（1）式，则有 b=A/8即：所产生的均匀分布序列的值域为（即：所产生的均匀分布序列的值域为（-A/8，A/8）现在学习的是第32页，共65页2、随机干扰的产生 随机干扰既要求概率统计特性，又要求频谱特性。一般认为随机干扰的频谱可通过实测或分析得到，即是已知的，故可通过下述方法产生随机干扰。第一步，产生白色随机过程（在一定频带内）。第二步，根据随机干扰的频谱设计成形滤波器。第三步，将白色随机过程作用于成形滤波器，得到有色随机过程。现在学习的是第33页，共65页7.3

13、Monte Carlo法 适用于线性或非线性随机系统仿真，使用时限制条件少，但仿真工作量大。通过对仿真所得数据进行处理得到系统状态变量和输出的统计特性。现在学习的是第34页，共65页蒙特卡罗方法的基本思想蒙特卡罗方法的基本思想l 二十世纪四十年代中期，由于科学技术的发二十世纪四十年代中期，由于科学技术的发展和电子计算机的发明，蒙特卡罗方法作为一种展和电子计算机的发明，蒙特卡罗方法作为一种独立的方法被提出来，并首先在核武器的试验与独立的方法被提出来，并首先在核武器的试验与研制中得到了应用。但其基本思想并非新颖，人研制中得到了应用。但其基本思想并非新颖，人们在生产实践和科学试验中就已发现，并加以利

14、们在生产实践和科学试验中就已发现，并加以利用。用。两个例子 例1.蒲丰氏问题 例2.射击问题（打靶游戏）基本思想现在学习的是第35页，共65页例1.蒲丰氏问题l 为了求得圆周率 值，在十九世纪后期，有很多人作了这样的试验：将长为2 l的一根针任意投到地面上，用针与一组相间距离为2a（la）的平行线相交的频率代替概率P，再利用准确的关系式：l 求出值 其中为投计次数，n为针与平行线相交次数。这就是古典概率论中著名的蒲丰氏问题。现在学习的是第36页，共65页l 一些人进行了实验，其结果列于下表：实验者年份投计次数的实验值沃尔弗(Wolf)185050003.1596斯密思(Smith)185532

15、043.1553福克斯(Fox)189411203.1419拉查里尼(Lazzarini)190134083.1415929现在学习的是第37页，共65页例2.射击问题（打靶游戏）l 设r表示射击运动员的弹着点到靶心的距离，(r)表示击中r处相应的得分数（环数），f(r)为该运动员的弹着点的分布密度函数，它反映运动员的射击水平。该运动员的射击成绩为 l 用概率语言来说，是随机变量(r)的数学期望，即现在学习的是第38页，共65页l 现假设该运动员进行了次射击，每次射击的弹着点依次为r1，r2，rN，则次得分g(r1)，g(r2)，g(rN)的算术平均值l 代表了该运动员的成绩。换言之，为积分的

16、估计值，或近似值。l 在该例中，用次试验所得成绩的算术平均值作为数学期望的估计值（积分近似值）。现在学习的是第39页，共65页基本思想 l 由以上两个例子可以看出，当所求问题的解是某个事件的概率，或者是某个随机变量的数学期望，或者是与概率、数学期望有关的量时，通过某种试验的方法，得出该事件发生的频率，或者该随机变量若干个具体观察值的算术平均值，通过它得到问题的解。这就是蒙特卡罗方法的基本思想。l当随机变量的取值仅为1或0时，它的数学期望就是某个事件的概率。或者说，某种事件的概率也是随机变量（仅取值为1或0）的数学期望。现在学习的是第40页，共65页l 因此，可以通俗地说，蒙特卡罗方法是用随机试

17、验的方法计算积分，即将所要计算的积分看作服从某种分布密度函数f(r)的随机变量(r)的数学期望 l 通过某种试验，得到个观察值r1，r2，rN（用概率语言来说，从分布密度函数f(r)中抽取个子样r1，r2，rN，），将相应的个随机变量的值g(r1)，g(r2)，g(rN)的算术平均值l 作为积分的估计值（近似值）。现在学习的是第41页，共65页l 为了得到具有一定精确度的近似解，所需试验的次数是很多的，通过人工方法作大量的试验相当困难，甚至是不可能的。因此，蒙特卡罗方法的基本思想虽然早已被人们提出，却很少被使用。本世纪四十年代以来，由于电子计算机的出现，使得人们可以通过电子计算机来模拟随机试验

18、过程，把巨大数目的随机试验交由计算机完成，使得蒙特卡罗方法得以广泛地应用，在现代化的科学技术中发挥应有的作用。现在学习的是第42页，共65页l一、蒙特卡罗法仿真步骤1、建立随机系统模型，产生随机参数序列。2、改变随机参数的值进行仿真。多采用数值积分法。一般，每个随机参数都要取不同的值仿真50次以上，以使结果可信。3、处理仿真结果，求状态变量、输出的统计特性。现在学习的是第43页，共65页l二、仿真结果处理设对某参数共仿真N次，每次计算M步。yi1：y11,y21,yM1 yi2：y12,y22,yM2 .yiN：y1n,y2n,yMN现在学习的是第44页，共65页1、均值：2、方差：根据仿真结

19、果求统计特性：均值、方差、概率密度及频谱。现在学习的是第45页，共65页3、概率密度：设yi为yi值域的分组区间，yi=(b a)/LNyij为介于区间(yij-0.5yi，yij+0.5yi)中yij的个数，N为yi中yij的总数，即仿真次数。则，概率密度为：pi(yij)=Nyij/Nyi4、频谱：其中：j为虚单位，r为谐波次数，Ar为幅值。现在学习的是第46页，共65页l例例 dotX(t)=(A+a)X(t)+(B+b)u(t)，y(t)=CX(t)初值：X0+X0X0为初值误差，正态分布，最大误差为A用Monte Carlo法法仿真，求初值误差X0对输出y(t)统计特性的影响。步骤1

20、：产生正态分布的随机数序列X0；用乘同余法产生12个值域为（-A/8，A/8）的均匀分布序列，然后合成即为正态分布的随机数序列X0。步骤2：取不同的X0进行仿真。采用数值积分法，仿真次数在50次以上。步骤3：数据处理，得出初值误差对输出的影响。现在学习的是第47页，共65页仿真流程图现在学习的是第48页，共65页3.蒙特卡罗法的特点蒙特卡罗法的特点优点(1)能够比较逼真地描述具有随机性质的事物的特点及物理实验过程。(2)受几何条件限制小。(3)收敛速度与问题的维数无关。(4)具有同时计算多个方案与多个未知量的能力。(5)误差容易确定。(6)程序结构简单易于实现。缺点(1)收敛速度慢。(2)误差

21、具有概率性。现在学习的是第49页，共65页1)能够比较逼真地描述具有随机性质的事物的特点及物理实验过程l 从这个意义上讲，蒙特卡罗方法可以部分代替物理实验，甚至可以得到物理实验难以得到的结果。用蒙特卡罗方法解决实际问题，可以直接从实际问题本身出发，而不从方程或数学表达式出发。它有直观、形象的特点。现在学习的是第50页，共65页2)受几何条件限制小l 在计算s维空间中的任一区域Ds上的积分l无论区域Ds的形状多么特殊，只要能给出描述Ds的几何特征的条件，就可以从Ds中均匀产生N个 l ，得到积分的近似值。l其中Ds为区域Ds的体积。这是数值方法难以实现的。另外，在具有随机性质的问题中，如考虑的系

22、统形状很复杂，难以用一般数值方法求解，而使用蒙特卡罗方法，不会有原则上的困难。现在学习的是第51页，共65页3)收敛速度与问题的维数无关l 由误差定义可知，在给定置信水平情况下，蒙特卡罗方法的收敛速度为，与问题本身的维数无关。维数的变化，只引起抽样时间及估计量计算时间的变化，不影响误差。也就是说，使用蒙特卡罗方法时，抽取的子样总数N与维数s无关。维数的增加，除了增加相应的计算量外，不影响问题的误差。这一特点，决定了蒙特卡罗方法对多维问题的适应性。而一般数值方法，比如计算定积分时，计算时间随维数的幂次方而增加，而且，由于分点数与维数的幂次方成正比，需占用相当数量的计算机内存，这些都是一般数值方法

23、计算高维积分时难以克服的问题。现在学习的是第52页，共65页4)具有同时计算多个方案与多具有同时计算多个方案与多个未知量的能力个未知量的能力l（1）对于那些需要计算多个方案的问题，使用蒙特卡罗方法有时不需要像常规方法那样逐个计算，而可以同时计算所有的方案，其全部计算量几乎与计算一个方案的计算量相当。例如，对于屏蔽层为均匀介质的平板几何，要计算若干种厚度的穿透概率时，只需计算最厚的一种情况，其他厚度的穿透概率在计算最厚一种情况时稍加处理便可同时得到。l（2）另外，使用蒙特卡罗方法还可以同时得到若干个所求量。例如，在模拟粒子过程中，可以同时得到不同区域的通量、能谱、角分布等，而不像常规方法那样，需

24、要逐一计算所求量。现在学习的是第53页，共65页5)误差容易确定误差容易确定l（1）对于一般计算方法，要给出计算结果与真值的误差并不是一件容易的事情，而蒙特卡罗方法则不然。根据蒙特卡罗方法的误差公式，可以在计算所求量的同时计算出误差。对干很复杂的蒙特卡罗方法计算问题，也是容易确定的。l（2）一般计算方法常存在着有效位数损失问题，而要解决这一问题有时相当困难，蒙特卡罗方法则不存在这一问题。现在学习的是第54页，共65页6)程序结构简单，易于实现程序结构简单，易于实现l 在计算机上进行蒙特卡罗方法计算时，程序结构简单，分块性强，易于实现。现在学习的是第55页，共65页7.4伴随系统仿真法l一、伴随

25、系统设系统的状态方程为dot X(t)=A(t)X(t)+B(t)u(t)y(t)=C(t)X(t)X(0)=X0，0=t=tf，A：n x n，B：n x m，C：r x n则其伴随系统为：现在学习的是第56页，共65页：伴随状态，n x 1：伴随输入，r x 1：伴随输出，m x 1：伴随系统时间自变量现在学习的是第57页，共65页l二、伴随系统结构图 已知原系统的结构图，求其伴随系统的结构图，需进行以下四步。1.将原系统的输入变为脉冲输入/白色随机输入2.将初始条件变为脉冲输入3.将原系统的t用 代替4.原系统信号流向倒转，节点变为相加点，相加点变为节点。现在学习的是第58页，共65页1

26、.将原系统的输入变为脉冲输入/白色随机输入将u（t）变为（t）（1）u(t)=1(t)时现在学习的是第59页，共65页(2)u(t)=Vt 时（3）u(t)=at2/2 时（4）u(t)为有色随机输入时x(t)为白色随机过程，H(jw)为成形滤波器Su()=|H(j )|2 Sx()现在学习的是第60页，共65页2、将初始条件转变为脉冲输入 阶跃输入：K(t)=a t+b K(t)=a()+b3、将原系统的 t 用 代替f(t)=1/a(tf-t)+b f(t)=1/a(t+b)原系统 伴随系统现在学习的是第61页，共65页4、原系统信号流向倒转 节点 相加点，相加点节点节点 相加点：相加点节

27、点：原系统伴随系统现在学习的是第62页，共65页l例、系统结构图转变初值：x1(0),x2(0)原系统结构图原系统状态方程：现在学习的是第63页，共65页伴随系统结构图伴随系统状态方程：现在学习的是第64页，共65页l三、伴随系统法仿真1、伴随系统的性质输入输出位置互换伴随系统在(t)(t=0)作用下各输出端的输出分别对应原系统各输入引起的输出在tF时刻的值。tF为系统运行结束时刻，即仿真结束时刻。设原系统脉冲响应为h(tO,tI)，tI为脉冲加入时刻，tO为输出观察时刻；伴随系统脉冲响应为h*(tO,tI)，tI为脉冲加入时刻，tO为输出观察时刻，则有h*(tF-tI，tF-tO)=h(tO,tI)当tO=tF时，该性质变为h*(tF-tI，0)=h(tF,tI)现在学习的是第65页，共65页