数学建模常用的代码

Dijkstra算法找最短路径代码

算法的核心就是从原点出发(原点可以是自己定义的任意一个点),以原点为圆心,半径从小到大,判断原点到半径上面的点的最短距离,这个距离可能是圆心r0->r1(半径较小)->r2(半径较大)或者是r0->r2(如果存在r0到r2这条路径的话)

例子:

某公司在六个城市c1, c2,,,, c6 中有分公司,从 ci到 cj 的直接航程票价记在

下述矩阵的 (i, j) 位置上。(∞ 表示无直接航路),请帮助该公司设计一张城市 c1 到其它城市间的票价最便宜的路线图。

符号含义:用矩阵 a[n,n](n 为顶点个数)存放各边权的邻接矩阵, 行向量 pb 、 index1、 index2 、d 分别用来存放 P 标号信息、标号顶点顺序、标号顶点索引、最短通路的值。

其中分量 index2(i) 存放始点到第i 点最短通路中第i 顶点前一顶点的序号;

d(i) 存放由始点到第i 点最短通路的值。

求第一个城市到其它城市的最短路径的 Matlab 程序如下:

(可以直接复制下方代码运行)

其中a(1,2)表示第一个点到第二个点的距离,以此类推,在实际应用中先把所有点直接的距离矩阵写出来,不连通的点用无穷大表示

代码

Matlab 复制代码
clc,clear all
a=zeros(6);
a(1,2)=50;a(1,4)=40;a(1,5)=25;a(1,6)=10;               
a(2,3)=15;a(2,4)=20;a(2,6)=25;
a(3,4)=10;a(3,5)=20;
a(4,5)=10;a(4,6)=25;
a(5,6)=55;
a=a+a'                                                  
a(find(a==0))=inf %将a=0的数全部替换为无强大               
pb(1:length(a))=0;pb(1)=1;  %当一个点已经求出到原点的最短距离时,其下标i对应的pb(i)赋1
index1=1; %存放存入S集合的顺序
index2=ones(1,length(a)); %存放始点到第i点最短通路中第i顶点前一顶点的序号
d(1:length(a))=inf;d(1)=0;  %存放由始点到第i点最短通路的值
temp=1;  %temp表示c1,算c1到其它点的最短路。
while sum(pb)<length(a)  %看是否所有的点都标记为P标号
tb=find(pb==0); %找到标号为0的所有点,即找到还没有存入S的点
d(tb)=min(d(tb),d(temp)+a(temp,tb));%计算标号为0的点的最短路,或者是从原点直接到这个点,又或者是原点经过r1,间接到达这个点
tmpb=find(d(tb)==min(d(tb)));  %求d[tb]序列最小值的下标
temp=tb(tmpb(1));%可能有多条路径同时到达最小值,却其中一个,temp也从原点变为下一个点
pb(temp)=1;%找到最小路径的表对应的pb(i)=1
index1=[index1,temp];  %存放存入S集合的顺序
temp2=find(d(index1)==d(temp)-a(temp,index1));
index2(temp)=index1(temp2(1)); %记录标号索引
end
d, index1, index2
相关推荐
优美的赫蒂21 小时前
理解欧拉公式
线性代数·算法·数学建模
人大博士的交易之路1 天前
龙虎榜——20250422
大数据·数学建模·数据挖掘·缠论·缠中说禅·涨停回马枪·龙虎榜
烟锁池塘柳02 天前
【数学建模】孤立森林算法:异常检测的高效利器
算法·数学建模
烟锁池塘柳03 天前
【数学建模】随机森林算法详解:原理、优缺点及应用
算法·随机森林·数学建模
嵌入式冰箱3 天前
2025妈妈杯数学建模D题完整分析论文
数学建模
烟锁池塘柳03 天前
齐次坐标系下的变换矩阵
线性代数·数学建模·矩阵
CC数学建模3 天前
第十七届“华中杯”大学生数学建模挑战赛题目C题就业状态分析与预测完整 思路 代码 模型 结果 分享
数学建模
Lucifer三思而后行3 天前
零基础玩转AI数学建模:从理论到实战
人工智能·数学建模
小陈爱建模4 天前
【已更新】2025华中杯B题数学建模网络挑战赛思路代码文章教学:校园共享单车的调度与维护问题
数学建模
CC数学建模4 天前
2025第十七届“华中杯”大学生数学建模挑战赛题目B 题 校园共享单车的调度与维护问题完整成品正文33页(不含附录)文章思路 模型 代码 结果分享
数学建模