MATLAB system identification系统辨识app的使用

系统辨识

      • 前言
      • 系统辨识
        • [第一步 选取时域数据到app](#第一步 选取时域数据到app)
        • [第二步 分割数据](#第二步 分割数据)
        • [第三步 设置传递函数的参数](#第三步 设置传递函数的参数)
        • [第四步 Estimate](#第四步 Estimate)
        • 第五步
      • 结束

前言

接上节:simulink-仿真以及PID参数整定

系统模型的辨识工作,在控制领域,一般用于开发控制器的先手工作。一般而言,设计控制器,会依据被控对象的数学模型。依据其数学模型,可以分析其各种特性,所以数学模型就显得很重要。

数学模型,通俗一点讲就是一个数学表达式,f(x, t),这个表达式有个特点就是,给(输入)x 赋值,那么这个表达式就会是一个关于时间 t 的函数,就会在时间轴上,随着时间变化而变化,呈现出的输出值,就会和实际的物理模型的输出近似。

为什么说是近似呢,实际生活中,物理模型的实际输出大多都会相对复杂,并不是有特别明显的规律,有的可能会随时间、环境以及各种情况变化而变化,那么还用一个固定的公式表示其输出,就不准确了。也有一些简单的物理模型,这样的模型就很容易,就直接可以用简单的数学公式表示其输出。

那么可能会有新手和我一样疑惑,既然有的不准确,那就用一些高级的方法去拟合得到实际物理对象的数学模型。有没有方法,也有,但是考虑到情况的复杂度,用统计学的方法去拟合,也不能覆盖全部的情况下的数学模型,所以随后的数学模型也只能做到近似。

那么问题又来了,既然是近似,横竖结果不能完全一致,那何必用那么复杂的方法去得到数学模型呢。所以就是这个想法,就有了高阶模型线性化处理,或者降阶数,能够用最简单的一阶模型去表示物理对象的输出,就不用二阶模型,就秉承这个原则,去简化数学模型。

那么什么情况下用一阶什么时候用二阶呢,有这么一段话,但是很绕口,简单说就是要给物理对象一个输入信号,可以是阶跃信号,看从一个值变化到另一个值,物理对象的实际输出结果的变化曲线,由变化曲线来判断。

这里说的是对一个系统不知道的情况下,采取的办法,这个办法又被通俗的叫作工程试验法,给一个激励信号,看输出动态响应,然后有输入和输出结果,就可以借助工具箱里的系统辨识去拟合系统的数学模型。

在开始说系统辨识的方法前,还要解释一下前边的两个问题:

  1. 既然最后得到的数学模型只是近似模型,那么设计的控制器能有用么?这个问题就是控制器需要解决的问题,最经典的PID控制器,就是包含了这点,它允许模型不准确。所以大多的控制器都具备一定的容错能力,本身控制器的作用就是降低误差,只要这个误差不是太大,就能救回来。那些严格绑定被控对象数学模型的控制器,就更加适用在准确的数学模型上。
  2. 准确的数学模型又从哪里来呢?最常见的运动的质点,物理中学过的 s = v0 x t + 1/2 x a x t^2 , 这个质点可以稍微放大一下,就变成遥控小车。那么遥控小车的运动学模型就是这个公式,输入加速度后,路程就和时间相关,随时间变化。但是实际小车没这么简单,小车还有质量,小车也不是手推的,需要提供电,由电机转动驱动小车,那么就涉及驱动力多大,能提供多大加速度,小车需要变速,那么要平滑变速,就需要加速度控制,那么就要有动力学公式,一个运动学,一个动力学,两个公式就是小车的数学模型,在不考虑什么复杂环境下,风阻,地面摩擦力降低这些情况,数学模型就是精准的,所以科研分析大多是先有推导的公式,然后才会有一些没法计算推导得到的系统进行扫频然后系统辨识 。

系统辨识

上边说了,辨识模型就需要激励信号输入,和响应信号输入,由输入、输出才能辨识。

常用的激励信号,是一个组合信号,阶跃、正弦、斜坡、啁啾信号:

然后会得到一个对应的输出:

这个输出肯定不是一个实际的物理对象的输出,一般情况下,实际输出和输入不在一个维度上,比如上边的小车模型,实际输出是路程,输入是加速度,在做控制器的时候,路程和加速度没法做加减,为什么作加减,闭环控制嘛,就得把输出送回来和输入目标值进行对比嘛,那么就得在同一个维度下进行了。为了方便辨识,我们把输出转换成输入相同的维度,这样就可以直接使用,最后得到的系统输出的值,就是转换后,这个在最后不要忘记了。

图片是为了展示效果,实际上要把数据导入matlab整理好成这样:

第一步 选取时域数据到app
第二步 分割数据

这一步将数据分割成测试集和验证集

然后就会得到这个状态:

第三步 设置传递函数的参数

最保守的就是第一个,或者 process models 比较直观,或者就是最后一个什么都不用设置

第四步 Estimate


第五步

在选择阶数前可以先用state space model 推荐最合适的阶数

然后就来回试一试,找到拟合度最好的一个模型,那个就是辨识得到的模型:

结束

打完收工,其实得到的模型只是一个辅助,应为采集到的数据到拟合,有很多需要注意的地方,最后得到的模型不一定理想。

相关推荐
无限大.几秒前
c语言200例 067
java·c语言·开发语言
余炜yw2 分钟前
【Java序列化器】Java 中常用序列化器的探索与实践
java·开发语言
篝火悟者3 分钟前
问题-python-运行报错-SyntaxError: Non-UTF-8 code starting with ‘\xd5‘ in file 汉字编码问题
开发语言·python
Death2006 分钟前
Qt 中的 QListWidget、QTreeWidget 和 QTableWidget:简化的数据展示控件
c语言·开发语言·c++·qt·c#
六点半8887 分钟前
【C++】速通涉及 “vector” 的经典OJ编程题
开发语言·c++·算法·青少年编程·推荐算法
惜.己7 分钟前
javaScript基础(8个案例+代码+效果图)
开发语言·前端·javascript·vscode·css3·html5
niu_sama11 分钟前
基于muduo库函数实现protobuf协议的通信
开发语言·qt
不写八个13 分钟前
Qt教程(001):Qt概述与安装
开发语言·qt
quaer19 分钟前
Open-Sora全面开源?
开发语言·算法·机器学习·matlab·矩阵
hakesashou23 分钟前
python如何比较字符串
linux·开发语言·python