用MATLAB来做智能小车的建模与仿真

2023-05-16

两种智能小车的构造简介

在市面上常见的两种智能小车都是基于轮式的,在某宝上面卖的最多的,各位在学生时代拿来应付课程设计和毕业设计用的各种小车分为两种。

1.后轮驱动,前轮阿克曼转向的,通常后轴通过有刷直流电机驱动,前轴通过舵机带动连杆转向。英文通常称为Car-like Model.

2.差速驱动,这种最为多见,可以做成轮式的,也可以做成履带式的。通过每个轮子都有电机,当然你可以说成为轮毂电机驱动。这种通常称为Tank-like Model.

通过在分析控制问题的时候,第一步就是需要对被控对象进行数学模型的建立。

这里我们不强调模型的复杂性,因为大多数都在低速情况行驶,所以可以只考虑运动学模。

如果您需要参加XXX智能车比赛,这类竞速类的项目,就必须考虑动力学模型,且需要通过大量的试验对模型进行参数估计和辨识。

本次只谈论简易的运动学模型,并且进行Simulink环境的模型搭建,最后通过一个仿真环境来验证模型的可用性!

需要的软件与硬件

1.电脑一台

2.装有Matlab/Simulink的软件本次只在仿真环境中建模,不需要小车的硬件和控制器硬件。

Car-like小车建模与仿真

1.模型推导

先给大家一个直观的印象,如在我的上一篇博客中的小车就属于Car-like 小车。您可以跳回上一篇文章进行查看小车图片。


那么对于小车研究运动学模型,是需要对其简化的,可以简化成如下图所示的模型。如果您不想看数学推导,可以直接看最后的推导结果。

2354e22e5300ae0f47545c5dfa7aa58a.png

ee81fcf1167e96c6299a5c1c3135c7bc.png

3d42da62b4e92ebab2748c2f8331cb6b.png

dfc391356b136234b3711f6910e22c74.png

291797cde2b9429af157416e3464de83.png

其中v是小车轴向速度,x,y代表小车后轴中心点的大地坐标。

2.Simulink模型搭建

首先,根据上节中的内容搭建仿真模型,具体模型如下图所示,如果您才入门Simulink,可以按照我的模型自己动手建一遍。

8933f19c432023fafa085356aae5a5c1.png

您会问我为什么会报错呢,原因很简单,我并没有在Baseworkspace或者数据字典定义参数的具体数值,所示编译器无法识别具体参数,为什么我这么做呢,因为会将这部分模型封装成子系统。

在子系统封装界面输入参数。

具体报错原因我点进去几个模块给您看一下。例如这个积分模块

0f3794fb0c53e3dfb5dca776c90a6cc5.png

这个模块输出值代表小车的大地坐标中的横坐标,在Initial Condition中是定义小车的初始值,这里通过参数定义小车的初始横坐标,采样时间用参数表示方便自定义采样时间,不用每次都打开模型输入。

其他方位角,纵坐标,轴距其实都是一样的。输出模块我用了Bus模块,也是方便我们查看自己感兴趣的输出信号。接下来封装子系统。全选界面,右键选择Create Subsystem from Selection。

如下图所示,可以看到,还是存在报错情况,不过没关系,稍后会进行解决。

车身系统目前就建立完了。

56e0364aca84e7bc6e57cf962e54eb95.png

然后我们将进行转向系统和驱动系统的建立,驱动系统的建立我就不用示意了,这时只需要把上次的模型拿过来,接入速度接口就可以了。

输入:期望车速。输出:实际速度。当然控制器不属于系统本身的模型,这里我为了后期方便仿真,就放在一起了。

10660d3993dae9743e76fd44b503a0cd.png

那么转向系统的模型如何建立呢?

小车使用的是一个舵机,舵机具有非常大的惯性,如果用过的同学肯定都知道,这里我使用一个一阶惯性环节来表示转向系统,一阶惯性环节的滞后这里我就不做详细的数学推导了,从bode图上可以看相位有90度的滞后。

这里就用一个开环控制,因为转向舵机我目前无法通过传感器进行反馈控制。时间常数tao的选取决定了系统滞后时间。

不过注意不要将惯性环节和延迟环节进行混淆,千万需要注意!其中Saturation模块是对小车的期望转角做一个限制的,限制最大转角,在模块中输入参数MaxSteerAngle,后期进行封装使用。

1d544fccd7fbef274ab72b4c050e5e43.png

如果您不确定滞后参数,可以加个Scope看一下滞后的时间,从而选定时间常数tao。这时我们将整个系统表示为如下图所示。

f59e7d270e257e6733f257a6c1819603.png

软件的架构一定要清晰,这是最基本的!

接下来创建一个整体的子系统。

80784e8a5b6e1780c817890273c7e603.png

在VehicleBody子系统里,参数没有定义还是处于报错状态。那么右键子系统,选择Mask->Create Mask,出现下图界面:

618432225eb17bb26580bf3bcb825d90.jpeg

先选中Parameters,再点击左侧的Edit。编辑成如下图所示:

b1780228353e68637872658b57422986.png

其中name对应您设置的参数名字,Prompt对应您在封装界面显示的名字。

并且您可以在Documentation条目中进行子系统的描述。最后点击Ok。

再次点击您的子系统,可见已经封装完成。

27650a9c031d120553a05736b1311b2c.png

当然您也可以右键子系统Mask->Add Icon Image选择您小车的图片,我也将自己的小车图片加上去了。

并且您可以在Library中将您的模型生成自定义库。

5471068c6d6881e0df17ceaa18084270.jpeg

3.Car-like小车模型仿真

本次的仿真环境我使用Mathworks培训视频提供的一个demo,在此demo上进行修改。

这个demo的需求是:通过规划好路径,设计纯追踪控制器,使小车能够良好的跟踪轨迹。

7e950c9c3dab98edaa8399f6304c5149.png

emo中分为三个部分:

  • 第一个部分是寻找预瞄点,通过对车辆反馈的当前位置和参考路径位置进行计算。

  • 第二部分就是控制器算法,具体算法这里我就不解释了,在CSDN搜索纯跟踪您可以看到关于此类车辆模型详细的控制器数学推导。

  • 第三部分就是车辆模型部分了。

从图可以看到这个demo使用的Powertrain Blockset工具箱中的3dof自行车动力学模型。

我们需要将此模型替换成我们刚刚建立好的模型用我们的模型进行仿真,其他参数一律不进行修改。

137090cb24f0838187bd6ec229d51d5c.jpeg

已经替换完成,设置小车参数以及车辆在坐标系初始位置。其中SignalProcess子系统里是信号的处理。

4fe6423a1bdd9f41e29275bd3b9e7d99.png

467ac11aa3745e6b9ac27b8372477966.png

ba6241a52bc75a44c25fe23e45fb9a86.png

这时点击仿真。红色点是小车的实际路径,参考路径是蓝色的,可以明显看到,小车在预瞄距离0.3米的情况下,跟踪的很好。

和官方的二自由度小车动力学模型结果几乎误差,这里我就没有把两者的误差进行对比了,大家可以看下轨迹效果。

说明基于运动学建立的模型在低速环境下使用是没有问题的。

e6d567289c7d92d504c0d64671ca5197.png

前轮转角也十分平稳,这里我用的离散模型,采用时间是5e3秒,官方提供的模型是连续的,只不过它在后面的信号处理用了RateTransition模块。

7e6939593e04319ffe43bd36c618baaa.png

Tank-like小车建模与仿真

1.模型推导

Tank-like小车模型较前者简单,并且在某宝的普及率是远远高于带前轮转向的小车的,因为价格相对便宜很多,通过用两块亚克力板和四个轮子加香蕉电机就可以完成组装。下图即为小车简化后的模型。

65430c2f914b732b9cb31b563bebb0c8.png

16144ff55d6dfd7e37ed75bad2d428b2.png

运动学关系就非常清晰了,可以通过速度和角速度求解出x与y的位置,也可以求解出左右的速度分别是多少。那么下一步就可以在Simulink中建模了。

2.Simulink模型搭建

建模也十分简单,将上面推导的公式代入,就可以得到一系列的输出了,这里我输出了8个值,分别是小车的x和y的位置信息,x和y的速度信息,方位角和角速度信息和两个车辆不同的前进速度!

b8478bce5f3c96e1afd5ff9c3f136152.png

注意在高亮的模块中,是需要和上个模型一样在模块中定义初始值和采样时间的,最后用总线信号输出。然后按照上一个模型类似步骤的简单封装,就可以做成如下图所示的样子。

e003bfb7e4b0b4f00fe0075b27db96e3.png

当您需要输出什么信号的时候,后接一个Bus Selector即可选择您想要的任何信号。

3.控制器推导

3d09c55ffbe966766d6ff811d1bd024e.png

2a659c94de5bc74309191913d6d54011.png

当然您可以发现我们需要的距离 x xx是在小车坐标系下的,而预瞄距离是在大地坐标系下的,所以在求出预瞄点坐标的时候,您需要将预瞄点坐标转换到小车的坐标系中,而坐标系按如下方法进行转换的:

(1)首先将预瞄点坐标平移到以小车坐标为原点的坐标系中。

(2)将平移后的坐标系旋转小车的方位角θ ,得到最终的小车坐标系。

4.仿真环境搭建

参考轨迹按照上个模型的仿真环境不变,计算预瞄点距离的算法子系统不变,我们需要改变的是控制器子系统和车辆模型子系统。

f648724797b7e3612f2ddab7676bc54c.png

1e4cba34e3a48d0620f2c6e0bff77a70.png

再看控制器如何设计的,最重要的就是坐标系的转换(蓝色区域),建议有兴趣的同学可以自己推算一遍,十分简单的数学几何关系。

这里就没有使用参考轨迹的方位角信息了!

4b327ce297803b5abbccbdc9c61f1608.png

整体模型架构如下:

4f3ca807b952538adfb38422df917f42.png

本次仿真设置的是预瞄距离20cm,车速0.5m/s,开始仿真。看一下期望轨迹。

32611d9823a468f1c927e61526eaa374.png

这是差速小车,当然您可以将左轮速度和右轮速度输出看一下。总体来说弄懂原理十分简单。参数的具体值就不对比了。

如果您有兴趣可以从头自己搭建一个您自己期望的参考轨迹,根据已经设计好的控制器进行仿真看一下效果。将x和y坐标接入xyGraph模块一样可以看到实际轨迹。

总结

本次进行了两种常见的小车模型推导,建模与仿真,其中Tank-like小车您也可以沿用成履带式的小车,在低速有约束的情况下都是模型都是通用的。

当您从某宝等电商平台购买回一台小车,我建议您从小车的构造本质进行入手,才能最其进行更好的控制,控制首先需要对被控对象进行深入的了解,无论是正向原则和数据驱动的方式建模,都需要对模型进行搭建,才能展开后续控制器的设计。

如果您使用的是Matlab2019b,那您可以在机器人工具箱中找到我今天做的两种小车模型和差速小车纯跟踪控制器,区别在于工具箱中的模块是基于面向对象编程的方法,没有用Simulink模块搭建,底层代码是不可见的。

您可以验证下本次教程的模型搭建是否和官方的模型有出入。

31991e275f61925f03c62fc91366c516.jpeg

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

用MATLAB来做智能小车的建模与仿真 的相关文章

  • 在Matlab中使用中心切片定理实现滤波反投影算法

    我正在研究一种使用中心切片定理的滤波反投影算法作为家庭作业 虽然我理解纸上的理论 但在 Matlab 中实现它时遇到了问题 我得到了一个可以遵循的框架 但我认为我可能误解了一个步骤 这是我所拥有的 function img sampleFB
  • 如何将向量标准化/非标准化到范围 [-1;1]

    我怎么能够正常化到范围的向量 1 1 我想使用函数norm 因为它会更快 也让我知道我该怎么做非规范化之后的向量正常化 norm对向量进行归一化 使其平方和为 1 如果要对向量进行归一化 使其所有元素都在 0 和 1 之间 则需要使用最小值
  • Matlab 中二维插值的函数形式

    我需要从二维数据数组构造一个插值函数 我需要返回实际函数的东西的原因是 我需要能够将函数作为我需要进行数值积分的表达式的一部分进行计算 因此 interp2 并没有解决这个问题 它不返回函数 我可以使用 TriScatteredInterp
  • 朴素分类器 matlab

    在 matlab 中测试朴素分类器时 即使我在相同的样本数据上进行训练和测试 我也会得到不同的结果 我想知道我的代码是否正确 是否有人可以帮助解释这是为什么 dimensionality reduction columns 6 U S V
  • 在 MATLAB 中高效地形成动态窗口

    有人可以帮助我提供一种有效的方法 或者帮助我执行提供的代码 以尽可能少的步骤获得相同的结果 我将感激你 我有一个原始数组 A 1 1 1 4 3 4 5 4 4 3 3 1 0 0 2 6 2 6 3 6 6 2 7 4 8 7 2 2 2
  • MATLAB 图像锐化 - 使用(1-高斯低通滤波器)的高斯高通滤波器

    我试图通过设计高斯高通滤波器来锐化图像 我想利用高通滤波器相当于单位矩阵减去低通滤波器的事实来做到这一点 所以我执行了以下操作 image imread Question3 Data Cats jpg read image H 1 fspe
  • 如何在 MATLAB 中创建带有分级标记颜色的散点图?

    我想在 MATLAB 中绘制一个简单的散点图 标记颜色从光谱的一端到另一端各不相同 例如红色 橙色 黄色 蓝色 紫色 我的数据比较了一段时间内河流的水量和水质 3 个简单的列 时间 水量 质量 我想绘制数量与质量的 x y 散点图 但颜色随
  • MATLAB 中元胞数组的左连接

    I ve 2 cellMATLAB 中的数组 例如 A jim 4 paul 5 sean 5 rose 1 第二个 B jim paul george bill sean rose 我想做一个 SQL 左连接 这样我就可以得到 B 中的所
  • 将输出从符号数学 (sym) 转换为浮点型

    我的问题类似于这个问题 https stackoverflow com questions 11114101 how to convert mupad symbol i sqrt 1 to i in matlab 11114959 1111
  • 如何将Matlab代码库与Android集成?

    我有一个算法和 MATLAB 中的一些其他代码 我想在我的 Android 应用程序中使用它 我怎样才能做到这一点 我可以从 MATLAB 制作 jar 文件以便与 Android 一起使用吗 我必须做点别的事吗 如果您还有其他产品 适用于
  • MATLAB 中的抗锯齿线和标记

    您好 我在 MATLAB 中有一张图像 我希望这条线是平滑的 看看从 0 4 到 0 8 的线 这太可怕了 当在图中使用 LineSmoothing on 运算符时 我得到了这个 我在线条上做得很好 但它也使标记变得平滑 而且它们太可怕了
  • Microsoft Visual C++ 2008 和 R2007b 的 Mex 类型

    我想对 vs2008 和 matlab2007b 使用 mex 类型 我尝试了下面的代码 include
  • 查找数组中元素之间的平均差异的有效方法

    希望标题不会让人困惑 通过例子来展示很简单 我有一个像这样的行向量 1 5 6 我想找到每个元素之间的平均差异 此示例中的差异为 4 和 1 因此平均值为 2 5 这是一个小例子 我的行向量可能非常大 我是 MatLab 新手 那么有没有一
  • 如何以编程方式指定 MATLAB 编辑器键绑定

    我想将键盘键绑定设置为Windows 默认设置我想在启动时使用startup m因为我希望在大量系统上设置此设置 首选项对话框中的等效设置是 MATLAB gt Keyboard gt Shortcuts gt Active Setting
  • 是否有一个函数可以检查矩阵是否对角占优(行占优)

    矩阵是对角占优 http en wikipedia org wiki Diagonally dominant matrix 按行 如果对角线处的值在绝对意义上大于该行中所有其他绝对值的总和 对于列也是如此 只是相反 matlab中有没有函数
  • 如何从绘图处理程序中绘图?

    我有绘图的处理程序或图形的处理程序 例子 h plot 1 0 2 10 xx get h xx DisplayName Annotation 1x1 handle Color 0 0 1 LineStyle LineWidth 0 500
  • 在 MATLAB 中创建共享库

    一位研究人员在 MATLAB 中创建了一个小型仿真 我们希望其他人也能使用它 我的计划是进行模拟 清理一些东西并将其变成一组函数 然后我打算将其编译成C库并使用SWIG https en wikipedia org wiki SWIG创建一
  • MATLAB 图中轴标签与轴之间的距离

    我正在使用 MATLAB 绘制一些数据 我想调整轴标签与轴本身之间的距离 但是 只需向标签的 位置 属性添加一点即可使标签移出图窗窗口 是否有 保证金 属性或类似的东西 在上图中 我想增加数字和标签 Time s 之间的距离 同时自动扩展数
  • 作为动画的八度情节点

    我有以下八度脚本 TOTAL POINTS 100 figure 1 for i 1 TOTAL POINTS randX rand 1 randY rand 1 scatter randX randY hold on endfor 当我运
  • 定义自定义 Mupad 程序的一般相对搜索路径

    假设我有一个 mupad 笔记本myMupadNotebook mn在路径上 C projectFolder ABC abc 它调用程序MyMupadProcedure mu它位于 C DEF GHI 现在我有一个 Matlab 脚本mai

随机推荐

  • VSCode的C/C++扩展功能

    VSCode的C C 43 43 扩展功能 1 在 Linux 上 使用 C 43 43 1 1 创建 Hello World1 2 探索 IntelliSense1 3 构造 helloworld cpp1 3 1 运行 build1 3
  • 从源码理解智能指针(二)—— shared_ptr、weak_ptr

    目录 计数器 Ref count Ref count del Ref count del alloc Ptr base Ptr base的成员变量 构造函数 赋值重载 获取引用计数 减少引用计数 Reset函数 Resetw函数 share
  • muduo源码学习(1):异步日志——日志消息的存储及输出

    目录 前言 日志存储的实现 日志输出的实现 总结 前言 muduo中的日志 xff0c 是诊断日志 用于将代码运行时的重要信息进行保存 xff0c 方便故障诊断和追踪 日志一般有两种 xff0c 一种是同步日志 xff0c 一种是异步日志
  • muduo源码学习(2):异步日志——异步日志的实现

    目录 什么是异步日志 异步日志的实现 前端与后端 前端与后端的交互 资源回收 后端与日志文件 滚动日志 自动flush缓冲区 开启异步日志功能 总结 在前文中分析了日志消息的存储和输出 xff0c 不过并没有涉及到异步日志 xff0c 下面
  • muduo异步日志——core dump后查找还未来得及写出的日志

    目录 前言 生成core文件 gdb调试Core文件 前言 通过异步日志的实现可以知道 xff0c 日志消息并不是生成后立刻就会写出 xff0c 而是先存放在前端缓冲区currentBuffer或者前端缓冲区队列buffers中 xff0c
  • C++知识积累:成员函数运算符重载与非成员函数运算符重载

    运算符重载 xff0c 是C 43 43 多态的表现形式之一 xff0c 可以通过对运算符进行重载来实现运算符特定的功能 运算符重载一般具有以下原则 xff1a xff08 1 xff09 不可重载不存在的运算符 xff0c 如重载 来表示
  • (二叉树)二叉树的最近公共祖先

    题目描述 给定一个二叉树 找到该树中两个指定节点的最近公共祖先 百度百科中最近公共祖先的定义为 xff1a 对于有根树 T 的两个结点 p q xff0c 最近公共祖先表示为一个结点 x xff0c 满足 x 是 p q 的祖先且 x 的深
  • 有符号数、无符号数理解

    大家都知道 xff0c 在C C 43 43 中 xff0c 对于w位编译器 xff0c 其有符号数表示的数值范围为 2 w 1 2 w 1 1 xff0c 无符号数表示的数值范围为0 2 w 1 xff0c 举个例子 xff0c 在16位
  • ​PCB的 “ 黑科技 ” ,应该是这个。

    大家好 xff0c 我是张巧龙 xff0c 前段时间炒的很火的折叠屏手机不知道大家还记得不 xff1f 折叠屏手机之所以这么具有 34 韧性 34 xff0c 全靠背后的柔性电路板 FlexiblePrintedCircuit xff0c
  • 指针数组、数组指针——用指针访问数组方法总结

    目录 1 数组元素的访问 2 通过指针访问数组 2 1 通过指针访问一维数组 2 2 通过指针访问二维数组 2 2 1 指向元素的指针 2 2 2 指向每一行的指针 xff08 指针数组方式 xff09 2 2 3 指向整个数组的指针 xf
  • C++知识积累:如何获取虚函数表以及虚函数地址

    如果一个类中存在虚函数的话 xff0c 那么编译器就会为这个类生成一个虚函数表 xff0c 这个虚函数表中按照个虚函数的声明顺序存放了各个虚函数的地址 xff0c 需要注意的是 xff0c 这个虚函数表并不存在于类中 xff0c 而对于这个
  • C++多线程:互斥锁

    目录 1 前言 2 互斥锁 2 1 互斥锁的特点 2 2 互斥锁的使用 2 2 std lock guard 3 死锁 3 1 死锁的含义 3 2 死锁的例子 3 3 死锁的解决方法 1 前言 比如说我们现在以一个list容器来模仿一个消息
  • Linux下MySQL中文显示问号乱码问题解决

    本文主要针对于Linux下MySQL插入中文数据显示问号的问题 网上一种普遍使用的方法是修改 etc my cnf文件 xff08 我的这个文件位于 etc mysql my cnf xff09 xff0c 修改步骤如下 xff1a 1 在
  • Linux下更改文件权限

    目录 查看文件权限 修改文件权限 查看文件权限 查看文件权限可以通过ls l命令查看 xff0c 如下所示 xff1a 如果只想查看某一个文件的权限 xff0c 可以使用grep xff0c 如下所示 xff1a 可以发现 xff0c 每一
  • 【蓝桥杯算法提高VIP-开灯游戏(两种超易理懂解法:暴力/位操作(切换位))(纯正C语言代码)】

    蓝桥杯算法提高VIP 开灯游戏 题目描述 有9盏灯与9个开关 xff0c 编号都是1 9 每个开关能控制若干盏灯 xff0c 按下一次会改变其控制的灯的状态 亮的变成不亮 xff0c 不亮变成亮的 具体如下 xff1a 第一个开关控制第二
  • C语言二叉查找树(图文详解)(超详细)

    二叉查找树 本人在第一次学习二叉树的时候 感觉很懵懵懂懂 勉强知道了二叉树的结构和查找方式 但要我自己去动手写的时候 可是难上加难 所以这里我用干货 43 实际例子的方式让你上手二叉树 这个例子几乎可以套用到所有链式结构问题 我们开始吧 文
  • STM32F4中断(Interrupt)详解

    STM32F4中断 Interrupt 详解 文章目录 STM32F4中断 Interrupt 详解一 中断是什么 二 STM32的中断体系2 1 STM32的中断分类2 2 STM32中断优先级2 2 1 抢占优先级2 2 2 响应优先级
  • 基于RK3399&ESP8285自动售货柜项目—ESP8266(8285)程序编写与烧录

    基于RK3399 amp ESP8285自动售货柜项目 ESP8266 8285 程序编写与烧录 本系列文章讲详细讲解该基于RK3399及ESP8285自动售货柜的完整实现方法 xff0c 从硬件连接到网络通信再到软件实现 xff0c 本产
  • Windows10安装OpenCV4.1.0+opencv_contrib

    Windows10安装OpenCV4 1 0 43 opencv contrib 文章目录 Windows10安装OpenCV4 1 0 43 opencv contrib一 Visual Studio 2015安装二 下载和安装OpenC
  • 用MATLAB来做智能小车的建模与仿真

    两种智能小车的构造简介 在市面上常见的两种智能小车都是基于轮式的 xff0c 在某宝上面卖的最多的 xff0c 各位在学生时代拿来应付课程设计和毕业设计用的各种小车分为两种 1 后轮驱动 xff0c 前轮阿克曼转向的 xff0c 通常后轴通