数学建模作业求解常微分方程和人口模型问题

数学建模作业求解常微分方程和人口模型问题
数学建模作业求解常微分方程和人口模型问题

实验报告

课程名称:数学建模

课题名称:求解常微分方程与人口模型

专业:信息与计算科学

姓名:胡家炜

班级: 123132

完成日期: 2016 年 6 月 10 日

一.求解微分方程的通解

(1). dsolve('2*x^2*y*Dy=y^2+1','x')

ans =

(exp(C3 - 1/x) - 1)^(1/2)

-(exp(C3 - 1/x) - 1)^(1/2)

i

-i

(2). dsolve('Dy=(y+x)/(y-x)','x')

ans =

x + 2^(1/2)*(x^2 + C12)^(1/2)

x - 2^(1/2)*(x^2 + C12)^(1/2)

(3). dsolve('Dy=cos(y/x)+y/x','x')

ans =

(pi*x)/2-x*log(-(exp(C25 + log(x)) - i) /(exp(C25 + log(x))*i - 1))*i (4). dsolve('(x*cos(y)+sin(2*y))*Dy=1','x')

ans =

-asin(x/2 + lambertw(0, -(C30*exp(- x/2 - 1))/2) + 1)

(5). dsolve('D2y+3*Dy-y=exp(x)*cos(2*x)','x')

ans =

C32*exp(x*(13^(1/2)/2 - 3/2)) + C33*exp(-x*(13^(1/2)/2 + 3/2)) + (13^(1/2)*exp(x*(13^(1/2)/2-3/2))*exp((5*x)/2(13^(1/2)*x)/2)*

(2*sin(2*x) - cos(2*x)*(13^(1/2)/2 - 5/2)))/(13*((13^(1/2)/2 - 5/2)^2 +4))-(13^(1/2)*exp(x*(13^(1/2)/2+3/2))*exp((5*x)/2

+(13^(1/2)*x)/2)*(2*sin(2*x)+cos(2*x)*(13^(1/2)/2+5/2)))

/(13*((13^(1/2)/2 + 5/2)^2 + 4))

(6)dsolve('D2y+4*y=x+1+sin(x)','x')

ans =

cos(2*x)*(cos(2*x)/4 - sin(2*x)/8 + sin(3*x)/12 - sin(x)/4 + (x*cos(2*x))/4 - 1/4) + sin(2*x)*(cos(2*x)/8 - cos(3*x)/12 + sin(2*x)/4 + cos(x)/4 + (x*sin(2*x))/4 + 1/8) + C35*cos(2*x) + C36*sin(2*x)

二.求初值问题的解

(1). dsolve('x^2+2*x*y-y^2+(y^2+2*x*y-x^2)*Dy=0','y(1)=1','x') ans =

(x*((- 4*x^2 + 4*x + 1)/x^2)^(1/2))/2 + 1/2

(2). dsolve('D2x+2*n*Dx+a^2*x=0','x(0)=x0','Dx(0)=V0')

ans =

(exp(-t*(n - (-(a + n)*(a - n))^(1/2)))*(V0 + n*x0 + x0*(-(a + n)*(a - n))^(1/2)))/(2*(-(a + n)*(a - n))^(1/2)) - (exp(-t*(n + (-(a + n)*(a - n))^(1/2)))*(V0 + n*x0 - x0*(-(a + n)*(a - n))^(1/2)))/(2*(-(a + n)*(a - n))^(1/2))

三.给出函数f(x)=sinx+cosx在x=0点的7阶taylor展开式以及在x=1处的5阶taylor展开式。(1). sym x;taylor(exp(x)*sin(x)+2^x*cos(x),7,0)

ans =

(log(2)^2/48 - log(2)^4/48 + log(2)^6/720 - 1/80)*x^6 + (log(2)/24 - log(2)^3/12 + log(2)^5/120 - 1/30)*x^5 + (log(2)^4/24 - log(2)^2/4 + 1/24)*x^4 + (log(2)^3/6 - log(2)/2 + 1/3)*x^3 + (log(2)^2/2 + 1/2)*x^2 + (log(2) + 1)*x + 1

(2). sym x;taylor(exp(x)*sin(x)+2^x*cos(x),5,1)

ans =

2*cos(1) + exp(1)*sin(1) - (x - 1)^2*(cos(1) - cos(1)*exp(1) + 2*log(2)*sin(1)-cos(1)*log(2)^2)+(x-1)^3*(sin(1)/3+(cos(1)* exp(1)) /3-cos(1)*log(2)-(exp(1)*sin(1))/3+(cos(1)*log(2)^3)/3 - log(2)^2* sin(1)) + (x - 1)^4*(cos(1)/12 - (exp(1)*sin(1))/6 + (log(2)*sin(1))/3 - (cos(1)*log(2)^2)/2 + (cos(1)*log(2)^4)/12 - (log(2)^3*sin(1))/3) + (x - 1)*(cos(1)*exp(1) - 2*sin(1) + 2*cos(1)*log(2) + exp(1)*sin(1))

四.判别下列级数的敛散性,若收敛求其和. (1). sym n;symsum(1/(2*n-1),n,1,inf)

ans =

Inf

因此不收敛

(2). sym n;symsum(tan(pi/(2*n*sqrt(n+1))),n,1,inf)

ans =

sum(tan(pi/(2*n*(n + 1)^(1/2))), n == 1..Inf)

limit(tan(pi/(2*n*(n + 1)^(1/2)))/(1/n^2),n,inf)

ans =

Inf

因此不收敛

(3). sym n;symsum((-1)*(n/2)*(1/(n*sqrt(n+1))),n,1,inf) ans =

-sum(1/(n + 1)^(1/2), n == 1..Inf)/2

limit(1/(n + 1)^(1/2)/(1/n^2),n,inf)

ans =

Inf

因此不收敛

(4). sym n;symsum((-1)^n*(1/(n*log(n))),n,1,inf)

ans =

sum((-1)^n/(n*log(n)), n == 1..Inf)

limit((-1)^n*(1/(n*log(n)))/(1/n^2),n,inf)

ans =

NaN

因此不收敛

五.求幂级数的和函数

syms x n;symsum((-1)^n*x^n/sqrt(n^2-n),n,2,inf)

ans =

sum(((-1)^n*x^n)/(n^2 - n)^(1/2), n == 2..Inf)

六.求函数项级数的和函数

syms x n;symsum((-1)^n*sin(pi/(2^n))*x^n,n,1,inf)

ans =

sum((-1)^n*x^n*sin(1/2^n*pi), n == 1..Inf)

七.人口模型

一、实验名称

建立适合的拟合模型,预测人口增长

二、实验目的

认识人口数量的变化规律,建立人口模型,运用专用的拟合函数polyfit求解,并作出较准确的预报

三、实验原理

对于情况较复杂的实际问题(因素不易化简,作用机理不详)可直接使用数据组建模 ,寻找简单的因果变量之间的数量关系,从而对未知的情形作预报。拟合模型的组建是

通过对有关变量的观测数据的观察、分析和选择恰当的数学表达方式得到的拟合模型组建的实质是数据拟合的精度和数学表达式简化程度间的一个折中。

四、实验题目

建立Logistic人口阻滞增长模型,利用表1中的数据分别根据从1954年、1963年、1980年到2005年三组总人口数据建立模型,进行预测我国未来50年的人口情况.并把预测结果与《国家人口发展战略研究报告》中提供的预测值进行分析比较。分析那个时间段数据预测的效果好?并结合中国实情分析原

因。

建立模型

阻滞增长模型(Logistic 模型)阻滞增长模型的原理:阻滞增长模型是考虑到自然资源、环境条件等因素对人口增长的阻滞作用,对指数增长模型的基本假设进行修改后得到的。阻滞作用体现在对人口增长率r 的影响上,使得r 随着人口数量x 的增加而下降。若将r 表示为x 的函数)(x r 。则它应是减函数。于是有: 0)0(,)(x x x x r dt dx == (1)

对)(x r 的一个最简单的假定是,设)(x r 为x 的线性函数,即

)0,0()(>>-=s r sx r x r (2)

设自然资源和环境条件所能容纳的最大人口数量m x ,当m x x =时人口不再

增长,即增长率0)(=m x r ,代入(2)式得

m x r s =,于是(2)式为 )1()(m x x r x r -= (3)

将(3)代入方程(1)得: ?????=-=0)0()1(x x x x rx dt dx m (4)

解得:

rt m m e x x x t x --+=

)1(

1)(0 (5)

模型求解

用Matlab 求解,程序如下:

t=1954:1:2005;

x=[60.2,61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129

.988,130.756];

x1=[60.2,61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80 .7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008 ,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,11 9.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,12 9.988];

x2=[61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83 ,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104. 357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85, 121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129.988 ,130.756];

dx=(x2-x1)/x2;

a=polyfit(x2,dx,1);

r=a(2),xm=-r/a(1)%求出xm和r

x0=61.5;

f=inline('xm./(1+(xm/x0-1)*exp(-r*(t-1954)))','t','xm','r','x0');%定义函数

plot(t,f(t,xm,r,x0),'-r',t,x,'+b');

title('1954-2005年实际人口与理论值的比较')

x2010=f(2010,xm,r,x0)

x2020=f(2020,xm,r,x0)

x2033=f(2033,xm,r,x0)

解得:x(m)=180.9516(千万),r=0.0327/(年),x(0)=61.5

人口增长模型的确定

题目:人口增长模型的确定 摘要 人口问题已成为当前世界上最普遍关注的问题之一,人口增长规律的发现以及人口增长的预测问题对一个国家制定长远的发展规划有着非常重要的意义。本文分别使用了马尔萨斯人口指数增长模型和阻滞增长模型,以美国1790-1980年间每隔10年的人口数量为依据,对接下来的每隔十年进行了预测五次人口数量。通过对比我们可以发现阻滞增长模型在预测准确度方面要明显优于原始的马尔萨斯人口指数增长模型。 关键词:人口增长;马尔萨斯人口指数增长模型;阻滞增长模型;人口预测

一、问题重述 1.1 问题背景 1790-1980年间美国每隔10年的人口记录如下表所示。 表1 人口记录表 1.2 问题提出 我们需要解决以下问题: 1.试用以上数据建立马尔萨斯(Malthus)人口指数增长模型,并对接下来的每隔十年预测五次人口数量,并查阅实际数据进行比对分析。 2.如果数据不相符,再对以上模型进行改进,寻找更为合适的模型进行预测,并对两次预测结果进行对比分析。 3.查阅资料找出中国人口与表1同时期的人口数量,用以上建立的两个模型进行人口预测与分析。 二、问题分析 首先,我们运用Matlab 软件绘制出1790到1980年的美国人口数据图,如图1。 17801800182018401860188019001920194019601980 050 100 150 200 250

图1 1790到1980年的美国人口数据图 从图表中我们可以清晰地看到人口数在1790—1980年是呈增长趋势的,而且我们很容易发现上述图表和我们学过指数函数的图表有很大的相似性,所以我们很自然想到建立指数模型。因此我们首先建立马尔萨斯模型,马尔萨斯生物总数增长定律指出:在孤立的生物群体中,生物总数N的变化率与生物总数成正比。 三、问题假设 为简化问题,我们做出如下假设: (1)在模型中预期的时间内,人口不会因发生大的自然灾害,突发事件或战争而受到大的影响; (2)所给出的数据具有代表性,能够反映普遍情况; (3)一段时间内我国人口死亡率不发生大的波动; (4)在查阅的资料与文献中,所得数据可信; (5)假设人口净增长率为常数。 四、变量说明 在此,对本文所使用的符号进行定义。 表2 变量说明 符号符号说明 N(0)起始年人口容纳量 N(t)t年后人口容纳量 t年份 r增长率 五、模型建立 5.1 问题一:马尔萨斯(Malthus)人口指数增长模型 设:t表示年份(起始年份t=0),r表示人口增长率,N(t)表示t年后的人口数量。 当考察一个国家或一个很大地区的人口时,N(t)是很大的整数。为了利用微积分这一数学工具,将N(t)视为连续、可微函数。记初始时刻(t=0)的人口为N(0),人口增长率为r,r是单位时间内N(t)的增量与N(t)的比例系数。根据r是常数的基本假设,于是N(t)满足如下的微分方程: dN(t)/dt=r*N(t) (5-1) 由这个线性常系数微分方程容易解出: N(t)=N(0)e rt(5-2) 表明人口将按指数规律无限增长(r>0)。将以t年为单位,上式表明,人口以e r为公

常微分方程在数学建模中的应用(免费版)

常微分方程在数学建模中的应用 这里介绍几个典型的用微分方程建立数学模型的例子. 一、人口预测模型 由于资源的有限性,当今世界各国都注意有计划地控制人口的增长,为了得到人口预测模型,必须首先搞清影响人口增长的因素,而影响人口增长的因素很多,如人口的自然出生率、人口的自然死亡率、人口的迁移、自然灾害、战争等诸多因素,如果一开始就把所有因素都考虑进去,则无从下手.因此,先把问题简化,建立比较粗糙的模型,再逐步修改,得到较完善的模型. 例1( 马尔萨斯 (Malthus ) 模型) 英国人口统计学家马尔萨斯(1766—1834)在担任牧师期间,查看了教堂100多年人口出生统计资料,发现人口出生率是一个常数,于1789年在《人口原理》一书中提出了闻名于世的马尔萨斯人口模型,他的基本假设是:在人口自然增长过程中,净相对增长(出生率与死亡率之差)是常数,即单位时间内人口的增长量与人口成正比,比例系数设为r ,在此假设下,推导并求解人口随时间变化的数学模型. 解 设时刻t 的人口为)(t N ,把)(t N 当作连续、可微函数处理(因人口总数很大,可近似地这样处理,此乃离散变量连续化处理),据马尔萨斯的假设,在t 到t t ?+时间段内,人口的增长量为 t t rN t N t t N ?=-?+)()()(, 并设0t t =时刻的人口为0N ,于是 ?????==. , 00)(d d N t N rN t N 这就是马尔萨斯人口模型,用分离变量法易求出其解为 )(00e )(t t r N t N -=, 此式表明人口以指数规律随时间无限增长. 模型检验:据估计1961年地球上的人口总数为9 1006.3?,而在以后7年中,人口总数以每年2%的速度增长,这样19610=t ,901006.3?=N ,02.0=r ,于是 ) 1961(02.09 e 1006.3)(-?=t t N . 这个公式非常准确地反映了在1700—1961年间世界人口总数.因为,这期间地球上的人 口大约每35年翻一番,而上式断定34.6年增加一倍(请读者证明这一点). 但是,后来人们以美国人口为例,用马尔萨斯模型计算结果与人口资料比较,却发现有很大的差异,尤其是在用此模型预测较遥远的未来地球人口总数时,发现更令人不可思议的问题,如按此模型计算,到2670年,地球上将有36 000亿人口.如果地球表面全是陆地(事实上,地球表面还有80%被水覆盖),我们也只得互相踩着肩膀站成两层了,这是非常荒谬的,因此,这一模型应该修改. 例2(逻辑Logistic 模型) 马尔萨斯模型为什么不能预测未来的人口呢?这主要是地

数学建模大作业

兰州交通大学 数学建模大作业 学院:机电工程学院 班级:车辆093 学号:200903812 姓名:刘键学号:200903813 姓名:杨海斌学号:200903814 姓名:彭福泰学号:200903815 姓名:程二永学号:200903816 姓名:屈辉

高速公路问题 1 实验案例 (2) 1.1 高速公路问题(简化) (2) 1.1.1 问题分析 (3) 1.1.2 变量说明 (3) 1.1.3 模型假设 (3) 1.1.4 模型建立 (3) 1.1.5 模型求解 (4) 1.1.6 求解模型的程序 (4) 1实验案例 1.1 高速公路问题(简化) A城和B城之间准备建一条高速公路,B城位于A城正南20公里和正东30公里交汇处,它们之间有东西走向连绵起伏的山脉。公路造价与地形特点有关,图4.2.4给出了整个地区的大致地貌情况,显示可分为三条沿东西方向的地形带。 你的任务是建立一个数学模型,在给定三种地形上每公里的建造费用的情况下,确定最便宜的路线。图中直线AB显然是路径最短的,但不一定最便宜。而路径ARSB过山地的路段最短,但是否是最好的路径呢? A B 图8.2 高速公路修建地段

1.1.1 问题分析 在建设高速公路时,总是希望建造费用最小。如果要建造的起点、终点在同一地貌 中,那么最佳路线则是两点间连接的线段,这样费用则最省。因此本问题是一个典型的最优化问题,以建造费用最小为目标,需要做出的决策则是确定在各个地貌交界处的汇合点。 1.1.2 变量说明 i x :在第i 个汇合点上的横坐标(以左下角为直角坐标原点),i =1,2,…,4;x 5=30(指目的地B 点的横坐标) x=[x 1,x 2,x 3,x 4]T l i :第i 段南北方向的长度(i =1,2, (5) S i :在第i 段上地所建公路的长度(i =1,2, (5) 由问题分析可知, () ()() () 2 542552 432442 322332212 222 1211x x l S x x l S x x l S x x l S x l S -+=-+=-+=-+=+= C 1:平原每公里的造价(单位:万元/公里) C 2:高地每公里的造价(单位:万元/公里) C 3:高山每公里的造价(单位:万元/公里) 1.1.3 模型假设 1、 假设在相同地貌中修建高速公路,建造费用与公路长度成正比; 2、 假设在相同地貌中修建高速公路在一条直线上。在理论上,可以使得建造费用最少, 当然实际中一般达不到。 1.1.4 模型建立 在A 城与B 城之间建造一条高速公路的问题可以转化为下面的非线性规划模型。优化目标是在A 城与B 城之间建造高速公路的费用。 () 4,3,2,1300. .)(min 5142332211=≤≤++++=i x t s S C S C S C S C S C x f i

数学建模作业——实验1

数学建模作业——实验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)>0 h(π)=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 ----=对岸状态。安全状态集合为 :

数学建模人口模型

摘要 以2010年11月1日零时为标准时点,中国大陆31个省、自治区、直辖市和现役军人的人口共13.397亿。13亿是一个忧虑的数字。13亿人要吃饭、要穿衣、要上学、要就业、要住房……,消费的需求乘以13亿,就是一个庞大的数目,而我国的耕地、水资源、森林以及矿产资源本来就稀缺,再除以13亿,就少得可怜。平均每人耕地面积只有1.4亩,水资源只相当于世界人均水平的1/4…….、 中国是世界上人口最多的发展中国家,人口多,底子薄,人均耕地少,人均占有资源相对不足,是我国的基本国情,人口问题一直是制约中国经济发展的首要因素。当前中国的人口存在着最为明显的三大特点:(1)人口基数大,人口数量的控制难度仍很大。(2)人口整体素质不高,特别是县域及以下农村人口素质普遍偏低。(3)人口结构不合理,城乡差别、地区差别和人口素质差别很大。 人口数量、质量和年龄分布直接影响一个地区的经济发展、资源配置、社会保障、社会稳定和城市活力。在我国现代化进程中,必须实现人口与经济、社会、资源、环境协调发展和可持续发展,进一步控制人口数量,提高人口质量,改善人口结构。对此,单纯的人口数量控制(如已实施多年的计划生育)不能体现人口规划的科学性。政府部门需要更详细、更系统的人口分析技术,为人口发展策略的制定提供指导和依据。 我国是世界第一人口大国,地球上每九个人中就有二个中国人,在20世纪的一段时间内我国人口的增长速度过快,如下表: 有效地控制人口的增长,不仅是使我国全面进入小康社会、到21世纪中叶建成富强民主文明的社会主义国家的需要,而且对于全人类社会的美好理想来说,也是我们义不容辞的责任。 长期以来,对人口年龄结构的研究仅限于粗线条的定性分析,只能预测年龄结构分布的大致范围,无法用于分析年龄结构的具体形态。随着对人口规划精准度要求的提高,通过数学方法来定量计算各种人口指数的方法日益受到重视,这就是人口控制和预测。 我国人口问题已积重难返,对我国人口进行准确的预测是制定合理的社会经济发展规划

常微分方程在数学建模中的应用.

微分方程应用 1 引言 常微分方程的形成与发展和很多学科有着密切的联系,例如力学、天文学、物理学等.数学的其他分支的快速发展,产生出很多新兴学科,这些新兴学科的产生都对常微分方程的发展有着深刻的影响,而且当前计算机的快速发展更是为常微分方程的应用及理论研究提供了非常有力的工具. 数学解决实际问题就必须建立模型,而数学建模就是把数学语言描述实际现象的过程.利用数学去解决各类实际问题时,建立数学模型是十分重要的一步,但是也是最困难的一步.建立数学模型的过程,是把错综复杂的实际问题简化、抽象为合理的数学结构的过程.要通过大量调查、收集相关数据资料,观察和研究实际对象的固有特征和内在规律,抓住问题的主要矛盾,建立起反映实际问题的数量关系,然后利用数学的理论和方法去分析和解决问题. 因此本文先简要介绍了如何建立微分方程模型,并通过具体的实例来简单地介绍了微分方程在数学建模中的应用. 2 数学模型简介 通常我们把现实问题的一个模拟称为模型.如交通图、地质图、航空模型和建筑模型等.利用字母、数学及其它数学符号建立起来的等式或不等式以及图表、图象、框图等来模拟现实的模型称为数学模型.数学模型在实际生活中经常碰到,如求不规则图形的面积,可建立定积分的数学模型,求变化率的问题可建立导数模型,统计学中抽样调查,买彩票中奖的概率问题等等.学会建立数学模型对解决实际生活问题会有很大的帮助. 建立数学模型是沟通摆在面前的实际问题与数学工具之间联系的一座必不可少的桥梁.随着科学技术的进步,特别是电子计算机技术的迅速发展,数学已经渗透到从自然科学技术到工农业生产建设,从经济生活到社会生活的各个领域.一般地说,当实际问题需要我们对所研究的现实对象提供分析、预报、决策、控制等方面的定量结果时,往往都离不开数学的应用,而建立数学模型则是这个过程的关键环节. 3 常微分方程模型 3.1 常微分方程的简介

数学建模作业

郑重声明: 本作业仅供参考,可能会有错误,请自己甄别。 应用运筹学作业 6.某工厂生产A,B,C,D四种产品,加工这些产品一般需要经刨、磨、钻、镗四道工序,每种产品在各工序加工时所需设备台时如表1-18所示,设每月工作25天,每天工作8小时,且该厂有刨床、磨床、钻床、镗床各一台。问:如何安排生产,才能使月利润最大?又如A,B,C,D四种产品,每月最大的销售量分别为300件、350件、200件和400件,则该问题的线性规划问题又该如何? 1234 四种产品的数量,则得目标函数: Max=(200?150)x1+(130?100)x2+(150?120)x3+(230?200)x4 =50x1+30x2+30x3+30x4 生产四种产品所用时间: (0.3+0.9+0.7+0.4)x1+(0.5+0.5+0.5+0.5)x2+(0.2+0.7+0.4+ 0.8)x3+(0.4+0.8+0.6+0.7)x4≤25×8 即:2.3x1+2.0x2+2.1x3+2.5x4≤200 又产品数量不可能为负,所以:x i≥0(i=1,2,3,4) 综上,该问题的线性规划模型如下: Max Z=50x1+30x2+30x3+30x4 S.T.{2.3x1+2.0x2+2.1x3+2.5x4≤200 x i≥0(i=1,2,3,4) 下求解目标函数的最优解: max=50*x1+30*x2+30*x3+30*x4; 2.3*x1+2.0*x2+2.1*x3+2.5*x4<200; Global optimal solution found. Objective value: 4347.826 Total solver iterations: 0 Variable Value Reduced Cost X1 86.95652 0.000000 X2 0.000000 13.47826 X3 0.000000 15.65217

matlab曲线拟合人口增长模型及其数量预测

实验目的 [1] 学习由实际问题去建立数学模型的全过程; [2] 训练综合应用数学模型、微分方程、函数拟合和预测的知识分析和解决实际问题; [3] 应用matlab 软件求解微分方程、作图、函数拟合等功能,设计matlab 程序来求解 其中的数学模型; [4] 提高论文写作、文字处理、排版等方面的能力; 通过完成该实验,学习和实践由简单到复杂,逐步求精的建模思想,学习如何建立反映人口增长规律的数学模型,学习在求解最小二乘拟合问题不收敛时,如何调整初值,变换函数和数据使优化迭代过程收敛。 应用实验(或综合实验) 一、实验内容 从1790—1980年间美国每隔10年的人口记录如表综2.1所示: 表综2.1 用以上数据检验马尔萨斯(Malthus)人口指数增长模型,根据检验结果进一步讨论马尔萨斯人口模型的改进,并利用至少两种模型来预测美国2010年的人口数量。 二、问题分析 1:Malthus 模型的基本假设是:人口的增长率为常数,记为 r 。记时刻t 的人口为x (t ),(即x (t )为模型的状态变量)且初始时刻的人口为x 0,于是得到如下微分方程: ?????==0 )0(d d x x rx t x 2:阻滞增长模型(或Logistic 模型) 由于资源、环境等因素对人口增长的阻滞作用,人 口增长到一定数量后,增长率会下降,假设人口的增长率为x 的减函数,如设r(x)=r(1-x/x m ),其中r 为固有增长率(x 很小时),x m 为人口容量(资源、环境能容纳的最大数量),于是得到如下微分方程: ?? ???=-=0)0()1(d d x x x x rx t x m

数学建模之微分方程建模与平衡点理论

微分方程 列微分方程常用的方法: (1)根据规律列方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律来建 立微分方程模型。 (2)微元分析法 利用已知的定理与规律寻找微元之间的关系式,与第一种方法不同的是对 微元而不是直接对函数及其导数应用规律。 (3)模拟近似法 在生物、经济等学科的实际问题中,许多现象的规律性不很清楚,即使有 所了解也是极其复杂的,建模时在不同的假设下去模拟实际的现象,建立能 近似反映问题的微分方程,然后从数学上求解或分析所建方程及其解的性 质,再去同实际情况对比,检验此模型能否刻画、模拟某些实际现象。 一、模型的建立与求解 1.1传染病模型 (1)基础模型 假设:t时刻病人人数() x t连续可微。每天每个病人有效接触(使病人治病的接触)的人数为λ,0 t=时有0x个病人。 +?病人人数增加 建模:t到t t

()()()x t t x t x t t λ+?-=? (1) 0,(0)dx x x x dt λ== (2) 解得: 0()t x t x e λ= (3) 所以,病人人数会随着t 的增加而无限增长,结论不符合实际。 (2)SI 模型 假设:1.疾病传播时期,总人数N 保持不变。人群分为两类,健康者占总人数的比例为s(t),病人占总人数的比例为i(t)。 2.每位病人每天平均有效接触λ人,λ为日接触率。有效接触后健康者变为病人。 依据:患病人数的变化率=Ni(t)(原患病人数)* λs(t)(每个病人每天使健康人变为病人的人数) 建模: di N Nsi dt λ= (4) 由于 ()()1s t i t += (5) 设t=0时刻病人所占的比例为0i ,则可建立Logistic 模型 0(1),(0)di i i i i dt λ=-= (6)

数学建模期末大作业

数学建模期末大作业论文 题目:A题美好的一天 组长:何曦(2014112739) 组员:李颖(2014112747)张楚良(2014112740) 班级:交通工程三班 指导老师:陈崇双

美好的一天 摘要 关键字:Dijkstra算法多目标规划有向赋权图 MATLAB SPSS

1 问题的重述 Hello!大家好,我是没头脑,住在西南宇宙大学巨偏远的新校区(节点22)。明天我一个外地同学来找我玩,TA叫不高兴,是个镁铝\帅锅,期待ing。我想陪TA在城里转转,当然是去些不怎么花钱的地方啦~~。目前想到的有林湾步行街(节点76)、郫郫公园(节点91),大川博物院(节点72)。交通嘛,只坐公交车好了,反正公交比较发达,你能想出来的路线都有车啊。另外,进城顺便办两件事,去老校区财务处一趟(节点50),还要去新东方(节点34)找我们宿舍老三,他抽奖中了两张电影票,我要霸占过来明晚吃了饭跟TA一起看。电影院嘛,TASHIWODE电影院(节点54)不错,比较便宜哈。我攒了很久的钱,订了明晚开心面馆(节点63)的烛光晚餐,额哈哈,为了TA,破费一下也是可以的哈。哦,对了,老三说了,他明天一整天都上课,只有中午休息的时候能接见我给我票。 我主要是想请教一下各位大神: 1)明天我应该怎么安排路线才能够让花在坐车上的时间最少? 2)考虑到可能堵车啊,TA比较没耐心啊,因为TA叫不高兴嘛。尤其是堵车啊,等车啊,这种事,万一影响了气氛就悲剧了。我感觉路口越密的地方越容易堵,如果考虑这个,又应该怎么安排路线呢? 3)我们城比较挫啊,连地图也没有,Z老师搞地图测绘的,他有地图,跟他要他不给,只给了我一个破表格(见附件,一个文件有两页啊),说“你自己画吧”。帮我画一张地图吧,最好能标明我们要去的那几个地方和比较省时的路线啊,拜托了~ 2 问题的分析 2.1 对问题一的分析 问题一要求安排路线使得坐车花费的时间最少。 对于问题一,假设公交车的速度维持不变,要使花费的时间最少,则将问题转化为对最短路径的求解。求解最短路径使用Dijkstra算法很容易进行求解,在运用MATLAB编程,得到最优的一条路径,则这条路径所对应的时间即为最少用时。 2.2 对问题二的分析 问题二要求在考虑堵车的情况下,路口越密越容易发生拥堵,安排路线是乘车时间最短。 对于问题二,在问题的基础上增加了附加因素,即公交车的速度会因道路的密集程度而发生改变,从而问题一建立的基本Dijkstra算法对于问题二就不再适用了,因此对问题一的基本Dijkstra算法进行改进,并结合蚁群算法的机理与特点,运用MATLAB求解出最短路径,保证了花费时间的最少性。 2.3 对问题三的分析 问题三要求根据提供的附件,画出一张地图,标明要去的那几个地方和比较省时的路线。 对于问题三,在问题一和问题二的基础上,根据求解的结果,运用SPSS软件画出地图。

数学建模作业

习 题 1 1. 请编写绘制以下图形的MA TLAB 命令,并展示绘得的图形. (1) 221x y +=、224x y +=分别是椭圆2241x y +=的内切圆和外切圆. (2) 指数函数x y e =和对数函数ln y x =的图像关于直线y=x 对称. (3) 黎曼函数 1, (0)(0,1) 0 , (0,1), 0,1 q x p q q x y x x x =>∈?=? ∈=?当为既约分数且当为无理数且或者 的图像(要求分母q 的最大值由键盘输入). 3. 两个人玩双骰子游戏,一个人掷骰子,另一个人打赌掷骰子者不能掷出所需点数,输赢的规则如下:如果第一次掷出3或11点,打赌者赢;如果第一次掷出2、7或12点,打赌者输;如果第一次掷出4、5、6、8、9或10点,记住这个点数,继续掷骰子,如果不能在掷出7点之前再次掷出该点数,则打赌者赢. 请模拟双骰子游戏,要求写出算法和程序,估计打赌者赢的概率. 你能从理论上计算出打赌者赢的精确概率吗?请问随着试验次数的增加,这些概率收敛吗?

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) 通过变量替换,可以将属于非线性模型的指数增长模型转化成线性模型,并用MA TLAB 函数polyfit 进行计算,请说明转化成线性模型的详细过程,然后写出程序,给出拟合参数和误差平方和的计算结果,并展示拟合效果图. (3) 请分析指数增长模型非线性拟合和线性化拟合的结果有何区别?原因是什么? (4) 如果用阻滞增长模型00 () 00()()e r t t Nx x t x N x --= +-模拟美国人口从1790年至2000年的变化过程,请用MA TLAB 统计工具箱的函数nlinfit 计算阻滞增长模型的以下三个数据拟合问题: (i) 取定0x =3.9,0t =1790,拟合待定参数r 和N ; (ii) 取定0t =1790,拟合待定参数0x 、r 和N ; (iii) 拟合待定参数0t 、0x 、r 和N . 要求写出程序,给出拟合参数和误差平方和的计算结果,并展示误差平方和最小的拟合效果图. 年份 1790 1800 1810 1820 1830 1840 1850 1860 1870 1880 1890

数学建模创新思维大作业

数学建模创新思维课大作业 一、使用MATLAB 求解一下问题,请贴出代码. 1. cos 1000x mx y e =,求''y >>clear >>clc >> syms x m; >> y=exp(x)*cos(m*x/1000); >> dfdx2=diff(y,x,2) dfdx2 = exp(x)*cos((m*x)/1000) - (m*exp(x)*sin((m*x)/1000))/500 - (m^2*exp(x)*cos((m*x)/1000))/1000000 >> L=simplify(dfdx2) L = -(exp(x)*(2000*m*sin((m*x)/1000) - 1000000*cos((m*x)/1000) + m^2*cos((m*x)/1000)))/1000000 2.计算22 1100x y e dxdy +?? >> clear >> clc; >> syms x y >> L=int(int(exp(x^2+y^2),x,0,1),y,0,1) L = (pi*erfi(1)^2)/4 3. 计算4 224x dx m x +? >> clear; >> syms x m; >> f=x^4/(m^2+4*x^2); >> intf=int(f,x) intf =

(m^3*atan((2*x)/m))/32 - (m^2*x)/16 + x^3/12 >> L=simplify(intf) L = (m^3*atan((2*x)/m))/32 - (m^2*x)/16 + x^3/12 4. (10)cos ,x y e mx y =求 >> clear; >> syms x m; >> y=exp(x)*cos(m*x); >> L=diff(y,x,10); >> L=simplify(L) L = -exp(x)*(10*m*sin(m*x) - cos(m*x) + 45*m^2*cos(m*x) - 210*m^4*cos(m*x) + 210*m^6*cos(m*x) - 45*m^8*cos(m*x) + m^10*cos(m*x) - 120*m^3*sin(m*x) + 252*m^5*sin(m*x) - 120*m^7*sin(m*x) + 10*m^9*sin(m*x)) 5. 0x =的泰勒展式(最高次幂为4). >> clear; >> syms m x; >> y=sqrt(m/1000.0+x); >> y1=taylor(y,x,'order',5); >> L=simplify(y1) L = (10^(1/2)*(m^4 + 500*m^3*x - 125000*m^2*x^2 + 62500000*m*x^3 - 39062500000*x^4))/(100*m^(7/2)) 6. Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4, )n n n x x x n --=+=用循环语句编程 给出该数列的前20项(要求将结果用向量的形式给出)。 >> x=[1,1]; >> for n=3:20

数学建模logistic人口增长模型

Logistic 人口发展模型 一、题目描述 建立Logistic 人口阻滞增长模型 ,利用表1中的数据分别根据从1954年、1963年、1980年到2005年三组总人口数据建立模型,进行预测我国未来50年的人口情况.并把预测结果与《国家人口发展战略研究报告》中提供的预测值进行分析比较。分析那个时间段数据预测 表1 各年份全国总人口数(单位:千万) 二、建立模型 阻滞增长模型(Logistic 模型)阻滞增长模型的原理:阻滞增长模型是考虑到自然资源、环境条件等因素对人口增长的阻滞作用,对指数增长模型的基本假设进行修改后得到的。阻滞作用体现在对人口增长率r 的影响上,使得r 随着人口数量x 的增加而下降。若将r 表示为x 的函数)(x r 。则它应是减函数。于是有: )0(,)(x x x x r dt dx == (1) 对)(x r 的一个最简单的假定是,设)(x r 为x 的线性函数,即 ) 0,0()(>>-=s r sx r x r (2) 设自然资源和环境条件所能容纳的最大人口数量m x ,当 m x x =时人口不再增长,即增 长率 )(=m x r ,代入(2)式得 m x r s = ,于是(2)式为

)1()(m x x r x r -= (3) 将(3)代入方程(1)得: ?? ???=-=0 )0() 1(x x x x rx dt dx m (4) 解得: rt m m e x x x t x --+= )1( 1)(0 (5) 三、模型求解 用Matlab 求解,程序如下: t=1954:1:2005; x=[60.2,61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129.988,130.756]; x1=[60.2,61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129.988]; x2=[61.5,62.8,64.6,66,67.2,66.2,65.9,67.3,69.1,70.4,72.5,74.5,76.3,78.5,80.7,83,85.2,87.1,89.2,90.9,92.4,93.7,95,96.259,97.5,98.705,100.1,101.654,103.008,104.357,105.851,107.5,109.3,111.026,112.704,114.333,115.823,117.171,118.517,119.85,121.121,122.389,123.626,124.761,125.786,126.743,127.627,128.453,129.227,129.988,130.756]; dx=(x2-x1)./x2; a=polyfit(x2,dx,1); r=a(2),xm=-r/a(1)%求出xm 和r x0=61.5; f=inline('xm./(1+(xm/x0-1)*exp(-r*(t-1954)))','t','xm','r','x0');%定义函数 plot(t,f(t,xm,r,x0),'-r',t,x,'+b'); title('1954-2005年实际人口与理论值的比较')

数学建模之微分方程建模与平衡点理论

微分方程 列微分方程常用的方法: (1)根据规律列方程 利用数学、力学、物理、化学等学科中的定理或经过实验检验的规律来建立微分方程模型。 (2)微元分析法 利用已知的定理与规律寻找微元之间的关系式,与第一种方法不同的是对微元而不是直接对函数及其导数应用规律。 (3)模拟近似法 在生物、经济等学科的实际问题中,许多现象的规律性不很清楚,即使有所了解也是极其复杂的,建模时在不同的假设下去模拟实际的现象,建立能近似反映问题的微分方程,然后从数学上求解或分析所建方程及其解的性质,再去同实际情况对比,检验此模型能否刻画、模拟某些实际现象。 一、模型的建立与求解 1.1传染病模型 (1)基础模型 假设:t 时刻病人人数()x t 连续可微。每天每个病人有效接触(使病人治病的接触)的人数为λ,0t =时有0x 个病人。 建模:t 到t t +?病人人数增加 ()()()x t t x t x t t λ+?-=?(1) 0,(0)dx x x x dt λ==(2) 解得: 0()t x t x e λ=(3) 所以,病人人数会随着t 的增加而无限增长,结论不符合实际。 (2)SI 模型

假设:1.疾病传播时期,总人数N 保持不变。人群分为两类,健康者占总人数的比例为s(t),病人占总人数的比例为i(t)。 2.每位病人每天平均有效接触λ人,λ为日接触率。有效接触后健康者变为病人。 依据:患病人数的变化率=Ni(t)(原患病人数)*λs(t)(每个病人每天使健康人变为病人的人数) 建模: di N Nsi dt λ=(4) 由于 ()()1s t i t +=(5) 设t=0时刻病人所占的比例为0i ,则可建立Logistic 模型 0(1),(0)di i i i i dt λ=-=(6) 解得: 01()111kt i t e i -= ??+- ??? (7) 用Matlab 绘制图1()~i t t ,图2 ~di i dt 图形如下, 结论:在不考虑治愈情况下 ①当12i = 时di dt 达到最大值m di dt ?? ???,这时101ln 1m t i λ-??=- ???

数学建模作业43950

题目: 某种电子系统由三种元件组成,为了使系统正常运转,每个元件都必须工作良好,如果一个或多个元件安装备用件将会提高系统的可靠性,已知系统运转的可靠性为各元件可靠性的乘积,而每一个元件的可靠性是备用元件函数,具体数值见下表。 若全部备用件费用限制为150元,重量限制为20公斤,问每个元件安装多少备用件可使系统可靠性达到极大值? 要求:①作出全局最优解 ②列出这个问题的整数规划模型

假设:系统在运转过程中相互间没有影响,并且系统在增加备用件后 可靠性可以相互叠加。 建模: 设原件1,2,3需要的备用件各为x,y,z,可靠性为p分别为xp,yp,zp,整 个设备的可靠性为p,则由题意可得到: p=xp*yp*zp; 2x+4y+6z<=20; 20x+30y+40z<=150; x,y,z均为整数; 求出适当的x,y,z使p的值最大。 运用穷举法,编写C++程序如下: #include void main() { using namespace std; int x=0,y=0,z=0;//备à?用??零¢?件t数oy目? double xp[6]={0.5,0.6,0.7,0.8,0.9,1},yp[4]={0.6,0.75,0.95,1},zp[3]={0.7,0.9,1}; double p=0,temp=0;//可¨|靠?性? int i=0,j=0,k=0; cout<<"x\ty\tz\tp\n"; for(i=0;i<6;i++) { y=0; for(j=0;j<4;j++) { z=0; for(k=0;k<3;k++) {if((x+2*y+3*z<=10)&&(2*x+3*y+4*z<=15)) {temp=p; p=xp[x]*yp[y]*zp[z]; cout<

人口增长数学模型

软件学院 人口增长模型数学建模报告 专业:软件工程 班级:卓越131班 学号:201370044120 学生姓名:郭俊成 指导教师:于志云 2015 年11 月12 日 题目:计划生育政策调整对人口数量、结构及其影响的研究

摘要 本论文针对2007年国家人口发展战略研究课题组发布的《国家人口发展战略研究报告》中关于“计划生育实施以来,全国少生了4亿多人,使世界60亿人口日推迟4年”的论述做了研究。论文根据计划生育实施之前1949-1980年的人口普查数据,使用最小二乘法拟合并建立灰色预测模型,利用数学软件,预测出了如果未实行计划生育现今中国人口的数量,从而对研究报告中“少生4亿”的结论产生质疑。 同时,本论文针对2006年全国老龄工作委员会发布的《中国人口老龄化发展趋势预测研究报告》中关于“2051年,中国老年人口规模将达到峰值4.37亿,老龄化水平基本稳定在31%左右”的论述做了研究,根据近几年的人口老龄化程度、老龄人口比重、老龄人口数量、死亡率的变化等诸多因素,建立阻滞增长模型(Logistic模型),预测40年到70年的老龄人口数量和老龄化率,验证了报告中的关于老龄人口数目持续增加、数目庞大、老龄化严重的预测。 论文基于近期的计划生育调整、“单独二孩”政策的逐步实施、城镇化所导致的人口迁移等现象,结合江苏省的实际情况,利用差分方程模型、LESLIE矩阵,分析新政策对江苏人口数量的影响。论文从出生率着手,重点研究了新政策对江苏省14岁以下儿童、60岁以上老人的影响,分析了儿童和老人数量的变化对人口结构、教育改革、养老的直接影响作用。 关键字 单独二孩、人口老龄化、Logistic 模型、差分方程模型、LESLIE模型 一、问题描述

数学建模——微分方程的应用

第八节 数学建模——微分方程的应用举例 微分方程在物理学、力学、经济学和管理科学等实际问题中具有广泛的应用,本节我们将集中讨论微分方程的实际应用,尤其是微分方程经济学中的应用. 读者可从中感受到应用数学建模的理论和方法解决实际问题的魅力. 分布图示 ★衰变问题 ★逻辑斯谛方程 ★价格调整问题 ★人才分配问题 内容要点: 一、衰变问题 镭、铀等放射性元素因不断放射出各种射线而逐渐减少其质量, 这种现象称为放射性物质的衰变. 根据实验得知, 衰变速度与现存物质的质量成正比, 求放射性元素在时刻t 的质量. 用x 表示该放射性物质在时刻t 的质量, 则 dt dx 表示x 在时刻t 的衰变速度, 于是“衰变速度与现存的质量成正比”可表示为 .kx dt dx -= (8.1) 这是一个以x 为未知函数的一阶方程, 它就是放射性元素衰变的数学模型, 其中0>k 是比例常数, 称为衰变常数, 因元素的不同而异. 方程右端的负号表示当时间t 增加时, 质量x 减少. 解方程(8.1)得通解.kt Ce x -=若已知当0t t =时, ,0x x =代入通解kt Ce x -=中可得,00kt e x C -= 则可得到方程(8.1)特解 ,)(00t t k e x x --= 它反映了某种放射性元素衰变的规律. 注: 物理学中, 我们称放射性物质从最初的质量到衰变为该质量自身的一半所花费的时间为半衰期, 不同物质的半衰期差别极大. 如铀的普通同位素( U 238)的半衰期约为50亿年;通常的镭( Ra 226)的半衰期是上述放射性物质的特征, 然而半衰期却不依赖于该物质的初始量, 一克Ra 226 衰变成半克所需要的时间与一吨Ra 226衰变成半吨所需要的时间同样都是1600年, 正是这种事实才构成了确定考古发现日期时使用的著名的碳-14测验的基础.

2015年数学建模作业题

数学模型课程期末大作业题 要求: 1)选题方式:共53题,每个同学做一题,你要做的题目编号是你的学号mod52所得的值+1。(例如:你的学号为119084157,则你要做的题为mod(119084157,52)+1=50)。 2)该类题目基本为优划问题,要求提交一篇完整格式的建模论文,文字使用小四号宋体,公式用word的公式编辑器编写,正文中不得出现程序以及程序冗长的输出结果,程序以附录形式附在论文的后面,若为规划求解必须用lingo 集合形式编程,其它可用Matlab或Mathmatica编写。 3)论文以纸质文档提交,同时要交一份文章和程序电子文档,由班长统一收上来,我要验证程序。 1、生产安排问题 某厂拥有4台磨床,2台立式钻床,3台卧式钻床,一台镗床和一台刨床,用以生产7种产品,记作p1至p7。工厂收益规定作产品售价减去原材料费用之余。每种产品单件的收益及所需各机床的加工工时(以小时计)列于下表(表1): 表 到6月底每种产品有存货50件。 工厂每周工作6天,每天2班,每班8小时。 不需要考虑排队等待加工的问题。 在工厂计划问题中,各台机床的停工维修不是规定了月份,而是选择最合

适的月份维修。除了磨床外,每月机床在这6个月中的一个月中必须停工维修;6个月中4台磨床只有2台需要维修。扩展工厂计划模型,以使可作上述灵活安排维修时间的决策。停工时间的这种灵活性价值若何? 注意,可假设每月仅有24个工作日。 5、生产计划 某厂有4台磨床,2台立钻,3台水平钻,1台镗床和1台刨床,用来生产7种产品,已知生产单位各种产品所需的有关设备台时以及它们的利润如表所示: 台镗床,4月—1台立钻,5月—1台磨床和1台立钻,6月—1台刨床和1台水平钻,被维修的设备在当月内不能安排生产。又知从1月到6月份市场对上述7种产品最大需求量如表所示: 量均不得超过100件。现在无库存,要求6月末各种产品各贮存50件。若该厂每月工作24天,每天两班,每班8小时,假定不考虑产品在各种设备上的加工顺序,要求: (a)该厂如何安排计划,使总利润最大; (b)在什么价格的条件下,该厂可考虑租用或购买有关的设备。 34、瓶颈机器上的任务排序 在工厂车间中,经常会出现整个车间的生产能力取决于一台机器的情况(例如,仅有一台的某型号机床,生产线上速度最慢的机器等)。这台机器就称为关键机器或瓶颈机器。此时很重要的一点就是尽可能地优化此机器将要处理的任务计划。

相关文档
最新文档