02 第二节 一阶微分方程
- 格式:doc
- 大小:348.50 KB
- 文档页数:8
一阶微分方程第二节 一阶微分方程一阶微分方程的一般形式为F (x ,y ,y ′)=0或y ′=f (x ,y ),其中F (x ,y ,y ′)是x ,y ,y ′的已知函数,f (x ,y )是x ,y 的已知函数.这一节只介绍几种较简单的一阶微分方程的解法.它们通过求积分就可以找到未知函数与自变量的函数关系,我们称这种求解微分方程的方法为初等积分法.一、 可分离变量的方程 形如xyd d =f (x )g (y ) (10-2-1)或M 1(x )M 2(y )d y =N 1(x )N 2(y )d x (10-2-2) 的一阶微分方程称为可分离变量方程.其中f (x ),g (y )及M 1(x ),M 2(y ),N 1(x )及N 2(y )均为已知连续函数.方程(10-2-1)的求解步骤如下: 先将方程(10-2-1)分离变量得将21y -作为分母时丢失了两个特解.故所求方程的通解为:arcsin y =x +C (C 为任意常数), 另外还有两个特解y =±1.例2 已知某商品的需求量x 对价格P 的弹性e =-3P 3,而市场对该商品的最大需求量为1(万件),求需求函数.解 需求量x 对价格P 的弹性e =pxx P d d . 依题意,得pxx P d d =-3P 3,于是xx d =-3P 2d P ,积分得ln x =-P 3+C 1,即x =C3P -e (C =1C -e ).由题设知P =0时,x =1,从而C =1.因此所求的需求函数为x =3P -e .例3 根据经验知道,某产品的净利润y 与广告支出x 之间有如下关系:xy d d =k (N -y ),其中k ,N 都是大于零的常数,且广告支出为零时,净利润为y 0,0<y 0<N ,求净利润函数y =y (x ),解 分离变量yN y -d =k d x ,两边同时积分得-ln |N -y |=kx +C 1 (C 1为任意常数), 因N -y >0,所以ln |N -y |=ln(N -y ),上式经整理得y =N -C e -kx (C =1C -e >0).将x =0,y =y 0代入上式得C =N -y 0,于是所求的利润函数为y =N -(N -y 0)e -kx .由题设可知xy d d >0,这表明y (x )是x 的单调递增函数;另一方面又有)(lim x y x ∞→=N ,即随着广告支出增加,净利润相应地增加,并逐渐趋向于y =N .因此,参数N 的经济意义是净利润的最大值.二、 齐次微分方程 1. 齐次微分方程 形如xy d d =⎪⎭⎫ ⎝⎛xy f (10-2-3) 的一阶微分方程,称为齐次微分方程,简称齐次方程.对于方程(10-2-3),通常可通过变量替换u =x y将方程化为可分离变量的方程来解.具体过程如下:令 u =x y(或y =ux ),其中u 是新的未知函数.对y =ux 两端关于x 求导,得xyd d =u +x xud d . 代入(10-2-3)得u +x xu d d =f (u ). 分离变量并积分得⎰-uu f u )(d =⎰x x d ,即Φ(u )=ln|x |+C (C 为任意常数),其中Φ(u )是⎰-u u f u )(d 的一个原函数,再将u =xy代入上式中,便得到方程(10-2-3)的通解Φ(x y)=ln|x|+C .上面的推导要求f (u )-u ≠0,如果f (u )-u =0,也就是⎪⎭⎫ ⎝⎛x y f =xy.这时,方程(10-2-3)为 x yd d =x y .这已是一个可分离变量的方程,不必作代换就可求出它的通解为y =Cx .例4 求微分方程xy xy d d =x 2+y 2满足条件y |x =e=2e 的解.解 原方程可化为x y d d = y x +xy,这是一个齐次方程.作代换u =x y,即y =ux ,则xyd d =u +x xud d . 代入前一方程得u +x x u d d =u 1+u 即 x x u d d =u1, 分离变量并积分得u 2=2ln |x |+2C (C 为任意常数),将u 替换为x y,便得原方程的通解:y 2=2x 2ln |x |+2Cx 2,再将初始条件代入通解得4e 2=2e 2·ln e +2C e 2,求得 C =1, 于是,所求的特解为y 2=2x 2(ln |x |+1).例 5 设甲、乙两种商品的价格分别为P 1,P 2,且价格P 1相对于P 2的弹性为21d d P P P P 12=1212PP P P +-,求价格P 1与P 2的函数关系.解 将所给方程整理为21d d P P =21212111P P P P P P +-.这是齐次方程.令u =21P P ,即P 1=uP 2,则21d d PP=u +P 22d d Pu ,代入上式得 u +P 22d d P u =uu+-11·u . 整理得⎪⎭⎫ ⎝⎛--211u u d u =222d P P.两边积分得u1-ln |u |=2ln |P 2|+C 1 (C 1为任意常数). 将u 替换为21P P ,便得方程的通解(注意到u >0,P 22>0)12P P e=CP 1P 2(C =1C e , C 为正数).2. 可化为齐次方程的微分方程形如xy d d =⎪⎪⎭⎫ ⎝⎛++++222111C y b x a C y b x a f (10-2-4)的微分方程,当C 1=C 2=0时,就是一个齐次方程.当C 1,C 2中至少有一个不为零时,尽管本身不是齐次方程,但经过适当的变量替换后,可化为齐次方程.下面分两种情况讨论:(1) 若a 1b 2-a 2b 1≠0,这时方程组⎩⎨⎧=++=++0,0222111C y b x a C y b x a有惟一解x =α,y =β.作变量替换⎩⎨⎧-=-=,,βαy v x u则222111C y b x a C y b x a ++++=222111)()()()(C v b u a C v b u a ++++++++βαβα=vb u a v b u a 22111++. 于是方程(10-2-4)化为u v d d =⎪⎪⎭⎫ ⎝⎛++v b u a v b u a f 22111.这是关于变量u 和v 的齐次方程.求出其通解后再换回原来的变量x 和y ,即得原方程的通解.(2) 若a 1b 2-a 2b 1=0,这时令21a a =21b b=λ,即有a 1=λa 2,b 1=λb 2. 方程(10-2-4)可写为xyd d =⎪⎪⎭⎫ ⎝⎛++++222122)(C y b x a C y b x a f λ.作变量替换t =a 2x +b 2y ,此时x t d d =a 2+b 2xyd d ,方程(10-2-4)化为xtd d =a 2+b 212()t C f t C λ++.这是关于变量t 和x 的可分离变量的方程.例6 求方程xyd d =51+++-x y x y 的解. 解 解方程组⎩⎨⎧=++=+-05,01x y x y得x =-2,y =-3.作变换x =u -2,y =v -3,原方程化为uv d d =uv u v +-. 这是一个齐次方程,按齐次方程的解法可求得它ln(u 2+v 2)+2arctan u v=C .再将u =x +2,v =y +3代入上式,便得原方程的通解为.ln [(x +2)2+(y +3)2]+2arctan 23++x y =C . 三、 一阶线性微分方程 形如y ′+P (x )y =Q (x ) (10-2-5)的方程叫做一阶线性微分方程.其中P (x ),Q (x )为x 的已知连续函数,Q (x )称为自由项.如果Q (x )≡0,方程(10-2-5)即为y ′+P (x )y =0. (10-2-6)该方程称为一阶齐次线性微分方程.而当Q (x ) ≠0时,方程(10-2-5)称为一阶非齐次线性微分方程.也称(10-2-6)为(10-2-5)所对应的齐次方程.注意这里所说的齐次方程与上段讨论的齐次方程是不同的.下面来讨论一阶非齐次线性方程(10-2-5)先考虑非齐次线性方程(10-2-5)所对应的齐次方程(10-2-6)的通解.显然y=0是它的一个解,当y≠0时分离变量得y d=-P(x)d x.y两边积分得ln|y|=⎰-xx(+C1,P d)即y=C⎰-x x P de)((C=±1C e).y=0也是方程(10-2-6)的解,这时在上式中取C=0即可.于是得到方程(10-2-6)的通解为y=C⎰-x x P de)((C为任意常数).(10-2-7)再利用“常数变易法”求非齐次线性方程(10-2-5)的通解.由于方程(10-2-5)与(10-2-6)的左端相同,右端不同,方程(10-2-5)的左端比方程(10-2-6)的左端多了一项Q(x),因此,我们猜想方程(10-2-5)的通解也具有(10-2-7)的形式,而其中的C不可能还是常数,而是x的某个函数C(x).于是,可设方程(10-2-5)的解为y=C(x)·⎰-x x P de)(,(10-2-8)其中C(x)是待定函数.将(10-2-8)代入方程(10-2-5),得[C (x ) ⎰-xx P d e )(]'+P (x )C (x )⎰-xx P d e )(=Q (x ).化简,得C '(x )=Q (x )⎰xx P d e )(.上式两端同时积分,得C (x )=⎰)(x Q⎰xx P d e )(d x +C (C 为任意常数).将上式代入(10-2-8)式,得非齐次线性方程(10-2-5)的通解 y =⎰-xx P d e)([⎰)(x Q⎰xx P d e )(d x +C ] (C 为任意常数). (10-2-9)这种将任意常数变成待定函数求解的方法,称为常数变易法.将通解(10-2-9)改写为y =C ⎰-xx P d e)(+⎰-xx P d e)(⎰⎰xx Q xx P d )e(d )(.不难看出: 通解由两部分构成,其中第一项是方程(10-2-5)所对应的齐次线性方程(10-2-6)的通解,第二项是方程(10-2-5)本身的一个特解[对应于通解(10-2-9)中C =0的特解].这并不偶然,这是线性方程解的结构的一个重要性质.例7 求方程xy ′+y =e x (x >0)的通解. 解 所给方程可化为y ′+xy=xx e . (10-2-10)先求得方程(10-2-10)对应的齐次线性方程的通解为y =xC , 再利用常数变易法,设方程(10-2-10)的解为y =x x C )(,代入方程(10-2-10)得22)()()(xx C x x C x C x +-'=xxe ,化简,得C '(x )=e x ,积分得C (x )=e x +C ,故得方程(10-2-10)的通解为y =x1 (e x+C )(C 为任意常数). 这也就是所求方程的通解.以上是按“常数变易法”的思路求解,本题也可直接利用通解公式(10-2-9)求解.但是,必须先将方程化为形如方程(10-2-5)的标准形式.这里,P (x )=x1,Q (x )=xx e ,代入公式(10-2-9),得方程的通解为y =⎥⎦⎤⎢⎣⎡+⎰⎰⎰-C x x x x x x xd e e ed d 11=x1(e x+C ). 例8 求方程y ′=3y x y +满足初始条件y (0)=1的特解.解 先求出所给方程的通解.这个方程乍一看不像一阶线性方程,但把它改写成yxd d -y1x =y 2, 则是以y 为自变量,x 为未知函数的一阶线性微分方程.利用通解公式(10-2-9)得 x =⎥⎥⎦⎤⎢⎢⎣⎡+⎰⎰-⎰C y y y yy yd eed 2d 11=[]⎰+-Cy y yyd ee 2ln ln =[]⎰+C y y y d =Cy +21y 3, 将初始条件y (0)=1代入上述通解中,得C =21-,故所求 方程的特解为x =21-y +21y 3. 例9 已知连续函数f (x )满足条件f (x )=t f x t d ⎰303)(+e 2x,求f (x ). 解 因原方程右端函数可导,所以f (x )可导.对方程两端同时求导,得f ′(x )=3f (x )+2e 2x .由一阶线性方程的通解公式,得 f (x )=()⎰+⎰-⎰Cx xx xd e ee d d 3232=e 3x (-2e -x +C )=-2e 2x +C e 3x .例10 设y =f (x )是第一象限内连接点A (0,1),B (1,0)的一段连续曲线,M (x ,y )为该曲线上任意一点,点C 为M 在x 轴上的投影,O 为坐标原点.若梯形OCMA 的面积与曲边三角形CBM 的面积之和为63x +31,求f (x )的表达式.图10-2解 参看图10-2,由题设得2x[1+f (x )]+⎰1)(xtt f d =63x +31, 求导,得21[1+f (x )]+21xf ′(x )-f (x )=22x ,即f ′(x )-x1f (x )=xx 12- (x ≠0).利用一阶线性微分方程的通解公式,得f (x )=⎪⎪⎭⎫ ⎝⎛+-⋅⎰⎰-⎰C x x x x x x xd e e d d 1211=e=x221d x x C x ⎛⎫-+ ⎪⎝⎭⎰=x 2+1+Cx .当x =0时,f (0)=1.说明上述解在x =0时有意义.将条件f (1)=0代入到通解中,得C =-2,于是有f (x )=x 2-2x +1.形如xyd d +P (x )y =Q (x )y a (α≠0,1)(10-2-11)的方程称为伯努利(Bernoulli)方程.它不是线性方程,但是经过适当的变量替换,可将它化成线性方程求解.事实上,只要将方程(10-2-11)两端除以y α,得y-αxy d d +P (x )y 1-α=Q (x ),即xy d d -1αα-11+P (x )y 1-α=Q (x ).若令y 1-α=z , 则上面这个方程为xz d d α-11+P (x )z =Q (x ). (10-2-12)这是一个线性方程.求出这个方程的通解后,用y 1-α替换z ,便得到伯努利方程的通解.例11 求方程y ′+y x x21- =21xy 的通解.解 这是α=21的伯努利方程.方程两边同时除以21y ,得21211y xxx y y -+d d =x .令z =y1-α=21211yy=-,则上面的方程化为xzd d +z x x )1(22-=2x. 这是一阶线性微分方程,其通解为z =⎪⎪⎭⎫ ⎝⎛+⋅⎰⎰-⎰-C x x x x x x x xd e ed d 221211212=⎥⎦⎤⎢⎣⎡+---C x x 43242)1(311=)1(311242x x C---.将21y 替换z ,得原方程的通解为y =2242)1(311⎥⎦⎤⎢⎣⎡---x x C (C 为任意常数).习题10-21. 求下列微分方程的通解或在给定的初始条件下的特解:(1) y ′=xy-+11; (2) xy d x +21x -d y =0;(3) (xy 2+x )d x +(y -x 2y )d y =0; (4) sin x cos 2y d x +cos 2x d y =0;(5)1,0110==+-+=x y y xyx y x d d ;(6) yy ′+x e y =0, y (1)=0; (7) y ′=e 2x -y , 00==x y .2. 物体冷却速度与该物质和周围介质的温差成正比,具有温度为T 0的物体放在保持常温为α的室内,求温度T 与时间t 的关系:3. 求下列微分方程的通解或在给定条件下的特解: (1) xy ′-y -22y x +=0;(2) y ′=x y +sin x y ;(3) 3xy 2d y =(2y 3-x 3)d x ; (4) x 2y ′+xy =y 2, y (1)=1; (5) xy ′=y (ln y -ln x ), y (1)=1; (6) (y -x +2)d x =(x +y +4)d y ; (7) (x +y )d x +(3x +3y -4)d y =0.4. 求下列微分方程的通解或在给定初始条件下的特解: (1) y ′-y =sin x ;(2) y ′-xn y =x n e x; (3) (x -2y )d y +d x =0; (4) (1+x sin y )y ′-cos y =0;(5) y ′-1+x y =(x +1)e x , y (0)=1; (6) y ′+2221212x x y x x +=+,y (0)=23; (7) y ′-y x 1=-x2ln x , y (1)=1; (8) y ′+2xy =(x sin x )·2x -e ,y (0)=1;(9) y ′=234xy y x +;(10) y ′=xy y x +331.5. 设函数f (x )在[1,+∞]上连续,若由曲线y =f (x ),直线x =1,x =t (t >1)与x 轴所围成的 平面图形绕x 轴旋转一周所成的旋转体的体积为V (t )=3π[t 2f (t )-f (1)]. 试求y =f (x )所满足的微分方程,并求该微分方程满足条件y (2)=92的特解.6. 设某生物群体的出生率为常数a ,由于拥挤及对食物的竞争的加剧等原因,死亡率与当时群体中的个体量成正比(比例系数为b >0).如果t =0时生物个体总数为x 0,求时刻t 时的生物个体的总数(注: 将生物群体中的个体量当做时间t 的连续可微变量看待).7. 已知f (x )=x t f xd ⎰⎪⎭⎫ ⎝⎛303+3x -3, 求f (x ).8. 已知某商品的成本C =C (x )随产量x 的增加而增加,其增长率为C ′(x )=xC x +++11, 且产量为零时,固定成本C (0)=C 0>0.求商品的生产成本函数C (x ).9. 某公司对某种电器设备的使用费用进行考察,结果发现,随该电路使用时间x 的延长,它的保养维修费会加倍增长,因而平均单位时间的使用费S 也在增加,即S 为x 的函数S =S (x ),其变化率为a xb S x b x S 21+-=d d ,其中a ,b 均为正常数.若当x =x 0时S =S 0,试问:使用时间为多少时,其平均单位时间的使用费S最高?。
一阶微分方程解法在数学的领域中,一阶微分方程是一个重要的研究对象,它在物理学、工程学、经济学等众多学科中都有着广泛的应用。
那么,什么是一阶微分方程呢?简单来说,一阶微分方程就是指方程中只含有一阶导数的微分方程。
一阶微分方程的一般形式可以表示为:$y' + P(x)y = Q(x)$,其中$y'$表示$y$对$x$的一阶导数,$P(x)$和$Q(x)$是关于$x$的已知函数。
接下来,我们就来探讨一下一阶微分方程的常见解法。
一、可分离变量的一阶微分方程如果一阶微分方程可以写成$g(y)y' = f(x)$的形式,那么我们就称它为可分离变量的一阶微分方程。
对于这种类型的方程,我们可以通过将变量分离,然后两边积分来求解。
具体的求解步骤如下:首先,将方程变形为$\frac{g(y)}{y'}= f(x)$。
然后,将两边分别积分:$\int \frac{g(y)}{y'}dx =\intf(x)dx$。
最后,经过积分运算,求出$y$的表达式。
例如,对于方程$y' = 2xy$,我们可以将其变形为$\frac{dy}{y} = 2xdx$,然后两边积分得到$\ln|y| = x^2 + C$,进而得到$y = Ce^{x^2}$(其中$C$为常数)。
二、一阶线性微分方程一阶线性微分方程是形如$y' + P(x)y = Q(x)$的方程。
对于这种类型的方程,我们可以使用积分因子法来求解。
首先,求出积分因子$\mu(x) = e^{\int P(x)dx}$。
然后,将原方程两边同时乘以积分因子$\mu(x)$,得到:$e^{\int P(x)dx}y' + P(x)e^{\int P(x)dx}y = Q(x)e^{\intP(x)dx}$此时,左边可以变形为$(ye^{\int P(x)dx})'$。
于是,原方程就变成了$(ye^{\int P(x)dx})'= Q(x)e^{\int P(x)dx}$。
第二节 一阶微分方程分布图示★ 可分离变量微分方程★ 例1 ★ 例2★ 例3★ 例4★ 例5 ★ 例6 ★ 例7★ 一阶线性微分方程及其解法★ 例7 ★ 例8★ 例9★ 例10★ 内容小结★ 课堂练习★ 习题6—2内容要点一、可分离变量的微分方程设有一阶微分方程),(y x F dxdy =,如果其右端函数能分解成)()(),(x g x f y x F =,即有)()(y g x f dxdy =. (2.1)则称方程(2.1)为可分离变量的微分方程,其中)(),(x g x f 都是连续函数. 根据这种方程的特点,我们可通过积分来求解. 求解可分离变量的方程的方法称为分离变量法.二、一阶线性微分方程形如)()(x Q y x P dxdy =+ (3.1)的方程称为一阶线性微分方程. 其中函数)(x P 、)(x Q 是某一区间I 上的连续函数. 当,0)(≡x Q 方程(3.1)成为0)(=+y x P dxdy (3.2)这个方程称为一阶齐次线性方程. 相应地,方程(3.1)称为一阶非齐次线性方程.方程(3.2)的通解.)(⎰-=dxx P Cey (3.3)其中C 为任意常数.求解一阶非齐次线性微分方程的常数变易法:即在求出对应齐次方程的通解(3.3)后,将通解中的常数C 变易为待定函数)(x u ,并设一阶非齐次方程通解为,)()(⎰-=dx x P e x u y一阶非齐次线性方程(3.1)的通解为[]⎰-⎰+=⎰dxx P dxx P eC dx ex Q y )()()( (3.5)例题选讲一阶线性微分方程例1(E01)求微分方程xy dxdy 2=的通解.解 分离变量得xdxydy 2=两端积分得⎰⎰=xdxy dy 212||ln C x y +=从而2112xC C xee e y ⋅±=±=+,记,1C e C ±=则得到题设方程的通解 .2x Ce y =例2(E02)求微分方程ydy dx y xydy dx +=+2的通解. 解 先合并dx 及dy 的各项, 得 dx y dy x y )1()1(2-=- 设 ,01,012≠-≠-x y 分离变量得 dxx dy y y 1112-=-两端积分 ⎰⎰-=-dxx dy yy1112得||ln |1|ln |1|ln 2112C x y +-=-于是 2212)1(1-±=-x C y记,21C C ±= 则得到题设方程的通解 .)1(122-=-x C y注:在用分离变量法解可分离变量的微分方程的过程中, 我们在假定0)(≠y g 的前提下, 用它除方程两边, 这样得到的通解, 不包含使0)(=y g 的特解. 但是, 有时如果我们扩大任意常数C 的取值范围, 则其失去的解仍包含在通解中. 如在例2中,我们得到的通解中应该0≠C ,但这样方程就失去特解1±=y ,而如果允许0=C ,则1±=y 仍包含在通解22)1(1-=-x C y 中.例3设一物体的温度为100℃,将其放置在空气温度为20℃的环境中冷却. 试求物体温度随时间t 的变化规律.解 设物体的温度T 与时间t 的函数关系为),(t T T =在上节的例1中我们已经建立了该问题的数学模型:⎪⎩⎪⎨⎧=--==100|)20(0t T T k dtdT )2()1( 其中)0(>k k 为比例常数. 下面来求上述初值问题的解. 分离变量, 得 ;20kdt T dT -=-两边积分,201⎰⎰-=-kdt dT T得 1|20|ln C kt T +-=-(其中1C 为任意常数), 即 kt kt C C kt Ce e e e T --+-=±=±=-1120(其中1C e C ±=).从而 ,20kt Ce T -+=再将条件(2)代入,得,8020100=-=C 于是,所求规律为 .8020kt e T -+=注:物体冷却的数学模型在多个领域有广泛的应用. 例如,警方破案时,法医要根据尸体当时的温度推断这个人的死亡时间,就可以利用这个模型来计算解决,等等.例4设降落伞从跳伞塔下落后, 所受空气阻力与速度成正比, 并设降落伞离开跳伞塔时)0(=t 速度为零, 求降落伞下落速度与时间的关系.解 设降落伞下落速度为),(t v 降落伞下落时,同时收到重力P 与阻力R 的作用. 降落伞所受外力为 kv mg F -=根据牛顿第二定律: αm F =,得到)(t v 满足微分方程kv mg dtdv m -= (1)初始条件 .00==t v将方程(1)分离变量得mdt kvmg dv =-两边积分得⎰⎰=-m dtkvmgdv1)ln(1C mt kv mg k+=--,即 ⎪⎭⎫ ⎝⎛+-=-1C m t k e kv mg 或 ⎪⎪⎭⎫⎝⎛+=--k e C Cekmg v kCtmk 1=-代入初始条件得 kmg C -=故所求特解为 ⎪⎪⎭⎫ ⎝⎛-=-t m k e k mg v 1. 例5(E03)在一次谋杀发生后,尸体的温度从原来的C 37按照牛顿冷却定律开始下降.假设两个小时后尸体温度变为C 35,并且假定周围空气的温度保持C 20不变,试求出尸体温度T 随时间t 的变化规律.又如果尸体被发现时的温度是C 30,时间是下午4点整,那么谋杀是何时发生的?解 根据物体冷却的数学模型,有⎪⎩⎪⎨⎧=>--=.37)0(,0),20(T k T k dt dT其中0>k 是常数.分离变量并求解得ktCeT -=-20,为求出k 值,根据两个小时后尸体温度为C 35这一条件,有2172035⋅-+=k e,求得063.0≈k ,于是温度函数为teT 063.01720-+=,将30=T 代入上式求解t ,有te063.01710-=,即得4.8≈t (小时).于是,可以判定谋杀发生在下午4点尸体被发现前的8.4小时,即8小时24分钟,所以谋杀是在上午7点36分发生的.例6(E04)饮酒量与事故风险率。
大量的研究所提供的数据表明,汽车司机发生事故的风险率R (百分比)与其血液中的酒精浓度b (百分比)有关。
使用两个有代表性的点(0,1%)和(0.14,20%)可用一个指数函数来近似这组数据。
假设风险率R 的变化率与血液酒精浓度b 的关系为kR dbdR =(1)设%10=R ,求满足方程的函数; (2)利用数据点()%2014.0=R ,求k ; (3)用求出的k 写出()b R ;(4)当血液酒精浓度是多少时发生事故的风险率为100%?四舍五入后精确到百分之一。
解 (1)因为kR dbdR =,将方程分离变量并积分得kbCeR c kb R kdbRdR =+==⎰⎰ln 由01.00=R 可得01.0=C()kbeb R 01.0=(2)用第二个点(0.14,0.2)计算数k 。
解方程()kbeb R 01.0=,即()14.020k e =取自然对数得ke14.0ln 20ln =k 14.020ln = 4.2114.0995732.214.020ln ≈==k(3)()b e b R 4.2101.0=(4)把()100=b R %代入3)中得be4.21100=100ln 4.21=b22.04.21605170.44.21100ln ≈==b按照这个模型,当血液酒精浓度达到0.22%时,事故的风险率是100%。
将xy z =代回,得所求通解为 .4)(2s i n 2C x xy xy +=-例7(E05)求微分方程xey dxdy -=+的通解.解 注意到.)(,1)(xex Q x P -==由一阶线性微分方程通解公式得:,e e dx -dx -⎪⎭⎫ ⎝⎛+⎰⋅⎰=⎰-C dx e y x故所求通解为.)(x e C x y -+= 例8(E06) 求方程x xy xy sin 1=+'的通解.解 ,s i n )(,1)(xx x Q xx P ==于是所求通解为⎪⎪⎭⎫ ⎝⎛+⎰⋅⎰=--⎰C dx e x x e y dx x dx x 11sin⎪⎭⎫⎝⎛+⋅=--⎰C dx exxe xxxx ln ln sin).cos (1C x x+-=例9 求下列微分方程满足所给初始条件的特解. ,0)ln (ln =-+dx x y xdy x .1==ex y解 将方程标准化为,1ln 1xy xx y =+'于是⎪⎪⎭⎫⎝⎛+=⎰⎰⎰-C dx exey xx dxx x dx ln ln 1⎪⎭⎫⎝⎛+=⎰-C dx exexx ln ln ln ln 1.ln21ln 12⎪⎭⎫⎝⎛+=C x x由初始条件,1==ex y得,21=C 故所求特解为.ln 1ln 21⎪⎭⎫ ⎝⎛+=x x y 例10(E07)在一个石油精炼厂,一个存储罐装8000L 的汽油,其中包含100g 的添加剂. 为冬季准备,每升含2g 添加剂的石油以40L/min 的速度注入存储罐. 充分混合的溶液以45L/min 的速度泵出. 在混合过程开始后20分钟罐中的添加剂有多少?解 令y 是在时刻t 罐中的添加剂的总量. 易知100)0(=y . 在时刻t 罐中的溶液的总量 ()()t t t V 5800045408000-=-+= 因此,添加剂流出的速率为()()()()tt y tt y t V t y 58000454558000-=⋅-=⋅溶液流出的速率添加剂流入的速率80402=⨯,得到微分方程 ty d td y 580004580--=即805800045=⋅-+y td t d y于是,所求通解为()()9580004558000451600101600080-+-=⎪⎪⎭⎫ ⎝⎛+⎰⋅⎰=---⎰t C t C dt e e y dt t dt t由100)0(=y 确定C ,得()()016000010160009=-+⨯-C ,8160010=C ,故初值问题的解是()()9816001600101016000-+-=t t y ,所以注入开始后20分钟时的添加剂总量是 ()()58.1512160020160010201016000)20(98≈-+⨯-=y g.注:液体溶液中(或散布在气体中)的一种化学品流入装有液体(或气体)的容器中,容器中可能还装有一定量的溶解了的该化学品. 把混合物搅拌均匀并以一个已知的速率流出容器. 在这个过程中,知道在任何时刻容器中的该化学品的浓度往往是重要的. 描述这个过程的微分方程用下列公式表示:容器中总量的变化率=化学品进入的速率—化学品离开的速率.课堂练习1. 求微分方程yx y y ydxdy sin 2sin cos cos -=的通解.2. 求线性微分方程23x y y =+'满足条件00==x y 的特解.雅各布.伯努利(Jacob Bermoulli ,1654~1705)伯努利瑞士数学、力学、天文学家。