数学建模及数学实验课程设计
数学建模与数学实验
课程设计报告
学院数理学院专业数学与应用数学班级数学二班学号131412211 学生姓名任娟娟指导教师周树克
2014年6月
工厂最优生产计划模型
【摘要】本文针对工厂利用两种原料生产三种商品制定最优生产计划的问题，
建立优化问题的线性规划模型。在求解中得到了在不同生产计划下收益最优化的各产品的产量安排策略、最大收益，以及最优化生产计划的灵敏度分析。
对于问题一，通过合理的假设，首先根据题中所给的条件找出工厂收益的决定条件，利用线性规划列出目标函数MAX。由题目中所得，工厂原料及价格的约束条件下运用lingo软件算出最优生产条件下最大收益为1920元，其次是不同产品的产量。
对于问题二，灵敏度分析是研究当目标函数的费用系数和约束右端项在什么范围变化时，最优基保持不变。对产品结构优化制定及调整提供了有效的帮助。根据问题一所给的数据，运用lingo软件做灵敏度分析。
关键词：最优化线性规划灵敏度分析 LINGO
一、问题重述
某工厂利用两种原料甲、乙生产A1、A2、A3三种产品。如果每月可供
应的原料数量（单位：t ），每万件产品所需各种原料的数量及每万件产品
的价格如下表所示：
（1）试制定每月和最优生产计划，使得总收益最大；
（2）对求得的最优生产计划进行灵敏度分析。 二、模型假设
（1）在产品加工时不考虑排队等待加工的问题。
（2）假设工厂的原材料足够多，不会出现原材料断货的情况。
（3）忽略生产设备对产品加工的影响。
（4）假设工厂的原材料得到充分利用，无原材料浪费的现象。
三、符号说明
Xij （i=1,2,；j=1,2,3；）表示两种原料分别生产出产品的数量（万件）；
Max 为最大总收益；
A1，A2，A3为三种产品。
四、模型分析
问题一分析：对于问题一的目标是制定每月和最优生产计划，求其最大生产
效益。由题中所给的条件找出工厂收益的决定条件，利用线性规划列出目标函数MAX 。由题目中所得，工厂原料工厂原料及价格的约束，列出约束条件。
问题二分析：研究当目标函数的费用系数和约束右端项在什么范围变化时，
最优基保持不变。通过软件数据进行分析。
五、模型建立与求解
问题一的求解：
建立模型：
题目的目标是寻求总利益最大化，而利润为两种原料生产的六种产品所获得
的利润之和。
设Xij （i=1,2,；j=1,2,3；）表示两种原料分别生产出产品的数量（万件）
则目标函数:max=12(x11+x21）+5(x12+x22)+4(x13+x23)
原料
每万件产品所需原料（t ） 每月原料供应量（t ） A1 A2 A3 甲
4 3 1 180 乙
2 6
3 200 价格（万元/万
件）
12 5 4
约束条件：
1）原料供应：4x11+3x12+x13<=180;
2x21+6x22+3x23<=200
2）非负约束：x11,x12,x13,x21,x22,x23>=0
所以模型为：
max=12(x11+x21）+5(x12+x22)+4(x13+x23)
S.t
200x x 6x 2180x x 34x 232221131211<=++<=++ 0x >=ij (i=1,2;j=1,2,3且为整数)}
模型求解：
model :
max =12*x11+12*x21+5*x12+5*x22+4*x13+4*x23;
4*x11+3*x12+x13<=180;
2*x21+6*x22+3*x23<=200;
End
计算结果：
Global optimal solution found.
Objective value: 1920.000
Infeasibilities: 0.000000
Total solver iterations: 0
Variable Value Reduced
Cost
X11 0.000000 4.000000
X21 100.0000 0.000000
X12 0.000000 7.000000
X22 0.000000 31.00000
X13 180.0000 0.000000
X23 0.000000 14.00000
Row Slack or Surplus Dual
Price
1 1920.000
1.000000
2 0.000000 4.000000
3 0.000000 6.000000
结论：从数据表明，这个线性规划的最优解为
x11=0,x12=0,x13=180,x21=100,x22=0,x23=0 ,最优值为1920.即这个工厂的最
优生产计划为：用甲原料生产A1，A2，A3产品数量分别为0万件，0万件，180万
件；用乙原料生产A1,A2,A3产品数量分别为100万件，0万件，0万件。
问题二的求解：
用lingo软件对模型进行灵敏度分析的结果如下：
Ranges in which the basis is unchanged:
Objective Coefficient Ranges
Current Allowable Allowable
Variable Coefficient Increase Decrease
X11 12.00000 4.000000 INFINITY
X21 12.00000 INFINITY 9.333333
X12 5.000000 7.000000 INFINITY
X22 5.000000 31.00000 INFINITY
X13 4.000000 INFINITY 1.000000
X23 4.000000 14.00000 INFINITY
Righthand Side Ranges
Row Current Allowable Allowable
RHS Increase Decrease
2 180.0000 INFINITY 180.0000
3 200.0000 INFINITY 200.0000
显然可以看出：在最优值不变的条件下目标函数系数允许变化的范围：x11
的系数为（12,12+4）=（12,16）；x12的系数为（5,5+7）=（5,12）；x13的系数
为（4-1,4）=（3,4）；x21的系数为（12-9.333333,12）=（2.666667,12）；x22
的系数为（5,5+31）=（5,36）；x23的系数为（4,4+14）=(4,18)。同样看出约束右端的限制数没有发生变化。由于目标函数的系数并不影响约束条件，所以最优解保持不变。
六、模型的优缺点
模型的优点：
（1）模型的适用性好，线性规划性比较好，能够随着市场的变化而做出相应的变动，从而得到更大的效益，具有更强的应用指导意义。
（2）模型的建立运用线性规划的方法，可理解性强，应用广泛。
（3）Lingo软件执行速度很快，易于输入，修改，求解，分析数学规划的问题。
模型的缺点：
（1）没有考虑到机床维修的费用对工厂总体效益的影响，与实际情况有出入。
（2）模型比较单一，并没有用更好的办法去进行相应的检验其最大收益，及最优生产计划。
七、模型的推广
本文的模型是一个典型的线性规划的模型，用来求解最大或最小目标函数极值问题。此问题有很多的推广应用价值。优化问题可以说是人们应用科学、工程设计、商业贸易等领域中常遇到的一类问题。这种数学建模的方法来处理优化问题，即建立和求解所谓的优化模型。虽然，由于建模时要适当做出简化，可能是结果不一定完全可行或达到实际上的困扰，但是它基于客观规律和数据，模型的建立与求解并不需要耗费太多的时间。如果在建模的基础上在赋予其现实的意义，就可以期望得到实际问题的一个圆满的结果。
八、参考文献
[1]赵静，但琦，数学建模与数学实验，北京，高等教育版社，2008.1
[2]姜启源，谢金星，叶俊，数学模型 [M]，北京：高等教育出版社，2003
数学建模-大学生就业问题
2010-2011第二学期 数学建模课程设计 2011年6月27日－7月1日 题目大学生就业问题 第 11 组组员1 组员2 组员3 组员4 姓名 学号 0808060217 0808060218 0808060219 0808060220 专业信计0802 信计0802 信计0802 信计0802 成绩
论文摘要 本文讨论了在新的形势下大学生的就业问题。20世纪90年代以来，我国出现了一种前所未有的现象，有着“天之骄子”美誉的大学生也开始面临失业问题。大学生就业难问题已受到普遍关注。大学生毕业失业群体正在不断扩大，已成为我国扩大社会就业，构建和谐稳定社会的急需解决的社会问题。 本文针对我国现有的国情，综合考虑了高校毕业生的就业率和高校招生规模的扩大之间的关系，建立了定量分析的微分方程模型，随后又建立了了离散正交曲线拟合模型对得出的结果进行了检验，并分析模型得出的结果得合理性。最终得到生源数量与失业率之间的拟合多项式和拟合曲线，并预测出了未来高校招生规模的变化趋势。 在找到大学生失业规律以后，本文还具体的对毕业生的性别、出生地对失业的影响做出了定量分析。 关键词：大学生就业微分方程模型多项式曲线拟合MATLAB软件 1、问题重述 大学生就业问题：如果我们将每年毕业的大学生中既没有找到工作又没有继续深造的情况视为失业，就可以用失业率来反映大学生就业的状况。下面的表中给出了某城市的大学生失业数占城市总失业人数的比率，比率的计算是按照国际劳工组织的定义，对16岁以上失业人员进行统计的结果。 表 1
请建立相应的模型对大学生就业状况进行分析找出其中的规律并讨论下面两个问题： (1)、就业中是否存在性别歧视； (2)、学生的出生对就业是否有影响。 2、模型假设 2.1在本次研究中做出以下假设： (1)、假设毕业生求职时竞争是公平的； (2)、假设考研等继续深造的毕业生属于已就业人群； (3)、假设每个毕业生都有就业或者继续深造的意图 (4)、假设就业率和失业率之和为1； (5)、假设本文搜集的数据全部真实可靠； 2.2 在定量分析性别、出生地对失业的影响时还要做以下假设： (1)、假设毕业生就业情况只受性别、出生地等因素的影响； (2)、假设具有上述同等条件的毕业生间就业机会相同 (3)、假设附件中的数据信息均合理； 3、问题分析 3.1 对问题的分析 若要分析新失业群体产生的主要原因，并就其重要性给出各种因素的排序，就需要对搜集的数据进行整理，并进行系统的分析，划分为不同的体系和矛盾，然后我们考虑用Logistic模型分析。 为了得到新失业群体对高校招生生源的影响和预测未来高校招生规模的变
数学建模课程设计报告范本
数学建模课程设计 报告 1 2020年4月19日
数学建模课程设计 题目： 学院： 专业： 班级： 姓名： 学号： 指导教师： 实验日期： 2 2020年4月19日
摘要 本文针对葡萄酒的质量分析与评价问题，以置信区间、优势矩阵、逐步回归分析等方法和方差分析理论为基础，首先分别构建了以评酒员和样酒为组别的方差数据序列，经过进行双向显著性检验，接着经过置信区间法处理的数据进行了方差分析，并确定可信的评价组别。然后以评酒员感官评价为主、葡萄酒的理化指标为辅，采用回归分析、聚类分析、判别分析法建立葡萄分级模型，继而使用相关系数矩阵确立葡萄酒与葡萄理化指标中具有较大相关性的指标，实现对葡萄理化指标的初步筛选，进行等级划分。再利用逐步回归的方法拟合酿葡萄酒理化指标与葡萄理化指标间一对多的函数关系得出二者之间的联系。最后经过上文函数关系，同时提取对香气与口感评分相关度较大的芳香物质，建立芳香物质与葡萄酒质量的函数关系，论证葡萄和葡萄酒的理化指标只在一定程度上对葡萄酒的质量有影响。 关键字：双向显著性检验；方差分析；置信区间；聚类分析；标准化； 1 2020年4月19日
一、问题重述 确定葡萄酒质量时一般是经过聘请一批有资质的评酒员进行品评。每个评酒员在对葡萄酒进行品尝后对其分类指标打分，然后求和得到其总分，从而确定葡萄酒的质量。酿酒葡萄的好坏与所酿葡萄酒的质量有直接的关系，葡萄酒和酿酒葡萄检测的一级理化指标会在一定程度上反映葡萄酒和葡萄的质量。附件1给出了某一年份一些葡萄酒的评价结果，附件2和附件3分别给出了该年份这些葡萄酒的和酿酒葡萄的成分数据。请尝试建立数学模型讨论下列问题： 1. 分析附件1中两组评酒员的评价结果有无显著性差异，哪一组结果更可信？ 2. 根据酿酒葡萄的一级理化指标和葡萄酒的质量对这些酿酒葡萄进行分级。 3. 分析酿酒葡萄与葡萄酒的理化指标之间的联系。 4．分析酿酒葡萄和葡萄酒的一级理化指标对葡萄酒质量的影响，并论证能否用葡萄和葡萄酒的一级理化指标来评价葡萄酒的质 2 2020年4月19日
数学建模实验答案-概率模型
数学建模实验答案-概率模型
实验10 概率模型（2学时） （第9章 概率模型） 1.（验证）报童的诀窍p302~304, 323（习题2） 关于每天报纸购进量的优化模型： 已知b 为每份报纸的购进价，a 为零售价，c 为退回价（a > b > c ），每天报纸的需求量为r 份的概率是f (r )(r =0,1,2,…)。 求每天购进量n 份，使日平均收入，即 1 ()[()()()]()()()n r r n G n a b r b c n r f r a b nf r ∞ ==+=----+ -∑∑ 达到最大。 视r 为连续变量，f (r )转化为概率密度函数p (r )，则所求n *满足 * ()n a b p r dr a c -= -? 已知b =, a =1, c =，r 服从均值μ=500（份），均方差σ=50（份）的正态分布。报童每天应购进多少份报纸才能使平均收入最高，这个最高收入是多少 [提示：normpdf, normcdf] 要求：
(1) 在同一图形窗口内绘制10 ()()n y n p r dr =?和2()a b y n a c -= -的图形，观察其交点。 [提示] 22 ()2()r p r μσ-- = ，0 ()()()n n p r dr p r dr p r dr -∞ -∞ =-?? ? ☆(1) 运行程序并给出结果： (2) 求方程0()n a b p r dr a c -= -?的根n *（四舍五入取整），并求G (n *)。
mu=500;sigma=50; a=1; b=; c=; r=n+1; while (a-b)*n*normpdf(r,mu,sigma)>1e-6 r=r+1; end r=n+1:r; G=sum((a-b)*n*normpdf(r,mu,sigma)); r=0:n; G=G+sum(((a-b)*r-(b-c)*(n-r)).*normpdf(r,mu,sigma)) ☆(2) 运行程序并给出结果： 2.（编程）轧钢中的浪费p307~310 设要轧制长l=的成品钢材，由粗轧设备等因素决定的粗轧冷却后钢材长度的均方差σ=，问这时钢材长度的均值m应调整到多少使浪费最少。 平均每得到一根成品材所需钢材的长度为 () () m J m P m = 其中， 2 2 () 2 ()(), () 2 x m l P m p x dx p xσ πσ - - ∞ == ? 求m使J(m)达到最小。 等价于求方程 () () z z z λ ? Φ =- 的根z*。 其中：
大学生就业问题数学模型
重庆交通大学学生实验报告 实验课程名称数学模型课程设计 开课实验室数学实验室 学院 XXX级 XXX 专业 1 班 开课时间 2013 至 2014 学年第 2 学期设计题目大学生就业问题
2013 年 12月 大学生就业问题 摘要：近年来，我国高校毕业生数量逐年增多，加之当前金融危机的影响，毕业生的就业形势受到前所未有的挑战，甚至出现了所谓“毕业即失业”的说法。因此大学生毕业后能否顺利就业，已成为全社会普遍关注的热点问题。大学生就业难不仅有社会原因，也有大学生自身的原因。如何解决大学生就业难的问题不仅关系到大学生的切身利益，更关系到社会的和谐稳定，需要政府、企业、高校和大学生共同的努力。本文从大学生自身，企业和社会三个大方面方面进行了分析和论述，从而总结出相关的结论及解决大学生就业难题的可行方法。 关键词大学生就业 Matlab 数据拟合 一、问题重述 据中国媒体援引人力和社会保障部的最新统计数据，二零一零年全国高校毕业生为630万人，比去年的611万多19万人，加上往届未能就业的，需要就业的毕业生数量很大，高校毕业生就业形势十分严峻。 随着九十年代末大学扩招和教育产业化政策推行以来，大学生人数的增幅远远超过经济增长所需要的人才增长，大学生就业不难才是怪事，"毕业即失业"成为中国大学生的普遍现象。 尽管如此，中国教育部决定继续扩大全日制专业学位硕士研究生招生规模，努力培养更多高层次、应用型人才。表面上看，研究生扩招能提高大学生学历层次，可以缓解就业难。但是，如果不清理高等教育积弊，扩招研究生来应对就业难将是饮鸩止渴，使就业矛盾更加突出。 现在大学生就业难的问题，是由许多原因造成的，既有社会原因，也有历史原因。 请用数学建模的方法从以下几个侧面探讨大学生就业问题： （1）利用网上大学生就业统计数据建立大学生就业供需预测模型，利用所建模型对2012年就业形势进行预测； （2）分析影响大学生就业的主要因素，建立就业竞争力评价模型，利用所建模型评估你的竞争力；
数学建模实验报告
数学建模实验报告
一、实验目的 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 0 1 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 c = 1 1 0 1 0 1 1 1 4、解决方法（MATLAB程序代码）：
n=10;r=10;d=1000; a=0; for l=1:d m=full(sparse(randint(1,r,[1,n]),1:r,1,n,r)); c=zeros(n,1); for i=1:n for j=1:r if m(i,j)==1 c(j)=1; break; end continue; end end s=0; for x=1:n if c(x)==1 s=s+1; end continue; end a=a+s; end a/d 5、实验结果 ans = 6.5150 那么，当楼高11层，乘坐10人时，电梯需停次数的数学期望为6.5150。 （二）题目二 1、问题：某厂生产甲乙两种口味的饮料,每百箱甲饮料需用原料6 千克,工人10名,可获利10万元;每百箱乙饮料需用原料5千 克,工人20名,可获利9万元.今工厂共有原料60千克,工人 150名,又由于其他条件所限甲饮料产量不超过8百箱.问如何 安排生产计划,即两种饮料各生产多少使获利最大.进一步讨 论: 1)若投资0.8万元可增加原料1千克,问应否作这项投资. 2)若每百箱甲饮料获利可增加1万元,问应否改变生产计划. 2、问题分析 （1）题目中共有3个约束条件，分别来自原料量、工人数与甲饮料产量的限制。 （2）目标函数是求获利最大时的生产分配，应用MATLAB时要转换
数学模型课程设计一
课程设计名称： 设计一：MATLAB 软件入门 指导教师： 张莉 课程设计时数： 8 课程设计设备：安装了Matlab 、C ++软件的计算机 课程设计日期： 实验地点： 第五教学楼北902 课程设计目的： 1. 熟悉MA TLAB 软件的用户环境； 2. 了解MA TLAB 软件的一般目的命令； 3. 掌握MA TLAB 数组操作与运算函数； 4. 掌握MATLAB 软件的基本绘图命令； 4. 掌握MA TLAB 语言的几种循环、条件和开关选择结构。 课程设计准备： 1. 在开始本实验之前，请回顾相关内容； 2. 需要一台准备安装Windows XP Professional 操作系统和装有数学软件的计算机。 课程设计内容及要求 要求：设计过程必须包括问题的简要叙述、问题分析、实验程序及注释、实验数据及结果分析和实验结论几个主要部分。 1. 采用向量构造符得到向量[1,4,7,,31] 。 //a=[1:3:31] 2. 随机产生一向量x ，求向量x 的最大值。 // a=rand(1,6) max(a) 3. 利用列向量(1,2,3,,6)T 建立一个范德蒙矩阵A ，并利用位于矩阵A 的奇数行偶数列的元素建立一个新的矩阵B ，须保持这些元素的相对位置不变。 4. 按水平和竖直方向分别合并下述两个矩阵： 100234110,5670018910A B ????????==???????????? 5. 当100n =时，求1121n i y i ==-∑的值。 6. 一个三位整数各位数字的立方和等于该数本身则称该数为水仙花数。输出全部水仙花数。 7. 求[1000,2000]之间第一个被17整除的整数。 8. 用MATLAB 绘制两条曲线，[0,2]x π∈，以10 π为步长，一条是正弦曲线，一条是余弦曲线，线宽为6个象素，正弦曲线为绿色，余弦曲线为红色，线型分别为实线和虚线，并给所绘的两条曲线增添图例，分别为“正弦曲线”和“余弦曲线”。
环境数模课程设计说明书
2016《环境数学模型》课程设计说明书 1.题目 活性污泥系统生化反应器中底物降解与微生物增长数学模型的建立 2.实验方法与结果 2.1.实验方法 2.1.1.工艺流程与反应器 本设计采用的工艺流程如下图所示： 图2-1 活性污泥系统工艺流程图 本设计工艺采用活性污泥法处理污水，工艺的主要反应器包括生化反应器和沉淀池。污水通过蠕动泵恒速加到生化反应器中，反应器内活性污泥和污水在机械搅拌设备和鼓风曝气设备的共同作用下充分接触，并在氧气充足的条件下进行反应。经处理后，污泥混液通过管道自流到沉淀池中，在里面实现泥水分离。分离后的水通过溢流堰从周边排出，直接被排放到下水道系统，沉淀下来的污泥则通过回流泵，全部被抽回进行回流。 系统运行过程中，进出水流量、进水质量、污水的停留时间、生化反应器的容积、机械搅拌设备转轴转速、鼓风曝气装置的曝气风量气速、污泥回流量等参数在系统运行的过程中都保持不变。待系统持续运行一周稳定后再取样进行分析。 实验的进水为实验室配置的污水，污水分别以葡萄糖、尿素、磷酸二氢钾为碳源、氮源和磷源，其中C：N：P=100：40：1（浓度比），TOC含量为200mg/L。生化反应器内污泥混液的容量为12L，污水停留时间为6h。系统运行时间为两周，第一周是调适阶段，第二周取样测试，测得的数据作为建模的原始数据。 表2-1 污水中各营养物质的含量 2.1.2.取样方法
每隔24h取一次样，通过虹吸管取样。每次取样时，先取进水和出水水样用于测水体的COD指标，其中进水直接取配得的污水溶液，出水取沉淀池上清液。取得的水样过膜除去水中的悬浮固体和微生物，保存在5ml玻璃消解管中，并在4℃下冷藏保存。 取完用于测COD的水样后，全开污泥回流泵，将沉淀池中的污泥全部抽回生化反应器（由于实验装置的原因，沉淀池排泥管易堵，污泥易积聚在沉淀池中，为更准确测定活性污泥的增长情况，在此实验中将泥完全抽回后再测定），待搅拌均匀后，取5ml污泥混液于干净、衡重的坩埚中，待用于测污泥混液的SS。 2.1. 3.分析方法 本实验一共分析进出水COD和污泥混液SS两个指标。其中COD采用《水质快速消解分光光度法》（HJ/T 399-2007）方法进行分析，SS采用《水质悬浮物的测定重量法》（GB 11901-89）方法进行分析。 准确取2ml经过膜处理的水样于5mlcod消解管中，以重铬酸钾为氧化剂，硫酸银-浓硫酸为催化剂，硫酸汞为抗氯离子干扰剂，按一定比例与水样混合均匀。将消解管放在COD 消解仪中，在150℃条件下消解2h。待经消解的溶液冷却后，以空白样为参比液，在COD 分析仪上读出待测水样的COD值，记录数据。 将装在已衡重称重的坩埚中的污泥混液放在烘箱中，在105℃温度下烘3h以上，保证污泥中的水分被充分除去。坩埚冷却后衡重称重，记录干污泥的质量，求得活性污泥的SS。 实验过程的所有样品都设置两个平行样，最后结果取平行样的算术平均值。 2.2.实验结果 2.2.1.实验数据 实验测得数据如下表： 表2-2 活性污泥系统水质分析结果 2.2.2.数据分析
数学建模实验报告
数学建模实验报告 实验一计算课本251页A矩阵的最大特征根和最大特征向量 1 实验目的 通过Wolfram Mathematica软件计算下列A矩阵的最大特征根和最大特征向量。 2 实验过程 本实验运用了Wolfram Mathematica软件计算，计算的代码如下：
3 实验结果分析 从代码的运行结果，可以得到最大特征根为5.07293，最大特征向量为 {{0.262281},{0.474395},{0.0544921},{0.0985336},{0.110298}},实验结果 与标准答案符合。
实验二求解食饵-捕食者模型方程的数值解 1实验目的 通过Wolfram Mathematica或MATLAB软件求解下列习题。 一个生物系统中有食饵和捕食者两种种群，设食饵的数量为x（t），捕食者为y（t），它们满足的方程组为x’(t)=(r-ay)x,y’(t)=-(d-bx)y,称该系统为食饵-捕食者模型。当r=1，d=0.5，a=0.1，b=0.02时，求满足初始条件x（0）=25，y（0）=2的方程的数值解。 2 实验过程 实验的代码如下 Wolfram Mathematica源代码： Clear[x,y] sol=NDSolve[{x'[t] (1-0.1y[t])x[t],y'[t] 0.02x[t]y[t]-0.5y[t],x[0 ] 25,y[0] 2},{x[t],y[t]},{t,0,100}] x[t_]=x[t]/.sol y[t_]=y[t]/.sol g1=Plot[x[t],{t,0,20},PlotStyle->RGBColor[1,0,0],PlotRange->{0,11 0}] g2=Plot[y[t],{t,0,20},PlotStyle->RGBColor[0,1,0],PlotRange->{0,40 }] g3=Plot[{x[t],y[t]},{t,0,20},PlotStyle→{RGBColor[1,0,0],RGBColor[ 0,1,0]},PlotRange->{0,110}] matlab源代码 function [ t,x ]=f ts=0:0.1:15; x0=[25,2]; [t,x]=ode45('shier',ts,x0); End function xdot=shier(t,x)
数学建模实验答案初等模型
实验02 初等模型（4学时） （第2章初等模型） 1.（编程）光盘的数据容量p23~27 表1 3种光盘的基本数据 CAV光盘：恒定角速度的光盘。 CLV光盘：恒定线速度的光盘。 R2=58 mm, R1=22.5 mm，d, ρ见表1。
CLV光盘的信息总长度(mm) L CLV 22 21 () R R d π- ≈ CLV光盘的信息容量(MB) C CLV = ρL CLV / (10^6) CLV光盘的影像时间(min) T CLV = C CLV / (0.62×60) CAV光盘的信息总长度(mm) L CAV 2 2 2 R d π≈ CAV光盘的信息容量(MB) C CAV = ρL CAV / (10^6) CAV光盘的影像时间(min ) T CAV = C CAV / (0.62×60) 1.1（验证、编程）模型求解 要求： ①（验证）分别计算出LCLV, CCLV和TCLV三个3行1列的列向量，仍后输出结果，并与P26的表2（教材）比较。 程序如下：
②（编程）对于LCAV, CCAV和TCAV，编写类似①的程序，并运行，结果与P26的表3（教材）比较。 ★要求①的程序的运行结果： ★要求②的程序及其运行结果：
1.2（编程）结果分析 信道长度LCLV 的精确计算：21 2R CLV R L d π=? 模型给出的是近似值：2221() CLV R R L L d π-= ≈ 相对误差为：CLV L L L δ-= 要求：
①取R2=58 mm, R1=22.5 mm，d, ρ见表1（题1）。 分别计算出LCLV, L和delta三个3行1列的列向量，仍后将它组合起来输出一个3行3列的结果。 ②结果与P26的表2和P27（教材）的结果比较。 [提示] 定积分计算用quad、quadl或trapz函数，注意要分别取d的元素来计算。要用数组d参与计算，可用quadv（用help查看其用法）。 ★编写的程序和运行结果： 程序：
数学建模课程设计汇本参考模板
2015-2016第1学期数学建模课程设计题目：医疗保障基金额度的分配 ： 学号： 班级： 时间：
摘要 随着人们生活水平的提高及社会制度的发展,医疗保险事业显得越来越重要,各企业也随之越来越注重员工的福利措施，医疗保障基金额度的分配也成为了人们的关注热点。扩大医疗保障受益人口也是政府和企业面临的难题，因而根据历史统计数据，合理的构造出拟合曲线，分析拟合函数的拟合程度，从而为基金的调配以及各种分配方案做方向上的指导。 本文针对A,B两个公司关于医疗保障基金额度的合理分配问题，根据两公司从1980-2003年统计的医疗费用支出数据，科学地运用了MATLAB软件并基于最小二乘法则进行了多项式曲线拟合，成功建立了医疗保障基金额度的分配模型。最后，对不同阶数的多项式拟合曲线的拟合程度进行了残差分析，并输出相关结果，得出拟合程度与多项式阶数的关联。 此问题建立在收集了大量数据的基础上，以及利用了MATLAB编程拟合曲线，使问题更加简单，清晰。该模型经过适当的改造，可以推广到股票预测，市场销售额统计等相关领域。
关键字：matlab，最小二乘多项式拟合,阶数，残差分析 一．问题重述 某集团下设两个子公司：子公司A、子公司B。各子公司财务分别独立核算。每个子公司都实施了对雇员的医疗保障计划，由各子公司自行承担雇员的全部医疗费用。过去的统计数据表明，每个子公司的雇员人数以及每一年龄段的雇员比例，在各年度都保持相对稳定。各子公司各年度的医疗费用支出见下表（附录1）。 试利用多项式数据拟合，得到每个公司医疗费用变化函数，并绘出标出原始数据的拟合函数曲线。需给出三种不同阶数的多项式数据拟合，并分析拟合曲线与原始数据的拟合程度。 二．模型假设 1．假设A,B两公司在1980年底才发放医疗保障基金。
数学建模与数学实验报告
数学建模与数学实验报告 指导教师__郑克龙___ 成绩____________ 组员1：班级______________ 姓名______________ 学号_____________ 组员2：班级______________ 姓名______________ 学号______________ 实验1.(1)绘制函数cos(tan())y x π=的图像，将其程序及图形粘贴在此。 >> x=-pi:0.01:pi; >> y=cos(tan(pi*x)); >> plot(x,y) -4 -3 -2 -1 1 2 3 4 -1-0.8-0.6-0.4-0.200.20.40.60.8 1 （2）用surf,mesh 命令绘制曲面2 2 2z x y =+，将其程序及图形粘贴在此。（注：图形注意拖放，不要太大）（20分） >> [x,y]=meshgrid([-2:0.1:2]); >> z=2*x.^2+y.^2; >> surf(x,y,z)
-2 2 >> mesh(x,y,z) -2 2 实验2. 1、某校60名学生的一次考试成绩如下:
93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55 1)计算均值、标准差、极差、偏度、峰度，画出直方图；2)检验分布的正态性；3)若检验符合正态分布，估计正态分布的参数并检验参数. （20分） 1） >> a=[93 75 83 93 91 85 84 82 77 76 77 95 94 89 91 88 86 83 96 81 79 97 78 75 67 69 68 84 83 81 75 66 85 70 94 84 83 82 80 78 74 73 76 70 86 76 90 89 71 66 86 73 80 94 79 78 77 63 53 55]; >> pjz=mean(a) pjz = 80.1000 >> bzhc=std(a) bzhc = 9.7106 >> jc=max(a)-min(a) jc = 44 >> bar(a)
数学建模实验报告
matlab 试验报告 姓名 学号 班级 问题：.（插值） 在某海域测得一些点(x,y)处的水深z 由下表给出，船的吃水深度为5英尺，在矩形区域（75，200）*（-50，150）里的哪些地方船要避免进入。 问题的分析和假设： 分析:本题利用插值法求出水深小于5英尺的区域，利用题中所给的数据，可以求出通过空间各点的三维曲面。随后，求出水深小于5英尺的范围。 基本假设：1表中的统计数据均真实可靠。 2矩形区域外的海域不对矩形海域造成影响。 符号规定：x ―――表示海域的横向位置 y ―――表示海域的纵向位置 z ―――表示海域的深度 建模： 1.输入插值基点数据。 2．在矩形区域（75,200）×（－50,150）作二维插值，运用三次插值法。 3．作海底曲面图。 4．作出水深小于5的海域范围，即z=5的等高线。 x y z 129 140 103.5 88 185.5 195 105 7.5 141.5 23 147 22.5 137.5 85.5 4 8 6 8 6 8 8 x y z 157.5 107.5 77 81 162 162 117.5 -6.5 -81 3 56.5 -66.5 84 -33.5 9 9 8 8 9 4 9
求解的Matlab程序代码： x=[129 140 103.5 88 185.5 195 105.5 157.5 107.5 77 81 162 162 117.5]; y=[7.5 141.5 23 147 22.5 137.5 85.5 -6.5 -81 3 56.5 -66.5 84 -33.5]; z=[-4 -8 -6 -8 -6 -8 -8 -9 -9 -8 -8 -9 -4 -9]; cx=75:0.5:200; cy=-50:0.5:150; cz=griddata(x,y,z,cx,cy','cubic'); meshz(cx,cy,cz),rotate3d xlabel('X'),ylabel('Y'),zlabel('Z') %pause figure(2),contour(cx,cy,cz,[-5 -5]);grid hold on plot(x,y,'+') xlabel('X'),ylabel('Y') 计算结果与问题分析讨论： 运行结果： Figure1:海底曲面图：
数学建模与数学实验习题
数学建模与数学实验课程总结与练习内容总结 第一章 1.简述数学建模的一般步骤。 2.简述数学建模的分类方法。 3.简述数学模型与建模过程的特点。 第二章 4.抢渡长江模型的前3问。 5.补充的输油管道优化设计。 6.非线性方程（组）求近似根方法。 第三章 7.层次结构模型的构造。 8.成对比较矩阵的一致性分析。 第五章 9.曲线拟合法与最小二乘法。 10 分段插值法。 第六章 11 指数模型及LOGISTIC模型的求解与性质。 12.VOLTERRA模型在相平面上求解及周期平均值。 13 差分方程（组）的平衡点及稳定性。 14 一阶差分方程求解。 15 养老保险模型。
16 金融公司支付基金的流动。 17 LESLLIE 模型。 18 泛函极值的欧拉方法。 19 最短路问题的邻接矩阵。 20 最优化问题的一般数学描述。 21 马尔科夫过程的平衡点。 22 零件的预防性更换。 练习集锦 1. 在层次分析法建模中，我们介绍了成对比较矩阵概念，已知矩阵P 是成对比较矩阵 31/52a b P c d e f ?? ??=?????? ，（1）确定矩阵P 的未知元素。 （2）求 P 模最大特征值。 （3）分析矩阵P 的一致性是否可以接受（随机一致性指标ＲＩ取0.58）。 2. 在层次分析法建模中，我们介绍了成对比较矩阵概念，已知矩阵P 是三阶成对比较矩阵 322P ? ???=?????? ，（1）将矩阵P 元素补全。 （2）求P 模最 大特征值。 （3）分析矩阵P 的一致性是否可以接受。 3.考虑下表数据
(1)用曲改直的思想确定经验公式形式。 （2）用最小二乘法确定经验公式系数。 4.. 考虑微分方程 (0.2)0.0001(0.4)0.00001dx x xy dt dy y xy dt εε?=--????=-++?? （1）在像平面上解此微分方程组。（2）计算0ε=时的周期平均值。（3）计算0.1ε=时，y 的周期平均值占总量的周期平均值的比例增加了多少？ 5考虑种群增长模型 '()(1/1000),(0)200x t kx x x =-= （1）求种群量增长最快的时刻。（2）根据下表数据估计参数k 值。 6. 布均匀，若环保部门及时发现并从某时刻起切断污染源，并更新湖水（此处更新指用新鲜水替换污染水），设湖水更新速率是 3 (m r s 单位：）。 （1） 试建立湖中污染物浓度随时间下降的数学模型？ 求出污染物浓度降为控制前的5%所需要的时间。 7. 假如保险公司请你帮他们设计一个险种:35岁起保,每月交费400元,60岁开始领取养老金,每月养老金标准为３６００元，请估算该保险费月利率为多少（保留到小数点后５位）？ 8. 某校共有学生40000人，平时均在学生食堂就餐。该校共有,,A B C 3 个学生食堂。经过近一年的统计观测发现：A 食堂分别有10%，25%的学生经常去B ，C 食堂就餐，B 食堂经常分别有15%，25%的同学去
数学建模课程设计
攀枝花学院 学生课程设计（论文） 题目：产品广告费用分配对销量及利润的影响模型学生姓名：梁忠 学号： 201210802007 所在院(系)：数学与计算机学院 专业：信息与计算科学 班级： 12信本1班 指导教师：马亮亮职称：讲师 2014年12 月19 日 攀枝花学院教务处制
攀枝花学院本科学生课程设计任务书 题目具有自身阻滞作用的食饵—捕食者模型 1、课程设计的目的 数学建模课程设计是让学生通过动手动脑解决实际问题，让学生学完《数学建模》课程后进行的一次全面的综合训练，是一个非常重要的教学环节。 2、课程设计的内容和要求（包括原始数据、技术要求、工作要求等） 根据指导教师所下达的课程设计题目和课程设计要求，在规定的时间内完成设计任务；撰写详细的课程设计论文一份。 3、主要参考文献 【1】姜启源，数学模型（第二版），高等教育出版社，北京。 【2】寿纪麟，数学建模——方法与范例，西安交大出版社。 【3】（美）JOHN A.QUELCH 等著吕—林等译，市场营销管理教程和案例, 北京大学出版社 2000。 【4】戴永良广告绩效评估，中国戏剧出版社，2001。 4、课程设计工作进度计划 序号时间（天）内容安排备注 1 2 分析设计准备周一至周二 2 4 编程调试阶段周三至周一 3 2 编写课程设计报告周二至周三 4 2 考核周四至周五 总计10（天） 指导教师（签字）日期年月日 教研室意见： 年月日 学生（签字）： 接受任务时间：2014 年12 月15 日
注：任务书由指导教师填写。 课程设计（论文）指导教师成绩评定表题目名称具有自身阻滞作用的食饵—捕食者模型 评分项目分 值 得 分 评价内涵 选题15% 01 能结合所学课程知识，有 一定的能力训练。符合选 题要求 5 遵守各项纪律，工作刻苦努力，具有良好的科学 工作态度。 02 工作量适中，难易度合理10 通过实验、试验、查阅文献、深入生产实践等渠 道获取与课程设计有关的材料。 能力水平35% 04 综合运用知识的能力10 能运用所学知识和技能去发现与解决实际问题， 能正确处理实验数据，能对课题进行理论分析， 得出有价值的结论。 05 应用文献的能力 5 能独立查阅相关文献和从事其他调研；能提出并 较好地论述课题的实施方案；有收集、加工各种 信息及获取新知识的能力。 06 设计（实验）能力，方案 的设计能力 5 能正确设计实验方案，独立进行装置安装、调试、 操作等实验工作，数据正确、可靠；研究思路清 晰、完整。 07 计算及计算机应用能力 5 具有较强的数据运算与处理能力；能运用计算机 进行资料搜集、加工、处理和辅助设计等。 08 对计算或实验结果的分析 能力（综合分析能力、技 术经济分析能力） 10 具有较强的数据收集、分析、处理、综合的能力。 成果质量45% 09 插图（或图纸）质量、篇 幅、设计（论文）规范化 程度 5 符合本专业相关规范或规定要求；规范化符合本 文件第五条要求。 10 设计说明书（论文）质量30 综述简练完整，有见解；立论正确，论述充分， 结论严谨合理；实验正确，分析处理科学。 11 创新10 对前人工作有改进或突破，或有独特见解。 成绩 指 导 教 师 评 语 指导教师签名：年月日
数学建模与数学实验试卷及答案
数学建模与数学实验试卷及答案 二、本题10分(写出程序和结果) 蚌埠学院2010—2011学年第二学期 2,x在 [-5 ，5] 区间内的最小值，并作图加以验证。求函数yxe,，,3《数学建模与数学实验》补考试卷答案 f1=inline('x.^2 +exp(-x)-3') 注意事项:1、适用班级:09数学与应用数学本科1,2班 2、本试卷共1页，附答题纸1页。满分100分。 x=fmin(f1,-5,5) 3、考查时间100分钟。 y=f1(x) 4、考查方式:开卷 fplot(f1,[-5,5]) 一、填空:(每空4分，共60分) x = 0.3517，y== -2.1728 123111,，，，， ,,,,三、本题15分(写出程序和结果) 1. 已知，，则A的秩为 3 ，A的特征值为 A,612B,234,,,, ,,,,,215531,，，，，360000xx，,,12,max2.5fxx,，求解:， stxx..250000，,,1212-1.9766 4.4883 + 0.7734i 4.4883 - 0.7734i ，若令 A([1,3],:)= B([2,3],:)，则,x,150001,A(2,:)= 6 1 2 ; 解: xxx，，,22,123,model: 2. 的解为 1.25 ,0.25 0.5 ; xxx，，,521,123max=2.5*x1+x2; ,242xxx，，,123,3*x1+x2<=60000; 装订线内不要答题 2*x1+x2<=50000; 3. 将1234521 分解成质因数乘积的命令为_factor(sym(‘1234521’)),
数学模型课程设计
数学模型课程设计
文档仅供参考，不当之处，请联系改正。 攀枝花学院 学生课程设计（论文） 题目：蔬菜的运输问题 学生姓名：孟蕾 学号： 1080 所在院(系)：数学与计算机学院 专业：信息与计算科学 班级：级信本 指导教师：李思霖 6 月 29 日 攀枝花学院教务处制
攀枝花学院本科学生课程设计任务书
课程设计（论文）指导教师成绩评定表
摘要 本文针对蔬菜的运输问题进行分析，针对蔬菜运输时所需要注意的蔬菜供应量，需求量，运输距离，运输补贴，短缺补偿等约束性条件，运用lingo编程的方法解决如何进行蔬菜运输来分别使各类要求的支出最少的问题。 问题一中，要求如果不考虑短缺补偿，只考虑运费补贴最少，请为该市设计最优蔬菜运输方案。我们将供货商和销售点需求分别编号a和b，数量是从1~8和1~35。从题中能够看出其约束条件，所有销售点从第 A基地获得的蔬菜数量应该等于该基地所 i 生产的蔬菜数量；所有基地给 B销售点提供的蔬菜数量要大于等 j 于0，而且应该小于或等于该点的需求量。 问题二中，增添了对短缺补缺的考虑，规定各蔬菜销售点的短缺量一律不超过需求量的30%，在同时考虑短缺补偿和运费补贴的情况下再次设计最有蔬菜方案。由题意即是要求总费用，具体步骤仍同问题一，需要变化的分别是总费用w的表示式和关于销售点需求的约束条件。w变为原运输补贴的公式再加上每个销售点每吨短缺蔬菜的数量乘上各个销售点不同的短缺补偿，短缺数量需要用各个销售点的需求减去所有基地供给给这个的销售点的蔬菜数量之和。 问题三中，要求增加任意两个基地的生产数量，使得不存在短缺情况出现，然后视运费补贴最小的情况来确定哪两个基地分
数学建模
长江学院课程设计报告 课程设计题目：海岛服务中心的建设问题 姓名1：学号： 姓名2：学号： 姓名3：学号： 专业：材料成型 班级：083115 指导教师：黄雯 2010年11 月01日
海岛服务中心建设 摘要 本论文主要讨论了如何选择海岛服务中心，并使得其工作效率高，经济效益也高，成本低，利润大。选址问题是一种极其重要的长期决策，它的好坏直接影响到服务方式，服务质量，服务效率，服务成本，及才生利润。因此能影响到利润和市场竞争里，决定了企业的命运，甚至影响到本地的经济发展，所以选址问题的研究有着企业和经济发展的重要意义。 “在海岛上建一个服务中心为居民提供各种服务”数学模型是通过服务中心的建立来探讨建在那里比较合适，使得人数多的居民点希望距离近且到各居民点的距离最小。这是海岛服务中心选择地址问题，使得服务中心起的作用效率最大化，即到每个居民点的总时间最短，或者说到每个居民点的距离总和最短，从而经济效益高。在考虑居民点与服务中心之间为直线道路连通的情况下：由于海岛上的居民点比较分散和各居民点的人数也不一样的影响，利用数学知识联系实际问题，作出相应的解答和处理。并运用lingo软件编程和处理相关数据，从而得到最优决策方案。 该问题是一个非线性规划问题，我们首先建立单位目标的优化模型，也即模型一。根据题意得到了模型一的目标函数通过lingo软件的计算，从而使得总距离最短。 经过本小组成员之间的思考和讨论，得出了另一个优化模型，即模型二。根据题意得到了模型二的目标函数通过lingo软件的计算，从而使得总时间最短，效益也为最高。 关键词：服务中心居民点最佳路径方案效率高选地址
数学建模与数学实验课后习题答案
P59 4.学校共1002名学生，237人住在A 宿舍，333人住在B 宿舍，432人住在C 宿舍。学生要组织一个10人的委员会，使用Q 值法分配各宿舍的委员数。 解:设P 表示人数，N 表示要分配的总席位数。i 表示各个宿舍（分别取A,B,C ），i p 表示i 宿舍现有住宿人数，i n 表示i 宿舍分配到的委员席位。 首先，我们先按比例分配委员席位。 A 宿舍为:A n = 365.21002 10237=? B 宿舍为:B n =323.31002 10333=? C 宿舍为:C n =311.4100210432=? 现已分完9人，剩1人用Q 值法分配。 5.93613 22372 =?=A Q 7.92404 33332 =?=B Q 2.93315 44322 =?=C Q 经比较可得，最后一席位应分给A 宿舍。 所以，总的席位分配应为：A 宿舍3个席位，B 宿舍3个席位，C 宿舍4个席位。
商人们怎样安全过河
由上题可求：4个商人，4个随从安全过河的方案。 解：用最多乘两人的船，无法安全过河。所以需要改乘最多三人乘坐的船。 如图所示，图中实线表示为从开始的岸边到河对岸，虚线表示从河对岸回来。商人只需要按照图中的步骤走，即可安全渡河。总共需要9步。
P60 液体在水平等直径的管内流动，设两点的压强差ΔP 与下列变量有关：管径d,ρ,v,l,μ,管壁粗糙度Δ，试求ΔP 的表达式 解：物理量之间的关系写为为()?=?,,,,,μρ?l v d p 。 各个物理量的量纲分别为 []32-=?MT L p ，[]L d =，[]M L 3-=ρ，[]1-=LT v ，[]L l =，[]11--=MT L μ，Δ是一个无量纲量。 ???? ??????-----=?0310100011110010021113173A 其中0=Ay 解得 ()T y 00012111---=， ()T y 00101102--=， ()T y 01003103--=， ()T y 10000004= 所以 l v d 2111---=ρπ，μρπ112--=v ，p v ?=--313ρπ，?=4π 因为()0,,,,,,=??p l v d f μρ与()0,,,4321=ππππF 是等价的，所以ΔP 的表达式为： ()213,ππψρv p =?
数学建模课程设计——优化问题
在手机普遍流行的今天，建设基站的问题分析对于运营商来说很有必要。本文针对现有的条件和题目的要求进行讨论。在建设此模型中，核心运用到了0-1整数规划模型，且运用lingo 软件求解。 对于问题一： 我们引入0-1变量，建立目标函数：覆盖人口最大数=所有被覆盖的社区人口之和，即max=15 1j j j p y =∑，根据题目要求建立约束条件，并用数学软件LINGO 对其模型求解，得到最优解。 对于问题二： 同样运用0-1整数规划模型，建立目标函数时，此处假设每个用户的正常资费相同，所以68%可以用减少人口来求最优值，故问题二的目标函数为：max=∑=15 1j j j k p 上述模型得到最优解结果如下： 关键字：基站; 0-1整数规划；lingo 软件
1 问题的重述．．．．．．．．．．．．．．．．．．．．．．．．．3 2 问题的分析．．．．．．．．．．．．．．．．．．．．．．．．．4 3 模型的假设与符号的说明．．．．．．．．．．．．．．．．．．．5 3.1模型的假设．．．．．．．．．．．．．．．．．．．．．． 5 3.2符号的说明．．．．．．．．．．．．．．．．．．．．．． 5 4 模型的建立及求解．．．．．．．．．．．．．．．．．．．．．. 5 4.1模型的建立．．．．．．．．．．．．．．．．．．．．．． 5 4.2 模型的求解．．．．．．．．．．．．．．．．．．．．．． 6 5 模型结果的分析．．．．．．．．．．．．．．．．．．．．．．．7 6 优化方向．．．．．．．．．．．．．．．．．．．．．．．．．．7 7 参考文献．．．．．．．．．．．．．．．．．．．．．．．．．．8 8、附录．．．．．．．．．．．．．．．．．．．．．．．．．．． 9