数学建模模型入门(二)
1.图论模型-Dijkstra算法
-
应用于求初始点到其他所有顶点的最短路径。
-
本质是一种标号法,给赋权图的每一个顶点记一个数,称为顶点的标号。(临时标号,称T标号。固定标号,称为P标号)T标号表示从始顶点到该标点的最短路长的上界,p标号则是从始顶点到该顶点的最短路长。
-
如图, 先为各个点命名,线段上的数字代表它们之间的距离,或者也可以理解成权重。按照图片我们需要的是v1到v11之间的最短路程。
-
对第一个点而言,到点2、3、4的距离分别为2、8、1。在图中用黄色的框里的数字标试出来路长。数字圈起来标为无穷大,是因为从第一个点没有直接到他们的路径。
-
容易看出第一个点所直接接触的到点四的路程最短,并把最短路程标出来。然后再可以得到到第二个点的距离最短。同理找到两点之间通过路长最短的,并把他们的路长总的标示出来。可以得到下面这个图,可以得到最短路长的路径。
-
通常这样一个一个看比较麻烦,而且一旦数据节点比较多就会更为繁琐易错。这样就引入了带权邻接矩阵,是表示顶点之间相邻关系的矩阵。
-
变换为矩阵形式为:
-
[0 2 8 1 Inf Inf Inf Inf Inf Inf Inf;
2 0 6 Inf 1 Inf Inf Inf Inf Inf Inf; 8 6 0 7 5 1 2 Inf Inf Inf Inf; 1 Inf 7 0 Inf Inf 9 Inf Inf Inf Inf; Inf 1 5 Inf 0 3 Inf 2 9 Inf Inf; Inf Inf 1 Inf 3 0 4 Inf 6 Inf Inf; Inf Inf 2 9 Inf 4 0 Inf 3 1 Inf; Inf Inf Inf Inf 2 Inf Inf 0 7 Inf 9; Inf Inf Inf Inf 9 6 3 7 0 1 2; Inf Inf Inf Inf Inf Inf 1 Inf 1 0 4; Inf Inf Inf Inf Inf Inf Inf 9 2 4 0;] -
再代入算法程序,即可快速得到最终结果。
2.图论模型-Floyd算法
- 可以看这样一个例子,带权邻接矩阵代码也如下。
function [D,path,min1,path1]=floyd(a,start,terminal)
D=a;n=size(D,1);path=zeros(n,n);
for i=1:n
for j=1:n
if D(i,j)~=inf
path(i,j)=j;
end,
end,
end
for k=1:n
for i=1:n
for j=1:n
if D(i,k)+D(k,j)<D(i,j)
D(i,j)=D(i,k)+D(k,j);
path(i,j)=path(i,k);
end,
end,
end,
end
if nargin==3
min1=D(start,terminal);
m(1)=start;
i=1;
path1=[ ];
while path(m(i),terminal)~=terminal
k=i+1;
m(k)=path(m(i),terminal);
i=i+1;
end
m(i+1)=terminal;
path1=m;
end
- 最后结果输出如图。第一遍是最廉价的航费表。由下面的路径矩阵可以得到,比如二到五的路径:先对应第二排第五列得到四,然后第四排第五列对应五,因此二到五的最廉价路径就是先从二到四再到五。去查上面对应的价格,第二行第五列对应30,第四行第五列对应10,一共花费40。
3.模拟退火模型
- 是模仿自然界推火现象得到的一种算法。利用概率突跳特征在解空间中随机寻找全局最优解。即,两两相比较时,给予一定概率跳出某个范围。
- 常于解决tsp问题。