数学建模作业
数学建模结课作业
一． 某旅游景点从山脚到山顶有一缆车索道，全长约1471m,高度 差为380m 。
采用循环单线修建，从下站到上站行经8个铁塔，将缆绳分为九段，各段的水平距离用i d 表示，高差用i h 表示，其数据见下表：每一段缆绳垂下来的最低点不低于两端铁塔最低塔顶悬挂绳处1m 。
要求：（1）折线法；（2）抛物线法，估计整个索道工程所用的缆绳总长度。
解：（一）折线法思路：考虑到实际中工程架线不能过紧，但又为了节省原料，我们采取求出最大折线和最小折线，对两者求取平均值，以得到对缆线总长度的估测。
由于八个铁塔分九段，因此此题分两部分考虑：（1） 第一段：直接求出发点到第一个铁塔的距离，即21211h d l +=（2） 第二到九段：建立坐标系，运用距离公式求取l 的长度。
设A （x -,1）,B(i d x -,1i h +)得：l =用此公式求最大最小值。
matlab 求解第一段syms h1 d1h1=50d1=220l1=sqrt(d1.^2+h1.^2)第二段求最小值clearl='sqrt((-x)^2+1)+sqrt((200-x)^2+(45+1)^2)' ezplot(l,[0,200]);[xmin,lmin]=fminbnd(l,0,200)得图形可得当x=4.2553时，取得最小值205.45由图形可得当x=200时取得最大值，即clearl='sqrt((-x)^2+1)+sqrt((200-x)^2+(45+1)^2)' ezplot(l,[0,200]);[xmin,lmin]=fminbnd(l,0,200)x=200;lmax=eval(l);l=(lmin+lmax)/2;得lmax=246.0025l=225.7254第三段到第九段算法与第二段相同，所以结果为第一段：l1 = 225.6103第二到九段分别为: 225.7254 ，163.5839 ，142.7476，120.6438，142.7476，163.5839，225.7254，248.5321总长为：1658.9m抛物线法思路：参照示意图，因为将绳的形状看做抛物线，为了方便研究，以抛物线的最低点为原点建立抛物线2y ax =，则每段绳的长度为l =，最后相加求总长。
数学建模作业
1.给出一个你所感兴趣的数学建模的实际问题。
（1）写出问题的实际背景。
（2）给出解答问题的建模与解答路径。
（3）解决什么样的问题。
答：（1）我们学校的教学楼中，教室的灯管的开关时间并没有一个明确的规定，这就造成了即使在大白天教室也开着灯的现象，浪费了很多的电力资源。
所以我们应该设计一个模型来对教室灯管的开关情况进行优化，以达到节省电力资源的同时又不影响同学们的正常学习。
（2）首先要统计出全校教学楼中共有多少个教室，以及每个教室的灯管的数量；其次要上网查资料，对西安一年四季的天气情况有一个初步的了解，分别统计出一年中雨天的比例，阴天的比例，和晴天的比例；最后查阅相关的资料，了解声控和光控开关的相关知识。
（3）通过建模来解决教室用电不合理的现象，即为学校节省了开销，也节约了电力资源，更可以通过这种潜移默化的形式，给同学们树立一个好榜样，使大家意识到节约用电的重要性。
2.找一本与本课有关的参考资料。
（1）你为何选择这一本书。
（2）这本资料对你的建模思想有什么启示作用。
（3）这本资料书对我么数学专业的学习有什么帮助。
答：我选择了《最优化方法》这本书。
（1）之所以选择这本书，首先是因为课堂上老师说数学建模里边，有很多问题都是要对某个问题进行优化的；其次是，随着科学技术的日益进步和生产经营的日益发展，最优化方法已成为现代管理科学的重要理论基础和不可缺少的方法，被人们广泛地应用到公共管理、经济管理、国防等各个领域，发挥着越来越重要的作用。
所以我选择了《最优化方法》这本书。
（2）这本书中主要是介绍线性规划问题的模型、求解及其应用――运输问题；以及动态规划的模型、求解、应用――资源分配问题。
其中的微分学中求极值、等式约束最优化问题、不等式约束最优化问题对数学建模都有很大的帮助。
用最优化方法解决实际问题，一般可经过下列步骤：①提出最优化问题，收集有关数据和资料；②建立最优化问题的数学模型,确定变量,列出目标函数和约束条件；③分析模型，选择合适的最优化方法；④求解，一般通过编制程序，用计算机求最优解；⑤最优解的检验和实施。
数学建模课后习题作业
【陈文滨】1、在稳定的椅子问题中，如设椅子的四脚连线呈长方形，结论如何？【模型假设】（1）椅子四条腿一样长，椅脚与地面接触处视为一点，四脚的连线呈长方形．（2）地面高度是连续变化的，沿任何方向都不会出现间断 (没有像台阶那样的情况)，即从数学的角度看，地面是连续曲面．这个假设相当于给出了椅子能放稳的必要条件．（3）椅子在任何位置至少有三只脚同时着地．为保证这一点，要求对于椅脚的间距和椅腿的长度而言，地面是相对平坦的．因为在地面上与椅脚间距和椅腿长度的尺寸大小相当的范围内，如果出现深沟或凸峰(即使是连续变化的)，此时三只脚是无法同时着地的。
【模型建立】在上述假设下，解决问题的关键在于选择合适的变量，把椅子四只脚同时着地表示出来．首先，引入合适的变量来表示椅子位置的挪动．生活经验告诉我们，要把椅子通过挪动放稳，通常有拖动或转动椅子两种办法，也就是数学上所说的平移与旋转变换．然而，平移椅子后问题的条件没有发生本质变化，所以用平移的办法是不能解决问题的．于是可尝试将椅子就地旋转，并试图在旋转过程中找到一种椅子能放稳的情形．注意到椅脚连线呈长方形，长方形是中心对称图形，绕它的对称中心旋转180度后，椅子仍在原地．把长方形绕它的对称中心O旋转，这可以表示椅子位置的改变。
于是，旋转角度θ这一变量就表示了椅子的位置．为此，在平面上建立直角坐标系来解决问题．如下图所示，设椅脚连线为长方形ABCD，以对角线AC所在的直线为x轴，对称中心O为原点，建立平面直角坐标系．椅子绕O点沿逆时针方向旋转角度θ后，长方形ABCD转至A1B1C1D1 的位置，这样就可以用旋转角θ（0≤θ≤π）表示出椅子绕点O旋转θ后的位置．其次，把椅脚是否着地用数学形式表示出来．我们知道，当椅脚与地面的竖直距离为零时，椅脚就着地了，而当这个距离大于零时，椅脚不着地．由于椅子在不同的位置是θ的函数，因此，椅脚与地面的竖直距离也是θ的函数．由于椅子有四只脚，因而椅脚与地面的竖直距离有四个，它们都是θ的函数．而由假设(3)可知，椅子在任何位置至少有三只脚同时着地，即这四个函数对于任意的θ，其函数值至少有三个同时为0．因此，只需引入两个距离函数即可．考虑到长方形ABCD是中心对称图形，绕其对称中心 O沿逆时针方向旋转180°后，长方形位置不变，但A,C和B,D对换了．因此，记A、B两脚与地面竖直距离之和为f（θ），C、D两脚与地面竖直距离之和为g（θ），其中θ∈[0，π]，从而将原问题数学化。
数学建模作业及答案
数学建模作业姓名：叶勃学号：班级：024121一：层次分析法1、 分别用和法、根法、特征根法编程求判断矩阵1261/2141/61/41A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦11/2433217551/41/711/21/31/31/52111/31/5311A ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦的特征根和特征向量（1）冪法求该矩阵的特征根和特征向量 程序为：#include<iostream> #include<math.h> using namespace std;#define n 3 //三阶矩阵#define N 20 #define err 0.0001 //幂法求特征值特征向量 void main(){cout<<"**********幂法求矩阵最大特征值及特征向量***********"<<endl; int i,j,k;double A[n][n],X[n],u,y[n],max;cout<<"请输入矩阵:\n"; for(i=0;i<n;i++) for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵 cout<<"请输入初始向量:\n"; for(i=0;i<n;i++)cin>>X[i]; //输入初始向量 k=1; u=0;while(1){ max=X[0]; for(i=0;i<n;i++) {if(max<X[i]) max=X[i]; //选择最大值 }for(i=0;i<n;i++)y[i]=X[i]/max; for(i=0;i<n;i++)X[i]=0;for(j=0;j<n;j++)X[i]+=A[i][j]*y[j]; //矩阵相乘}if(fabs(max-u)<err){cout<<"A的特征值是 :"<<endl; cout<<max<<endl; cout<<"A的特征向量为:"<<endl; for(i=0;i<n;i++) cout<<X[i]/(X[0]+X[1]+X[2])<<" ";cout<<endl;break;}else{if(k<N) {k=k+1;u=max;} else {cout<<"运行错误\n";break;}}} }程序结果为：（2）和法求矩阵最大特征值及特征向量程序为：#include<stdio.h>#include<iostream>#include<math.h> using namespace std;#define n 3 //三阶矩阵#define N 20void main(){int i,j,k;double A[n][n],w[n],M[n],u[n],W[n][n],max;cout<<"********和法求矩阵的特征根及特征向量*******"<<endl;cout<<"请输入矩阵:\n";for(i=0;i<n;i++)for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵 //计算每一列的元素和M[0]=0;M[1]=0;M[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){M[i]+=A[j][i];}//将每一列向量归一化for(i=0;i<n;i++)for(j=0;j<n;j++){W[j][i]=A[j][i]/M[i];}//输出按列归一化之后的矩阵Wcout<<"按列归一化后的矩阵为："<<endl;for(i=0;i<n;i++)for(j=0;j<n;j++){cout<<W[i][j]<<" ";if(j==2)cout<<endl;} //求特征向量w[0]=0;w[1]=0;w[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){w[i]+=W[i][j];}cout<<"特征向量为:"<<endl; for(i=0;i<n;i++){u[i]=w[i]/(w[0]+w[1]+w[2]);cout<<u[i]<<" "<<endl;}//求最大特征值max=0;for(i=0;i<n;i++){w[i] = 0;for(j=0;j<n;j++){w[i] += A[i][j]*u[j];}}for(i = 0;i < n;i++){max += w[i]/u[i];}cout<<"最大特征根为："<<endl;cout<<max/n<<endl; }运行结果为：（3）根法求矩阵最大特征值及特征向量：程序为：#include<stdio.h>#include<iostream>#include<math.h>using namespace std;#define n 3 //三阶矩阵#define N 20void main(){int i,j;double A[n][n],w[n],M[n],u[n],W[n][n],max;cout<<"********根法求矩阵的特征根及特征向量*******"<<endl; cout<<"请输入矩阵:\n";for(i=0;i<n;i++)for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵//计算每一列的元素和M[0]=0;M[1]=0;M[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){M[i]+=A[j][i];}//将每一列向量归一化for(i=0;i<n;i++)for(j=0;j<n;j++){W[j][i]=A[j][i]/M[i];}//输出按列归一化之后的矩阵Wcout<<"按列归一化后的矩阵为："<<endl;for(i=0;i<n;i++)for(j=0;j<n;j++){cout<<W[i][j]<<" ";if(j==2)cout<<endl;}//求特征向量//w[0]=A[0][0];w[1]=A[0][1];w[2]=A[0][2];w[0]=1;w[1]=1;w[2]=1;for(i=0;i<n;i++){for(j=0;j<n;j++){w[i]=w[i]*W[i][j];}w[i]=pow(w[i], 1.0/3);}cout<<"特征向量为:"<<endl;for(i=0;i<n;i++){u[i]=w[i]/(w[0]+w[1]+w[2]);cout<<u[i]<<" "<<endl;}//求最大特征值max=0;for(i=0;i<n;i++){w[i] = 0;for(j=0;j<n;j++){w[i] += A[i][j]*u[j];}}for(i = 0;i < n;i++){max += w[i]/u[i];}cout<<"最大特征值为:"<<endl; cout<<max/n;}运行结果为：2、编程验证n阶随机性一致性指标RI：运行结果：3、考虑景色、费用、居住、饮食、旅途五项准则，从桂林、黄山、北戴河三个旅游景点选择最佳的旅游地。
数学建模一周作业题目
对作业题目的说明1. 本次数学建模周一共提供十五道题目供大家选择。
每支队伍（2-3人/队）必须从以下题目中任意选取一题（只须选择一道），并完成一篇论文，对论文的具体要求参阅《论文格式规范》。
2. 题目标注为“A ”的为有一定难度的题目，指导老师会根据题目的难度对论文最后的评分进行调整。
（一）乒乓球赛问题 (A)A 、B 两乒乓球队进行一场五局三胜制的乒乓球赛，两队各派3名选手上场，并各有3种选手的出场顺序（分别记为123,,ααα 和123,,βββ）。
根据过去的比赛记录，可以预测出如果A 队以i α次序出场而B 队以j β次序出场，则打满5局A 队可胜ija 局。
由此得矩阵()ij R a =如下：123123214034531R βββααα⎛⎫=⎪ ⎪ ⎪⎝⎭（1） 根据矩阵R 能看出哪一队的实力较强吗？（2） 如果两队都采取稳妥的方案，比赛会出现什么结果？ （3） 如果你是A 队的教练，你会采取何种出场顺序？（4） 比赛为五战三胜制，但矩阵R 中的元素却是在打满五局的情况下得到的，这样的数据处理和预测方式有何优缺点？（二）野兔生长问题在某地区野兔的数量在连续十年的统计数量(单位十万)如下：分析该数据，得出野兔的生长规律。
并指出在哪些年内野兔的增长有异常现象，预测T=10 时野兔的数量。
（三）停车场的设计问题在New England的一个镇上，有一位于街角处面积100 200平方英尺的停车场，场主请你代为设计停车车位的安排方式，即设计在场地上划线的方案。
容易理解，如果将汽车按照与停车线构成直角的方向，一辆紧挨一辆地排列成行，则可以在停车场内塞进最大数量的汽车，但是对于那些缺乏经验的司机来说，按照这种方式停靠车辆是有困难的，它可能造成昂贵的保险费用支出。
为了减少因停车造成意外损失的可能性，场主可能不得不雇佣一些技术熟练的司机专门停车；另一方面，如果从通道进入停车位有一个足够大的转弯半径，那么，看来大多数的司机都可以毫无困难地一次停车到位。
数学建模作业题+答案
数学建模MATLAB 语言及应用上机作业11. 在matlab 中建立一个矩阵135792468101234501234A ⎡⎤⎢⎥⎢⎥=⎢⎥-----⎢⎥⎣⎦答案：A = [1,3,5,7,9;2,4,6,8,10;-1,-2,-3,-4,-5;0,1,2,3,4]2. 试着利用matlab 求解出下列方程的解（线性代数22页例14）123412423412342583692254760x x x x x x x x x x x x x x +-+=⎧⎪--=⎪⎨-+=-⎪⎪+-+=⎩ 答案：A=[2 ,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6]; B=[8;9;-5;0]; X=A\B 或A=[2,1,-5,1;1,-3,0,-6;0,2,-1,2;1,4,-7,6] b=[8,9,-5,0]' X=inv(A)*b3. 生成一个5阶服从标准正态分布的随机方阵，并计算出其行列式的值，逆矩阵以及转置矩阵。
答案:A=randn(5) det(A) inv(A) A'4. 利用matlab 求解出110430002A -⎡⎤⎢⎥=-⎢⎥⎢⎥⎣⎦的特征值和特征向量。
答案：A=[-1,1,0;-4,3,0;0,0,2] [V,D]=eig(A)5.画出衰减振荡曲线3sin3t y et -=在[0,4]π上的图像。
要求，画线颜色调整为黑色，画布底面为白色。
（在实际中，很多打印机时黑白的，因此大多数作图要考虑黑白打印机的效果。
） 给出恰当的x ，y 坐标轴标题，图像x 轴的最大值为4π。
6. 生成一个0-1分布的具有10个元素的随机向量，试着编写程序挑选出向量中大于0.5的元素。
数学建模和Matlab 上机作业2（2016-9-20）跟老师做（不用整合进作业中）：上机演示讲解：函数，递归的两个例子的写法。
附：1. Fibonacci Sequence （斐波那契数列）在数学上，费波那西数列是以递归的方法来定义： F1= 1；F2= 1；F （n ）=F （n-1）+F （n-2） 2. 阶乘举例：数学描述：n!=1×2×……×n ；计算机描述：n!=n*(n-1)!自己做（需要整合进作业中，提交到系统中）：1. 写一个m 文件完成分值百分制到5分制的转换（即输入一个百分制，转换后输出一个5级对应的得分，联系条件控制语句）。
数学建模作业(1)
数学建模作业(1)
数模
数模
1.学校共学校共1000名学生，235人住在宿名学生，人住在A宿名学生人住在人住B宿舍人住在C宿舍舍，333人住宿舍，432人住在宿舍人住宿舍，人住在宿舍.学生们要组织一个10人的委员会人的委员会，学生们要组织一个人的委员会，试用下列办法分配各宿舍的委员数：列办法分配各宿舍的委员数：(1)按比例分配取整数的名额后，剩下的名按比例分配取整数的名额后，按比例分配取整数的名额后额按惯例分给小数部分较大者。
额按惯例分给小数部分较大者。
(2)用Q值方法。
值方法。
用值方法
数模
如果委员会从10人增至人如果委员会从人增至15人，用以上人增至2种方法再分配名额。
将2种方法两次分配种方法再分配名额。
种方法再分配名额种方法两次分配的结果列表比较。
的结果列表比较。
(3)你能提出其它的方法吗？用你的方你能提出其它的方法吗？你能提出其它的方法吗法分配上面的名额。
法分配上面的名额。
数模
2.考察模拟水下爆炸的比例模型.爆炸物质量m,在距爆炸点距离r处设置仪器，接收到的冲击波压强为p,记大气初始压强p0,水的密度ρ,水的体积弹性模量k,用量纲分析法已经得到
p0ρrp=p0(,)km3
数模
设模拟实验与现场的p0,ρ,k相同，而爆炸物模型的质量为原模型的1/1000.为了使实验中接收到与现场相同的压强p,问实验时应如何设置接收冲击波的仪器，即求实验仪器与爆炸点之间的距离是现场的多少倍？
p0,ρ,k。
数学建模作业及答案
数学建模作业姓名：叶勃学号：班级：024121一：层次分析法1、 分别用和法、根法、特征根法编程求判断矩阵1261/2141/61/41A ⎡⎤⎢⎥=⎢⎥⎢⎥⎣⎦11/2433217551/41/711/21/31/31/52111/31/5311A ⎡⎤⎢⎥⎢⎥⎢⎥=⎢⎥⎢⎥⎢⎥⎣⎦的特征根和特征向量（1）冪法求该矩阵的特征根和特征向量 程序为：#include<iostream> #include<math.h> using namespace std;#define n 3 //三阶矩阵#define N 20 #define err 0.0001 //幂法求特征值特征向量 void main(){cout<<"**********幂法求矩阵最大特征值及特征向量***********"<<endl; int i,j,k;double A[n][n],X[n],u,y[n],max;cout<<"请输入矩阵:\n"; for(i=0;i<n;i++) for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵 cout<<"请输入初始向量:\n"; for(i=0;i<n;i++)cin>>X[i]; //输入初始向量 k=1; u=0;while(1){ max=X[0]; for(i=0;i<n;i++) {if(max<X[i]) max=X[i]; //选择最大值 }for(i=0;i<n;i++)y[i]=X[i]/max; for(i=0;i<n;i++)X[i]=0;for(j=0;j<n;j++)X[i]+=A[i][j]*y[j]; //矩阵相乘}if(fabs(max-u)<err){cout<<"A的特征值是 :"<<endl; cout<<max<<endl; cout<<"A的特征向量为:"<<endl; for(i=0;i<n;i++) cout<<X[i]/(X[0]+X[1]+X[2])<<" ";cout<<endl;break;}else{if(k<N) {k=k+1;u=max;} else {cout<<"运行错误\n";break;}}} }程序结果为：（2）和法求矩阵最大特征值及特征向量程序为：#include<stdio.h>#include<iostream>#include<math.h> using namespace std;#define n 3 //三阶矩阵#define N 20void main(){int i,j,k;double A[n][n],w[n],M[n],u[n],W[n][n],max;cout<<"********和法求矩阵的特征根及特征向量*******"<<endl;cout<<"请输入矩阵:\n";for(i=0;i<n;i++)for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵 //计算每一列的元素和M[0]=0;M[1]=0;M[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){M[i]+=A[j][i];}//将每一列向量归一化for(i=0;i<n;i++)for(j=0;j<n;j++){W[j][i]=A[j][i]/M[i];}//输出按列归一化之后的矩阵Wcout<<"按列归一化后的矩阵为："<<endl;for(i=0;i<n;i++)for(j=0;j<n;j++){cout<<W[i][j]<<" ";if(j==2)cout<<endl;} //求特征向量w[0]=0;w[1]=0;w[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){w[i]+=W[i][j];}cout<<"特征向量为:"<<endl; for(i=0;i<n;i++){u[i]=w[i]/(w[0]+w[1]+w[2]);cout<<u[i]<<" "<<endl;}//求最大特征值max=0;for(i=0;i<n;i++){w[i] = 0;for(j=0;j<n;j++){w[i] += A[i][j]*u[j];}}for(i = 0;i < n;i++){max += w[i]/u[i];}cout<<"最大特征根为："<<endl;cout<<max/n<<endl; }运行结果为：（3）根法求矩阵最大特征值及特征向量：程序为：#include<stdio.h>#include<iostream>#include<math.h>using namespace std;#define n 3 //三阶矩阵#define N 20void main(){int i,j;double A[n][n],w[n],M[n],u[n],W[n][n],max;cout<<"********根法求矩阵的特征根及特征向量*******"<<endl; cout<<"请输入矩阵:\n";for(i=0;i<n;i++)for(j=0;j<n;j++)cin>>A[i][j]; //输入矩阵//计算每一列的元素和M[0]=0;M[1]=0;M[2]=0;for(i=0;i<n;i++)for(j=0;j<n;j++){M[i]+=A[j][i];}//将每一列向量归一化for(i=0;i<n;i++)for(j=0;j<n;j++){W[j][i]=A[j][i]/M[i];}//输出按列归一化之后的矩阵Wcout<<"按列归一化后的矩阵为："<<endl;for(i=0;i<n;i++)for(j=0;j<n;j++){cout<<W[i][j]<<" ";if(j==2)cout<<endl;}//求特征向量//w[0]=A[0][0];w[1]=A[0][1];w[2]=A[0][2];w[0]=1;w[1]=1;w[2]=1;for(i=0;i<n;i++){for(j=0;j<n;j++){w[i]=w[i]*W[i][j];}w[i]=pow(w[i], 1.0/3);}cout<<"特征向量为:"<<endl;for(i=0;i<n;i++){u[i]=w[i]/(w[0]+w[1]+w[2]);cout<<u[i]<<" "<<endl;}//求最大特征值max=0;for(i=0;i<n;i++){w[i] = 0;for(j=0;j<n;j++){w[i] += A[i][j]*u[j];}}for(i = 0;i < n;i++){max += w[i]/u[i];}cout<<"最大特征值为:"<<endl; cout<<max/n;}运行结果为：2、编程验证n阶随机性一致性指标RI：运行结果：3、考虑景色、费用、居住、饮食、旅途五项准则，从桂林、黄山、北戴河三个旅游景点选择最佳的旅游地。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
数学建模作业：成靖学号：1408030311班级：计科1403班日期：2015.12.301.某班准备从5名游泳队员中选4人组成接力队，参加学校的4×100m混合泳接力比赛，5名队员4种泳姿的百米平均成绩如下表所示，问应如何选拔队员组成接力队？如果最近队员丁的蛙泳成绩有较大的退步，只有1′15"2；而队员戊经过艰苦训练自由泳成绩有所进步，达到57"5，组成接力队的方案是否应该调整？名队员4种泳姿的百米平均成绩ij若参选择队员i加泳姿j 的比赛，记x ij=1, 否则记x ij=0目标函数:即min=66.8*x11+75.6*x12+87*x13+58.6*x14+57.2*x21+66*x22+66.4*x23+53*x24 +78*x31+67.8*x32+84.6*x33+59.4*x34+70*x41+74.2*x42+69.6*x43+57.2*x44+ 67.4*x51+71*x52+83.8*x53+62.4*x54;约束条件: x11+x12+x13+x14<=1;x21+x22+x23+x24<=1;x31+x32+x33+x34<=1;x41+x42+x43+x44<=1;x51+x52+x53+x54<=1;x11+x21+x31+x41+x51=1;x12+x22+x32+x42+x52=1;x13+x23+x33+x43+x53=1;x14+x24+x34+x44+x54=1;∑∑===4151j iijijxcZMinlingo模型程序和运行结果因此，最优解为x14=1,x21=1,x32=1,x43=1,其余变量为0 成绩为253.2(秒)=4′13"2即：甲~ 自由泳、乙~ 蝶泳、丙~ 仰泳、丁~ 蛙泳.(2).若丁的蛙泳成绩退步为1′15"2=75.2(秒)，戊的自由泳成绩进步为57"5=57.5(秒)，则目标函数：min=66.8*x11+75.6*x12+87*x13+58.6*x14+57.2*x21+66*x22+66.4*x23+53*x24+78*x31+67.8*x32+84.6*x 33+59.4*x34+70*x41+74.2*x42+75.2*x43+57.2*x44+67.4*x51+71*x52+83.8*x53+57.5*x54;约束条件：x11+x12+x13+x14<=1;x21+x22+x23+x24<=1;x31+x32+x33+x34<=1;x41+x42+x43+x44<=1;x51+x52+x53+x54<=1;x11+x21+x31+x41+x51=1;x12+x22+x32+x42+x52=1;x13+x23+x33+x43+x53=1;x14+x24+x34+x44+x54=1lingo模型程序和运行结果因此，最优解为x21=1,x32=1,x43=1,x54=1 ,其余变量为0;成绩为257.7(秒)= 4′17"7 ,新方案:乙~ 蝶泳、丙~ 仰泳、丁~ 蛙泳、戊~ 自由泳。
2.某工厂用A1，A2两台机床加工B1，B2，B3三种不同零件，已知在一个生产周期A1只能工作80机时，A2只能工作100机时。
一个生产周期加工B1为70件，B2为50件，B3为20件。
两台机床加工每个零件的时间和加工每个零件的成本，分别如下所示加工每个零件时间表（单位:机时/个）加工每个零件成本表(单位：元/个)问怎样安排两台车床一个周期的加工任务，才能使加工成本最低？解：设在A1机床上加工零件B1、B2、B3的数量分别为x1、x2、x3，在A2机床上加工零件B1、B2、B3的数量分别为x4、x5、x6,可建立以下线性规划模型：目标函数：min=2*x1+3*x2+5*x3+3*x4+3*x5+6*x6约束条件：x1,x2,x3,x4,x5,x6均为整数x1+2*x2+3*x3<=80x1+x2+3*x3<=100x1+x4=70x2+x5=50x3+x6=20lingo模型程序和运行结果最优解为x1=70,x2=0,x3=3,x4=0,x5=50,x6=17;最低成本价为407元。
即：在A1机床上加工零件B1、B2、B3的数量分别为70、0、3，在A2机床上加工零件B1、B2、B3的数量分别为0、50、17。
3.某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券以及其信用等级、到期年限、收益如下表所示按照规定，市政证券的收益可以免税，其他证券的收益需按50%的税率纳税。
此外还有以下限制：（1）政府及代办机构的证券总共至少要购进400万元；（2）所购证券的平均信用等级不超过1.4（信用等级数字越小，信用程度越高）；（3）所购证券的平均到期年限不超过5年。
（1）若该经理有1000万元资金，应如何投资?（2）如果能够以2.75%的利率借到不超过100万元资金，该经理应如何操作?（3）在1000万元资金情况下，若证券A的税前收益增加为4.5%，投资应否改变?若证券C的税前收益减少为4.8%，投资应否改变?解：设投资证劵A、证劵B、证劵C、证劵D、证劵E的金额分别为：X1、X2、X3、X4、X5（百万元），投资之后获得的总收益为Y百万元。
(1).建立如下的线性规划模型：目标函数：maxY=0.043*X1+(0.054*0.5)*X2+(0.05*0.5)*X3+(0.044*0.5)*X4+0.045*X5 约束条件: X2+X3+X4>=4X1+X2+X3+X4+X5<=10(2*X1+2*X2+X3+X4+5*X5)/( X1+X2+X3+X4+X5)<=1.4(9*X1+15*X2+4*X3+3*X4+2*X5)/( X1+X2+X3+X4+X5)<=5整理化简可得：Max Y=0.043*X1+0.027*X2+0.025*X3+0.022*X4+0.045*X5；X2+X3+X4>=4；X1+X2+X3+X4+X5<=10；6*X1+6*X2-4*X3-X4+36*X5<=0；4*X1+10*X2-X3-2*X4-3*X5<=0；lingo模型程序和运行结果因此，最优解为Y=0.298，X1=2.182，X3=7.364，X5=0.454最优解方案不投资证劵B和证劵D，投资证劵A为218.2万元，投资证劵C为736.4万元，投资证劵E为45.4万元;总收益为29.8万元。
（2）.由问题（1）得:投资金额每增加100万元，收益可增加2.98万元，而借贷100万元所要支付的利息是2.75万元，比2.98万元少，因此应该借贷这100万元去投资。
目标函数仍为：Max Y=0.043*X1+0.027*X2+0.025*X3+0.022*X4+0.045*X5；X2+X3+X4>=4；X1+X2+X3+X4+X5<=11；6*X1+6*X2-4*X3-X4+36*X5<=0；4*X1+10*X2-X3-2*X4-3*X5<=0；lingo模型程序和运行结果因此，最优解为: X1=2.40,X3=8.10,X5=0.50,Y=0.328;即应投资证劵A 240万元，证劵C 810万元，证劵E 50万元。
此时收益总额为32.8万元，再减去所要支付的利息2.75万元，还剩30.05万元，比问题(1)中的收益总额29.8万元还要多，这也证明了借贷100万元来投资是明智的。
（3). 问题(1)的灵敏度分析可得下图：则在最优解不变的条件下目标函数系数允许的变化围：X1的系数为（0.043-0.013，0.043+0.0035），即（0.030，0.0465）；X3的系数为（0.025-0.0006，0.025+0.017），即（0.02494，0.042）；当证劵A的税前收益增加为4.5%时，其在目标函数中的系数为0.045，在最优解不变的条件下目标函数系数允许的变化围，因此投资方案不应该改变。
当证劵C的税前收益减少为4.8%时，其在目标函数中的系数为0.024，不在最优解不变的条件下目标函数系数允许的变化围，因此只有改变投资方案，才能使银行经理获得最大收益值。
4.某医院负责人每日至少需要下表数量的护士。
每班的护士在值班开始时向病房报到，连续工作8小时，医院领导为满足每班所需要的护士数，最少需要用多少护士？解：设在i班刚加入工作的人数分别为x1、x2、x3、x4、x5、x6；目标函数为：min=x1+x2+x3+x4+x5+x6;约束条件：x1,x2,x3,x4,x5,x6均为整数x1+x2>=70x2+x3>=60x3+x4>=50x4+x5>=20x5+x6>=30x6+x1>=60lingo模型程序和运行结果因此，最优解为：x1=60、x2=10、x3=50、x4=0、x5=30、x6=0；最少需要护士150人。
5．某海岛上有12个主要的居民点，每个居民点的位置（用平面坐标x,y表示，距离单位：km）和居住的人数R如表下表所示，现在准备在岛上建一个服务中心为居民提供各种服务，那么服务中心应该建在何处？解：设第i 个居民点的位置（x i ，x j ），居住的人数为R i ，i=1,2,3,4,5,6,7,8,9,10,11,12;设服务中心的位置为（a ，b ）,无约束条件；服务中心应该让所有的人都方便，因此目标函数为min=∑R i 12i =1√(x i −a )2 +(x j −b )2lingo 模型程序和运行结果因此，服务中心应该建的位置是（3.19，3.20）第十一个小岛。
6.某厂向用户提供发动机，合同规定，第一、二、三季度末分别交货40台、60台、80台。
每季度的生产费用为()2bx ax x f +=(元)，其中x 是该季生产的发动机台数，若交货后有剩余，可用于下季度交货，但需支付存储费，每台每季度c 元。
已知工厂每季度最大生产能力为100台，第一季度开始无存货，设a=50，b=0.2，c=4，问工厂应如何安排生产计划，才能既满足合同有使总费用最低？讨论a 、b 、c 、变化对计划的影响，并作出合理的解释。
解：(1).设工厂第一季度生产x1台发动机，第二季度生产x2台发动机，第三季度生产x3台发动机。
目标函数：min=50*x1+0.2*x1^2+50*x2+0.2*x2^2+50*x3+0.2*x3^2+4*(x1-40)+4*(x1+x2-100); 约束条件：x 1,x2,x3均为整数 x1<=100;x2<=100; x3<=100; x1>=40;x1+x2>=100; x1+x2+x3>=180;lingo 模型程序和运行结果因此，最优解为：x1=50,x2=60,x3=70;即：工厂第一季度生产50台发动机，第二季度生产60台发动机，第三季度生产70台发动机。