系统辨识的matlab实现方法(手把手).pdfVIP

最近在做一个项目的方案设计，应各位老总的要求，只有系统框 图和器件选型可不行，为了凸显方案设计的高大上，必须上理论分析， 炫一下“技术富”，至于具体有多大实际指导意义，那就不得而知了！ 本人也是网上一顿百度，再加几日探索，现在对用 matlab 实现系统辨 识有了一些初步的浅薄的经验，在此略做一小节。 必须要指出的是，本文研究对象是经典控制论理最简单最常用的 线性时不变的 siso 系统，而且是 2 的哦，也就是具有如下形式的传 递函数： 1 G(s)  2 2 T s  2 Ts  1 本文要做的就是，对于有这样传递函数的一个系统，要辨识得到 其中的未知数 T， ξ!!这可是控制系统设计分析的基础哦，没有系统 模型，啥理论、算法都是白扯，在实际工程中非常重要哦！ 经过总结研究，在得到系统阶跃响应实验数据之后（当然如果是其 他响应，也有办法可以辨识，在此还是只讨论最简单的阶跃响应实验 曲线，谁让你我是菜鸟呢），利用 matlab 至少可以有两种方法实现实 现（目前我只会两种，呵呵）！ 一、函数法 二、GUI 系统辨识工具箱 下面分别作详细介绍 ！ 一、 函数法 看官别着急，先来做一段分析（请看下面两排红*之间部分），这段 分析是网上找来的，看看活跃一下脑细胞吧，如果不研读一下，对于 下面 matlab 程序，恐怕真的就是一头雾水咯！ ******************************************************* ************************ 1 G(s)可以分解为：G(s)  T 2 (s   )(s   ) 1 2 1  2  1      1 T 其中， 1  2  2      1 T  、 都是实数且均大于零。 1 2 则有： 1    T  ，  1 2   2   1 2 1 2 传递函数进一步化为：   G(s)  1 2 (s   )(s   ) 1 2 因此，辨识传递函数就转化为求解 、 。 1 2 当输入为单位阶跃函数时，对上式进行拉普拉斯反变换，得系统 时域下的单位阶跃响应为： 