《基于matlab在连续性系统中的仿真研究 计算机科学和技术专业》由会员分享，可在线阅读，更多相关《基于matlab在连续性系统中的仿真研究 计算机科学和技术专业（23页珍藏版）》请在金锄头文库上搜索。
1、基于matlab在连续性系统中的仿真研究摘要：从数学建模到仿真系统的搭建，再到加进控制环节进行实时控制，最后得出结果的过程中，参考了大量的资料，通过对比整合，设计出了适合自己的一套实验方法：本研究的核心在于构建相应的倒立摆数学模型，首席按通过牛顿-欧拉法完成相关模型的基本搭建，之后通过动态系统空间状态方程来确定其系数矩阵的值，完成上述工作之后，通过MATLAB编程的方式，将构建的数学模型转化成相应的传递函数，从而获得传递函数模型。此后便可在此基础之上，通过Simulink软件通过PID控制算法实现相应的仿真工作，先用连续系统的设计方法设计出模拟控制器，然后在满足一定条件下，对其进行离散化处理，
2、(采用加零阶保持器的Z变换法)形成数字控制器。接着进行PID参数整定，利用试凑法,根据PID控制器各部分组件和性能之间的联系，在某个初始化参数的条件下重复进行试凑工作，最终得到较好的结果。研究过程中，系统的控制非常稳定，性能较好。关键词：Matlab；仿真；PID；连续系统目录1绪论11.1连续性系统简介11.2连续性系统基本处理方法11.2.1数值积分法11.2.2替换法21.2.3离散相似法31.3MATLAB在连续性系统的应用52 MATLAB基础62.1MATLAB语言介绍62.2MATLAB工具箱62.3MATLAB对连续性系统的仿真性能分析73直线一级倒立摆系统设计73.1倒立摆系
3、统概述73.2倒立摆系统工作原理83.3倒立摆系统数学模型建立93.3.1确定系统输入输出量及中间变量：103.3.2受力分析，列写运动方程：103.3.3系统非线性方程的线性化：113.3.4零初始条件下的拉氏变换：113.3.5代入参数求解传递函数：114基于MATLAB的一级倒立摆系统仿真结果及分析124.1校正前系统性能仿真分析124.1.1稳定性分析：124.1.2阶跃响应分析：124.1.3频率特性分析：134.2PID控制设计校正装置144.2.1计算开环增益K：144.2.2计算校正后的增益交界频率c：164.2.3计算参数确定传递函数：164.3 控制实验仿真分析165总结和
4、展望21参考文献2221绪论1.1连续性系统简介所谓连续系统，其实就是体系的状态会在时间迁移的条件下产生光滑且连续的改变的一个系统。其也包含那些因离散采样而引起非连续改变的体系。这类系统通常能够通过微分方程组进行表示。如果方程之中的相关系数皆是常数，那么可将之称作定常系统，如果相关系数会因时间的改变发生改变，那么可将之称作。除此之外，相关的数学模型也可分为多个类别，除了连续模型（即通过微分方程进行表示）之外，还有离散模型（通过差分方程进行表示）以及二者的混合模型。如果某个体系能够通过非线性的常微分方程组进行表示，那么其就属于非线性系统。其中的各个方程不仅含有状态变量及其各阶导数的线性组合，也含
5、有非线性项。其非线性特性主要表现在：（1）线性系统同时满足均匀性与叠加性。线性方程必定有解，并且其解为零输人响应和零状态响应之和。非线性系统则不具有这些特征，即对非线性系统，通常叠加性原理不成立，非线性方程通常没有解析表达形式的解。（2）给定一个线性系统，对于不同的初始状态，它只有一种类型的运动。但对于非线性定常自由系统则不然，初始状态不同，它对应的运动类型可以不一样。非线性系统有三种类型的运动，即稳定运动、不稳定运动和周期运动。1.2连续性系统基本处理方法1.2.1数值积分法如果通过相关的状态方程来表示被模拟的体系，其对应的一阶微分方程以及初始值可以通过下式表示： （1-1）在这之中，Y代表
6、维度为n的状态向量，而F则代表维度为n的向量函数。对于上式，当时其对应的连续解可以写成：（1-2）设 ，令（1-3）则有：也就是说，（1-4）如果准确解为近似值，是准确积分值的近似值，则式（1-4）就是式（1-2）的近似公式。换句话说，连续系统的数值解就转化为求两个相邻时刻之间的积分。综上所述，（1-1）的数值求解策略本质上就是求得该式的实际解在相应的离散时间点之中的近似值，对于处在相邻时刻的两个点，其间距，通常被命名为步长，一般令，保持不变。1.2.2替换法以上使用的通过数值积分对连续体系进行模拟的策略尽管理论上较为完善、求解精度也较为不错，然而过程较为繁琐，需要进行大量的计算，一般仅将之应
7、用于对于效率没有过高要求的模拟工作之中。如果对于模拟工作有着实时性的需要，要求计算速度快，以便能在一个采样周期内完成全部计算任务，这就必须使用其他的高效率的求解方法。通过数值积分策略实现对连续体系的模拟，本质上已对其实施了离散化的操作，只是在进行离散化的各个步骤之中都会使用连续体系的模型，每一步都是离散之后再进行求解。因此，可以先将相关连续模型实施离散化操作，获得等价的离散模型，在这之后的各个计算工作中可以不必再涉及最初的连续模型，而是能够直接使用该离散模型。通过对连续模型简化得到相应的离散模型，能够很好地适应计算机上的计算环境，不但能够应用在连续体系的数字模拟工作中，同时也能够通过计算机实现
8、相应的数字控制器。替换法的核心思路为：对某个传递函数G（s），寻找将s域映射至z域的策略，从而实现属于S域的变量到z平面之间的映射工作，从而能够将和连续体系相关的G（s）转化为离散的G（z）。之后可以通过G（z）进行相应的z反变换，得到对应的通过差分方程进行表示的离散模型，从而实现高效率的计算工作。根据z变换理论，s域到z域的最基本的映射关系是或 如果按这一映射关系直接代入G（s），得到的G（z）是相当复杂的，不便于算法实现，所以往往借助于Z变换的基本映射关系或作一些简化和近似处理。1.2.3离散相似法离散相似法能够把某个连续的体系实施离散化操作，并获得与之对应的差分方程形式的离散模型。得到对
9、应的离散模型的策略包括两类： （1）对传递函数实施相应的离散化操作，从而获得离散化的传递函数这类模型通常叫做“频域离散相似模型”；（2）对状态方程实施相应的离散化操作这类模型通常叫做“时域离散相似模型3”；为了实现连续体系的数字模拟工作，首先可向体系之中添加相应的虚拟采样器以及保持器，具体的架构可以参考下图： 图1-1 连续系统离散化结构图附注：上图之中给出的采样开关以及保持器实际上是不存在的，而是为了将(1-5)式离散化而虚构的。然后通过z变换实现对应体系脉冲传递函数的求解，之后可在此基础之上获得体系所对应的差分方程。 根据上图不难得到相应的脉冲传递函数模型，其可通过下式进行表示: （1-5
10、）在这之中，Gh(s)代表保持器的传递函数，并且可根据保持器的种类获得相应的G(z)，具体可参考下表： 表1-1不同保持器的G(z)保持器的传递函数Gh(s)脉冲传递函数G(z)零阶： 一阶： 三角形： 对于某个连续体系，其状态方程可以表示成下式：（1-6）如果在体系的输入以及输出部分都添加相应的采样控制开关，在输入部分额外添加保持器以确保相关信号能够良好地还原，那么此时可获得下图所给出的结构： 图1-2采样控制系统结构图若对方程(1-6)式两边进行拉普拉斯变换，得： 即： 以左乘上式的两边可得 ：（1-7）考虑到状态转移矩阵：（1-8）故对(1-7)式反变换可得：（1-9）此为(1-5)式的
11、连续解，由此可推导出系统的离散解。 根据上式，n及n+1两个相连的采样瞬间，有： （1-10）（1-11）（1-12）将(1-10)式减去(1-10)式后乘以，得： 将(1-11)式右边积分进行变量代换，即令： （1-13）则得：（1-14）但由图1-2可知：若系统采用零阶保持器时，则两个采样点之间输入量可看做常数，即u(nT+t)=u(nT)，这样(1-14)式可写为： 式中 ：1.3MATLAB在连续性系统的应用瑞典的著名高校Lunds大学曾经开展过时间为3个月的极地研究工作，并在这个过程中大量运用到了MATLAB和相关工具箱，其主要用于考察北极的放射性物质和当地环境之间的作用关系。相关研
12、究者通过MATLAB来完成数据处理以及分析工作，同时通过提供的神经网络工具箱实现了对北冰洋之上漂浮木材的年轮图像识别工作，从而确定相关浮木的来源。除此之外，Forsmark核电站也通过MATLAB来完成对核反应堆的功率优化工作，相关的研究者能够在各个堆芯之中得到海量的数据，并在此基础上求解各燃料棒以及控制设备坐标的最优值，从而获得最高的功率。这个数据处理过程的规模非常大，且复杂度也极高，涉及到超过1700个结点的处理工作，为能精简这个流程， Forsmark通过MATLAB开发得到了相应的GUI交互界面，其能够让毫无经验的人员进行数据的处理和分析工作。不仅如此，Forsmark还进一步通过MA
13、TLAB对一系列失效和扰动的场景进行建模处理。一旦反应堆出现了扰动现象，应当对相关信息进行详实的分析，从而实现引发扰动因素的定位工作，利用MATLAB的种种优势，上述工作所耗费的时间已经从原本的7天降低到15分钟。Calspan先进技术中心曾对相关的飞行器进行仿真测试，并有效提高了其性能。相关的模拟软件工作与一系列并行的浮点DSP中，硬件部分则使用的dSPACE。在地面阶段，首先通过Simulink模块实现相关飞行器及其控制系统的建模以及模拟工作，之后通过Workshop产生对应的C语言代码，并传输到飞行器的DSP之中。进行试飞之时，试飞员能够通过DSP模块对该系统进行考察，飞行过程之中涉及到
14、的一系列参数也能够实时调整，同时也可下载到硬件之中进行测试，从而有效降低了成本，提高了效率，能够在控制系统原型的实际构建之前便将相关试验全部完成。2 MATLAB基础2.1MATLAB语言介绍MATLAB诞生初期，其仅仅是用于为数据分析工具Eispack与Linpack提供接口，其最早通过Fortran实现，之后通过C实现。1975年左右，M公司将其推出面向市场，并且将其功能逐渐丰富，开发了文字图像的处理、数值的运算以及符号的解译等功能，将其编程语言升级为面向对象语言，并以此为用户界面的设计语言，从而将其打造成了一个能够在多个领域、多个学科都很好地发挥功能的软件。其版本也逐渐的更新换代，其有效
15、综合了各种数据处理功能，包括矩阵求解、数据可视化、非线性体系的分析等等，目前在工程领域以及科研工作之中都得到了普遍的应用，其主要包括下面几个部分： 2.2MATLAB工具箱现阶段，MATLAB已经拥有了超过37个工具箱，同时其主要包括两个类别，分别是功能类以及领域类。后者通常依赖于较高的特定领域专业知识，主要有信号处理、自控、逻辑分析、图象处理、通讯、偏微分方程等等很多细化的工具箱，可应用的领域涵盖了化学、数学、工程还有经济等等很多的领域范围。功能型的工具箱主要是增强其图形文字处理能力、模型建立和仿真功能、数值符号运算解议以及和软硬件系统的通信等一系列作用。使用者可以直接对相关工具箱的源代码进行自行定制，以满足特定的需求。3直线一级倒立系统设计3.1倒立摆系统概述倒立摆系统是用于进行控制理论分析工作的一个非常经典的设备，其有着非常精简的构造，十分低廉的成本，同时相关参数以及构造可以灵活地更改。该系统处于非稳态，涉及变量多、阶次高、耦合高。通过该系统，可以有效体现出控制流程里面的可能问题，包括稳健性、可定位性等等，该系统能够用于对相关的控制理论进行检验。在各类控制理论的基础
