常微分方程4.5

《常微分方程4.5》由会员分享，可在线阅读，更多相关《常微分方程4.5（27页珍藏版）》请在装配图网上搜索。

1、4.3 高阶方程的降阶和幂级数解法一般的高阶方程没有普遍适用的解法.处理问题的原则是降阶.4.3.1可降阶的一些方程的类型n阶方程的一般形式为 F(t,x,x,x(n)=0.这里介绍几类可降阶的方程：不显含未知函数x，或一般地，不含x及x的各个较低阶导数的方程F(t,x(k),x(k+1),x(n)=0.对此类方程，令x(k)=y,即可将方程降k阶例x(n)=f(t).解可用逐次积分法求通解显然，这个解满足初始条件x(t0)=cn,x(t0)=cn-1,x(n-1)(t0)=c1.为计算方便，可把关于t的n次求积分的表达式,化为对参数u的一次积分n=2时，右端看作vou平面上的二重积分，积分区

2、域为图中三角形交换积分次序，得n=3时，类似地对F(t,x(n)=0若能解出x(n)，则按例方法求解，否则引入参数u，将其改写为参数形式t=(u),x(n)=(u).例2解设x=p，则原方程等价于例解设x()y,则所以不显含自变量t的方程 F(x,x,x(n)=0.令x=y，则原方程可化为以y为未知函数，x为自变量的n-1阶方程事实上，有例解设x=y，则 代入原方程得于是例求数学摆的运动方程满足初值条件：t=0时,的解.解设则初值条件代入得考虑摆从最大的正偏离角=0到最大的负偏离角=-0间的第一次摆动的情况此时0,所以应讨论下述方程积分上式,并注意初值条件是t=2t0时=-0.方程(3)的解在

3、2t0 t 4t0上适用,4t0 t 6t0上摆的运动又由方程(1)描述,.注1 不能用初等函数表示,是一个椭圆积分椭圆积分(形如 或 的积分,最先在求椭圆弧长时遇到此类积分).注2 显然,上述讨论的非线性情形比线性化了的情形要复杂.齐线性方程如果知道(1)的k个线性无关的解x1,x2,xk,则可通过一系列同类型的 变换,使方程降k阶,且新方程是n-k阶齐线性的.具体做法如下:作变换x=xky,则所以(1)被化为因xk是(1)的解,故的y系数为0.因此,再作变换z=y,并在xk0的区间上将y(n)的系数变为1,得方程方程(1),(2)的解的关系是所以(2)有k-1个线性无关的解事实上,显然zi

4、是(2)的解.其次,假设有关系式 a 1 z 1+a 2z 2+a k-1z k-10,即因为x1,x2,xk线性无关,所以 a 1=a 2=a k-1=a k=0,因此z1,z2,zk-1线性无关.继续上述做法.最后将原方程降k阶.特别地,设 是的解,作变换x=x1(t)y,可得c2=1,c3=0得另一特解 线性无关.事实上,上述结果也可直接运用Liouville公式得到例 6 已知 的特解 求通解.解 令 则代入原方程,得例 6 第二宇宙速度的计算 设法计算发射人造卫星的最小初速(第二宇宙速度)使卫星摆脱地球引力绕太阳运行,成为一颗人造行星.先建立物体垂直上抛运动的微分方程.以M,m表示地

5、球与物体的质量(物体的体积,质相对于地球很小,可看作质点).由牛顿万有引力定律(在Kepler关于行星绕太阳运行的三定律的基础上,Newton通过研究两体问题而导出了万有引力定律,而万有引力定律包含了开普勒的三定律).地球作用于物体的引力Fr为物体重心到地球中心的距离.再根据牛二定律(运动物体所受力总和=质量x加速度)初始条件:t=0,r=R(=6.37x106m),(初速).设 则 方程化为注意到r=R时,v=v0所以要使卫星摆脱地球引力,须足够大,所以要使上式成立,必须因此,最小初速将 K=6.672x10-23(千米)3/克秒2,M=5.977x1027克,R=6.37x106米=6.3

6、7x103千米,代入得第二宇宙速度第二宇宙速度 v0=11.2x103米米/秒秒.如果已知重力加速度g，而K，M未知，则由地球表面r=R的重力加速度g9.81米米/秒秒2可可得得米米/秒秒.若要计算使卫星摆脱太阳引力的初速，则将改为太阳的质量1.983x1033克，地球的半径R改为地球离太阳的距离1.495x108千米，得，要摆脱太阳引力的初速(即相对太阳的速度)超过又地球绕太阳速度是29.76千米秒，所以地球上卫星发射后的速度应超过42.11-29.76=12.35千米秒(摆脱地球引力后的速度)设v0为地球上发射的初速，那么即，第三宇宙速度是第三宇宙速度是16.716.7千米秒千米秒.若计算第一宇宙速度第一宇宙速度(7.9(7.9千米秒千米秒)，可类似方法建立方程组使运行轨道为圆的最小初速即为第一宇宙速度