第六十七章 ：两秒！

第六十七章 ：两秒！ (第2/2页)

“你原来用协调参数处理的时候，每一条经过西直门的线路都要和其他所有经过西直门的线路做配对协调。”
　　
　　说着，他提起笔，在稿纸图层边重写了一个算式。
　　
　　“假设西直门有n条线路经过，协调参数的数量是O(n²)。”
　　
　　“那么在共享层里，这n条线路各自向西直门提交一个需求值，共享层只需要处理一个变量。”
　　
　　“也就是说，西直门的枢纽容量分配，那么复杂度从O(n²)降到了O(n)。”
　　
　　“更重要的是，这个O(n)的过程是有解析解的，即【L(x，λ)=i=1∑n(x−di)2+λ(x−C)】”
　　
　　“其中d_i是第i条线路提交的需求值，x是共享层分配给该枢纽的总流量，C是枢纽容量上限，λ是拉格朗日乘子。对x求导，令导数为零.....”
　　
　　许志远看着那个矩阵表达式，沉默了几秒。
　　
　　然后他把韩川的稿纸拉过来，对照着上面的公式，开始修改自己的MATLAB代码。
　　
　　屏幕上的代码一行行流畅地刷过，不到半个小时左右的时间，许志远把代码粘进主程序里，然后打开题目附件里的站点数据文件，配好输入参数，按下了运行键。
　　
　　屏幕上没有立刻跳出结果，MATLAB的命令窗口里，光标一闪一闪地，正在跑迭代。
　　
　　等了大概两秒的时间，结果出来了！
　　
　　两秒！
　　
　　看着屏幕上的结果，许志远下意识地以为是程序出buG了。
　　
　　这速度，是不是也太快了？
　　
　　他滚动鼠标，仔细地检查了一下迭代日志，从头到尾看了一遍。
　　
　　第一轮迭代，共享层收到来自三层共一百二十条线路的需求数据，汇总成三十七个枢纽的容量分配方案。
　　
　　第二轮迭代，需求变化量比第一轮下降了百分之七十。
　　
　　第三轮，再降百分之五十。
　　
　　到第五轮的时候，各层的线路方案和共享层的枢纽分配方案已经基本一致，偏差量降到了初始值的千分之一以下。
　　
　　“嘶！”
　　
　　“这怎么可能！？”
　　
　　“五轮收敛，每一轮的计算时间不到一秒！”
　　
　　“这不是超算啊？！”
　　
　　看着屏幕上的结果，许志远的喉结滚动了一下。
　　
　　他参加过这届建模大赛，这道题拿国一的队伍有十七支，而其中从建模到计算出结果，用时最短的也花费了足足六个小时的时间。
　　
　　即便是单纯看模型计算的速度，最短也耗时半个多小时。
　　
　　而现在，他用自己的笔记本电脑做测算，仅仅是两秒的时间，就完成了线路规划。
　　
　　虽然说测试用的数据量远比不上正式建模大赛的数据库，但别忘了，他跑模型用的设备，也远比不上建模的计算机啊。
　　
　　两秒钟！
　　
　　这速度，如果放到07年的建模大赛上，用爆杀全场来形容都太保守了。
　　
　　真要说，这种级别的产品，理论上已经不再是单纯的竞赛建模了，它具备了商业化的可能性！