《数学模型》实验指导书
实验项目与学时分配表
实验项目一：预测模型
一、实验目的和要求：
熟悉科学计算软件MATLAB的图形功能，会用软件画图，并进行数据模拟。掌握数据预测方法。
二、实验内容：
某乡镇企业2010-2016年的生产利润如下表；
试预测2017年和2018年的利润。
三、过程：
１．利用MATLAB软件或其它绘图软件，对所给数据画出散点图；
２．根据散点图，分析合适的函数，并试探（画图作对比）；
３．确定函数类型，作数据拟合，确定函数中的参数；
４．作误差分析；
5. 考虑用cftool进行数据拟合，并作出结果分析.
实验项目二：初等模型
一、实验目的和要求：
掌握线性方程组建模，并会用它解决一些实际问题；熟悉科学计算软件MATLAB求线性方程组的命令。
二、实验内容：
问题一：某城市有下图所示的交通图, 每条道路都是单行线, 需要调查每条道路每小时的车流量. 图中的数字表示该条路段的车流数. 如果每个交叉路口进入和离开的车数相等,
图4 某城市单行线车流量
(1)建立确定每条道路流量的线性方程组；
(2)使用MATLAB求线性方程组；
(3)分析哪些流量数据是多余的；
(4)为了唯一确定未知流量, 需要增添哪几条道路的流量统计；
问题二：某地有一座煤矿, 一个发电厂和一条铁路. 经成本核算, 每生产价值1元钱的煤需消耗0.3元的电; 为了把这1元钱的煤运出去需花费0.2元的运费; 每生产1元的电需0.6元的煤作燃料; 为了运行电厂的辅助设备需消耗本身0.1元的电, 还需要花费0.1元的运费; 作为铁路局, 每提供1元运费的运输需消耗0.5元的煤, 辅助设备要消耗0.1元的电. 现煤矿接到外地6万元煤的订货, 电厂有10万元电的外地需求, 问: 煤矿和电厂各生产多少才能满足需求
过程：
（1）建立确定的线性方程组；
（2）使用MATLAB求线性方程组；
实验项目三：简单优化模型
一、实验目的和要求：
学习函数极值的相关知识，熟悉科学计算软件MATLAB求极值的方法。
二、实验内容：
问题一：一栋楼房的后面是一个很大的花园。在花园中紧靠着楼房有一个温室，温室伸
入花园宽2m ，高3m ，温室正上方是楼房的窗台。清洁工打扫窗台周围，他得用梯子越过温室，一头放在花园中，一头靠在楼房的墙上。因为温室是不能承受梯子压力的，所以梯子太短是不行的。现清洁工只有一架7m 长的梯子，你认为它能达到要求吗？能满足要求的梯子的最小长度为多少？
过程：
1、 设温室宽为a ，高为b ，梯子倾斜的角度为x ，当梯子与温室顶端 A 处恰 好接触时，梯子的长度L 只与x 有关。试写出函数L （x ） 及其定义域。
3、 在 Matlab 环境，先用命令 clear x 清除x 的值，再定义函数L(x) ，并求导。
4、 将a 、b 赋值，画出L(x) 的图形。注意自变量x 的范围选取。
5、 求驻点，即求方程()0L x '=的根，有什么命令求根？并计算函数在驻点的值。驻点唯一吗？
6、 观测图形，选取初始点，用fminbnd 直接求L(x)的极小值。并与（5）的结果比较。
7、 取a=2,b=2.8,重新运行程序，结果如何？
问题二：在某医院走廊拐角处，垂直相交的两通道宽度分别是1m 与1.5m, 病床宽为0.80m ，问病床至多为多长才能被推过此拐角？ 过程：
1 建立数学模型；
2 求解数学模型；
3 改动模型中一些数据，再求解，观测结果。
实验项目四：数学规划问题
一、实验目的和要求：
熟悉数学规划软件LINGO 的运用，理解数学规划模型及其应用，建模过程和求解方法。重点是模型的约束条件的建立和结果的分析。
二、实验内容：
1、某银行经理计划用一笔资金进行有价证券的投资，可供购进的证券以及信用等级、到期年限、收益如表2所示。按照规定，市政证券的收益可以免税，其他证券的收益需按50%的税率纳税。此外还有以下限制：
（1）政府及代办机构的证券总共至少要够进400万元；
（2）所购证券的平均信用等级不超过1.4（信用等级数字越小，信用程度越高）； （3）所购证券的平均年限不超过5年；
表2证券以及信用等级、到期年限、收益
2)如果能够以2.75%的利率借到不超过100万元资金，该经理应如何操作？
3)在1000万元资金情况下，若证券A的税前收益增加为4.5%，投资应否改变？若证券C 的税前收益减少为4.8%，投资应否改变？
实验项目五：微分方程模型
一、实验目的和要求：
理解一阶、二阶微分法在建模过程中的应用，熟悉利用MATLAB软件求解微分方程的方法。注意模型的普遍性和模型的广泛性。
二、实验内容：
问题一：一个半球体状的雪堆，其体积V的融化速率与半球面面积S成正比，比例系数K>0. 假设在融化过程中雪堆始终保持半球体状，已知初始半径为r0的雪堆在开始融化的3小时内，融化了其原体积的7/8，问该雪堆全部融化需要多少时间？
过程：
１．分析雪堆的融化过程；
２．建立雪堆融化的微分方程模型；
３．利用所给数据，确定参数；
４．确定初始条件，求解方程（模型）．
５．扩展讨论：雪堆形状不同时的建模和求解方法（供参考，不作要求）
问题二：现有一只兔子、一匹狼，兔子位于狼的正西100米处，假设兔子与狼同时发现对方并一起起跑，兔子往正北60米处的巢穴跑，而狼在追兔子。已知兔子、狼是匀速跑且狼的速度是兔子的两倍。问兔子能否安全回到巢穴？
过程：（1）建立狼的运动轨迹微分模型。
（2）画出兔子与狼的运动轨迹图形。
（3）用解析方法求解，问兔子能否安全回到巢穴？
（4）用数值方法求解，问兔子能否安全回到巢穴？
实验项目六：差分方程模型
一、实验目的和要求：
理解差分法在建模过程中的应用，熟悉利用MATLAB软件求解差分方程的方法。注意模型的普遍性和模型的广泛性。
二、实验内容：
1、作为房产公司的代理人，你要迅速准确地回答客房各方面的问题。现有人看中了贵公司一套建筑面积为S（m2），单价为P（元/ m2）的房子。他计划首付30%，其余70%用20年按揭贷款（贷款年利率r）。请你提供下列信息：房屋总价格、首付款额、月付还款额。当S=120m2，P=5200元/ m2，r=5.58% 时上述三个值。
要求：先求出房屋总价格、首付款额、月付还款额三者的符号解；再求出当S=120m2，P=5200元/ m2，r=5.58% 时三者的数值解。
过程：（1）给出模型假设及建立相应的差分方程；
（2）利用递推公式法求解差分方程的符号解；
（3）利用Matlab求解差分方程的符号解；
（4）求出当S=120m2，P=5200元/ m2，r=5.58% 时三者的数值解；
（5）分析按月还款与按年还款哪种对贷款者更有利？
2、一年生植物春季发芽，夏天开花，秋季产种，没有腐烂，风干，被人为掠取的那些种子可以活过冬天，其中一部分能在第2年春季发芽，然后开花，产种，其中的另一部分虽未能发芽，但如又能活过一个冬天，则其中一部分可在第三年春季发芽，然后开花，产种，如此继续，一年生植物只能活1年，而近似的认为，种子最多可以活过两个冬天，试建立数学模型研究这种植物数量变化的规律，及它能一直繁殖下去的条件。
过程：（1）记一棵植物秋季产种的平均数为c,种子能活过一个冬天的(1岁种子)比例为b,活过一个冬天没有发芽又活过一个冬天的（2岁种子）比例仍为b,1岁种子发芽率a1，2岁种子发芽率a2。建立相应的差分方程研究这种植物数量变化的规律；
（2）利用递推公式法求解差分方程的符号解；
（3）利用Matlab求解差分方程的符号解；
（4）考察当X0=100,a1=0.5,a2=0.25,c=10,b=0.18~0.20这种植物能一直繁殖的条件
实验项目七：层次分析法模型
一、实验目的和要求：
理解层次分析法，熟悉层次分析法在实际问题中的应用。注意层次分析法建模的几个基本步骤，重点是成对比较矩阵的建立和一致性检验。
二、实验内容：
假设你是一位应届毕业生，现有P1、P2、P3等三个就业单位可供你选择。结合你的实际情况，建立一个优选模型，作出你的最优选择。
P1：广东某计算机软件公司，从事软件编程工作；
P2：上海某国际贸易公司，从事报关工作；
P3：武汉某机械制造公司，从事生产管理工作。
过程：
１．分析自己选择就业岗位所要考虑的因素；
２．构建目标层、准则层、方案层的层次结构模型；
３．利用１－９尺度，依据个人的认同，构造各成对比较矩阵；
４．对各成对比较矩阵进行一致性检验，不通过，应作修改，直到全部通过；
５．求权向量和组合权向量，并作组合权向量的一致性检验，不通过，应作修改（主要是对一致性比率较大的成对比较矩阵作调整），直到全部通过．
６．根据组合权向量，选取权向量最大的方案，作为决策方案．
实验项目八：概率统计模型
一、实验目的和要求：
理解概率统计模型的应用，了解蒙特卡洛算法来进行模拟随机数。
二、实验内容：
某报童以每份0.03元的价格买进报纸,以0.05元的价格出售. 根据长期统计,报纸每天的销售量及百分率为
已知当天销售不出去的报纸，将以每份0.02元的价格退还报社.试用模拟方法确定报童每天买进报纸数量,使报童的平均总收入为最大?
过程：（1）给出模型假设；
（2）让时间从1开始变化到365产生销售量的随机数；
（3）计算当天的收入，计算累计利润；
（4）通过比较得出最优订货量。
数学建模实验答案-概率模型
实验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*。 其中：
数学建模实验报告
一、实验目的 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时要转换
《数学建模与数学实验》本科教学日历 数学建模部分 开设课程课程名称数学建模课程编号0701107 施教单位理学院 课内学时 总课时36 课程性质公共基础讲授课时28 修读要求选修实践课时8 选用教材教材名称数学建模教程出版社名称高等教育出版社 出版时间 及版次 2011年出版，第一版印刷时间2011年 其他情况 教学安排 班次授课对象及人数任教教员（指导教员）姓名及职称数学建模A 各专业本科学员 吴孟达教授 段晓君教授 毛紫阳讲师 王丹讲师 数学建模B 各专业本科学员 吴孟达教授 段晓君教授 毛紫阳讲师 王丹讲师 课次节 次 授课内容 教学 方法 采用现代化教学手段(课时) 多媒体电教双语网络实验 1 1 （1）什么是数学建模？数学建模的一般概念 （2）几个数学建模问题 讲授 1 2 （1）数学建模的一般步骤 （2）敏感问题调查案例 讲授 1 2 3 （1）行走步长问题 （2）雨中行走淋雨量最小问题 （3）道路是越多越通畅吗? 讲授 1 4 （1）有奖销售的抽奖策略问题 （2）“非诚勿扰”女生最佳选择问题 （3）网络文章流行度预测和招聘匹配 讲授 1 3 5 （1）线性规划模型基本概念 （2）整数规划模型 （3）0-1规划模型 讲授 1 6 （1）非线性规划 （2）多目标规划 讲授 1 4 7 （1）最短路算法 （2）最小生成树算法 讲授 1 8 （1）最大流算法 （2）PageRank算法 讲授 1 5 9 规划模型上机实践实践 1
课次节 次 授课内容 教学 方法 采用现代化教学手段(课时) 多媒体电教双语网络实验10 图论模型上机实践实践 1 6 11 （1）博弈模型基本概念 （2）Nash平衡和Pareto最优 （3）博弈论案例 讲授 1 12 （1）贝叶斯纳什均衡 （2）拍卖模型 讲授 1 7 13 社会选择理论中的选举问题数学模型-阿罗不可能定理讲授 1 14 越野长袍团体赛排名规则公平性问题讲授 1 8 15 军事作战模型-Lanchester作战模型讲授 1 16 自动化车床管理模型讲授 1 9 17 （1）“边际效应”基本概念 （2）实物交换模型，最佳消费模型、报童售报问题 讲授 1 18 （1）价格弹性模型 （2）合作效益的Shapley值分配模型 讲授 1 10 19 （1）聚类分析基本概念 （2）常用聚类算法 讲授 1 20 （1）方差分析基本概念 （2）单因素方差分析 （3）双因素方差分析 讲授 1 11 21 （1）主成分分析基本概念 （2）因子分析 讲授 1 22 （1）一元回归分析 （2）多元回归分析 （3）多元回归模型的检验与优化 讲授 1 12 23 聚类分析和方差分析上机实践实践 1 24 主成分分析和多元回归分析上机实践实践 1 13 25 （1）遗传算法基本思想 （2）算法步骤 讲授 1 26 遗传算法计算实例讲授 1 14 27 （1）模拟退火算法基本思想 （2）算法步骤 讲授 1 28 模拟退火算法计算实例讲授 1 15 29 （1）蚁群算法基本思想 （2）算法步骤 讲授 1 30 （1）数学建模中的计算机仿真 （2）不可召回的秘书招聘问题 （3）车灯光源优化设计 （4）生命游戏 讲授 1 16 31 遗传算法上机实践实践 1 32 模拟退火算法上机实践实践 1
实验二 SolidWorks 草绘特征和放置特征操作（一） 一、 实验目的 1. 掌握基本零件建模的一般步骤和方法 2. 掌握SolidWorks 草绘特征：拉伸凸台、拉伸切除、旋转凸台、旋转切除、扫描、 放样的操作方法。 3. 掌握放置（应用）特征：钻孔特征、倒角特征、圆角特征、抽壳特征、拔模斜度特 征、筋的操作方法 二、 实验内容 完成下列下列零件造型 三、 实验步骤 1. 连接件设计 完成如图 1 (1) (2) 2 所示。 图 1连接件 图 2草图 (3) 单击【拉伸凸台/ 框内选择【两侧对称】选项，在【深度】文本框内输入“54mm ”，单击【确定】按钮，如图 3所示。 图 3 “拉伸”特征 (4) 120°”，然后 在第二参考中选择图形的一条下边线。单击【确定】按钮，建立新基准面，如
错误！未找到引用源。所示。 (5) 1，选择“反转法线” 1，单选择 4所示。 图4草图 图4建立基准面 底面边线
(6) 单击【拉伸凸台/ 列表框内选择【给定深度】选项，在【深度】文本框内输入“12mm”，单击【确定】按钮，如图5所示。 图5“拉伸”特征 (7)选取基体上表面，单击【草图绘制】进入草图绘制，使用中心线工具在 上表面的中心位置绘制直线，注意不要捕捉到表面边线，如图6所示。 图 6 中心线 (8) 内输入“8mm”，在图形区域选择中心线，在属性管理器中选中【添加尺寸】、【选择链】、【双向】和【顶端加盖】复选框，选中【圆弧】单选按钮，单击【确定】按钮，标注尺寸，完成草图，如图7所示。 运用“等距实体”绘制草图 (8) -拉伸】属性管理器，在【终止条件】下拉 列表框内选择【完全贯穿】选项，单击【确定】按钮，如图8所示。
数学建模实验报告 实验一计算课本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查看其用法）。 ★编写的程序和运行结果： 程序：
数学建模与数学实验课程总结与练习内容总结 第一章 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%的同学去
数学建模习题指导 第一章 初等模型 讨论与思考 讨论题1 大小包装问题 在超市购物时你注意到大包装商品比小包装商品便宜这种现象吗？比如洁银牙膏50g 装的每支1.50元，120g 装的每支3.00元，二者单位重量的价格比是1.2:1，试用比例方法构造模型解释这种现象。 （1）分析商品价格C 与商品重量w 的关系。 （2）给出单位重量价格c 与w 的关系，并解释其实际意义。 提示： 决定商品价格的主要因素：生产成本、包装成本、其他成本。 单价随重量增加而减少 单价的减少随重量增加逐渐降低 思考题2 划艇比赛的成绩 赛艇是一种靠浆手划桨前进的小船，分单人艇、双人艇、四人艇、八人艇四种。各种艇虽大小不同，但形状相似。T.A.McMahon 比较了各种赛艇1964—1970年四次2000m 比赛的最好成绩(包括1964年和1968年两次奥运会和两次世界锦标赛)，见下表。建立数学模型解释比赛成绩与浆手数量之间的关系。 各种艇的比赛成绩与规格 γβα++=3 2w w C w w c γβα++=-3 123 431w w c γβ--='-3 2943 4w w c γβ+=''-
第二章 线性代数模型 森林管理问题 森林中的树木每年都要有一批砍伐出售。为了使这片森林不被耗尽且每年都有所收获，每当砍伐一棵树时，应该就地补种一棵幼苗，使森林树木的总数保持不变。被出售的树木，其价值取决于树木的高度。开始时森林中的树木有着不同的高度。我们希望能找到一个方案，在维持收获的前提下，如何砍伐树木，才能使被砍伐的树木获得最大的经济价值。 思考： 试解释为什么模型中求解得到的 为每周平均销售量会略小于模型假设中给出的1。 练习： 将钢琴销售的存贮策略修改为：当周末库存量为0或1时订购，使下周初的库存 达到3架；否则，不订购。建立马氏链模型，计算稳态下失去销售机会的概率和每周的平均销售量。 2.将钢琴销售的存贮策略修改为：当周末库存量为0时订购本周销售量加2架；否则，不订购。建立马氏链模型，计算稳态下失去销售机会的概率和每周的平均销售量。 第三章 优化模型 讨论题 1）最优下料问题 用已知尺寸的矩形板材加工半径一定的圆盘。给出几种加工排列方法，比较出最优下料方案。 2）广告促销竞争问题 甲乙两公司通过广告竞争销售商品，广告费分别为 x 和 y 。设甲乙公司商品的售量在两公司总售量中所占份额是它们的广告费在总广告费中所占份额的函数 又设公司的收入与售量成正比，从收入中扣除广告费后即为公司的利润。试构造模型的图形，并讨论甲公司怎样确定广告费才能使利润最大。 （1）令 （2）写出甲公司的利润表达式 对一定的 y ，使 p (x ) 最大的 x 的最优值应满足什么关系。用图解法确定这个最优值。 练习1 三个家具商店购买办公桌：A 需要30张，B 需要50张，C 需要45张。这些办公桌由两个工厂供应：工厂1生产70张，工厂2生产80张。下表给出了工厂和商店的距离（单位公里） ， 857.0=n R ) (),(y x y f y x x f ++的示意图。。画出则)()()(,t f t f t f y x x t 11=-++= 。 )(t p
数学建模与数学实验报告 指导教师__郑克龙___ 成绩____________ 组员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)
数学建模与数学实验试卷及答案 二、本题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’)),
质量-弹簧-阻尼系统实验教学指导书 北京理工大学机械与车辆学院 2016.3
实验一：单自由度系统数学建模及仿真 1 实验目的 （1）熟悉单自由度质量-弹簧-阻尼系统并进行数学建模； （2）了解MATLAB 软件编程，学习编写系统的仿真代码； （3）进行单自由度系统的仿真动态响应分析。 2 实验原理 单自由度质量-弹簧-阻尼系统，如上图所示。由一个质量为m 的滑块、一个 刚度系数为k 的弹簧和一个阻尼系数为c 的阻尼器组成。系统输入：作用在滑块上的力f (t )。系统输出：滑块的位移x (t )。 建立力学平衡方程： m x c x kx f ??? ++= 变化为二阶系统标准形式： 22f x x x m ζωω?? ? ++= 其中：ω是固有频率，ζ是阻尼比。 ω= 2c m ζω= = 2.1 欠阻尼(ζ<1)情况下，输入f (t )和非零初始状态的响应： ()()sin()))] t t x t t d e ζωττζωττ +∞ --=? -= -+-?
2.2 欠阻尼(ζ<1)情况下，输入f(t)=f0*cos(ω0*t) 和非零初始状态的的响应： 022 3 00 22222 00 222222 2 ()cos(arctan()) 2f [(0)]cos() [()(2)] sin( t t x t t x e k e ζω ζω ζωω ω ωω ζωω ωωζωω - ? - =- - ++ -+ +) 输出振幅和输入振幅的比值：A= 3 动力学仿真 根据数学模型，使用龙格库塔方法ODE45求解，任意输入下响应结果。 仿真代码见附件 4 实验 4.1 固有频率和阻尼实验 （1）将实验台设置为单自由度质量-弹簧-阻尼系统。 （2）关闭电控箱开关。点击setup菜单，选择Control Algorithm，设置选择Continuous Time Control，Ts=0.0042，然后OK。 （3）点击Command菜单，选择Trajectory，选取step,进入set-up,选取Open Loop Step 设置(0)counts, dwell time=3000ms,(1)rep, 然后OK。此步是为了使控制器得到一段时间的数据，并不会驱动电机运动。 （4）点击Data菜单，选择Data Acquisition,设置选取Encoder#1 ，然后OK离开；从Utility菜单中选择Zero Position使编码器归零。 （5）从Command菜单中选择Execute，用手将质量块1移动到2.5cm左右的位置（注意不要使质量块碰触移动限位开关），点击Run, 大约1秒后，放开手使其自由震荡，在数据上传后点击OK。 （6）点击Plotting菜单，选择Setup Plot，选取Encoder #1 Position；然后点击Plotting 菜单，选择Plot Data，则将显示质量块1的自由振动响应曲线。 （7）在得到的自由振动响应曲线图上，选择n个连续的振幅明显的振动周期，计算出这段振动的时间t，由n/t即可得到系统的频率，将Hz转化为rad/sec即为系统的振动频率ω。
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:海底曲面图：
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 =?
数学建模实验指导书 数学建模实验项目一 初等模型 一、 实验目的与意义： 1、练习初等问题的建模过程；熟悉数学建模步骤 2、练习Matlab 基本编程命令； 二、 实验要求： 1、较能熟练应用Matlab 基本命令和函数； 2、注重问题分析与模型建立，了解建模小论文的写作过程； 3、提高Matlab 的编程应用技能。 三、 实验学时数：4学时 四、 实验类别：综合性 五、 实验内容与步骤： 练习： 基本命令 ：循环、绘图、方程（组）求解 作业： 1、某大学青年教师从31岁开始建立自己的养老基金，他把已有的积蓄10000元也一次性地存入，已知月利率为0.001（以复利计），每月存入700元，试问当他60岁退休时，他的退休基金有多少？又若他退休后每月要从银行提取1000元，试问多少年后他的基金将用完？ 2、试对公平席位分配问题进行编程求解。 3、编程求解差分方程的阻滞增长模型1(1)k k k x bx x +=-，分别令b 从1.8逐渐增加，考察序列k x 收敛、2倍周期收敛、4倍周期收敛……，直至一片混乱的情况，试以b 为横坐标，收敛点为纵坐标作图。（与7.3节图8比较）。 数学建模实验项目二 数学规划 一、实验目的与意义： 1、认识数学规划的建模过程； 2、认识数学规划的各种形式和解法。 二、实验要求： 1、熟练应用Matlab 、lindo 、lingo 求解工具箱求解数学规划； 2、掌握建立数学规划的方法和步骤； 3、提高Matlab 、lindo 、lingo 的编程应用技能。 三、实验学时数：4学时 四、实验类别：综合性 五、实验内容与步骤： 练习： 1、奶制品生产销售计划问题的再讨论。 2、自来水输送问题。 3、货机装运问题。 4、选课策略问题。 5、第四章 习题4的模型求解及灵敏度分析。 6、第四章 习题6的模型求解及灵敏度分析。 作业： 1、市场上有n 种资产i s （i=1,2……n ）可以选择，现用数额为M 的相当大的资金作一个时期
10级信息《数学建模与数学实验(实践)》任务书 一、设计目的 通过《数学建模与数学实验(实践)》实践环节，掌握本门课程的众多数学建模方法和原理，并通过编写C语言或matlab程序，掌握各种基本算法在计算机中的具体表达方法，并逐一了解它们的优劣、稳定性以及收敛性。在熟练掌握C 语言或matlab语言编程的基础上，编写算法和稳定性均佳、通用性强、可读性好，输入输出方便的程序，以解决实际中的一些科学计算问题。 二、设计教学内容 1线性规划（掌握线性规划的模型、算法以及Matlab 实现）。整数线性规划（掌握整数线性规划形式和解法）。 2微分方程建模（掌握根据规律建立微分方程模型及解法；微分方程模型的Matlab 实现）。 3最短路问题（掌握最短路问题及算法，了解利用最短路问题解决实际问题）。 行遍性问题（了解行遍性问题，掌握其TSP算法）。 4回归分析（掌握一元线性回归和多元线性回归，掌握回归的Matlab实现）。 5计算机模拟（掌握Monte-carlo方法、了解随机数的产生；能够用Monte-carlo 解决实际问题）。 6插值与拟合（了解数据拟合基本原理，掌握用利用Matlab工具箱解决曲线拟合问题）。 三、设计时间 2012—2013学年第1学期：第16周共计一周 目录 一、10级信息《数学建模与数学实验(实践)》任务书 (1) 二、饭店餐桌的布局问题 (3) 摘要 (3)
问题重述 (3) 模型假设 (3) 模型分析 (4) 模型的建立和求解 (4) 模型推广 (9) 参考文献 (9) 三、白酒配比销售问题 (10) 摘要 (10) 问题重述 (11) 问题分析 (12) 模型假设 (12) 符号及变量说明 (12) 模型的建立与求解 (13) 模型的检验 (18) 模型的评价与推广 (19) 附录 (21) 饭店餐桌的布局问题 摘要 饭店餐桌的布局对于一个饭店有着很重要的作用。本文讨论的就是饭店餐桌的布局问题，根据实际需求及规定建立模型，同时考虑餐桌的类型及规格，尤其是餐桌的摆放技巧，保证使饭店能容纳的人数达到最大。根据所需餐桌的数量
《控制工程基础》实验指导书常熟理工学院机械工程学院 2009.9
目录 1.MATLAB时域分析实验 (2) 2.MATLAB频域分析实验 (4) 3.Matlab校正环节仿真实验 (8) 4.附录：Matlab基础知识 (14)
实验1 MATLAB 时域分析实验 一、实验目的 1． 利用MATLAB 进行时域分析和仿真。 要求：(1)计算连续系统的时域响应（单位脉冲输入，单位阶跃输入，任意输入）。 2．掌握Matlab 系统分析函数impulse 、step 、lsim 、roots 、pzmap 的应用。 二、实验内容 1．已知某高阶系统的传递函数为 ()265432 220501584223309240100 s s G s s s s s s s ++=++++++，试求该系统的单位脉冲响应、单位阶跃响应、单位速度响应和单位加速度响应。 MATLAB 计算程序 num=[2 20 50]; den=[1 15 84 223 309 240 100]; t= (0: 0.1: 20); figure (1); impulse (num,den,t); %Impulse Response figure (2); step(num,den,t);%Step Response figure (3); u1=(t); %Ramp.Input hold on; plot(t,u1); lsim(num,den,u1,t); %Ramp. Response gtext(‘t’); figure (4); u2=(t.*t/2);%Acce.Input u2=(0.5*(t.*t)) hold on; plot(t,u2); lsim(num,den,u2,t);%Acce. Response
数学实验报告 班级：应用物理091 学号：0902052016 姓名：李海瑞
2011.11.26 探索实验二非线性迭代 一，实验指导书解读 迭代是数学研究中的一个非常重要的工具，通过函数或向量函数由初始结点生成迭代结 点列，也可通过函数或向量函数由初值（向量）生成迭代数列或向量列。本实验主要分四步：首先利用蛛网图和迭代数列研究不动点的类型； 其次通过蛛网图和迭代数列研究Logistic映射，探索周期点的性质、认识混沌现象； 第三通过迭代数列或向量列求解方程（组）而寻求有效的求解方法； 最后，利用结点迭代探索分形的性质。 二,试验计划 2.1.迭代序列与不动点 对函数的迭代过程，我们可以用几何图象来直观地显示它——“蜘蛛网”。运行下列 Mathematica程序： Clear[f] f[x_] := (25*x - 85)/(x + 3); （实验时需改变函数） Solve[f[x]==x , x] （求出函数的不动点） g1=Plot[f[x], {x, -10, 20}, PlotStyle -> RGBColor[1, 0, 0], DisplayFunction -> Identity]; g2=Plot[x, {x, -10, 10}, PlotStyle -> RGBColor[0, 1, 0], DisplayFunction -> Identity]; x0=5.5; r = {}; r0=Graphics[{RGBColor[0, 0, 1], Line[{{x0, 0}, {x0, x0}}]}]; For[i = 1, i <= 100, i++, r=Append[r, Graphics[{RGBColor[0, 0, 1], Line[{{x0, x0}, {x0, f[x0]}, {f[x0], f[x0]}}] }]]; x0=f[x0] ]; Show[g1, g2, r, r0, PlotRange -> {-1, 20}, （PlotRange控制图形上下范围） DisplayFunction -> $DisplayFunction] x[0]=x0; x[i_]:=f[x[i-1]]; （定义序列） t=Table[x[i],{i,1,10}]//N ListPlot[t] （散点图） 如果只需迭代n次产生相应的序列，用下列Mathematica程序：
实验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 =0.75, a =1, c =0.6，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=0.75; c=0.6; 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 =2.0m的成品钢材，由粗轧设备等因素决定的粗轧冷却后钢材长度的均方差σ=0.2m，问这时钢材长度的均值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*。 其中：