模糊pid控制算法matlab_matlab模糊控制器

模糊pid控制算法matlab_matlab模糊控制器上一篇写了模糊自整定PID的理论,这篇来做MATLAB仿真。目录补充内容:如何计算临界稳定下的开环增益Ku和震荡周期TuMATLAB进行模糊PID仿真1、准备工作2、模糊控制器的设计补充内容:如何计算临界稳定下的开环增益Ku和震荡周期Tu学过控制工程或者相关理论的同学应该比较了解,判断系统稳定性的条件一般用到劳斯表(劳斯判据)。而PID控制和模糊PI…

大家好,又见面了,我是你们的朋友全栈君。

上一篇写了模糊自整定PID的理论,这篇来做MATLAB仿真。

目录

补充内容:如何计算临界稳定下的开环增益 Ku 和震荡周期 Tu

MATLAB进行模糊PID仿真

1、准备工作

2、模糊控制器的设计


前置说明:由于本人长期在外地出差,还没有时间来做本文中模型的Debug, 也未对每一步深入研究。请各位先参考评论区自行检查,我后续会更新本文。十分抱歉。

补充内容:如何计算临界稳定下的开环增益 Ku 和震荡周期 Tu

    学过控制工程或者相关理论的同学应该比较了解,判断系统稳定性的条件一般用到劳斯表(劳斯判据)。而PID控制和模糊PID控制极大地依赖系统传递函数的建立,因此如果对于系统复杂,难以建立模型的,还是需要考虑一下。

(1)通过自己对所需要控制的系统进行建模,拉普拉斯变换得到传递函数Gc(s),然后求得其特征方程(形式如下):

模糊pid控制算法matlab_matlab模糊控制器

 其闭环系统特征方程可以表示为:

模糊pid控制算法matlab_matlab模糊控制器

特征方程的计算就是要求闭环系统传递函数的分母为0,如果是开环系统,就将传递函数的分子和分母加起来作为特征方程。

(2)将系统特征方程中的s用(-jw)代替,然后令实部和虚部分别等于0,解出特征根。特征根 w 就是震荡频率,K 为开环增益。

则,周期 Td = 2PI/w 。这样就可以用来求取PID参数的论域了。

 以上下划线理论有问题,验证后发现并不是这样计算。

修改内容:

修改PID三个参数的模糊论域为

P = [0 0.1];

I = [0.0.05];

D = [0.0.1];

MATLAB进行模糊PID仿真

1、准备工作

    首先需要选取传递函数,设系统传递函数方程如下:

模糊pid控制算法matlab_matlab模糊控制器

纯属为了计算方便。解得 K=8,w^2=3。

先看使用 Simulink 自带的仿真结果,其PID参数整定的情况:

模糊pid控制算法matlab_matlab模糊控制器

模糊pid控制算法matlab_matlab模糊控制器

接下来设计模糊自整定PID。

2、模糊控制器的设计

步骤:

(1)打开 MATLAB ,选择一个工作目录,这里一般选择一个可写的磁盘即可。为什么要可写的呢,因为C盘有些目录下你的用户权限是不可写的,所以就算你设计的没有问题也没法运行。

模糊pid控制算法matlab_matlab模糊控制器

(2)命令行输入fuzzy,打开模糊规则编辑器。

模糊pid控制算法matlab_matlab模糊控制器

(3)Edit -> Add variables -> Input;Edit -> Add variables -> Output;Edit -> Add variables -> Output。单击选中,分别对 input1 、input2、output1、output2、output3 进行隶属度编辑。对于输入进行命名,分别为误差 E 和误差变化率 EC。输出分别命名为P、I、D,其他不变。

模糊pid控制算法matlab_matlab模糊控制器

(4)双击 Input1 进入隶属度编辑器,Edit -> Add MFS 。添加隶属度函数,这里采用三角形隶属度函数,添加4个,实际应该对每一个都添加4个,并进行编辑。

模糊pid控制算法matlab_matlab模糊控制器

(5)其中 E 和 EC 的论域(Range)取 [-6 6],然后对每个隶属度曲线进行编辑。这里有两个小问题(或者说是小技巧)。

  a. 论域的划分,抛开模糊PID理论来说,其实论域的划分很大程度上取决于个人的方法,灵活性很大。你可以认为身高180的人很高、也可以认为较高,甚至可以认为较矮,只是隶属度不同。这些都取决于你的系统,但是常规是均匀划分的,例如这篇博客。然后我自己想了一种方法(适用于7隶属度):半区间宽度=(论域总长度/12)*2,也就是说一个大三角由两个半区间组成。

  b. 还有就是隶属度函数命名的时候,最好是有顺序的,单击隶属度函数曲线,设计顺序:mf1->NL;mf2->NM;mf3->NS;mf4->ZE;mf5->PS;mf6->PM;mf7->PL。这样设置的好处是后面添加模糊规则的时候很方便。 

模糊pid控制算法matlab_matlab模糊控制器

(6)设计PID三个参数的输出隶属度

   按照上一篇的介绍,三个参数的隶属度选择范围如下,近似的话建议扩大论域区间:

Kp = [0.32Ku 0.6Ku] = [2.56 4.80] \approx [2,5]

Ki = [0.32Tu 0.47Tu] = [1.1608 1.7050] \approx [1 2]     // 这里用[1 2]可能存在结果1的情况,结果2用的[1 3]

Kd = [2Tu 5Tu] = [7.2552 18.1380] \approx [7 19]

按照之前的方法进行划分,获得隶属度函数分布区间:

                                                                               Kp 的论域划分

隶属度 NL NM NS ZE PS PM PL
区间 [1.5 2 2.5] [2 2.5 3] [2.5 3 3.5] [3 3.5 4] [3.5 4 4.5] [4 4.5 5] [4.5 5 5.5]

                                                                              Ki 的论域划分

隶属度 NL NM NS ZE PS PM PL
区间 [0.833 1 1.167] [1 1.167 1.333] [1.167 1.333 1.5] [1.333 1.5 1.667] [1.5 1.667 1.833] [1.667 1.833 2] [1.833 2 2.167]

                                                                             Kd 的论域划分

隶属度 NL NM NS ZE PS PM PL
区间 [5 7 9] [7 9 11] [9 11 13] [11 13 15] [13 15 17] [15 17 19] [17 19 21]

直接复制上去吧!!!!

模糊pid控制算法matlab_matlab模糊控制器

其实这里还有很多不严谨的地方,比如隶属度论域的计算方法。

(7)设计模糊规则,这里的设计规则选取的是一篇文献中的规则表,我主要是考虑方便,实际上,这种模糊规则表是决定模糊PID效果的关键之一(还有就是隶属度函数的确定)。如果这里设计的不好,最后出来的PID阶跃响应可能还不如手工调整PID,虽然这种方式对于在线调整是不太可能实现的。调整规则如下,参考文献 [2]。

双击下图圈中方框:

模糊pid控制算法matlab_matlab模糊控制器

调整规则: 

模糊pid控制算法matlab_matlab模糊控制器

模糊pid控制算法matlab_matlab模糊控制器

模糊pid控制算法matlab_matlab模糊控制器

第一个隶属度图中的圈内换成PL,第二个图换成ZE(其实这里不确定,我换的ZE)。

根据隶属度规则表挨个添加进 Rule Editor 即可

模糊pid控制算法matlab_matlab模糊控制器

(8)导入工作空间

关闭规则编辑器。单击 file-> export -> To file… 命名为Fuzzy_PID.fis -> 保存。

单击 File -> Export -> To workspace -> OK。

模糊pid控制算法matlab_matlab模糊控制器

我们会发现在MATLLAB的工作区(Workspace)有来了一个1*1 struct的值。

(9)Simulink中进行仿真。

    在命令行里输入simulink,或者在MATLAB主页点击Simulink,打开Simulink工具箱。新建一个空白Blank。

模糊pid控制算法matlab_matlab模糊控制器

在库里面找到Fuzzy Logic Controlller模块,拖到Blank里面。相关的还需要step(阶跃信号),sum(输入输出反馈),PID(一个完整的控制算法块,也可以自己写),mux(用于整合图形),scope(显示结果)。

(10)然后一步一步设置:

双击sum,将sum的配置改为下图:

模糊pid控制算法matlab_matlab模糊控制器

双击PID Controller,使source为external,也就是PID接收外部参数。

模糊pid控制算法matlab_matlab模糊控制器

PID改完之后是这样的:

模糊pid控制算法matlab_matlab模糊控制器

双击传递函数(Transfer Fun),修改如下:

模糊pid控制算法matlab_matlab模糊控制器

如果显示 num(s)/den(s) 的形式,可以把它放大,和你放大窗口是一样的。

设置Fuzzy controller,同样双击。

模糊pid控制算法matlab_matlab模糊控制器

这个名字要和你工作空间的一样才行,有些同学虽然改了,但是还是不成功,可能就是你没能导入到工作空间里面。自寻解决办法(readfis()函数)。

模糊pid控制算法matlab_matlab模糊控制器

 以下为错误

模糊pid控制算法matlab_matlab模糊控制器

结果1:

模糊pid控制算法matlab_matlab模糊控制器

结果2:Ki 的区间从[1,2] 修改为[1,3].

模糊pid控制算法matlab_matlab模糊控制器

模糊pid控制算法matlab_matlab模糊控制器

目前结果是这样,我再调一调看看。

不妥之处恳请广大朋友指出。

增加:评论区有位指出积分时间常数和PID系数Ki的问题,这里贴个PID的方程:

典型PID控制器的传递函数为:

G_{c}(s)=K_{p}(1+\frac{1}{T_{i}s}+T_{d}s)

而离散PID控制方程可以表示为:

U(n)=K_{p}[e(n)+\frac{T}{T_{i}}\sum e(i)+\frac{T_{d}}{T}(e(n)-e(n-1)))]

源程序!源程序!GitHub – NidhoghostX/Fuzzy-PID-on-Matlab

CSDN 要收费下载,资源转到Github了。

——————————————————————————————————————————————

更新:

二维模糊PID的matlab仿真(官网教程):Fuzzy PID Control with Type-2 FIS- MATLAB & Simulink

现在在忙,没有时间研究,今后可能会做一下。

参考文献

[1] https://sci-hub.tw/10.1115/1.2899060 (最早经典的PID参数整定文章,大部分文献都会用到)

[2] S. krishna, S. Vasu,Fuzzy PID based adaptive control on industrial robot system,Materials Today: Proceedings,Volume 5, Issue 5, Part 2,2018,Pages 13055-13060.
 

版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/135038.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • Python 根据AIC准则定义向前逐步回归进行变量筛选(二)

    Python 根据AIC准则定义向前逐步回归进行变量筛选(二)Python根据AIC准则定义向前逐步回归进行变量筛选(二)AIC简介AIC即赤池值,是衡量模型拟合优良性和模型复杂性的一种标准,在建立多元线性回归模型时,变量过多,且有不显著的变量时,可以使用AIC准则结合逐步回归进行变量筛选。AICD数学表达式如下:AIC=2p+n(log(SSE/n))AIC=2p+n(log(SSE/n))AIC=2p+n(log(SSE/n))其中,ppp…

    2022年5月24日
    181
  • java垃圾回收机制原理_java垃圾回收的缺点

    java垃圾回收机制原理_java垃圾回收的缺点  最近做一个ETL的项目模块,经常由于查询数据量比较大用消息中间件MQ时引起了内存溢出的报错。做完后没事研究了一下JVM和垃圾回收的相关知识点。一:垃圾回收机制的意义java 语言中一个显著的特点就是引入了java回收机制,是c++程序员最头疼的内存管理的问题迎刃而解,它使得java程序员在编写程序的时候不在考虑内存管理。由于有个垃圾回收机制,java中的额对象不在有“作用域”的概念,只有…

    2022年10月13日
    0
  • 清理电脑垃圾的命令代码

    清理电脑垃圾的命令代码@echooffecho正在清除系统垃圾文件,请稍等……del/f/s/q%systemdrive%\*.tmpdel/f/s/q%systemdrive%\*._mpdel/f/s/q%systemdrive%\*.logdel/f/s/q%systemdrive%\*.giddel/f/s/q%system

    2022年6月18日
    327
  • [Python]多线程详解

    [Python]多线程详解最近在研究pythonthread的一些用法,有看到一篇关于thread的介绍文章,写的很好很详细,看完后转载收藏一下,供温故知新。谢谢作者大大转载于:python多线程详解(超详细)impo

    2022年7月5日
    21
  • Raid0、Raid1、Raid0+1、Raid5

    Raid0、Raid1、Raid0+1、Raid5Raid0:最少需要两块盘,没用冗余数据,不做备份,任何一块磁盘损坏都无法运行。n块磁盘(同类型)的阵列理论上读写速度是单块磁盘的n倍(实际达不到),风险性也是单一n倍(实际更高),是磁盘阵列中存储性能最好的。适用于安全性不高,要求比较高性能的图形工作站或者个人站。Raid1:至少需要两块盘,磁盘数量是2的n倍,每一块磁盘要有对应的备份盘,利用率是50%,只要有一对磁盘没有损坏就可以正常使用…

    2022年7月15日
    13
  • J1939协议之通俗易懂—-简介

    J1939协议之通俗易懂—-简介J1939简介 J1939协议简介J1939协议是由美国汽车工程师协会(SAE)(SAE协会简介)定义的一组标准。J1939标准用于卡车、公共汽车和移动液压等重型车辆。在许多方面,J1939标准类似于旧版J1708和J1587标准,但J1939标准协议建立在CAN(控制器区域网络,ISO11898)上。物理层(J1939/11)描述了针对客车的电气接口。数据链路层描述…

    2022年5月1日
    315

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号