差分方程x_n_1_x_n_e_r_n_1_x_n_k_正解的渐近性

差分方程x_n_1_x_n_e_r_n_1_x_n_k_正解的渐近性
差分方程x_n_1_x_n_e_r_n_1_x_n_k_正解的渐近性

差分方程x n +1=

x n e r n (1-x n -k )

正解的渐近性

刘玉记

(岳阳师院数学系,岳阳,414000)

摘 要

本文研究时滞Log istic 方程

x n +1=x n e r n

(1-x

n-k

)

,n =0,1,2, (*)

其中r n 是非负实数列, (0,1],k 为非负整数,获得了保证方程(*)的每一正解趋于1的一些充分条件,推

广和改进了文[1]的结果.

关键词 全局吸引性,Logistic 方程,充分条件主题分类 AMS(1991) 39A10,39A11中图分类号:O 175

1 引言和主要结果

本文考虑Logistic 方程

x n+1=x n e r n

(1-x n-k )

n =0,1,2, (1)

其中{r n }是非负实数列, [0,1],k 是非负整数.当 =1,r n !r 时,文[1]证明了当r (k +1)?1时,(1)的每一正解x n 满足lim x n =1.

本文我们证明下面的定理.定理1 设

#?

n=0

-n r n <

?,0? <1

(2) lim sup

#n i =n-k

-k-2-i

r i <1+

12

-1

(1+1k +1)?2

(3)则(1)的每一正解x n 满足lim x n =1.

(4)

定理2 设

#?

n=k

r n =

?, =1(5)lim sup

#n

i=n -k

r i <

32+12(k +1)

(6)

则(1)的每一正解x n 满足(4)式.

推论1 设

#?

n=k r n =

?,lim sup r n <2,则

第17卷第3期2000年9月 经 济 数 学MATHE MATICS ECONOMICS

Vo1 17 No.3

Sep.2000

湖南省教委科研基金资助:99C12

收稿日期:1998-1-09

x n+1=x n e r n

(1-x n )

(7)

的每一正解x n 满足(4)式.

定理的证明依赖于第2节给出的引理.显然,我们的结果改进了文[1]的结果.方程(1)的初始条件

x i =a i %0,i =-k,-k +1, ,0,a 0>0,

(8)

易知(1)与(8)有唯一正解x n .这里说x n 是(1)的解指的是数列{x n }(n %-k ),x n 满足(1)式.本文约定j ?i -1时,#j

i

r n =

0,由于(1)有实际背景,因此,本文结果对差分方程理论及实际

应用均有一定意义.

2 引 理

引理1 设f (x )=ux -vx 2,u >0,v >0,则存在0

max 0?x ?1

f (x )?max {u -v,r }.证明 u %2v 时,则f &(x )=u -2vx %0(x [0,1])故f (x )?f (1)=u -v.当u <2v 时,f (x )?f (u

2v )=u 24v

令ln x n =y n ,(1)式化为

y n+1-2y n =r n (1-e y n-k )

(9)引理2 设y n 是(9)的振动解,y N 2>0是任一极大项,N 2>2k +1,满足

y n+1- y n ?cr n ,n ?N 2

(10)

则存在实数 [N 2-k -1,N 2-k ],使得当n ?N 2-1时

y n+1- y n ?cr n N 2-n [#N 2-k-1

i =n-k

N 2-k-1-i

r i -r N 2-k-1(N 2-k - )](11)

证明 由于y N 2是极大项,故y N 2-1?y N 2%y N 2+1,从而 y N 1-1? y N 2?y N 2% y N 2% y N 2+1,于是

0?y N 2- y N 2-1=r N 2-1(1-e y N 2-1-k )

从而y N 2-k-1?0,同理y N 2-k %0.由中值定理,存在 [N 2-k -1,N 2-k ],使得

y N 2-k-1+(y N 2-k - y N 2-k-1( -N 2+k +1)=0(12)

由(10)式,当n ?N 2-k -1时

y N 2-k-1- N 2

-k-n-1

y n ?c

#

N 2

-k-2

i=n

r i N 2

-k-2-i

而结合(12)式立得

-

N 2-k-n

y n ?(y N 2-k - y N 2-

k -1)( -N 2+k +1)+c

#

N 2

-k-2

i=n

N 2-k-1-i

r i

从而n ?N 2-1时

- N 2-n

y n-k

?-cr N 2-k-1(N 2-k - )+c

#

N 2

-k-1

i=n -k

N 2

-k -1-i

r i

又由(9)式,y n+1- y n ?-r n y n-k ,结合前一式立得(11)式

?67?第3期

刘玉记:差分方程x n +1=x n e r n (1-x n-k )

正解的渐近性

引理3 设y n 满足引理2中条件,又存在实数a 1>0使得

#n

i=n-k

-k-2-i

r i ?a 1(n %M )(13)

则存在实数0

,使得

y N 2?c max {a 1-12(1+1k +1

) -1

,r }

(14)

证明分两种情形讨论.

情形1 d =#

N 2

-1

n=N 2

-k

r n +r N 2-k-1(N 2-k - )?1.由(10)及(11)式

y N 2= k

y N 2-k +

#

N 2

-1

i=N 2

-k

N 2-i -1

(y i+1- y i )

= k (y N 2-k - y N 2-k-1(N 2-k - )+

#

N 2

-1

n=N 2

-k

(y n+1- y n )

N 2-1-n

?

k

cr N 2-k-1

[

#

N 2

-k-1

i=N 2

-2k-1

N 2-k-1-i

r i -r N 2-k-1(N 2-k - )](N 2-k - )

+c

#

N 2

-1n=N 2

-k

N 2

-1-n

r n

N 2

-k-n+k [

#

N 2

-k-1

i=n -k

N 2-k -1-i

r i -r N 2-k-1(N 2-k - )]

=c -1

r N 2-k-1[

#

N 2

-k-1

i=N 2

-2k-1

N 2-k-1-i

r i -r N 2-k-1(N 2-k - ](N 2-k - )

+c

#

N 2

-1

N =N 2

-k

-1

r n [

#

N 2

-k-1

i =n-k

N 2-k-1-i

r i -r N 2-k-1(N 2-k - )]

=c[r N 2-k -1

#

N 2

-k-1

i=N 2

-2k-1

N 2

-k-2-i

r i - -1r 2

N 2-k-1(N 2-k - )](N 2-k - )

+c[

#

N 2

-1N =N 2

-k

r n ([

#

N 2

-k-1

i=n-k

N 2-k-1-i

r i - -1

r N 2-k-1(N 2-

k - ))]

?c[r N 2-K-1a 1- -1r 2

-N 2-k-1(N 2-k - )](N k -k - )+

c[

#

N 2

-1

n=N 2

-k

r n (a 1-

#

n

n=N 2

-k

N 2

-k-2-i

r i - -1

r N 2-k-1(N 2-k - ))]

=a 1cd -c[ -1r 2

N 2

-k-1(N 2-k - )2

+

#

N 2

-1

n =N 2

-k

r n

#

n

i=N 2

-k

N 2-k-2-i

r i

+ -1

r N 2-k-1(N 2-k - )

#

N 2

-1

i=N 2

-k

r i ]

因i %N 2-k ,所以N 2-k -2-i ?k -k -1=-1<0,从而有

y N 2?a 1cd -c[

r 2N 2-k-1(N 2-k - )2

+

#

N 2

-1

i =N 2

-k

r n

#

n

i =N 2

-k

r i

+r N 2-k-1(N 2-k - )

#

N 2

-1

i=N 2

-k

r i )

-1

?68? 经

济数学 第17卷

由#N 2

-1

i =N 2-k r n #n

i =N 2-k r i =12(#N 2

-1

n=N 2-k r n )2+12#N 2

-1

N=N 2

-k r 2n 及#m

n=1

r 2

n %1m (#m

n=1

r n )可知,

y N 2

?a 1cd -12c -1

[(#N 2

-1

n=N 2

-k

r n +r N 2-k -1(N 2-

k - )]2+

(

#

N 2

-1

n=N 2

-k

r n )2

+

r 2N 2-k-1(N 2-

k - )2

+(

#

N 2

-1

n=N 2

-k

r n )2

+

#

N 2

-1

n=N 2

-k

r 2n ]

?c[a 1d -12d 2 -1-12(#N 2

-1

n =N 2-k

r 2n +r 2

N 2-k-1(N 2-k - )2]?c[a 1d -12

-1d 2

-121k +1(#N 2

-1

n=N 2-k r n +r N 2-k-1(N 2-

k - ))2]

=c(a 1d -12(1+1k +1) -1d 2

)由引理1立得,存在0

) -1

,使(14)式成立.

情形2 d >1.由于d =#

N 2

-1

n =N 2

-k-1

r n +r N 2-

k -1(N 2-k - -1)>1,而N 2-k - -1<

0,故存在

N 3 [N 2-k ,N 2],使得

#N 2

-1

n=N

3

r n

?1,

#

N 2

-1

n=N 3

-1

r n >1

这样就存在 (N 3-1,N 3),使得

#N 2

-1

n=N

3

r n +

r N 3-1(N 3- )=1

(15) -1

(#N 2

-1

n =N

3

r n +

r N 3-1(N 3- ))%1

(16)

我们断言

y N 2? k

y N 2-k +

#

N 3

-2

n=N 2

-k

N 2

-1-n

(y n+1- y n )+ N 2

-N

3

( -N 3+1)(y N 3- y N 3-1)

+ N 2

-N

3

(N 3- )(y N 3- y N 3-1)+

#

N 2

-1

n=N

3

N 2-1-n

(y n+1- y n )

(17)

事实上,若N 2-k +2?N 3?N 2-1,则(1)式显然成立.而结合N 3 [N 2-k ,N 2],只需说明N 3=N 2,N 2-k,N 2-k +1时,(17)成立.

N 3=N 2时,(17)式右边化为

k

y N 2-k +

#

N 2

-1

n=N 2

-k

N 2-1-n (y n+1- y n )+ N 2-N 2(y N 2- y N 2-1)

=y N 2- y N 2-1+

#

N 2

-1

n=N 2

-k

N 2

-1-n

(y n+1- y n )+ k

y N 2-k .

?

69?第3期

刘玉记:差分方程x n +1=x n e r n (1-x n-k )

正解的渐近性

又当k %2时,(17)成立.当k =1时,(17)式右边化为 y N 2-1+y N 2- y N 2-1,(17)成立.当k =0时(17)式右端化为y N 2- y N 2-1+y N 2.由于y N 2是极大项,y N 2- y N 2-1%0,故(17)成立.

N 3=N 2-k 时,(17)右边化为

k

y N 2-k +

#

N 2

-1

n=N 2

-k

N 2-1-n (y n+1- y n )+ k

(y N 2-k - y N 2-k-1).

当k %1时,化为y N 2+ k

(y N 2-k - y N 2-k-1).由于y N 2是极大项,y N 2-k-1

?0,y N 2-

k

%0,(17)

成立.而k =0时,化为y N 2+y N 2- y N 2-1.由y N 2是极大项立知(17)成立.

N 3=N 2-k +1时,(17)右边化为

k

y N 2-k +

#

N 2

-1

n=N 2

-k+1

N 2-1-n (y n+1- y n )+ k -1(y N 2-k+1- y N 2-k )

当k %2时,化为 k

y N 2-k +y N 2- k-2 y N 2-k+1+ k -1

(y N

2

-k+1-

y N 2-K ,(17)成立,而k =1时,

化为 y N 2-1+y N 2- y N 2-1,(17)式成立,k =0不可能(因N 3 [N 2-k ,N 2]).从而(17)有

y N 2? k

(y N 2-k - y N 2-k-1)(N 2-k - )+

#

N 2

-2

n=N 2

-k

N 2-1-n

(y n+1- y n ) + N 2

-N

3

( -N 3+1)(y N 3- y N 3-1)+ N 2

-N 3

(n 3- )(y N 3- N 3

-1) +

#N 2

-1

n=N

3

N 2

-1-n

(y n+1- y n )

?c[ k

r N 2-K-1(N 2-k - )+

#

N 3

-2

n=N 2

-k

N 2-1-n r n +( -N 3+1) N 2-N 3

r N 3-1]

+c[

#

N 2

-1

n=N

3

N 2

-1-n

r n

N 2

-k -n+k (

#

N 2

-k-1

i=n -k

N 2

-k-1-i

r i +r N 2-k-1(N 2-k - ))]

+c N 2

-N

3

(N 3- )

r N 3-1

N 2

-k-N 3

+k+1[

#

N 2

-k-1

i=N 3

-k-1

!2

-k-1-i

r i -r N 2-

k-1(

N 2-k - )]

=c[ k

r N 2-k-1(N 2-k - )+

#

N 3

-2

n=N 2

-k

N 2-1-n r n +( -N 3+1) N 2-N

3r N 3

-1] +c

-1

#N 2

-1n=N

3

r n [#

N 2

-k-1

i=n -k

N 2-k-1-i

r i -r N 2-k-1(N 2-k - )]

+c -1

(N 3- )r N 3-1[#

N 2

-K-1

I =N 3

-k-1

N 2-k -1-i

r i -r N 2-k-1(N 2-K - )]

现在应用(16)式,并整理得

y N 2?c{

-1

#N 2

-1

n=N

3

r n [#

N 2

-k-1

i=n -k

N 2-k-1-i r i -r N 2-k-1(N 2-k - )]}

+c{[ k

r N 2-k-1(N 2-k - )+

#

N 3

-2

n=N 2

-k

N 1-1-n

r n +( -N 3+1)r N 3-1](

[

-1

#N 2

-1

n=N

3

r n +

-1

r N 3-1(N 3- )]}

?70? 经

济数学 第17卷

+c -1

(N 3- )r N 3-1[

#

N 2

-K-1

i=N 3

-k-1

N 2-k -1-i

r i -r N 2-k-1(N 2-k - )]

=c{

-1

#N 2

-1

n=N

3

r n [#

N 2

-k-1

i=n-k

N 2

-k-1-i

r i -r N 2-k-1(N 2-

k - )]

+ k

r N 2-k -1(N 2-k - )+

#

N 3

-1

n=N 2

-k

N 2

-n-1

r n +( -N 3)r N 3-1

+ -1

(N 3- )r N 3-1[#

N 2

-k-1

i=N 3

-k-1

N 2-k -1-i

r i -r N 2-k-1(N 2-k - )

+ k

r N 2-k -1(N 2-k - )+

#

N 3

-1

n=N 2

-k

N 2-1-n

r n +( -N 3)r N 3-1]}

注意 k

?1, N 2

-1-i

? N 2

-k -1-i

上式得出

y N 2?c

-1

#N 2

-1

n=N

3

r n [#N 3

-1

i=n -k

N 2

-k-1-i

r i -(N 3- )r N 3-1]

+c -1

(N 3- )r N 3-1[

#

N 3

-1

i=N 3

-1-k

N 2-k-1-i

r i -(N 3- )N 3-1]

=c

#N 2

-1

n=N

3

r n [#

N 3

-1

i=n -k

N 2

-k -2-i

r i - -1

r N 3-1(N 3- )]

+c(N 3- )r N 3-1[

#

N 3

-1

i=N 3

-k-1

N 2-k-2-i r i - -1

(N 3- )r N 3-1]?c

#N 2

-1

n=N

3

r n [a 1-#n

i=N

3

N 2

-k-2-i

r i - -1(N 3- )r N 3-1]

+c(N 3- )r N 3-1[a 1- -1

(N 3- )r N 3-1]

=c(a 1-

#N 2

-1

n=N

3

r n #n

i=N

3

N 2-k-2-i

r i - -1

(N 3- )r N 3-1

#N 2

-1

n =N

3

r n - -1

(N 3- )2r 2N 3-1)

?c(a 1-12 -1(#N 2

-1

n=N 3r n )2-12 -1#N 2

-1

n=N 3r 2n - -1

(N 3- )r N 3-1#N 2

-1

n=N

3

r n - -1

(N 3- )2r 2N 3

-1)=c(a 1-12 -1(#N 2

-1

n=N 3r n +(N 3- )(r N 3-1)2-12 -1(#N 2

-1

n=N

3r 2n +(N 3- )2r 2N 3-1)))

=c(a 1-12 -1-12 -1((N 3- )2r 2N 3-1+(#N 2

-1

n=N

3

r 2

n )))

?c(a 1-12 -1-12 -11N 2-N 3+1(#N 2

-1

n=N

3

r n +(N 3- )r 2N 3-1)2

)

?c(a 1-

12(1+1k +1

) -1

.综合情形1与情形2,存在0

) -1

使得

?

71?第3期

刘玉记:差分方程x n +1=x n e r n (1-x n-k )

正解的渐近性

y N 2?cmax {a 1-

12(1+1k +1

) -1

,r }3 定理证明

定理1的证明:

设x n 是(1)有任一正解,令ln x n =y n ,(1)化为 y n+1- y n =r n (1-e y n-k )

(19)

我们证明在条件(2)(3)(8)下,(19)的每一解y n 满足lim n )?

y n =0.情形1 y n 是最终正解.

y n+1-y n

?0

从而y n 最终单调减少,而y n 最终正,故y *=lim y n %0,若y *>0,则存在N,当n %N 时

y n >y *

2

>0,n %N 对(19)式

y n+1- y n ?-r n y n-k

(20)

对(20)式求和,

y n+1-

n+1-N-k

y N+k ?-

#

n

i=N+k

y

*

2 n-i

r i =-

n+1

2

#

n

i=N+k

-i

r i

由条件(3),令n )?,得y *?0,矛盾.lim y n =0.

情形2 y n 最终为负时,-y n 最终为正.

-y n+1-(-y n )?-y n+1- (-y n )=r n (e y n-k -1)?0

-y n 单调减少,y

*

=lim (-y n )%0,若y

*

>0,则r n =y n+1- y n

1-e y

n-l limsup r n =y *( -1)

1-e *

<0

与r n 是非负数列矛盾,故lim y n =0.

情形3 若y n 是振动解.这时,由条件(2),存在?>0使得

1<#+??1+12 -1(1+1

k +1

)

#=limsup

#n

i=N-k

-k-2-i

r i

且存在N ,当n >N 时

#n

i=n -k

-k-2-i

r i ?#+?(21)

先证y n 有上界,由x n 是(1)的正解,故

x n+1=x n e r n (1-x n-k )?x n e r

n

设ln x n =y n ,化为

y n+1- y n ?r n .

现令y n *+

k+1

是y n 的任一极大项,则0?y n *+

k

?y n *+k+1%y n *+k+2,由于

?72? 经

济数学 第17卷

y n *+k+1-

y n *+k =r n *+k (1-e y n *

y n *+k+1-y n *+k

?y n *+k+1- y n *+

k

知y n *?0.同理,y n *+1%0.从n *到n *+k 求和得

y n

*

+k +1- k+1

y n

*

?

#n *

+k

i =n

*

n *

+k-i

r i

y n *+

k+1

? k+1y n *

+

#

n *

+k

i=n

*

-i

r i

? k+2

(#+?)(由(21)式).

故y n 有上界.再证y n 有下界.由y n 有上界,设y n ?M(M 是常数),故由(19)式

y n+1- y n %r n (1-e M )

设y n *+

k +1

是y n 的任一极小项,则0%y n *+k %y n *+

k+1

?y n *+k+2.又由

y n *+k+1- y n *+k =r n *+k (1-e y n *)

0%y n *+

k +1-

y n *+

k

%y n *+k+1- y n *+

k

故y n *%0.同理,y n *+1?0.对y n+1- y n %r n (1-e M ).从n *到n *+k 求和得

y n

*

+k+1- k+1

y n

*

%(1-e M

)

#

n *

+k

i=n *

n *

+k-i

r i %(1-e M ) k+2

(#+?)

由于y n *+

k+1

是任一极小项,故y n 有下界.从而y n 有界.现令:

u =limsup y n , v =liminf y n

(22)

易知:-?N ,当n >N 1时,

?1=v -%

y n+1- y n ?r n (1-e ?

1)(23)y n+1- y n %r n (1-e ?2)

(24)

由于y n 振动,可选单增整数列{n i }n i >N

*

+2k +1,y n i >0,y n *>y n *-1,lim i )?

y n i =u 由(23)及条件(2),经引理3

y n i ?(1-e ?1)max {#+?-12

-1

(1+1k +1),r }

由于#+?-12 -1(1-1k +1

)<1,r <1,故y n i ?1-e ?1.令i )?,得u ?1-e ?

1.又令%

)0,得

u ?1-e v

(25)

同理,考虑-y n 时,可得

v %1-e u (26)

若v ?0,则v <0,由(25)(26)式知:1-e v %ln (1-v ),f (v)=1-e v -ln (1-v)%0,v <0,又有f &(v)=-e v

+11-v =1-(1-v )e v

1-v

,令g (v)=1-(1-v)e v ,g (0)=0,g &(v)=e v

-(1-v)e v =ve v

<0,故g (v )>0,(v <0),于是f &(v)>0,从而f (v)是严格单增函数,则当v <0时,f (v)<0,即

?

73?第3期

刘玉记:差分方程x n +1=x n e r n (1-x n-k )

正解的渐近性

1-e v -ln (1-v )<0这与前一式矛盾,故v =0,从而推出u =0,证毕.

定理2的证明与定理1的证明相似,略.

注:运用与本文相类似方法,可以得出差分方程x n+1=x n

f (x n-k ,n)的正解的渐近性,限于篇幅,略.

参 考 文 献

[1] Kocic V.Lj.and Ladas G.,Global attracti vity in nonli near delay difference eq uations,Pr oc.Amer.Math.Soc ,115

(1992),1083-1088.[2] Zhang B.G.and Gopalsmy K.,Global attractivi ty in delay log i stic equati on with variable parameters,Ma th.Pr oc,

Camb.phil.Soc ,107(1990),579-590.

GLOBAL ATTRACTIVITY IN LOGISTIC EQUATION

Liu Yuji

(Yue y an g Tea che r s College,Yueyang,414000)

Abstract This paper s tudies the global attractivity of positive solution of logistic equation

x n+1=x n e r n

(1-x n-k )

,n=0,1,2,3, (*)

where r n is a sequence of nonegative real number. (0,1].k %0is integer,Follow theorems are proved.

Th1.If 0< <1and

limsup m,n

#n

i=n-k

m-k-2-i

r i <1+

12

-1

(1+1k+1)?2

#?

i=1

-i

r i

Then every positive solution of (*)tends to 1ans n )?.

Th2If =1and limsup #n

i=n-k

r i <

32+

12(k+1)

#?

i=1

r i =?

Then every positive solution of (*)tends to 1as n )?..Keywords Global Attractivity,Logistic equation..

?74? 经

济数学 第17卷

(完整版)大连理工大学高等数值分析抛物型方程有限差分法

抛物型方程有限差分法 1. 简单差分法 考虑一维模型热传导方程 (1.1) )(22x f x u a t u +??=??,T t ≤<0 其中a 为常数。)(x f 是给定的连续函数。(1.1)的定解问题分两类: 第一,初值问题(Cauchy 问题):求足够光滑的函数()t x u ,,满足方程(1.1)和初始条件: (1.2) ()()x x u ?=0,, ∞<<∞-x 第二,初边值问题(也称混合问题):求足够光滑的函数()t x u ,,满足方程(1.1)和初始条件: ()13.1 ()()x x u ?=0,, l x l <<- 及边值条件 ()23.1 ()()0,,0==t l u t u , T t ≤≤0 假定()x f 和()x ?在相应的区域光滑,并且于()0,0,()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。

现在考虑边值问题(1.1),(1.3)的差分逼近 取 N l h = 为空间步长,M T = τ为时间步长,其中N ,M 是 自然数, jh x x j ==, ()N j ,,1,0Λ=; τ k y y k ==, ()M k ,,1,0Λ= 将矩形域G {}T t l x ≤≤≤≤=0;0分割成矩形网格。其中 ()j i y x ,表 示网格节点; h G 表示网格内点(位于开矩形G 中的网格节点)的集合; h G 表示位于闭矩形G 中的网格节点的集合; h Γ表示h G -h G 网格边界点的集合。 k j u 表示定义在网点()k i t x ,处的待求近似解,N j ≤≤0,M k ≤≤0。 注意到在节点()k i t x ,处的微商和差商之间的下列关系 ((,)k j k j u u x t t t ????≡ ? ????): ()() ()ττ O t u t x u t x u k j k j k j +??? ????=-+,,1 ()() ()2112,,ττ O t u t x u t x u k j k j k j +??? ????=--+ ()()()h O x u h t x u t x u k j k j k j +??? ????=-+,,1 ()() ()h O x u h t x u t x u k j k j k j +??? ????=--,,1 ()() ()2112,,h O x u h t x u t x u k j k j k j +??? ????=--+ ()()() ()2 222 11,,2,h O x u h t x u t x u t x u k j k j k j k j +???? ????=+--+ 可得到以下几种最简差分格式

差分方程的解法

1、常系数线性差分方程的解 方程( 8)其中为常数,称方程(8)为常系数线性方程。 又称方程(9) 为方程(8)对应的齐次方程。 如果(9)有形如的解,带入方程中可得: (10) 称方程(10)为方程(8)、(9)的特征方程。 显然,如果能求出(10)的根,则可以得到(9)的解。 基本结果如下: (1)若(10)有k个不同的实根,则(9)有通解: , (2)若(10)有m重根,则通解中有构成项: (3)若(10)有一对单复根,令:,,则(9)的通解中有构成项: (4)若有m 重复根:,,则(9)的通项中有成项:

综上所述,由于方程(10)恰有k 个根,从而构成方程 (9)的通解中必有k个独立的任意常数。通解可记为: 如果能得到方程(8)的一个特解:,则(8)必有通解: + (11) (1)的特解可通过待定系数法来确定。 例如:如果为n 的多项式,则当b不是特征根时,可设成形如形式的特解,其中为m次多项式;如果b是r重根时,可设特解:,将其代入(8)中确定出系数即可。 2、差分方程的z变换解法 对差分方程两边关于取Z变换,利用的Z 变换F(z)来表示出的Z变换,然后通过解代数方程求出F(z),并把F(z)在z=0的解析圆环域中展开成洛朗级数,其系数就是所要求的 例1设差分方程,求 解:解法1:特征方程为,有根: 故:为方程的解。 由条件得: 解法2:设F(z)=Z(),方程两边取变换可得:

由条件得 由F(z)在中解析,有 所以, 3、二阶线性差分方程组 设,,形成向量方程组 (12)则 (13)(13)即为(12)的解。 为了具体求出解(13),需要求出,这可以用高等代数的方法计算。常用的方法有: (1)如果A为正规矩阵,则A必可相似于对角矩阵,对角线上的元素就是A的特征值,相似变换矩阵由A的特征向量构成:。 (2)将A 分解成为列向量,则有 从而,

【毕业设计(论文)】二维热传导方程有限差分法的MATLAB实现

第1章前言 1.1问题背景 在史策教授的《一维热传导方程有限差分法的MATLAB实现》和曹刚教授的《一维偏微分方程的基本解》中,对偏微分方程的解得MATLAB实现问题进行过研究,但只停留在一维中,而实际中二维和三维的应用更加广泛。诸如粒子扩散或神经细胞的动作电位。也可以作为某些金融现象的模型,诸如布莱克-斯科尔斯模型与Ornstein-uhlenbeck过程。热方程及其非线性的推广形式也被应用与影响分析。 在科学和技术发展过程中,科学的理论和科学的实验一直是两种重要的科学方法和手段。虽然这两种科学方法都有十分重要的作用,但是一些研究对象往往由于他们的特性(例如太大或太小,太快或太慢)不能精确的用理论描述或用实验手段来实现。自从计算机出现和发展以来,模拟那些不容易观察到的现象,得到实际应用所需要的数值结果,解释各种现象的规律和基本性质。 科学计算在各门自然科学和技术科学与工程科学中其越来越大的作用,在很多重要领域中成为不可缺少的重要工具。而科学与工程计算中最重要的内容就是求解科学研究和工程技术中出现的各种各样的偏微分方程或方程组。 解偏微分方程已经成为科学与工程计算的核心内容,包括一些大型的计算和很多已经成为常规的计算。为什么它在当代能发挥这样大的作用呢?第一是计算机本身有了很大的发展;第二是数值求解方程的计算法有了很大的发展,这两者对人们计算能力的发展都是十分重要的。 1.2问题现状 近三十年来,解偏微分方程的理论和方法有了很大的发展,而且在各个学科技术的领域中应用也愈来愈广泛,在我国,偏微分方程数值解法作为一门课程,不但在计算数学专业,而且也在其他理工科专业的研究生的大学生中开设。同时,求解热传导方程的数值算法也取得巨大进展,特别是有限差分法方面,此算法的特点是在内边界处设计不同于整体的格式,将全局的隐式计算化为局部的分段隐式计算。而且精度上更好。 目前,在欧美各国MATLAB的使用十分普及。在大学的数学、工程和科学系科,MATLAB

常微分方程和偏微分方程的数值解法教学大纲

上海交通大学致远学院 《常微分方程和偏微分方程的数值解法》教学大纲 一、课程基本信息 课程名称(中文):常微分方程和偏微分方程的数值解法 课程名称(英文):Numerical Methods for Ordinary and Partial Differential Equations 课程代码:MA300 学分 / 学时:4学分 / 68学时 适用专业:致远学院与数学系相关专业 先修课程:偏微分方程,数值分析 后续课程:相关课程 开课单位:理学院数学系计算与运筹教研室 Office hours: 每周二19:00—21:00,地点:数学楼1204 二、课程性质和任务 本课程是致远学院和数学系应用数学和计算数学方向的一门重要专业基础课程,其主要任务是通过数学建模、算法设计、理论分析和上机实算“四位一体”的教学方法,使学生掌握常微分方程与偏微分方程数值解的基本方法、基本原理和基本理论,进一步提升同学们利用计算机解决实际问题的能力。在常微分方程部分,将着重介绍常微分方程初值问题的单步法,含各类Euler方法和Runge-Kutta方法,以及线性多步法。将简介常微分方程组和高阶常微分方程的数值方法。在偏微分方程部分,将系统介绍求解椭圆、双曲、抛物型方程的差分方法的构造方法和理论分析技巧,对于椭圆型方程的边值问题将介绍相应变分原理与有限元方法。将在课堂上实时演示讲授的核心算法的计算效果,以强调其直观效果与应用性。本课程重视实践环节建设,学生要做一定数量的大作业。 三、教学内容和基本要求 第一部分:常微分方程数值解法 1 引论 1.1回顾:一阶常微分方程初值问题及解的存在唯一性定理

【文献综述】热传导方程差分格式的收敛性和稳定性

文献综述 信息与计算科学 热传导方程差分格式的收敛性和稳定性在实际研究物理问题过程中, 往往能给出问题相应的数学表达式, 但是由于实际物理问题的复杂性, 它的解却一般不容易求出. 由此计算物理应运而生, 计算物理是以计算机为工具, 应用数学的方法解决物理问题的一门应用性学科, 是物理、数学和计算机三者结合的交叉性学科. 它产生于二战期间美国对核武器的研究, 伴随着计算机的发展而发展. 计算物理的目的不仅仅是计算, 而是要通过计算来解释和发现新的物理规律. 这一点它与传统的实验物理和理论物理并无差别, 所不同的只是使用的工具和方法. 计算物理早已与实验物理和理论物理形成三足鼎立之势, 甚至有人提出它将成为现代物理大厦的“栋梁”. 在一个物理问题中一个数值解往往比一个式子更直观, 更有价值. 在实际求解方程时, 除了一些特殊的情况下可以方便地求得其精确解外, 在一般情况下, 当方程或定解条件具有比较复杂的形式, 或求解区域具有比较复杂的形状时, 往往求不到, 或不易求到其精确解. 这就需要我们去寻找方程的近似解, 特别是数值近似解, 简称数值解. 这里主要研究的是热传导方程. 有限差分法是微分方程和积分微分方程数值解的方法. 其基本思想是把连续的定解区域用有限个离散点构成的网格来代替, 这些离散点称作网格的节点;把连续定解区域上的连续变量的函数用在网格上定义的离散变量函数来近似;把原方程和定解条件中的微商用差商来近似, 积分用积分和来近似, 于是原微分方程和定解条件就近似地代之以代数方程组, 即有限差分方程组, 解此方程组就可以得到原问题在离散点上的近似解. 然后再利用插值方法便可以从离散解得到定解问题在整个区域上的近似解. 热传导的差分法是求解热传导方程的重要方法之一. 对于差分格式的的求解, 我们首先要关注差分格式的收敛性和稳定性. 对于一个微分方程建立的各种差分格式, 为了有实用意义, 一个基本要求是它们能够任意逼近微分方程, 即相容性要求. 一个差分格式是否有用, 就要看差分方程的精确解能否任意逼近微分方程的解, 即收敛性的概念. 此外, 还有一个重要的概念必须考虑, 即差分格式的稳定性. 因为差分格式的计

差分方程的解法分析及MATLAB实现(程序)

差分方程的解法分析及MATLAB 实现(程序) 摘自:张登奇,彭仕玉.差分方程的解法分析及其MATLAB 实现[J]. 湖南理工学院学报.2014(03) 引言 线性常系数差分方程是描述线性时不变离散时间系统的数学模型,求解差分方程是分析离散时间系统的重要内容.在《信号与系统》课程中介绍的求解方法主要有迭代法、时域经典法、双零法和变换域 法[1]. 1 迭代法 例1 已知离散系统的差分方程为)1(3 1)()2(81)1(43)(-+=-+--n x n x n y n y n y ,激励信号为)()4 3()(n u n x n =,初始状态为21)2(4)1(=-=-y y ,.求系统响应. 根据激励信号和初始状态,手工依次迭代可算出24 59)1(,25)0(==y y . 利用MATLAB 中的filter 函数实现迭代过程的m 程序如下: clc;clear;format compact; a=[1,-3/4,1/8],b=[1,1/3,0], %输入差分方程系数向量,不足补0对齐 n=0:10;xn=(3/4).^n, %输入激励信号 zx=[0,0],zy=[4,12], %输入初始状态 zi=filtic(b,a,zy,zx),%计算等效初始条件 [yn,zf]=filter(b,a,xn,zi),%迭代计算输出和后段等效初始条件 2 时域经典法 用时域经典法求解差分方程:先求齐次解;再将激励信号代入方程右端化简得自由项,根据自由项形 式求特解;然后根据边界条件求完全解[3].用时域经典法求解例1的基本步骤如下. (1)求齐次解.特征方程为081432=+-αα,可算出4 1 , 2121==αα.高阶特征根可用MATLAB 的roots 函数计算.齐次解为. 0 , )4 1()21()(21≥+=n C C n y n n h (2)求方程的特解.将)()4 3()(n u n x n =代入差分方程右端得自由项为 ?????≥?==-?+-1,)4 3(9130 ,1)1()43(31)()43(1n n n u n u n n n 当1≥n 时,特解可设为n p D n y )4 3()(=,代入差分方程求得213=D . (3)利用边界条件求完全解.当n =0时迭代求出25)0(=y ,当n ≥1时,完全解的形式为 ,)4 3(213 )41()21()(21n n n C C n y ?++=选择求完全解系数的边界条件可参考文[4]选)1(),0(-y y .根据边界条件求得35,31721=-=C C .注意完全解的表达式只适于特解成立的n 取值范围,其他点要用 )(n δ及其延迟表示,如果其值符合表达式则可合并处理.差分方程的完全解为

一维导热方程 有限差分法 matlab实现

第五次作业(前三题写在作业纸上) 一、用有限差分方法求解一维非定常热传导方程,初始条件和边界条件见说明.pdf 文件,热扩散系数α=const , 22T T t x α??=?? 1. 用Tylaor 展开法推导出FTCS 格式的差分方程 2. 讨论该方程的相容性和稳定性,并说明稳定性要求对求解差分方程的影响。 3. 说明该方程的类型和定解条件,如何在程序中实现这些定解条件。 4. 编写M 文件求解上述方程,并用适当的文字对程序做出说明。(部分由网络搜索得到,添加,修改后得到。) function rechuandaopde %以下所用数据,除了t 的范围我根据题目要求取到了20000,其余均从pdf 中得来 a=0.00001;%a 的取值 xspan=[0 1];%x 的取值范围 tspan=[0 20000];%t 的取值范围 ngrid=[100 10];%分割的份数,前面的是t 轴的,后面的是x 轴的 f=@(x)0;%初值 g1=@(t)100;%边界条件一 g2=@(t)100;%边界条件二 [T,x,t]=pdesolution(a,f,g1,g2,xspan,tspan,ngrid);%计算所调用的函数 [x,t]=meshgrid(x,t); mesh(x,t,T);%画图,并且把坐标轴名称改为x ,t ,T xlabel('x') ylabel('t') zlabel('T') T%输出温度矩阵 dt=tspan(2)/ngrid(1);%t 步长 h3000=3000/dt;

h9000=9000/dt; h15000=15000/dt;%3000,9000,15000下,温度分别在T矩阵的哪些行T3000=T(h3000,:) T9000=T(h9000,:) T15000=T(h15000,:)%输出三个时间下的温度分布 %不再对三个时间下的温度-长度曲线画图,其图像就是三维图的截面 %稳定性讨论,傅里叶级数法 dx=xspan(2)/ngrid(2);%x步长 sta=4*a*dt/(dx^2)*(sin(pi/2))^2; if sta>0,sta<2 fprintf('\n%s\n','有稳定性') else fprintf('\n%s\n','没有稳定性') error end %真实值计算 [xe,te,Te]=truesolution(a,f,g1,g2,xspan,tspan,ngrid); [xe,te]=meshgrid(xe,te); mesh(xe,te,Te);%画图,并且把坐标轴名称改为xe,te,Te xlabel('xe') ylabel('te') zlabel('Te') Te%输出温度矩阵 %误差计算 jmax=1/dx+1;%网格点数 [rms]=wuchajisuan(T,Te,jmax) rms%输出误差

时间序列分析讲义 第01章 差分方程

第一章 差分方程 差分方程是连续时间情形下微分方程的特例。差分方程及其求解是时间序列方法的基础,也是分析时间序列动态属性的基本方法。经济时间序列或者金融时间序列方法主要处理具有随机项的差分方程的求解问题,因此,确定性差分方程理论是我们首先需要了解的重要内容。 §1.1 一阶差分方程 假设利用变量t y 表示随着时间变量t 变化的某种事件的属性或者结构,则t y 便是在时间t 可以观测到的数据。假设t y 受到前期取值1-t y 和其他外生变量t w 的影响,并满足下述方程: t t t w y y ++=-110φφ (1.1) 在上述方程当中,由于t y 仅线性地依赖前一个时间间隔自身的取值1-t y ,因此称具有这种结构的方程为一阶线性差分方程。如果变量t w 是确定性变量,则此方程是确定性差分方程;如果变量t w 是随机变量,则此方程是随机差分方程。在下面的分析中,我们假设t w 是确定性变量。 例1.1 货币需求函数 假设实际货币余额、实际收入、银行储蓄利率和商业票据利率的对数变量分别表示为t m 、t I 、bt r 和ct r ,则可以估计出美国货币需求函数为: ct bt t t t r r I m m 019.0045.019.072.027.01--++=- 上述方程便是关于t m 的一阶线性差分方程。可以通过此方程的求解和结构分析,判断其他外生变量变化对货币需求的动态影响。 1.1.1 差分方程求解:递归替代法 差分方程求解就是将方程变量表示为外生变量及其初值的函数形式,可以通过以前的数据计算出方程变量的当前值。 由于方程结构对于每一个时间点都是成立的,因此可以将(1.1)表示为多个方程: 0=t :01100w y y ++=-φφ 1=t :10101w y y ++=φφ t t =:t t t w y y ++=-110φφ 依次进行叠代可以得到: 1011211010110101)()1()(w w y w w y y ++++=++++=--φφφφφφφφ 0111122113121102)1(w w w y y φφφφφφφ++++++=- i t i i t t i i t w y y ∑∑=-=++=0 111 1 0φφφφ (1.2) 上述表达式(1.2)便是差分方程(1.1)的解,可以通过代入方程进行验证。上述通过叠代将 t y 表示为前期变量和初始值的形式,从中可以看出t y 对这些变量取值的依赖性和动态变化 过程。 1.1. 2. 差分方程的动态分析:动态乘子(dynamic multiplier) 在差分方程的解当中,可以分析外生变量,例如0w 的变化对t 阶段以后的t y 的影响。假设初始值1-y 和t w w ,,1 不受到影响,则有:

常微分方程边值问题的数值解法

第8章 常微分方程边值问题的数值解法 引 言 第7章介绍了求解常微分方程初值问题的常用的数值方法;本章将介绍常微分方程的边值问题的数值方法。 只含边界条件(boundary-value condition)作为定解条件的常微分方程求解问题称为常微分方程的边值问题(boundary-value problem). 为简明起见,我们以二阶边值问题为 则边值问题(8.1.1)有唯一解。 推论 若线性边值问题 ()()()()()(),, (),()y x p x y x q x y x f x a x b y a y b αβ'''=++≤≤?? ==? (8.1.2) 满足 (1) (),()p x q x 和()f x 在[,]a b 上连续; (2) 在[,]a b 上, ()0q x >, 则边值问题(8.1.1)有唯一解。 求边值问题的近似解,有三类基本方法: (1) 差分法(difference method),也就是用差商代替微分方程及边界条件中的导数,最终化为代数方程求解; (2) 有限元法(finite element method);

(3) 把边值问题转化为初值问题,然后用求初值问题的方法求解。 差分法 8.2.1 一类特殊类型二阶线性常微分方程的边值问题的差分法 设二阶线性常微分方程的边值问题为 (8.2.1)(8.2.2) ()()()(),,(),(), y x q x y x f x a x b y a y b αβ''-=<

差分方程求解

例题:已知差分方程51 (2)(1)()(+1)+0.5()66 x k x k x k r k r k +-++=,其中r (k )=1,k ≥0,x (0)=1, x (1)=2。 (1) 试由迭代法求其全解的前5项; (2) 分别由古典法求其零输入解、零状态解,以及全解; (3) 用Z 变换法求解差分方程。 解:注:解题过程中出现的下标“zi ”和“zs ”分别表示零输入条件和零状态条件。 1. 迭代法 题目中给出的条件仅仅是零输入初始条件,进行迭代求解时的初始条件应该是全解初始条件。 (1) 零输入初始条件 本题已给出零输入时的两个初始条件x zi (0)=1,x zi (1)=2。 (2) 零状态初始条件 取k =-2时,则51 (0)(1)(2)(1)0.5(2)66x x x r r --+-=-+-,得x zs (0)=0; 取k =-1 时,则51 (1)(0)(1)(0)0.5(1)66 x x x r r -+-=+-,求得x zs (1)=1。 (3) 全解初始条件 x (0)= x zi (0)+ x zs (0)=1; x (1)= x zi (1)+ x zs (1)=3。 (4) 根据求出的全解x (0)和x (1),利用迭代法求解 取k =0时,则51(2)(1)(0)(1)0.5(0)66x x x r r -+=+,求得23(2)6x =; 取k =1时,则51(3)(2)(1)(2)0.5(1)66x x x r r -+=+,求得151 (3)36x =; 取k =2时,则51(4)(3)(2)(3)0.5(2)66x x x r r -+=+,求得941 (4)216 x =。 2. 古典法 (1) 零输入解 令输入为零,则得齐次方程 51 (2)(1)()066 x k x k x k +-++= (a) 根据差分方程定义的算子()()n d x k x k n =+,可得它的特征方程251 066 d d -+= 求得特征根为: 112d = ,21 3 d =

差分方程方法

第四章 差分方程方法 在实际中,许多问题所研究的变量都是离散的形式,所建立的数学模型也是离散的,譬如,像政治、经济和社会等领域中的实际问题。有些时候,即使所建立的数学模型是连续形式,例如像常见的微分方程模型、积分方程模型等等,但是,往往都需要用计算机求数值解。这就需要将连续变量在一定条件下进行离散化,从而将连续型模型转化为离散型模型,因此,最后都归结为求解离散形式的差分方程解的问题。关于差分方程理论和求解方法在数学建模和解决实际问题的过程中起着重要作用。 下面就不同类型的差分方程进行讨论。所谓的差分方程是指:对于一个数列{}n x ,把数列中的前1+n 项()n i x i ,2,1,0=关联起来所得到的方程。 4.1常系数线性差分方程 4.1.1 常系数线性齐次差分方程 常系数线性齐次差分方程的一般形式为 02211=+?+++---k n k n n n x a x a x a x (4.1) 其中k 为差分方程的阶数,()k i a i ,,2,1 =为差分方程的系数,且()n k a k ≤≠0。对应的代数方程 02 211=++++--k k k k a a a λ λλ (4.2) 称为差分方程的(4.1)的特征方程,其特征方程的根称为特征根。 常系数线性齐次差分方程的解主要是由相应的特征根的不同情况有不同的形式。下面分别就特征根为单根、重根和复根的情况给出差分方程解的形式。 1. 特征根为单根 设差分方程(4.1)有k 个单特征根 k λλλλ,,,,321 ,则差分方程(4.1)的通解为 n k k n n n c c c x λλλ+++= 2211, 其中k c c c ,,,21 为任意常数,且当给定初始条件 () 0 i i x λ= ()k i ,,2,1 = (4.3) 时,可以唯一确定一个特解。 2. 特征根为重根 设差分方程(4.1)有l 个相异的特征根()k l l ≤≤1,,,,321λλλλ 重数分别为 l m m m ,,,21 且k m l i i =∑=1 则差分方程(4.1)的通解为

差分方程的解法

第三节 差分方程常用解法与性质分析 1、常系数线性差分方程的解 方程)(...110n b x a x a x a n k k n k n =+++-++ ( 8) 其中k a a a ,...,,10为常数,称方程(8)为常系数线性方程。 又称方程0...110=+++-++n k k n k n x a x a x a (9) 为方程(8)对应的齐次方程。 如果(9)有形如 n n x λ=的解,带入方程中可得: 0 ...1110=++++--k k k k a a a a λλλ (10) 称方程(10)为方程(8)、(9)的特征方程。 显然,如果能求出(10)的根,则可以得到(9)的解。 基本结果如下: (1) 若(10)有k 个不同的实根,则(9)有通解: n k k n n n c c c x λλλ+++=...2211, (2) 若(10)有m 重根λ,则通解中有构成项: n m m n c n c c λ )...(121----+++

(3)若(10)有一对单复根 βαλi ±=,令:?ρλi e ±=, αβ?βαρarctan ,22=+=,则(9)的通解中有构成项: n c n c n n ?ρ?ρsin cos 21--+ (4) 若有m 重复根:βαλi ±=,φρλi e ±=,则(9)的通项中有成 项: n n c n c c n n c n c c n m m m m n m m ?ρ?ρsin )...(cos )...(1221121---++---+++++++ 综上所述,由于方程(10)恰有k 个根,从而构成方程 (9)的通解中必有k 个独立的任意常数。通解可记为:-n x 如果能得到方程(8)的一个特解:*n x ,则(8)必有通解: =n x -n x +* n x (11) (1) 的特解可通过待定系数法来确定。 例如:如果)(),()(n p n p b n b m m n =为n 的多项式,则当b 不是特征 根时,可设成形如)(n q b m n 形式的特解,其中)(n q m 为m 次多项式;如 果b 是r 重根时,可设特解:r n n b )(n q m ,将其代入(8)中确定出系 数即可。

常微分方程数值解法

i.常微分方程初值问题数值解法 常微分方程初值问题的真解可以看成是从给定初始点出发的一条连续曲线。差分法是常微分方程初值问题的主要数值解法,其目的是得到若干个离散点来逼近这条解曲线。有两个基本途径。一个是用离散点上的差商近似替代微商。另一个是先对微分方程积分得到积分方程,再利用离散点作数值积分。 i.1 常微分方程差分法 考虑常微分方程初值问题:求函数()u t 满足 (,), 0du f t u t T dt =<≤ (i.1a ) 0(0)u u = (i.1b) 其中(,)f t u 是定义在区域G : 0t T ≤≤, u <∞上的连续函数,0u 和T 是给定的常数。我们假设(,)f t u 对u 满足Lipschitz 条件,即存在常数L 使得 121212(,)(,), [0,]; ,(,)f t u f t u L u u t T u u -≤-?∈∈-∞∞ (i.2) 这一条件保证了(i.1)的解是适定的,即存在,唯一,而且连续依赖于初值0u 。 通常情况下,(i.1)的精确解不可能用简单的解析表达式给出,只能求近似解。本章讨论常微分方程最常用的近似数值解法-差分方法。先来讨论最简单的Euler 法。为此,首先将求解区域[0,]T 离散化为若干个离散点: 0110N N t t t t T -=<< <<= (i.3) 其中n t hn =,0h >称为步长。 在微积分课程中我们熟知,微商(即导数)是差商的极限。反过来,差商就是微商的近似。在0t t =处,在(i.1a )中用向前差商 10()()u t u t h -代替微商du dt ,便得 10000()()(,())u t u t hf t u t ε=++ 如果忽略误差项0ε,再换个记号,用i u 代替()i u t 便得到 1000(,)u u hf t u -= 一般地,我们有 1Euler (,), 0,1, ,1n n n n u u hf t u n N +=+=-方法: (i.4) 从(i.1b) 给出的初始值0u 出发,由上式可以依次算出1,,N t t 上的差分解1,,N u u 。

常微分方程差分解法、入门、多解法

毕业论文 题目抛物型方程的差分解法学院数学科学学院 专业信息与计算科学 班级计算0802 学生王丹丹 学号20080901045 指导教师王宣欣 二〇一二年五月二十五日

摘要 偏微分方程的数值解法在数值分析中占有重要的地位,很多科学技术问题的数值计算包括了偏微分方程的数值解问题【1】。近三十多年来,数值解法的理论和方法都有了很大的发展,而且在各个科学技术的领域中应用也愈来愈广泛。本文的研究主要集中在依赖于时间的问题,借助于简单的常系数扩散方程,介绍抛物型方程的差分解法。本文以基本概念和基本方法为主,同时结合算例实现算法。 第一部分介绍偏微分方程及差分解法的基本概念,引入本文的研究对象——常系 数扩散方程: 2 2 ,,0 u u a x R t t x ?? =∈>?? 第二部分介绍上述方程的几种差分格式及每种格式的相容性、收敛性与稳定性。 第三部分通过算例检验每种差分格式的可行性。 关键词:偏微分方程;抛物型;差分格式;收敛性;稳定性;算例

ABSTRACT The numerical solution of partial differential equation holds an important role in numerical analysis .Many problems of compution in the field of science and techology include the numerical solution of partial differential equation. For more than 30 years, the theory and method of the numerical computation made a great development and its applications in various fields of science and technology are more and more widely. This paper focuses on the problems based on time. I will use object-constant diffusion equation to introduces the finite difference method of parabolic equation. This paper mainly focus on the basic concept ,basic method and simple numerical example. The first part of this paper introduces partial differential equations and basic concepts of finite difference method.I will introduce the object-constant diffusion equation for the first time. 2 2 ,,0 u u a x R t t x ?? =∈>?? The second part of this paper introduces several difference schemes of the above equation and their compatibility ,convergence and stability. The third part tests the accuracy of each scheme. Key words:partial differential equation;parabolic;difference scheme;convergence;stability;application

常微分方程的差分方法

第五章 常微分方程的差分方法 一、教学目标及基本要求 通过对本节课的学习,使学生掌握常微分方程、常微分方程方程组的数值解法。 二、教学内容及学时分配 本节课主要介绍常微分方程的数值解法。具体内容如下: 讲授内容:欧拉公式、改进的欧拉公式。 三、教学重点难点 1.教学重点:改进的欧拉公式、龙格库塔方法、收敛性与稳定性。 2. 教学难点:收敛性与稳定性。 四、教学中应注意的问题 多媒体课堂教学为主。适当提问,加深学生对概念的理解。 五、正文 基于数值积分的求解公式:欧拉公式、改进的欧拉公式 引 言 1.主要考虑如下的一阶常微分方程初值问题的求解: 00()(,)()y x f x y y x y '=??=? 微分方程的解就是求一个函数y=y(x),该函数满足微分方程并且符合初值条件。 2. 例如微分方程: xy'-2y=4x ;初始条件: y(1)=-3。 于是可得一阶常微分方程的初始问题 24(1)3y y x y ?'=+???=-?。 显然函数y(x)=x 2-4x 满足以上条件,因而是该初始问题的微分方程的解。

3. 但是,只有一些特殊类型的微分方程问题能够得到用解析表达式表示的函数解,而大量的微分方程问题很难得到其解析解,有的甚至无法用解析表达式来表示。因此,只能依赖于数值方法去获得微分方程的数值解。 4. 微分方程的数值解: 设微分方程问题的解y(x)的存在区间是[a,b],初始点x 0=a ,将[a,b]进行划分得一系列节点x 0 , x 1 ,...,x n ,其中a= x 0< x 1<…< x n =b 。y(x)的解析表达式不容易得到或根本无法得到,我们用数值方法求得y(x)在每个节点x k 的近似值y(x k ),即 y≈y(x k ),这样y 0 , y 1 ,...,y n 称为微分方程的数值解。 如果计算y n 时,只利用y n-1,称这种方法为单步法;如果在计算y n 时不仅利用y n-1,而且还要利用y n-2, y n-3,…, y n-r ,则称这种方法为r 步方法,也称多步法。 §5.1 欧拉方法 §5.1.1 欧拉格式 方程()(,)n n n y x f x y '=中,1()()()n n n y x y x y x h +-'≈ 1()()(,())n n n n y x y x hf x y x +≈+?1(,)n n n n y y hf x y +=+ 称为解一阶常微分方程初值问题的欧拉公式,也称显示欧拉公式。 欧拉公式的几何意义非常明显,因为微分方程的解在xoy 平面上表示一族积分曲线。用欧拉公式求数值解的几何意义如图: 容易验证,该折线各个顶点的纵坐标(1,2...)n y n =就是欧拉公式算得的近似值解,所以,欧拉方法又称为折线法。 算例:P98

差分方程

1 设一阶采样离散控制系统的差分方程为 ()()()1c k bc k r k +-= 已知输入信号()k r k a =,初始条件为()00c =,求系统的输出响应()c k 。 解:对差分方程两边进行Z 变换,得 ()()()()0zC z zc bC z R z --= ()k z R z Z a z a ??==??- 代入初始条件()00c =,得: ()()() z C z z a z b = --= 1 z z a b z a z b ??-??---?? 查表得 ()()1 k k c k a b a b = -- 2. 求解差分方程 ()()()()2413x k x k x k k δ+-++= 已知()0x k =,0k ≤, ()1,000k k k δ=?=? ≠?, 解:对差分方程两端作z 变换,得 ()()()()()()2 2 014031z X z z x zx zX z zx X z ----+=???? 已知x (0)=0,将k =-1代入差分方程得 x (1) = 0 将x (0)=0,x (1) = 0代入z 变换式,得: ()()() 2 11 43 31X z z z z z = = -+-- ()() () 1 1 2 2 1 3 lim 1lim 343 43 k k z z z z x k z z z z z z --→→=-+--+-+ =1 0.50.53k --+? 3. 求差分方程 ()()()2 1.510.50f k f k f k -+-+=的解。已知初始条件为()0.5f T -=-, ()20.75f T -=。

有限差分法解偏微分方程

有限差分法解偏微分方程综述 绪论 有限元方法最早应用于结构力学,后来随着计算机的发展慢慢用于流体力学的数值模拟。在有限元方法中,把计算域离散剖分为有限个互不重叠且相互连接的单元,在每个单元内选择基函数,用单元基函数的线形组合来逼近单元中的真解,整个计算域上总体的基函数可以看为由每个单元基函数组成的,则整个计算域内的解可以看作是由所有单元上的近似解构成。 有限差分方法(FDM)是计算机数值模拟最早采用的方法,至今仍被广泛运用。该方法将求解域划分为差分网格,用有限个网格节点代替连续的求解域。有限差分法以Taylor 级数展开等方法,把控制方程中的导数用网格节点上的函数值的差商代替进行离散,从而建立以网格节点上的值为未知数的代数方程组。该方法是一种直接将微分问题变为代数问题的近似数值解法,数学概念直观,表达简单,是发展较早且比较成熟的数值方法。 对于有限差分格式,从格式的精度来划分,有一阶格式、二阶格式和高阶格式。 从差分的空间形式来考虑,可分为中心格式和逆风格式。 考虑时间因子的影响,差分格式还可以分为显格式、隐格式、显隐交替格式等。 目前常见的差分格式,主要是上述几种形式的组合,不同的组合构成不同的差分格式。差分方法主要适用于有结构网格,网格的步长一般根据实际地形的情况和柯朗稳定条件来决定。构造差分的方法有多种形式, 目前主要采用的是泰勒级数展开方法。其基本的差分表达式主要有三种形式:一阶向前差分、一阶向后差分、一阶中心差分和二阶中心差分等,其中前两种格式为一阶计算精度,后两种格式为二阶计算精度。通过对时间和空间这几种不同差分格式的组合,可以组合成不同的差分计算格式。 有限元方法的基础是变分原理和加权余量法,其基本求解思想是把计算域划分为有限个互不重叠的单元,在每个单元内,选择一些合适的节点作为求解函数的插值点,将微分方程中的变量改写成由各变量或其导数的节点值与所选用的插值函数组成的线性表达式,借助于变分原理或加权余量法,将微分方程离散求解。采用不同的权函数和插值函数形式,便构成不同的有限元方法。 有限差分法求解偏微分方程 在采用数值计算方法求解偏微分方程时,若将每一处导数由有限差分近似公式替代,从而把求解偏微分方程的问题转换成求解代数方程的问题,即所谓的有限差分法。有限差分法求解偏微分方程的步骤如下: 1、区域离散化,即把所给偏微分方程的求解区域细分成由有限个格点组成的网格; 2、近似替代,即采用有限差分公式替代每一个格点的导数; 3、逼近求解。换而言之,这一过程可以看作是用一个插值多项式及其微分来代替偏微分方程的解的过程 有限差分法的应用 抛物型方程的差分方法 1. 简单差分法

抛物型方程有限差分法

抛物型方程有限差分法 1. 简单差分法 考虑一维模型热传导方程 (1.1) )(22x f x u a t u +??=??,T t ≤<0 其中a 为常数。)(x f 是给定的连续函数。(1.1)的定解问题分两类: 第一,初值问题(Cauchy 问题):求足够光滑的函数()t x u ,,满足方程(1.1)和初始条件: (1.2) ()()x x u ?=0,, ∞<<∞-x 第二,初边值问题(也称混合问题):求足够光滑的函数()t x u ,,满足方程(1.1)和初始条件: ()13.1 ()()x x u ?=0,, l x l <<- 及边值条件 ()23.1 ()()0,,0==t l u t u , T t ≤≤0 假定()x f 和()x ?在相应的区域光滑,并且于()0,0,()0,l 两点满足相容条件,则上述问题有唯一的充分光滑的解。 现在考虑边值问题(1.1),(1.3)的差分逼近 取 N l h = 为空间步长,M T =τ为时间步长,其中N ,M 是自然数, jh x x j ==, ()N j ,,1,0 =; τk y y k ==, ()M k ,,1,0 = 将矩形域G {}T t l x ≤≤≤≤=0;0分割成矩形网格。其中 ()j i y x ,表示网格节点; h G 表示网格内点(位于开矩形G 中的网格节点)的集合; h G 表示位于闭矩形G 中的网格节点的集合; h Γ表示h G -h G 网格边界点的集合。 k j u 表示定义在网点()k i t x ,处的待求近似解,N j ≤≤0,M k ≤≤0。

注意到在节点()k i t x ,处的微商和差商之间的下列关系((,)k j k j u u x t t t ???? ≡ ?????) : 可得到以下几种最简差分格式 (一) 向前差分格式 ()24.1 ()j j j x u ??==0, k u 0=k N u =0 其中1,,1,0-=N j ,1,,1,0-=M k 。取2h a r τ = 为网比,则进一步有 ()14.1' 1+k j u =k j ru 1++()r 21-k j u +k j ru 1-+j f τ 此差分格式是按层计算:首先,令0=k ,得到 1j u =01+j ru +()r 21-0j u +0 1-j ru +j f τ 于是,利用初值() j j j x u ??==0和边值k u 0=k N u =0,可算出第一层的1 j u , 1,,1,0-=N j 。再由()14.1'取1=k ,可利用1j u 和k u 0=k N u =0算出2j u , 1,,1,0-=N j 。如此下去,即可逐层算出所有k j u (1,,1,0-=N j , 1,,1,0-=M k ) 。 由于第()1+k 层值可以通过第()k 层值直接得到,如此的格式称为显格式。并 视k j u 为()k j t x u ,的近似值。 若记 () T k N k k k u u u 1 21,,,-= u ,()()()()T N x x x 121,,,-=???? ,()()()()T N x f x f x f 121,,,-=τττ f 则显格式()14.1' 可写成向量形式 其中 若记 那末截断误差 (1.5) ()=u R k j () ()[]k j k j h Lu t x u L -,1=()ττO t x t u r k j +??? ? ??????? ??--)~,~(2112122=()2 h O +τ。 其中(,)j k x t 是矩形11+-<

相关文档
最新文档