四(微分方程模型)训练题2
- 格式:doc
- 大小:146.00 KB
- 文档页数:3
微分方程练习题及答案1、已知微分方程组⎪⎪⎩⎪⎪⎨⎧=+-=--05303y x dtdy y x dt dx (1)利用matlab 软件求此方程组在初始条件1|,2|11====t t y x 下的特 解,并画出解函数()y f x =的图形。
(2)利用matlab 软件分别用 ode23、ode45 求微分方程初值问题的数值解(近似解),求解区间为[0,3]t ∈,并作图来比较两种求解器之间的差异。
2、已知微分方程组⎪⎪⎩⎪⎪⎨⎧=-+=++03023y x dtdy y x dt dx (1)利用matlab 软件求此方程组在初始条件2|,1|00====t t y x 下的特 解,并画出解函数()y f x =的图形。
(2)利用matlab 软件分别用 ode23、ode45 求微分方程初值问题的数值解(近似解),求解区间为[0,2]t ∈,并作图来比较两种求解器之间的差异。
1、参考答案:(1)程序代码:syms x y t[x,y]=dsolve('Dx-x-3*y=0','Dy-3*x+5*y=0','x(1)=2','y(1)=1','t') ezplot(x,y,[0,3]);(2)程序代码:M函数文件verderpol.m:function xprime=verderpol(t,x)xprime=[x(1)+3*x(2); 3*x(1)-5*x(2)];在程序中调用此函数:clear;y0=[2;1];[t,x]=ode45('verderpol',[0,3],y0); plot(x(:,1),x(:,2),'r-'); hold onclear;y0=[2;1];[t,x]=ode23('verderpol',[0,3],y0); plot(x(:,1),x(:,2),'b-');2、参考答案:(1)程序代码:syms x y t[x,y]=dsolve('Dx+3*x+2*y=0','Dy+x-3*y=0','x(0)=1','y(0)=2','t') ezplot(x,y,[0,2]);(2)程序代码:M函数文件verderpol.m:function xprime=verderpol(t,x)xprime=[-3*x(1)-2*x(2); 3*x(2)-x(1)];在程序中调用此函数:clear;y0=[1;2];[t,x]=ode45('verderpol',[0,2],y0); plot(x(:,1),x(:,2),'r-'); hold onclear;y0=[1;2];[t,x]=ode23('verderpol',[0,2],y0); plot(x(:,1),x(:,2),'b-');。
第十章 微分方程习题一.填空题:(33)1-1-40、微分方程4233''4''')'(x y x y y =++的阶数是. 1-2-41、微分方程0'2'2=+-xy yy xy 的阶数是. 1-3-42、微分方程0d d d d 22=++s x sx s 的阶数是.1-4-43、x y y y y sin 5''10'''4)()4(=-+-的阶数是 . 1-5-44、微分方程xy x y2d d =满足条件1|'0==x y 的特解是. 1-6-45、微分方程0d d =+y x y的通解是.1-7-46、方程y e y x='的通解是 . 1-8-47、方程y y y ln '=的通解是 . 1-9-48、方程04'4''=+-y y y 的通解是 . 1-10-49、方程04'4''=+-y y y 的通解是 . 1-11-50、方程013'4''=+-y y y 的通解是 .1-12-51、已知特征方程的两个特征根,3,221-==r r 则二阶常系数齐次微分方程为1-13-52、微分方程xe y =''的通解为 .1-14-53、微分方程x ey xsin ''2-=的通解为 .1-15-54、若0d ),(dx ),(=+y y x Q y x P 是全微分方程, 则Q P ,应满足 .1-16-55、与积分方程xy x f y x x d ),(0⎰=等价的微分方程初值问题是 .1-17-56、方程0d )2(d )(22=-++y xy x x y xy 化为齐次方程是 . 1-18-57、通解为21221,(C C e C e C y xx +=为任意常数)的微分方程为 .1-19-58、方程yx e y -=2'满足条件00==x y 的特解是 .1-19-59、方程0dy 1dx 2=-+x xy 化为可分离变量方程是 1-20-60、方程xy y 2'=的通解是1-21-61、方程x yxy x y x y d d d d 22=+化为齐次方程是1-22-62、若t y ωcos =是微分方程09''=+y y 的解, 则=ω .1-23-63、若ktCe Q =满足Qdt dQ03.0-=, 则=k .1-24-64、y y 2'=的解是1-25-65、某城市现有人口50(万), 设人口的增长率与当时的人口数x (万)和x -1000的积成正比,则该城市人口)(t x 所满足的微分方程为1-26-66、圆222r y x =+满足的微分方程是1-27-67、axae y =满足的微分方程是1-28-68、一阶线性微分方程)()(d dyx Q y x P x =+的通解是 .1-29-69、已知特征方程的两个根3,221-==r r , 则二阶常系数线性齐次微分方程为 .1-30-70、方程25x y =是微分方程y xy 2'=的解.1-31-71、二阶常系数非齐次微分方程的结构为其一个特解与之和.1-32-72、二阶常系数齐次线性微分方程0'''=++qy py y 对应的特征方程有两个不等实根,则其通解为.1-33-73、将微分方程0)2()(22=---dy xy x dx y xy 写成齐次微分方程的标准形式为二.选择题:(29)2-1-56、微分方程y x2dx dy=的通解是 ( )A.2x y =B. 25x y =C. 2Cx y = D.Cx y =2-2-57、微分方程0dy 1dx 2=-+x xy 的通解是 ( ) A.21x ey -= B.21x Cey -= C.x C y arcsin = D. 21x C y -=2-3-58、下列方程中是全微分方程的是 ( )A.0dy dx )(2=--x y x B. 0dy dx =-x y C. 0dy )(1dx )1(=-++xy y xy D.0dy dx )(22=++xy y x 2-4-59、下列函数组中,线性无关的是 ( )A.x x e e 32,B.x x 2sin ,2cosC. x x x sin cos ,2sinD.2ln ,ln x x2-5-60、方程03'2''=--y y y 的通解是 ( )A.x x e C e C y 321--+=B. xx e C e C y 321+= C. x x e C e C y 321-+= D. x x e C e C y 321+=-2-6-61、方程0''=+y y 的通解是 ( )A.x C y sin =B.x C y cos =C.x C x y cos sin +=D.x C x C y cos sin 21+=2-7-62、下列方程中是可分离变量的方程是 ( )A. xy y x -=33dx dyB.0dy 2dx )3(2=++xy y e x C. 234dx dy xy y x += D.y x xy y 321dx dy ++= 2-8-63、微分方程0cot '=-x y y 的通解是 ( ) A. x C y cos = B. x C y sin = C. x C y tan = D. x C y csc = 2-9-64、已知微分方程0''=+-p y y 的通解为)(212x C C e y x +=,则p 的值是 ( )A.1B.0C.21D.412-10-65、微分方程02'=-y y 的通解是 ( )A.C x y +=2sinB.C e y x +=24C.x Ce y 2=D. xCe y =2-11-66、方程xy 2dx dy=的通解是 ( )A.C e x +2B.Cxe+2C. 2Cx eD. 2)(C x e +2-12-67、xe y -=''的通解为=y ( ) A.x e -- B. xe - C. 21C x C e x ++- D.21C x C e x ++-- 2-13-68、微分方程x e 21dx dy -=满足10-==x y 的特解为 ( )A.1221+-=-x e y B. 3221-=-x e y C. C ey x +-=-212 D.212121--=-xe y2-14-69、微分方程0ydy -dx 3=x 的通解是 ( )A.C y x =-2422B. C y x =+2422C. 02422=-y xD. 12422=+y x2-15-70、微分方程0ydy -dx 3=x 的通解是 ( )A.222=+y xB. 933=+y xC. 133=+y x D. 13333=+y x2-16-71、过点,0()2-的曲线,使其上每一点的切线斜率都比这点纵坐标大5的曲线方程是( )A.32-=x yB. 52+=x yC.53-=x e yD.5-=x Ce y 2-17-72、齐次方程x yxy tandx dy =化为可分离变量的方程, 应作变换 ( ) A. 2ux y = B. 22x u y = C. ux y = D.33x u y = 2-18-73、设方程)()('x Q y x P y =+有两个不同的解21,y y ,若21y y βα+也是方程的解,则( )A.βα=B. 0=+βαC. 1=+βαD. βα,为任意常数2-19-74、方程dx 2dx dy y x x =+的通解是 ( )A.x Cx y +=2B. x x C y +=2sinC. C x y +=2cosD.C x y +=22-20-75、下面各微分方程中为一阶线性方程的是 ( )A.x y xy =+2'B .x xy y sin '=+C .x yy =' D .xy y -=2' 2-21-76、曲线上任一点P 的切线均与OP 垂直的曲线方程是 ( )A.y x y -=' B. y x y =' C. x y y -=' D. x y y ='2-22-77、方程2)3(,0'==+y y y 的解是 ( )A.x e y -=32B. x e y --=32C. 32-=x e yD. 32--=x e y2-23-78、微分方程x y y ln '=的通解是 ( )A.x x e y ln =B. x x Ce y ln =C. x x x e y -=lnD. xx x Ce y -=ln 2-24-79、下列哪个不是方程y y 4''=的解 ( )A. x e y 22=B. x e y 2=C. x e y 2-=D. x e y 2=2-25-80、方程0sin '''653)4(=-+++y y y y x xy y 的阶是 ( ) A. 6 B. 5 C. 4 D. 32-26-81、如果一条曲线在它任意一点的切线斜率等于y x2-,则这条曲线是( )A. 椭圆B. 抛物线C. 双曲线D. 圆2-27-82、下列可分离变量的方程是 ( )A. xy y x dx dy-=33 B. 02)3(2=++xydy dx y e xC. xy yx dx dy += D.y x xy y dx dy 321++= 2-28-83、微分方程0cot '=-x y y 的通解是 ( ) A. x C y cos = B. x C y sin = C. x C y tan = D. x C y csc =2-29-84、已知微分方程0''=+-p y y 的通解为)(212x C C e y x+=,则p 的值( )A. 1B. 0C. 21D. 41三.计算题:(59)3-1-52、0d tan sec d tan sec 22=+y x y x y x 3-2-53、0ln '=-y y xy3-3-54、0d sec )2(d tan 32=-+y y e x y e x x 3-4-55、y x y y x x y 22222')1(=-+- 3-5-56、y xe y e x dx dy +-=-3-6-57、0)1()1(=-++xdy y ydx x3-7-58、x x y y y x d sin cos d sin cos =,4|0π==x y3-8-59、0)0(,02')1(22==+-y xy y x3-9-60、1)(,ln 2'==e y x y y3-10-61、x x y y y x d sin cos d sin cos =, 4|0π==x y3-11-62、0y)dx -(x dy )(=++y x3-12-63、)ln (ln dx d x y y yx-=3-13-64、0)2(22=+-dy x dx xy y3-14-65、x y x y xy tan'=-3-15-66、x y x y x y xy ++=-ln)('3-16-67、dx dy xy dx dy x y =+223-17-68、x y y x y +=', 2|1==x y3-18-69、x y x y y +=', e y e x ==|3-19-70、2|,'122=-=-=x y y x y xy3-20-71、x x y x y sin 1'=+, 1|==πx y3-21-72、x e x y x y 43'=-3-22-73、342'x xy y =-3-23-74、x y x y ln 11'=-3-24-75、x e y x x y x 21'=-+ 3-25-76、x x y y sec tan '=-, 0|0==x y3-26-77、x x y x y sin 1'=+, 1|==πx y3-27-78、22112'x y x xy +=+-, 0|0==x y3-28-79、x xy xy ln '=-, e y e x ==|3-29-80、22d dyxxe xy x -+=3-30-81、)sin (cos d dy2x x y y x -=+ 3-31-82、5d dyxy y x =- 3-32-83、02d dy4=++xy xy x 3-33-84、4)21(3131d dy y x y x -=+3-34-85、xy xy x 2d dy 2-= 3-35-86、x y y +='''3-36-87、01)'(''2=++y yy 3-37-88、01''3=+y y3-38-89、y y 3''=, 1|0==x y , 2|'0==x y3-39-90、223''yy =, 1|3==x y , 1|'3==x y3-40-91、02''=+y y 3-41-92、013'4''=++y y y3-42-93、0'2''=+-y y y 3-43-94、04'5''=+-y y y3-44-95、04'3''=--y y y , 0|0==x y , 5|'0-==x y 3-45-96、029'4''=++y y y , 0|0==x y , 15|'0==x y 3-46-97、0'4''4=++y y y , 2|0==x y , 0|'0==x y 3-47-98、0'4''4=++y y y , 2|0==x y , 0|'0==x y 3-48-99、013'4''=+-y y y , 0|0==x y , 3|'0==x y 3-49-100、04'4''=+-y y y , 0|0==x y , 1|'0==x y3-50-101、xe y y y 2'''2=-+ 3-51-102、x e y y xcos ''+=+ 3-52-103、xe x y y y 3)1(9'6''+=+-3-53-104、'''22xy y y e --= 3-54-105、123'2''+=--x y y y3-55-106、''sin20y y x ++=, 1|==πx y , 1|==πx y 3-56-107、52'3''=+-y y y , 1|0==x y , 2|'0==x y 3-57-108、xe y y y 29'10''=+-,76|0==x y ,733|'0==x y 3-58-109、xxe y y 4''=-, 0|0==x y , 1|'0==x y 3-59-110、x xe y y y 26'5''=+-四.应用解答题:(14)4-1-9、一曲线通过点)3,2(, 它在两坐标轴间的任一切线段均被切点所平分, 求这曲线方程.4-2-10、已知⎰--=+xx x y t t y t t 03231d )(12, 求函数)(x y 4-3-13、求一曲线, 这曲线通过原点, 并且它在点),(y x 处的切线斜率等于y x =2.4-4-14、试求x y =''的经过点)1;0(M 且在此点与直线12+=xy 相切的积分曲线.4-5-15、设某曲线,它上面的任一点的切线与两坐标轴所围成的三角形面积总等于2,求这条曲线的方程所满足的微分方程.4-6-16、已知某曲线经过点)1,1(, 它的切线在纵轴上的截距等于切点的横坐标,求它的方程.4-7-17、设可导函数)(x ϕ满足⎰+=+xx t t t x x 01d sin )(2cos )(ϕϕ, 求)(x ϕ.4-8-10、已知某商品需求量Q 对价格p 的弹性为22p Ep EQ-=, 最大需求量为1000=Q , 求需求函数)(p f Q =.4-9-11、设质量为m 的物体在高空中静止下落, 空气对物体运动的阻力与速度成正比. 求物体下落的数率v 与时间t 的关系, 再求物体下落距离与时间t 的关系4-10-12、在串联电路中, 设有电阻R, 电感L 和交流电动势t E E ωsin 0=, 在时刻0=t 时接通电路, 求电流i 与时间t 的关系(0E , ω为常数).4-11-13、如图, 位于坐标原点的我舰向位于x 轴上)0,1(A 点处的敌舰发射制导鱼雷, 鱼雷始终对准敌舰, 设敌舰以常数0v 沿平行与y 轴的直线行驰, 又设鱼雷的速度为02v , 求鱼雷的航行曲线方程.4-12-14、根据经验可知, 某产品的纯利润L 与广告支出x 有如下关系)(ddL L A k x -=,(其中0,0>>A k ), 若不做广告, 即0=x 时纯利润为0L , 且A L <<00, 试求纯利润L 与广告费x 之间的函数关系.4-13-15、在宏观经济研究中, 知道某地区的国民收入y , 国民储蓄S 和投资I 均是时间t 的函数, 且在任一时刻t , 储蓄)(t S 为国民收入)(t y 的101, 投资额)(t I 是国民收入增长率t d dy 的31. 设0=t 时国民收入为5(亿元),假定在时刻t 的储蓄全部用于投资,试求国民收入函数.4-14-16、试建立描述市场价格形成的动态过程的数学模型.五.证明题:(2)5-1-18、设),(1x y )(2x y 是二阶齐次线性方程0)(')(''=++y x q y x p y 的两个解,令)()(')(')()(')(')()()(21212121x y x y x y x y x y x y x y x y x w -==证明: )(x w 满足方程0)('=+w x p w5-2-19、设1y , 2y , 3y 是线性方程)()(d dyx Q y x P x =+的3个相异特解,证明1213y y y y --为一常数.部分应用题答案487.在串联电路中, 设有电阻R, 电感L 和交流电动势t E E ωsin 0=, 在时刻0=t 时接通电路, 求电流i 与时间t 的关系(0E , ω为常数).解. 设)(t i i =, 由回路电压定律t E dt di LRi ωsin 0=+, 即t L E L R dt di ωsin 0=+∴⎰+⎰⎰=-]sin [)(0C dt te L E e t i t dt LR L Rω=⎰+-]sin [0C dt te L E ett L R LR ω=)cos sin (2220t L t R L R E Cet LR ωωωω-++-将0|0==t i 代入通解得2220L R LE C ωω+=∴)cos sin ()(2220t L t R Le L R E t i t LR ωωωωω-++=-488.设质量为m 的物体在高空中静止下落, 空气对物体运动的阻力与速度成正比. 求物体下落的数率v 与时间t 的关系, 再求物体下落距离与时间t 的关系 解:.物体重力为mg w =, 阻力为kv R -=, 其中g 是重力加速度, k 是比例系数.由牛顿第二定律得kv mg dt dv m-=,从而得线性方程gv m k dt dv =+, 0|0==t v ∴⎰--+=+⎰⎰=t mkdt dt Ce g k mC dt ge e v kmmk][, 将0|0==t v 代入通解得g k m C -= ∴)1(tm k e g k m v --=, 再积分得122C ge k m gt k m S t m k++=-,将0|0==t S 代入求得gk m C 221-=∴)1(22-+=-t m ke g k m gt k m S489.如图, 位于坐标原点的我舰向位于x 轴上)0,1(A 点处的敌舰发射制导鱼雷, 鱼雷始终对准敌舰, 设敌舰以常数0v 沿平行与y 轴的直线行驰, 又设鱼雷的速度为02v , 求鱼雷的航行曲线方程.解:设鱼雷的航行曲线方程为)(x y y =, 在时刻t , 鱼雷的坐标巍巍),(y x P , 敌舰的坐标为),1(0t v Q .因鱼雷始终对准敌舰, 故x y t v y --=1'0, 又弧OP 的长度为⎰=-xtv dx y 0022'1,从以上两式消去t v 0得''121''')1(2y y y y x -+=--, 即2'121'')1(y y x +=-根据题意, 初始条件为0)0(=y , 0)0('=y令p y =', 原方程化为2121')1(p p x +=-, 它是可分离变量得方程, 解得21)1(112--=++x C p p , 即21)1('1'12--=++x C y y将0)0('=y 代入上式得11=C , 故21)1('1'2--=++x y y而21)1(''1'1'122--=-+=++x y y y y , 得2121)1()1(21'x x y -+-=-积分得22321)1(31)1(C x x y +-+--=, 将0)0(=y 代入上式得322=C , 所以鱼雷的航行曲线为32)1(31)1(2321+-+--=x x y 490.根据经验可知, 某产品的纯利润L 与广告支出x 有如下关系)(ddL L A k x -=,(其中0,0>>A k ), 若不做广告, 即0=x 时纯利润为0L , 且A L <<00, 试求纯利润L 与广告费x 之间的函数关系.解:依题意得)(L A k dx dL-=,00|L L x ==, 解可分离变量得微分方程, 得通解 kx Ce A L -+=, 将00|L L x ==代入通解, 得A L C -=0, 所以纯利润L 与广告费x 之间的函数关系为kxe A L A x L --+=)()(.491.在宏观经济研究中, 知道某地区的国民收入y , 国民储蓄S 和投资I 均是时间t 的函数, 且在任一时刻t , 储蓄)(t S 为国民收入)(t y 的101, 投资额)(t I 是国民收入增长率t d dy 的31.设0=t 时国民收入为5(亿元), 假定在时刻t 的储蓄全部用于投资,试求国民收入函数.解:依题意:y S 101=, dt dy I ⋅=31, 解之得通解t Ce y 103=, 将5|0==t y 代入通解得5=C , 所以国民收入函数为t e y 1035=492.试建立描述市场价格形成的动态过程的数学模型.解:设在某一时刻t , 商品的价格为)(t p , 因供需差价, 促使价格变动. 对新的价格,又有新的供需差, 如此不断地调节价格, 就构成了市场价格形成的动态过程.假设价格)(t p 的变化率dt dp与需求和供给之差成正比. 记需求函数为),(r p f , 供给函数为)(p g , 其中r 为参数. 于是得微分方程)](),([p g r p f k dt dp-=,0)0(p p =, 其中0p 为0=t 时商品的价格, k 为正常数.若需求供给函数均为线性函数, b kp r p f +-=),(, d cp p g +=)(, 则方程为)()(d b k p c k k dt dp-++=,0)0(p p =, 其中d c b k ,,,均为正常数, 其解为 c k db ec kd b p t p t c k k +-++--=+-)(0)()(下面对所得结果进行讨论:(1) 设p 为静态均衡价格, 则应满足0)(),(=-p g r p f , 即d p c b p k +=+-,则c k db p +-=, 从而价格函数p e p p t p c k k +-=+-)(0)()(,取极限: p t p t =∞→)(lim .它表明: 市场价格逐步趋于均衡价格. 若初始价格p p =0 , 则动态价格就维持在均衡价格p 上, 整个动态过程就变为静态过程.(2) 由于t c k k e c k k p p dt dp )(0)()(+-+-=, 所以当p p >0时, 0<dt dp , )(t p 单调下降向p靠拢, 这说明: 初始价格高于均衡价格时,动态价格会逐渐降低, 逐渐接近均衡价格; 而当初始价格低于均衡价格时, 动态价格会逐渐增高, 逐渐接近均衡价格.。
广州大学2010年数学建模暑期培训答题卡
题目:
以前美国原子能委员会将放射性核废料装在封闭的圆桶里,扔到水深为91米的海底。
生态学家和科学家担心这种做法不安全,而原子能委员会向他们保证,圆桶绝不会破漏。
经过周密的试验,证明了圆桶的封闭性确实很好,但工程师们提出疑问,圆桶是否会因与海底相撞而发生破裂而原子能委员会却说决不会。
于是,工程师们通过大量试验发现:当圆桶在海水中的速度超过12.2米时,圆桶就会因碰撞破裂。
现在的问题是,圆桶到达91米的海底时的速度是不是超过了12.2米的问题。
现在来计算这一问题。
其中圆桶的质量(重力)为)(456.239N W =,海水的浮力为)/(94.10253
m N ,圆桶的体积为3
208.0m V =,
圆桶下沉过程所受
的阻力为D ,工程师们做了大量试验得出结论:阻力与圆桶的方位无关,而与下沉的速度成正比,比例系数12.0=k . 要求给出计算机仿真,拟题人:钟老师
答案: 海中物体下落问题 问题摘要 问题重述
对于以前美国原子能委员会将放射性核废料装在封闭的圆桶里,扔到水深为91米的海底很多科学家对其安全问题持有怀疑,之前工程师已对圆桶经过周密的试验,证明了圆桶的封闭性确实很好,但却发现当圆桶在海水中的速度超过12.2米时,圆桶就会因碰撞破裂, 而我们要解决的问题是圆桶在到达91米的海底时的速度会不会超过了12.2米的问题。
问题假设 1. 假设圆桶在下落到海底之前不受任何因素干扰; 2. 假设圆桶内的东西不会在下落时发生爆炸;
3. 假设圆桶在下落过程中受到的重力加速度不会因为深度的改变而改变;
4. 假设圆桶在下落时是垂直下落的;
5. 假设圆桶在下落开始时速度为零,海面距离为零; 符号说明 1.W 圆桶受到的重力;
2. F 圆桶受到海水的浮力;
3.M 圆桶的质量;
4. s 圆桶下落过程中与海面的距离;
5. t 圆桶在下落过程中于海平面距离s 时的时间力;
6. t v 圆桶在t 时刻下落时的速度;
7. k 阻力与速度的比例系数;
8. t a 圆桶在t 时刻下落过程中的加速度; 9. V 圆桶的体积;
10. t D 圆桶在t 时刻下落时受到的阻; 11. g 为重力加速度;
问题分析
圆桶在下沉过程中受到3个力的作用的作用:重力W 垂直向下,受到海水的浮力F=V )/(94.10253m N 垂直向上,受到的阻力t D 且与圆桶在水中的方位无关,与下落速度
t v 成正比,比例系数为k ;速度t v 积分与t 时间比就是t 时刻的加速度t a ,t 时刻的距离s 积
分与t 时间积分的比就是速度t v ,根据牛顿第二定理可得出距离s 与t v 的二介的关系。
模型的建立
根据上述可知:t D =t kv ; (1)
下落过程时垂直下落, 在t 时刻是下落到与海平面距离s 处,可得速度与距离的函数为
()v s ,g 2
9.8()m s =;根据牛顿定律可得:
M=
W g
; (2);
t v =
d s d t
; (3); t a =
d v d t
; (4 );
t t M a W F D =--; (5);
综上(1)(2)(3)(4)(5)的到微分方程为(6)(7): 2
2d s ds M W F k
dt
dt
=--; (6);
(0)0,(0)0s v ==; (7); 模型的求解
根据(2)(3)(4)(6)(7)进行化为一介微分方程:
v dv M ds W F kv
=--; (8);
两边积分:
v dv M ds W F kv
=
--⎰
⎰
; (9)
; 求的解为:
2
1()()
ln(
)v s W F W F kv s s M
k k W F --+=--;
即2
()()
ln(
)g v s W F W F kv s s W
k
k
W F
--+=
-
- ; (10)
S=91,即(91)v 通过MA TLAB 求的(91)13.64v ≈m s ; 根据(10)的函数关系式大概描述如下图: m
v s
s/m
结果分析:图像和结果所显示,圆桶在到达海底91米时,圆桶的速度已达到了和超过12.2m s ,因此他们的做法是不安全的。
模型分析
对于动态变量可通过建立微分方程,建立模型,可以求出答案,但没有用到计算。