广度优先搜索（BFS）及其matlab代码

爱学习的一一一

已于 2022-07-13 15:48:33 修改

阅读量3.9k

点赞数 7

分类专栏：路由算法介绍文章标签： matlab 宽度优先算法

于 2022-05-03 10:20:30 首次发布

本文链接：https://blog.csdn.net/weixin_42105848/article/details/124551006

版权

路由算法介绍专栏收录该内容

5 篇文章 0 订阅

订阅专栏

文章目录

前言
一、广度优先搜索是什么？
二、广度优先搜索的基本思路即其matlab代码
- 1、广度优先搜索的基本思路
- 2、matlab代码
总结

前言

广度优先搜索常被拿来与深度优先搜索进行比较，二者都是图形搜索算法但是寻找手段不同，接下来让我们一起进行学习吧~

一、广度优先搜索是什么？

广度优先搜索，全称为 $B r e a d t h f i r s t S e a r c h$ ，简称为 $B F S$ ，又称宽度优先搜索。它是从初始结点开始,应用产生式规则和控制策略生成第一层结点,同时检查目标结点是否在这些生成的结点中。若没有,再用产生式规则将所有第一层结点逐一拓展,得到第二层结点,并逐一检査第二层结点是否包含目标结点。若没有,再用产生式规则拓展第二层结点。如此依次拓展,检査下去,直至发现目标结点为止。如果拓展完所有结点,都没有发现目标结点,则问题无解。
以上介绍参考自以下网址：

https://blog.csdn.net/ZhuRanCheng/article/details/115588272

二、广度优先搜索的基本思路即其matlab代码

1、广度优先搜索的基本思路

广度优先搜索使用队列（queue）来实现，整个过程也可以看做一个倒立的树形：
（1）把根节点放到队列的末尾；
（2）每次从队列的头部取出一个元素，查看这个元素所有的下级元素，把它们放到队列末尾。并把这个元素记为它下一级元素的前驱；
（3）找到所要找的元素时结束程序；
（4）如果遍历整个树还没有找到，结束程序。
以上基本思路参考自百度百科：

https://baike.baidu.com/item/%E5%AE%BD%E5%BA%A6%E4%BC%98%E5%85%88%E6%90%9C%E7%B4%A2/5224802?fromtitle=%E5%B9%BF%E5%BA%A6%E4%BC%98%E5%85%88%E6%90%9C%E7%B4%A2&fromid=2148012&fr=aladdin

2、matlab代码

广度优先搜索：BFS.m

clear all;close all;clc
%初始化邻接压缩表
b=[1 2;1 3;1 4;2 4;
   2 5;3 6;4 6;4 7];

m=max(b(:));                %压缩表中最大值就是邻接矩阵的宽与高
A=compresstable2matrix(b);  %从邻接压缩表构造图的矩阵表示
netplot(A,1)                %形象表示

head=1;             %队列头
tail=1;             %队列尾，开始队列为空，tail==head
queue(head)=1;      %向头中加入图第一个节点
head=head+1;        %队列扩展

flag=1;             %标记某个节点是否访问过了
re=[];              %最终结果
while tail~=head    %判断队列是否为空
    i=queue(tail);  %取队尾节点
    for j=1:m
        if A(i,j)==1 && isempty(find(flag==j,1))    %如果节点相连并且没有访问过
            queue(head)=j;                          %新节点入列
            head=head+1;                            %扩展队列
            flag=[flag j];                          %对新节点进行标记
            re=[re;i j];                            %将边存入结果
        end
    end
    tail=tail+1;            
end

A=compresstable2matrix(re);
figure;
netplot(A,1)

构造图的矩阵表示：compresstable2matrix.m

function A=compresstable2matrix(b)
    [n ~]=size(b);
    m=max(b(:));
    A=zeros(m,m);

    for i=1:n
        A(b(i,1),b(i,2))=1;
        A(b(i,2),b(i,1))=1;
    end
end

画图函数：netplot.m（矩阵生成无向网络图）

%函数名netplot
%使用方法输入请help netplot
%无返回值%函数只能处理无向图
%作者：tiandsp
%最后修改：2012.12.26
function netplot(A,flag)
    %调用方法输入netplot(A,flag)，无返回值
    %A为邻接矩阵或关联矩阵
    %flag=1时处理邻接矩阵
    %flag=2时处理关联矩阵
    %函数只能处理无向图
    if flag==1      %邻接矩阵表示无向图
        ND_netplot(A);
        return;
    end
    
    if flag==2      %关联矩阵表示无向图
        [m n]=size(A);      %关联矩阵变邻接矩阵
        W=zeros(m,m);
        for i=1:n
            a=find(A(:,i)~=0);
            W(a(1),a(2))=1;
            W(a(2),a(1))=1;
        end
        ND_netplot(W);
        return;
    end
           
    function ND_netplot(A)
        [n n]=size(A);
        w=floor(sqrt(n));       
        h=floor(n/w);        
        x=[];
        y=[];
        for i=1:h           %使产生的随机点有其范围，使显示分布的更广
            for j=1:w
                x=[x 10*rand(1)+(j-1)*10];
                y=[y 10*rand(1)+(i-1)*10];
            end
        end
        ed=n-h*w;
        for i=1:ed
           x=[x 10*rand(1)+(i-1)*10]; 
           y=[y 10*rand(1)+h*10];
        end
        plot(x,y,'r*');    
        
        title('网络拓扑图'); 
        for i=1:n
            for j=i:n
                if A(i,j)~=0
                    c=num2str(A(i,j));                      %将A中的权值转化为字符型              
                    text((x(i)+x(j))/2,(y(i)+y(j))/2,c,'Fontsize',10);  %显示边的权值
                    line([x(i) x(j)],[y(i) y(j)]);      %连线
                end
                text(x(i),y(i),num2str(i),'Fontsize',14,'color','r');   %显示点的序号            
                hold on;
            end
        end  
    end
    
end