数学建模课程设计优化问题.docx
《数学建模课程设计优化问题.docx》由会员分享，可在线阅读，更多相关《数学建模课程设计优化问题.docx（18页珍藏版）》请在冰豆网上搜索。
数学建模课程设计优化问题
摘要
在手机普遍流行的今天，建设基站的问题分析对于运营商来说很有必要。
本文针对现有的条件和题目的要求进行讨论。
在建设此模型中，核心运用到了0-1整数规划模型，且运用lingo软件求解。
对于问题一：
我们引入0-1变量，建立目标函数：
覆盖人口最大数=所有被覆盖的社区人口之和，即max=
，根据题目要求建立约束条件，并用数学软件LINGO对其模型求解，得到最优解。
对于问题二：
同样运用0-1整数规划模型，建立目标函数时，此处假设每个用户的正常资费相同，所以68%可以用减少人口来求最优值，故问题二的目标函数为：
max=
上述模型得到最优解结果如下：
研究问题
建中继站位置
所需费用
最优值
问题一
2、4、6、7
45百万
覆盖中人口数109.5千人
问题二
2、4、6、7
45百万
获得资费83.74a
关键字：
基站;0-1整数规划；lingo软件
1问题的重述．．．．．．．．．．．．．．．．．．．．．．．．．3
2问题的分析．．．．．．．．．．．．．．．．．．．．．．．．．4
3模型的假设与符号的说明．．．．．．．．．．．．．．．．．．．5
3.1模型的假设．．．．．．．．．．．．．．．．．．．．．．5
3.2符号的说明．．．．．．．．．．．．．．．．．．．．．．5
4模型的建立及求解．．．．．．．．．．．．．．．．．．．．．.5
4.1模型的建立．．．．．．．．．．．．．．．．．．．．．．5
4.2模型的求解．．．．．．．．．．．．．．．．．．．．．．6
5模型结果的分析．．．．．．．．．．．．．．．．．．．．．．．7
6优化方向．．．．．．．．．．．．．．．．．．．．．．．．．．7
7参考文献．．．．．．．．．．．．．．．．．．．．．．．．．．8
8、附录．．．．．．．．．．．．．．．．．．．．．．．．．．．9
1、问题的重述
某手机运营商准备在一个目前尚未覆盖的区域开展业务，计划投资5000万元来建设基站。
该区域由15个社区组成，有7个位置可以建设基站，每个基站只能覆盖有限个社区。
图1是该区域的示意图，每个社区简化为一个多边形，每个可以建设基站的位置已用黑点标出。
由于地理位置等各种条件的不同，每个位置建设基站的费用也不同，且覆盖范围也不同。
表1中列出了每个位置建设基站的费用以及能够覆盖的社区，表2列出了每个社区的人口数。
表1每个位置建设基站的费用及所能覆盖的社区
位置
1
2
3
4
5
6
7
费用（百万元）
9.5
7
19
14
17.5
13
11
覆盖社区
1,2,4
2,3,5
4,7,8,
10
5,6,8,9
8,9,12
7,10,11,12,15
12,13,
14,15
表2每个社区的人口数量
社区
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
人口（千人）
2
4
13
6
9
4
7．5
12.5
10
11
6
14
9
3.5
6
问题一：
在不超过5000万建设费用的情况下，在何处建设基站，能够覆盖尽可能多的人口；
问题二：
考虑到基站出现故障维修的时候可能会出现所覆盖的社区信号中断等问题，为此对通讯资费进行了调整，规定，仅有一个基站信号覆盖的小区通讯资费按正常资费的68%收取，有两个或两个以上基站信号覆盖的小区的通讯资费按正常收取，针对于5000万元的预算，应该如何建设基站，才能够使得资费的收入达到最大。
2、问题的分析
手机是通过在地面上建立了大量的无线基站来传递信号,达到通话目的。
若某手机运营商准备在一个目前尚未覆盖的区域开展业务，则需要考虑基站的覆盖能力，即某基站覆盖的那些社区以及社区的人数等问题，在此基础上建立基站网络，最大程度上服务于小区的居民。
根据题目条件，为了更好地分析问题，我们将基站对于小区的覆盖情况用下表来描述。
表3每个基站所能覆盖的社区
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
1
O
O
O
2
O
O
O
3
O
O
O
O
4
O
O
O
O
5
O
O
O
6
O
O
O
O
O
7
O
O
O
O
考虑到有的小区仅仅只有一个基站覆盖，因此要想实现所有社区的全面覆盖，有些基站是不能缺少的。
例如，1号、3 号、6 号、11 号、13号、14号社区均只可能有一个基站覆盖，那么为这些社区服务的基站是必不可少的。
因此，基站1号、2号、4号、6号、7号必须要设。
建设这些基站的费用
9.5+7+14+13+11=54.5>50;此时，仅仅必须建设的基站的费用已经不能满足要求。
因此，要想在实现不超过5000万建设费用的情况下实现对所有社区的覆盖是不可能的。
针对问题一：
建立0-1整数规划，通过对题目条件和问题的挖掘，列写出规模型中的目标函数和约束条件。
运用数学软件lingo求解，得到合理的基站建设方案。
针对问题二：
在满足基站建设成本不超过5000 万元的情况下，确定一个合理的基站建设方案，使得运营商的资费收入最高。
问题关键在于确定每一个社区用哪几个社区覆盖，然后计算根据题目中的“仅有一个基站信号覆盖的小区通讯资费按正常资费的68%收取，有两个或两个以上基站信号覆盖的小区的通讯资费按正常收取”的原则，可以列写出关于资费收入的函数表达式。
运用数学软件lingo最终把满足条件的基站建设方案解出，最终确定出最理想的基站建设方案。
3、模型的假设与符号的说明
3.1模型的假设
（1）若某社区处在某一基站覆盖范围内,则该社区中的人口全部被该基站覆盖；
（2）各社区的手机使用率相同；
（3）每位手机使用者的通讯资费相同；
（4）该区域只存在这一种通信网络；
（5）每个基站覆盖且仅覆盖图1所列出的覆盖区域；
（6）通讯信号不受地形地貌,气候变化等因素影响；
（7）社区人口保持不变；
（8）不考虑手机漫游等情况；
（9）每个基站位置最多只建一个基站。
3.2符号的说明
表示第i个基站建设情况（i=1,2,..7），当
=1时，表示第i个基站要被建设；当
=0时表示第i个基站不要被建设
表示第j个社区被覆盖情况（j=1,2,...15）,当
=1时，表示第j个社区被覆盖；当
=0时表示第j个社区未被覆盖
表示第j个社区的人口数（j=1,2,...15）
表示第i个基站被建设所需的费用（i=1,2,...7）
表示第j个社区被覆盖情况（j=1,2,...15）,当j=i，表示第j个社区被多个基站覆盖；当
=0.68时，表示第j个社区被1个基站覆盖；当
=0时表示第j个社区未被覆盖
4、模型的建立及求解
4.1模型的建立
问题一：
设
（i=1，2,...7表示7个中继站）表述每一个基站的建设情况。
引入0-1变量，即
=1,表示第i个基站要建立
0，表示第i个基站不建立
在此模型的建立过程中，由于同一个社区可能有多个基站覆盖，如果覆盖同一社区的基站都要建设时，那么基站覆盖的人口就会被重复计算。
故我们将目标转移到社区上，每个社区的被覆盖情况只有两种，要么被覆盖要么不被覆盖我们也引入0-1变量，即
=1,表示第j个社区被覆盖
0，表示第j个社区不被覆盖
这样就可避免了对同一社区人口的重复计算。
本问题的目标是使得基站覆盖的人口尽量多。
根据表1、2、3我们可以得到目标函数：
max=
由于考虑到1号、3 号、6 号、11 号、13号、14号社区均只可能有一个基站覆盖，这里我们让
代替
（即第j个社区只被第i个基站覆盖），则目标函数：
max=2*x1+4*（y2）+13*x2+6*（y4）+9*（y5）+4*x4+7.5*（y7）+12.5*（y8）+10*（y9）+11*（y10）+6*x6+14*（y12）+9*x7+3.5*x7+6*（y15）;
要求建设基站的费用不超过5000 万元
故约束条件：
（9.5*x1+7*x2+19*x3+14*x4+17.5*x5+13*x6+11*x7）<=50;
问题二：
题中考虑到基站出现故障维修的时候可能会出现所覆盖的社区信号中断等问题，为此对通讯资费进行了调整，规定，仅有一个基站信号覆盖的小区通讯资费按正常资费的68%收取，有两个或两个以上基站信号覆盖的小区的通讯资费按正常收取，为此，我们需要得到新的模型来进行求解，因为假设每个用户的正常资费相同，所以68%可以用减少人口来求最优值，与问题一类似，考虑到1号、3 号、6 号、11 号、13号、14号社区均只可能有一个基站覆盖，这里我们让
代替
（即第j个社区只被第i个基站覆盖），故问题二的目标函数为：
max=2*x1+4*（y2）+13*x2+6*（y4）+9*（y5）+4*x4+7.5*（y7）+12.5*（y8）+10*（y9）+11*（y10）+6*x6+14*（y12）+9*x7+3.5*x7+6*（y15）;
题目要求建设中继站的费用不超过5000 万元
故约束条件：
（9.5*x1+7*x2+19*x3+14*x4+17.5*x5+13*x6+11*x7）<=50;
在此方案下，获得的资费为：
s=2*x1*（k1）+4*（y2）*（k2）+13*x2*（k3）+6*（y4）*（k4）+9*（y5）*（k5）+4*x4*（k6）+7.5*（y7）*（k7）+12.5*（y8）*（k8）+10*（y9）*（k9）+11*（y10）*（k10）+6*x6*（k11）+14*（y12）*（k12）+9*x7*（k13）+3.5*x7*（k13）+6*（y15）*（k15）;
4.2模型的求解
问题一：
根据附录中的程序一利用LINGO求解得到最佳的方案如下表4所示：
表4
基站
1
2
3
4
5
6
7
建设情况
不建设
建设
不建设
建设
不建设
建设
建设
此方案所需费用为45百万元，覆盖人口为109.5千人。
问题二：
根据附录中的程序二利用LINGO求解得到最佳的方案如下表5所示：
表5
基站
1
2
3
4
5
6
7
建设情况
不建设
建设
不建设
建设
不建设
建设
建设
此方案所需要的费用为45百万元，获得资费83.74a（a为标准的资费常数）。
5、结果分析
对于问题一，要求在基站建设成本不超过50百万元的情况下，确定一个合理的基站建设方案，使得覆盖的人口尽可能的多。
所以我们根据题意建立了0-1规划模型，运用LONGO软件对规划模型求解，得到在2,4,6,7号位置建设基站时，覆盖人口最多为109.5千人，同时建设基站的费用为45百万元，满足约束条件中的费用不超过50百万的要求。
对于问题二，要求的是在满足基站建设成本不超过5000万元预算条件下，怎样建设基站，使得运营商的资费收入最高。
根据题目中“仅有一个基站信号覆盖的小区人均通讯资费按正常资费的68%收取，而有两个或两个以上站信号覆盖的小区人均的通讯资费按正常收取”的要求，我们运用了0-1规划方法，并且用lingo数学软件得出最大资费收益为S=83.74a。
6、优化方向
该模型巧妙的解决了相邻信号站重复覆盖的人口数的问题，使得LINGO求解方便，缺点是当数据量更大时计算会比较复杂，所以可以考虑用MATLAB编程求解，列出基站和小区的关系矩阵。
并且考虑问题时我们只考虑了两个重要的因素，因此，对于本问题的延伸，可更改规划目标，并加入更多的约束条件，如：
通过研究得出地区信号覆盖层数对信号质量的影响，继而影响用户数量及收费标准，在通过各种方法将对这些因素进行定量分析，建立合理的基站最大覆盖模型。
以最大收益为目标函数。
新问题的规划方法可以再上述模型为框架的基础上修改而得。
7、参考文献
[1].胡运权编著《运筹学教程》清华大学出版社2007.04 第三版；
[2].蒋启源编著《数学模型》高等教育出版社2003.08 第三版；
[3].吴礼斌，李柏年数学实验与建模[M]，北京：
国防工业出版社，2007年；
[4] 王兵团数学建模基础[M]，北京：
北京交通大学出版社，2004 年；
[5] 胡守信，李柏年基于MATLAB 的数学试验[M]，北京：
科学出版社，2004年；
[6]李明月移动通讯基站建设问题
2012.12.17/2015.07.02
附录：
程序一：
问题一
model:
max=2*x1+4*（y2）+13*x2+6*（y4）+9*（y5）+4*x4+7.5*（y7）+12.5*（y8）+10*（y9）+11*（y10）+6*x6+14*（y12）+9*x7+3.5*x7+6*（y15）;
（9.5*x1+7*x2+19*x3+14*x4+17.5*x5+13*x6+11*x7）<=50;
Y2=@if（x1+x2#eq#0,0,1）;
Y4=@if（x1+x3#eq#0,0,1）;
Y5=@if（x2+x4#eq#0,0,1）;
Y7=@if（x3+x6#eq#0,0,1）;
Y8=@if（x3+x4+x5#eq#0,0,1）;
Y9=@if（x4+x5#eq#0,0,1）;
Y10=@if（x3+x6#eq#0,0,1）;
Y12=@if（x5+x6+x7#eq#0,0,1）;
Y15=@if（x6+x7#eq#0,0,1）;
@bin（x1）;@bin（x2）;@bin（x3）;@bin（x4）;@bin（x5）;@bin（x6）;@bin（x7）;
end
运行结果：
Localoptimalsolutionfound.
Objectivevalue:
109.5000
Extendedsolversteps:
3
Totalsolveriterations:
185
VariableValueReducedCost
X10.000000-2.000000
Y21.0000000.000000
X21.000000-13.00000
Y40.0000000.000000
Y51.0000000.000000
X41.000000-4.000000
Y71.0000000.000000
Y81.0000000.000000
Y91.0000000.000000
Y101.0000000.000000
X61.000000-6.000000
Y121.0000000.000000
X71.000000-12.50000
Y151.0000000.000000
X30.0000000.000000
X50.0000000.000000
RowSlackorSurplusDualPrice
1109.50001.000000
25.0000000.000000
30.0000004.000000
40.0000006.000000
50.0000009.000000
60.0000007.500000
70.00000012.50000
80.00000010.00000
90.00000011.00000
100.00000014.00000
110.0000006.000000
程序二：
问题二
model:
max=2*x1+4*（y2）+13*x2+6*（y4）+9*（y5）+4*x4+7.5*（y7）+12.5*（y8）+10*（y9）+11*（y10）+6*x6+14*（y12）+9*x7+3.5*x7+6*（y15）;
（9.5*x1+7*x2+19*x3+14*x4+17.5*x5+13*x6+11*x7）<=50;
y2=@if（x1+x2#eq#0,0,1）;
y4=@if（x1+x3#eq#0,0,1）;
y5=@if（x2+x4#eq#0,0,1）;
y7=@if（x3+x6#eq#0,0,1）;
y8=@if（x3+x4+x5#eq#0,0,1）;
y9=@if（x4+x5#eq#0,0,1）;
y10=@if（x3+x6#eq#0,0,1）;
y12=@if（x5+x6+x7#eq#0,0,1）;
y15=@if（x6+x7#eq#0,0,1）;
k1=@if（x1#eq#1,0.68,0）;
k2=@if（x1+x2#eq#1,0.68,1）;
k3=@if（x2#eq#1,0.68,1）;
k4=@if（x1+x3#eq#1,0.68,0）;
k5=@if（x4+x2#eq#1,0.68,1）;
k6=@if（x4#eq#1,0.68,1）;
k7=@if（x3+x6#eq#1,0.68,1）;
k8=@if（x3+x4+x5#eq#1,0.68,1）;
k9=@if（x4+x5#eq#1,0.68,1）;
k10=@if（x3+x6#eq#1,0.68,1）;
k11=@if（x6#eq#1,0.68,1）;
k12=@if（x5+x6+x7#eq#1,0.68,1）;
k13=@if（x7#eq#1,0.68,1）;
k14=@if（x7#eq#1,0.68,1）;
k15=@if（x6+x7#eq#1,0.68,1）;
s=2*x1*（k1）+4*（y2）*（k2）+13*x2*（k3）+6*（y4）*（k4）+9*（y5）*（k5）+4*x4*（k6）+7.5*（y7）*（k7）+12.5*（y8）*（k8）+10*（y9）*（k9）+11*（y10）*（k10）+6*x6*（k11）+14*（y12）*（k12）+9*x7*（k13）+3.5*x7*（k13）+6*（y15）*（k15）;
@bin（x1）;@bin（x2）;@bin（x3）;@bin（x4）;@bin（x5）;@bin（x6）;@bin（x7）;
end
运行结果：
Localoptimalsolutionfound.
Objectivevalue:
109.5000
Extendedsolversteps:
0
Totalsolveriterations:
115
VariableValueReducedCost
X10.000000-2.000000
Y21.0000000.000000
X21.000000-13.00000
Y40.0000000.000000
Y51.0000000.000000
X41.000000-4.000000
Y71.0000000.000000
Y81.0000000.000000
Y91.0000000.000000
Y101.0000000.000000
X61.000000-6.000000
Y121.0000000.000000
X71.000000-12.50000
Y151.0000000.000000
X30.0000000.000000
X50.0000000.000000
K10.0000000.000000
K20.68000000.000000
K30.68000000.000000
K40.0000000.000000
K51.0000000.000000
K60.68000000.000000
K70.68000000.000000
K80.68000000.000000
K90.68000000.000000
K100.68000000.000000
K110.68000000.000000
K121.0000000.000000
K130.68000000.000000
K140.68000000.000000
K151.0000000.000000
S83.740000.000000
RowSlackorSurplusDualPrice
1109.50001.000000
25.0000000.000000
30.0000004.000000
40.0000006.000000
50.0000009.000000
60.0000007.500000
70.00000012.50000
80.00000010.00000
90.00000011.00000
100.00000014.00000
110.0000006.000000
120.0000000.000000
130.0000000.000000
140.0000000.000000
150.0000000.000000
160.0000000.000000
170.0000000.000000
180.0000000.000000
190.0000000.000000
200.0000000.000000
210.0000000.000000
220.0000000.000000
230.0000000.000000
240.0000000.000000
250.0000000.000000
260.0000000.000000
270.0000000.000000
东华理工大学
学年课程设计报告评分表
学生姓名：
张仕纬学号：
201320400223班级：
1324002
学生姓名：
揭小兰学号：
201320400215班级：
1324002
学生姓名：
林绵庄学号：
201320400109班级：
1324001
课程设计题目：
移动通讯基站建设问题
项目内容
满分
实评
选
题
能结合所学课程知识、有一定的能力训练。
符合选题要求
（3人一题）
5
工作量适中，难易度合理
10
能
力
水
平
能熟练应用所学知识，有一定查阅文献及运用文献资料能力
10
理论依据充分，数据准确，公式推导正确
10
能应用计算机软件进行编程、资料搜
升级会员 