[算法]-最短路径算法总结「建议收藏」

[算法]-最短路径算法总结「建议收藏」Dijkstra最短路径算法按路径长度的递增次序,逐步产生最短路径的贪心算法基本思想:首先求出长度最短的一条最短路径,再参照它求出长度次短的一条最短路径,依次类推,直到从顶点v到其它各顶点的最短路径全部求出为止。时间复杂度为O(n2)算法流程:首先选定源点1,建立邻接矩阵C[5][5],初始化三个数组分别为D[n],P[n],S[n],分别用来存储从源点到对应点的最短距离和最短路…

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

Dijkstra最短路径算法

按路径长度的递增次序,逐步产生最短路径的贪心算法

基本思想:首先求出长度最短的一条最短路径,再参照它求出长度次短的一条最短路径,依次类推,直到从顶点v 到其它各顶点的最短路径全部求出为止。

时间复杂度为O(n2)

算法流程:
在这里插入图片描述

  1. 首先选定源点1,建立邻接矩阵C[5] [5],初始化三个数组分别为D[n],P[n],S[n],分别用来存储从源点到对应点的最短距离和最短路径上到该序列节点的前序节点,以及是否已经找到最短路径。
  2. 开始遍历源点的出边,将邻接矩阵的第一行放入数组D,找出距离最小的节点序号2,将数组P中的P[2]置为1,因为2的前序节点为1。
  3. 以上一步骤找到的节点为起点,继续遍历其邻接点(此处为2),若C[1][2]+C[2][m]<D[m] 则将其替换进数组D中,并将数组P中的P[m]置为2,因为m最短路径上的前序节点为2,节点2的邻接点全部遍历完成后,再从数组D中找出值最小且未被选中过的节点
  4. 重复以上步骤3,直到所有点都已经加入到数组S中。
    在这里插入图片描述

程序实现:

public class Dijikstra extends Strategy {

  // open表
  Map<Integer, Node> open = new HashMap<>();
  // close表
  Map<Integer, Node> close = new HashMap<>();
  // 动作列表
  int[][] motion = get_motion();

  @Override
  public List<Location> nextstep(Location start, Location work) {
    Node startnode = new Node(start.getX(), start.getY(), 0, -1);
    Node worknode = new Node(work.getX(), work.getY(), 0, -1);

    // 将起点放入open表
    open.put(cal_grid_index(startnode), startnode);
    // 向起点的上下左右四个方向扩展
    while (true) {

      // 与Astar算法不同 A*算法含有启发式函数 可以保证更快找到目标点
      // 但A*算法有可能永远无法找到目标点
      // Dijikstra算法虽然较慢 会遍历全部方向的点 但一定可以找到一条到目标点的路径

      // 找到map中节点cost最小的元素
      int temp_cost = Integer.MAX_VALUE;
      int current_id = 0;
      Node current = null;
      for (Map.Entry<Integer, Node> entry : open.entrySet()) {
        int current_cost =
            entry.getValue().cost + cal_Manhattan_distance(entry.getValue(), worknode);
        if (current_cost < temp_cost) {
          current_id = entry.getKey();
          temp_cost = current_cost;
          current = entry.getValue();
        }
      }
      // 判断是否找到目标点
      if (current.x == work.getX() && current.y == work.getY()) {
        System.out.println("找到目标点!");
        worknode.pind = current.pind;
        worknode.cost = current.cost;
        break;
      }
      // 移除open表中的current节点
      for (Iterator<Integer> iterator = open.keySet().iterator(); iterator.hasNext();) {
        Integer key = iterator.next();
        if (key == current_id) {
          iterator.remove();
        }
      }

      // 将current节点放入close表中
      close.put(current_id, current);
      for (int i = 0; i < motion.length; i++) {
        Node node = new Node(current.x + motion[i][0], current.y + motion[i][1],
            current.cost + motion[i][2], current_id);
        int node_id = cal_grid_index(node);
        // 该节点不能被扩展
        if (!verify_node(node)) {
          continue;
        }
        // 该节点已被确认过
        if (close.containsKey(node_id)) {
          continue;
        }
        // 将新扩展节点放入open表中
        if (!open.containsKey(node_id)) {
          open.put(node_id, node);
        } else {
          // 更新open表
          if (open.get(node_id).cost > node.cost) {
            open.put(node_id, node);
          }
        }
      }
    }
    open.clear();
    return cal_final_path(worknode, close);
  }

  /**
   * 根据终点回溯计算最短路径
   * 
   * @param worknode
   * @param close
   * @return
   */
  private List<Location> cal_final_path(Node worknode, Map<Integer, Node> close) {
    List<Location> ans = new ArrayList<>();
    int pind = worknode.pind;
    int i = 0;
    while (pind != -1) {
      // System.out.println(i);
      Node node = close.get(pind);
      Location location = new Location(node.x, node.y);
      ans.add(location);
      pind = node.pind;
    }
    return ans;
  }

  /**
   * 计算曼哈顿距离
   * 
   * @param now
   * @param end
   * @return
   */
  private int cal_Manhattan_distance(Node now, Node end) {
    return Math.abs(now.x - end.x) + Math.abs(now.y - end.y);
  }

  /**
   * 判断节点是否合理 1. 是否超出范围 2. 是否遇到货柜
   * 
   * @param node
   * @return
   */
  private boolean verify_node(Node node) {
    Container[][] map = ContainerMap.getMap();
    if (node.x < 0 || node.x >= ContainerMap.N || node.y < 0 || node.y >= ContainerMap.N) {
      return false;
    }
    if (map[node.x][node.y] != null) {
      return false;
    }
    return true;
  }

  /**
   * 计算栅格地图中的线性坐标
   * 
   * @param node
   * @return
   */
  private int cal_grid_index(Node node) {
    return node.y * ContainerMap.N + node.x;
  }

  /**
   * 定义动作及其代价 motion = {dx, dy, cost}
   * 
   * @return
   */
  private int[][] get_motion() {
    int[][] motion = {
  
  {1, 0, 1}, {0, 1, 1}, {-1, 0, 1}, {0, -1, 1}};
    return motion;
  }
}

A*算法

利用A*算法在拓扑地图上实现的路径规划:
关键:代价函数

  • 对于任意节点n
    • g(n)=从树根到n的代价
    • h*(n)=从n到目标节点的优化路径的代价
    • f*(n)=g(n) + h*(n)是节点n的代价
  • 估计h*(n)
    • 使用任何方法去估计h*(n), 用h(n)表示h*(n)的估计
    • h(n) <= h*(n)总为真
    • f(n)=g(n)+h(n) <= g(n)+h*(n)=f*(n)定义为n的代价

在这里插入图片描述
在这里插入图片描述
在这里插入图片描述
在这里插入图片描述

利用A*算法在栅格地图上实现的路径规划:
A*算法的核心思想是通过设置代价函数对未遍历点进行筛选,使搜索方向更偏向目标点,缩小搜索范围,进而提升搜索速度。公式f(n)=g(n)+h(n)是路径成本表达式。其中,f(n)为起点P—遍历中间点n—目标点g的路径总成本;g(n)为起点P—遍历中间点n之间的实际路径成本;h(n)为遍历中间点n—目标点Q的预估路径成本。寻找到最优解的条件是f(n)尽可能接近起点P到目标点Q的真实路径成本,此时搜寻到的点都是最短路径的过程点。由于遍历中间点n的g(n)值己知,因此只能提升h(n)的估计准确性,根据AGV在作业环境中只有“上下左右”四个行走方向,所以选取 h(n) = |Xn – XQ|+|Yn – YQ| (即曼哈顿距离),
算法流程描述如下:
(1)新建两个列表用作路径数据点存储:open list和close list,open list初始默认包含路径起点;
(2)计算与起点相邻的节点的f(n),按数值大小进行排序,f(n)最小的节点m转移到close list作为当前最优路径过程点,其余点留在open list中做备选点。
(3)若m为目标点,则停止搜索,并输出最优路径结果,否则继续;
(4)计算遍历中间栅格点m的所有相邻栅格点的f(n),若相邻节点n未加入open list,及时添加,并令g(n) = g(m) + d(m,n),此时n的父节点为m;若相邻点n在open list中,当g(n) > g(m) + d(m,n)时,更新g(n) = g(m) + d(m,n),设置n的父节点为m;
(5)遍历完所有备选点后重复步骤(2)-(4),直到目标点Q出现在open list中,搜索结束,按进入closed list的先后顺序输出中的节点栅格序列,该序列即为最优路径。

程序实现:

public class Astar extends Strategy {

  // open表
  Map<Integer, Node> open = new HashMap<>();
  // close表
  Map<Integer, Node> close = new HashMap<>();
  // 动作列表
  int[][] motion = get_motion();

  @Override
  public List<Location> nextstep(Location start, Location work) {
    Node startnode = new Node(start.getX(), start.getY(), 0, -1);
    Node worknode = new Node(work.getX(), work.getY(), 0, -1);

    // 将起点放入open表
    open.put(cal_grid_index(startnode), startnode);
    // 向起点的上下左右四个方向扩展
    while (true) {
      if (open.size() == 0) {
        System.out.println("open表空!查找失败!");
        break;
      }
      // 找到map中节点cost最小的元素
      int temp_cost = Integer.MAX_VALUE;
      int current_id = 0;
      Node current = null;
      for (Map.Entry<Integer, Node> entry : open.entrySet()) {
        int current_cost =
            entry.getValue().cost + cal_Manhattan_distance(entry.getValue(), worknode);
        if (current_cost < temp_cost) {
          current_id = entry.getKey();
          temp_cost = current_cost;
          current = entry.getValue();
        }
      }
      // 判断是否找到目标点
      if (current.x == work.getX() && current.y == work.getY()) {
        System.out.println("找到目标点!");
        worknode.pind = current.pind;
        worknode.cost = current.cost;
        break;
      }
      // 移除open表中的current节点
      for (Iterator<Integer> iterator = open.keySet().iterator(); iterator.hasNext();) {
        Integer key = iterator.next();
        if (key == current_id) {
          iterator.remove();
        }
      }

      // 将current节点放入close表中
      close.put(current_id, current);
      for (int i = 0; i < motion.length; i++) {
        Node node = new Node(current.x + motion[i][0], current.y + motion[i][1],
            current.cost + motion[i][2], current_id);
        int node_id = cal_grid_index(node);

        // 该节点不能被扩展则查找下一个节点
        if (!verify_node(node)) {
          continue;
        }
        // 该节点已被确认过则查找下一个节点
        if (close.containsKey(node_id)) {
          continue;
        }

        // 将新扩展节点放入open表中
        if (!open.containsKey(node_id)) {
          open.put(node_id, node);
        } else {
          // 更新节点的代价
          if (open.get(node_id).cost > node.cost) {
            open.put(node_id, node);
          }
        }
      }
    }
    open.clear();
    return cal_final_path(worknode, close);
  }

  /**
   * 根据终点回溯计算最短路径
   * 
   * @param worknode
   * @param close
   * @return
   */
  private List<Location> cal_final_path(Node worknode, Map<Integer, Node> close) {
    List<Location> ans = new ArrayList<>();
    int pind = worknode.pind;
    int i = 0;
    while (pind != -1) {
      // System.out.println(i);
      Node node = close.get(pind);
      Location location = new Location(node.x, node.y);
      ans.add(location);
      pind = node.pind;
    }
    return ans;
  }

  /**
   * 计算曼哈顿距离
   * 
   * @param now
   * @param end
   * @return
   */
  private int cal_Manhattan_distance(Node now, Node end) {
    return Math.abs(now.x - end.x) + Math.abs(now.y - end.y);
  }

  /**
   * 判断节点是否合理 1. 是否超出范围 2. 是否遇到货柜
   * 
   * @param node
   * @return
   */
  private boolean verify_node(Node node) {
    Container[][] map = ContainerMap.getMap();
    if (node.x < 0 || node.x >= ContainerMap.N || node.y < 0 || node.y >= ContainerMap.N) {
      return false;
    }
    if (map[node.x][node.y] != null) {
      return false;
    }
    return true;
  }

  /**
   * 计算栅格地图中的线性坐标
   * 
   * @param node
   * @return
   */
  private int cal_grid_index(Node node) {
    return node.y * ContainerMap.N + node.x;
  }

  /**
   * 定义动作及其代价 motion = {dx, dy, cost}
   * 
   * @return
   */
  private int[][] get_motion() {
    int[][] motion = {
  
  {1, 0, 1}, {0, 1, 1}, {-1, 0, 1}, {0, -1, 1}};
    return motion;
  }
}
版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

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

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


相关推荐

  • 做了6年的Java,java简历包装项目经验[通俗易懂]

    高频问题1.上一家公司,你为什么会离职?公司很好,但是公司调整了业务,接下来的发展路线和自己的目标不一致,所以要换工作工作太清闲,学不到知识,我不怕累,就是想多锻炼自己,想找具有挑战力工作公司的管理制度不也是很完善,没有晋升机会,我比较想进步,找一个更好的平台我想去优秀的公司,让自己变得更好2.为什么来我们这里?对原单位充满感恩,这是我人生中非常重要的经历,我认同原单位领导和文化两份工作的本质是一致的,行业和工作性质都有紧密联系选择一份新的工作不代表背叛过去,发扬原公司魅力,同时为了

    2022年4月11日
    340
  • VirtualBox安装Debian6的方法和步骤(详细)

    VirtualBox安装Debian6的方法和步骤(详细)下面是用VirtualBox安装Debian6的方法和步骤l新建一个文件夹,用于存放虚拟硬盘,如Debianl打开VirtualBox,点击新建l输入虚拟机名称,Debian_6l给虚拟机分配

    2022年7月2日
    25
  • 两地 三中心

    两地 三中心1、两地三中心同城双中心+异地灾备中心,“两地三中心”的灾备模式,方案兼具高可用性和灾难备份的能力。同城双中心是指在同城或邻近城市建立两个可独立承担关键系统运行的数据中心,双中心具备基本等同的业务处理能力并通过高速链路实时同步数据,日常情况下可同时分担业务及管理系统的运行,并可切换运行;灾难情况下可在基本不丢失数据的情况下进行灾备应急切换,保持业务连续运行。与异地灾备模式相比较,同城双中心具有投资成本低、建设速度快、运维管理相对简单、可靠性更高等优点。异地灾备中心是指在异地的城市建立一.

    2022年6月30日
    29
  • intellij idea 2021 激活码【2021.10最新】

    (intellij idea 2021 激活码)2021最新分享一个能用的的激活码出来,希望能帮到需要激活的朋友。目前这个是能用的,但是用的人多了之后也会失效,会不定时更新的,大家持续关注此网站~IntelliJ2021最新激活注册码,破解教程可免费永久激活,亲测有效,下面是详细链接哦~https://javaforall.net/100143.html…

    2022年3月27日
    28
  • pandas缺失值填充_python缺失值处理 fillna

    pandas缺失值填充_python缺失值处理 fillna约定:importpandasaspdimportnumpyasnpfromnumpyimportnanasNaN填充缺失数据fillna()是最主要的处理方式了。df1=pd.DataFrame([[1,2,3],[NaN,NaN,2],[NaN,NaN,NaN],[8,8,NaN]])df1代码结果:…

    2022年8月12日
    10
  • TOF相机-非扫描三维成像

    TOF相机-非扫描三维成像飞行时间法(TOF)3D成像,是通过给目标连续发送光脉冲,然后用传感器接收从物体返回的光,通过探测光脉冲的飞行(往返)时间来得到目标物距离。这种技术跟3D激光传感器原理基本类似,只不过3D激光传感器是逐点扫描,而TOF相机则是同时得到整幅图像的深度信息;TOF相机与普通机器视觉成像过程也有类似之处,都是由光源、光学部件、传感器、控制电路以及处理电路等几部

    2022年5月26日
    41

发表回复

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

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