新手进阶建模（7）图示法图论

2336426014

1 o2 t( F  r& d; ^% X2 ^
       关于图示法，百度给出的定义是：图示法是用曲线或图形表示数据之间的关系，从图形中能直观地反映出数据变化的趋势，如递增性或递减性，是否具有周期性变化规律等。在图上作进一步处理可以获得更多信息，如 最大值、最小值，做出切线，求出曲线下包围的面积等。但是图形的缺点为不能进行数学分析。工程测试中，多采用直角坐标系绘制测量数据的图形，也可采用对数坐标系、极坐标系等坐标系来描述。在直角坐标系中描绘曲线时，应该使曲线通过尽可能多的数据，曲线以外的数据则尽可能靠近曲线，并且曲线两侧数据点数目要大致相等，最后得到一条平滑曲线。
3 E+ n4 R* O. |7 o/ A6 T; i      我自己的定义是：图示法就是用方块加箭头来表示元素之间的关系（具体啥关系在箭头上加文字表达就好）。
: T& Y3 _% a+ ~4 H% r* D7 r      建模中用图的好处有很多，我自己经验感觉的话，主要就是方便评卷人阅读，能一目了然我们思路（前提是图作的好看和整齐），不会因为论文看着没意思毙掉。另一个就是方便后面的论文排版。可以提前准备好论文各个部分内容的版式。
3 F- ~0 Z9 p; j   
1 f4 k( {0 q. V8 j. W     图论与图示法我感觉有点同根的意思，基本的思想都是表示两个事物的特定的联系，只不过图论后来发展成了一门单独的学科。
      建模中遇到指派问题（通俗讲就是其群人如何从一堆鞋子中找到适合自己的鞋子）时候，图论就会被排上用场，需要建模者对矩阵运算和集合知识有一定的基础（会matlab运算矩阵也可以）。其优点是通过矩阵的变换，找到我们想要的最佳指派方案（找鞋子步骤）或者步骤。相比于编程序让计算机挨个试，这种办法计算更快。
     举个简单例子：
       某公司在六个城市（c1，c2...c6）中有都分公司，从c(i) 到 c(j) 的直接航程票价记在下面的矩阵 ,（i排j列表示从ci到cj的价格），请帮助该公司设计一张城市 1 c 到其它任意城市的最便宜路线。
8 Q: I4 d+ O" D! J3 G  q( G- k（矩阵发现进不来这个位置，所以放附件图片了）
用上面矩阵存放各边权，行向量 pb、 1 index 、 2 index 、 d 分别用来存放P 标号信息、标号顶点顺序、标号顶点索引、最短通路的值。
5 C, s: k3 e9 Q      pb(i)=0表示该点未标号，pb(i)=1表示已经标号
: x7 C4 d8 T. ^& ^    index(i) 存放始点到第i 点最短通路中第i 顶点前一顶点的序号；  d(i) 存放由始点到第i 点最短通路的值。
9 S1 `9 t! E3 `# F9 N; Y) }7 T# A求解程序如下：
      
# [, S+ B/ r! m. A+ a; @3 tclc,clear a=zeros(6);
( I' j% Z) q- `( g5 r2 d9 sa(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;
9 m1 m# c- |! h/ K6 Ga=a+a';
a(find(a==0))=inf;
pb(1:length(a))=0;pb(1)=1;index1=1;index2=ones(1,length(a));
* j6 ~% ?7 S- n7 d7 r7 w! wd(1:length(a))=inf;d(1)=0;temp=1;
while sum(pb)<length(a)   
         tb=find(pb==0);   
# Q3 ]9 n7 j& m         d(tb)=min(d(tb),d(temp)+a(temp,tb));   
; D5 p& N* [4 D3 Y5 X          tmpb=find(d(tb)==min(d(tb)));   
: D' c9 z0 B, ^% C$ X         temp=tb(tmpb(1));   
% i6 ?$ T! X6 P2 H  ?7 }" n          pb(temp)=1;   
) ~6 R8 E+ a8 d          index1=[index1,temp];    temp2=find(d(index1)==d(temp)-a(temp,index1));    index2(temp)=index1(temp2(1));
end
! j& W- ^) A. A6 @# s* x( A6 `1 j3 Q

: F0 R5 T  I% H* D& w& L更多图论内容比如迪克斯屈拉算法，Floyd算法之类的。见附件
  l2 [# s. g0 L4 c   

1 }8 g6 p  g7 B5 H; x8 E+ z2 U

李江杰

好好好88888