数学建模与数学实验

1、数学建模与数学实验数学建模与数学实验总复习总复习西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验一、一、 绪论绪论 数学模型数学模型就是对于一个特定的对象为了一个特定目标，根据特就是对于一个特定的对象为了一个特定目标，根据特有的内在规律，做出一些必要的简化假设，运用适当的数学工具有的内在规律，做出一些必要的简化假设，运用适当的数学工具，得到的一个数学结构。数学结构可以是数学公式，算法、表格，得到的一个数学结构。数学结构可以是数学公式，算法、表格、图示等。、图示等。数学建模就是建立数学模型数学建模就是建立数学模型，建立数学模型的过程就是数学建，建立数学模型的过程就是数学建模的

2、过程。数学建模是一种数学的思考方法，是运用数学的语言模的过程。数学建模是一种数学的思考方法，是运用数学的语言和方法，通过抽象、简化建立能近似刻划并和方法，通过抽象、简化建立能近似刻划并 解决解决 实际问题的一实际问题的一种强有力的数学手段。种强有力的数学手段。数学建模需要一下几方面的能力：数学建模需要一下几方面的能力：数学知识的应用能力数学知识的应用能力、计计算机的运用能力算机的运用能力、论文的写作能力论文的写作能力。 模型的逼真性和可行性模型的逼真性和可行性：希望模型尽可能逼近研究对象，越：希望模型尽可能逼近研究对象，越逼真的模型越复杂，数学上有时难以处理，实用上不一定可行，逼真的模型越复杂

3、，数学上有时难以处理，实用上不一定可行，通常要在逼真性与可行性之间做出折衷与抉择；通常要在逼真性与可行性之间做出折衷与抉择； 模型的渐进性模型的渐进性：复杂的实际问题的建模通常不可能一次完成，：复杂的实际问题的建模通常不可能一次完成，包括由简到繁，也包括由繁到简，以获得满意的数学模型包括由简到繁，也包括由繁到简，以获得满意的数学模型； 模型的可转移性模型的可转移性：模型是现实对象的抽象化、理想化的产物，：模型是现实对象的抽象化、理想化的产物，不为对象所属领域所独有，具有应用的极端广泛性；不为对象所属领域所独有，具有应用的极端广泛性； 模型的非预测性模型的非预测性：建模本身是事先没有答案的，在建

4、模过程：建模本身是事先没有答案的，在建模过程中有时会伴随新的数学方法或数学概念产生；中有时会伴随新的数学方法或数学概念产生； 除此之外，建模还具有条理性、局限性等。除此之外，建模还具有条理性、局限性等。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学1 1依据对实际问题的了解程度可以分为依据对实际问题的了解程度可以分为：白箱模型、灰箱：白箱模型、灰箱模型和黑箱模型。模型和黑箱模型。2 2依据模型中的变量特征可以分为依据模型中的变量特征可以分为：连续型模型、离散型：连续型模型、离散型模型或者确定性模型、随机性模型等。模型或者确定性模型、随机性模型等。3 3依据建模中所用的数学

5、方法可以分为依据建模中所用的数学方法可以分为：初等数学模型、：初等数学模型、微分方程模型、优化模型、统计分析模型、控制论模型等。微分方程模型、优化模型、统计分析模型、控制论模型等。还有其他的不同分类方法。还有其他的不同分类方法。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模的步骤：数学建模的步骤：1 1模型的准备模型的准备：了解问题的实际背景，明确建模的目的，搜：了解问题的实际背景，明确建模的目的，搜集建模所需的各种信息（如现象、数据等），弄清研究对象的集建模所需的各种信息（如现象、数据等），弄清研究对象的特征、机理等，由此初步确定用哪一类模型；特征、机理等，由此

6、初步确定用哪一类模型； 2 2问题的假设问题的假设：由于实际问题比较复杂，直接建模比较困难：由于实际问题比较复杂，直接建模比较困难，所以合理的假设可以简化建模的难度，但也不能过于简化，否，所以合理的假设可以简化建模的难度，但也不能过于简化，否则所得到的模型与实际问题会有比较大的差异。合理的假设要符则所得到的模型与实际问题会有比较大的差异。合理的假设要符合客观事实，又要有一定的依据。由于考虑问题的视点不同，所合客观事实，又要有一定的依据。由于考虑问题的视点不同，所作的简化假设不同，因而对同一问题会得到不同的数学模型。作作的简化假设不同，因而对同一问题会得到不同的数学模型。作假设的依据，首先是出于

7、对问题内在规律的认识，其次是来自对假设的依据，首先是出于对问题内在规律的认识，其次是来自对数据、现象的分析，或者是二者的结合。数据、现象的分析，或者是二者的结合。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 3 3模型的建立模型的建立：主要说明建模的思路与依据。充分表述自己：主要说明建模的思路与依据。充分表述自己的思路与所用到的依据，表述要简洁、清晰，论据要充分，推的思路与所用到的依据，表述要简洁、清晰，论据要充分，推理及运算要正确。理及运算要正确。 4 4模型的求解模型的求解：对于已经建立的模型给出一个正确的解答，：对于已经建立的模型给出一个正确的解答，应用所学过的知

8、识，也可以借助于计算机，并应用相应的数学应用所学过的知识，也可以借助于计算机，并应用相应的数学软件等。软件等。 5 5模型分析模型分析：对模型求解进行数学上的分析，主要是进行误：对模型求解进行数学上的分析，主要是进行误差分析、稳定性分析及灵敏性分析等。差分析、稳定性分析及灵敏性分析等。6 6模型的检验与推广模型的检验与推广：把数学上的分析结果：把数学上的分析结果“翻译翻译”回到实回到实际问题，并用实际的现象、数据与之比较，检验模型的合理性际问题，并用实际的现象、数据与之比较，检验模型的合理性和适用性。这是建模成败的关键。模型检验的结果如果不符合和适用性。这是建模成败的关键。模型检验的结果如果不

9、符合或者部分不符合实际，应该对模型简化假设部分再做进一步的或者部分不符合实际，应该对模型简化假设部分再做进一步的修改、补充，重新建立模型。针对已经建立的模型指出其优点修改、补充，重新建立模型。针对已经建立的模型指出其优点与不足，同时将模型改进、推广以适应更大的应用范围。与不足，同时将模型改进、推广以适应更大的应用范围。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验二、二、 MatlabMatlabplotplot函数格式函数格式西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验数学建模与

10、数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学例3 图数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学3. 3. 图形标注说明图形标注说明 为了对图形进行清晰的说明，可以使用为了对图形进行清晰的说明，可以使用MATLABMATLAB标注函数将标标注函数将标题、坐标轴标记、网格线及文字注释加注到图形上。相关函数如题、坐标轴标记、网格线及文字注释加注到图形上。相关函数如下表所示。下表所示。其中，其中，axisaxis的用法还有：的用法还有： axis(xmin xmax ymin ymax) axi

11、s(xmin xmax ymin ymax) 用行向量中给出的用行向量中给出的值设定坐标轴的最大和最小值。值设定坐标轴的最大和最小值。 如如axis (-2 2 0 5)axis (-2 2 0 5)数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 axis(equal) axis(equal) 将两坐标轴设为相等将两坐标轴设为相等 axis on(off) axis on(off) 显示和关闭坐标轴的标记、标志显示和关闭坐标轴的标记、标志 axis auto axis auto 将坐标轴设置返回自动缺省值将坐标轴设置返回自动缺省值数学建模与数学实验数学建模与数学实验西安建

12、筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学plot3 函数说明函数说明 plot3是函数是函数plot的三维扩展。使用格式与的三维扩展。使用格式与plot相似，二维图形相似，二维图形的所有基本特性对三维图形全都适用，只是增加了一个维数而的所有基本特性对三维图形全都适用，只是增加了一个维数而已。已。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学定义三维坐标轴大小定义三维坐标轴大小axis(

13、xmin xmax ymin ymax zmin zmax ) grid on(off) 绘制三维网格绘制三维网格 text(x,y,z,string) 三维图形标注三维图形标注 子图和多窗口也可以用到三维图形中子图和多窗口也可以用到三维图形中例例1：用用plot3(x,y,z)格式，绘制参数方程格式，绘制参数方程 x=sin(t)，y=cos(t)，z=t。2. plot3 函数举例函数举例数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 函数函数meshgrid将给定的区域按一定的方式划分成平面网格，将

14、给定的区域按一定的方式划分成平面网格，该平面网格可以用来绘制三维曲面给定该平面网格可以用来绘制三维曲面给定x轴的范围，给定轴的范围，给定y轴的范轴的范围，生成围，生成xoy平面内的该区域划分平面网格，给出网格点坐标矩平面内的该区域划分平面网格，给出网格点坐标矩阵（同维），便于继续计算阵（同维），便于继续计算xoy平面内的点对应的平面内的点对应的z坐标。坐标。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学MATLAB还提供了不少特殊的三维图

15、形函数，如下表所示。还提供了不少特殊的三维图形函数，如下表所示。西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验三、三、 数学规划模型数学规划模型线性规划的一般形式如下：线性规划的一般形式如下：(1.3)数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学可以用向量矩阵形式表示上述优化问题。可以用向量矩阵形式表示上述优化问题。于是上述优化问题可表示成如下形式：于是上述优化问题可表示成如下形式：数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 在在MATLAB优化工具箱中，求解线性规划的命令为优化工具箱中，求解线性规划的命令为linprog

16、。对于形如。对于形如(1.4)形式的线性规划，形式的线性规划，linprog的调用的调用格式为：格式为：x,fval=linprog(c,A,B,Aeq,Beq,lb,ub)其中输出变量其中输出变量x为最优解，为最优解，fval为最优值。为最优值。 数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 例例6. 某车间有甲、乙两台机床，可用于加工三种工件。假定某车间有甲、乙两台机床，可用于加工三种工件。假定这两台车床的可用台时数分别为这两台车床的可用台时数分别为800和和900，三种工件的数量分别，三种工件的数

17、量分别为为400、600和和500，且已知用三种不同车床加工单位数量不同工，且已知用三种不同车床加工单位数量不同工件所需的台时数和加工费用如下表。问怎样分配车床的加工任务件所需的台时数和加工费用如下表。问怎样分配车床的加工任务，才能既满足加工工件的要求，又使加工费用最低？，才能既满足加工工件的要求，又使加工费用最低？数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学MATLAB求解程序如下：求解程序如下： c = 13 9 10 11 12 8; A = 0.4 1.1 1 0 0 0; 0 0 0 0.5

18、1.2 1.3; b = 800; 900; Aeq=eye(3) eye(3);beq=400 600 500; vb = zeros(6,1); x,fval = linprog(c,A,b,Aeq,beq,vb)数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 在最优化模型中，若目标函数或约束条件中至少有一个为非在最优化模型中，若目标函数或约束条件中至少有一个为非线性的，则称这类模型为非线性规划问题。非线性规划的一般形线性的，则称这类模型为非线性规划问题。非线性规划的一般形式如下：式如下：数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 求解无约束

19、优化问题常用的求解无约束优化问题常用的MATLAB函数是函数是fminunc，调用，调用格式如下格式如下: x,fval,exitflag=fminunc(FUN, x0)其中其中FUN为目标函数，为目标函数，x0为初始值，为初始值，x为最优解，为最优解，fval为最优值为最优值，exitflag为算法终止标志：若为算法终止标志：若exitflag=1，则输出的最优结果可，则输出的最优结果可靠；若靠；若exitflag=0，则输出的最优结果不可靠。，则输出的最优结果不可靠。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学x=linspace(0,20);y=sin(x)./

20、x;plot(x,y);xlabel(x);ylabel(y);title(sinc(x);grid onsinc曲线如右图所示曲线如右图所示：数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学下面建立下面建立sinc函数文件：函数文件： function y=sinc(x) y=sin(x)./x;初值初值x0分别取分别取2和和8，求最小值的命令为，求最小值的命令为 x1 f1 e1=fminunc(sinc,2); x2 f2 e2=fminunc(sinc,8);输出最优解输出最优解x1 =4.4934, x2 =10.9041, 最优值为最优值为f1 =-0.2172

21、, f2 =-0.0913。由此可以看出：。由此可以看出：fminunc只能寻找到初值只能寻找到初值x0附近的局部附近的局部最优解。最优解。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 将非线性规划将非线性规划(1.5)的约束按线性和非线分开表示，则它可的约束按线性和非线分开表示，则它可表示如下的标准型：表示如下的标准型：min( )s.t.( )0( )0eqeqfcceqAbAblbubxxxxxx数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 用用MATLAB求解上述问题，基本步骤分三步：求解上述问题，基本步骤分三步： 1 首先建立首先建立M

22、文件文件fun.m,用来定义目标函数用来定义目标函数F(X): function f=fun(X); f=F(X);02. 若约束条件中有非线性约束若约束条件中有非线性约束:C(X) 0或或Ceq(X)=0,则建立则建立M文件文件nonlcon.m定义函数定义函数C(X)与与Ceq(X): function C,Ceq=nonlcon(X) C=; Ceq= ;数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学3 建立主程序建立主程序.求解非线性规划的函数是求解非线性规划的函数是fmincon,命令的基本格式命令的基本格式如下：如下： (1) x=fmincon(fun,X0

23、,A,b) (2) x=fmincon(fun,X0,A,b,Aeq,beq) (3) x=fmincon(fun,X0,A,b, Aeq,beq,lb,lu) (4) x=fmincon(fun,X0,A,b,Aeq,beq, lb,lu,nonlcon) (5)x=fmincon(fun,X0,A,b,Aeq,beq, lb,lu,nonlcon,options) (6) x,fval= fmincon() (7) x,fval,exitflag= fmincon() (8)x,fval,exitflag,output= fmincon()输出极值点M文件迭代的初值变量上下限参数说明数学建

24、模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学因此，可用因此，可用fmincon函数求解函数求解(1.6)，讲义常用调用格式为：，讲义常用调用格式为： x,fval,exitflag = fmincon(fun,x0,A,b,Aeq,beq,lb,ub,nonlcon)数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学先将目标函数( )f x表示为表示为M文件：文件： function f = fun2(x) f = -x(1)*x(2)*x(3); x0 = 10;10;10; A=-1 -2 -2;1 2 2;b=0;72; x,f = fmincon(f

25、un2,x0,A,b)结果为结果为x =24.0000 12.0000 12.0000，f =-3.4560e+003。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学 多目标优化是在给定的约束范围内求多个目标的最值，一般多目标优化是在给定的约束范围内求多个目标的最值，一般形式如下：形式如下：12min( )( ),( ),.,( )s.t.( )0,1,2,.,( )0,1,2,.,nijFfffgimhjlxxxxxx数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验

26、数学建模与数学实验西安建筑科技大学西安建筑科技大学1212312412min3s.t. 333044162120,1,2,3,4ixxxxxxxxxxxi(1) 求解程序为：求解程序为：c=-3 -1 0 0;A=2 -1 0 0;b=12;Aeq=3 3 1 0;4 -4 0 1;beq=30;16;lb=zeros(4,1);x=linprog(c,A,b,Aeq,beq,lb)课后题：课后题：数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学1221212212121212min(42421)s.t.01.50100 xexxx xxxxx xxxx x解：解：对目标函

27、数建立对目标函数建立M文件：文件：function f=ex2(x)f=exp(x(1)*(4*x(1)2+2*x(2)2+4*x(1)*x(2)+2*x(2)+1);对非线性约束建立对非线性约束建立M文件：文件：function c ceq=ex2con(x)c=1.5+x(1)*x(2)-x(1)-x(2);-x(1)*x(2)-10;ceq=;主程序为主程序为 Aeq=1 1;beq=0;x0=1 -1;x fval exitflag=fmincon(ex2,x0,Aeq,beq,ex2con)数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学2、营养学家指出、营养学家

28、指出，成人良好的日常饮食应至少提供，成人良好的日常饮食应至少提供0.075kg的碳的碳水化合物，水化合物，0.06kg的蛋白质，的蛋白质，0.06kg的脂肪。的脂肪。1kg食物食物A含有含有0.105kg碳水化合物，碳水化合物，0.07kg的蛋白质，的蛋白质，0.14kg的脂肪，花的脂肪，花 费费28元；而元；而1kg食物食物B含有含有0.105kg碳水化合物，碳水化合物，0.14kg的蛋白质，的蛋白质，0.07kg的脂肪，花费的脂肪，花费21元。为了满足营养学家指出的日常饮食要元。为了满足营养学家指出的日常饮食要求，同时使花费最低，需要同时食用食物求，同时使花费最低，需要同时食用食物A和食物

29、和食物B多少多少kg？ 数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学1212121212min 2821s.t. 0.1050.1050.0750.070.140.060.140.070.060,0 xxxxxxxxxx求解程序为求解程序为c=28 21; A=-0.105 0.105;0.07 0.14;0.14 0.07;b=-0.075;0.06;0.06;x=linprog(c, A,b,0;0,);最优结果为：最优结果为：x=0.1429 0.5714.数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学3. 某公司生产甲、乙、丙三种产品。今已

30、知上述三种产品的单位某公司生产甲、乙、丙三种产品。今已知上述三种产品的单位产品原材料消耗定额分别为产品原材料消耗定额分别为4kg、4kg和和5kg；三种产品的单位产；三种产品的单位产品所需工时分别为品所需工时分别为6台时、台时、3台时和台时和6台时。由于生产该三种产品台时。由于生产该三种产品的原材料和工时的供应有一定限制，原材料每天只能供应的原材料和工时的供应有一定限制，原材料每天只能供应180kg，工时每天只有工时每天只有150台时。另外，三种产品的利润分别为台时。另外，三种产品的利润分别为400元每件、元每件、250元每件和元每件和300元每件。试建立能获得最大利润的优化模型，并元每件。试

31、建立能获得最大利润的优化模型，并进行求解。进行求解。数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学123123123123max 400250300s.t. 4451806361500,0,0 xxxxxxxxxxxx求解程序为求解程序为c=-400 250 300;A=4 4 5;6 3 6;b=180;150;x f=linprog(c,A,b,0 0 0,); -f西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 MATLAB优化工具箱提供了求解优化工具箱提供了求解0-1线性规划的函数命令线性规划的函数命令bintprog。考虑如下形式的。考虑如

32、下形式的0-1线性规划：线性规划：Tmins.t.1,2,.,eqeqjzc xAxbA xbxjn0,1， (3)bintprog调用格式为：调用格式为： x,z = bintprog(c, A, b, Aeq, beq)四、四、 整数规划模型整数规划模型西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验例例4. 求解求解12341234341224min9564s.t. 635291000,1,1,2,3,4izxxxxxxxxxxxxxxxi 解解：c = -9 5 6 4; A = 6 3 5 2; 0 0 1 1; -1 0 1 0; 0 -1 0 1; b = 9

33、 1 0 0; x z = bintprog(c,A,b)西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 下面给出使用分支界定法求解下面给出使用分支界定法求解整数线性规划的整数线性规划的MATLAB函函数程序。数程序。 function x,fval,status = intprog(c,A,b,Aeq,beq,lb,ub) global I e e= 0.00001;I = 1:length(c) if nargin 7, ub = ; if nargin 6, lb = ; if nargin 5, beq = ; if nargin 4, Aeq = ; end,

34、end, end, end西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验options = optimset(display,off);x0,fval0,exitflag = linprog(c,A,b,Aeq,beq,lb,ub,options);if exitflag 0 disp(无整数解无整数解); x = x0;fval = fval0; status = exitflag; return;else bound = inf; x,fval,status = branchbound (c,A,b,x0,fval0,bound,Aeq,beq,lb,ub);end西

35、安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验%以下为子函数以下为子函数function x_new,fval_new,status,bound_new = branchbound(c,A,b,x,fval,bound,Aeq,beq,lb,ub)global I eoptions = optimset(display,off);x0,fval0,status0=linprog(c,A,b,Aeq,beq,lb,ub,options);if status0 = bound x_new = x; fval_new = fval; bound_new = bound; stat

36、us = status0; return;end西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验int_index = find(abs(x0(I) - round(x0(I) e);if isempty(int_index) x_new(I) = round(x0(I);fval_new = fval0; bound_new = fval0;status = 1; return;endn = I(int_index(1);addA = zeros(1,length(c);addA(n) = 1;A = A;addA;b = b;floor(x(n); x1,fval1,s

37、tatus1,bound1= branchbound(c,A,b,x0,fval0,bound,Aeq,beq,lb,ub);西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验A(end,:) = ;b(end,:) = ;status = status1;if status1 0 & bound1 0 & bound2 bound status = status2; x_new = x2;fval_new = fval2; bound_new = bound2;end西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验例例5. 求解求解1221212min3

38、s.t.3.132234285,zxxxxxx x为非负整数解：clear c=1 3; A=0 -1;-22 -34;b=-3.13; 285; lb=0;0; x,fval,status = intprog(c,A,b,lb);x,fval,status = intprog(c,A,b,Aeq,beq,lb,ub) 西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验解：c=3 -2 5; A=1 2 -1;1 4 -1;1 1 0;0 4 1; b=2 4 3 6; x z=bintprog(c,A,b)课后题：课后题：西安建筑科技大学西安建筑科技大学数学建模与数学实验数

39、学建模与数学实验解：解：clearc=-1 1 -4;A=1 1 2;1 1 -1;-1 1 1;b=9 2 4;lb=0 0 0;x,fval,status = intprog(c,A,b,lb);fval=-fval;西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验2. 有一份中文说明书，需翻译成英、日、德、俄、法五种语言。有一份中文说明书，需翻译成英、日、德、俄、法五种语言。现有甲、乙、丙、丁、戊五个人，他们用各种语言翻译所用时现有甲、乙、丙、丁、戊五个人，他们用各种语言翻译所用时间如下表，问如何指派时间最少间如下表，问如何指派时间最少 ?西安建筑科技大学西安建筑科技

40、大学数学建模与数学实验数学建模与数学实验该指派问题的优化模型为该指派问题的优化模型为西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验123456789101112131415161718192021222324251611162127121722381318234914192min3868987441022226109739375510s.t.111zxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx45101520251234567891011121314151617181920212223242511111110,1,1,2,25

41、ixxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxiL西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验程序：程序：Data=3 8 2 10 3;8 7 2 9 7;6 4 2 7 5; 8 4 2 3 5;9 10 6 9 10;c=Data(:);Aeq=zeros(10,25);for i=1:5Aeq(i,i:5:25)=1;Aeq(i+5,1:5+(i-1)*5)=1;endbeq=ones(10,1);x fval exitflag=bintprog(c,Aeq,beq);reshape(x,5,5)西安建筑科技大学西安建筑科技大学数学建模与数学实

42、验数学建模与数学实验微分方程的解析解微分方程的解析解 在在MATLAB中，求解微分方程中，求解微分方程(组组)的函数是的函数是desolve，调用，调用格式如下：格式如下： r = dsolve(eq1,eq2,., cond1,cond2,. , v)或或 r = dsolve(eq1,eq2,.,cond1,cond2,. , v)微分方程模型微分方程模型西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 解：命令为解：命令为 r=dsolve(Dy=a*x,x) 结果为结果为r=1/2*a*x2+C1。解：命令为解：命令为 r=dsolve(Dx=x+sin(t),x(

43、0)=1,t) 结果为结果为r =-1/2*cos(t)-1/2*sin(t)+3/2*exp(t)。 其中其中eq1,eq2,用来表示常微分方程用来表示常微分方程(组组)，cond1, cond2,表表示初始或边界条件，示初始或边界条件，v表示自变量，缺省时默认为表示自变量，缺省时默认为t。西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 解：命令为解：命令为 r=dsolve(D2y-Dy2/y=0,y(0)=1,Dy(0)=2, x) 结果为结果为r =exp(2*x)。解：命令为解：命令为x,y=dsolve(Dx+5*x+y=exp(t),Dy-x-3*y=exp

44、(2*t),t)西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 解含初始条件的常微分方程的解含初始条件的常微分方程的MATLAB格式格式 t,Y = solver(odefun,tspan,y0)其中，其中，odefun表示以函数文件存储的微分方程表示以函数文件存储的微分方程(不含初始条件不含初始条件)，tspan为二维向量，用来表示区间的起点和终点，为二维向量，用来表示区间的起点和终点，y0表示自变量表示自变量取取tspan区间起点时的初始条件。区间起点时的初始条件。solver指指ode45、ode23、 ode113、ode15s、ode23s、ode23t和和od

45、e23tb七种求解方法中的一种。七种求解方法中的一种。西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 例例8. 求下列微分方程求下列微分方程1232133121230.51(0)0,(0)1,(0)1yy yyy yyy yyyy 解：先建立函数文件来表示微分方程：解：先建立函数文件来表示微分方程： function dy = rigid(t,y) dy = zeros(3,1); % dy应为列向量应为列向量 dy(1) = y(2) * y(3); dy(2) = -y(1) * y(3); dy(3

46、) = -0.51 * y(1) * y(2);西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验再在脚本函数或命令窗口中输入再在脚本函数或命令窗口中输入：t,Y = ode45(rigid,0 12,0 1 1);plot(t,Y(:,1),ro-,t,Y(:,2),g*:,t,Y(:,3),bx-)xlabel(t,FontSize,15)legend(y_1,y_2,y_3,3)运行后得运行后得右右图图.024681012-1-0.8-0.6-0.4-0.200.20.40.60.81ty1y2y3西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 例

47、例9. 求下列微分方程求下列微分方程2221000(1)0(0)0,(0)1d xdxxxdtdtxx在在t0,3000上的数值解，并绘出曲线。上的数值解，并绘出曲线。 解：先将上述高阶微分方程转化为下列微分方程组：解：先将上述高阶微分方程转化为下列微分方程组：1222121121000(1)(0)0,(0)1yyyyyyyy 西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验建立函数文件来表示微分方程：建立函数文件来表示微分方程：function dy = vdp1000(t,y)dy = zeros(2,1); % dy应为列向量应为列向量dy(1) = y(2);dy(

48、2) = 1000*(1 - y(1)2)*y(2) - y(1);在脚本函数或命令窗口中输入在脚本函数或命令窗口中输入：t,Y = ode15s(vdp1000,0 3000,0 1); % ode45运行时间较长运行时间较长plot(t,Y(:,1),r.-)xlabel(t,FontSize,15)运行结果如运行结果如右图右图。050010001500200025003000-2-1.5-1-0.500.511.522.5t西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验微分方程建模案例微分方程建模案例人口预测模型人口预测模型 在当前世界，人口问题已成为人们所面对的最

49、为严峻的问在当前世界，人口问题已成为人们所面对的最为严峻的问题之一。较大的人口基数给自然生态、人类健康和人口素质带题之一。较大的人口基数给自然生态、人类健康和人口素质带来的严重的影响。控制和预测人口数量是当前人口问题的核心来的严重的影响。控制和预测人口数量是当前人口问题的核心问题。表问题。表1给出了陕西省给出了陕西省1955年至年至2010年每间隔年每间隔5年的人口数据年的人口数据，根据这些数据建立人口统计数据模型，并预测，根据这些数据建立人口统计数据模型，并预测2015年和年和2020年陕西省的总人口。年陕西省的总人口。西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安

50、建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验上述线性方程组的最小二乘解为：上述线性方程组的最小二乘解为：从而从而 西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验r下面给出了求解下面给出了求解r的的MATLAB程序：程序：pop=1711.8 1954.2 2144.3 2427 2692.1 2831.4 3001.7 3316 35143644 3720 3735;u=mean(log(pop(2:end)/pop(1)./(1:length(pop)-1);r=exp(u)-1;绘

51、制实际人口与预测人口对比图、相对误差图的程序为：绘制实际人口与预测人口对比图、相对误差图的程序为：figure(1)plot(0:11*5+1955,pop,ro-)popp=zeros(size(pop);popp(1)=pop(1);for i=2:length(pop) popp(i)=popp(i-1)*(1+r);end西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验hold onplot(0:11*5+1955,popp,gs:)xlabel(年份年份),ylabel(人口人口 (单位：万单位：万)legend(实际人口实际人口,预测人口预测人口,2)axis(

52、1955 2010 1500 5500)grid onfigure(2)plot(0:11*5+1955,100*abs(popp-pop)./pop,ro-)xlabel(年份年份)ylabel(相对误差相对误差 (%)axis(1955 2010 0 40)grid on西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验结果如图结果如图196019701980199020002010150020002500300035004000450050005500年份人口 (单位：万)实际人口预测人口196019701980199020002010010203040年份相对误差 (

53、%)图图7 几何模型的实际与预测人口几何模型的实际与预测人口 图图8 几何模型的相对误差几何模型的相对误差西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验 从图从图7和图和图8可以看出，在前期，相对误差较小；在后期相可以看出，在前期，相对误差较小；在后期相对误差非常大，对误差非常大， 2010年的预测人口的相对误差高达年的预测人口的相对误差高达35.73%。因。因此，此模型不适合来预测此，此模型不适合来预测2015年和年和2020年的人口。年的人口。西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验即即因为因为所以可建立含初始条件的微分方程所以可建立含初始

54、条件的微分方程西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验求解上述微分方程得：求解上述微分方程得： 上述模型即为指数增长模型，也称为马尔萨斯人口模型，它上述模型即为指数增长模型，也称为马尔萨斯人口模型，它最早由英国人口学家马尔萨斯在最早由英国人口学家马尔萨斯在1798年提出。为了根据拟合方法年提出。为了根据拟合方法求解参数求解参数r，对上述公式两边取对数：，对上述公式两边取对数：西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验MATLAB求解程序为：求解程序为：t=0:11;y=log(pop);pop0=pop(1);r=t(y-log(pop0)

55、;下面分析指数增长模型的性能，程序如下：下面分析指数增长模型的性能，程序如下：tt=1955:5:2010;plot(tt,pop,o-,tt,pop0*exp(r*t),gs-,MarkerSize,10)xlabel(年份年份) ylabel(人口人口 (单位：万单位：万)legend(实际人口实际人口,预测人口预测人口,2)axis(1955 2010 1500 4500)grid on西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验figure(2)plot(tt,100*abs(pop-pop0*exp(r*t)./ pop,ro-,MarkerSize,10)x

56、label(年份年份) ylabel(相对误差相对误差 (%)axis(1955 2010 0 20)grid on 实验结果如图实验结果如图9和图和图10所示。从图所示。从图9可以看出：在多数时可以看出：在多数时刻，预测模型都有稍大的误差；从图刻，预测模型都有稍大的误差；从图10可以看出：指数模型可以看出：指数模型的相对误差相对稳定，最大值为的相对误差相对稳定，最大值为16.35%，比几何模型的相对，比几何模型的相对误差低很多。误差低很多。西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验1960197019801990200020101500200025003000350

57、040004500年份人口 (单位：万)实际人口预测人口19601970198019902000201005101520年份相对误差 (%)图图9 指数模型的实际与预测人口指数模型的实际与预测人口图图10 指数模型的相对误差指数模型的相对误差西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验于是建立下列微分方程：于是建立下列微分方程：西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验则则西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验xm/2dx/dtxOxmtx(t)x

58、mx0 xm/2O西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验附：附：x=0:0.01:1;dx=0.1*(1-x).*x;plot(x,dx,-)axis(0 1.1 0 0.04)hold onplot(0.5*ones(1,26),0:0.001:0.025,r:)text(0.5-0.03,-0.002,x_m/2)text(1-0.03,-0.002,x_m)text(1.1,-0.002,x)text(0-0.01,0-0.001,O)text(0-0.1,0.04,dx/dt)%=t=0:0.05:30;r=0.2237;x0=1711;xm=4310;x=xm./(1+(xm/x0-1)*exp(-r*t);plot(t,x)axis(0 32 0 4400)box offhold onplot(t,xm*ones(length(t),r:)西安建筑科技大学西安建筑科技大学数学建模与数学实验数学建模与数学实验function x=fun(beta,t)x0=1711.8;x=beta(1)./(1+(beta(1)/x0-1)*exp(-beta(2)*t);建立建立M