【数模】用matlab程序求出有向图的所有有向圈（有向环）并可视化

Dream of Grass

于 2022-06-13 00:17:19 发布

阅读量2.1k

点赞数 7

分类专栏：数学建模文章标签： matlab 图像处理矩阵图论

本文链接：https://blog.csdn.net/dream_of_grass/article/details/125252414

版权

数学建模专栏收录该内容

12 篇文章 13 订阅

订阅专栏

这篇博客展示了如何使用MATLAB程序来绘制有向图，并利用allcycles函数找出图中的所有循环。代码首先创建了一个有向图，然后通过plot函数显示，接着调用allcycles函数获取所有循环，并使用highlight高亮显示了前五个循环。这是一个简单的MATLAB图形操作实例，适合初学者在15分钟内完成。

摘要由CSDN通过智能技术生成

用matlab程序求出下图的所有有向圈。

在这里插入图片描述

首先在matlab中把图画出来

% 首先在matlab中把图画出来
s = [1 2 2 3 3 3 4 4 5];
t = [2 5 4 3 1 2 3 5 2];
G = digraph(s,t)
p = plot(G);

在这里插入图片描述

然后再通过allcycles函数找到圈

[cycles,edgecycles] = allcycles(G);

最后可视化出来

p = plot(G);
highlight(p,'Edges',edgecycles{1},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{2},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{3},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{4},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{5},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)

第一个圈

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-Mzs0tROL-1655050403540)(D:/typora%E5%9B%BE%E7%89%87/image-20220613001024751.png)]

第二个圈

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-zQj0QatG-1655050403541)(D:/typora%E5%9B%BE%E7%89%87/image-20220613001111536.png)]

第三个圈

在这里插入图片描述

第四个圈

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-8bM32QfZ-1655050403543)(D:/typora%E5%9B%BE%E7%89%87/image-20220613001136209.png)]

第五个圈

[外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传(img-aaCO1Tei-1655050403544)(D:/typora%E5%9B%BE%E7%89%87/image-20220613001151485.png)]

完整代码如下

% 首先在matlab中把图画出来
s = [1 2 2 3 3 3 4 4 5];
t = [2 5 4 3 1 2 3 5 2];
G = digraph(s,t)
p = plot(G);
[cycles,edgecycles] = allcycles(G);
p = plot(G);
highlight(p,'Edges',edgecycles{1},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{2},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{3},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{4},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)
p = plot(G);
highlight(p,'Edges',edgecycles{5},'EdgeColor','r','LineWidth',1.5,'NodeColor','r','MarkerSize',6)

谢谢大家的观看，这题不难，15min搞定。