基于模糊控制算法的倒立摆控制系统matlab仿真

news/2024/9/18 3:53:56/文章来源:https://www.cnblogs.com/hardware/p/18357970

1.课题概述

        基于模糊控制算法的倒立摆控制系统,模糊规则,模糊控制器等通过MATLAB编程实现,通过模糊控制器对小车倒立摆平衡系统进行控制,输出倒立摆从不稳定到稳定的动画过程,最后输出小车,倒立摆的收敛过程。

 

2.系统仿真结果

 

 

 

3.核心程序与模型

版本:MATLAB2022a

 

for ij=Ts:Ts:Tend% 计算小车位置误差  ex  = xp - X(1);            % 计算小车位置误差变化量dex = Xt - X(1);      % 计算摆角度误差  eq  = qp - X(3);     % 计算摆角度误差变化量 deq = Qt - X(3);   % 模糊控制(u_x:小车位置,u_q:摆角度)tmp1    = func_check(ex/xn, dex/xdn);tmp2    = func_check(eq/qn, deq/qdn);Fr_tmp1 = func_fuzzy_rule(tmp1);Fr_tmp2 = func_fuzzy_rule(tmp2);Fr_dtmp1= func_defuzzy(Fr_tmp1);Fr_dtmp2= func_defuzzy(Fr_tmp2);u_x     = Fr_dtmp1 * Un*1;u_q     = Fr_dtmp2 * Un*2;Uset(ic,:) = [u_x, u_q];% 保存控制量  % 计算外力F  F       = -u_x + u_q;% 更新上一次的位置和角度  Xt      = X(1);Qt      = X(3);% 使用ode45求解摆的下一个状态 [T, X_next] = ode45(@func_pendulum, [0, Ts], X);% 更新当前状态X  X       = X_next(end,:)';% 保存当前状态到X_Fuzzy  Xf(ic,:)= X';% 保存当前时间到time_Fuzzy Tf(ic)  = ij;% 保存当前外力F到F_save  Fset(ic)= F;% 更新计数器  ic      = ic + 1;   
endfigure
plot(Time_result, X_result(:,1:2),'linewidth',2)
grid on
xlabel('Time [s]')
legend('小车位置[m]','小车速度[m/s]')figure
plot(Time_result, X_result(:,3:4)*180/pi,'linewidth',2)
grid on
xlabel('Time [s]')
legend('摆角', '摆锤速度')figure
plot(Time_result, Fset,'b','linewidth',2)
hold on
plot(Time_result, Uset(:,1),'r','linewidth',2)
hold on
plot(Time_result, Uset(:,2),'m','linewidth',2)
grid on
xlabel('Time [s]')
ylabel('力[N]')
legend('输入F', 'U_x', 'U_q')
0006

  

 

4.系统原理简介

       倒立摆是一个经典的控制问题,其目标是保持摆在垂直位置。由于其非线性、不稳定特性,传统控制方法往往难以实现理想效果。模糊控制作为一种先进的控制策略,能够处理不确定性和非线性问题,因此在倒立摆控制中具有显著优势。

 

4.1. 模糊控制算法原理

        模糊控制是基于模糊集合理论、模糊语言变量及模糊逻辑推理的知识表示和推理方式。它模仿人的模糊思维方式和决策过程,不需要建立精确的数学模型。其核心组成包括:

 

模糊化:将输入量转化为模糊量,用隶属度函数表示。

μA(x)

 

其中,A 为模糊集合,x 为具体值,μA(x) 表示 x 对 A 的隶属度。

 

规则库:根据专家经验或数据,建立模糊规则,如“IF-THEN”规则。

 

模糊推理:基于规则库和当前模糊输入,进行推理得到模糊输出。

 

去模糊化:将模糊输出转化为精确量,常见的去模糊化方法有最大隶属度法、重心法等。

 

4.2. 倒立摆控制系统设计

        对于倒立摆,我们定义其角度为 θ,角速度为 ω。控制目标是使得 θ 接近 0。设计模糊控制器如下:

 

输入:e=θ−θd e = θ - θ_de=θ−θd(角度误差)和 ec=ω−ωd ec = ω - ω_dec=ω−ωd(角速度误差)。

 

输出:u(控制力)。

 

模糊化:e 和 ec 的论域为 [-3, 3],u 的论域为 [-1, 1]。选择三角形或高斯型隶属度函数。

 

规则库:根据经验和试验,建立如下规则(仅举例):

 

如果 e 是负大且 ec 是负大,则 u 是正大。

如果 e 是零且 ec 是零,则 u 是零。

如果 e 是正大且 ec 是正大,则 u 是负大。

 

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.hqwc.cn/news/783513.html

如若内容造成侵权/违法违规/事实不符,请联系编程知识网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

第46届金砖国家世界技能大赛 内存取证样题一

第46届金砖国家世界技能大赛 内存取证样题一题目:从内存中获取到用户admin的密码并且破解密码,以Flag{admin,password} 形式提交(密码为 6 位); 获取当前系统 ip 地址及主机名,以 Flag{ip:主机名}形式提交; 获取当前系统浏览器搜索过的关键词,作为 Flag 提交; 当前系统中…

TS 项目如何从 ESlint8 升级至 9并集成 Prettier

eslint 8 到 9 属于破坏性更新(Break Change),因此导致 eslint 8 配置方式无法直接使用(可以使用兼容包,但这不是本文的主题)。 其实大家最关心的就是从 eslint 8 到 9 之后的写法,而与 eslint 息息相关的多种配置插件也需要大量的变更,因此本文核心就是:使用 eslint9…

Spring Cloud Gateway 实现简单自定义过滤器

背景 Spring Cloud Gateway 是 Spring Cloud 退出的第二代网关框架,我们可以用它来实现 反向代理,路由转发,权限校验等功能,这里介绍一个它的基础功能,通过 Filter 机制实现一个简单的 HTTP 接口处理。 从总体上来看 Spring Cloud Gateway 提供的过滤器可以分为两类,一种…

【Java手写RPC框架系列-1】—— 基础知识准备:RPC+Netty

代码随想录知识星球介绍 https://articles.zsxq.com/id_m76jd72243bi.html 基于Netty手写实现RPC https://www.cnblogs.com/mic112/p/15565795.html项目背景与介绍RPC:远程过程调用协议:客户端在不知道调用细节的情况下,调用存在于远程计算机上的某个对象,就像调用本地应用…

非线性规划的经典例题--选址问题

本章会介绍如何利用非线性规划解决选址问题,这个问题是文章线性规划在数学建模中的两道例题中第二道投料问题的第二小题,本章为基于这道题的基础上进行介绍,建议读者返回去看一看 目录一、问题提出二、问题分析三、模型建立四、代码实现1.输入目标函数2.输入线性约束 一、问…

了解10X文库组成

转的文章 版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。 原文链接:https://blog.csdn.net/LittleComputerRobot/article/details/137207083R1: 26 表示10X barcode 的 16bp碱基 + 10bp UMI; i7: 8表示 8bp 样本index序列 Read…

Java数组07:稀疏数组

1. 线性结构线性结构是最常用的数据结构,其特点是数据元素之间存在一对一的线性关系。 线性结构有两种不同存储结构,即顺序存储结构和链式存储结构。 顺序存储的线性表称为顺序表,顺序表中的存储元素是连续的,即在内存中是连续的,例如数组。 链式存储的线性表称为链表,链表…

【Django开发】django美多商城项目完整开发4.0第2篇:项目准备【附代码文档】

本教程的知识点为:美多商城 项目准备 项目准备 配置 1. 修改settings/dev.py 文件中的路径信息 2. INSTALLED_APPS 3. 数据库 用户部分 图片验证码 1. 后端接口设计: 视图原型 2. 具体视图实现 用户部分 使用Celery完成发送短信 判断帐号是否存在 1. 判断用户名是否存在 后端…

IntelliJ IDEA 单元测试插件 TestMe

IntelliJ IDEA 添加单元测试插件 TestMe

Redis 使用场景

Redis 使用场景前三种使用广泛,另加队列用

Datawhale X 魔搭 AI夏令营 AIGC方向 Task2

代码逐行解析 先记录一下baseline代码通义千问生成的逐行解析: # 安装 Data-Juicer 和 DiffSynth-Studio !pip install simple-aesthetics-predictor # 安装simple-aesthetics-predictor !pip install -v -e data-juicer # 安装data-juicer !pip uninstall pytorch-lightning …

[Paper Reading] Multiple View Geometry Transformers for 3D Human Pose Estimation

Multiple View Geometry Transformers for 3D Human Pose Estimation link 时间:CVPR2024 机构:University of Toronto && Southeast University && Microsoft Research Asia TL;DR 提出一种基于Transformer端到端3D Human Pose Estimation方法MVGFormer,核…