系统辨识 课程论文

系统辨识 课程论文
系统辨识 课程论文

课程论文

系统辨识原理及其在飞行器中的应用

目录

摘要............................................................................................................................... II 第一章引言.. (1)

1.1系统辨识简介 (1)

1.2系统辨识原理 (1)

1.3系统建模 (2)

第二章系统辨识发展 (3)

2.1系统辨识的过去 (3)

2.2系统辨识的现在 (3)

第三章系统辨识在飞行器中的应用 (5)

3.1飞行器参数辨识及其进展 (5)

3.2极大似然估计 (5)

3.3云模型优化 (6)

3.4广义卡尔曼滤波法 (6)

3.5傅里叶变换法 (7)

第四章结束语 (8)

参考文献 (9)

摘要

系统辨识在现代飞行控制系统设计中扮演越来越重要的角色,飞行器模型的在线更新使得人们可以采用更智能的控制方法。飞行控制系统的设计首先需要建立系统的数学模型,由于建模不可避免地存在误差,通过在线系统辨识,可以获得更准确的参数、跟踪参数的变化、检测故障的发生,进而控制系统可以做出相应的调整,保持系统的稳定性和预期性能。然而,从带有噪声和干扰的数据中辨识出有用的信息是非常困难的。对于不稳定飞行器,通常只能在闭环稳定控制情形下进行辨识,这将导致辨识所需的激励信息因反馈作用而减弱。基于计算精度和速度的考虑,在线辨识方法通常以递推方式进行,主要分为时域和频域两大类方法。

本文主要介绍了系统辨识原理及发展,系统建模的方法,并且阐述了系统辨识在飞行器中的应用,着重说明了极大似然估计、云模型优化、广义卡尔曼滤波法、傅里叶变换法。

第一章引言

1.1系统辨识简介

辨识、状态估计和控制理论是现代控制理论三个互相渗透的领域。辨识和状态估计离不开控制理论的支持,控制理论的应用又几乎不能没有辨识和状态估计技术。随着控制过程复杂性的提高,控制理论的应用日益广泛,但其实际应用不能脱离被控对象的数学模型。然而在大多数情况下,被控对象的数学模型是不知道的,或者在正常运行期间模型的参数可能发生变化,因此利用控制理论去解决实际问题时,首先需要建立被控对象的数学模型,因此建模在控制器设计中有着广泛的应用,是设计控制器首要解决的问题。系统辨识正是适应这一需要而形成的,他是现代控制理论中一个很活跃的分支。

系统辨识是建模的一种方法,不同的学科领域,对应着不同的数学模型。从某种意义上来说,不同学科的发展过程就是建立他的数学模型的过程。L.A.Zadeh 给出辨识这样的定义:“辨识就是在输入和输出数据的基础上,从一组给定的模型类中,确定一个与所测系统等价的模型。”当然按照这个定义,寻找一个与实际过程完全等价的模型无疑是非常困难的。而从实用性观点出发,对模型的要求并非如此苛刻,为此,对辨识又有一些实用性的定义。比如,P.E.ykhoff给出辨识的定义为:“辨识问题可以归结为用一个模型来表示客观系统(或将要构造的系统)本质特征的一种演算,并用这个模型把对客观系统的理解表示成有用的形式。”总而言之,辨识的实质就是从一组模型类中选择一个模型,按照某种准则,使之能最好地拟合所关心的实际过程的静态或动态特性。

1.2系统辨识原理

系统辨识是在已知或测得系统输入和输出数据的基础上,从一组给定的模型类中,确定一个与所测系统等价的模型,简单点说,就是根据系统的输入输出时间函数来确定描述系统行为的数学模型。系统辨识要素为:

数据:指系统过程的输入数据和输出数据,它是辨识的基础;

模型类:指各种已知的系统过程模型集合,它是辨识时寻找模型的范围;

等价准则:指系统行为相似性、效用等同性的识别标准,它是辨识优化的目标。

辨识的实质就是按某种准则,从一组已知模型类中选择一个模型,使之能最

好地拟合实际过程的动态特性。观测数据含有噪声,因此辨识建模实际上是一种实验统计的方法,所获得的模型只是与实际过程的外特性等价的一种近似描述。

1.3系统建模

建模目的:

(1)估计具有特定物理意义的参数:有些表征系统行为的重要参数是难以直接测量的,例如在生理、生态、环境、经济等系统中就常有这种情况。这就需要通过能观测到的输入输出数据,用辨识的方法去估计那些参数。

(2)仿真:仿真的核心是要建立一个能模仿真实系统行为的模型。用于系统分析的仿真模型不仅要求能真实反映系统的特性,还要求设计参数能正确地符合它本身的物理意义。

(3)预测:这是辨识的一个重要应用方面,其目的是用系统的可测量的输入和输出去预测系统输出的未来的演变。预测模型辨识的等价准则主要是使预测误差平方和最小,对模型的结构及参数则很少再有其他要求。

(4)控制器设计:建立数学模型的最终目的还是设计一个好的控制器。

建模基础是参数设计,必须合理选择输入信号、采样时间、辨识时间、开环或闭环辨识、离线或在线辨识等参数或方式。而模型结构辨识是建模的前提,必须明确模型的基本构型,如动态或静态、离散或连续、线性或非线性等模式。线性离散模型是最简单的辨识模型。线性离散模型仅含单个或数个变量,可表示成另外一些变量在时间或空间的离散点上的线性组合。

随着系统辨识新理论、新方法的不断发展,以及与其他学科的相互渗透交融,系统辨识与建模理论日趋成熟,被广泛应用于国民经济和国防建设的各个工程技术领域,包括航天、航空、机器人领域,各种工业过程及社会经济与管理等内容。系统辨识与建模理论已成为众多先进控制技术(如神经网络、模糊控制及遗传算法等)的研究基础,同时,这些技术反过来又推动了该理论体系、内容以及技术的快速更新。

第二章系统辨识发展

近二十年来,系统辨识获得了长足的发展,已经成为控制理论的一个十分活跃而又重要的分支。从线性现象和线性系统的研究过渡到非线性现象和非线性系统的研究是科学发展的必然结果,这不仅是对科学家们一种新的挑战,而且也是人类社会向更高级形式演化的一种必然。随着智能控制理论、遗传算法理论等的不断成熟,逐渐形成了形式多样的现代的系统辨识方法,并且已在实际问题应用中取得了较好的使用效果。

2.1系统辨识的过去

系统辨识方法有经典和现代之分。经典的系统辨识方法包括阶跃响应法、脉冲响应法、频率响应法、相关分析法、谱分析法、最小二乘法和极大似然法等。其中最小二乘法是一种经典的和最基本的,也是应用最广泛的方法。但是,最小二乘估计是非一致的,是有偏差的,所以为了克服他的缺陷,而形成了一些以最小二乘法为基础的诸如广义最小二乘法的系统辨识方法。

虽然最小二乘法能适应很多的系统,但现实中的实际系统很多都是具有不确定性的复杂系统。对于这类系统,经典的辨识建模方法难以得到较好的控制结果,或者说,经典的辨识建模方法存在一定的缺陷。比如:最小二乘法要求输入信号已知,并且必须具有较明显的变化,然而,在某些动态系统中,系统的输入条件往往无法保证,此时就需要寻求新的方法了。

2.2系统辨识的现在

复杂系统分析的难点在于线性和非线性特性的描述,如何建立有效的模型进行真实地评价和高效地优化,进而开展系统全寿命周期的可行性研究、设计开发、测试完善等任务一直是比较难以突破的节点。随着系统辨识的发展,非线性系统建模已从用线性模型逼近发展到用非线性模型逼近的阶段,具体的现代的系统辨识方法有:集员系统辨识法、多层递阶系统辨识法、神经网络系统辨识法、遗传算法系统辨识法、模糊逻辑系统辨识法以及小波网络系统辨识法。

近十年来,随着小波分析理论的发展与成熟,小波网络作为一种有突出特点的前向神经网络受到越来越多的关注和重视。小波网络具有相对有效和简洁的建

模方法(平移和伸缩小波母小波),能够构成框架、紧框架,甚至正交基,构造效率高,收敛速度快,并能解决一般的“维数灾”问题,逼近单变量函数的渐进最优逼近器已经被大量应用于系统辨识中。在系统辨识中,尤其在非线性系统辨识中的应用潜力越来越大,为不确定的复杂的非线性系统辨识提供了一种新的有效途径,其具有良好的应用前景。

图1系统辨识原理图

第三章系统辨识在飞行器中的应用

3.1飞行器参数辨识及其进展

飞行器参数辨识是动力学系统辨识研究中起步早,取得成绩最为显著的一个分支,其研究的侧重点是对飞行器气动参数(或流体动力参数)进行辨识,包括气动力数学模型结构及其中未知参数的辨识。这是一门介于空气(流体)动力学,飞行(航行)力学,弹性力学和现代控制理论间的交叉学科。自二十世纪20年代初Warner等的先导性工作以来,已有70多年的研究历史。

早期的气动参数辨识是通过飞行器稳定飞行数据得到飞行特性,而今已能采用随意大幅值机动飞行试验数据辨识其气动,控制导数。40一50年代,通过飞行试验数据确定气动参数还主要是采用频率法、回归技术、时间矢量法和模拟匹配技术。随着系统辨识学科的形成和在飞行试验数据分析中的应用,使飞行器气动参数辨识研究于70年代后得到了飞速发展,并在飞行器设计中起到越来越重要的作用。被应用于验证和校正飞行器气动参数风洞实验和理论计算结果;为飞行器控制系统设计和改进提供基本依据;鉴定定型飞行器的飞行品质;研究高性能飞行器的飞行品质;对失事飞行器事故进行分析;研制飞行器自适应控制系统等。飞行器气动或流体动力参数辨识的主要内容有以下八个方面:模型辨识、参数估计、数据预处理和相容性检验、试验设计与最佳输入选择、弹性与非定常效应、频域辨识、闭环辨识和系统验证。

3.2极大似然估计

在飞行器气动参数辨识研究中,由于极大似然估计具有良好的收敛特性而被认为是最有效的估计方法之一。极大似然估计的主要问题是似然函数的优化计算,常用的算法是修正牛顿一拉夫森方法。它需要预先确定气动模型结构以推导灵敏度方程。进而积分求解灵敏度方程获得灵敏度矩阵。这对于由非线性连续一离散系统描述的飞行器状态方程和观测方程,其积分计算灵敏度方程数目是n*p个(n 是状态维数,P是参数维数)。当待估参数较多时, 灵敏度方程可达上百个,无论是方程推导还是积分计算都是相当繁杂的。为此,Murphy提出了通过曲面拟合近似计算灵敏度的思想川,避开了灵敏度方程的推导与计算。但由于参数初始摄动量选择、灵敏度初值计算与递推计算等问题没有解决,影响了此算法的实际应

用。

3.3云模型优化

极大似然估计方法在飞行器参数辨识中得到了广泛应用,该方法需要预先推导灵敏度方程,进而求解灵敏度矩阵,在应用过程中比较繁杂,且容易陷入局部最优。提出一种基于云模型优化的飞行器参数辨识算法,根据极大似然估计原理,利用云模型的优化理论对极大似然函数进行优化,从而得到待辨识参数值。该算法不必推导灵敏度矩阵,对初值要求不高,应用便捷,且保留了云模型优化的特点,收敛速度较快、不易陷入局部最优。

云模型是李德毅院士提出的一种定性定量转化模型,已经在智能控制、模糊评测等多个领域得到应用。云模型在知识表达时具有不确定中带有确定性、稳定之中又有变化的特点。根据其转化过程中的优良特性,结合遗传算法的基本思想,产生的基于云模型的全局最优化算法易于实现,不存在编码问题,具有精度高、收敛速度快等优点。

3.4广义卡尔曼滤波法

其基本思想乃是使观测量出现概率最高的模型是最优模型,采用滤波算法。具体算法是:对给定的汉选模型集与输入,逐个或同时采用广义卡尔曼滤波迸行增广状态、新息和协方差的估计,并用来计算观测量出现的概率,定义此概率为该模型到该观测点为止所获得的可洁度水平。可信度水平为最高者为最优模型。滤波过程中对可信度小的模型要及时剔除;可信度明显优于其他模型者,可及时定为最优模型。

还有提出以残差协方差和重复测量的频率响应曲线协方差之比为判据,以自由度的残差平方和为判据等等。采用F检验与其他判据进行建模仿真试验,结果是几种判据相差不多。也有人得出不同结论,认为一般子集回归的F检验法不能保证建立的模型是正确的,并给出了F检验法所得模型很不正确的仿真试验例子,作者提出采用特定区间中的正交—正则多项式做为候选模型项,并采用特征值分解法或矩阵变换将实测数据变换到主元基系,然后采用子集回归法确定模型的方法。其仿真计算结果表明,这种算法有较好的稳健性,即使模型数据不完备或有较大数值计算误差也能得到较好的结果。

3.5傅里叶变换法

在建立飞行器系统模型结构的基础上,利用频域递推傅里叶变换及最小二乘方法,实现对气动及控制偏导数的在线辨识。

Delft两步法:其基本思想是先辨识出系统状态,通过迭代扩展卡尔曼滤波器,消除传感器噪声、偏差以及风的影响,确定真实的飞行器状态;然后采用递推最小二乘方法辨识出气动参数、推力、力矩以及气动导数。通过辨识得到了准确的对象模型,使得对模型误差极其敏感的自适应动态逆控制方法可以应用于实际,该技术已在试验研究飞机上得到了初步应用。Delft两步法解决了时域辨识方法不能很好处理噪声的难题。该方法辨识飞行器结构上的故障非常有效,但对于执行机构的故障辨识仍需发展有效的算法。并且,由于时域最小二乘方法的局限性,参数辨识的误差特别是收敛过程中的偏差不能可靠地给出。

傅里叶变换回归辨识方法:其基本思想是对飞行器模态进行充分的激励,采用有限递推傅里叶变换RFT得到系统模型参数的最小二乘解。频域辨识方法假定参数辨识模型是线性的且结构形式在辨识过程中固定不变。对于非线性参数模型,该方法仍然适用,需要先用时变线性模型进行等效转换。

第四章结束语

通过研究有关的文献可以看出,辨识研究要回答的问题是:模型辨识的趋势是什么,系统辨识模型的最终用途在那里。许多文献表明系统辨识已被用于研究和教学,或者是实际用于工程,比如仿真数据库,控制器设计,飞行器仿真,用于加快证明过程,或是为飞行器进行包线扩展等。涉及技术既包括时域也包括频域,应用范围从测量数据校准到空气动力学的数据库建立,从固定翼的飞行器、航天飞机、直升机到无人驾驶飞行器。其中许多研究人员指出如何从飞行数据辨识出高精度模型是当今面对的一个巨大挑战。

实践证明,应用的多样化和普及水平在二十年前是很难想象的。可以明确地讲,在过去的十年里,无论应用的数量,还是建模的复杂度上都有明显的提高。因此,企图用一本手册推导出高精度模型是困难的,也不现实,但是可以归纳出一定的规范和指导准则。

此外,实时参数估计将作为当前和未来研究的另一领域重新引起重视。一些理论将聚焦于模型的不确定性估计。其它方面进展将从实际需求中产生,例如追求更好的空气数据测定方法。不久的将来,另一个进展有可能发生在系统辨识和计算流体动力学的综合集成方法上。这将要求有新的模型概念和相关联的测量技术以及实验技术来进行参数估计。我们将面临更苛刻和更吸引人的建模问题。但可以相信,系统辨识还将提供最根本的解决问题方案来揭示建模过程中的缺陷,降低开发风险和提高飞行安全保证。

参考文献

[1]崔平远.递推计算灵敏度的极大似然估计算法实现.自动化学报,1995,21(1)

82一86.

[2]付斌,李道国,王慕快.云模型研究的回顾与展望[J].计算机应用研究,2011,28(2):420-426.

[3]李德毅,刘常昱.论正态云模型的普适性[J].中国工程科学,2004,6(8):28-35.

[4]蔡金狮.飞行器系统辨识学[M].北京:国防工业出版社,2003:135-139.

[5]Suk J,Lee Y,Kim S,et al.System identification and stability evaluation of an unmanned aerial vehicle from automated flight tests[J].KSME International Journal,2003,17(5):654-667.

[6]Morelli E A.Flight test maneuvers for efficient aerodynamic modeling,AIAA-2011-6672[R].Reston:AIAA,2011.

[7]王文正,蔡金狮.行器动参数的集员辨识[J].宇航学报,1998,19(2):31-36.

[8]Jwo D J,Wang S H.Adaptive fuzzy strong trackingextended kalman filtering for GPS navigation[J].IEEESensors Journal, 2007,7(5):778-789.

[9]于黎明,王占林,裘丽华.飞机操纵系统建模及降阶仿真研究.计算机仿真,2000,17(3):15-186.

[10]吴建德.基于频域辨识的微小型无人直升机的建模与控制研究.杭州:浙江大学博士学位论文,2007.

系统辨识大作业1201张青

《系统辨识》大作业 学号:12051124 班级:自动化1班 姓名:张青 信息与控制工程学院自动化系 2015-07-11

第一题 模仿index2,搭建对象,由相关分析法,获得脉冲响应序列?()g k ,由? ()g k ,参照讲义, 获得系统的脉冲传递函数()G z 和传递函数()G s ;应用最小二乘辨识,获得脉冲响应序列? ()g k ;同图显示两种方法的辨识效果图;应用相关最小二乘法,拟合对象的差分方程模型;构建时变对象,用最小二乘法和带遗忘因子的最小二乘法,(可以用辨识工具箱) 辨识模型的参数,比较两种方法的辨识效果差异; 答:根据index2搭建结构框图: 相关分析法:利用结构框图得到UY 和tout 其中的U 就是题目中要求得出的M 序列,根据结构框图可知序列的周期是 1512124=-=-=n p N 。 在command window 中输入下列指令,既可以得到脉冲相应序列()g k :

aa=5;NNPP=15;ts=2; RR=ones(15)+eye(15); for i=15:-1:1 UU(16-i,:)=UY(16+i:30+i,1)'; end YY=[UY(31:45,2)]; GG=RR*UU*YY/[aa*aa*(NNPP+1)*ts]; plot(0:2:29,GG) hold on stem(0:2:29,GG,'filled') Grid;title('脉冲序列g(τ)') 最小二乘法建模的响应序列 由于是二阶水箱系统,可以假设系统的传递函数为2 21101)(s a s a s b b s G +++= ,已知)(τg ,求2110,,,a a b b

系统辨识实验1实验报告

实验报告 --实验1.基于matlab的4阶系统辨识实验 课程:系统辨识 题目:基于matlab的4阶系统辨识实验 作者: 专业:自动化 学号:11351014 目录 实验报告 (1) 1.引言 (2) 2.实验方法和步骤 (2) 3.实验数据和结果 (2) 4.实验分析 (4)

1、 引言 系统辨识是研究如何确定系统的数学模型及其参数的理论。而模型化是进行系统分析、仿真、设计、预测、控制和决策的前提和基础。 本次实验利用matlab 工具对一个简单的4阶系统进行辨识,以此熟悉系统辨识的基本步骤,和matlab 里的一些系统辨识常用工具箱和函数。 这次实验所采取的基本方法是对系统输入两个特定的激励信号,分别反映系统的动态特性和稳态特性。通过对输入和输出两个系统信号的比较,来验证系统的正确性。 2、 实验方法和步骤 2.1 实验方法 利用matlab 对一个系统进行辨识,选取的输入信号必须能够反映系统的动态和稳态两个方面的特性,才能更好地确定系统的参数。本次实验采取了两种输入信号,为反映动态特性,第一个选的是正弦扫频信号,由下面公式产生: 选定频率范围 ,w(t)是时间t 的线性函数,具有扫频性质,可以反映系统的动态特性。 为反映稳态特性,选的输入信号是阶跃信号。以上的到两组数据,利用matlab 的merge()函数,对两组数据融合,然后用matlab 系统辨识工具箱中的基于子空间方法的状态空间模型辨识函数n4sid()来对系统进行辨识 2.2 实验步骤 (1)建立一个4阶的线性系统,作为被辨识的系统,传递函数为 3243211548765 ()125410865 s s s G s s s s s -+-+=++++ (2)产生扫频信号u1和阶跃信号u2 (3)u1、u2作为输入对系统进行激励,分别产生输出y1和y2 (4)画出稳态测试输入信号u1-t 的曲线,和y1-t 的曲线 画出动态测试输入信号u2-t 的曲线,和y2-t 的曲线 (5)使用merge()函数对u1-y1数据和u2-y2数据进行融合,并使用n4sid()函数对系统进行辨识。 (6)画出原系统和辨识出的系统的零极点图,画出原系统和辨识出的系统的阶跃响应特性曲线,通过对比,验证辨识出的系统的准确性。 3、 实验数据和结果 (1) 分别以扫频正弦函数、阶跃函数作为系统的激励,得到的输出:

系统辨识之经典辨识法

系统辨识作业一 学院信息科学与工程学院专业控制科学与工程 班级控制二班 姓名 学号

2018 年 11 月 系统辨识 所谓辨识就是通过测取研究对象在认为输入作用的输出响应,或正常运行时 的输入输出数据记录,加以必要的数据处理和数学计算,估计出对象的数学模型。 辨识的内容主要包括四个方面: ①实验设计; ②模型结构辨识; ③模型参数辨识; ④模型检验。 辨识的一般步骤:根据辨识目的,利用先验知识,初步确定模型结构;采集 数据;然后进行模型参数和结构辨识;最终验证获得的最终模型。 根据辨识方法所涉及的模型形式来说,辨识方法可以分为两类:一类是非参 数模型辨识方法,另一类是参数模型辨识方法。 其中,非参数模型辨识方法又称为经典的辨识方法,它主要获得的是模型是 非参数模型。在假定过程是线性的前提下,不必事先确定模型的具体结构,广泛 适用于一些复杂的过程。经典辨识方法有很多,其中包括阶跃响应法、脉冲响应法、相关分析法和普分析法等等,本次实验所采用的辨识方法为阶跃响应法和脉 冲响应法。 1.阶跃响应法 阶跃响应法是一种常用非参数模型辨识方法。常用的方法有近似法、半对数法、切线法、两点法和面积法等。本次作业采用面积法求传递函数。 1.1面积法 ① 当系统的传递函数无零点时,即系统传递函数如下: G(S) = + ?11?1+?+ 1+1 (1-1) 系统的传递函数与微分方程存在着一一对应的关系,因此,可以通过求取 微分方程的系数来辨识系统的传递函数。在求得系统的放大倍数K后,要得到无 因次阶跃响应y(t)(设τ=0),其中y(t)用下式描述: () ?1 () (1-2) 面积法原则上可以求出n为任意阶的个系数。以n为3为例。有: 3() 2() () {| →∞ =| →∞ =| →∞ = 0 (1-3) ()| →∞ = 1

最优估计大作业1.

最优估计大作业 姓名:李海宝 学号:S314040186 导师:刘胜 专业:控制科学与工程

模糊逻辑卡尔曼滤波器在智能AUV导航系统中的自适应调 整 摘要 本论文基于全球定位系统(GPS)和几个惯性导航系统(INS)传感器描述了对于自主水下航行器(AUV)应用的一种智能导航系统的执行过程。本论文建议将简单卡尔曼滤波器(SKF)和扩展卡尔曼滤波器(EKF)一前一后地用于融合INS 传感器的数据并将它们与GPS数据结合到一起。传感器噪声特性里潜在的变化会引起SKF和EKF的初始统计假定的调整,本论文针对这一问题着重突出了模糊逻辑方法的使用。当这种算法包含实际传感器噪特性的时候,SKF和EKF只能维持他们的稳定性和性能,因此我们认为这种自适应机制同SKF与EKF一样有必要。此外,在提高导航系统的可靠性融合过程期间,故障检测和信号恢复算法也需在此要讨论。本论文建议的这种算法用于使真实的实验数据生效,这些数据都是从Plymouth大学和Cranfield大学所做的一系列AUV实验(运行低成本的锤头式AUV)中获得的。 关键词:自主水下航行器;导航;传感器融合;卡尔曼滤波器;扩展卡尔曼滤波器;模糊逻辑 1.引言 对于以科学、军事、商业为目的应用,如海洋勘察、搜索未爆弹药和电缆跟踪检查,AUV的发展需要相应导航系统的发展。这样的系统提供航行器位置和姿态的数据是很有必要的。在这样的系统中对精度的要求是最重要的:错误的位置和姿态数据会导致收集数据的一个毫无意义的解释,或者甚至AUV的一个灾难性故障。 越来越多来自整个世界的研究团队正利用INS和GPS来研发组合导航系统。然而,他们的工作中几乎都没有明确几个INS传感器融合的本质要求,这些传感器用于确保用户保持精度或甚至用来防止在与GPS融合之前导航系统这部分的完全失败。例如,金赛和惠特科姆(2003)使用一个切换机制来防止INS的完全失败。虽然这个方法简单易行,但是可能不适合用于维持一个确定的精度等级。 出于多传感器数据融合和集成的目的,几种估计方法在过去就已经被使用过。为此,SKF/EKF和它们的变形在过去就已经是流行的方法,并且一直到现在都对开发算法感兴趣。然而,在设计SKF/EKF过程中,一个显著的困难经常会被

系统辨识实验报告30288

一、相关分析法 (1)实验原理 图1 实验原理图 本实验的原理图如图1。过程传递函数()G s 中12120,8.3, 6.2K T Sec T Sec ===;输入变量()u k ,输出变量()z k ,噪声服从2(0,)v N σ,0()g k 为过程的脉冲响应理论 值,?()g k 为过程脉冲响应估计值,()g k 为过程脉冲响应估计误差。 过程输入()u k 采用M 序列,其输出数据加白噪声()v k 得到输出数据()z k 。利 用相关分析法估计出过程的脉冲响应值?()g k ,并与过程脉冲响应理论值0()g k 比较,得到过程脉冲响应估计误差值()g k 。 M 序列阶次选择说明:首先粗略估计系统的过渡过程时间T S (通过简单阶跃响应)、截止频率f M (给系统施加不同周期的正弦信号或方波信号,观察输出)。本次为验证试验,已知系统模型,经计算Hz T T f M 14.01 2 1≈= ,s T S 30≈。根据式M f t 3 .0≤ ?及式S T t N ≥?-)1(,则t ?取值为1,此时31≥N ,由于t ?与N 选择时要求完全覆盖,则选择六阶M 移位寄存器,即N =63。

(2)编程说明 图2 程序流程图 (3)分步说明 ① 生成M 序列: M 序列的循环周期63126=-=N ,时钟节拍1t Sec ?=,幅度1a =,移位寄存器中第5、6位的内容按“模二相加”,反馈到第一位作为输入。其中初始数据设为{1,0,1,0,0,0}。程序如下:

② 生成白噪声序列: 程序如下: ③ 过程仿真得到输出数据: 如图2所示的过程传递函数串联,可以写成形如1212 11 ()1/1/K G s TT s T s T = ++, 其中112 K K TT = 。 图2 过程仿真方框图 程序如下: ④ 计算脉冲响应估计值:

系统辨识作业2

系统辨识作业 学院: 专业: 姓名: 学号: 日期:

系统辨识作业: 以下图为仿真对象 图中,v(k)为服从N(0,1)正态分布的不相关随即噪声,输入信号采用循环周期Np>500的逆M 序列,幅值为1,选择辨识模型为: )()2()1()2()1()(2121k v k u b k u b k z a k z a k z +-+-=-+-+ 加权因子1)(=Λk ,数据长度L=500,初始条件取I P 610)0(= ,????????? ???=001.0001.0001.0)0(? θ 要求:(1)采用一次完成最小二乘法对系统进行辨识,给出数据u(k)和z(k), 及L H ,L Z 和θ 和)?(θ J 的值。 (2)采用递推最小二乘法进行辨识,要给出参数收敛曲线以及新息)(~k Z ,残差)(k ε,准则函数)(k J 随着递推次数K 的变化曲线。 (3)对仿真对象和辨识出的模型进行阶跃响应对比分析以检验辨识结果的实效。 1、一次完成法对系统进行辨识: 估计L T L L T L LS Z H H H 1)(?-=θ ,其中 []2121,,,b b a a LS =θ ????? ? ??????=L L Z Z Z Z 21 ????????????------------=????????? ???=)2()1()2()1()0()1()0()1()1()0()1() 0()()2()1(L u L u L z L z u u z z u u z z L h h h H L 一次完成算法对系统辨识的Matlab 程序见附录: 部分输入、输出数据如下,全部的输入输出数据用图1.1所示 输入数据u(k)=Columns 1 through 16 0 0 1 1 1 1 0 0 0 0 1 0 0 1 1 0

系统辨识与自适应控制硕士研究生必修课程考核

《系统辨识与自适应控制》硕士研究生必修课程考核(检测技术与自动化装置专业)2003.5. 22 可下载自https://www.360docs.net/doc/c710045064.html,/xuan/leader/mrj/ 学生姓名:考核成绩: 一、笔试部分 (占课程成绩的 80% ) 考试形式:笔试开卷 答卷要求:笔答,可以参阅书籍,要求简明扼要,不得大段抄教材,不得相互抄袭 试题: 1 简述系统辨识的基本概念(概念、定义和主要步骤)(10分) 2 简述相关辨识的基本原理和基于二进制伪随机序列的相关辩识方法(原理、 框图、特点)。(10分) 3 简述离散线性动态(SI / SO)过程参数估计最小二乘方法(LS法)的主要 内容和优缺点。带遗忘因子递推最小二乘估计(RLS法)的计算步骤和主要递推算式的物理意义(10分) 4 简述什么是时间序列?时间序列建模如何消除恒定趋势、线性趋势和季节性 的影响?(10分) 5 何谓闭环系统的可辨识性问题,它有那些主要结论?(10分) 6 何谓时间离散动态分数时滞过程?“分数时滞”对过程模型的零点和极点有 什么影响?(10分) 7 简述什么是自适应控制,什么是模型参考自适应控制(MRAC)?,试举一例说明MRAC的设计方法(10分)。 8 请设计以下过程( yr = 0 ) y(k) -1.6y(k-1)+0.8y(k-2) = u(k-2)- 0.5u(k-3)+ε(k)+1.5ε(k-1)+0.9ε(k-2) 的最小方差控制器(MVC)和广义最小方差控制器(GMVC), 并分析他们的主要性能。(10分) 二、上机报告RLS仿真(占课程成绩的 20%) 交卷时间:6月9日下午

系统辨识大作业论文Use

中南大学 系统辨识大作业 学院:信息科学与工程学院 专业:控制科学与工程 学生姓名:龚晓辉 学号:134611066 指导老师:韩华教授 完成时间:2014年6月

基于随机逼近算法的系统辨识设计 龚晓辉1, 2 1. 中南大学信息科学与工程学院,长沙410083 2. 轨道交通安全运行控制与通信研究所, 长沙410083 E-mail: csugxh@https://www.360docs.net/doc/c710045064.html, 摘要:本文对系统辨识的基本原理和要素进行了详细阐述,介绍和分析了系统辨识中常用的最小二乘算法,极大似然法,神经网络算法和随机逼近算法。随机逼近算法只需利用输入输出的观测来辨识系统参数,在实际中有重要运用。本文对随机逼近算法进行了详细说明。同时,针对一个三阶系统设计了KW随机逼近算法进行了参数辨识,并且和递推最小二乘法进行了对比。实验证明在实际辨识过程中两种算法各有优缺点。 关键词: 系统辨识, 随机逼近法, 递推最小二乘法 1.引言 在我们所学的线性系统理论中,都是在系统模型已知的情况来设计控制率,使系统达到稳定性,准确性和快速性的要求。然而,在实际系统中,对象的模型往往是未知的。而且,非线性是普遍存在的,线性系统只是对非线性系统的一种近似。因此,了解对象准确的模型,对设计控制器及其重要。在一些实际对象中,如导弹,化学过程,生物规律,药物反应,以及社会经济等,这些对象使用机理分析法比较困难,但是通过使用辨识技术可以建立系统精确的模型,确定最优控制率[1]。如今,系统辨识技术已经在航空航天,海洋工程,生物学等各个领域获得了广泛运用。 2.系统辨识的基本思想与常用方法 辨识的目的是为了获得对象模型。对象的模型有多种表现形式,它包括直觉模型,图表模型,数学模型,解析模型,程序模型和语言模型。这些模型之间可以相互转换。我们在建立系统模型时,需要遵循目的性,实在性,可辨识性,悭吝性的基本原则。目的性指的是建模的目的要明确,实在性指的是模型的物理概念要明确。可辨识性指的是模型结构合理,输入信号持续激励,数据量充足。悭吝性指的是被辨识参数的个数要尽量少。 辨识对象模型要遵循上面的基本原则。它是将对象看成一个黑箱。从含有噪声的输入输出数据中,按照一个准则,运用辨识理论,从一组给定的模型中,确定一个与所测系统等价的模型,是现代控制理论的一个分支。系统辨识由数据、模型类和准则三要素组成。数据是由观测实体而得,它不是唯一的,受观测时间、观测目的、观测手段等影响。模型类就是模型结构,它也不是唯一的,受辨识目的、辨识方法等影响。而准则是辨识的优化目标,用来衡量模型接近实际系统的标准。它也不是唯一的,受辨识目的、辨识方法的影响。由于存在多种数据拟合

系统辨识实验报告

南京理工大学 电加热炉动态特性辨识实验报告 作者: 张志鹏(94)学号:813001010014 实验时间2013年11月24日 组员: 刘心刚(63)李昊(88)倪镭(90) 任课老师:郭毓教授 2013 年 11 月

1.熟悉对实际控制系统的辨识与参数估计,并利用所得模型进行控制仿真,进而控制实际系统。 2.掌握实际工程中常用的辨识方法,如LS,RLS,RLES等。 二、实验平台: 嵌入式温度控制系统主要由嵌入式温度控制器、立式RGL-9076A 型温箱、NETGEAR 无线路由器和24V 开关电源等组成。系统电气连接如图1 所示。系 统采用CS(客户端—服务器)模式实现了一对一的服务器、客户端的数据通信。 嵌入式控制系统软软硬件运行平台. 硬件:PC 机、嵌入式温度控制器、NETGEAR 无线路由器等。 软件:Windows XP、Microsoft Visual C++ 6.0、Matlab 2007a 等。 图1 实验硬件平台

1.设置硬件。根据实验手册上的连接方式,确认硬件连接是否正确。根据使用手册进行IP设置、系统参数设置,直至软件可以实时显示温度曲线。 2.达到稳态。我们首先采用81V的加热电压加热使系统尽快到达某一较稳定温度。使用3S的采样周期进行采样温度信号。当温箱实际温度达到135度左右时,温度变化曲线几乎持平,我们认定此时温箱系统处于稳态。 3.加入辨识信号。这里选选取M序列进行辨识,在试验阶段我们组做了一组数据:选取M序列幅值为+20,-20,,辨识信号的采样周期为40s。加入辨识信号后继续进行数据采集。 4.数据处理、辨识系统模型。 5.分析辨识结果得出结论。 四、辨识算法及过程 经过分析研究,确定使用计算残差平方和的RELS方法验证模型的阶次及延时并辨识系统模型参数。 1、确定系统的延迟d

现代控制理论课程报告

现代控制理论课程总结 学习心得 从经典控制论发展到现代控制论,是人类对控制技术认识上的一次飞跃。现代控制论是用状态空间方法表示,概念抽象,不易掌握。对于《现代控制理论》这门课程,在刚拿到课本的时候,没上张老师的课之前,咋一看,会认为开课的内容会是上学期学的控制理论基础的累赘或者简单的重复,更甚至我还以为是线性代数的复现呢!根本没有和现代控制论联系到一起。但后面随着老师讲课的风格的深入浅出,循循善诱,发现和自己想象的恰恰相反,张老师以她特有的讲课风格,精心准备的ppt 课件,向我们展示了现代控制理论发展过程,以及该掌握内容的方方面面,个人觉得,我们不仅掌握了现代控制理论的理论知识,更重要的是学会了掌握这门知识的严谨的逻辑思维和科学的学习方法,对以后学习其他知识及在工作上的需要大有裨益,总之学习了这门课让我受益匪浅。 由于我们学习这门课的课时不是很多,并结合我们学生学习的需求及所要掌握的课程深入程度,张老师根据我们教学安排需要,我们这学期学习的内容主要有:1.绪论;2.控制系统的状态表达式;3.控制系统状态表达式的解;4.线性系统的能空性和能观性;5.线性定常系统的综合。而状态变量和状态空间表达式、状态转移矩阵、系统的能控性与能观性以及线性定常系统的综合是本门课程的主要学习内容。当然学习的内容还包括老师根据多年教学经验及对该学科的研究的一些深入见解。 在现代科学技术飞速发展中,伴随着学科的高度分化和高度综合,各学科之间相互交叉、相互渗透,出现了横向科学。作为跨接于自然科学和社会科学的具有横向科学特点的现代控制理论已成为我国理工科大学高年级的必修课。 经典控制理论的特点 经典控制理论以拉氏变换为数学工具,以单输入-单输出的线性定常系统为主要的研究对象。将描述系统的微分方程或差分方程变换到复数域中,得到系统的传递函数,并以此作为基础在频率域中对系统进行分析和设计,确定控制器的结构和参数。通常是采用反馈控制,构成所谓闭环控制系统。经典控制理论具有明显的局限性,突出的是难以有效地应用于时变系统、多变量系统,也难以揭示系统更为深刻的特性。当把这种理论推广到更为复杂的系统时,经典控制理论就显得无能为力了,这是因为它的以下几个特点所决定。 1.经典控制理论只限于研究线性定常系统,即使对最简单的非线性系统也是无法处理的;这就从本质上忽略了系统结构的内在特性,也不能处理输入和输出皆大于1的系统。实际上,大多数工程对象都是多输入-多输出系统,尽管人们做了很多尝试,但是,用经典控制理论设计这类系统都没有得到满意的结果;2.经典控制理论采用试探法设计系统。即根据经验选用合适的、简单的、工程上易于实现的控制器,然后对系统进行分析,直至找到满意的结果为止。虽然这种设计方法具有实用等很多完整,从而促使现代控制理论的发展:对经典理论的精确化、数学化及理论化。优点,但是,在推理上却是不能令人满意的,效果也

系统辨识与自适应控制论文

XXXXXXXXXX 系统辨识与自适应控制课程论文 题目:自适应控制综述与应用 课程名称:系统辨识与自适应控制 院系:自动化学院 专业:自动化 班级:自动化102 姓名: XXXXXX 学号: XXXXXXXXX 课程论文成绩: 任课教师: XXXXX 2013年 11 月 15 日

自适应控制综述与应用 一.前言 对于系统辨识与自适应控制这门课,前部分主要讲了系统辨识的经典方法(阶跃响应法、频率响应法、相关分析法)与现代方法(最小二乘法、随机逼近法、极大似然法、预报误差法)。对于系统辨识,简单的说就是数学建模,建立黑箱系统的输入输出关系;而其主要分为结构辨识(n)与参数辨识(a、b)这两个任务。 由于在课上刘老师对系统辨识部分讲的比较详细,在此不再赘述,下面讨论自适应控制部分的相关内容。 对于自适应控制的概念,我觉得具备以下特点的控制系统,可以称为自适应控制系统: 1、在线进行系统结构和参数辨识或系统性能指标的度量,以便得到系统当前状态的改变情况。 2、按一定的规律确定当前的控制策略。 3、在线修改控制器的参数或可调系统的输入信号。 二.自适应控制综述 1.常规控制系统与自适应控制系统比较 (1)控制器结构不同 在传统的控制理论与控制工程中,常规控制系统的结构主要由控制器、控制对象以及反馈控制回路组成。 而自适应控制系统主要由控制器、控制对象、自适应器及反馈控制回路和自适应控制回路组成。 (2)适用的对象与条件不同 传统的控制理论与控制工程中,当对象是线性定常、并且完全已知的时候,才能进行分析和控制器设计。无论采用频域方法,还是状态空间方法,对象一定是已知的。这类方法称为基于完全模型的方法。在模型能够精确地描述实际对象时,基于完全模型的控制方法可以进行各种分析、综合,并得到可靠、精确和满意的控制效果。 然而,有一些实际被控系统的数学模型是很难事先通过机理建模或离线系统辨识来确知的,或者它们的数学模型的某些参数或结构是处于变化之中的.对于这类事先难以确定数学模型的系统,通过事先整定好控制器参数的常规控制往往难以对付。 面对上述系统特性未知或经常处于变化之中而无法完全事先确定的情况,如何设计一个满意的控制系统,使得能主动适应这些特性未知或变化的情况,这就 是自适应控制所要研究解决的问题.自适应控制的基本思想是:在控制系统的运行过程中,系统本身不断地测量被控系统的状态、性能和参数,从而“认识”或“掌握”系统当前的运行指标并与期望的指标相比较,进而作出决策,来改变控制器的结构、参数或根据自适应规律来改变控制作用,以保证系统运行在某种意义下的最优或次优状态。按这种思想建立起来的控制系统就称为自适应控制系统。

闭环系统辨识 报告

闭环系统辨识 气动参数辨识在导弹研发中的作用 气动力参数辨识是飞行器系统辨识中发展最为成熟的一个领域。对于导弹而言,采用系统辨识技术从飞行试验数据获取导弹空气动力特性,已经成为导弹研制和评估程序的重要组成部分。导弹气动参数辨识的作用主要体现在以下几个方面: (1)验证气动力数值计算和风洞试验结果。如前所述,数值计算和风洞试验各有其优点,也各有其局限性,必须通过飞行试验进行验证。如果飞行试验气动参数辨识结果与数值计算和风洞试验结果一致,则说明数值计算和风洞试验结果是正确的;如果不一致,就要找出产生不一致的原因,通过相关性分析,将地面试验结果换算到真实飞行状态下。 (2)为导弹系统仿真提供准确的气动参数。在导弹打靶仿真中,控制系统的执行元件、旋转台、控制系统、目标源等都可以采用实物,但导弹所受外作用力,特别是空气动力是飞行状态参数的函数,无法用实物实现,应代之以数学模型。该数学模型是否正确决定了系统仿真的置信度,因此,采用系统辨识技术,辨识出导弹的外作用力数学模型,特别是气动力数学模型,是导弹系统仿真技术的关键环节之一。 (3)为导弹飞行控制系统设计提供准确的气动参数。控制律设计取决于导弹的气动特性。如果控制律设计所依赖的气动数据误差过大,可能会导致控制失效;如果气动数据误差带很大,为了满足控制系统鲁棒性要求,或者控制精度降低,或者对指令的响应时间加长。利用飞行试验气动参数辨识结果,经过相关性分析给出的导弹气动特性,其可信度可望显著提高,用于飞行控制律设计,可以大大提高控制系统的性能。 (4)自适应控制。自适应控制系统能根据系统的状态和环境参数变化,自动调节控制系统的相应系数,以达到最佳控制状态。系统实时辨识是自适应控制系统的重要组成部分。对于导弹,机动性与导弹的静稳定裕度和动压关系很大,实

自适应控制大作业

自适应控制结课作业 班级: 组员: 2016年1月

目录 1 遗忘因子递推最小二乘法 (1) 1.1最小二乘理论 (1) 1.2带遗忘因子的递推最小二乘法 (1) 1.2.1白噪声与白噪声序列 (1) 1.2.2遗忘因子递推最小二乘法 (2) 2.2仿真实例 (3) 2 广义最小方差自校正控制 (5) 2.1广义最小方差自校正控制 (5) 2.2仿真实例 (6) 3 参考模型自适应控制 (9) 3.1参考模型自适应控制 (9) 3.2仿真实例 (12) 3.2.1数值积分 (12) 3.2.2仿真结果 (12) 参考文献 (16)

1 遗忘因子递推最小二乘法 1.1最小二乘理论 最小二乘最早的想法是高斯在1795年预测行星和彗星运动轨道时提出来的,“未知量的最大可能的值是这样一个数值,它使各次实际观测和计算值之间的差值的平方乘以度量其精确度的数值以后的和为最小”。这一估计方法原理简单,不需要随机变量的任何统计特性,目前已经成为动态系统辨识的主要手段。最小二乘辨识方法使其能得到一个在最小方差意义上与实验数据最好拟合的数学模型。由最小二乘法获得的估计在一定条件下有最佳的统计特性,即统计结果是无偏的、一致的和有效的。 1.2带遗忘因子的递推最小二乘法 1.2.1白噪声与白噪声序列 系统辨识中所用到的数据通常含有噪声。从工程实际出发,这种噪声往往可以视为具有理想谱密度的平稳随机过程。白噪声是一种最简单的随机过程,是由一系列不相关的随机变量组成的理想化随机过程。白噪声的数学描述如下:如果随机过程()t ξ均值为0,自相关函数为2()σδτ,即 2()()R ξτσδτ= 式中,()δτ为单位脉冲函数(亦称为Dirac 函数),即 ,0 ()0,0τδττ∞=?=? ≠?,且-()1d δττ∞ ∞ =? 则称该随机过程为白噪声,其离散形式是白噪声序列。 如果随机序列{}()V k 均值为零,且两两互不相关,即对应的相关函数为: 2,0 ()[()()]0,0v n R n E v k v k n n σ?==+=?=? 则这种随机序列称为白噪声序列。其谱密度函数为常数2(2)σπ。白噪声序列的功率在π-到π的全频段内均匀分布。 建立系统的数学模型时,如果模型结构正确,则模型参数辨识的精度将直接依赖于输入信号,因此合理选用辨识输入信号是保证能否获得理想的辨识结果的

系统辨识报告

系统辨识实验报告

实验一 最小二乘法 1 最小二乘算法 1.1 基本原理 系统模型 )()()()()(11k n k u z B k z z A +=-- a a n n z a z a z a z A ----++++= 221111)( b b n n z b z b z b z B ----+++= 22111)( 最小二乘格式 )()()(k n k h k z T +=θ [][] ?????=------=T n n T b a b a b b a a n k u k u n k z k z k h 11)()1()()1()(θ 对于L k ,,2,1 =,构成线性方程组 L L L n H z +=θ 式中, []T L L z z z z )()2()1( = []T L L n n n n )()2()1( = ? ????? ???? ??--------------= ??????????????=)()1()()1()2()1()2()1()1() 0() 1()0()()2()1(b a b a b a T T T L n L u L u n L z L z n u u n z z n u u n z z L h h h H 参数估计值为 ()L T L L T L LS z H H H 1 ?-=θ 1.2 Matlab 编程 % 基本最小二乘法LS clear;clc A=ones(5,1);B=ones(4,1);%A 为首1多项式,B 中体现时滞(d=1) na=length(A)-1;nb=length(B); load dryer2

系统辨识

作业1 如图1.1所示一阶系统,系统传递函数为G(s)=1/(0.1s+1),如果采用M序列作为输入信号进行系统辨识,采用5级移位寄存器产生M序列作为输入信号,取M序列的时钟脉冲△=15ms,a=2辨识该系统的脉冲响应。并说明取5级移位寄存器合理与否。 图1.1 一阶RC系统 答: 1.解题步骤 1.初始化参数,设置模型参数,设置产生M序列的各个关键参数; 2.利用产生伪随机二进制序列信号的函数getPRBS产生M序列,并作为 系统输入; 3.通过系统模型,产生系统输出,并将输入输出画在同一图中; 4.计算系统输入输出相关函数R xy; 5.计算系统脉冲估计值ghat和系统真实脉冲输出g 2.程序清单 主程序 clc; close all; clear all; %% Initialization R = 100e3; % system initialization resistance=100k ohm C = 1e-6; % capacitance=1uf tc = R*C; % Time Constant % generate M-sequence n=5; a=2; % Level of the PRBS

del = 15e-3; % clock pulse period N=2^n-1; % Period of M sequence total=2*N; % Generate m-sequence using the 'getPRBS' function Out = getPRBS(n,a,del,total); % Generate response y(t) of the system s = tf('s'); G = 1/(tc*s+1); tf = total*del; tim = 0:del:tf-del; y = lsim(G, Out, tim); %plot input and output of the system figure stairs(tim,Out); axis([0 1.0 -2.5 2.5]); hold on plot(tim,y,'r'); hold off % Compute Rxy(i*del) sum = 0.0; Rxy = []; iDel_vec=[]; for i=1:N tau=i-1; iDel_vec=[iDel_vec;tau*del]; for j=1:N sum=sum+sign(Out(j))*y(j+tau); end Rxy_i = (a/N)*sum; sum=0.0; Rxy = [Rxy; tau Rxy_i]; end % Compute ghat & g ind = length(Rxy); C = -Rxy(ind, 2); S = (N+1)*a^2*del/N; Rxy_iDel = Rxy(:,2); ghat=(Rxy_iDel+ C )/S; ghat(1)=2*ghat(1); g = 10*exp(-10.*iDel_vec);

2003版系统辨识最小二乘法大作业

西北工业大学系统辩识大作业 题目:最小二乘法系统辨识

一、 问题重述: 用递推最小二乘法、加权最小二乘法、遗忘因子法、增广最小二乘法、广义最小二乘法、辅助变量法辨识如下模型的参数 离散化有 z^4 - 3.935 z^3 + 5.806 z^2 - 3.807 z + 0.9362 ---------------------------------------------- = z^4 - 3.808 z^3 + 5.434 z^2 - 3.445 z + 0.8187 噪声的成形滤波器 离散化有 4.004e-010 z^3 + 4.232e-009 z^2 + 4.066e-009 z + 3.551e-010 ----------------------------------------------------------------------------- = z^4 - 3.808 z^3 + 5.434 z^2 - 3.445 z + 0.8187 采样时间0.01s 要求:1.用Matlab 写出程序代码; 2.画出实际模型和辨识得到模型的误差曲线; 3.画出递推算法迭代时各辨识参数的变化曲线; 最小二乘法: 在系统辨识领域中 ,最小二乘法是一种得到广泛应用的估计方法 ,可用于动态 ,静态 , 线性 ,非线性系统。在使用最小二乘法进行参数估计时 ,为了实现实时控制 ,必须优化成参数递推算法 ,即最小二乘递推算法。这种辨识方法主要用于在线辨识。MATLAB 是一套高性能数字计算和可视化软件 ,它集成概念设计 ,算法开发 ,建模仿真 ,实时实现于一体 ,构成了一个使用方便、界面友好的用户环境 ,其强大的扩展功能为各领域的应用提供了基础。对 4324326.51411.5320120232320 Y s s s s G U s s s s ++++== ++++432 120120232320 E N W s s s s == ++++

最优控制实验报告

实验报告 课程名称:现代控制工程与理论实验课题:最优控制 学号:12014001070 姓名:陈龙 授课老师:施心陵

最优控制 一、最优控制理论中心问题: 给定一个控制系统(已建立的被控对象的数学模型),选择一个容许的控制律,使被控对象按预定要求运行,并使给定的某一性能指标达到极小值(或极大值) 二、最优控制动态规划法 对离散型控制系统更为有效,而且得出的是综合控制函数。这种方法来源于多决策过程,并由贝尔曼首先提出,故称贝尔曼动态规划。 最优性原理:在一个多级决策问题中的最优决策具有这样的性质,不管初始级、初始状态和初始决策是什么,当把其中任何一级和状态做为初始级和初始状态时,余下的决策对此仍是最优决策 三、线性二次型性能指标的最优控制 用最大值原理求最优控制,求出的最优控制通常是时间的函数,这样的控制为开环控制当用开环控制时,在控制过程中不允许有任何干扰,这样才能使系统以最优状态运行。在实际问题中,干扰不可能没有,因此工程上总希望应用闭环控制,即控制函数表示成时间和状态的函数。 求解这样的问题一般来说是很困难的。但对一类线性的且指标是

二次型的动态系统,却得了完全的解决。不但理论比较完善,数学处理简单,而且在工际中又容易实现,因而在工程中有着广泛的应用。 一.实验目的 1.熟悉Matlab的仿真及运行环境; 2.掌握系统最优控制的设计方法; 3.验证最优控制的效果。 二.实验原理 对于一个给定的系统,实现系统的稳定有很多途径,所以我们需要一个评价的指标,使系统在该指标下达到最优。如果给定指标为线性二次型,那么我们就可以利用MATLAB快速的计算卡尔曼增益。 三.实验器材 PC机一台,Matlab仿真平台。 四.实验步骤 例题1 (P269)考虑液压激振系统简化后的传递函数方框图如下,其中K a为系统前馈增益,K f为系统反馈增益,w h为阻尼固有频率。(如图5-5所示) 将系统传递函数变为状态方程的形式如下: ,

系统辨识

系 统 辨 识 作 业 系统辨识作业: ?已知某系统为单输入/单输出系统,其测量噪声为有色噪声,分布未知。 现给出一个实验样本(如下表所示),求该系统模型。 说明: 可采用GLS ,ELS ,IV 等,要定阶,要比较仅用RLS 的计算结果 一、问题分析 在估计模型参数时需要已知模型的阶数,但是由于本系统模型阶数也是未知的,所以本系统需要先由输入/输出数据通过辩识得出系统的阶数。然后根据辨识的系统阶数再分析求解系统模型。 二、模型阶数的辨识 按照品质指标“残差平方总和”定阶,如高阶系统模型相应的系数为零,则可退化成相应的低阶系统即低阶模型可视为高阶模型的特例。理论上高阶模型的精度不低于低阶模型,但是考虑到计算机的舍入误差的影响,过高的阶数亦能引起模型精度的下降。一般说低阶模型描述粗糙,高阶模型精度高,但是代价亦大。根据逼近的观点,定阶往往是考虑多种因素的折衷。定阶一般是按照假设——检验的步骤进行的,检验过程中往往带有主观成分。 一般说来低阶模型描述粗糙,高阶模型精度高。残差平方总和J(n)是模型阶数的函数 在不同的模型阶数的假设下,参数估计得到的J(n)值亦不同。定阶的最简单办法是直接用J(n)。设模型阶数的“真值”为n 0 ,当n < n 0 时随着n 的增加,J(n)值将明显的下降;而当n ≥ n 0 时随着n 的增加,J(n)值变化将不显著。因此,由J(n)曲线随着n 的增加最后一次陡峭下降的n 值定做n 的估计值。用数理统计的检验方法,判断n 的增加使得J(n)值改善是否明显。 讨论如下 (1).当n=1时程序如下: clear u=zeros(100,1);%构造输入矩阵 z=zeros(100,1);%构造输出矩阵 u=[-0.93249 0.34935 0.76165 -0.9964 -0.38894 -0.12288 0.021565 -0.49555 -0.61624 -1.912 0.22207 -0.31231 -0.17866 -1.8356 -0.26472 1.7642 -1.0418 1.1146 -2.0856 0.8152 1.5094 -0.5822 0.61097 0.35521 2.5907 1.5843 -0.9603 -0.27341 0.39947 0.17493 -1.7451 0.8112 1.2645 1.5682 0.63959 -0.47757 0.99697 0.058774 -0.16174 -1.2928 -0.04722 0.73182 -0.19644 0.091783 -1.1908 -0.90716 0.85388 0.33836 0.74074 0.54181 0.15676 -0.50569 -0.17521 1.3255 -2.488 0.50261 -1.1533 0.36407 0.65283 -0.05983 ∑=-=N k T K k y n J 12 ) )(()(θ?

系统辨识大作业加学习心得

论文 系统辨识 姿态角控制 1.系统辨识概述 辨识、状态估计和控制理论是现代控制理论三个相互渗透的领域。辨识和状态估计离不开控制理论的支持,控制理论的应用又几乎不能没有辨识和状态估计技术。随着控制过程复杂性的提高,控制理论的应用日益广泛,但其实际应用不能脱离被控对象的数学模型。然而在大多数情况下,被控对象的数学模型是不知道的,或者在正常运行期间模型的参数可能发生变化,因此利用控制理论去解决实际问题时,首先需要建立被控对象的数学模型。系统辨识正是适应这一需要而形成的,他是现代控制理论中一个很活跃的分支。社会科学和自然科学领域已经投入相当多的人力去观察、研究有关的系统辨识问题。 系统辨识是建模的一种方法,不同的学科领域,对应着不同的数学模型。从某种意义上来说,不同学科的发展过程就是建立他的数学模型的过程。辨识问题可以归结为用一个模型来表示可观系统(或将要改造的系统)本质特征的一种演算,并用这个模型吧对客观系统的理解表示成有用的形式。当然可以刻有另外的描述,辨识有三个要素:

数据,模型类和准则。辨识就是按照一个准则在一组模型类中选择一个与数据拟合得最好的模型。总而言之,辨识的实质就是从一组模型类中选择一个模型,按照某种准则,使之能最好地拟合所关心的实际过程的静态或动态特性。 通过辨识建立数学模型的目的是估计表征系统行为的重要参数,建立一个能模仿真实系统行为的模型,用当前可测量的系统的输入和输出预测系统输出的未来演变,以及设计控制器。对系统进行分析的主要问题是根据输入时间函数和系统的特性来确定输出信号。对系统进行控制的主要问题是根据系统的特性设计控制输入,使输出满足预先规定的要求。而系统辨识所研究的问题恰好是这些问题的逆问题。通常,预先给定一个模型类{}M(即给定一类已知结构的模型),一类输入信号u和等价准则(,)JLyyM(一般情况下,J是误差函数,是过程输出y和模型输出yM的一个泛函);然后选择是误差函数J达到最小的模型,作为辨识所要求的结果。系统辨识包括两个方面:结构辨识和参数估计。在实际的辨识过程中,随着使用的方法不同,结构辨识和参数估计这两个方面并不是截然分开的,而是可以交织在一起进行的 一、控制对象 本文采用了控制不同电机转速组合的方法,对四轴旋翼蝶形飞行器进行姿态控制,使四旋翼蝶形飞行器在不同姿态下飞行时具有较好的性能。为了实现四轴旋翼蝶形飞行器的飞行控制,对飞行的控制系统进行了初步的设计,并给出了设计流程。同时利用matlab对四轴旋翼

matlab实验报告

专业仿真课程设计题目: 学院: 专业班级: 学号: 学生姓名: 指导教师: 设计时间:

专业仿真课程设计题目 主要研究内容: 从所拍摄的多个目标物中检测三角形物,给出三角形物几何中心、三个边长以及边长的方向、面积。 设计要求: (1)提交能够实现题目要求、并通过演示验收的可执行文件。 (2)提交课程设计报告(包括程序清单)。 (3)通过答辩,答辩成绩满分20分,其中个人设计部分10分,非个人设计部分10分。 (4)软件设计要求:有一个人机交互界面,模块化设计,在模块之间通过BMP文件或者文本文件传送数据,可以查看中间结果。 (5)5个人一组,组长协调分工,每个组员一定要有具体任务,以便考核。预期达到的目标: 1、能够通过相关文献查阅、文献综述和总结,给出问题求解的多种可行方案。 2、能够综合运用测控技术与仪器专业理论和技术手段,设计实验方案、分析实验结果,得出有效的结论。 3、能够借助MATLAB仿真软件,进一步掌握高等数学、复变函数与积分变换等相关数学和自然科学知识以及测控技术与仪器专业的基本理论知识,能够结合本专业“自动控制原理”、“数字信号处理”、“误差理论”等相关课程,采用MATLAB软件对复杂工程问题建立模型并进行预测与模拟; 4、能够与团队中其他学科成员合作开展工作,能够与其他队员很好地沟通和交流意见,能够通过口头或书面方式表达自己的设计思路,具有一定的表达能力和人际交往能力。

目录 第一章课程设计相关知识综述 1.1 MATLAB相关知识叙述 1.1.1 MATLAB基本知识介绍 1.1.2 MATLAB的优势特点 1.1.3 MATLAB的发展历程 1.2 MATLAB工具箱与函数 1.2.1 MATLAB图像处理工具箱 1.2.2 课程设计所用图像处理函数介绍第二章课程设计内容和要求 2.1 课程设计主要研究内容 2.2 课程设计要求 2.3 课程设计预期目标 第三章设计过程 3.1 设计方案 3.2 设计步骤及流程图 3.3 程序清单及相关注释 3.4 实验结果分析 3.5 结论 第四章团队情况 第五章总结 第六章参考文献

相关文档
最新文档