六章分方程模型

第六章第六章 微分方程模型微分方程模型一、经济增长模型一、经济增长模型 发展经济、提高生产力主要有以下手段：增加投资、发展经济、提高生产力主要有以下手段：增加投资、增加劳动力、技术革新增加劳动力、技术革新.本节的模型将首先建立产值与资金、劳动力之间的关本节的模型将首先建立产值与资金、劳动力之间的关系，然后再研究资金与劳动力的最佳分配，使投资效益系，然后再研究资金与劳动力的最佳分配，使投资效益最大，最后讨论如何调节资金与劳动力的增长率，使劳最大，最后讨论如何调节资金与劳动力的增长率，使劳动生产率得到有效的增长动生产率得到有效的增长.1.Douglas生产函数生产函数 用用 分别表示某一地区或部门在时分别表示某一地区或部门在时刻刻 的产值、资金和劳动力，相互的关系为的产值、资金和劳动力，相互的关系为 ,Q tK tL tt ,Q tF K tL t其中其中 为待定函数，对于固定的时刻为待定函数，对于固定的时刻 上面的函数简上面的函数简记为记为F,t,.QF K L 现来探讨函数现来探讨函数 的具体表达式，引入记号的具体表达式，引入记号F,QKzyLL分别表示每个劳动力的产值和投资。有如下的假设：分别表示每个劳动力的产值和投资。有如下的假设：随着随着 的增加而增加，但增长速度递减。从而可以假设的增加而增加，但增长速度递减。从而可以假设为为yz01,.zcy0.511.522.530.20.40.60.811.21.4函数函数 满足上面的要求，常数满足上面的要求，常数 可看成是可看成是技术的作用。将上面的形式代入到技术的作用。将上面的形式代入到中，即：中，即：g yy0c 1.KQzLcy LcLcK LL由由式知函数式知函数 有如下的性质：有如下的性质：Q2222,0,0.QQQQKLKL式的具体意义是：产值是资金和劳动力的递增函数，式的具体意义是：产值是资金和劳动力的递增函数，但增长率逐渐下降（即加速度为负数）但增长率逐渐下降（即加速度为负数）.（见前图）（见前图）记记 表示单位资金创造的产值（又称为对表示单位资金创造的产值（又称为对,kkQQQK,LLQQQL资金的边际产值）；资金的边际产值）；表示单位劳动力能创造表示单位劳动力能创造的产值（又称为对劳动力的边际产值）的产值（又称为对劳动力的边际产值）.则从则从式得到式得到,1.kLKLKQLQKQLQQQQ 的具体意义是：的具体意义是：是资金在产值中占有的份额，是资金在产值中占有的份额，是劳动力在产值中占有的份额。所以是劳动力在产值中占有的份额。所以 的大小直接反映的大小直接反映了资金、劳动力二者对于创造产值的轻重关系。了资金、劳动力二者对于创造产值的轻重关系。1 式是经济学中著名的式是经济学中著名的Cobb-Douglas生产函数，更生产函数，更一般的形式是一般的形式是,0,1.QcK L 2.资金与劳动力的最佳分配资金与劳动力的最佳分配 本段根据本段根据式讨论，如何分配资金与劳动力，使生产式讨论，如何分配资金与劳动力，使生产创造的效益达到最大。创造的效益达到最大。假定资金来自贷款，利率为假定资金来自贷款，利率为 每个劳动力都要支付工每个劳动力都要支付工资资 因而总效益为因而总效益为,r,w则相应的问题转化为资金与劳动力的分配比例则相应的问题转化为资金与劳动力的分配比例 （即每个劳动力占有的资金），使效益（即每个劳动力占有的资金），使效益 为最大。为最大。/K LS.SQrKwL此问题由微分法可得到（为一个求极值的问题）。在此问题由微分法可得到（为一个求极值的问题）。在式两边对式两边对 求导，并令其为零，则有求导，并令其为零，则有.KLQrQw再由再由式；可得到式；可得到此即为资金与劳动力的最佳分配。此即为资金与劳动力的最佳分配。.1KwLr,K L0,0.KKLLSQrSQw例例 取取则由关系则由关系1000000,100,800,0.004,KLwr,1KwLr得得0.05,0.0476.1 3.劳动生产率增长的条件劳动生产率增长的条件 衡量经济增长的指标：总产值衡量经济增长的指标：总产值 和每个劳动力的产值和每个劳动力的产值 这个模型讨论的就是这个模型讨论的就是 满满足什么条件才能使足什么条件才能使 保持增长？保持增长？Q /.z tQ tL t ,K tL t ,Q tz t 假设假设1.投资增长率与产值成正比，比例系数投资增长率与产值成正比，比例系数0.2.劳动力的相对增长率为常数劳动力的相对增长率为常数.注：这两个条件的数学表达式分别为：注：这两个条件的数学表达式分别为：.0.dKQdt.dLLdt方程方程的解为的解为 0.tL tL e将将 代入代入 得到得到1,zcyQcK L,dKQdt,dKc Lydt又因又因 再由再由式，两边求导式，两边求导,得得,KLy.dKdyLLydtdt比较上面两个式子，就有比较上面两个式子，就有.dyyc ydt此方程为此方程为Bernoulli方程，其解为方程，其解为 1110011.tKcy teK 其中其中00.tdKKdt ,Q tz t 以下根据以下根据式来研究式来研究 保持增长的条件。保持增长的条件。1)增长增长 Q t 由条件由条件 及及 得得0,dQdt,QcLy1dQdycL yc Lydtdt2111,cLycyy其中的其中的 以以代入，可知条件代入，可知条件 等价于等价于0dQdt10011.1tKeK 注意到上式右端大于注意到上式右端大于 所以当所以当 （即劳动力不减少）（即劳动力不减少）时，上式总成立；而当时，上式总成立；而当 时，上式成立的条件是时，上式成立的条件是1,00此说明在劳动力减少的条件下，产值只能在短时间内增此说明在劳动力减少的条件下，产值只能在短时间内增长，同时注意到，若长，同时注意到，若001ln11.1KtK 00111,KK则不存在这样的时间。则不存在这样的时间。2)增长增长 z t 由条件得由条件得 再由再由 知知0,dzdt,zcy00.dzdydtdt由于：由于：.dyyc ydt所以当所以当 时，该条件成立，而当时，该条件成立，而当 时，时，00100010,tKdyedtK 此条件等价于此条件等价于00.KK此条件的意义是此条件的意义是:劳动力增长率小于初始投资增长率劳动力增长率小于初始投资增长率.二、香烟过滤嘴的作用二、香烟过滤嘴的作用 问题问题 烟草公司普遍地在香烟尾部装上一截过滤嘴，烟草公司普遍地在香烟尾部装上一截过滤嘴，但是过滤嘴的作用有多大，和使用的材料有什么直接的但是过滤嘴的作用有多大，和使用的材料有什么直接的关系？在这一段中，我们建立相应的数学模型来具体讨关系？在这一段中，我们建立相应的数学模型来具体讨论这个问题。论这个问题。分析分析 吸烟时毒物吸入人体的过程大致为：由于毒物基本上吸烟时毒物吸入人体的过程大致为：由于毒物基本上均匀地分布在烟草中，吸烟时点燃的烟草大部分化为烟均匀地分布在烟草中，吸烟时点燃的烟草大部分化为烟雾，毒物由烟雾携带一部分直接进入空气，另一部分沿雾，毒物由烟雾携带一部分直接进入空气，另一部分沿香烟穿行，在穿行过程中又部分地被未点燃的烟草和过香烟穿行，在穿行过程中又部分地被未点燃的烟草和过滤嘴吸收而沉淀下来，剩下的进入人体。被烟草吸收而滤嘴吸收而沉淀下来，剩下的进入人体。被烟草吸收而沉淀下来的那一部分毒物，当香烟燃烧到那里的时候又沉淀下来的那一部分毒物，当香烟燃烧到那里的时候又通过烟雾部分进入空气，部分沿香烟穿行，这个过程一通过烟雾部分进入空气，部分沿香烟穿行，这个过程一直到香烟燃烧到过滤嘴处为止。直到香烟燃烧到过滤嘴处为止。在整个过程中，原来分布在烟草中的毒物除了进入空在整个过程中，原来分布在烟草中的毒物除了进入空气和被过滤嘴吸收的一部分外，剩下的全部被人体吸入气和被过滤嘴吸收的一部分外，剩下的全部被人体吸入.模型假设模型假设 1.烟草和过滤嘴的长度分别是烟草和过滤嘴的长度分别是 ，毒素，毒素 （毫克）（毫克）均匀地分布在烟草中，密度为均匀地分布在烟草中，密度为1,l lM01/.wM l 2.点燃处的毒物随烟雾进入空气和沿烟草穿行的数量点燃处的毒物随烟雾进入空气和沿烟草穿行的数量比例是比例是:,1.3.未点燃的烟草和过滤嘴对烟雾穿行的毒物的吸收率未点燃的烟草和过滤嘴对烟雾穿行的毒物的吸收率分别是常数分别是常数,.b 4.烟草沿香烟穿行的速度是常数烟草沿香烟穿行的速度是常数 香烟的燃烧速度香烟的燃烧速度是常数是常数 且且 ,v.u.vu 将一支烟吸完后毒物进入人体的总量记为将一支烟吸完后毒物进入人体的总量记为 常识告常识告诉我们，量诉我们，量 与过滤嘴的吸收率、烟草中毒物的初始含与过滤嘴的吸收率、烟草中毒物的初始含量有关，因此降低量有关，因此降低 的含量及降低烟雾在香烟总的穿的含量及降低烟雾在香烟总的穿行速度都是降低吸入量的有效方法。行速度都是降低吸入量的有效方法。,QQQ 模型建立模型建立xxx x1ll q xq xx 设在时刻设在时刻 及及 处点燃香烟，为此建立相应处点燃香烟，为此建立相应的坐标系统：吸入毒物量的坐标系统：吸入毒物量 由毒物穿过香烟的流量确由毒物穿过香烟的流量确定，后者又与毒物在烟草中的密度有关，为此定义函定，后者又与毒物在烟草中的密度有关，为此定义函数：数：0t 0 x Q 毒物流量毒物流量 表示在时表示在时刻刻 时单位时间通过香烟截面时单位时间通过香烟截面 处的毒物量；处的毒物量；,q x ttx0 xl 毒物密度毒物密度 表示时刻表示时刻 截面截面 处单位长度烟草处单位长度烟草中的毒物含量。中的毒物含量。,w x ttx 如果知道了流量函数如果知道了流量函数 ，吸入的毒物流量，吸入的毒物流量 就就是是 处的流量在吸一支烟时间内的总和。即处的流量在吸一支烟时间内的总和。即,q x tQxl下面分下面分4个过程来计算个过程来计算 值。值。Q10,/.TQq x t dtTlu 1.在在 瞬间由烟雾携带的毒物在单位时间内通过瞬间由烟雾携带的毒物在单位时间内通过处的数量处的数量 由假设由假设4 知在香烟点燃处知在香烟点燃处 静止不动。静止不动。0t x,0,q x,vu0 x 记记 考察考察 所对应的一段所对应的一段香烟，毒物通过香烟，毒物通过 和和 处的流量分别是处的流量分别是 和和 则根据守恒定律，这两个流量之差等于这则根据守恒定律，这两个流量之差等于这一段未点燃的烟草或过滤嘴的毒物的吸收量，由假设一段未点燃的烟草或过滤嘴的毒物的吸收量，由假设2及及4，有，有,0,q xq x,x xx xxx q x,q xx 11 0,bq xxlq xq xxq xlxl 其中其中 是烟雾穿过长度为是烟雾穿过长度为 的这一段香烟所需的这一段香烟所需要的时间，令要的时间，令 则有微分方程：则有微分方程：xvx0,101 0lim .xbq xxlq xxq xdqvbdxxq xlxlv 在在 处点燃的烟草单位时间内放出的毒物量记作处点燃的烟草单位时间内放出的毒物量记作0 x 根据假设根据假设1、3、4，则可得微分方程，则可得微分方程2的初始条件为的初始条件为求解微分方程求解微分方程。首先由方程及初始条件。首先由方程及初始条件 00,dqbq xdxvqaH 方程的通解为方程的通解为 再由初始条件得再由初始条件得从而得到方程的解为从而得到方程的解为,bxvq xCe0.CaH 0000,.qaHHuw 0,0.bxvqq xaH exl对于方程及初始条件对于方程及初始条件 110,blvdqq xdxvq laH e 同样可得方程的解为同样可得方程的解为 1101,.x lblvvq xaH eelxl即方程的解可以表达为即方程的解可以表达为 110101 0,.bxvx lblvvaH exlq xaH eelxl 2.在香烟燃烧的任何时刻在香烟燃烧的任何时刻 求毒物在单位时间内通求毒物在单位时间内通过过 的数量的数量,txl,.q l tt 因为在时刻因为在时刻 时，香烟燃至时，香烟燃至 处，记此时点燃的处，记此时点燃的烟草在单位时间放出的毒物量为烟草在单位时间放出的毒物量为 则则xut,H t仿仿1可得可得,.H tuw ut t 1111 0,.b x utvb lutx lvvaH t exlq x taH t eelxl 3.确定确定,w x t 由于在吸烟过程中，未点燃的烟草不断得吸收烟雾中由于在吸烟过程中，未点燃的烟草不断得吸收烟雾中的毒物，所以毒物在烟草中的密度的毒物，所以毒物在烟草中的密度 由初始值由初始值逐渐增加。考察烟草截面逐渐增加。考察烟草截面 处在处在 时间内毒物密度的增时间内毒物密度的增量量 根据守恒定律它等于在单根据守恒定律它等于在单位长度烟雾中的毒物被吸收的部分，由假设位长度烟雾中的毒物被吸收的部分，由假设4，则有，则有,w x t0wtx,w x xtw x t,q x tw x ttw x tbtv 令令 并将并将、代入得方程代入得方程0,t 0,0.b x utvwabuw ut t etvw xw此方程为偏微分方程，该方程的解为此方程为偏微分方程，该方程的解为00,1.,1bxbutabutvvvabutvaw x tweeeaww ut taea其中其中 （由假设（由假设2）1.aa 4.计算计算Q 将将代入代入中，得中，得120,.blblbutabutvvvvauwq l teeeea最后将最后将代入代入式作定积分得：式作定积分得：21100,1.blabllvvauwQq l t dteea b 为便于对结果的分析，将上式写成为便于对结果的分析，将上式写成1211,allvveQaMea blv记记 则则式可写成式可写成 11,ra blerrvr 2.lvQaMer式即是我们需要的最终，该式表达了吸入量式即是我们需要的最终，该式表达了吸入量 与与 等因素之间的关系。等因素之间的关系。Q,a M21,l v b l 结果分析结果分析 1.与烟草含毒物量与烟草含毒物量 、毒物随烟雾沿香烟穿行比例、毒物随烟雾沿香烟穿行比例 成正比。可以设想为将毒物成正比。可以设想为将毒物 集中在集中在 处，则吸处，则吸入量为入量为QMaMxl.aM 2.因子因子 体现了过滤嘴对减少毒物进入人体的作体现了过滤嘴对减少毒物进入人体的作用，提高过滤嘴的长度用，提高过滤嘴的长度 及吸收率及吸收率 降低穿行速度降低穿行速度也可降低吸入量。也可降低吸入量。2lve2,l.v r 3.因子因子 表示的是未点燃的烟草对毒物的吸收而起表示的是未点燃的烟草对毒物的吸收而起到的减少到的减少 的作用。的作用。Q211.2lva blQaMev 根据实际数据有根据实际数据有 则则11,a blrv 4.比较比较 为了比较过滤嘴的作用，取两支香烟作比较，两支香为了比较过滤嘴的作用，取两支香烟作比较，两支香烟的长度均为烟的长度均为 一支带过滤嘴，长度为一支带过滤嘴，长度为 吸入量分别吸入量分别为为 则有则有.l2.l12,.Q Q212.b lvQeQ可见只要可见只要 过滤嘴就体现了相应的作用。过滤嘴就体现了相应的作用。b三、烟雾的扩散和消失三、烟雾的扩散和消失 问题的提出问题的提出 一颗炮弹在平原上爆炸，放出的烟雾以爆炸点为中心一颗炮弹在平原上爆炸，放出的烟雾以爆炸点为中心向四周迅速扩散，形成一个近似于圆形的不透光区域。向四周迅速扩散，形成一个近似于圆形的不透光区域。由于这个这个区域逐渐增大，其边界逐渐明亮起来，不由于这个这个区域逐渐增大，其边界逐渐明亮起来，不透光区域逐渐变小，最后烟雾完全消失。本节建立一个透光区域逐渐变小，最后烟雾完全消失。本节建立一个相应的数学模型来描述烟雾的扩散和消失的过程，分析相应的数学模型来描述烟雾的扩散和消失的过程，分析消失的时间与哪些因素有关。消失的时间与哪些因素有关。问题的分析问题的分析 爆炸引起的烟雾传播可以看成是无穷空间由瞬时点源爆炸引起的烟雾传播可以看成是无穷空间由瞬时点源导致的扩散过程。能够用二阶抛物型偏微分方程描述烟导致的扩散过程。能够用二阶抛物型偏微分方程描述烟雾浓度的变化规律。雾浓度的变化规律。整个建模过程应当包括能够：烟雾浓度的变化规律；整个建模过程应当包括能够：烟雾浓度的变化规律；穿过烟雾的光的强度的变化规律，仪器辨别亮暗的灵敏穿过烟雾的光的强度的变化规律，仪器辨别亮暗的灵敏度的描述；不透光区域边界的变化过程等。度的描述；不透光区域边界的变化过程等。模型假设模型假设 1.炮弹的爆炸看作是在空中某一点向四周等强度地瞬炮弹的爆炸看作是在空中某一点向四周等强度地瞬时释放烟雾，烟雾在空间扩散，不计风力和大地的影时释放烟雾，烟雾在空间扩散，不计风力和大地的影响；响；2.烟雾的传播遵从扩散定律：即单位时间通过法向面烟雾的传播遵从扩散定律：即单位时间通过法向面积的流量与它的浓度成正比；积的流量与它的浓度成正比；3.光线穿过烟雾时其强度由于烟雾的吸收而减少光线穿过烟雾时其强度由于烟雾的吸收而减少,单位单位距离上光强的相对减少量与烟雾浓度成正比距离上光强的相对减少量与烟雾浓度成正比,没有烟雾的没有烟雾的大气对光线的吸收作用忽略不记；大气对光线的吸收作用忽略不记；4.在烟雾的扩散过程中，不穿过烟雾直接进入观测仪在烟雾的扩散过程中，不穿过烟雾直接进入观测仪器的标准光强器的标准光强 保持不变，对于穿过烟雾而进入仪器的保持不变，对于穿过烟雾而进入仪器的光强光强 观测结果只有亮暗之分，仅当观测结果只有亮暗之分，仅当 时观测时观测结果为量。结果为量。称为仪器的灵敏度。称为仪器的灵敏度。0I,I0/1I I 1 建模建模 将爆炸时刻记为将爆炸时刻记为 爆炸点设为坐标原点。在时刻爆炸点设为坐标原点。在时刻 时空间中任一点时空间中任一点 的烟雾浓度记为的烟雾浓度记为由假设由假设2，单位时间通过单位法向面积的流量，单位时间通过单位法向面积的流量0,t t,x y z,.C x y z tgrad.qkC 其中其中 是扩散系数，负号表示由浓度高向浓度低的地方是扩散系数，负号表示由浓度高向浓度低的地方扩散。设空间区域扩散。设空间区域 的体积为的体积为 包围包围 的曲面为的曲面为 的外法向量为的外法向量为 则在时间区间则在时间区间 内，通过内，通过 的的k,V,S,nS,t tt 的流量为的流量为1.tttSQq nd dt 而而 的烟雾增量为的烟雾增量为2,.QC x y z tC x y z ttdV 由质量守恒定律，知由质量守恒定律，知12.QQ由高斯公式，有由高斯公式，有.Sq nddivqdV 由积分中值定理得：由积分中值定理得：222222div grad,CCCCkCktxyz0,.tx y z 再由假设再由假设1，初始条件为在坐标原点的点源函数，记为，初始条件为在坐标原点的点源函数，记为,0,.C x y zQx y z其中其中 表示烟雾总量，表示烟雾总量，是单位强度的点源函是单位强度的点源函数。数。Q,x y z 方程方程的解为的解为22243/2,.4xyzktQC x y z tek t方程的解说明：在任何时刻方程的解说明：在任何时刻 烟雾浓度函数烟雾浓度函数 的等值面的等值面为球面为球面 并且当并且当 时，时，,tC2222.xyzRR 0.C 2.穿过烟雾的光强的变化规律穿过烟雾的光强的变化规律 考察沿一定方向穿过烟雾的光线，此方向的长度坐标考察沿一定方向穿过烟雾的光线，此方向的长度坐标记为记为 烟雾浓度为烟雾浓度为 光强为光强为 由假设由假设3，则有，则有,l,I l,C l .dIaC l I ldl 是烟雾对光线的吸收系数，光线未进入烟雾时是烟雾对光线的吸收系数，光线未进入烟雾时的强度记作的强度记作 即即a0ll0,I 00.II l则方程则方程在条件在条件下的解为下的解为 00.llC s dsI lI e 3.仪器的灵敏度与不透光区域的边界仪器的灵敏度与不透光区域的边界 注意到，烟雾浓度在空间中是连续变化的，穿过烟雾注意到，烟雾浓度在空间中是连续变化的，穿过烟雾而进入仪器的光强也是连续变化的，之所以会观察到烟而进入仪器的光强也是连续变化的，之所以会观察到烟雾扩散时不透光区域的边界有一个先大后小、最终消失雾扩散时不透光区域的边界有一个先大后小、最终消失的过程，是由于仪器的观测结果只有亮暗之分。亮暗分的过程，是由于仪器的观测结果只有亮暗之分。亮暗分界线由灵敏度界线由灵敏度 决定。由假设决定。由假设4仅当仅当01,II 观测结果为亮。观测结果为亮。由假设由假设3，光强，光强 穿过没有烟雾的大气时其衰减可以穿过没有烟雾的大气时其衰减可以忽略，因而不必对它与直接进入仪器的标准光强加以区忽略，因而不必对它与直接进入仪器的标准光强加以区分，从而分，从而式中的光强可以由式中的光强可以由式进行计算。式进行计算。0I 为方便起见，取沿着为方便起见，取沿着 轴的光线，不妨设光源（太轴的光线，不妨设光源（太阳）在阳）在 处而仪器在处而仪器在 处，则处，则可以写可以写成成zz z ,1.aC x y z t dse 因为因为 的等值面为球面，所以仪器观测到的的等值面为球面，所以仪器观测到的投影在投影在 平面上的不透光区域的边界是圆周，记作平面上的不透光区域的边界是圆周，记作,C x y z txy222.xyr其中圆半径由其中圆半径由式确定式确定.4.不透光区域的边界的变化规律不透光区域的边界的变化规律 对对式取对数，并利用一次近似式取对数，并利用一次近似则则式可化为式可化为ln 11xx x11,ln.1C x y z t dzaa于是不透光区域的边界由于是不透光区域的边界由,C x y z t dza确定，将确定，将式中的式中的 代入上式并进行积分代入上式并进行积分,并并,C x y z t利用公式利用公式 可得可得2,xaedxa224,4xyQekta以以 代入上式代入上式,有有222.xyr 4ln.4aQr tktk t 结果分析结果分析 由由式可大致画出函数式可大致画出函数 的曲线的曲线,并且能算出当并且能算出当 r t14aQttke时不透光的区域的半径达到时不透光的区域的半径达到最大最大 当当.mrt r t1tmrO2t24aQttk时时,即此时烟雾完全消散即此时烟雾完全消散.0.r ,式表明式表明 与烟雾施放量与烟雾施放量 和烟雾对光线的吸和烟雾对光线的吸收系数收系数 成正比成正比,与扩散系数与扩散系数 和系统灵敏度和系统灵敏度 成反比成反比.12,t tQak 再从最后两式可以得到再从最后两式可以得到2112.7.ttet此说明当知道此说明当知道 后后,可预测烟雾完全消失的时刻可预测烟雾完全消失的时刻1t2.t四、万有引力的发现四、万有引力的发现 历史背景历史背景 15世纪下半叶开始世纪下半叶开始,欧洲商品经济的繁荣促进了航海欧洲商品经济的繁荣促进了航海事业的发展事业的发展.哥伦布新大陆的发现哥伦布新大陆的发现,麦哲伦的环球远航麦哲伦的环球远航,引起了社会的普遍关注引起了社会的普遍关注.当时的远洋航船的方位全考星当时的远洋航船的方位全考星球的位置来确定球的位置来确定.在强大的社会需要推动下在强大的社会需要推动下,天文观测天文观测的精确程度不断提高的精确程度不断提高.在大量的实际观测数据面前在大量的实际观测数据面前,一一直处于天文学统治地位的直处于天文学统治地位的“地心说地心说”开始摇动了开始摇动了.波兰天文学家哥白尼在天文观测的基础上波兰天文学家哥白尼在天文观测的基础上,冲破宗教冲破宗教统治和统治和“地心说地心说”的束缚的束缚,提出提出“日心说日心说”.这是天文这是天文学乃学乃整个科学的一大革命整个科学的一大革命.但是由于历史条件和科学水平的但是由于历史条件和科学水平的限制限制,哥白尼的理论还有些缺陷哥白尼的理论还有些缺陷.他接受了圆周运动是他接受了圆周运动是最完美的天体运动形式的概念最完美的天体运动形式的概念,认为行星绕太阳的运行认为行星绕太阳的运行轨道是圆形的轨道是圆形的.意大利物理学家伽利略不仅用观察方法证明了哥白尼意大利物理学家伽利略不仅用观察方法证明了哥白尼的学说的学说,而且用实验方法发现了落体定律和惯性原理而且用实验方法发现了落体定律和惯性原理,揭示了物体在不受阻扰时作匀速直线运动的规律揭示了物体在不受阻扰时作匀速直线运动的规律.德国天文学家、数学家开普勒在第谷德国天文学家、数学家开普勒在第谷.布拉赫对于行布拉赫对于行星运动大量观测资料的基础上用数学方法研究发现星运动大量观测资料的基础上用数学方法研究发现,火火星的实际位置与按哥白尼理论计算的位置相差星的实际位置与按哥白尼理论计算的位置相差8弧分弧分.经过对观测数据长期深入的分析经过对观测数据长期深入的分析,开普勒终于归纳出著开普勒终于归纳出著名的所谓行星运动三大定律名的所谓行星运动三大定律:即即 各颗行星分别在不同的椭圆轨道上绕太阳运行各颗行星分别在不同的椭圆轨道上绕太阳运行,太阳太阳位于这些椭圆的一个焦点上位于这些椭圆的一个焦点上;每颗行星运行过程中单位时间内太阳每颗行星运行过程中单位时间内太阳行星向径扫过行星向径扫过的面积是一个常数的面积是一个常数;各颗行星运行的周期的平方与其轨道长半轴的各颗行星运行的周期的平方与其轨道长半轴的3次方次方成正比成正比.在伽利略、开普勒的基础上在伽利略、开普勒的基础上,17,18世纪许多科学家致世纪许多科学家致力于行星沿椭圆轨道运行时受力状况的研究力于行星沿椭圆轨道运行时受力状况的研究.从开普勒从开普勒定律中可以看出定律中可以看出,行星运行速度是变化的行星运行速度是变化的.而在当时没而在当时没有计算变速运动的动力学方法有计算变速运动的动力学方法.英国物理学家胡克和荷英国物理学家胡克和荷兰物理学家惠更斯等人虽然都取得了一些成果兰物理学家惠更斯等人虽然都取得了一些成果,但终未但终未得到有关引力的定律得到有关引力的定律.英国物理学家、数学家牛顿认为一切运动都有其力学英国物理学家、数学家牛顿认为一切运动都有其力学原因原因.开普勒三定律的背后必定有力学定律起作用开普勒三定律的背后必定有力学定律起作用.他他在研究变速运动中发明了微积分（流数法）在研究变速运动中发明了微积分（流数法）,又以微积又以微积分为工具在开普勒三定律和牛顿第二定律的基础上用演分为工具在开普勒三定律和牛顿第二定律的基础上用演绎的得到所谓的万有引力定律绎的得到所谓的万有引力定律.模型假设模型假设 开普勒三定律和牛顿第二定律为真开普勒三定律和牛顿第二定律为真;对任意一颗行星的椭圆运行轨道建立极坐标对任意一颗行星的椭圆运行轨道建立极坐标 并并以太阳为坐标原点以太阳为坐标原点,以椭圆长半轴方向为极轴正向以椭圆长半轴方向为极轴正向.,.rrOrp 1.轨道方程轨道方程.1cospre2222,1.bpbaea其中其中 为椭圆的长、短半轴为椭圆的长、短半轴,为离心率为离心率.,a be 2.单位时间向径单位时间向径 扫过的面积是常数扫过的面积是常数 即即r.A21.2rA 3.行星运行周期行星运行周期 满足满足T23.Ta其中其中 是绝对常数是绝对常数,与行星无关与行星无关.4.行星运行时所承受的作用力行星运行时所承受的作用力 等于行星加速度等于行星加速度 和和fr质量质量 的乘积的乘积.即即m.fmr 模型建立模型建立 万有引力定律研究的是力万有引力定律研究的是力 的大小和方向的大小和方向.由由知知,它它取决于行星的运动加速度取决于行星的运动加速度 由此必须研究行星的位置由此必须研究行星的位置和速度和速度 为此先引入基向量为此先引入基向量f.rr.rcossin,sincos.ruijuij rOrpruu 向径向径 可表示为可表示为r.rrru000,.tKdyyc yydtL作变换：作变换：则原方程变形为则原方程变形为1,zy11.zzc 此方程为一阶线性微分方程，代入求解公式得此方程为一阶线性微分方程，代入求解公式得11111.tttcczeeCC e 代入代入 则方程的通解为则方程的通解为1,zy111.tcyC e 再由初始条件，得再由初始条件，得11001100,KKccCCLL 1110011.tKcy teK