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

发布时间:2023-07-27  

两种智能小车的构造简介

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


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

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

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

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

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

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


需要的软件与硬件

1.电脑一台

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

Car-like小车建模与仿真

1.模型推导

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


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

98b1f622-ee6e-11ed-90ce-dac502259ad0.png

98c72998-ee6e-11ed-90ce-dac502259ad0.png

98df9488-ee6e-11ed-90ce-dac502259ad0.png

98f5f98a-ee6e-11ed-90ce-dac502259ad0.png

9922523c-ee6e-11ed-90ce-dac502259ad0.png

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

2.Simulink模型搭建

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

9946eac0-ee6e-11ed-90ce-dac502259ad0.png

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

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

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

995b6860-ee6e-11ed-90ce-dac502259ad0.png

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

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

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

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

997aece4-ee6e-11ed-90ce-dac502259ad0.png

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

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

9992d160-ee6e-11ed-90ce-dac502259ad0.png

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

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

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

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

99adc2a4-ee6e-11ed-90ce-dac502259ad0.png

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

99c0b60c-ee6e-11ed-90ce-dac502259ad0.png

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

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

99d374e0-ee6e-11ed-90ce-dac502259ad0.png

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

99ebbc62-ee6e-11ed-90ce-dac502259ad0.png

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

9a17b5ce-ee6e-11ed-90ce-dac502259ad0.png

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

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

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

9a29cde0-ee6e-11ed-90ce-dac502259ad0.png

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

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

9a45a45c-ee6e-11ed-90ce-dac502259ad0.jpg

3.Car-like小车模型仿真

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

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

9a642d82-ee6e-11ed-90ce-dac502259ad0.png

emo中分为三个部分:

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

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

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

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

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

9a79e4ce-ee6e-11ed-90ce-dac502259ad0.png

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

9a975bd0-ee6e-11ed-90ce-dac502259ad0.png

9aabf1bc-ee6e-11ed-90ce-dac502259ad0.png

9ac6aa48-ee6e-11ed-90ce-dac502259ad0.png

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

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

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

9ada09da-ee6e-11ed-90ce-dac502259ad0.png

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

9aedc420-ee6e-11ed-90ce-dac502259ad0.png

Tank-like小车建模与仿真

1.模型推导

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

9b0940ce-ee6e-11ed-90ce-dac502259ad0.png

9b24c51a-ee6e-11ed-90ce-dac502259ad0.png

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

2.Simulink模型搭建

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

9b3c278c-ee6e-11ed-90ce-dac502259ad0.png

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

9b544678-ee6e-11ed-90ce-dac502259ad0.png

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

3.控制器推导

9b6db8e2-ee6e-11ed-90ce-dac502259ad0.png

9b82f31a-ee6e-11ed-90ce-dac502259ad0.png

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

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

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

4.仿真环境搭建

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

9ba12a42-ee6e-11ed-90ce-dac502259ad0.png

9bbebe4a-ee6e-11ed-90ce-dac502259ad0.png

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

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

9be355b6-ee6e-11ed-90ce-dac502259ad0.png

整体模型架构如下:

9bf64b76-ee6e-11ed-90ce-dac502259ad0.png

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

9c0fee1e-ee6e-11ed-90ce-dac502259ad0.png

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

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

总结

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

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

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

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

9c235526-ee6e-11ed-90ce-dac502259ad0.png



文章来源于:电子工程世界    原文链接
本站所有转载文章系出于传递更多信息之目的,且明确注明来源,不希望被转载的媒体或个人可与我们联系,我们将立即进行删除处理。

我们与500+贴片厂合作,完美满足客户的定制需求。为品牌提供定制化的推广方案、专属产品特色页,多渠道推广,SEM/SEO精准营销以及与公众号的联合推广...详细>>

利用葫芦芯平台的卓越技术服务和新产品推广能力,原厂代理能轻松打入消费物联网(IOT)、信息与通信(ICT)、汽车及新能源汽车、工业自动化及工业物联网、装备及功率电子...详细>>

充分利用其强大的电子元器件采购流量,创新性地为这些物料提供了一个全新的窗口。我们的高效数字营销技术,不仅可以助你轻松识别与连接到需求方,更能够极大地提高“闲置物料”的处理能力,通过葫芦芯平台...详细>>

我们的目标很明确:构建一个全方位的半导体产业生态系统。成为一家全球领先的半导体互联网生态公司。目前,我们已成功打造了智能汽车、智能家居、大健康医疗、机器人和材料等五大生态领域。更为重要的是...详细>>

我们深知加工与定制类服务商的价值和重要性,因此,我们倾力为您提供最顶尖的营销资源。在我们的平台上,您可以直接接触到100万的研发工程师和采购工程师,以及10万的活跃客户群体...详细>>

凭借我们强大的专业流量和尖端的互联网数字营销技术,我们承诺为原厂提供免费的产品资料推广服务。无论是最新的资讯、技术动态还是创新产品,都可以通过我们的平台迅速传达给目标客户...详细>>

我们不止于将线索转化为潜在客户。葫芦芯平台致力于形成业务闭环,从引流、宣传到最终销售,全程跟进,确保每一个potential lead都得到妥善处理,从而大幅提高转化率。不仅如此...详细>>