《微积分模型》ppt课件

1、第五章第五章 微积分模型微积分模型例1：不允许缺货的存储模型设某厂消费假设干种产品，在轮换消费不同的产品时因改换设备要付消费预备费与产品数量无关，同一的产量大于需求时因占用仓库要付存储费。知某一产品日需求量为100件，消费预备费5000元，存储费每件每日1元，假设消费才干远大于需求，并且不允许出现缺货，试安排该产品的消费方案，即多少天消费一次消费周期，每次产量多少，可使总费用最小。 解：模型假设 1每天的需求量为常数 2每次消费预备费a，每天每件产品存储费 为b3.当存储量为零时，Q件产品立刻消费出来 模型建立: 设存储量为q(t)，当t=0时消费Q件，即q(0)= Q，以需求速率r递减，q(

2、T)=0，如图AQqrOTt一个周期内的存储费用为 ，其中积分等于图中三角形A的面积 ，一个周期内的预备费用a，那么总费用 Tdttqb 0 2QT2QTbac22rTba rTQ 那么每天的平均费用为2rTbTaTcc模型求解: 求T使c最小，那么braT2 有 barQ2最小的总费用为 abrc2例2报童的战略报童每天早上购进报纸零售，晚上将剩余的报纸退报答社，每份报纸的购进价为b元，零售价为a元，退回价为c元，cbx，那么x份全部售出，售出一份赚 元，退回一份赔 元，所以有 xgxr rxba cb xrxrrxfbarfrxcbrbaxg01问题为求x使 最大 xg由于需求量r和购进量

3、x都很大，将r看成延续的量更便于分析与计算，那么上式转化为 xxdrrxfbadrrfrxcbrbaxg0 xxdrrfcbxdrrfrcbrba00 xdrrfbax计算 的导数有 xg xfcbxdrrfcbxxfcbbadxxdgx0 xfbaxdrrfbax xxdrrfcbdrrfba0令 ，有 0dxxdg cbbadrrfbadrrfcbxx0由 上式可以化为 10drrf cabadrrfcbx0上式阐明报童应购进的报纸数量应该使卖不完与卖完的概率之比恰好等于卖出一份赚的钱与退回一份赔的钱之比。三级火箭的设计是最优的三级火箭的设计是最优的 发射卫星普通用三级火箭，而不是一级、发

4、射卫星普通用三级火箭，而不是一级、二级或四级，这是为什么呢？下面讨论运载二级或四级，这是为什么呢？下面讨论运载火箭的数学模型来阐明这个问题。火箭的数学模型来阐明这个问题。一火箭速度的计算一火箭速度的计算设火箭在设火箭在t时辰的质量为时辰的质量为 ，速度为，速度为 ，在在 时辰的质量为时辰的质量为 ，质量的变，质量的变化量为化量为 tm tvttttm totdtdmtmttm 1 质量 单调减少， ，减少的质量是火箭喷出的气体的质量。设火箭喷出的气体相对于火箭的速度为常数u，那么气体相对于地球的速度为 。由动量守恒定律 tm0dtdm utv ttvttmtvtm toutvtdtdm 2 将

5、1中 代入2，两边同除以 ，令 有ttmt0t udttdmdttdvtm 3 解方程得 tmmlnuvtv00 4 其中 是火箭的初始质量， 是t=0时的速度。阐明火箭的速度依赖于气体相对于火箭的速度u以及t时辰和t=0时的质量比。0m0v二火箭质量的计算 设火箭质量有三部分组成：载重质量 、燃料的质量 、构造质量 。 首先思索一切燃料全都耗尽，只剩下的质量为 ，pmfmsmspmm 有 spmmmlnuv0 5普通来说构造质量 在中 有一定比例，设为smsfmm psfsmmkmmkm0将上式代入5得pmkkmmlnuv100当载重质量 时，火箭的最大速度为0pmklnuv1知目前 ，假设

6、取 ，有 。这样不能用于卫星的发射。必需求改良。s/kmu310.k s/kmv7三改良的火箭质量采用多级发射，不断丢弃无用的部分就会提高效率。设在t到 时间内丢弃的质量看成1包括丢弃的构造质量和熄灭的燃料质量，那么熄灭掉的燃料质量为 ，由动量守恒定律ttk1 ttvttmtvtm utvtdtdmkttvdtdmk1将1中 代入上式，两边同除以 ttm ，令 有 t0t dttdmkdttdvtm 1解得 tmmlnuktv01 火箭的最终速度为 pmmlnuktv01那么理想火箭的最终速度可以到达s/km.510四火箭的最优级数 设 为第i级的燃料加构造质量， 为构造质量，那么燃料质量为

7、，并假设k和u都一样。 下面分析三级火箭的情形 初始质量为imikmimk13210mmmmmp当第一级燃料用完时，剩余质量为321mmkmmp由4火箭的速度为32101mmkmmmlnuvp此时丢弃一级外壳，启动二级，当第二级燃料用完时，速度323212mkmmmmmlnuvvpp当第三级燃料用完时，速度3323lnkmmmmuvvpp记 3201mmmmap3322mmmmmappppmmma33111111332211akaakaakalnuv由上式可知当 时， 最大，此时 321aaauvkckaaa1321uvexpc3最大的载重质量为 301kckmmp取 ， ， ，那么s/km.

8、v510s/kmu310.k 770pmm其他情况可类似给出，当u、v、k与前一样，载重1吨，那么火箭级数与质量之间的关系为级数12345质量149776560 飞越北极飞越北极 CMCM2000C今年今年6月，扬子晚报月，扬子晚报发布音讯：发布音讯：“中美航线下月可飞越北极，北中美航线下月可飞越北极，北京至底特律可节省京至底特律可节省4小时，摘要如下：小时，摘要如下： 7月月1日起，加拿大和俄罗斯将允许民航班日起，加拿大和俄罗斯将允许民航班机飞越北极，此改动可大幅度缩短北美与亚机飞越北极，此改动可大幅度缩短北美与亚洲间的飞行时间，旅客可直接从休斯敦，丹洲间的飞行时间，旅客可直接从休斯敦，丹佛

9、及明尼阿波利斯直飞北京等地。据加拿大佛及明尼阿波利斯直飞北京等地。据加拿大空中交通控制局估计，如飞越北极，底特律空中交通控制局估计，如飞越北极，底特律至北京的飞行时间可节省至北京的飞行时间可节省4个小时。由于不个小时。由于不需中途降落加油，实践节省的时间不止此数。需中途降落加油，实践节省的时间不止此数。问题重述：飞机飞行高度约为问题重述：飞机飞行高度约为10公里，飞行速公里，飞行速度约为每小时度约为每小时980公里；从北京至底特律原来的公里；从北京至底特律原来的航线飞经以下航线飞经以下10处：处：A1 (北纬北纬31度，东经度，东经122度度)； A2 (北纬北纬36 度，度，东经东经140度

10、度)；A3 (北纬北纬 53度，西经度，西经165度度)； A4 (北纬北纬62度，度，西经西经150度度); A5 (北纬北纬 59度，西经度，西经140度度)； A6 (北纬北纬 55度，度，西经度西经度)；A7 (北纬北纬 50度，西经度，西经130度度)； A8 (北纬北纬 47度，度，西经西经125度度)；A9 (北纬 47度，西经122度)； A10 (北纬 42度，西经87度)。请对“北京至底特律的飞行时间可节省4小时“从数学上作出一个合理的解释，分两种情况讨论：1.设地球是半径为6371千米的球体；2.设地球是一旋转椭球体，赤道半径为6378千米，子午线短半轴为6357千米。

11、模型的假设模型的假设1不思索地球的自转。不思索地球的自转。2飞机每经相邻两地的航程，均以曲面上两飞机每经相邻两地的航程，均以曲面上两点间最短间隔进展计算。点间最短间隔进展计算。3飞机飞行中途不需降落加油，同时忽略升飞机飞行中途不需降落加油，同时忽略升降时间。降时间。4开辟新航线后开辟新航线后,飞机由北京经过北极上空飞机由北京经过北极上空直飞底特律。直飞底特律。数据与符号的阐明数据与符号的阐明在以下计算中，北京的坐标为在以下计算中，北京的坐标为(北纬北纬40度，东度，东经经116度度)， 底特律的坐标为(北纬 43度，西经83度)。 球面(或椭球面)上的点Ai的直角坐标。 球面(或椭球面)上的点

12、Ai的经度与纬度。 t 飞机的飞行时间。 飞机飞行所节省的时间。 R 地球半径。 H 飞机飞行高度。 旋转椭球体的长半轴与短半轴。模型的建立与求解模型 地球是一个半径为R的均匀球体时的情况：先在曲面上建立直角坐标系，以地心为坐标原点O，以赤道平面为xoy平面，以0度经线即本初子午线圈所在的平面为xoz平面.那么,我们就可以iiizyx,ii,tba,写出球面的参数方程如下：coscosRx sincosRy sinRz 20，0球面上恣意两点不是一条直径的两个端点之间的最小间隔就是过这两点的大圆即经过球心的圆的劣弧长，根据这一方法，我们就可以确定恣意两点之间的飞行最短航线，我们可以知道:过A、

13、B两点的大圆的劣弧长即为两点的最短间隔，A、B两点的坐标分别为11111sin,cossin,coshRhRcoehRA点坐标B点坐标22222sin,cossin,coshRhRcoehR从A到B的飞行路程 OBOAOBOAhRlABarccos从A到B的飞行时间 980AnABlt球面上各点向量可表示为:iiiiiikjihROAsincossincoscos111111sincossincoscosiiiiiikjihROA那么的出详细结果单位：公里，“北京与A1之间的间隔是：1113.191576； A1与A2之间的间隔是：1758.78953；A2与A3之间的间隔是：4624.407

14、995；A3与A4之间的间隔是：1339.08267；A4与A5之间的间隔是：641.162553； A5与A6之间的间隔是：538.594244；A6与A7之间的间隔是：651.536106；A7与A8之间的间隔是：497.569491；A8与A9之间的间隔是：227.843664；A9与A10之间的间隔是：2810.85917；A10与“底特律之间的间隔是：331.924046；“北京直达“底特律的间隔是：10684.861107；“北京原到达“底特律的间隔是：14535 沿飞行原道路飞行总时间t=14.831600(小时)；飞机从北京直达底特律的航行时间t1=10.902912(小时)；

15、飞行节省的时间t=3.928681(小时)。地球是一旋转椭球体情况：地球是一旋转椭球体情况：假定椭球面上恣意两点假定椭球面上恣意两点A、B之间的测地线在经过之间的测地线在经过A、B两点的一个平面内两点的一个平面内.因此，我们可以写出过直线因此，我们可以写出过直线AB的平面族含一个参数的平面族含一个参数.下面就要求出过直线下面就要求出过直线AB的平面与椭球面的交线。可以利用搜索的方法求出的平面与椭球面的交线。可以利用搜索的方法求出交线上交线上A、B两点之间的劣弧长，然后求出弧长的最两点之间的劣弧长，然后求出弧长的最小值。小值。 采用近似计算的方法，运用数学类推的思想，在球体的根底上进展类推,算出

16、旋转椭球面上两点之间的最短间隔。 假设地球是一个旋转椭球体赤道半径为6378km，子午线短半轴为6357km此时A、B是地球上空距地面10km的恣意两点，坐标分别为 和 ，在x轴上任取一点 ，那么过A、B两点的椭圆曲线方程为111,zyxA222,zyxB0 , 0 ,mC1106357106387022222222111zyxzymxzymxzymx由此方程组可以解的： mxzzmxyy,那么A、B两点之间的弧长 211xxxxdxzyml求出 的最小值，即为AB的曲面最短道路长 ml由于上面的含参椭圆积分没有显式解析表达式，所以我们对模型作如下修正：首先我们察看从北京究竟特律中途所经城市的

17、纬度关系，不难发现，各城市的纬度之差都不会超越10,可以作这样的近似处置， 与 之间的曲面最短间隔 可由过 的平面上的弧 在赤道上的投影弧的长度来近似计算 。当每两个相邻城市之间的最短间隔算出来以后，就可以得到： iAjAijLjiAAjiAA原航线的总长度为:14458.4公里；飞行时间为:14.7535小时；飞越北极直达的航线长为：10605.8公里；直达时间为：10.8222小时；节省的时间为：3.9313小时。 模型的改良模型的改良前面我们利用空间曲线积分的思想求出了旋转椭球前面我们利用空间曲线积分的思想求出了旋转椭球面上的弧线长面上的弧线长,计算所得的结果精度较高计算所得的结果精度较

18、高,且结果比较且结果比较符合实践情况符合实践情况.但是对弧长的积分推导计算量大但是对弧长的积分推导计算量大,表达表达式很复杂式很复杂.我们可采用我们可采用“空间紧缩比率法更直观地空间紧缩比率法更直观地计算出大致的结果计算出大致的结果,详细改良方案如下详细改良方案如下:轴向紧缩比率法：轴向紧缩比率法：当我们假设地球是旋转椭球体时，由题知长半轴当我们假设地球是旋转椭球体时，由题知长半轴a=6371公里公里R(R为地球半径为地球半径)，而短半轴，而短半轴b=6357公公里，与里，与R有差距有差距.假设忽视假设忽视R与与a的微小差距，的微小差距， 我们假设椭球是由球体沿短半轴b紧缩而得到的，在球面上的

19、原来各点A0、A1、A2、A11紧缩后在椭球面上，且两点间的纬度差越大,对紧缩后的弧长影响也越大，同一纬度上的点，紧缩后弧长变化很小，由标题所给的条件可知：A0、A1、A2、A11各相邻点的纬度差分别为：8.9、5、7、9、3、4、5、3、0、5、0.4，北京与底特律的纬度差为3，差值都很小。由此可见这里用紧缩法是行之有效的。紧缩比率为定值 997806. 063716367hRhbk假设地球为球体时，已求出了航线的长度,在将这一结果乘以k，那么可得到为椭球体时的航线长，下面是利用这一方法得出的结果:分段飞行时：总航线的长度L=14534.96114公里，总的飞行时间t=14.799小时；从北

20、京直达底特律的航线长度L1=10657.6公里，飞行时间t1=10.872522小时，那么缩短的路程 公里，节省的时间 小时。36.3877l9265. 3t 变分法变分法简介 设函数 知，具有二阶延续偏导数，求 函数，使zyxF, xyy 21,xxdxxyxyxFyJ到达极大小值，上式表示的泛函称为最简泛函 1.端点固定的情形 设 ，那么对于满足 的二阶可导函数 和绝对值充分小的恣意常数 ，有 minyJ 021xx x 0yJyJ而 yJyJ yJdxyyxFxx21,在 处展开 0 yJdxFFyyxFxxyy212, 221xxyydxFF从而有 021xxyydxFF又 212121xxyxxyxxydxdxdFFdxF21xxydxdxdF所以有 021xxyydxdxdFF由 的恣意性，得 x0.,yyxFdxdyy