数学建模作业——实验1

合集下载

数学建模实验报告

数学建模实验报告

湖南城市学院数学与计算科学学院《数学建模》实验报告专业:学号:姓名:指导教师:成绩:年月日目录实验一 初等模型........................................................................ 错误!未定义书签。

实验二 优化模型........................................................................ 错误!未定义书签。

实验三 微分方程模型................................................................ 错误!未定义书签。

实验四 稳定性模型.................................................................... 错误!未定义书签。

实验五 差分方程模型................................................................ 错误!未定义书签。

实验六 离散模型........................................................................ 错误!未定义书签。

实验七 数据处理........................................................................ 错误!未定义书签。

实验八 回归分析模型................................................................ 错误!未定义书签。

实验一 初等模型实验目的:掌握数学建模的基本步骤,会用初等数学知识分析和解决实际问题。

实验内容:A 、B 两题选作一题,撰写实验报告,包括问题分析、模型假设、模型构建、模型求解和结果分析与解释五个步骤。

数学建模实验报告

数学建模实验报告

数学建模实验报告1.流⽔问题问题描述:⼀如下图所⽰的容器装满⽔,上底⾯半径为r=1m,⾼度为H=5m,在下地⾯有⼀⾯积为B0.001m2的⼩圆孔,现在让⽔从⼩孔流出,问⽔什么时候能流完?解题分析:这个问题我们可以采⽤计算机模拟,⼩孔处的⽔流速度为V=sqrt[2*g*h],单位时间从⼩孔流出的⽔的体积为V*B,再根据⼏何关系,求出⽔⾯的⾼度H,时间按每秒步进,记录点(H,t)并画出过⽔⾯⾼度随时间的变化图,当⽔⾯⾼度⼩于0.001m 时,可以近似认为⽔流完了。

程序代码:Methamatic程序代码:运⾏结果:(5)结果分析:计算机仿真可以很直观的表现出所求量之间的关系,从图中我们可以很⽅便的求出要求的值。

但在实际编写程序中,由于是初次接触methamatic 语⾔,对其并不是很熟悉,加上个⼈能⼒有限,所以结果可能不太精确,还请见谅。

2.库存问题问题描述某企业对于某种材料的⽉需求量为随机变量,具有如下表概率分布:每次订货费为500元,每⽉每吨保管费为50元,每⽉每吨货物缺货费为1500元,每吨材料的购价为1000元。

该企业欲采⽤周期性盘点的),(S s 策略来控制库存量,求最佳的s ,S 值。

(注:),(S s 策略指的是若发现存货量少于s 时⽴即订货,将存货补充到S ,使得经济效益最佳。

)问题分析:⽤10000个⽉进⾏模拟,随机产⽣每个⽉需求量的概率,利⽤计算机编程,将各种S 和s 的取值都遍历⼀遍,把每种S,s的组合对应的每⽉花费保存在数组cost数组⾥,并计算出平均⽉花费average,并⽤类answer来记录,最终求出对应的S和s。

程序代码:C++程序代码:#include#include#include#include#define Monthnumber 10000int Need(float x){int ned = 0;//求每个⽉的需求量if(x < 0.05)ned = 50;else if(x < 0.15)ned = 60;else if(x < 0.30)ned = 70;else if(x < 0.55)ned = 80;else if(x < 0.75)ned = 90;else if(x < 0.85)ned = 100;else if(x < 0.95)ned = 110;else ned = 120;return ned;}class A{public:int pS;int ps;float aver;};int main(){A answer;answer.aver=10000000;//int cost[Monthnumber+1]={0}; float average=0;int i;float x;int store[Monthnumber];//srand((int)time(0));for(int n=6;n<=12;n++){// int n=11;int S=10*n;for(int k=5;k{// int k=5;int s=k*10;average=0;int cost[Monthnumber+1]={0};for(i=1;i<=Monthnumber;i++){store[i-1]=S;srand(time(0));x=(float)rand()/RAND_MAX; //产⽣随机数//cout<<" "<//cout<int need=Need(x);if(need>=store[i-1]){cost[i]= 1000*S + (need - store[i-1])*1500 + 500;store[i]=S;}else if(need>=store[i-1]-s){cost[i]=1000*(need+S-store[i-1]) + 50*(store[i-1]-need) + 500; store[i]=S;}else{cost[i]=(store[i-1]-need)*50;store[i]=store[i-1]-need;}average=cost[i]+average;}average=average/Monthnumber;cout<<"n="<cout<<"花费最少时s应该为:"<cout<<"平均每⽉最少花费为:"<}运⾏结果:结果分析:⽤计算机模拟的结果和⽤数学分析的结果有⼀定的差异,由于计算机模拟时采⽤的是随机模型⽽我⽤time函数和rand函数产⽣真随机数,所以在每次的结果上会有所差异,但对于⼀般的⽣产要求亦可以满。

数学建模活动(1)-练习题

数学建模活动(1)-练习题
数学建模活动(1)
课后练习
请同学们仿照上述过程开展一次建立模型解决实际问题的活动,可以继续研究不同室温下泡制一杯最佳口感茶水所需的时间,也可以从下列选题中选择一个:
1、应在炒菜之前多长时间将冰箱里的肉拿出来解冻?
2、根据某一同学的身高和体重,判断该同学是否超重.
3、用微波炉或电磁炉烧一壶开水,找到最省电的功率设定方法.
4估计阅读一本书所需要的时间.
也可以根据自己的兴趣成3—5人的研究小组,每位同学参加其中一个小组。在小组内,要确定一个课题负责人,使每位成员都有明确的分工.拟定研究课题、确定研究方案、规划研究步骤、编制研究手册,然后在班里进行一次开题报告。

数学建模实验报告

数学建模实验报告

《数学建模实验》实验报告学院名称数学与信息学院专业名称提交日期课程教师实验一:数学规划模型AMPL求解实验内容1. 用AMPL求解下列问题并作灵敏度分析:一奶制品加工厂用牛奶生产A1和A2两种奶制品,1桶牛奶可以在甲类设备上用12小时加工成3公斤A1或者在乙类设备上用8小时加工成4公斤A2,且都能全部售出,且每公斤A1获利24元,每公斤A2获利16元。

先加工厂每天能得到50桶牛奶的供应,每天工人总的劳动时间为480小时,并且甲类设备每天至多加工100公斤A1,乙类设备的加工能力没有限制,试为该厂制定一个计划,使每天的获利最大。

(1)建立模型文件:milk.modset Products ordered;param Time{i in Products }>0;param Quan{i in Products}>0;param Profit{i in Products}>0;var x{i in Products}>=0;maximize profit: sum{i in Products} Profit [i]* Quan [i]*x[i];subject to raw: sum{i in Products}x[i] <=50;subject to time:sum{i in Products}Time[i]*x[i]<=480;subject to capacity: Quan[first(Products)]*x[first(Products)]<=100;(2)建立数据文件milk.datset Products:=A1 A2;param Time:=A1 12 A2 8;param Quan:=A1 3 A2 4;param Profit:=A1 24 A2 16;(3) 建立批处理文件milk.runmodel milk.mod;data milk.dat;option solver cplex;solve;display x;(4)运行运行结果:CPLEX 11.0.0: optimal solution; objective 33602 dual simplex iterations (1 in phase I)x [*] :=A1 20A2 30;(5)灵敏度分析:model milk.mod;data milk.dat;option solver cplex;option cplex_options 'sensitivity';solve;display x;display x.rc, x.down, x.up;display raw, time, capacity;display raw.down, raw.up,raw.current, raw.slack;得到结果:【灵敏度分析】: x.rc x.down x.up:=A1 -3.55271e-15 64 96A2 0 48 72;raw = 48time = 2capacity = 0raw.down = 43.3333raw.up = 60raw.current = 50raw.slack = 0某公司有6个建筑工地,位置坐标为(a i, b i)(单位:公里),水泥日用量d i (单位:吨)1) 现有j j j吨,制定每天的供应计划,即从A, B两料场分别向各工地运送多少吨水泥,使总的吨公里数最小。

数学建模实验报告

数学建模实验报告

数学建模实验报告一、实验目的1.通过具体的题目实例, 使学生理解数学建模的基本思想和方法, 掌握数学建模分析和解决的基本过程。

2、培养学生主动探索、努力进取的的学风, 增强学生的应用意识和创新能力, 为今后从事科研工作打下初步的基础。

二、实验题目(一)题目一1.题目: 电梯问题有r个人在一楼进入电梯, 楼上有n层。

设每个乘客在任何一层楼出电梯的概率相同, 试建立一个概率模型, 求直到电梯中的乘客下完时, 电梯需停次数的数学期望。

2.问题分析(1)由于每位乘客在任何一层楼出电梯的概率相同, 且各种可能的情况众多且复杂, 难于推导。

所以选择采用计算机模拟的方法, 求得近似结果。

(2)通过增加试验次数, 使近似解越来越接近真实情况。

3.模型建立建立一个n*r的二维随机矩阵, 该矩阵每列元素中只有一个为1, 其余都为0, 这代表每个乘客在对应的楼层下电梯(因为每个乘客只会在某一层下, 故没列只有一个1)。

而每行中1的个数代表在该楼层下的乘客的人数。

再建立一个有n个元素的一位数组, 数组中只有0和1,其中1代表该层有人下, 0代表该层没人下。

例如:给定n=8;r=6(楼8层, 乘了6个人),则建立的二维随机矩阵及与之相关的应建立的一维数组为:m =0 0 1 0 0 01 0 0 0 0 00 0 0 0 0 00 1 0 0 0 00 0 0 0 0 00 0 0 0 0 10 0 0 0 1 00 0 0 1 0 0c = 1 1 0 1 0 1 1 14.解决方法(MATLAB程序代码):n=10;r=10;d=1000;a=0;for l=1:dm=full(sparse(randint(1,r,[1,n]),1:r,1,n,r));c=zeros(n,1);for i=1:nfor j=1:rif m(i,j)==1c(j)=1;break;endcontinue;endends=0;for x=1:nif c(x)==1s=s+1;endcontinue;enda=a+s;enda/d5.实验结果ans = 6.5150 那么, 当楼高11层, 乘坐10人时, 电梯需停次数的数学期望为6.5150。

数学建模 -实验报告1

数学建模 -实验报告1
推导出了动力学方程
������������⁄������������ = ������������(1 − (������ + ������)) − ������1������∗������,
(4 − 3)
������������∗⁄������������ = −������1������∗������ + ������2������
二、 问题分析
建立肿瘤细胞增长模型时,我们可以从自由增长模型开始分析,引进 Logistic 阻滞增长模型,构成肿瘤细胞增长初步框架。再者肿瘤细胞不同于普 通细胞,其生长受到人体自身免疫系统的制约。于是综合考虑正常细胞转化,癌 细胞增殖,癌细胞死亡,癌细胞被效应细胞消除等情况,建立动力学方程。并对 模型进行适当简化求解。在放射治疗方案的设计中,我们可以引入放射生物学中 广泛接受的 LQ 模型对问题进行分析,由于放疗对人体伤害相当大,因此我们采 取分次逐次放疗的方式进行治疗。我们具体分两种情形进行讨论,一是在总剂量 一定的条件下,不同的分次剂量组合对生物效应的影响;二是在产生相同生物效 应的情况下,分析最优的分次剂量组合。
易算出癌细胞转入活动期已有 300 多天,故如何在早期发现癌症是攻克癌症的关键之一 (2)手术治疗常不能割去所有癌细胞,故有时需进行放射疗法。射线强度太小无法杀
死癌细胞,太强病人身体又吃不消且会使病人免疫功能下降。一次照射不可能杀死全部癌细 胞,请设计一个可行的治疗方案(医生认为当体内癌细胞数小于 100000 个时即可凭借体内 免疫系统杀灭)。
进一步简化,根据(4-4),(4-5)式可知,效应细胞������∗和复合物������有出有进.假 设出入保持平衡,则有
������ + ������∗ = C (C 为常数)

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

数学建模第一次实验报告

数学建模第一次实验报告

数学建模第一次实验报告问题描述:2.某动物从食物中每天得到2500卡(1卡=4.18焦)的热量,其中1200卡用于基本的新陈代谢,每天每kg的体重要再消耗16卡,假如它每增加1kg体重需要10000卡的热量,问该动物体重将怎样变化?解:设动物体重为m。

令动物每日消耗热量等于获取的热量,可求得最大体重。

此时,2500=1200+16mm=81.25kg根据生物学知识可知,没有动物的出生时的体重会大于成年后的体重,即m≤81.25kg。

又设每天体重的变化量为dm,2500=1200+16m+10000dm/dtt=625In(16m-1300)m=1/16*(e^(t/625)+1300)3.37支球队进行冠军争夺赛,每轮比赛中出场的两只球队中的胜者及轮空着进入下一轮,直至比赛结束,问共需进行多少场比赛?解:各球队的胜负对比赛场次无影响,忽略。

将其拆解为一轮一轮的比赛分析:第一轮:37÷2=18……1——>19第二轮:19÷2=9 ……1——>10第三轮:10÷2=5 ——>5第四轮: 5÷2=2 ……1——>3第五轮: 3÷2=1 ……1——>2第六轮: 2÷2=1 ——>1所以,一共进行六轮比赛,其场数为:18+9+5+2+1+1=36(场)答:一共36场比赛。

4.1条河宽1km,两岸各有一个城镇A与B,A与B的直线距离为4km。

今需铺设一条电缆连接A与B,已知地下电缆的修建费是2万元/km,水下电缆的修建费是4万元/km,假设两岸为平行的直线,问应该如何架设电缆方可以使总建设费用最少?AB D C如图,设C 点。

易知A-D-B 为铺设路线,设AD 长为a ,BD 长为b ,总花费为m 。

其中1<a<4,0<b<15所以,m=4a+2b (1)(15-b )2+12=a 2 (2)所以, m=4215216b b +-+2b 求得m 最小值即可。

数学建模作业一:汽车刹车距离

数学建模作业一:汽车刹车距离

汽车刹车距离一、问题描写司机在碰到突发紧迫情形时都邑刹车,从司机决议刹车开端到汽车停滞行驶的距离为刹车距离,车速越快,刹车距离越长.那么刹车距离与车速之间具有什么样的关系呢?二、问题剖析汽车的刹车距离有反响距离和刹车距离两部分构成,反响距离指的是司机看到须要刹车的情形到汽车制动器开端起感化汽车行使的距离,刹车距离指的是制动器开端起感化到汽车完整停滞的距离.反响距离有反响时光和车速决议,反响时光取决于司机小我状态(敏锐.机灵等)和制动体系的敏锐性,因为很难对反响时光进行差别,是以,平日以为反响时光为常数,并且在这段时光内车速不变.刹车距离与制动感化力.车重.车速以及路面状态等身分有关系.由能量守恒制动力所做的功等于汽车动能的转变.设计制动器的一个合理原则是,最大制动力大体上与汽车的质量成正比,汽车的减速度根本上是常数.路面状态可以为是固定的.三、问题求解1、模子假设依据上述剖析,可作如下假设:①刹车距离d等于反响距离1d和制动距离2d之和;②反响距离1d 与车速v 成正比,且比例系数为反响时光t;③刹车时应用最大制动力F,F 作的功等于汽车动能的转变,且F 与车质量m 成正比;④人的反响时光t 为一个常数; ⑤在反响时光内车速v 不变 ; ⑥路面状态是固定的;⑦汽车的减速度a 根本上是一个常数. 2、 模子树立由上述假设,可得: ⑴tv d =2;⑵2221mv Fd =,而ma F =,则2221v ad =.所以22kv d =. 综上,刹车距离的模子为2kv tv d +=. 3.参数估量可用我国某机构供给的刹车距离现实不雅察数据来拟合未知参数t 和k.转化单位后得:车速(公里/小时) 20 40 60 80 100 120 140现实刹车距离(米)118.0用Mathematica 进行拟合,代码如下: Clear[x,v,d];x={{20/3.6,6.5},{40/3.6,17.8},{60/3.6,33.6},{80/3.6,57.1},{100/3.6,83.4},{120/3.6,118},{140/3.6,153.5}}; d=Fit[x,{v,v^2},v];Print["d=",d];Plot[d,{v,0,200/3.6}] 成果: 4. 成果剖析将拟合成果与现实成果比较:(代码) Clear[v,d];d=0.65218*v/3.6+0.0852792*(v/3.6)^2;For[v=20,v<=140,v=v+20,Print["速度为",v,"km/h 时刹车距离为",d]] 成果:车速(公里/小时) 20 40 60 80 100 120 140 现实刹车距离(米) 盘算刹车距离(米)盘算刹车距离与现实刹车距离基底细当.综上,反响时光t 约等于0.6522秒,刹车时减速度约等于2/62/1s m k ≈.刹车距离与车速的关系知足:208528.06522.0d v v +=.。

1数学建模实验-圆周率的计算

1数学建模实验-圆周率的计算

19-23
设计方案
在正方形 0< x <1, 0< y<1 上随机的投大量的点,那么 落在四分之一园内的点数 数m与在正方形内的点数n 之比m/n应为这两部分图形 面积之比=π/4,故 π=4 m/n 计算机模拟:产生区间[0,1]上数目为n的一组 随机数(x,y),计算满足x2+y2<1的点数m
Ä 1630年,最后一位用古典方法求π的人 年 最后一位用古典方法求π 格林伯格也只求到了π的第 位小数 格林伯格也只求到了 的第39位小数 的第
8-23
分析方法
从十七世纪中叶起, 从十七世纪中叶起,人们开始用更先进的 分析方法来求π的近似值,其中应用的主 分析方法来求π的近似值, 要工具是收敛的无穷乘积和无穷级数, 要工具是收敛的无穷乘积和无穷级数,在 本节中我们将介绍一些用此类方法求π 本节中我们将介绍一些用此类方法求π近 似值的实例。 似值的实例。
>>ans=
3.1415925796063512110
>> vpa(4*symsum((-1)^(k-1)/(2*k-1)*(1/2)^(2*k1),1,20)+4*symsum((-1)^(k-1)/(2*k-1)*(1/3)^(2*k-1),1,20),20)
>>ans=
3.1415926535897574098
16-23
方法2 利用数值积分
1 设 y(x) = 1+ x2
1 A = 4∫ dx = π 2 01 x +
1
将[0,1]区间 n 等分,取 xk=k/n, yk= 1/ (1+ k2) (1+x
2 梯形法⇒ A = [2( y1 + y2 +L+ yn−1) + y0 + yn ] n

线性代数数学建模案例1

线性代数数学建模案例1
网络分析要解决的问题是:在部分信息(如 网络的输入量)已知的情况下,确定每一分支中 的流量。
案例1 交通网络流量分析问题
城市道路网中每条道路、每个交叉 路口的车流量调查,是分析、评价及改 善城市交通状况的基础。根据实际车流 量信息可以设计流量控制方案,必要时 设置单行线,以免大量车辆长时间拥堵。
下图为某城市的局部单行示意图
【模型假设】假设不考虑价格变动等其他因素.
【模型建立】设煤矿, 电厂, 铁路分别产出x元, y元, z元刚好满足需求. 则有下表
产出(1元)
产出



煤0
0.6 0.5
x
分配 0.6y + 0.5z
订单 60000
消 电 0.3 0.1 0.1
y

0.3x + 0.1y + 0.1z 100000
几条道路的流量统计? (3) 当x4 = 350时, 确定x1, x2, x3的值. (4) 若x4 = 200, 则单行线应该如何改动才合
理? 。
【模型假设】: (1) 每条道路都是单行线 (2) 每个交叉路口进入和离开的车辆数目相等.
【模型建立】 根据图3和上述假设, 在①, ②, ③, ④ 四个路口进出车辆数目分别满足:
【模型分析】
(1) 由(A, b)的行最简形可见, 上述方程组中的最
后一个方程是多余的. 这意味着最后一个方程中的
数据“300”x可1 以x4不1用00统计.
(2)由

x2

x4

600
可得
x3 x4 300
x2 x1 500

x3

x1

200

数学建模实验报告1

数学建模实验报告1

数学建模实验报告1桂林电⼦科技⼤学2017-2018学年第1学期数学建模⼀、实验⽬的1. 熟悉MATLAB 软件的⽤户环境;2. 了解MATLAB 软件的⼀般命令;3. 掌握MATLAB 向量、数组、矩阵操作与运算函数;4. 掌握MATLAB 软件的基本绘图命令;5. 掌握MATLAB 语⾔的⼏种循环、条件和开关选择结构及其编程规范。

⼆、实验内容1. MATLAB 软件的矩阵输⼊和操作2. ⽤MA TLAB 语⾔编写命令M ⽂件和函数M ⽂件3. 直接使⽤MATLAB 软件进⾏作图练习;三、实验任务1. 有⼀个4×5的矩阵,编程求出其元素最⼤值及其所在的位置。

Jm.m ⽂件代码: clear;a=input('请输⼊⼀个4*5矩阵'); max=a(1,1); maxi=0; maxj=0; for i=1:4 for j=1:5if a(i,j)>max max=a(i,j); maxi=i; maxj=j;end end endfprintf('最⼤值为:%d 位置:o%d %d \n',max,maxi,maxj); 实验结果:2. 有⼀函数f(x,y)=x 2+sin xy+2y,写⼀程序,输⼊⾃变量的值,输出函数值。

Jm_5.m ⽂件代码: function f=Jm_5(x,y) f=x.^2+sin(x*y)+2*y;实验结果:3.⽤surf,mesh绘制曲⾯z=2x2+y2。

Jm5.m代码:x=-3:0.1:3;y=1:0.1:5;[X,Y]=meshgrid(x,y);Z=2*X.^2+Y.^2;subplot(1,2,1);surf(X,Y,Z);title('surf(x,y)');subplot(1,2,2);mesh(X,Y,Z);title('mesh(x,y)');实验结果:4.在同⼀平⾯的两个窗⼝中分别画出⼼形线和马鞍⾯。

数学建模作业实验1数学建模入门

数学建模作业实验1数学建模入门

数学建模作业(实验1数学建模入门)基本实验1.贷款问题小王夫妇计划贷款20万元购买一套房子,他们打算用20年的时间还清贷款。

目前,银行的利率是0.6%/月。

他们采用等额还款的方式(即每月的还款额相同)偿还贷款。

(1)在上述条件下,小王夫妇每月的还款额是多少?共计付了多少利息?(2)在贷款满5年后,他们认为他们有经济能力还完余下的款额,打算提前还贷,那么他们在第6年初,应一次付给银行多少钱,才能将余下全部的贷款还清?(3)如果在第6年初,银行的贷款利率由0.6%/月调到0.8%/月,他们仍然采用等额还款的方式,在余下的15年内将贷款还清,那么在第6年后,每月的还款额应是多少?(4)某借贷公司的广告称,对于贷款期在20年以上的客户,他们帮你提前三年还清贷款。

但条件是:(i)每半个月付款一次,但付款额不增加,即一次付款额是原付给银行还款额的1/2;(ii)因为增加必要的档案、文书等管理工作,因此要预付给借贷公司贷款总额10%的佣金。

试分析,小王夫妇是否要请这家借贷公司帮助还款。

解答解:设A k为第k个月的欠款额,r为月利率,x为每个月的还款额,则A0=200000,且第k个月的欠款额=第k-1个月的欠款额 月利率+第k-1个月的欠款额-每月的还款额,即:A k=A k-1(1+r)-x,k=1,2,…,N,由以上差分方程可得:第1个月的欠款额:A 1=A 0(1+r )-x第2个月的欠款额:A 2=A 1(1+r )-x=A 0(1+r)2-x[(1+r)+1]……第k 个月的欠款额:A k =A k -1(1+r )-x=A 0(1+r )k -x[(1+r)k-1+…+(1+r)+1] =0(1)1(1)(1)1k kr A r x r +-+-+- 贷款总月数为n ,也就是说,第n 个月的欠款额为0,即A n =0,令n=k ,导出00(1)(1)1nn A r r x A r r +=>⋅+-,说明每个月的还款额一定大于贷款额×月利率。

数学建模实验项目

数学建模实验项目

数学建模实验项⽬数学建模实验指导书数学建模实验项⽬⼀养⽼基⾦问题⼀、实验⽬的与意义:1、练习初等问题的建模过程;2、练习Matlab基本编程命令;⼆、实验要求:3、较能熟练应⽤Matlab基本命令和函数;4、注重问题分析与模型建⽴,了解建模⼩论⽂的写作过程;5、提⾼Matlab的编程应⽤技能。

三、实验学时数:2学时四、实验类别:综合性五、实验内容与步骤:(1.必做,2、3选⼀)1.某⼤学青年教师从31岁开始建⽴⾃⼰的养⽼基⾦,他把已有的积蓄10000元也⼀次性地存⼊,已知⽉利率为0.001(以复利计),每⽉存⼊700元,试问当他60岁退休时,他的退休基⾦有多少?⼜若,他退休后每⽉要从银⾏提取1000元,试问多少年后他的基⾦将⽤完?2.贷款助学问题。

3贷款购房问题。

⾃⼰调查设计具体情况数学建模实验项⽬⼆梯⼦问题⼀、实验⽬的与意义:1、进⼀步熟悉数学建模步骤;2、练习Matlab优化⼯具箱函数;3、进⼀步熟悉最优化模型的求解过程。

⼆、实验要求:1、较能熟练应⽤Matlab⼯具箱去求解常规的最优化模型;2、注重问题分析与模型建⽴,熟悉建模⼩论⽂的写作过程;3、提⾼Matlab的编程应⽤技能。

三、实验学时数:2学时四、实验类别:综合性五、实验内容与步骤:⼀幢楼房的后⾯是⼀个很⼤的花园。

在花园中紧靠着楼房建有⼀个温室,温室⾼10英尺,延伸进花园7英尺。

清洁⼯要打扫温室上⽅的楼房的窗户。

他只有借助于梯⼦,⼀头放在花园中,⼀头靠在楼房的墙上,攀援上去进⾏⼯作。

他只有⼀架20⽶长的梯⼦,你认为他能否成功?能满⾜要求的梯⼦的最⼩长度是多少?步骤:1.先进⾏问题分析,明确问题;2.建⽴模型,并运⽤Matlab函数求解;3.对结果进⾏分析说明;4.设计程序画出图形,对问题进⾏直观的分析和了解(主要⽤画线函数plot,line)5.写⼀篇建模⼩论⽂。

数学建模实验项⽬三确定肥猪的最佳销售时机⼀、实验⽬的与意义:1、认识微分法的建模过程;2、认识微分⽅程的数值解法。

数学建模作业一:汽车刹车距离

数学建模作业一:汽车刹车距离

汽车刹车距离一、 问题描述司机在遇到突发紧急情况时都会刹车,从司机决定刹车开始到汽车停止行驶的距离为刹车距离,车速越快,刹车距离越长。

那么刹车距离与车速之间具有什么样的关系呢?二、 问题分析汽车的刹车距离有反应距离和刹车距离两部分组成,反应距离指的是司机看到需要刹车的情况到汽车制动器开始起作用汽车行使的距离,刹车距离指的是制动器开始起作用到汽车完全停止的距离。

反应距离有反应时间和车速决定,反应时间取决于司机个人状况(灵敏、机警等)和制动系统的灵敏性,由于很难对反应时间进行区别,因此,通常认为反应时间为常数,而且在这段时间内车速不变。

刹车距离与制动作用力、车重、车速以及路面状况等因素有关系。

由能量守恒制动力所做的功等于汽车动能的改变。

设计制动器的一个合理原则是,最大制动力大体上与汽车的质量成正比,汽车的减速度基本上是常数。

路面状况可认为是固定的。

三、 问题求解1、 模型假设根据上述分析,可作如下假设:①刹车距离d 等于反应距离1d 和制动距离2d 之和;②反应距离1d 与车速v 成正比,且比例系数为反应时间t ;③刹车时使用最大制动力F ,F 作的功等于汽车动能的改变,且F 与车质量m 成正比; ④人的反应时间t 为一个常数;⑤在反应时间内车速v 不变 ;⑥路面状况是固定的;⑦汽车的减速度a 基本上是一个常数。

2、 模型建立由上述假设,可得:⑴tv d =2; ⑵2221mv Fd =,而ma F =,则2221v ad =。

所以22kv d =。

综上,刹车距离的模型为2kv tv d +=。

3、 参数估计可用我国某机构提供的刹车距离实际观察数据来拟合未知参数t 和k 。

转化单位后得:车速(公里/小时)20 40 60 80 100 120 140实际刹车距离(米) 6.5 17.8 33.6 57.1 83.4 118.0 153.5用Mathematica进行拟合,代码如下:Clear[x,v,d];x={{20/3.6,6.5},{40/3.6,17.8},{60/3.6,33.6},{80/3.6,57.1},{100/3.6,83.4},{120/ 3.6,118},{140/3.6,153.5}};d=Fit[x,{v,v^2},v];Print["d=",d];Plot[d,{v,0,200/3.6}]结果:4、结果分析将拟合结果与实际结果对比:(代码)Clear[v,d];d=0.65218*v/3.6+0.0852792*(v/3.6)^2;For[v=20,v<=140,v=v+20,Print["速度为",v,"km/h时刹车距离为",d]]结果:车速(公里/小时)20 40 60 80 100 120 140实际刹车距离(米) 6.5 17.8 33.6 57.1 83.4 118.0 153.5计算刹车距离(米) 6.2 17.8 34.6 56.6 83.9 116.5 154.3计算刹车距离与实际刹车距离基本相当。

(完整版)数学建模作业1长方形椅子能在不平的地面上放稳吗?

(完整版)数学建模作业1长方形椅子能在不平的地面上放稳吗?

注意到椅脚连线呈长方形,长方形是中心对称图形,绕它的对称中心旋转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,π],从而将原问题数学化。

数学模型:已知f(θ)和g(θ)是θ的非负连续函数,对任意θ,f(θ)•g(θ)=0,证明:存在θ0∈[0,π],使得f(θ0)=g(θ0)=0成立。

五、模型求解(显示模型的求解方法、步骤及运算程序、结果)如果f(0)=g(0)=0,那么结论成立。

如果f(0)与g(0)不同时为零,不妨设f(0)>0,g(0)=0。

此时,将长方形ABCD。

数学建模实验报告1,层次分析法

数学建模实验报告1,层次分析法

数学建模实验报告一、实验要求柴静的纪录片《穹顶之下》从独立媒体人的角度调查了席卷全国多个省份的雾霾的成因,提出解决的方法有:关停重污染的钢铁厂、提高汽柴油品质、淘汰排放不达标汽车、提高洗煤率等,请仔细观看该纪录片,根据雾霾的成因,选择你认为治理雾霾确实可行的几个方案,并用AHP方法给出这几个主要方案的重要性排序。

二、前期准备1、理解层次分析法(AHP)的原理、作用,掌握其使用方法。

2、观看两遍柴静所拍摄的纪录片《穹顶之下》,选出我认为可较为有效地治理雾霾的几个方法,初步确定各方法的有效性(即权重)。

3、初步拟定三个方案,每个方案中各个治理方法的权重不同。

三、思路&分析1、根据纪录片《穹顶之下》和个人的经验判断给出各个记录雾霾的方法对于治理雾霾的判断矩阵,以及三个不同方案对于五大措施的判断矩阵。

2、了解了AHP的原理后,不难发现MATLAB在其中的作用主要是将判断矩阵转化为因素的权重矩阵。

当然矩阵要通过一致性检验,得到的权重才足够可靠。

3、分别得到准则层对目标层、方案层对准则层的权重之后,进行层次总排序及一致性检验。

得到组合权向量(方案层对目标层)即可确定适用方案。

四、实验过程1、确定层次结构2、构造判断矩阵(1)五大措施对于治理雾霾(准则层对目标层)的判断矩阵(2)三个方案对于五大措施(方案层对准则层)的判断矩阵3、层次单排序及一致性检验该部分在MATLAB中实现,每次进行一致性检验和权向量计算时,步骤相同,输入、输出参数一致。

(虽然输入的矩阵阶数可能不同,但可以不把矩阵阶数作为参数输入,而通过[n,n]=size(A)来算得阶数。

)因此考虑将这个部分定义为一个函数judge,输入一个矩阵A,打印一致性检验结果和权向量计算结果,并返回权向量、一致性指标CI、平均随机一致性指标RI。

将此脚本存为judge.m,在另一脚本ahp.m中调用。

代码如下:调试通过后,下面便用此函数进行一致性检验及权向量计算:(1)准则层对目标层(A矩阵)(2)方案层对准则层(BB矩阵)代码:结果:注:实际实验时,一开始构造的五个矩阵中有两个没有通过一致性检验。

数学建模作业(一)1

数学建模作业(一)1

第一题: 某班共45人,要去离校7.7千米的风景区旅游。

学校派了一辆可坐12人的校车接送。

为了尽快又同时到达目的地,校车分段分批接送学生。

已知校车速度为每小时70千米,学生步行的速度为每小时5千米。

如果上午七点出发,问最快什么时候全班同时到达目的地?(班长作为联系人要始终跟车)
第二题:某人为了锻炼身体,每天早晨坚持晨跑30分钟, 其中从A到B为800米上坡路,从B到C为1000米平路。

问在30分钟内跑完1800米,怎样安排跑步计划,才能使锻炼效果最佳?(即总疲劳程度伟为最低)
第三题:一辆小汽车与一辆大卡车在一段狭路上相遇,只有倒车才能继续通行。

如果小汽车的速度为大卡车的3倍,两车倒车的速度是各自正常速度的1/5,在这段狭路上,小汽车需倒车的路程是大卡车需倒车路程的4倍。

那么,为了使后通过狭路的那辆车尽早地通过这段狭路,问怎样倒车较为合理?
第四题:某人在一家公司工作,目前年薪为1万元。

老板说,现在有两种方案可供选择:第一种,每一年加1000元;第二种,每半年加300元。

试问:
(1)如果你在该公司工作5年,用哪一种方案收入高?
(2)如果你在该公司工作5年,将第二种方案中的每半年加300元改为a元时,那一种方案收入高?
(3)如果你在该公司工作n年,用哪一种方案收入高?
第五题:一个直角走廊宽为1.5米,有一辆转动灵活的平板水平推车,宽为1米,长为2.2米,问能否将其推过直角走廊?说明理由。

相关主题
  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

数学建模作业——实验1学院:软件学院姓名:学号:班级:软件工程2015级 GCT班邮箱:电话:日期:2016年5月10日基本实验1.椅子放平问题依照1.2.1节中的“椅子问题”的方法,将假设中的“四腿长相同并且四脚连线呈正方形”,改为“四腿长相同并且四脚连线呈长方形”,其余假设不变,问椅子还能放平吗?如果能,请证明;如果不能,请举出相应的例子。

答:能放平,证明如下:如上图,以椅子的中心点建立坐标,O为原点,A、B、C、D为椅子四脚的初始位置,通过旋转椅子到A’、B’、C’、D’,旋转的角度为α,记A、B两脚,C、D两脚距离地面的距离为f(α)和g(α),由于椅子的四脚在任何位置至少有3脚着地,且f(α)、g(α)是α的连续函数,则f(α)和g(α)至少有一个的值为0,即f(α)g(α)=0,f(α)≥ 0,g(α)≥0,若f(0)>0,g(0)=0,则一定存在α’∈(0,π),使得f(α’)=g(α’)=0令α=π(即椅子旋转180°,AB 边与CD 边互换),则f(π)=0,g(π)>0定义h(α)= f(α)-g(α),得到h(0)=f(0)-g(0)>0h(π)=f(π)-g(π) <0根据连续函数的零点定理,则存在α’∈( 0,π),使得h(α’)= f(α’)-g(α’)=0结合条件f(α’)g(α’)=0,从而得到f(α’)=g(α’)=0,即四脚着地,椅子放平。

2. 过河问题依照1.2.2节中的“商人安全过河”的方法,完成下面的智力游戏:人带着猫、鸡、米过河,船除需要人划之外,至多能载猫、鸡、米之一,而当人不在场时,猫要吃鸡、鸡要吃米,试设计一个安全过河的方案,并使渡河的次数尽量的少。

答: 用i =1,2,3,4分别代表人,猫,鸡,米。

1=i x 在此岸,0=i x 在对岸,()4321,,,x x x x s =此岸状态,()43211,1,1,1x x x x D ----=对岸状态。

安全状态集合为 :乘船方案,记作()4321,,,u u u u U =,当i 在船上时记1=i u ,否则记0=i u ,允许决策集合为()()()(){}0,0,0,1,1,0,0,1,0,1,0,1,0,0,1,1=U因为乘船k 为奇数时船从此案驶向彼岸,k 为偶数时船由彼岸驶向此岸,所以状态k s 随决策k U 变化的规律为()kU kk s k s 11-+=+, 设计安全过河方案归结为求决策,,,,21U U U U n ∈ ,使状态Ss k ∈按转移律由初始状态()1,1,1,11=s 经有限步n 到达状态()0,0,0,01=+n s 。

根据题设条件:影响安全渡河的元素是猫、鸡、米,这3个元素中取2个元素的组合一共有=3种,分别为猫+鸡,猫+米,鸡+米。

其中“猫+鸡”和“鸡+米”组合不安全,而“猫+米”的组合是安全的。

第一次渡河需带3个元素中的1个元素,另外2个元素留在岸上,而留在岸上的3种组合中只有“猫+米”的组合是安全的,可见第一次渡河只能带鸡,即安全方案只有U=(1,0,1,0),第二次将米或者猫带到对岸,把鸡带回,第三次将猫或者米带到对岸,第四次将鸡带过河,至此,猫、鸡、米均安全过河。

具体有2种执行方案如下:上述方案直观表示如下:3.购房贷款问题(续)在1.2.3节“购房贷款”的问题中,我们讨论了小王夫妇借贷还贷的方式。

现进一步讨论此问题。

某借贷公司的广告称,对于贷款期在20年以上的客户(利率仍按06%/月计算),他们帮你提前3年还清贷款。

但条件如下:(1)每半个月付款一次,但付款额不增加,即一次付款额是原付给银行还款额的1/2;(2)因为增加必要的档案、文书等管理工作,因此,要预付给借贷公司贷款总额10%的佣金。

试分析,小王夫妇是否要请这家借贷公司帮助还款。

答:设每月还款额为x,利率r,贷款额A0,总还款月数为N,总还款额为X。

若贷款20万,20年期,则:x===1574.7元X=x×N=1574.7×240=377928元若请这家借贷公司还款,则还款总额X1=x(N-3×12)+ 10%A0=1574.7×204+20000=341238.8元但小王夫妇预付给借贷公司20000元的佣金,如果把这20000元作为首付,则小王夫妇只需贷款18万,则总还款额:X2=X(1-10%)=340135.2元还款月数:N2=X2÷x=340135.2÷1574.7=216月=18年X1>X2,小王夫妇不需用借贷公司还款。

4.冷却定律按照Newton冷却定律,温度为T的物体在温度为T0(T0<T)的环境中冷却的速度与温差T-T0成正比。

用此定律建立相应的微分方程模型。

如果空气的温度是20℃,且沸腾的水在20分钟内冷却到60℃,那么水温降低到30℃需要多长时间?答:首先,牛顿冷却定律为温度为T(t)的物体在温度的环境中冷却的速度与温度差成正比。

所以,得出微分方程,K为比例常数。

任意时刻t,物体的温度为,C为常数根据已知条件,T0=20℃,记t=0时刻,初始水温T(0)=100℃,20分钟后水温T(20)=60℃。

则:求解函数得,k=-0.034657,C=80,即当水温降低到30℃时,t=-ln((30-20)÷80)÷0.034657=60分钟。

5.锻炼想象力、洞察力和判断力的问题(只简单回答出理由即可)(1)某人早8时从山下旅店出发沿一条山路上山,下午5时到达山顶并留宿,次日8时沿同一山路下山,下午5时回到酒店。

该人比在两天中的同一时刻经过路径中的同一地点,为什么?答:题中人记作A,假定有另一个人B,第二天复制A第一天的上山过程(即任何时刻的足迹均相同),那么问题相当于A、B两人分别从一段路的两端同时出发,相向行走,同时到达另一端。

则A、B二人必然会在路途中相遇,相遇点即为A在两天中的同一时刻经过路径中的同一地点。

也可用下图加以说明,两条曲线分别为上山和下山的曲线,两天都是同一时刻出发,同一时刻到达,无论曲线如何变化,总会有一个相交点,相交点即为为A在两天中的同一时刻经过路径中的同一地点。

(2)甲乙两站之间有汽车相通,每隔10分钟甲乙两站相互发一趟车,但发车时刻不一定相同,甲乙两站之间有一中间站丙,某人每天在随机时刻到达丙站,并搭乘最先经过丙站的那趟车,结果发现100天中约有90天到达甲站,大约有10天到达乙站。

问开往甲乙两站的汽车经过丙站的时刻是如何安排的?答:由题意可知,搭乘乙站开往甲站的车概率是90%,搭乘甲站开往乙站的车概率是10%。

说明90%的时间等来的是乙站开往甲站的车,10%的时间等来的是甲站开往乙站的车,而发车时间间隔是10分钟,进而说明:10分钟内有9分钟等来的是乙站开往甲站的车,有1分钟等来的是甲站开往乙站的车。

因此可以得出结论:乙站开往甲站的车到站后1分钟,甲站开往乙站的车到站。

开往甲乙两站的汽车经过丙站的时刻可表示如下(H0为乙站开往甲站的首班车时刻):H0乙甲H0+1 甲乙H0+10 乙甲H0+11 甲乙H0+15 乙甲H0+16 甲乙…… ……(3)张先生家住在A市,在B市工作,每天下班后他乘城际火车于18:00抵达A市火车站,他妻子驾车至火车站接他回家。

一日他提前下班,乘早一班火车于17:30抵达A市火车站,随即步行回家,他妻子像往常一样驾车前来,在半路相遇将他接回家。

到家时张先生发现比往常提前了10分钟,问张先生步行了多长时间?答:张先生从火车站比往常早出发30分钟,早到家10分钟,则其步行路段用时比往常多了20分钟。

他妻子从家出发时间不变,早到家10分钟,则其少行驶了10分钟的路程,即张先生步行路段的往返路程,那么张先生步行路段的单程行车时间为5分钟。

张先生步行的时间应为比往常多用的时间和步行路段的单程行车时间之和,即25分钟。

(4)一男孩和一女孩分别在距家2公里和1公里且方向相反的两所学校上学,每天同时放学后分别以每小时4公里和每小时2公里的速度步行回家。

一小狗以每小时6公里的速度由男孩处奔向女孩,又从女孩处奔向男孩,如此往返直至回到家中。

问小狗奔波了多少路程。

如果男孩和女孩上学时,小狗也往返奔波在他们中间,问他们到达学校时小狗在何处?答:由于他们到家的用时都是2÷4=0.5小时,所以小狗奔波的路程为6×0.5=3公里。

上学过程是放学过程的逆向过程,所以在他们到达学校时小狗的位置为放学时的起始位置,即男孩所在的学校。

加分实验(公平投票问题)某部门推出一专项基金目的在于培养优秀人才,根据评比结果确定资助的额度。

许多单位的优秀者都申请了该基金,于是该基金的委员会聘请了数名专家,按照如下规则进行评比。

1.为了公平性,评委对本单位选手不给分;2.每位评委对每位参与申请的人(除本单位选手外)都必须打分,且不打相同的分;3.评委打分方法为给参加申请的人排序,根据优劣分别记1分、2分、…以此类推。

4.评判结束后,求出各选手的平均分,按平均分从低到高排序,依次确定本次评比的名次,即平均分最低者获得资助最高,依次类推。

本次基金申请中,甲所在单位有一名评委,这位评委将不参加对选手甲的评判,其他选手没有类似情况,评审结束后选手甲觉得这种评比规则对他不公平。

问选手甲的抱怨是否有道理?若不公平,能否做出修正来解决选手甲的抱怨?答:令评审个数为n,根据评分规则,在一般情况下评委给出的评分分别为,,…,所以该选手的平均得分为=由于回避规则,甲的单位不得为甲评分,则甲的得分与一般选手__________________________________________________ __________________________________________________ 有差别,应为=二者之差为由上式可知:只有当每个评委的打分均相等时,两者之差才恒等于0,其余情况则不恒等于0。

根据题意,每个评委的对每个选手的打分是不等的,所以此规则对于选手甲确实有不公平之处。

若在尽可能公平的情况下制定得分规则,可判断出甲所得到的公平平均分应介于之间。

可以用的几何平均值来定义,即K=。

相关文档
最新文档