路径规划算法:基于黄金正弦优化的路径规划算法- 附代码

路径规划算法:基于黄金正弦优化的路径规划算法- 附代码


摘要:本文主要介绍利用智能优化算法黄金正弦算法来进行路径规划。

1.算法原理

黄金正弦算法原理请参考:https://blog.csdn.net/u011835903/article/details/111699194

1.1 环境设定

在移动机器人的路径优化中,每个优化算法的解代表机器人的一条运动路径。优化算法会通过优化计算在众多路径中找出一条最优路径。
优化算法的设定必须和机器人运动环境模型相对应。不失一般性,假设在用栅格法对机器人运动环境建模后得出的结果是 m×n 的矩形区域,坐标值从 1 开始,如图1 。其中坐标原点栅格代表机器人的初始位置,坐标 (m,n)对应的栅格代表机器人的移动目标位置。优化算法设定的一个重要内容是确定优化算法的数学表达形式,在这里这个问题转化为用一个向量表示机器人的移动路径。经过分析发现,尽管栅格法建立的模型对空间进行了离散化,但本质上机器人的移动路径依然是连续的。

在这里插入图片描述

图1.栅格地图

1.2 约束条件

对于机器人的路径优化来说,其运动路径必须局限在栅格空间内,即搜索不能越过栅格的矩形边界。此外,还应受障碍物的限制,即机器人的运动轨迹不能穿过存在障碍物的栅格区域。

1.3 适应度函数

在本文的建模方法中,本文路径规划目标是路径长度最短。路径的长度可以表示为:

L ( P a t h ) = ∑ i = 0 n − 1 ( x l i + 1 − x l i ) 2 + ( y l i + 1 − y l i ) 2 (1) L(Path) = \sum_{i=0}^{n-1}\sqrt{(xl_{i+1} - xl_i)^2 + (yl_{i+1} - yl_{i})^2}\tag{1} L(Path)=i=0n1(xli+1xli)2+(yli+1yli)2 (1)
其中(x,y)是路径中间点的坐标

利用黄金正弦算法对上式进行寻优,找到最短路径。黄金正弦算法参数设定如下:

%% 黄金正弦算法参数设置
dim=length(noLM);%维度,即为非障碍物个数。
numLM0=round((EndPoint(1)-StartPoint(1))/4);%每次迭代选取的的中间路径点个数,可调
lb=0;%下边界
ub=1;%上边界
Max_iteration = 100;%最大迭代次数
SearchAgents_no = 30;%种群数量
fobj = @(x)fun(x,noS,noE,numLM0,net);%适应度函数

2.算法结果

在这里插入图片描述

3.MATLAB代码

本程序中,支持1.地图任意创建保存。2.其实点任意更改。

4.参考文献

[1]罗阳阳,彭晓燕.基于改进PSO的四轮移动机器人全局路径规划[J].计算机仿真,2020,37(07):373-379.

[2]鲁丹. 粒子群算法在移动机器人路径规划中的应用研究[D].武汉科技大学,2009.

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论
以下是基于趋势拟合算法路径规划代码的示例,其中使用了Python的numpy和matplotlib库: ```python import numpy as np import matplotlib.pyplot as plt # 定义趋势拟合函数 def fit_trend(x, y, degree=1): coeffs = np.polyfit(x, y, degree) trend = np.poly1d(coeffs) return trend # 定义路径规划函数 def plan_path(x, y, degree=1, step_size=0.1): # 进行趋势拟合,得到趋势函数 trend = fit_trend(x, y, degree) # 计算路径长度 path_length = 0 for i in range(1, len(x)): dx = x[i] - x[i-1] dy = trend(x[i]) - trend(x[i-1]) path_length += np.sqrt(dx**2 + dy**2) # 计算路径点数 num_points = int(path_length / step_size) # 生成路径点 path_x = np.zeros(num_points) path_y = np.zeros(num_points) path_x[0] = x[0] path_y[0] = trend(x[0]) dist = 0 for i in range(1, num_points): dx = x[-1] - path_x[i-1] dy = trend(x[-1]) - path_y[i-1] dist += np.sqrt(dx**2 + dy**2) target_dist = dist / (path_length / len(x)) idx = np.argmin(np.abs(np.cumsum(np.sqrt(np.diff(x)**2 + np.diff(trend(x))**2)) - target_dist)) path_x[i] = x[idx+1] path_y[i] = trend(path_x[i]) return path_x, path_y # 测试代码 x = np.linspace(0, 10, 100) y = np.sin(x) + np.random.normal(0, 0.1, 100) path_x, path_y = plan_path(x, y) plt.plot(x, y, 'b.') plt.plot(path_x, path_y, 'r-') plt.show() ``` 这段代码中,首先定义了一个`fit_trend`函数来进行趋势拟合,返回一个多项式函数对象。然后,定义了一个`plan_path`函数来进行路径规划。该函数首先利用`fit_trend`函数得到趋势函数,然后根据步长计算路径点数,最后生成路径点。在生成路径点时,利用`np.argmin`函数和`np.cumsum`函数来计算距离,以确定下一个路径点的位置。最后,测试代码生成一个带有随机噪声的正弦曲线,然后利用`plan_path`函数对其进行路径规划,并将原始数据和路径点绘制在图中。

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

智能算法研学社(Jack旭)

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值