1、机场选址问题摘 要针对机场选址问题,文章共建立了三个模型用以解决该类问题。为了计算出任意两城市之间的距离,我们利用公式(1)将利用题目中所给的大地坐标得出了任意两点之间的距离,见附录2。对于问题1,我们主要利用0-1变量法,从而对问题进行了简化。我们设了第i个城市是否建支线机场的以及第i个城市是否是以第j个支线机场为最近机场的。然后将任意两点之间的距离与该城市的总人数之积,再乘以0-1变量,最后得出每一个所有城市到最近机场的距离与该城市人口的乘积,然后利用LINGO进行编写程序,进行最优化求解,最后得出的结果见表1和表2,各大城市以及支线机场的分布见图2。 对于问题2, 该问题是属于多目标规划
2、的问题,目标一是居民距离最近机场的距离最短,目标二是每个机场覆盖人口数尽可能相等。我们在第一题的基础上,又假设了一些正、负偏差变量,对多个目标函数设立优先级,把目标函数转化为约束条件,进而求得满足题目要求的结果。对于问题3, 我们分析到影响客流量的因素是GDP跟居民人数,所以通过所搜集的资料分析我们给予这两个因素以不同的权重。然后同样采取问题2中所给的反求机场覆盖的方法,求的各个机场所覆盖的客流量,再让其在平均客流量水平上下浮动。通过LINGO程序的运行得到的六个机场的坐标见表6,六个机场的分布见图7。 针对论文的实际情况,对论文的优缺点做了评价,文章最后还给出了其他的改进方向,以用于指导实际
3、应用。关键词:选址问题;多目标规划;LINGO;0-1变量法;加权1问题的重述近年来,随着我国经济社会的迅猛发展,公共交通基础设施日趋需要进一步完善与提高。支线机场作为我国交通运输体系的有机组成部分,对促进欠发达地区经济社会的发展具有基础性的作用。现某区域有30个城市,本区域计划在未来的五年里拟建6个支线机场。任务1,确定6个支线机场的所在城市,建立居民到最近机场之间的平均距离最小的数学模型。任务2,在任务一基础上,确定6个支线机场的所在城市,建立使得每个支线机场所覆盖的居民人数尽可能均衡的数学模型。任务3,在任务一基础上,根据近一年每个城市的GDP情况,确定6个支线机场的所在城市,建立使得每
4、个支线机场的客流量尽量均衡的数学模型。2问题的分析21问题1题目要求是建立居民到最近机场之间的平均距离最小的数学模型,该问题其实就是利用的0-1变量建立的模型。首先我们设两个0-1变量,一个是控制某个城市是否为支线机场的,一个是控制某个城市的最近机场是哪一个的。针对于上述两个0-1变量,我们分别设立了约束条件。同时又为了满足问题所要求的使局面平均距离最小,我们将某一个城市到离它最近的机场的距离与该城市的人口乘积作为目标函数,在LINGO软件中,通过设立一约束条件,最后将目标函数进行最优化求解。22问题2该问题可以归结为多元目标线性规划的问题,所以我们在第一问的基础上又增加了一个目标函数,最后利
5、用加权的方法将两个目标函数转化成了一个目标函数,将另一个目标函数作为约束条件。同时我们又引入了正负偏差变量,通过控制该变量达到覆盖居民人数均衡以及居民到城市之间的平均距离尽量小。23问题3该问题要求的是客流量尽量均衡,经过分析可以知道,城市的GDP越高,说明该城市经济越繁荣,货币流通越快,从而反映出客流量越大。另一方面城市越大、人口越多,也在一定程度上反映出了该城市客流量越大。基于上述两点,我们对GDP跟城市人口分别给予了不同的权重来反映其对客流量的影响大小。按照第二问的方法,我们依然利用多元目标线性规划的只是进行求解。通过LINGO编写程序,最中求得可行解。3模型的假设与符号说明3模型的假设
6、(1)各个城市的人口在某个较长的时间段内是不进行流动的,基本保持不变。(2)两城市之间的距离都按照直线来计算,不存在弯曲线段的情况。(3)各个城市都满足支线机场的建设条件,不存在某个城市不能建设支线机场的情况。(4)假设各个支线机场是建在各个城市里面的,也就是说,我们将在这30个城市里面选取6个城市建设支线机场。(5)各城市在未来一段时间内发展水平基本不变。3符号说明 符号符号说明若第j个城市建设支线机场,则第i个城市到第j个城市距离最近时,其余的各种情况,。第i个城市到第j个城市之间的距离,其中,。第i个城市的总人口,其中。第i个城市被选为建设支线机场时,;否则, 。按优先顺序k极小化的现实
7、目标或约束偏差变量的线性函数,其中。、正偏差变量,其中。、负偏差变量,其中。第i个城市的GDP,。问题一中求得的目标函数的最小值。30个城市总的GDP分配到六个支线机场的平均GDP。4模型的准备首先我们将30座城市在大地坐标(经纬度)下的位置用MATLAB软件画出以下图形(源程序见附录1):图1 大地坐标系下的各大城市的位置题目条件所给的数据是经纬度,显然是不能进行距离计算的,首先我们从网上查取了一个公式,用于计算地球上任意两点之间的距离。所以,我们就利用该公式计算出了任意两个城市之间的距离。公式如下: (1) (2)其中代表的是地球上两点的经度,、代表的是地球上两点的纬度,R代表的是地球的平
8、均半径,。最后求的任意两点之间的距离见附录2。5模型的建立与求解51问题1的模型建立与求解我们将每个城市到离该城市最近的那个支线机场的距离与该城市的总人数之积做为目标函数,求当该目标函数最小时,支线机场所建立的城市。所以该模型的目标函数可以写作如下: (3)定义作为0-1变量,所以应该满足: (4) (5)为了使建设机场的个数为6个,还需要满足以下条件: (6)为了是一座城市只能到一个机场的距离最短,还需要满足以下条件: (7)为了使满足地i个城市是否是以第j个城市中的机场的距离作为最短距离时,所以,还应满足以下条件: (8)所以可得出该模型的目标函数为:约束条件: 我们利用LINGO软件进行
9、编程(源程序见附录3,运算结果见附录4),得到如下结果: 表1 建设支线机场的城市城市编号城市说明1、5、7、11、20、23在这些城市里面建立支线机场可以使居民到支线机场的平均距离最短。也就是说,只有在以上表中的几个城市建设机场才能使得居民离最近支线机场的距离最小,即应把机场建设在城市1、5、7、11、20、23六处。以下是以某个支线机场作为最近的机场的城市表2 以某支线机场作为最近的机场的城市编号城市编号到该支线机场最近的城市编号11、10、1754、5、6、8、9、1972、3、7、121111、13、14、15、162018、20、21、22、24、25、26、27、28、29、302
10、323现在我们将利用MATLAB将各个城市以及支线机场所在的位置以图像的形式表现出来(源程序见附录5),图形如下:图2 各大城市以及支线机场的位置52问题2的模型建立与求解我们将各个正负偏差变量的之和作为目标函数,将求解该目标函数的最小值作为问题的目的,所以目标函数可以写成下面的形式: (9)定义作为变量,所以应该满足: (10) (11)为了使建设机场的个数为6个,还需要满足以下条件: (12)为了是一座城市只能到一个机场的距离最短,还需要满足以下条件: (13)为了使满足地i个城市是否是以第j个城市中的机场的距离作为最短距离时,所以,还应满足以下条件: (14)为了使覆盖的居民人数达到平衡
11、,所以有: (15)在上式中,1655是30个城市的总人口平均分配到6个支线机场的人数。为了使覆盖居民到支线机场的平均距离较小,还需要满足以下情况: (16)所以该模型的目标函数为:约束条件为: 我们利用LINGO软件进行编程(源程序见附录6,运算结果见附录7),得到如下结果:表3 建设支线机场的城市城市编号城市说明3、8、9、13、15、28在这些城市里面建立支线机场可以使居民到支线机场的平均距离最短。 以下是以某个支线机场作为最近的机场的城市表4 以某支线机场作为最近的机场的城市编号城市编号到该支线机场最近的城市编号33、7、12、2086、8、14、17、18、2192、4、5、9、10
12、、26131、13、29、301511、15、16、24、252819、22、23、27、28现在我们将利用MATLAB将各个城市以及支线机场所在的位置以图像的形式表现出来(源程序见附录8),图形如所示下:图3 各大城市以及支线机场的位置53问题3的模型建立与求解对于该问题,我们同样利用的是多元线性规划的问题。写出多了目标函数,将这些目标函数转换成为但目标函数,最后利用LINGO进行求解。其中我们的目标函数是将各个正负偏移量取得最小值时,选取的某些城市作为支线机场的建设城市,从而作为最优解。所以该问题的目标函数如下: (17)定义作为0-1变量,所以应该满足: (18) (19)为了使建设机场
13、的个数为6个,还需要满足以下条件: (20)为了使一座城市只能到一个支线机场的距离最短,还需要满足以下条件: (21)为了使满足地i个城市是否是以第j个城市中的机场的距离作为最短距离时,所以,还应满足以下条件: (22)为了使覆盖的居民人数达到平衡,所以有: (23)在上式中,1655是30个城市的总人口平均分配到6个支线机场的人数。为了使覆盖居民到支线机场的平均距离较小,还需要满足以下情况: (24)为了使各个机场的客运量达到均衡,我们不禁用其覆盖的军民人数来表示,而且还用该支线机场所覆盖的城市的GDP之和均衡来表示各个支线机场的客流量均衡,所以约束条件见下式: (25)所以该模型的目标函数
14、见下式:约束条件为:我们利用LINGO软件进行编程(源程序见附录9,运算结果见附录10),得到如下结果:表5 建设支线机场的城市城市编号城市说明3、5、8、9、13、15在这些城市里面建立支线机场可以使居民到支线机场的平均距离最短。以下是以某个支线机场作为最近的机场的城市表6 以某支线机场作为最近的机场的城市编号城市编号到该支线机场最近的城市编号31、3、2052、5、7、2386、8、17、22、24、25、2894、9、12、14、181310、11、13、191515、16、21、26、27、29、30现在我们将利用MATLAB将各个城市以及支线机场所在的位置以图像的形式表现出来(源程序
15、见附录11),图形如下:图4 各大城市以及支线机场的位置6模型的推广与改进方向61模型的改进我们在建立模型的时候假设任意两城市之间的距离是直线的,也就是按照最短距离来算的,可是实际上并不是这样的。如果能知道两城市之间的真实距离,这样算出来的结果应该会更加准确些。62模型的推广本模型不仅适用于此题目中飞机场的选址,而且可应用于实际生活中多种优化问题,能够很好地解决实际问题,给出优化方案起到了减少资源能源的投入,提高了收益的作用。比如移动、联通信号塔选址覆盖问题(尽可能少的信号塔覆盖全部地区并满足客户需求),网通铁通线路布局(减少了电缆线路),银行选址,派出所驻地(及时的处理问题、服务人民),物流
16、公司配货站点设置(快速高效的配送货物),相邻几个村庄建设学校选址等问题。7模型的优缺点71模型的优点(1)该模型比较简单,运用了0-1变量和多目标线性规划,整体来说就是属于一个最优化的问题。我们建立的模型,能够运用较简单的程序就能将结果解答出来,简捷、方便、易懂。(2)该模型适应性广,可适用于许多问题的选址问题。72模型的缺点(1)我们假设各个城市都能建设支线机场,忽略了不能建设机场的情况。(2)模型所给的最优解只单纯性的给出了所在城市,而城市所在地区幅员辽阔,所以答案略显得笼统。参考文献1 J.A.邦迪图论及其应用M.北京:科学出版社,19842 韩中庚长江水质综合评价与预测的数学模型J工程
17、数学学报,20053 张宏伟kng08.0及其在环境系统优化中的应用M天津:天津大学出版社,20054 袁新生LINGO和EXCEL在数学建模中的作用M北京:科学出版社,20075 韩中庚数学建模方法及其应用(第二版)M北京:高等教育出版社,2009附 录附录1:画各城市在经纬度坐标系下的位置图程序(MATLAB) x=data3(1:30,2); y=data3(1:30,3); cftool hold on legend(城市)附录2:任意两座城市之间的距离表表7 任意两座城市之间的距离10146.2087527.2021319.5801330.579265.5566290.1473214
18、6.20870452.8953399.5192451.0169406.705144.03573527.2021452.89530461.9687604.6329662.4375425.90954319.5801399.5192461.96879.49E-05147.0755245.1018512.16365330.579451.0169604.6329147.07559.49E-05132.061583.98076265.5566406.705662.4375245.1018132.0610549.3757290.1473144.0357425.9095512.1636583.9807549.
19、37508322.1059434.9494570.4819112.042735.03323153.9619564.36429382.0044467.9079504.59668.94566140.9123263.8845580.995210106.027121.2429421.2866278.3842336.6706312.6631247.318211230.7495313.8892750.9839512.1978464.44345.7461426.168412301.3644222.7853638.7191609.0307631.6312552.4608218.431613431.511654
20、3.0116957.3559630.9566531.3655399.8005662.464714213.5075332.4318739.7067448.4028383.7803259.5661462.465315338.0798424.812861.6865596.7951527.6315399.0967531.803216310.9952367.8226816.9819603.565556.3732435.7018458.735117150.8514144.7333597.2657469.1772476.8148393.8177242.956418450.5393549.9555977.67
21、97672.9524580.5681448.5217660.778319399.0649544.2358858.9645443.0354309.6535201.3211686.925320506.7538560.60151012.185786.1422719.9278590.8671636.542321556.7958609.94251061.854833.963764.5274634.5581683.135522528.5275589.24891039.45799.467727.0272596.5862669.450723870.0234918.33651371.2051139.979105
22、8.702926.775978.436324621.0185701.22141144.005860.1246767.7907635.8289793.425625492.6097551.96321002.3767.4209698.7519569.2037632.703626646.4403667.06731116.981946.5458893.4887767.3262709.67427543.2451597.75371049.385819.8059750.4162620.5047672.571928599.0317634.98591087.693889.6416828.9302700.80696
23、92.575729428.5377479.8684931.4416715.3912657.6154531.4196558.218630658.5439703.33041156.204939.1613869.2498738.9001765.96838910111213141322.1059382.0044106.027230.7495301.3644431.5116213.50752434.9494467.9079121.2429313.8892222.7853543.0116332.43183570.4819504.596421.2866750.9839638.7191957.3559739.
24、70674112.042768.94566278.3842512.1978609.0307630.9566448.4028535.03323140.9123336.6706464.44631.6312531.3655383.78036153.9619263.8845312.6631345.7461552.4608399.8005259.56617564.3642580.9952247.3182426.1684218.4316662.4647462.465380110.9692317.7541472.2642623.337553.2906395.33149110.96920346.6834562
25、.5093675.3566661.4524492.418310317.7541346.68340334.4388335.5164536.6083318.674511472.2642562.5093334.43889.49E-05295.6673236.965891.9305412623.337675.3566335.5164295.66730521.8714367.632313553.2906661.4524536.6083236.9658521.87140.000134218.017514395.3314492.4183318.674591.93054367.6323218.01750.00
26、013415541.4334640.4299443.2496111.4915377.7094145.6242148.433316564.6577654.7271409.90292.47302288.1646239.6622176.989717471.0135532.5704215.6171183.2393158.6585419.8684229.527218600.7262707.1916556.5516236.5741503.911360.0657240.821519340.4724450.4702478.0227351.9951628.9406273.6213264.844620733.71
27、7831.7732606.8222278.8287437.8628250.7407339.385321779.2573878.5509656.9706328.4046480.9564277.647386.243522742.4258842.7858630.3214298.4273473.4299237.4302351.1342231076.9371180.451970.0621640.9429765.2642534.5851691.938124788.4059894.9956726.2713393.0356605.805237.7907420.562925713.1048812.0709593
28、.9085263.1121438.4273226.4582319.720226904.3961997.1976737.1299434.7804492.6558435.1421509.716427765.0983864.3661643.6835314.5747471.9155265.7957372.061128841.6099937.5706694.6222378.0484480.7581357.418446.328429669.0514763.7197527.3669203.9999364.9369237.9016273.928730884.3325983.8589757.0125432.00
29、84555.6791370.2764491.5565151617181920211338.0798310.9952150.8514450.5393399.0649506.7538556.79582424.812367.8226144.7333549.9555544.2358560.6015609.94253861.6865816.9819597.2657977.6797858.96451012.1851061.8544596.7951603.565469.1772672.9524443.0354786.1422833.9635527.6315556.3732476.8148580.568130
30、9.6535719.9278764.52746399.0967435.7018393.8177448.5217201.3211590.8671634.55817531.8032458.7351242.9564660.7783686.9253636.5423683.13558541.4334564.6577471.0135600.7262340.4724733.717779.25739640.4299654.7271532.5704707.1916450.4702831.7732878.550910443.2496409.902215.6171556.5516478.0227606.822265
31、6.970611111.491592.47302183.2393236.5741351.9951278.8287328.404612377.7094288.1646158.6585503.9113628.9406437.8628480.956413145.6242239.6622419.868460.0657273.6213250.7407277.64714148.4333176.9897229.5272240.8215264.8446339.3853386.243515094.4048289.7937129.3947344.8858192.3545238.16111694.40480.000
32、134224.5387215.9593418.9134197.0116247.152717289.7937224.53870418.2559480.026415.9109465.212218129.3947215.9593418.25590332.8155192.5851217.720719344.8858418.9134480.026332.81550513.4402547.013120192.3545197.0116415.9109192.5851513.4402050.1489621238.1611247.1527465.2122217.7207547.013150.1489602220
33、2.7009222.4977444.9029177.5981506.607837.686240.4246423544.0877560.1601774.0008480.0777803.7818363.3654313.259424283.022335.8087559.746187.8067504.6906175.7938150.779325171.7777185.431407.606168.5347489.266524.3823566.5523426369.5956342.9849527.1691375.7903699.5421186.171160.113427223.978234.0087453
34、.0499206.0371534.238937.2017114.1848328301.7102289.1376491.576297.7916624.1345111.681381.3153529135.8495117.6339335.2004189.2728477.321180.76234130.502330343.4481347.6351559.029310.5419643.2302153.4296105.3146222324252627281528.5275870.0234621.0185492.6097646.4403543.2451599.03172589.2489918.3365701
35、.2214551.9632667.0673597.7537634.985931039.451371.2051144.0051002.31116.9811049.3851087.6934799.4671139.979860.1246767.4209946.5458819.8059889.64165727.02721058.702767.7907698.7519893.4887750.4162828.93026596.5862926.775635.8289569.2037767.3262620.5047700.80697669.4507978.4363793.4256632.7036709.674
36、672.5719692.57578742.42581076.937788.4059713.1048904.3961765.0983841.60999842.78581180.451894.9956812.0709997.1976864.3661937.570610630.3214970.0621726.2713593.9085737.1299643.6835694.622211298.4273640.9429393.0356263.1121434.7804314.5747378.048412473.4299765.2642605.805438.4273492.6558471.9155480.7
37、58113237.4302534.5851237.7907226.4582435.1421265.7957357.41814351.1342691.9381420.5629319.7202509.7164372.0611446.328415202.7009544.0877283.022171.7777369.5956223.978301.710216222.4977560.1601335.8087185.431342.9849234.0087289.137617444.9029774.0008559.746407.606527.1691453.0499491.57618177.5981480.
38、0777187.8067168.5347375.7903206.0371297.791619506.6078803.7818504.6906489.2665699.5421534.2389624.13452037.6862363.3654175.793824.38235186.17137.20171111.68132140.42464313.2594150.779366.55234160.113414.1848381.31535220342.9691138.107937.29709198.882728.54578120.399123342.96910299.106377.8321281.1582326.7793285.992824138.1079299.1069.49E-05168.9968286.3377151.7908215.23392537.29709377.8321168.99680.000134210.45452.38587135.312126198.8827281.1582286.3377210.4540170.350178.800472728.54578326.7793151.790852.38587170.3501091.85902