论文部分内容阅读
【摘要】本文主要介绍了黎卡提(Riccati)微分方程运用特解来求解方程通解的几种特殊情形,即通过特解个数的变化来求解出方程的通解。
【关键词】黎卡提(Riccati);微分方程;伯努利方程;特解;通解
【中图号】S777【文献标示码】A【文章编号】1005-1074(2009)01-0146-02
1引言
若一阶微分方程可写成如下形式:
dydx=P(x)+Q(x)y+R(x)y2(1)
其中P(x),Q(x),R(x)在区间I上连续,且R(x)≠0,则方程(1)叫做黎卡提(Riccati)微分方程。
刘维尔在1846年证明出了少数情形外,它的解不能用初等函数的有限次积分以及有限次代数运算得到。但作为黎卡提(Riccati)微分方程的特殊情形,包含着以下几种可求解的情形。
当P(x),Q(x),R(x)为常数时,方程(1)可化为变量分离方程
dyay2+by+c=dx
这里a,b,c是常数。
当R(x)=0时,方程(1)化为线形方程
dydx=P(x)+Q(x)y
当P(x)=0时,方程(1)化为伯努利方程
dydx=Q(x)y+R(x)y2
黎卡提(Riccati)微分方程的一个重要特点就在于,当已知方程的特解时,就可以比较容易地求出其通解。
本文主要介绍了在已知特解的情形下来求解黎卡提(Riccati)微分方程的通解。
2主要定理的证明及应用
定理1 已知方程(1)的一个特解y1,则方程(1)可化为伯努利方程来求解,其通解为
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
证明 由题设知y1为特解
作变换y=y1+u,(其中u是关于自变量x的新的未知函数)
则有y'=y'1+u',代入(1)中得
y'1+u'=P(x)+Q(x)(y1+u)+R(x)(y1+u)2
整理即得
u'=R(x)u2+[Q(x)+2R(x)y1]u-y'1+R(x)y12+Q(x)y1+P(2)
而y1为(1)的特解,有
y'1=R(x)y12+Q(x)y1+P(x)
故(2)式整理得
u'=R(x)u2+[Q(x)+2R(x)y1]u (3)
(3)式是关于u的一伯努利方程,可解。
令v=1u,代入(3)式中化简即得
dvdx=[Q(x)+2R(x)y1]v+R(x)
对上式运用积分公式解的通解为
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
例1 求解下列黎卡提(Riccati)微分方程
dydx+x-2x-x2y+1x2-x3y2=0
已知y1=x为其特解。
解 由于y1=x是方程的特解,根据上面定理1,其通解公式为
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
故得到
1y-x=e∫(x-2x-x2+21x2-x3x)dx[∫1x2-x3e-∫(x-2x-x2+21x2-x3x)dxdx+c]化简得
(y-x)(cx+1)=x(x-1)为方程的通解。
定理2 若已知方程(1)的两个特解y1,y2(y1≠y2),则方程(1)仅用一次积分就能求解,其通解为
1y-y1=e∫(Q+2Ry1)dxc+1y2-y1.
证明 已知两特解为y1,y2(y1≠y2),设y=y1+u,
由定理1的证明过程可得到
v=1y-y1
v'=[Q(x)+2R(x)y1]v+R(x)(4)
又由假设得1y1-y2是方程(4)的特解。
一方面因为此方程是线性的,所以方程(1)的通解只要由v'-[Q(x)+2R(x)y1]v=0的通解加上一个特解即可。故方程(1)仅用一次积分即可求解。
而v'-[Q(x)+2R(x)y1]v=0的通解为
v=e∫(Q+2Ry1)dxc
即1y-y1=e∫(Q+2Ry1)dxc
所以方程(1)的通解为
1y-y1=e∫(Q+2Ry1)dxc+1y2-y1.
例2 求解下列黎卡提(Riccati)微分方程
dydx+x-2x-x2y+1x2-x3y2=0
已知y1=x,y2=0为其特解.
解 由题意y1=x,y2=0为方程的特解,根据上面定理2,其通解公式为
1y-y1=e∫(Q+2Ry1)dxc+1y2-y1.
即得1y-x=e∫(x-2x-x2+21x2-x3x)dxc-1x.
化简即有
1y-x=cx-1-1x
所以y(x-1)=c(y-x)x为方程的解.
定理3 若已知方程(1)的3个互异特解y1,y2,y3则不需要积分计算就能求得其通解,其通解为
y3-y1y3-y2=y-y1y-y2C.
其中C为任一常数.
证明 由定理1中黎卡提(Riccati)方程之解法所示的通解
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
令右边e∫(Q+2Ry1)dx=h(x),∫Re-∫(Q+2Ry1)dxdx=g(x)
即1y-y1=h(x)g(x)+h(x)c
所以 y=h(x)y1c+h(x)g(x)y1+1h(x)c+h(x)g(x)
又令h(x)y1=f1(x),h(x)g(x)y1+1=f2(x),h(x)=f3(x),h(x)g(x)=f4(x),
故 y=cf1(x)+f2(x)cf3(x)+f4(x) (a)
乃是关于常数C的一次分数函数。
又设c1,c2,c3为对应于上式特解y1,y2,y3
故有 y=cif1(x)+f2(x)cif3(x)+f4(x)(i=1,2,3) (b)
由(a)和(b)式消去f1(x),f2(x),f3(x),f4(x)得到
y3-y1y3-y2/y3-yy3-y=c3-c1c3-c2/c3-cc3-c
而其中c,c1,c2,c3均为任一常数,故令
c3-c1c3-c2/c3-cc3-c=C (C为任一常数)
对y求解即得到方程(1)的通解
y3-y1y3-y2=y-y1y-y2C.
例3 求解下列黎卡提(Riccati)微分方程
dydx+x-2x-x2y+1x2-x3y2=0
已知y1=0,y2=x,y3=x3为其特解。
解 由已知条件得y1=0,y2=x,y3=x3为方程的特解,根据定理3,不用积分便可求得方程的通解,其通解公式为
y3-y1y3-y2=y-y1y-y2C.
以代入y1=0,y2=x,y3=x3并对y求解得到
x2-0x2-x=y-0y-xC
即得(y-x)x=y(x-1)C为方程的通解。
4参考文献
[1]叶彦谦常微分方程讲义[M]北京:人民教育出版社,1979:45-57
[2]矢野健太郎微分方程的基本原理及习题详解[M]北京:晓圆出版社,1994:42-45
[3]丁同仁,李承治.常微分方程教程[M]北京:高等教育出版社,2002:61-68
[4]汤正谊微分方程解题分析[M]江苏:江苏科学技术出版社,1993:106-107
[5]王建锋浅谈黎卡提方程的求解[J]数学理论与应用,2002,22(3):107-109
注:本文中所涉及到的图表、注解、公式等内容请以PDF格式阅读原文。
【关键词】黎卡提(Riccati);微分方程;伯努利方程;特解;通解
【中图号】S777【文献标示码】A【文章编号】1005-1074(2009)01-0146-02
1引言
若一阶微分方程可写成如下形式:
dydx=P(x)+Q(x)y+R(x)y2(1)
其中P(x),Q(x),R(x)在区间I上连续,且R(x)≠0,则方程(1)叫做黎卡提(Riccati)微分方程。
刘维尔在1846年证明出了少数情形外,它的解不能用初等函数的有限次积分以及有限次代数运算得到。但作为黎卡提(Riccati)微分方程的特殊情形,包含着以下几种可求解的情形。
当P(x),Q(x),R(x)为常数时,方程(1)可化为变量分离方程
dyay2+by+c=dx
这里a,b,c是常数。
当R(x)=0时,方程(1)化为线形方程
dydx=P(x)+Q(x)y
当P(x)=0时,方程(1)化为伯努利方程
dydx=Q(x)y+R(x)y2
黎卡提(Riccati)微分方程的一个重要特点就在于,当已知方程的特解时,就可以比较容易地求出其通解。
本文主要介绍了在已知特解的情形下来求解黎卡提(Riccati)微分方程的通解。
2主要定理的证明及应用
定理1 已知方程(1)的一个特解y1,则方程(1)可化为伯努利方程来求解,其通解为
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
证明 由题设知y1为特解
作变换y=y1+u,(其中u是关于自变量x的新的未知函数)
则有y'=y'1+u',代入(1)中得
y'1+u'=P(x)+Q(x)(y1+u)+R(x)(y1+u)2
整理即得
u'=R(x)u2+[Q(x)+2R(x)y1]u-y'1+R(x)y12+Q(x)y1+P(2)
而y1为(1)的特解,有
y'1=R(x)y12+Q(x)y1+P(x)
故(2)式整理得
u'=R(x)u2+[Q(x)+2R(x)y1]u (3)
(3)式是关于u的一伯努利方程,可解。
令v=1u,代入(3)式中化简即得
dvdx=[Q(x)+2R(x)y1]v+R(x)
对上式运用积分公式解的通解为
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
例1 求解下列黎卡提(Riccati)微分方程
dydx+x-2x-x2y+1x2-x3y2=0
已知y1=x为其特解。
解 由于y1=x是方程的特解,根据上面定理1,其通解公式为
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
故得到
1y-x=e∫(x-2x-x2+21x2-x3x)dx[∫1x2-x3e-∫(x-2x-x2+21x2-x3x)dxdx+c]化简得
(y-x)(cx+1)=x(x-1)为方程的通解。
定理2 若已知方程(1)的两个特解y1,y2(y1≠y2),则方程(1)仅用一次积分就能求解,其通解为
1y-y1=e∫(Q+2Ry1)dxc+1y2-y1.
证明 已知两特解为y1,y2(y1≠y2),设y=y1+u,
由定理1的证明过程可得到
v=1y-y1
v'=[Q(x)+2R(x)y1]v+R(x)(4)
又由假设得1y1-y2是方程(4)的特解。
一方面因为此方程是线性的,所以方程(1)的通解只要由v'-[Q(x)+2R(x)y1]v=0的通解加上一个特解即可。故方程(1)仅用一次积分即可求解。
而v'-[Q(x)+2R(x)y1]v=0的通解为
v=e∫(Q+2Ry1)dxc
即1y-y1=e∫(Q+2Ry1)dxc
所以方程(1)的通解为
1y-y1=e∫(Q+2Ry1)dxc+1y2-y1.
例2 求解下列黎卡提(Riccati)微分方程
dydx+x-2x-x2y+1x2-x3y2=0
已知y1=x,y2=0为其特解.
解 由题意y1=x,y2=0为方程的特解,根据上面定理2,其通解公式为
1y-y1=e∫(Q+2Ry1)dxc+1y2-y1.
即得1y-x=e∫(x-2x-x2+21x2-x3x)dxc-1x.
化简即有
1y-x=cx-1-1x
所以y(x-1)=c(y-x)x为方程的解.
定理3 若已知方程(1)的3个互异特解y1,y2,y3则不需要积分计算就能求得其通解,其通解为
y3-y1y3-y2=y-y1y-y2C.
其中C为任一常数.
证明 由定理1中黎卡提(Riccati)方程之解法所示的通解
1y-y1=e∫(Q+2Ry1)dx〔∫Re-∫(Q+2Ry1)dxdx+c〕.
令右边e∫(Q+2Ry1)dx=h(x),∫Re-∫(Q+2Ry1)dxdx=g(x)
即1y-y1=h(x)g(x)+h(x)c
所以 y=h(x)y1c+h(x)g(x)y1+1h(x)c+h(x)g(x)
又令h(x)y1=f1(x),h(x)g(x)y1+1=f2(x),h(x)=f3(x),h(x)g(x)=f4(x),
故 y=cf1(x)+f2(x)cf3(x)+f4(x) (a)
乃是关于常数C的一次分数函数。
又设c1,c2,c3为对应于上式特解y1,y2,y3
故有 y=cif1(x)+f2(x)cif3(x)+f4(x)(i=1,2,3) (b)
由(a)和(b)式消去f1(x),f2(x),f3(x),f4(x)得到
y3-y1y3-y2/y3-yy3-y=c3-c1c3-c2/c3-cc3-c
而其中c,c1,c2,c3均为任一常数,故令
c3-c1c3-c2/c3-cc3-c=C (C为任一常数)
对y求解即得到方程(1)的通解
y3-y1y3-y2=y-y1y-y2C.
例3 求解下列黎卡提(Riccati)微分方程
dydx+x-2x-x2y+1x2-x3y2=0
已知y1=0,y2=x,y3=x3为其特解。
解 由已知条件得y1=0,y2=x,y3=x3为方程的特解,根据定理3,不用积分便可求得方程的通解,其通解公式为
y3-y1y3-y2=y-y1y-y2C.
以代入y1=0,y2=x,y3=x3并对y求解得到
x2-0x2-x=y-0y-xC
即得(y-x)x=y(x-1)C为方程的通解。
4参考文献
[1]叶彦谦常微分方程讲义[M]北京:人民教育出版社,1979:45-57
[2]矢野健太郎微分方程的基本原理及习题详解[M]北京:晓圆出版社,1994:42-45
[3]丁同仁,李承治.常微分方程教程[M]北京:高等教育出版社,2002:61-68
[4]汤正谊微分方程解题分析[M]江苏:江苏科学技术出版社,1993:106-107
[5]王建锋浅谈黎卡提方程的求解[J]数学理论与应用,2002,22(3):107-109
注:本文中所涉及到的图表、注解、公式等内容请以PDF格式阅读原文。