数学建模第四篇典型案例分析.ppt

第四篇第四篇 典型案例分析典型案例分析§1 投篮的出手角度投篮的出手角度§2 水塔流量估计水塔流量估计§3 钢管订购和运输钢管订购和运输§1.1 §1.1 问题的提出问题的提出L L H0Ph从罚球点投篮示意图从罚球点投篮示意图§1 §1 投篮的出手角度投篮的出手角度问题问题 ☞☞不考虑篮球和篮框大小不考虑篮球和篮框大小,讨论球心命中框心的条件讨论球心命中框心的条件 ☞☞考虑篮球和篮框大小考虑篮球和篮框大小,讨论球心命中框心且入框条件讨论球心命中框心且入框条件 ☞☞保证球入框保证球入框,出手角度和出手速度允许的最大偏差出手角度和出手速度允许的最大偏差 ☞☞考虑空气阻力的影响考虑空气阻力的影响 d0△xD§1.2 §1.2 问题的分析问题的分析0dD篮球入篮球入框框球心偏前球心偏前 ☞☞不考虑篮球和篮框大小不考虑篮球和篮框大小,讨论球心命中框心的条件讨论球心命中框心的条件 ☞☞考虑篮球和篮框大小考虑篮球和篮框大小,讨论球心命中框心且入框条件讨论球心命中框心且入框条件 ☞☞保证球入框保证球入框,出手角度和出手速度允许的最大偏差出手角度和出手速度允许的最大偏差 ☞☞考虑空气阻力的影响考虑空气阻力的影响 §1.3 §1.3 基本模型基本模型★ ★ 不考虑篮球和篮框大小不考虑篮球和篮框大小, ,不考虑空气阻力的影响不考虑空气阻力的影响 则球心运动轨迹为抛物线则球心运动轨迹为抛物线 球可视为质点（球心）的斜抛运动，其运动方程球可视为质点（球心）的斜抛运动，其运动方程 以以 代入上式，就得到球心命中代入上式，就得到球心命中框框心的条件心的条件且有且有最小出手速度最小出手速度★ ★ 球入篮框时的入射角度球入篮框时的入射角度 计算公式计算公式★ ★ 考虑篮球和篮框的大小，考虑篮球和篮框的大小，球心命中球心命中框框心且球入心且球入框框的条件为的条件为 代入得代入得 将将0dD§1.4 §1.4 出手角度和出手速度最大偏差估计出手角度和出手速度最大偏差估计球入球入框框时球心偏前时球心偏前(偏后偏后)的最大距离的最大距离 因为因为★ ★ 出手角度最大偏差估计出手角度最大偏差估计★ ★ 出手速度最大偏差估计出手速度最大偏差估计的相对偏差为的相对偏差为 只考只考虑水平方向的阻力，且阻力与速度成正比水平方向的阻力，且阻力与速度成正比§1.5 §1.5 空气阻力的影响空气阻力的影响不考虑篮球和篮不考虑篮球和篮框框大小，确定球心命中大小，确定球心命中框框心的条件心的条件水平方向运动水平方向运动解得解得h(m)v(m/s)(度度)1.87.678952.60121.97.598552.01812.07.518651.42902.17.439250.8344§1.6 §1.6 算法实现和计算结果算法实现和计算结果★ ★ 对不同出手高度的最小出手速度和相不同出手高度的最小出手速度和相应的出手角度的出手角度v（（m/s））h (m) (度度) (度度) (度度) (度度)8.01.81.92.02.162.4099 42.7925 53.8763 20.921363.1147 40.9188 55.8206 20143163.7281 39.1300 57.4941 19.647864.2670 37.4017 5809615 19.36988.51.81.92.02.167.6975 37.5049 62.1726 12.625068.0288 36.0075 63.1884 12.775368.3367 34.5214 64.1179 13.024068.6244 33.0444 64.9279 13.35839.01.81.92.02.171.0697 34.1327 67.1426 7.655071.2749 3207614 67.7974 8.1663 71.4700 31.3881 68.4098 8.732171.6561 30.0127 68.9840 9.3472★★对不同出手速度和出手高度的出手角度和入射角度不同出手速度和出手高度的出手角度和入射角度h(m) (度度) v(m/s)1.862.4099 8.067.6975 8.571.0697 9.0-0.7562 0.0528-0.5603 0.0694-0.4570 0.08031.2261 0.65970.8276 0.81670.6431 0.89252.063.7281 8.068.3367 8.571.4700 9.0-0.7100 0.0601-0.5411 0.0734-0.4463 0.08321.1140 0.75110.7918 0.86400.6244 0.9243★ ★ 出手角度和出手速度最大偏差出手角度和出手速度最大偏差★ ★ 空气阻力的影响空气阻力的影响v(m/s））h (m)8.01.81.92.02.160.7869 43.5424 61.6100 41.569362.3017 39.715662.9012 37.943362.4099 42.7925 63.1147 40.9188 63.7281 39.1300 64.2670 37.4017 8.51.81.92.02.166.5719 37.790566.9244 36.287067.2505 34.798267.5541 33.320967.6975 37.5049 68.0288 36.0075 68.3367 34.5214 68.6244 33.0444 9.01.81.92.02.170.1198 34.273670.3328 32.908770.5352 31.542870.7279 30.175671.0697 34.1327 71.2749 3207614 71.4700 31.3881 71.6561 30.0127 对不同出手速度和出手高度的出手角度和入射角度对不同出手速度和出手高度的出手角度和入射角度§1.7 §1.7 结果分析结果分析1．最小出手速度和出手角度（表．最小出手速度和出手角度（表4-1-1））2．出手速度和出手高度对出手角度的影响（表．出手速度和出手高度对出手角度的影响（表4-1-2））3．出手角度和出手速度的允许偏差（表．出手角度和出手速度的允许偏差（表4-1-3））4．空气阻力的影响．空气阻力的影响思考思考 若出手高度和角度固定，考察阻力对出手速若出手高度和角度固定，考察阻力对出手速度的影响。

度的影响§2 §2 水塔流量估计水塔流量估计§2.1 §2.1 问题的提出问题的提出时刻（刻（h））水位水位(cm) 0 0.92 1.84 2.95 3.87 4.98 5.90 7.01 7.93 8.97968 948 931 913 898 881 869 852 839 822时刻（刻（h））水位水位(cm)9.98 10.92 10.95 12.03 12.95 13.88 14.98 15.90 16.83 // // 1082 1050 1021 994 965 941 918时刻（刻（h））水位水位(cm)17.93 19.04 19.96 20.84 22.014 22.96 23.88 24.99 25.91 892 866 843 822 // // 1059 1035 1018已知一天水位测量记录已知一天水位测量记录.圆柱形水塔圆柱形水塔 :高高12.2、直径、直径17.4米米 水位降至约水位降至约8.2米升到约米升到约10.8米时米时, ,水泵工作水泵工作. .水泵每天供水一两次，每次约两小时水泵每天供水一两次，每次约两小时. . 估计任何时刻流水量、一天总用水量估计任何时刻流水量、一天总用水量. . §2.2 §2.2 问题的分析问题的分析流量是单位时间流出的水的体积流量是单位时间流出的水的体积.(2)拟合水位拟合水位~时间函数、求导数，得连续时间流量时间函数、求导数，得连续时间流量.★★水泵不工作时段流量计算：水位对时间的变化率水泵不工作时段流量计算：水位对时间的变化率.★★水泵供水时段的流量计算水泵供水时段的流量计算：：(1) 数值微分、拟合数值微分、拟合;★★模型检验：模型检验：t=0-8.97,水位下降水位下降968-822=146,这一时段用水量这一时段用水量=146*水塔的截面积水塔的截面积.§2.3 §2.3 模型假设模型假设1．流量只取决于水位差，与水位本身无关．流量只取决于水位差，与水位本身无关.2．水泵第．水泵第1次供水：次供水：t=9-11; 第第2次供水次供水:t=20.8-23.3．水泵工作时单位时间的供水量大致为常数．水泵工作时单位时间的供水量大致为常数.4．流量是时间的连续函数．流量是时间的连续函数.5流量与水泵是否工作无关流量与水泵是否工作无关.6．流量．流量:单位时间流出的水的高度单位时间流出的水的高度*水塔截面积水塔截面积.1．拟合水位．拟合水位~时间函数时间函数.§2.4 §2.4 流量估计流量估计3．一天总用水量的估计．一天总用水量的估计.2．确定流量．确定流量~时间函数时间函数.§2.5 §2.5 算法设计与编程算法设计与编程1．拟合第．拟合第1、、2时段的水位，并导出流量时段的水位，并导出流量.2．拟合供水时段的流量．拟合供水时段的流量.3．一天总用水量的估计．一天总用水量的估计.4．流量及总用水量的检验．流量及总用水量的检验.Watertower.mn=(3,4) 的拟合流量曲线的拟合流量曲线n=(5,6)的拟合流量曲线的拟合流量曲线§2.6 §2.6 计算结果计算结果各时段和一天总用水量各时段和一天总用水量 及两个供水时段水泵的功率及两个供水时段水泵的功率(n1,n2)y1 y2 y12 y3 y p1 p2 (3,4)146.18 258.10 48.50 78.501263.4154.25 143.59 (5,6)146.52 257.76 46.13 76.301252.5153.06 142.67§2.7 §2.7 分析及改进分析及改进数据拟合、数值积分精度足够数据拟合、数值积分精度足够.流量曲线与原始记录基本上相吻合，零点到十点钟流流量曲线与原始记录基本上相吻合，零点到十点钟流量很低，十点到下午三点是用水高峰量很低，十点到下午三点是用水高峰.供水时段用供水时段用3次曲线拟合次曲线拟合4点不够好点不够好.全天平均流量全天平均流量22(cm/h).一天总用水量：一天总用水量：§3 §3 钢管订购和运输钢管订购和运输要铺设一条要铺设一条3.1 3.1 问题的提出问题的提出 经筛选后可以生产这种主管道钢管的钢厂有经筛选后可以生产这种主管道钢管的钢厂有的输送天然气的主管道的输送天然气的主管道, 如如图一图一所示。

所示A1325801010312012427010881070627030202030450104301750606194205201680480300220210420500600306195202720690520170690462160320160110290115011001200A2A3A4A5A6A7A8A9A10A11A12A13A14A15S1S2S3S4S5S6S7图一一铁路铁路钢管厂钢管厂 公路公路铁路铁路钢管厂钢管厂 公路公路天然气管道天然气管道火车站火车站里程（单位里程（单位km））1 12 23 34 45 56 67 780800 080800 010100 00 020200 00 020200 00 020200 00 030300 00 016160 015155 515155 516160 015155 515150 016160 01km主管道钢管称为主管道钢管称为1单位钢管单位钢管.一钢厂如果承担制造这种钢管一钢厂如果承担制造这种钢管,至少生产至少生产500个单位个单位.钢厂钢厂 在指定期限内能生产该钢管最大数量为在指定期限内能生产该钢管最大数量为 个单位个单位.钢管出厂销价钢管出厂销价1单位钢管为单位钢管为 万元万元.1 1单位钢管的铁路运价如表：单位钢管的铁路运价如表：里程里程(km)(km)≤300≤300301301～～350350 351351～～400400 401401～～450450 451451～～500500运价运价( (万元万元) )20202323262629293232里程里程(km)(km)501501～～600600601601～～700700 701701～～800800 801801～～900900901901～～10001000运价运价( (万元万元) )373744445050555560601000km1000km以上每增加以上每增加1 1至至100km100km运价增加运价增加5 5万元万元. . 公路运输费用为公路运输费用为1 1单位钢管每公里单位钢管每公里0.10.1万元（不足整公里万元（不足整公里部分按整公里计算）部分按整公里计算）. . 钢管可由铁路、公路运往铺设地点钢管可由铁路、公路运往铺设地点 和和管道全线管道全线. .☞☞ (1) (1)请制定一个主管道钢管的订购和运输计请制定一个主管道钢管的订购和运输计划，使总费用最小（给出总费用划，使总费用最小（给出总费用).).☞☞ (2) (2)请就请就(1)(1)的模型分析：哪个钢厂钢管的销的模型分析：哪个钢厂钢管的销价的变化对购运计划和总费用影响最大，哪个价的变化对购运计划和总费用影响最大，哪个钢厂钢管的产量的上限的变化对购运计划和总钢厂钢管的产量的上限的变化对购运计划和总费用的影响最大，并给出相应的数字结果费用的影响最大，并给出相应的数字结果. .☞☞ (3) (3)如果要铺设的管道不是一条线，而是一如果要铺设的管道不是一条线，而是一个树形图，铁路、公路和管道构成网络，请就个树形图，铁路、公路和管道构成网络，请就这种更一般的情形给出一种解决办法，并对这种更一般的情形给出一种解决办法，并对图图二二按按(1)(1)的要求给出模型和结果的要求给出模型和结果. . 问题问题 A13258010103120124270108810706270302020450104301750606194205201680480300220210420500600306195202720690520170690462160320160110290115011001200A19130190260100A2A3A4A5A6A7A8A9A10A11A12A13A14A15S1S2S3S4S5S6S7A16A17A18A20(A21)图二二303.2.1 3.2.1 模型的假设模型的假设3.2 3.2 模型的假设与符号约定模型的假设与符号约定铺设钢管各线段运输费以铺设钢管各线段运输费以1 1公里为单位公里为单位(1(1公里内不计公里内不计).). 钢管的铺设情况以从某一头向另一头铺设为准钢管的铺设情况以从某一头向另一头铺设为准. .承担生产任务的钢管厂，供应量至少为承担生产任务的钢管厂，供应量至少为500500个单位个单位. .3.2.2 3.2.2 符号约定符号约定第第i 段长度段长度. .第第i 段内的运费段内的运费. . i 段用段用j 厂单位钢管最小费用厂单位钢管最小费用. .i 段用段用j工厂钢管的数量工厂钢管的数量. .j工厂单位钢管生产价格工厂单位钢管生产价格. .j工厂生产钢管的上限工厂生产钢管的上限. .☞☞ (3) (3)要铺设一个树形管道网络图要铺设一个树形管道网络图, ,解决办法解决办法. .☞☞ (1) (1)制定钢管的订购和运输计划，使总费用最小制定钢管的订购和运输计划，使总费用最小. .铺设钢管的总费用包括三部分铺设钢管的总费用包括三部分 3.3 3.3 问题的分析问题的分析 ？？ ？？ ？？ ？？ ？？钢管生产成本钢管生产成本 运到各段的费用运到各段的费用 铺设时的运输费用铺设时的运输费用 问题：求各厂单位钢管运到各段的最小费用、最优运量问题：求各厂单位钢管运到各段的最小费用、最优运量☞☞ (2) (2)哪个厂销价的变化（产量上限的变化）对购运哪个厂销价的变化（产量上限的变化）对购运计划和总费用影响最大计划和总费用影响最大. .处理问题的方法与问题处理问题的方法与问题(1)(1)类似类似. .将各厂的最大生产上限和生产价格逐一进行变动，在将各厂的最大生产上限和生产价格逐一进行变动，在一厂变动时，其他厂不变一厂变动时，其他厂不变. .3.4 3.4 模型的建立与处理模型的建立与处理目标函数是总费用目标函数是总费用: :从工厂到各段的运费及其材料费之和从工厂到各段的运费及其材料费之和: :铺设时的运费铺设时的运费0.10.1为单位钢管的每公里运费为单位钢管的每公里运费其中其中为一定值为一定值设单位生产价格与单位运费相加之和为单位运费设单位生产价格与单位运费相加之和为单位运费则则目目标函数的限定条件函数的限定条件1)1)各厂运往各段的各厂运往各段的钢管量之和管量之和为此段此段铺设钢管的管的总量量. .2)2)各路段所需各厂的钢管量之和应小于此厂的最大生各路段所需各厂的钢管量之和应小于此厂的最大生产量，大于其最小生产量产量，大于其最小生产量. .3.4.1 3.4.1 问题一问题一的计算：的计算：的计算：的计算：的计算方法：观察法、枚举法和动态规划法的计算方法：观察法、枚举法和动态规划法 采用观察法和枚举法采用观察法和枚举法, ,比较得出工厂到各路段的单比较得出工厂到各路段的单位产品的最优费用位产品的最优费用单位位供供应量量 800 800 1000 500 854 717 500利用利用MATLABMATLAB解线性规划，得解线性规划，得考虑各厂均供应考虑各厂均供应各厂供应量各厂供应量得得由于由于生产量全部供应给了生产量全部供应给了为使全局更优为使全局更优计算可得计算可得各厂供应量各厂供应量单位位供供应量量 800 800 1000 500 767 1304 0同样可知同样可知 供应给了供应给了然而然而 的材料价格运输情况相对的材料价格运输情况相对 更高更高即使使用即使使用 代替代替 其供应量不会有影响其供应量不会有影响为使全局最优再次考虑为使全局最优再次考虑计算可得计算可得各厂供应量各厂供应量单位位供供应量量800 800 1000 0 1201 1370 0各厂供应量（即各厂供应量（即 （（ ）） 代表路段，代表路段， 代表供应厂）分配如下：代表供应厂）分配如下：供供应厂厂路段路段A1----A2 0 104 0 0 0 0 0A2----A3 0 125.7778 67.4444 0 107.7778 0 0A3----A4 272.5 184.6111 126.2778 0 166.6111 0 0A4----A5 128.5 184.6111 126.2778 0 166.6111 0 0A5----A6 194 0 0 0 0 0 0A6----A7 205 0 0 0 0 0 0A7----A8 0 201 0 0 0 0 0A8----A9 0 0 680 0 0 0 0A9----A10 0 0 0 0 240 240 0A10---A11 0 0 0 0 300 0 0A11---A12 0 0 0 0 220 0 0A12---A13 0 0 0 0 0 210 0A13---A14 0 0 0 0 0 420 0A14---A15 0 0 0 0 0 500 0 分别将各厂的生产价格与生产上分别将各厂的生产价格与生产上限变化一个单位，看他们对总费用的影响。

限变化一个单位，看他们对总费用的影响S变化化为1时10335300000P变化化为1时8008001000500854717500厂的销价变化对总费用的影响大厂的销价变化对总费用的影响大 3.4.2 3.4.2 问题二问题二厂的产量上限变化对购运计划总费用影响最大厂的产量上限变化对购运计划总费用影响最大评注评注程序计算可得：程序计算可得： =1170626.1+122387.2=1293013.3=1170626.1+122387.2=1293013.3 对于更一般的情况，如树形图，考虑情况可与对于更一般的情况，如树形图，考虑情况可与（（1 1）类似 在厂的取舍过程中亦同具体模型在厂的取舍过程中亦同具体模型与求解与求解 的最佳方案一致，从而再次重复的最佳方案一致，从而再次重复 的求的求解过程，可得：解过程，可得： =1296491.1; =128973.8. =1296491.1; =128973.8.3.4.3 3.4.3 问题三问题三 各厂供应量各厂供应量单位位供供应量量800 800 1000 0 1673 1630 0各厂供应量（即各厂供应量（即 （（ ）） 代表路段，代表路段， 代表供应厂）分配如下：代表供应厂）分配如下： 供供应厂厂路段路段A1----A2 0 104 0 0 0 0 0A2----A3 0 125.7778 53.4444 0 121.7778 0 0A3----A4 272.5 184.6111 22.2778 0 180.6111 0 0A4----A5 128.5 184.6111 22.2778 0 180.6111 0 0A5----A6 194 0 0 0 0 0 0A6----A7 205 0 0 0 0 0 0A7----A8 0 201 0 0 0 0 0A8----A9 0 0 680 0 0 0 0A9----A10 0 0 0 0 480 0 0A10---A11 0 0 0 0 300 0 0A11---A12 0 0 0 0 220 0 0A12---A13 0 0 0 0 0 210 0A13---A14 0 0 0 0 0 420 0A14---A15 0 0 0 0 0 500 0A15---A16 0 0 42 0 0 0 0A16---A17 0 0 0 0 0 10 0A17---A18 0 0 0 0 0 130 0A18---A19 0 0 0 0 190 0 0A19---A20 0 0 0 0 0 260 0A20---A21 0 0 0 0 0 100 03.5 3.5 模型的优点与改进模型的优点与改进 改进改进：在实际情况中，可考虑从某一段从两头：在实际情况中，可考虑从某一段从两头对铺，使模型更优化，对铺，使模型更优化， 对于供应量较多，运费和单对于供应量较多，运费和单价较低情况下应就近考虑价较低情况下应就近考虑(如：如： ) 。

在无在无软件帮助的情况下，可考虑一些不可能的情况，使软件帮助的情况下，可考虑一些不可能的情况，使计算规模缩小计算规模缩小 优点优点：本模型适用范围较广，在辅助软件的求：本模型适用范围较广，在辅助软件的求解下，结果易得出，管道数量的变化对算法几乎没解下，结果易得出，管道数量的变化对算法几乎没有影响。