数学建模作业题

合集下载

数学建模作业练习

数学建模作业练习

优化作业(1)1.(本题只写模型不求解)某工厂向用户提供发动机,按合同规定,其交货数量和日期是:第一季度末交40台,第二季度末交60台,第三季度末交80台。

工厂的最大生产能力为每季度100台,每季度的生产费用是22.050)(x x x f +=元,其中x 为该季度生产发动机的台数。

若工厂生产得多,多余的发动机可移到下季度向用户交货,这样,工厂就需要支付存储费用,每台发动机每季度的存储费用为4元。

问该厂每季度生产多少台发动机,才能既满足交货合同,又使工厂所花费的费用最少(假定第一季度开始时发动机无存货)?2.(本题只写模型不求解)某市为方便小学生上学,拟在新建的8个居民小区821,,,A A A 增设若干所小学,经过论证知备选校址有621,,,B B B ,它们能够覆盖的居民小区如下表所列,试建立一个数学模型,确定出最小个数的建校地址,使其能覆盖所有的居民小区。

备选校址B 1 B 2 B 3 B 4 B 5 B 6 覆盖小区 A 1,A 5,A 7 A 1,A 2,A 5,A 8 A 1,A 3,A 5 A 2,A 4,A 8 A 3,A 6 A 4,A 6,A 83.写出下面LINGO 程序所对应的完整数学模型。

SETS: HANG/1..3/:B; LIE/1..4/:X,C; XISHU(HANG,LIE):A;ENDSETSDATA:A= 1 2 3 12 5 1 23 1 6 -2;B=4 5 7;C=1 3 4 5;ENDDATAmin=@sum(LIE(I):C(I)*X(I));@FOR(HANG(I):@SUM(LIE(J):A(I,J)*X(J))>B(I));4.根据下面LINGO 程序的集合段和模型段写出其所对应的数学模型。

SETS: HANG/1..3/:A;LIE/1..4/:B;XISHU(HANG,LIE):C,X;ENDSETSmin=@sum(XISHU(I,J):C(I,J)*X(I,J));@FOR(HANG(I):@SUM(LIE(J):X(I,J))=A(I));@FOR(LIE(J):@SUM(HANG(I):X(I,J))=B(J));5.某校篮球队准备从十名预备队员中选择五名作为正式队员,队员的各种情况如下表:队员号码身高(厘米)技术分位置1 185 8.6 中锋2 186 9 中锋3 193 8.4 中锋4 190 9.5 中锋5 182 9.1 前锋6 184 9 前锋7 188 8.1 前锋8 186 7.8 后卫9 190 8.2 后卫10 192 9.2 后卫队员的挑选要满足下面条件:(1)至少补充一名前锋。

数学建模作业题+答案

数学建模作业题+答案

数学建模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级对应的得分,联系条件控制语句)。

简单数学建模100例

简单数学建模100例

“学”以致用-----简单数学建模应用问题100例数学教学过程中学习了一个数学公式后,需要做大量的应用题,通过训练来加深理解所学公式。

但是在生活中又有多少实际问题是可以直接套用公式的呢?理想状态下的公式直接运用,在生产及生活中的实例是少之又少。

为此学生总感到学了数学没有什么实际用处,所以对学习数学少有兴趣。

数学建模的引入对培养学生利用数学方法分析、解决实际问题的能力开辟了一条有效的途径,让中职学生从中体会到数学是来源于生活并应用于生活的.数学建模是一种思维方式,它是一个动态的过程,通过此过程可以将一个实际的问题,经过模型准备、模型假设、模型构成、模型解析、模型检验与应用等五个具体步骤,转变为可以用数学方法(公式)来解决的,在理想状态下的数学问题,上述的整个流程统称为数学建模如果想解决某个实际问题(也许它和数学没有直接的关系),可以按下面流程对问题进行数学建模。

一.模型准备先了解该问题的实际背景和建模目的,尽量弄清要建模的问题属于哪一类学科的问题,可能需要用到哪些知识,然后学习或复习有关的知识,为接下来的数学建模做准备.由于人们所掌握的专业知识是有限的,而实际问题往往是多样和复杂的,模型准备对做好数学建模问题是非常重要的.二.模型假设有了模型准备的基础,要想把实际问题变为数学问题还要对其进行必要合理的简化和假设.明确了建模目的又掌握了相关资料,再去除一些次要因素.以主要矛盾为主来对该实际问题进行适当的简化并提出一些合理的假设。

模型假设不太可能一蹴而就,可以在模型的不断修改中得到逐步完善.三.模型构成在模型假设的基础上,选择适当的数学工具并根据已知的知识和搜集的信息来描述变量之间的关系或其他数学结构(如数学公式、定理、算法等).做模型构成时可以使用各种各样的数学理论和方法,但要注意的是在保证精度的条件下尽量用简单的数学方法是建模时要遵循的一个原则.四.模型解析在模型构成中建立的数学模型可以采用解方程、推理、图解、计算机模拟、定理证明等各种传统的和现代的数学方法对其进行求解,其中有些可以借助于计算机软件来做这些工作。

数学建模作业(1)

数学建模作业(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。

数学建模小作业例题

数学建模小作业例题

数学建模小作业例题1. 在冷却过程中,物体的温度在任何时刻变化的速率大致正比于它的温度与周围介质温度之差,这一结论称为牛顿冷却定律,该定律同样用于加热过程。

一个煮硬了的鸡蛋有98℃,将它放在18℃的水池里,5分钟后,鸡蛋的温度为38℃,假定没有感到水变热,问鸡蛋达到20℃,还需多长时间?解:题意没有感到水变热,即池水中水温不变。

设:鸡蛋的温度为T,温度变化率就是dT/dt 其中t为时间,水的温度为T1,则鸡蛋与水温差为T-T1由题意有:T- T1=kdT/dt (其中k为比例常数) (1)方程(1)化为:dt=kdT/(T- T1)(2)对(2)两边同时积分之后并整理一下就得到:t=k*ln(T- T1)+C则k*ln(98-18)+ C=05=k*ln(38-18)+ct1=k*ln(20-18)+c-[k*ln(38-18)+c]=8.3(min)所以,还需8.3(min)。

2. 报童每天清晨从报社购进报纸零售,晚上将没有卖完的报纸退回。

设每份报纸的购进价为,零售价为,退回价为,应该自然地假设。

这就是说,报童售出一份报纸赚,退回一份报纸赔。

报童如果每天购进的报纸太少,不够卖的,会少赚钱;如果购进太多,卖不完,将要赔钱。

请你为报童筹划一下,他应该如何确定每天购进报纸的数量,以获得最大的收入。

解:设:报纸具有时效性每份报纸进价b元,卖出价a元,卖不完退回份报纸c元。

设每日的订购量为n,如果订购的多了,报纸剩下会造成浪费,甚至陪钱。

订的少了,报纸不够卖,又会少赚钱。

为了获得最大效益,现在要确定最优订购量n。

n的意义。

n是每天购进报纸的数量,确定n一方面可以使报童长期以内拥有一个稳定的收入,另一方面也可以让报社确定每日的印刷量,避免纸张浪费。

所以,笔者认为n的意义是双重的。

本题就是让我们根据a、b、c及r来确定每日进购数n。

基本假设1、假设报童现在要与报社签定一个长期的订购合同,所以要确定每日的订购量n。

数学建模全部作业

数学建模全部作业

一、图论(组合优化)和排列论实验解:设cij表示i年开始到j-1年结束购车的总消费,则有:C12=2.5+0.3-2.0=0.8,C13=2.5+0.3+0.5-1.6=1.7,C14=2.5+0.3+0.5+0.8-1.3=2 .8,C15=2.5+0.3+0.5+0.8+1.2-1.1=4.2,C23=2.6+0.3-2.0=0.9,C24=2.6+0.3+0. 5-1.6=1.8,C25=2.6+0.3+0.5+0.8-1.3=2.9,C34=2.8+0.3-2.0=1.1,C35=2.8+0.3 +0.5-1.6=2,C45=3.1+0.3-2.0=1.4;建模如下:sets:nodes/1..5/;arcs(nodes, nodes)|&1 #lt# &2: c, x;endsetsdata:c = 0.8 1.7 2.8 4.20.9 1.8 2.91.12.01.4;enddatan = @size(nodes);min = @sum(arcs: c * x);@for(nodes(i)| i #ne# 1 #and# i #ne# n:@sum(arcs(i,j): x(i,j)) = @sum(arcs(j,i): x(j,i)));@sum(arcs(i,j)| i #eq# 1 : x(i,j)) = 1;LINGO运行如下:Global optimal solution found.Objective value: 3.700000Total solver iterations: 0Variable Value Reduced CostX( 1, 2) 1.000000 0.000000X( 2, 5) 1.000000 0.000000由计算结果分析可得,其最短路径为1->2->5,最小花费为3.7万元。

即:该单位应该在第一年购买新设备,年末卖掉设备;第二年初更换新设备,一直用到第四年年末,再卖出。

数学建模作业及答案

数学建模作业及答案

数学建模作业姓名:叶勃学号:班级: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、马青公式 π=16arctan1/5-4arctan1/239 这个公式由英国天文学教授约翰·马青于 1706 年发现。他利用这个公式计
算到了 100 位的圆周率。马青公式每计算一项可以得到 1.4 位的十进制精度。因 为它的计算过程中被乘数和被除数都不大于长整数,所以可以很容易地在计算机 上编程实现。
练习题 6:兄弟三人戴帽子问题 解放前,在一个村子里住着聪明的三兄弟,他们除恶杀了财主的儿子,犯了人命案。县太爷有意想免他们
一死,决意出一个难题测测他们是否真的聪明,如果他们能在一个时辰内回答出来,就免他们一死,否则就被 处死。题目如下:
兄弟三人站成一路纵队(老三选择了站在最前面,他后面是老二,老大站在了最后面 ),并分别被蒙住了眼 睛,县太爷说我这里有两顶黑帽子和三顶红帽子,接着分别给他们头上各带了一顶帽子,然后又分别把被蒙住 的眼睛解开。
还有很多类似于马青公式的反正切公式。在所有这些公式中,马青公式似乎是 最快的了。虽然如此,如果要计算更多的位数,比如几千万位,马青公式就力不 从心了。
2、拉马努金公式 1914 年,印度天才数学家拉马努金在他的论文里发表了一系列共 14 条圆周
率的计算公式。这个公式每计算一项可以得到 8 位的十进制精度。1985 年 Gosper 用这个公式计算到了圆周率的 17,500,000 位。
此时,老大只可以看见老三和老二头上的帽子,老二只可以看见老三头上的帽子,老三看不见帽子。 只有一个时辰的时间,看谁能说出自己头上帽子的颜色,第一句声音有效。现在开始! (县太爷有多少种带帽子的方案,那一种最难?你能回答吗?)
解答:
县太爷一共有 7 种戴帽子方案:
1 黑黑红 2 黑红黑 3 黑红红 4 红红红 5 红红黑 6 红黑红 7 红黑黑
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数学建模作业题
习题1第4题. 根据表1.14的数据,完成下列数据拟合问题:
(1) 如果用指数增长模型0()0()e r t t x t x -=模拟美国人口从1790年至2000年的变化过程,请用MATLAB 统计工具箱的函数nlinfit 计算指数增长模型的以下三个数据拟合问题:
(i) 取定0x =3.9,0t =1790,拟合待定参数r ;
(ii) 取定0t =1790,拟合待定参数0x 和r ; (iii) 拟合待定参数0t 、0x 和r .
要求写出程序,给出拟合参数和误差平方和的计算结果,并展示误差平方和最小的拟合效果图.
(2) 通过变量替换,可以将属于非线性模型的指数增长模型转化成线性模型,并用MATLAB 函数polyfit 进行计算,请说明转化成线性模型的详细过程,然后写出程序,给出拟合参数和误差平方和的计算结果,并展示拟合效果图.
(3) 请分析指数增长模型非线性拟合和线性化拟合的结果有何区别?原因是什么?
(4) 如果用阻滞增长模型00
()
00()()e r t t Nx x t x N x --=+-模拟美国人口从1790年至2000年的
变化过程,请用MATLAB 统计工具箱的函数nlinfit 计算阻滞增长模型的以下三个数据拟合问题:
(i) 取定0x =3.9,0t =1790,拟合待定参数r 和N ; (ii) 取定0t =1790,拟合待定参数0x 、r 和N ;
(iii) 拟合待定参数0t 、0x 、r 和N .
要求写出程序,给出拟合参数和误差平方和的计算结果,并展示误差平方和最小的拟合效果图.
习题2第1题. 继续考虑第2.2节“汽车刹车距离”案例,请问“两秒准则”和“一车长度准则”一样吗?“两秒准则”是否足够安全?对于安全车距,你有没有更好的建议?
习题2第2题. 一盘录像带,从头转到尾,时间用了184分钟,录像机计数器读数从0000变到6061. 表2.5是观测得到的计数器读数,图2.7是录像机计数器工作原理示意图. 请问当计数器读数为4580时,剩下的一段录像带还能否录下一小时的节目?
习题3第4题. 某成功人士向学院捐献20万元设立优秀本科生奖学金,学院领导打算将这笔捐款以整存整取一年定期的形式存入银行,第二年一到期就支取,取出一部分作为当年的奖学金,剩下的继续以整存整取一年定期的形式存入银行……请你研究这个问题,并向学院领导写一份报告.
习题3第5题. 有一位老人60岁时将养老金10万元以整存零取方式(指本金一次存入,分次支取本金的一种储蓄)存入,从第一个月开始每月支取1000元,银行每月初按月利率0.3%把上月结余额孳生的利息自动存入养老金. 请你计算老人多少岁时将把养老金用完?如果想用到80岁,问60岁时应存入多少钱?
习题4第3题. 继续考虑第3.4.2小节“酵母培养物的增长”案例,建立微分方程模型,模拟酵母培养物的增长.
习题6第2题. 13名儿童参加了一项睡眠时间(分钟)与年龄(岁)关系的调查,表6.18中的睡眠时间是根据连续3天记录的每天睡眠时间的平均值得到的. 请建立和求解回归模型,解释得到的结果,给出10岁儿童的平均睡眠时间及预测区间.
习题6第3题. 水的沸点与大气压强有密切关系,表6.19中包含了17次试验中所测得的水的沸点(华氏温度)和大气压强(水银英寸),请建立回归模型估计沸点和压强之间的关系,并给出当沸点为201.5F 时压强的预测值及预测区间.
习题7第2题. 某配件厂为装配线生产若干种部件. 每次轮换生产不同的部件时,因更换设备要付生产准备费(与生产数量无关). 同一部件的产量大于需求时,因积压资金、占用仓库要付库存费. 今已知某一部件的日需求量100件,生产准备费5000元,库存费每日每件1元. 如果生产能力远大于需求,并且不允许出现缺货,请制定最优生产计划.
习题7第3题. 某商场把销售所剩的空纸皮箱压缩并打成包准备回收,每天能产生5包,在商场后院存放的费用是每包每天10元. 另一家公司负责将这些纸包运送到回收站,要收取固定费用1000元租装卸车,外加运输费每包100元. 请制定运送纸包到回收站的最优策略.。

相关文档
最新文档