第七章 行波法和积分变换法

第七章  行波法和积分变换法
第七章  行波法和积分变换法

第七章 行波法和积分变换法

7.2 基础训练

7.2.1 例题分析

一、行波法

例1 求解下列初值问题

()

2

1

000

,0(,)|cos ,

(,)|.

tt xx t t t u a u x t u x t x u x t e -==?-=-∞<<∞>??

==??

解 本题中1()cos ,(),x x x e ?ψ-== 直接应用D ’Alembert 公式(7.2)有

[]1

11(,)cos()cos()22cos()cos x at x at u x t x at x at e d a t

at x e

ξ+--=++-+=+?

例2 试求解有阻尼的波动方程的初值问题

()

220020(1)(,)|(),

(,)|()

(2)

tt xx t t t t v a v v v x v x t x v x t x εε?ψ==-++=-∞<<∞==

解 问题中的泛定方程与无界弦的自由振动相比,多了阻尼项,故不能直接运用D ’Alembert 公式求解。但对于阻力的作用,常常可表示为其解中带一个随时间呈指数衰减的因子,故可令

()(,)(,)

0t

v x t e

u x t ββ-=> (3)

其中β为一待定参数,于是

()2,2,

t

t t t tt tt t xx xx u v e

u v e u u u v e u t ββββββ---???=-=-+= ????

代入阻尼振动方程(1),得

()()

222

220tt xx t u a u u u εβεεββ-+-+-+=

显然,若取βε=,则上式变成标准波动方程

20tt xx u a u -= (4) 将变换(3)代入(2),得 0

(,0)(,0)()v x e u x x ε?-?==

0(,0)(,)|()t

t t d v x e u x t x dt

εψ-?=??=

=??

(,0)(,0)()(,0)()()t u x v x x u x x x ?ψεψ==?

?=+?

(5)

于是求解定解问题(1)--(2)就转化为求解定界问题(4)--(5)。由D ’Alembert 公式(7.2)有

()()()()11(,)22x at

x at u x t x at x at d a ??ψξε?ξξ+-=++-++?????

???? 代入(3)式得远定界问题的解为 ()()()()11(,)22x at

t t

x at

u x t x at x at d e ae εε??ψξε?ξξ+-=

++-++?????????

注意:在求解本例题时我们引入了变量代换(,)(,)t v x t e u x t β-=,虽然是出于物理上的考虑而引进了衰减因子t

e

β-,但由于指数函数的特殊性(即,无论求多少次导数,函数的形式不

变,只是增加了一个常数)而使方程得到了化简,这一点,作为解决问题的方法,对于处理

其他定解问题也提供了一定的启示。

例3 证明球面波问题

()200

,0,0(1)|()

(2)|()(3)

tt xx yy zz t t t u a u u u r t u r u r ?ψ==?=++>>??

=??=?? 的解为

()()()()()1(,)22r at

r at r at r at r at r at u r t d r

ar

??ξψξξ+---+++=

+

? (4) 分析 本例虽然为三维问题,但由初始条件可以看出问题是球对称的,即如果用球坐标表示,未知函数u 只依赖于r ,而与θ?和无关,所以可以化为一维问题来处理。

证明 因为在球对称的情况下有

()22

2

2221u u u ru r r r r r

????=+=??? 令v ru =,此时方程(1)化为 2

,0,0tt rr v a v r t =>> (5)

其通解为

()()12v f r at f r at =++- 即(1)式的通解为

()()

12f r at f r at u r

++-= (6)

将(6)代入初始条件(2)、(3)得

()()

()()()()12''

12(7)

(8)

f r f r r r

af r f r r

r

?ψ+?=???-?=??

将(8)式两边对r 积分得

()()()1201r

f r f r d C a

ξψξξ-=+? (9) 联立(7)、(9)两式,求解得

()()()()()()10201(10)

2221(11)

222r r r r C

f r d a r r C

f r d a ?ξψξξ?ξψξξ?=++???

?=--??

??

将(10)式中的r 换为r at +,(11)式中的r 换为r at -后一并代入通解(6),即得到题目

所要证明的结论。

例4 均匀气体的初始振动区域是一个半径为R 的球,所有气体质点的初始速度为零,而初始浓聚度在求内是常数0u ,在球外是零。求任意位置M 处的浓聚度。

解 本问题的丁解问题为

2000

(),|0,|0tt xx yy zz t t t u a u u u u r R u r R u ==?=++?

=??≥??

?=?

其中r 为扰动区域中心到球外任意一点的距离。显然,此例是球对称的。我们用两种方法求解。方法1,利用上题结论(4);方法2,利用Poisson 公式。

仿照上例,此例的初始条件可以写为

()()0,;

0.0,u r R

r r r R ?ψ

方法1 (1) 若M 在球外,则因为()()(),,0.r R r at R r at r at ?≥+≥+=-所以即而 可能大于、等于或小于R 。故由上题结论(4),此时的解为

()()

(,)2,,20,,r at r at u r t r

r at u r at R r R at r R r r at R at r R at r R ?--=

-?-<-<<+?=??-≥<->+?

即即或

(2)若M 在球内,即r

R r r a t u r a t R

?+<<-+=-<当即时,此时更有

所以, ()0.r at u ?-=于是由上题结论(4),有问题的解为 ()()00

0(,)2r at u r at u u r t u r

-++=

=

()()002,0.,r at R at R r r at r at R r at u ??+≥≥-+=-<-=当即时,而当时,此时,0(0),1R r at R at R r at at R r -<-<<+≥≤-由有而的情况已经讨论过, 故由上题结论(4)有

()()()000(,)22r at u r at r at u u r t r

r

-++?-=

=

()()030,,0r at R r at r at r at R r at ??-≥-=+=+>+=当时,更有所以,此

时由,r at R at R r at R r -≤->+>+有即当时,有

(,)0u r t =。 方法2 用Poisson 公式求解。

当,M

at r R at r R r at R S -<<+-<即时,球面与初始扰动区域有相交部分,此时由定解条件有()0M u ?=,而()0M ψ≡,于是由Poisson 公式有

()()()()()02

'2000222000200sin 1

()44111cos 12221222M

at S M at u u M dS d d a t at a t at u r a t u t t t t art r at u u r at t ar r

πθ?θ?θππθ???

?==????????

?????

?+-=-=-???? ??????

?????---???==-????? 而当r at R ->时,由于球面M

at S 与初始扰动区域无相交部分,即()()

'',M at M M S ?ψ在上

均为零,从而由Poisson 公式有此时的解

(,)0u M t =。 例5 求解下列初值问题

200

()||tt xx yy zz t t t u a u u u u yz u xz ==?=++?

=??=?

解 本题(),(),M yz M xz ?ψ== 故有

()()'''2'

''

2(sin sin )(cos )

sin sin cos sin cos sin (sin cos )(cos )

sin cos cos sin cos cos M y z y r z r yz zr yr r M x z x r z r xz zr xr r ?θ?θθ?θθθ?

ψθ?θθ?θθθ?

==++=+++==++=+++

注意到at r =,从而有

()

'22

sin sin cos sin cos sin sin M at

S M dS yz zr yr r r d d at

π

π

?θ?θθθ?θθ???=+++????

?

?

()

'22

sin cos cos sin cos cos sin M

at

S M dS xz zr xr r r d d at

π

π

ψθ?θθθ?θθ???=+++?????

?

而有三角函数的周期性和正交性有

220

00

sin cos 0,

cos 0,

sin cos 0d d d d π

π

π

π

????θθθθθ====?

??

?

所以

()

'02sin 44M

at

S M dS yzr d yzr yzat at

π

?πθθππ===??

?

()'

4M

at S M dS ayz t at

?π?

=??? ()

'0

2sin 4M

at

S M dS xzr d xzat at

π

ψπθθπ==??

?

故由Poisson 公式有

()()()''1(,)4M M at at S S M M u M t dS dS a t at at tx y z

?ψπ???

??

=+?????=+????

例6 求解下列定解问题

()()()()

22

,,0(,;0),,(,;0),tt xx yy t v a v v c v

x y t v x y x y v x y x y ?ψ?=++-∞<<∞>??==??

其中c 为已知常数。

解 该方程不是标准的波动方程,为使之变成标准的波动方程,以便可以用Poisson 公式求解,受本节例题2的启发,我们令

()(),,;,;c z a

u x y z t e v x y t = 求导后代入原问题得

()()()()2,,0(,,;0),,(,,;0),tt xx yy zz c c

z z

a a t u a u u u x y t u x y z e x y u x y z e x y ?ψ?=++-∞<<∞>?

??==?

于是由Poisson 公式,有

()()

()()'''''',;,,;,,14M M

at at c

z a

c c z z a a c z a

S S v x y t e u x y z t e x y e x y e

dS dS a t at at

?ψπ--=?????=+???????

???? 例7 求解下列定解问题

()()

()()22

,,,0(,;0),

,,(,;0)0tt xx yy t

u a u u x y t u x y x x y x y u x y ?=+-∞<<∞>??=+-∞<<∞??=??

解 由二维Poisson 公式有

()

200

,;cos ,sin 12at u x y t x y d d a t πρθ?

π++?=???

其中ρθ=

为极坐标变量。由初始条件有

()()()

()()()()()()

2

2222

2

32cos ,sin cos cos sin 2cos cos sin cos 2sin cos cos cos sin cos x y x x y x x y x x y x y x x ?ρθρθρθρθρθρθρθρθθρθθθ

ρθθθ++=++++=+++++++++++

由三角函数的周期性、正交性、倍角公式等可得,3

2cos ,sin ,cos sin ,cos ,cos sin θθθθθθθ

在[]0,2π上积分均为零,而 220

cos d π

θθπ=?

故有

()

()()()20

32

233

cos ,sin 232

233

at at

at

x y d d x

x y x y x x y at x y a t π

?ρθρθρθ?

ππππ++=+++=+++??

??

于是

()()()()()

2332221

2,;23233u x y t x x y at x y a t a t x x y a t x y πππ???=

+++?????=+++

二、积分变换

例8 求解弦振动方程的初值问题

???

??∞<<∞-=∞<<∞-=>∞<<∞-=x x u x x x u t x u a u t

xx tt 0)0,()()0,(0,2?

解:视t 为参数,对上述方程和初始条件的两端均进行傅氏变换,并记

)()]([);,()],([ω??ω==x F t u t x u F

则有

???

?

?

????==-=0)0,()()0,(2

222ωω?ωωt u u u a dt u d 这是带参数ω的常微分方程的初值问题,解之得

t a t u ωω?ωcos )(),(=

于是由逆变换公式,得定解问题的解为

)]()([2

1

])[(41cos )(21]cos )([)],([),()

()(11at x at x d e e d te a t a F t u F t x u at x j at x j x j -++=+====??∞∞--+∞∞

---??ωωπωωωπωω?ωωωω 例9 求无界杆的热传导问题

???∞<<∞-=>∞<<∞-+=x x x u t x t x f u a u xx t )

()0,(0

,),(2?

解:对对方程和初始条件两端关于x 分别进行傅氏变换,并记

),()],([t u t x u F ω=

),()],([);()]([t f t x f F x F ωω??==

则有

??

???=+-=)()0,()

,(2

222ω?ωωωu t f u a dt

u d 这是带参数ω关于变量t 的常微分方程的初值问题,解之得

ττωω?ωτωωd e

f e

t u t a t

t

a )

(0

2222),()(),(---?+=

于是,由逆变换公式得问题的解应为

τ

τ?ττωω?ωτω

ωτω

ωd e F F x f F F e F F x F F d e f F e F t u F t x u t

t a t a t a t

t a ??-------------?+?=???

???+==0)

(1111)

(0

111)]}([)],([{)]}([)]([{),(])([)]

,([),(2

2

2

22

2

2

2

故由性质卷积定理有

ττ?τω

ωd e F x f e F x t x u t

t a t a ?-----*+*=0

)

(11)(),()()(),(2

22

2

t

a x t a t a x

j t

a t

a e

t

a xd e d x x e d e

e

e

F 222

22

2222240

1

21cos 1

)sin (cos 2121

)(-∞

-∞

--∞

----=

=+=

=?

?

?

πωωπ

ωωωπ

ωπωωωωω

这里利用了积分公式 )0(21cos 40

2

2

>=-∞

-?

a a

e

bxdx e a b ax

π

于是

??

?

?∞

----∞

-------+

=

-*

+*

=t t a x t

a x t

t a x t a x d d e t f a d e t

a d e t a x f e t

a x t x u 0)

(4)(4)(0

)(4422

22

2

2

),(21)(21)

(21

),(21)(),(τξτ

τξπ

ξξ?πτ

τπτπ?τξξτ

由此例看到,用傅氏变换解方程时不必像分离变量法那样区分齐次方程和非齐次方程,都是按同样的步骤解.

例10 求解半无界弦的振动问题

?

??

??∞<≤==≥==>∞<<=∞

→)

0(0)0,(,0)0,()0(0),(lim ),(),0()0,0(2x x u x u t t x u t f t u t x u a u t x xx tt 解:对方程两边关于变量t 作拉氏变换,并记

?∞

-==0

),()],([),(dt e t x u t x u L p x U pt

2

22

2

)

,()0,()0,(),(dx

p x U d a x u x pu p x U p t =-- 代入初始条件,得

0),(2

2

22=-p x U a p dx U d (1) 再对边界条件关于变量t 作拉氏变换,并记)]([)(t f L p F =,则有

?????==∞

→0),(lim )

(),0(p x U p F p U x (2)

常微分方程(1)的通解为

x a

p x a

p

e

p C e

p C p x U )()(),(21+=-

代入边界条件(2),得

)()(,0)(12p F p C p C ==

故 )(),(p F e p x U a

x p ?=-

而由位移定理有

??

?

?????? ??-=-a x t f L p F e

a

x p

)(

所以

???????

??? ?

?-<

=??????????????? ??-==--a x t a x t f a x t a x t f L L p x U L t x u 0)],([),(11

例11 求解长为l 的均匀细杆的热传导问题

??

?

??<≤=≥==><<=)

0()0,()0(),(,0),0()0,0(012l x u x u t u t l u t u t l x u a u x xx t 解:对方程和边界条件(关于变量t )进行拉氏变换,记),()],([p x U t x u L =,并考虑到初始条件,则得

02

0222=+-a u U a p dx U d 0),0(=p U x

p

u p l U 1

),(=

方程的通解为

x a

p ch p C x a p sh p C p u p x U )()(),(210

++=

由边界条件定出)(),(21p C p C ,便得

l

a

p

ch a x

p ch

p u u p u p x U 010),(-+= 变换公式p

L 1

)1(=

,则 111=???

? ??-p L 及

∑∞=------+=?????

?

??

????14)12(12222)12(cos 12)1(41k t l k a k e l x k k l a p pch x a p ch L πππ

∑∞

=------+

==14)12(112

2

222)12(cos 12)1(4

)],([),(k t

l k a k

e l x

k k u p x U L t x u πππ

例12 在传输线的一端输入电压信号)(t f ,初始条件均为零,求解传输线上电压的变化.

解:这是个半无界问题,定解条件如下:

????

???==+∞

<>=>∞<<=-+++===0

00)()0,0(0)(000t t t x xx tt t u u u x t f u t x u LCu u RC LG RGu 内且在 将方程和边界条件施以关于t 的Laplace 变换,并考虑初始条件,得到

0])([2

22

=-+++dx d u LCp p RC LG RG (1) )(0

p f u x == (1)

方程(1)的通解为

Ax Ax De Be p x u -+=),( (3) 其中

LC RC LG RG LC RC LG p LC RG

p RC LG LCp A 4)(2)(22

2+-

+??

? ??

++=+++= 在实际问题中,一个很重要的情形是LGRC RC LG 4)(2=+,则

RG p LC LC

RC

LG p LC A +=++

=2

其次,∞≠∞

→u x lim ,有0=B ,故

x

RG p LC Ax De De p x u )(

),(+--==

再由边界条件(2),得

x

RG p LC e p p x )(

)(),(+-=

通过反演求),(t x u ,则由延迟定理有

?????>-<--=?=----)

0(0

)0()(])([),(1

x LC t x LC t e

x LC t f e

e p

f L t x u x

RG x

RG px

LC

7.2.2 习题

1 确定下列初值问题的解:

(1) 02=-xx tt u a u ,0)0,(=x u ,1)0,(=x u t

(2) 02=-xx tt u a u ,x x u sin )0,(=,2)0,(x x u t = (3) 02=-xx tt u a u ,3)0,(x x u =,x x u t =)0,( (4) 02=-xx tt u a u ,x x u cos )0,(=,1)0,(-=e x u t

2 求解无界弦的自由振动,设初始位移为)(x ?,初始速度为)('x a ?-。

3 求方程 22u x y x y ?=??

满足边界条件

201|,|cos y x u x u y ====

的解。

4 证明定解问题

2s i n s i n 2c o s s i n s i n 0|()|()

x x x y

y y

y

y x y y x u xu xu xu u x u x ?φ==+--=

==

的解为

sin sin (sin )(sin )1

(,)()2

2x x y

x x y

x x y x x y u x y d ??φξξ-++--+++-=

+

?

5 利用泊松求解下列定解问题

???

??∞<<-∞+=∞<<-∞=++===),,(|),,(0

|)(20

02z y x yz x u z y x u u u u a u t t t zz yy xx tt 6 证明傅里叶变换的卷积定理

1121

2[()()]()*()

F F w F w

f t f t -= 其中

1111221211()[()],()[()]()*()()()f t F F w f t F F w f t f t f f t d ξξξ

--∞

-∞

===-?

7

证明:22

2

214[]x a w t a t

F e -

--=

8 求上半平面内静电场的电位,即求解下列定解问题:

222000

(1)

|()(2)lim 0(3)

y x y u y u f x u =+→∞

??=>??

=??=??

9 用积分变换法解下列定解问题:

222200,0|()|()t x u u

x t t x u x u

x x

?φ==???-∞<<∞>=?????

=????=??? 10 用积分变换法求解下列问题

200

10,0|1|1

x y u

x y x y u y u ==??=>>?????

=+??=???

11 用积分变换法解下列定解问题:

2

22

00

0l 10,0||0,|t x x u u x l t a t x u u u

u u x

===???<<>=?????=????==???

7.2.3 解答与提示

1 解:(1) 这里0)(=x ?,1)(=x ψ,由一维无界波动方程的D ’Alembert 公式,可得

t d a t x u at

x at

x ==?+-ξ21),(

(2) 这里x x sin )(=?,2

)(x x =ψ,由D ’Alembert 公式可得

?+-+-++=at x at x d a

at x at x t x u ξξ2

21)]sin()[sin(21),( )3(3

sin sin 2

22t a x t at x ++=

(3) 这里3

)(x x =?,x x =)(ψ,由D ’Alembert 公式可得

?+-+-++=at

x at x d a

at x at x t x u ξξ21])()[(21),(33 axt xt a x 232

2

3

++=

(4) 这里x x cos )(=?,1

)(-=e x ψ,由D ’Alembert 公式可得

?

+--+

-++=at

x at

x d a

e at x at x t x u ξ2)]cos()[cos(21),(1

e

t at x +=cos cos 2。解:其定解问题为:

??

?-==>∞<<-∞=-)

(')0,(),()0,()

0,(02x a x u x x u t x u a u t xx tt ?? 故由D ’Alembert 公式可得

?+--+-++=at

x at

x d a a at x at x t x u ξξ???)]('[21)]()([21),(

)

()]

()([2

1

)]()([21at x at x at x at x at x -=--+--++=????? 3 解: 方程两边对y 积分,得

22201

()2

y u x ydy x y f x x ?==+?? 其中()f x 任意函数,上式再对x 积分,得

22320

11

()()()26

x

u x y dx f x x y f x g y =

+=++?

其中()g x 任意函数。

由2

0|y u x ==,得 21()f x x c =+ 由1|cos x u y ==,得 2

21()cos 6

g y y y c =-+ 于是问题的解为

322211

cos 66

u x y x y y c =

++-+ 4 证明:

令sin sin x x y

x x y ξη=-+??

=+-?

222222(1cos )(1cos )(1cos )(1cos )(1cos )sin (1cos )(1cos )(1cos )sin x y xx u u u u u x x x x u u u u u

u u x x x x

u u

u x x x x

ξηξηξηξηξξηξξηηη

??????=

+=-++????????=-

???????=--+++????????

?????++-++-????????

2222

2

2

2

2

22(1cos )(1cos )yy xy u u u u u u u u x x ξη

ξξηξηη??=+??????????=--++-??????????????

代入定解问题得到

22

2

2cos sin sin (22cos )0xx xy yy y u u xu xu xu x ξη

?+--=+=??

由于 2

1cos 0x +≠故20u

ξη

?=??

故方程的通解为:

(,)()()(s i n )(s i n u x y f g f x x y g x x y

ξη=+=-+++- 由边界条件

sin sin sin |()()()

(1)

|[()()]|()()()

(2)

y x y y x y x u f x g x x u f g f x g x x ?ξηφ====+=''''=+=-=

对(2)积分得到

()()[()()]()(3)x x

f x

g x f x g x dx x dx

?''-=

-=??

解(1)(3)得到

00

()()()2()()()2

x

x

x x dx

f x x x dx

g x ?φ?φ?+?=

???-?=??? 于是,有

sin 00sin 1

()(sin )(sin )()2

1

()(sin )(sin )()2

x x y x x y f f x x y x x y d g f x x y x x y d ξ?φξξ

η?φξξ

-++-=-+=-++=+-=+-+??

代入通解表达式得

sin sin (sin )(sin )1

(,)()2

2x x y

x x y

x x y x x y u x y d ??φξξ-++--+++-=

+

?

5 解:此处yz x M M +==2

)(,0)(ψ?,故有

)cos )(sin sin ()cos sin (')'(22θ?θ?θψr z r y r x yz x M ++++=+=

?

θθ?θθ

?θ?θsin cos sin cos sin cos cos sin 2cos sin 2

2222r zr yr xr r yz x ++++++=

注意到r at =,从而有

0)

'(=??

M at

S dS at

M ?

?

θθ?θθ?θθ

?θ?θψπ

π

d d r r zr yr xr r yz x dS at

M M at

S sin ]sin cos sin cos sin cos cos sin 2cos sin [)

'(2222220

++++++=?

?

??

有三角函数的周期性和正交性可知

)2

1(4)

'(3

22yzt t a t x a dS at

M M at

S ++

=??

πψ 故由Poisson 公式得

])'()'([41),(????+=

M at

M at S S dS at

M dS at M a t M u ψ?π

yat t a t x ++=322

3

1

6 证明 由定义式并交换积分次序(由于()12f f -∞∞和都是在,上绝对可积的,故积分次序可交换),有

()()()()12121

122i x i x e dx f f x d f d f x e dx ωωξξξξξξπ

π

---∞

-∞

-∞

-∞

-=

-?

?

?

?

令,,x t dx dt ξ-== 上式成为

()()()()()()()

12121211122222i t i i t

f d f t e dt f e d f t e dt F F ξω

ωξωξξπξξ

ππππωω∞∞∞

∞-+---∞-∞-∞

-∞??=????=????

?

得证。

7 证明

(

)2222

22

22

2

2141

1

cos sin 221cos 2a t a t i x a t x a t a t

F e e e d e x i x d e xd ωωωωωωωωωπ

π

ωωπ

-----∞

-∞

-∞

--∞

??=

=

+??=

=

??

?

8 解 ① 记 -i r

-F[u(x,y)]=u(x,y)e dx u(,y)ωω-

=?

-i r

-F[(x)]=

f(x)e

dx f ()ωω-

=?

对上述定解问题中的各项施以x 为变量Fourier 的变换,则得到常微分方程的定解问题:

222

y u u(,y)0

y u(,0)f ()

u(,y)|0d d ωωωωω--

---

→+∞??-=???

=??→????

(4)(5)(6) ② 求解变量y 的常微分方程(4)得

y

-y u(,y)A()e

B()e ωωωωω-

=+ (7)

其中,A()ω和B()ω为含参变量ω的任意常数,由边界条件(5)和(6)所确定。 将(7)式代入(6)式得到

当0ω>时A()0ω=,当0ω<时B()0ω=,所以

-||y

u(,y)c()e

ωωω-

= (8)

其中,c()ω为含参变量ω的任意函数,可由(5)式决定。 将(8)式代入边界条件(5)得c()f ()ωω-

=,所以

-||y u(,y)f ()e

ωωω-

-

=

③ 求u(,y)ω-

的逆变换,即

-1

u(x,y)=F [u(,y)]ω--1

||y -1-1||y =F [f ()e ]F F[f (x)*F (e )]ωωω-

--= (9)

-1||y ||y i x 1F [e ]e e d 2ωωωωπ

---∞

=

?

0(y-ix)(y+ix)0-22

1e d e d 2111y 1

2y ix y+ix x y ωωωωπ

π

π∞-∞??

=

+??????=+=??-+??

??

代入(9)式并应用卷积的定义得原定解问题(1)(3)-的解为

22y

f ()

u x y d (x )y ξξπξ∞

-∞-+?(,)=

9 解 作傅里叶变换

F[u (x ,t )]u (,t )F[(x )]()F[(

ω??ωφψω=== 得

222

2

t 0u (i )u(,t)(,t)t u(,t)()u

|i ()t d d d d ωωωωω?ωωψ

ω=?==-???

=?

??=??

解方程得到 j t -j t 12u(

,t)C e C e ωωω=+ 由初始条件得到

121

2C C (),

C -C ()

?ωψω+== 所以 1211C [()()],C [()()].22

?ωψω?ωψω=

+=- 从而 j t -j t 11u(

,t)[()()]e [()()]e 22

ωωω?ωψω?ωψω=++- 两端作傅里叶积分逆变换得到

11j t -j t 1j t -j t 11

u(x,t)F [u(

,t)]F [()e ()e ]F [()e ()e ]22

ωωωωω?ωψω?ωψω---==++-

[][]00x x 1j t 1-j t x x x t x t

111

(x t)(x t)F F ()d e F F ()d e 22

211

(x t)(x t)()d 22ωω??ψξξψξξ??ψξξ--+-????????=

++-+-????????????????=++-+???

10 解 作Laplace 变换:根据 [()][()](0)L f t pL f t f '=-

01

[

()]()|[1]y u u u L p L L y x x x p =????=-==???? 01|0(

)y u

u pL x

x p =??==??因为

所以 即 2

11u u p x p x p ???=?=?? (1) 对0|1x u y ==+,作Laplace 变换有

0211

(,)|x u x p p p

==

+ (2) 由(1)得到 2()x

u g p p

=

+ 由(2)得到 221()g p p p

=

+

2211(,)x u x p p p p

=

++因此

, 再对(,)u x p 作逆Laplace 变换得到

(,)1u x y xy y =++

11 解 设_

u(x,t)u(x,p)?,则定解问题变为

02

21

x x 0x x l u p u=-(0x l)a a u

u |0,u |p

==?-<

?==??

这方程的通解是

12u u(x,p)C C p =

++ 而满足边界条件的特解则是

0u u(x,p)p =现用展开定理求其原函数,为此注意到上式右端第二项是p 的单值函数,而且它的奇点

p 0=和2

n (2n 1)a p p (n 1,2,3,)2l π-??

==-=????

全为单极点,所以

n

1010

p t p t 0n 1p 0p p (u -u (u -u u(x,t)u Res Res ∞===????

??????

=++??????∑

2

(2n 1)a

t n

2l 100n 1

4(u u )

1(2n 1)x

u (1)cos e 2n 12l

πππ

-??∞

-????

=--=+

--∑

第三章行波法与积分变换法教学提纲

第三章行波法与积分变换法 」 分离变量法,它是求解有限区域内定解问题常用的一种方法。 J 行波法,是一种针对无界域的一维波动方程的求解方法。 」 积分变换法,一个无界域上不受方程类型限制的方法。 作如下代换; X at, X at 利用复合函数求导法则可得 同理可得 2 a 2(£ 代入(1)可得 =0o u(x,t) F( ) G( ) F(X at) G(X at) 这里F,G 为二阶连续可微的函数。再由初始条件可知 F(X ) G(X ) (X ), aF (X ) aG (X ) (X ). X 2 u -2 )(」 2 2」 2 u ~2 先对求积分,再对 求积分,可得u(X,t)d 的一般形式 § 3.1 一维波动方程的达朗贝尔 (D 'alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: 2 u 下 u 2 2 u a 2 , X (X), u 0, (1) (X ),- (2) 2 ■4), (3)

由(3)第二式积分可得 1 X F(x) G(x) - 0 (t)dt C , a 0 利用(3)第一式可得 所以,我们有 1 1 x at u(x,t) [ (x at) (x at)] (t)dt 2 2a x at 此式称为无限弦长自由振动的达朗贝尔公式。 二、 特征方程、特征线及其应用 考虑一般的二阶偏微分方程 AU xx 2BU xy CU yy DU x EU y Fu 0 称下常微分方程为其特征方程 A(dy)2 2Bdxdy C(dx)2 0。 由前面讨论知道,直线x at 常数为波动方程对应特征方程的积分曲线, 称为特征线。已知,左行波F(x at)在特征线x at G 上取值为常数值F(CJ , 右行波G(x at)在特征线x at C 2上取值为常数值G(C 2),且这两个值随着特 征线的移动而变化,实际上,波是沿着特征线方向传播的。称变换( 2)为特征 变换,因此行波法又称特征线法。 注:此方法可以推广的其他类型的问题。 三、 公式的物理意义 由 U(x,t) F (x at) G(x at) 其中F(x at)表示一个沿x 轴负方向传播的行波,G(x at)表示一个沿x 轴正方 向传播的行波。达朗贝尔公式表明:弦上的任意扰动总是以行波形式分别向两个 方向传播出去,其传播速度为a 。因此此法称为行波法。 四、 依赖区间、决定区域、影响区域 F(x) 1 2(X ) 2a (t)dt G(x) (x) 1 x 2a o (t)dt (4)

第三章行波法(2)

补充:(习题2.1 ) 10.求解无限长理想传输线上电压和电流的传播情况,设始电压分布为kx cos A ,初始电流分布为 kx A L C cos 。 解:(1)电压的传播情况: 传输线方程:02=-xx tt v a v ,式中LC a 1 2 =。 初始条件: ?? ? ? ?==--=-======)(sin )sin ()1(1 )(cos 00 0x kx aAk kx Ak L C C j C v x kx A v t x t t ψ? 由达朗伯公式有: ?+-+-++=at x at x d a at x at x t x v ξξψ??)(21)]()([21),( )](cos )(cos [21at x k A at x k A -++=?+-+at x at x d k aA a ξξsin 21 )](cos )(cos [21at x k A at x k A -++=)](cos )(cos [2 at x k at x k A -++-+ )(cos at x k A -= (2)电流的传播情况: 传输线方程:02=-xx tt j a j ,式中LC a 1 2 = ,初始条件: ??? ??? ?==-======)(sin 1 )(cos 000 x kx L Ak v L j x kx A L C j t x t t t ψ? 应用一维无界空间解达朗伯公式: )](cos )(cos [21),(at x k A L C at x k A L C t x j -++=?+-+at x at x d k L Ak a ξξsin 21 )](cos )([cos 2at x k at x k L C A -++= )]cos()(cos [2at x at x k A L LC -++-+ )(cos at x k A L C -= 11.在G/C=R/L 条件下求无限长传输线上的电报方程的通解。 解:关于j 和v 的电报方程为

第三章 行波法与积分变换法

第三章 行波法与积分变换法 分离变量法,它是求解有限区域内定解问题常用的一种方法。 行波法,是一种针对无界域的一维波动方程的求解方法。 积分变换法,一个无界域上不受方程类型限制的方法。 §3.1 一维波动方程的达朗贝尔(D ’alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: .- ),(u ),(u 0, ,- ,0t 02 2 222+∞<<∞==>+∞<<∞??=??==x x x t x x u a t u t t ψ? (1) 作如下代换; ? ? ?-=+=at x at x ηξ, (2) 利用复合函数求导法则可得 22 2 2 2 22 2))((,ηηξξηξηξη ξηηξξ??+???+??=??+????+??=????+??=????+????=??u u u u u x u u u x u x u x u 同理可得 ),2(2222222 2ηηξξ??+???-??=??u u u a t u 代入(1)可得 η ξ???u 2=0。 先对η求积分,再对ξ求积分,可得),(t x u d 的一般形式 )()()()(),(at x G at x F G F t x u -++=+=ηξ 这里G F ,为二阶连续可微的函数。再由初始条件可知 ). ()()(),()()(' ' x x aG x aF x x G x F ψ?=-=+ (3)

由(3)第二式积分可得 C dt t a x G x F x += -?0)(1)()(ψ, 利用(3)第一式可得 . 2 )(21)(21)(, 2 )(21)(21)(00C dt t a x x G C dt t a x x F x x --=++=??ψ?ψ? 所以,我们有 ?+-+-++=at x at x dt t a at x at x t x u )(21)]()([21),(ψ?? (4) 此式称为无限弦长自由振动的达朗贝尔公式。 二、特征方程、特征线及其应用 考虑一般的二阶偏微分方程 02=+++++Fu Eu Du Cu Bu Au y x yy xy xx 称下常微分方程为其特征方程 0)(2)(22=+-dx C Bdxdy dy A 。 由前面讨论知道,直线常数=±at x 为波动方程对应特征方程的积分曲线,称为特征线。已知,左行波)(at x F +在特征线1C at x =+上取值为常数值)(1C F ,右行波)(at x G -在特征线2C at x =-上取值为常数值)(2C G ,且这两个值随着特征线的移动而变化,实际上,波是沿着特征线方向传播的。称变换(2)为特征变换,因此行波法又称特征线法。 注:此方法可以推广的其他类型的问题。 三、公式的物理意义 由 )()(),(at x G at x F t x u -++= 其中)(at x F +表示一个沿x 轴负方向传播的行波, )(at x G -表示一个沿x 轴正方向传播的行波。达朗贝尔公式表明:弦上的任意扰动总是以行波形式分别向两个 方向传播出去,其传播速度为a 。因此此法称为行波法。 四、依赖区间、决定区域、影响区域 由方程的解(4)可以看出,解在(x,t )点的数值由x 轴上区间[x-at,x+at]

第三章-行波法与积分变换法Word版

第三章 行波法与积分变换法 分离变量法,它是求解有限区域内定解问题常用的一种方法。 行波法,是一种针对无界域的一维波动方程的求解方法。 积分变换法,一个无界域上不受方程类型限制的方法。 §3.1 一维波动方程的达朗贝尔(D ’alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: .- ),(u ),(u 0, ,- ,0t 02 2 222+∞<<∞==>+∞<<∞??=??==x x x t x x u a t u t t ψ? (1) 作如下代换; ? ? ?-=+=at x at x ηξ, (2) 利用复合函数求导法则可得 22 2 2 2 22 2))((,ηηξξηξηξη ξηηξξ??+???+??=??+????+??=????+??=????+????=??u u u u u x u u u x u x u x u 同理可得 ),2(2 2222222ηηξξ ??+???-??=??u u u a t u 代入(1)可得 η ξ???u 2=0。 先对η求积分,再对ξ求积分,可得),(t x u d 的一般形式 )()()()(),(at x G at x F G F t x u -++=+=ηξ 这里G F ,为二阶连续可微的函数。再由初始条件可知

). ()()(),()()(' ' x x aG x aF x x G x F ψ?=-=+ (3) 由(3)第二式积分可得 C dt t a x G x F x += -?0)(1)()(ψ, 利用(3)第一式可得 .2 )(21)(21)(,2 )(21)(21)(00C dt t a x x G C dt t a x x F x x --=++=??ψ?ψ? 所以,我们有 ?+-+-++=at x at x dt t a at x at x t x u )(21)]()([21),(ψ?? (4) 此式称为无限弦长自由振动的达朗贝尔公式。 二、特征方程、特征线及其应用 考虑一般的二阶偏微分方程 02=+++++Fu Eu Du Cu Bu Au y x yy xy xx 称下常微分方程为其特征方程 0)(2)(22=+-dx C Bdxdy dy A 。 由前面讨论知道,直线常数=±at x 为波动方程对应特征方程的积分曲线,称为特征线。已知,左行波)(at x F +在特征线1C at x =+上取值为常数值)(1C F ,右行波)(at x G -在特征线2C at x =-上取值为常数值)(2C G ,且这两个值随着特征线的移动而变化,实际上,波是沿着特征线方向传播的。称变换(2)为特征变换,因此行波法又称特征线法。 注:此方法可以推广的其他类型的问题。 三、公式的物理意义 由 )()(),(at x G at x F t x u -++= 其中)(at x F +表示一个沿x 轴负方向传播的行波, )(at x G -表示一个沿x 轴正方向传播的行波。达朗贝尔公式表明:弦上的任意扰动总是以行波形式分别向两个 方向传播出去,其传播速度为a 。因此此法称为行波法。

数学物理方法之行波法与达朗贝尔公式

数学物理方法 泰山医学院 于承斌 cbyu@https://www.360docs.net/doc/cd14790986.html,

第十四章行波法与达朗贝尔公式 14.1 二阶线性偏微分方程的通解 对于给定的偏微分方程,一般不能简单的确定通解, 但对简单的标准形式的方程或一个标准形式进一步化简后, 有的可以得到通解。

例14.1.1 求偏微分方程的通解为:板书讲解P280 例14.1.2 求偏微分方程的通解为:板书讲解P281 14.2 二阶线性偏微分方程的行波解 通解法中有一种特殊的解法――行波法, 即以自变量的线性组合作变量代换,进行求解的一种方法,它对波动方程 类型的求解十分有效.

1.简单的含实系数的二阶线性偏微分方程 为了方便起见,我们首先讨论如下的含实常系数的简单二阶线性偏微分方程 xx xy yy au bu cu ++=(14.2.1) 方程中的系数 ,,a b c 为实常数. ,,a b c (,)x y (说明:这里我们用了小写字母 表示它是实常数,而不是 的函数)

假设方程的行波解具有下列形式 (,)() u x y F y x λ=+代入方程即得 2 ()()()0 a F y x b F y x cF y x λλλλλ′′′′′′+++++=需要求方程的非零解,故2 a b c λλ++=(14.2.2) ''()0 F y x λ+≠上述方程变为

(i) 2 40 b a c ?=?>12(,)()() u x y F y x G y x λλ=+++(14.2.3) 2 40 b a c ?=?=(ii) 122b a λλ==? 对应于抛物型方程,式(14.2.2)有相等的实根11(,)()() u x y F y x xG y x λλ=+++(14.2.4) 对应于双曲型方程,式(14.2.2)有两个不同的实根12 ,λλ

数学物理方法

数学物理方法课程教学大纲 一、课程说明 (一)课程名称:数学物理方法 所属专业:物理、应用物理专业 课程性质:数学、物理学 学分:5 (二)课程简介、目标与任务 这门课主要讲授物理中常用的数学方法,主要内容包括线性空间和线性算符、复变函数、积分变换和δ-函数、数学物理方程和特殊函数等,适当介绍近年来的新发展、新应用。本门课程是物理系学生建立物理直观的数学基础,其中很多内容是为后续物理课程如量子力学、电动力学等服务,是其必需的数学基础。 这门课中的一些数学手段将在今后的基础研究和工程应用中发挥重要的作用,往往构成了相应领域的数学基础。一般来讲,因为同样的方程有同样的解,掌握和运用这些数学方法所体现的物理内容将更深入,更本质。 (三)先修课程要求,与先修课与后续相关课程之间的逻辑关系和内容衔接 本课程以普通物理、高等数学和部分线性代数知识为基础,为后继的基础课程和专业课程研究有关的数学问题作准备,也为今后工作中遇到的数学物理问题求解提供基础。 (四)教材:《数学物理方法》杨孔庆编 参考书:1. 《数学物理方法》柯朗、希尔伯特著 2. 《特殊函数概论》王竹溪、郭敦仁编著 3. 《物理中的数学方法》李政道著 4. 《数学物理方法》梁昆淼编 5. 《数学物理方法》郭敦仁编 6. 《数学物理方法》吴崇试编 二、课程内容与安排 第一部分线性空间及线性算子 第一章R3空间的向量分析 第一节向量的概念 第二节R3空间的向量代数

第三节R3空间的向量分析 第四节R3空间的向量分析的一些重要公式 第二章R3空间曲线坐标系中的向量分析 第一节R3空间中的曲线坐标系 第二节曲线坐标系中的度量 第三节曲线坐标系中标量场梯度的表达式 第四节曲线坐标系中向量场散度的表达式 第五节曲线坐标系中向量场旋度的表达式 第六节曲线坐标系中Laplace(拉普拉斯)算符▽2的表达式第三章线性空间 第一节线性空间的定义 第二节线性空间的内积 第三节Hilbert(希尔伯特)空间 第四节线性算符 第五节线性算符的本征值和本征向量 第二部分复变函数 第四章复变函数的概念 第一节映射 第二节复数 第三节复变函数 第五章解析函数 第一节复变函数的导数 第二节复变函数的解析性 第三节复势 第四节解析函数变换 第六章复变函数积分 第一节复变函数的积分 第二节Cauchy(柯西)积分定理 第三节Cauchy(柯西)积分公式 第四节解析函数高阶导数的积分表达式 第七章复变函数的级数展开

数学物理方法习题解答(完整版)

数学物理方法习题解答 一、复变函数部分习题解答 第一章习题解答 1、证明Re z 在z 平面上处处不可导。 证明:令Re z u iv =+。Re z x =,,0u x v ∴==。 1u x ?=?,0v y ?=?, u v x y ??≠??。 于是u 与v 在z 平面上处处不满足C -R 条件, 所以Re z 在z 平面上处处不可导。 2、试证()2 f z z = 仅在原点有导数。 证明:令()f z u iv =+。()2 2222,0f z z x y u x y v ==+ ∴ =+=。 2,2u u x y x y ??= =??。v v x y ?? ==0 ??。 所以除原点以外,,u v 不满足C -R 条件。而 ,,u u v v x y x y ???? , ????在原点连续,且满足C -R 条件,所以()f z 在原点可微。 ()00 00 00x x y y u v v u f i i x x y y ====???????? '=+=-= ? ?????????。 或:()()()2 * 00 0lim lim lim 0z z x y z f z x i y z ?→?→?=?=?'==?=?-?=?。 2 2 ***0* 00lim lim lim()0z z z z z z z zz z z z z z z z z =?→?→?→+?+?+??==+??→???。 【当0,i z z re θ≠?=,*2i z e z θ-?=?与趋向有关,则上式中**1z z z z ??==??】

3、设333322 ()z 0 ()z=0 0x y i x y f z x y ?+++≠? =+??? ,证明()z f 在原点满足C -R 条件,但不可微。 证明:令()()(),,f z u x y iv x y =+,则 ()332222 22 ,=0 0x y x y u x y x y x y ?-+≠? =+?+??, 332222 22 (,)=0 0x y x y v x y x y x y ?++≠? =+?+?? 。 3 300(,0)(0,0)(0,0)lim lim 1x x x u x u x u x x →→-===, 3300(0,)(0,0)(0,0)lim lim 1y y x u y u y u y y →→--===-; 3300(,0)(0,0)(0,0)lim lim 1x x x v x v x v x x →→-===, 3300(0,)(0,0)(0,0)lim lim 1y y x v y v y v y y →→-===。 (0,0)(0,0),(0,0)(0,0)x y y x u v u v ∴ = =- ()f z ∴ 在原点上满足C -R 条件。 但33332200()(0)() lim lim ()()z z f z f x y i x y z x y x iy →→--++=++。 令y 沿y kx =趋于0,则 333333434322222 0()1(1)1(1) lim ()()(1)(1)(1)z x y i x y k i k k k k i k k k x y x iy k ik k →-++-++-++++-+==+++++ 依赖于k ,()f z ∴在原点不可导。 4、若复变函数()z f 在区域D 上解析并满足下列条件之一,证明其在区域D 上

《高等数学》第四册(数学物理方法)

第一章 复数与复变函数(1) 1.计算 )(1)2; i i i i i -- = -- =-()122(12)(34)(2)5102122. ; 345(34)(34)59165 5 i i i i i i i i i i i i +-++--+++ = + =- =- --+-+5 5 51(3). ; (1)(2)(3) (13)(3) 102i i i i i i i = = = ------ 4 2 2 2 (4).(1)[(1)](2)4; i i i -=-=-=- 1 1 22 ())]a b a b i =+= 1 1 2 2 24s sin )]()(co s sin ); 2 2 i a b i θθθθ=+=++ 3. 设 1z = 2;z i = 试用三角形式表示12z z 及1 2z z 。 解: 121co s sin ;(co s sin ); 4 4 2 6 6 z i z i ππππ=+= + 121155[co s( )sin ( )](co s sin ); 2 4 6 4 6 2 12 12 z z i i π π π π ππ= + ++ = + 12 2[co s( )sin ( )]2(co s sin ); 4 6 4 6 12 12 z i i z ππππππ=- +- =+ 11.设123,,z z z 三点适合条件1230z z z ++=及1231; z z z ===试证明123,,z z z 是一个内接于单位圆 z =1 的正三角形的顶点。 证明:1230;z z ++=z 123231;312;;z z z z z z z z z ∴=--=--=-- 122331;z z z z z z ∴-=-=-123 ,,z z z ∴所组成的三角形为正三角形。 1231z z z === 123 ,,z z z ∴为以z 为圆心,1为半径的圆上的三点。 即123z ,z ,z 是内接于单位圆的正三角形。

数学物理方程第三章行波法与积分变换法word版

第三章 行波法与积分变换法 (第十三讲) 分离变量法,它是求解有限区域内定解问题常用的一种方法。 行波法,是一种针对无界域的一维波动方程的求解方法。 积分变换法,一个无界域上不受方程类型限制的方法。 §3.1 一维波动方程的达朗贝尔(D ’alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: .- ),(u ),(u 0, ,- ,0t 02 2 222+∞<<∞==>+∞<<∞??=??==x x x t x x u a t u t t ψ? (1) 作如下代换; ? ? ?-=+=at x at x ηξ, (2) 利用复合函数求导法则可得 22 2 2 2 22 2))((,ηηξξηξηξη ξηηξξ??+???+??=??+????+??=????+??=????+????=??u u u u u x u u u x u x u x u 同理可得 ),2(2 2222222ηηξξ ??+???-??=??u u u a t u 代入(1)可得 η ξ???u 2=0。 先对η求积分,再对ξ求积分,可得),(t x u d 的一般形式 )()()()(),(at x G at x F G F t x u -++=+=ηξ 这里G F ,为二阶连续可微的函数。再由初始条件可知 ). ()()(),()()(' ' x x aG x aF x x G x F ψ?=-=+ (3)

由(3)第二式积分可得 C dt t a x G x F x +=-?0 )(1)()(ψ, 利用(3)第一式可得 . 2 )(21)(21)(, 2 )(21)(21)(00C dt t a x x G C dt t a x x F x x --=++=??ψ?ψ? 所以,我们有 ?+-+-++=at x at x dt t a at x at x t x u )(21)]()([21),(ψ?? (4) 此式称为无限弦长自由振动的达朗贝尔公式。 例 求解柯西问题: ?????+∞≤≤-∞==+∞≤≤-∞>=-+==.,0,3,,0,03202 x u x u x y u u u y y y yy xy xx 解:其特征方程为 0)(32)(22=--dx dxdy dy 由此可得特征线方程为 d y x c y x =+=-3 因此作变换 ?? ?+=-=y x y x μξ, 3 从而可得 η ξ???u 2=0 从而有 )()3(),(y x G y x F y x u ++-= 由初始条件可得 )()3(3)()3(' ' 2=+-=+x G x F x x G x F 所以有 C x G x F =-)(3)3(,

行波法和达朗贝尔公式

行波法与达朗贝尔公式 我们已经熟悉常微分方程的常规解法:先不考虑任何附加条件,从方程本身求出通解,通解中含有任意常数(积分常数),然后利用附加条件确定这些常数。偏微分方程能否仿照这种办法求解呢? (一)达朗贝尔公式 试研究均匀弦的横振动方程(7-1-6)、均匀杆的纵振动方程(7-1-9)、理想传输线方程(7-1-14),它们具有同一形式 即 (7-4-1) (1)通解 方程(7-4-1)的形式提示我们作代换 (7-4-2) 因为在这个代换下, 方程(7-4-1)就成为 。但为了以后的书写便利,把代换(7-4-2) 修改为 ,0 2 2 2 22 =??? ????-??u x a t . 0 =??? ????-????? ????+?? u x a t x a t , ),(ηξηξ-=+=t a x , x a t x x t t ??+??=????+????=??ξξξ, ??? ????-?? -=?? ??+????=?? x a t x x t t ηηη0 ) /(2 =???u ηξ

即 在此代换下,方程(7-4-1)化为 (7-4-3) 就很容易求解了。 先对 积分,得 (7-4-4) 其中 是任意函数。再对 积分,就得到通解 (7-4-5) 其中 和 都是任意函数。 式(7-4-5)就是偏微分方程(7-4-1)的通解。不同于常微分方程的情况,式中出现任意函数而不是任意常数。 通解(7-4-5)具有鲜明物理意义。以 而论,改用以速度 沿 正方向移动的坐标轴 ,则新旧坐标和时间之间的关系为 而 与时间 T 无关。这就是说,函数的图象在动坐标系中保持不变,亦即是随着动 坐标系以速度 沿 正方向移动的行波。同理,是以速度 沿 负方向移动的行波。 这样,偏微分方程(7-4-1)描写以速度 向两方向传播的行波。 ?????? ? -=+=),(21),(21ηξηξa t x ?? ?-=+=., at x at x ηξ, 0 2 =???η ξu η ) ( ξξ f u =??f ξ ), ()( ) ()()()(21212at x f at x f f f f d f u -++=+=+= ? ηξηξξ1 f 2 f ) (2at x f -a x X ?? ?=-=,,t T at x X ), ()(22X f at x f =-a x ) (1at x f +a x a

数学物理方程第三章行波法与积分变换法

第三章 行波法与积分变换法 (第十三讲 ) 分离变量法,它是求解有限区域内定解问题常用的一种方法。 行波法,是一种针对无界域的一维波动方程的求解方法。 积分变换法,一个无界域上不受方程类型限制的方法。 §3.1 一维波动方程的达朗贝尔(D ’alembert )公式 一、达朗贝尔公式 考察如下Cauchy 问题: .- ),(u ),(u 0, ,- ,0t 02 2 222+∞<<∞==>+∞<<∞??=??==x x x t x x u a t u t t ψ? (1) 作如下代换; ?? ?-=+=at x at x ηξ, (2) 利用复合函数求导法则可得 222222 22))((,ηηξξ ηξηξη ξηηξξ??+???+??=??+????+??=????+??=????+????=??u u u u u x u u u x u x u x u 同理可得 ),2(2222222 2ηηξξ ??+???-??=??u u u a t u 代入(1)可得 η ξ???u 2=0。 先对求积分,再对求积分,可得),(t x u d 的一般形式 )()()()(),(at x G at x F G F t x u -++=+=ηξ 这里G F ,为二阶连续可微的函数。再由初始条件可知 ). ()()(),()()(' ' x x aG x aF x x G x F ψ?=-=+ (3) 由(3)第二式积分可得 C dt t a x G x F x +=-?0 )(1)()(ψ, 利用(3)第一式可得

.2 )(21)(21)(,2 )(21)(21)(00C dt t a x x G C dt t a x x F x x --=++=??ψ?ψ? 所以,我们有 ?+-+-++=at x at x dt t a at x at x t x u )(21)]()([21),(ψ?? (4) 此式称为无限弦长自由振动的达朗贝尔公式。 例 求解柯西问题: ?????+∞≤≤-∞==+∞≤≤-∞>=-+==.,0,3,,0,03202 x u x u x y u u u y y y yy xy xx 解:其特征方程为 0)(32)(22=--dx dxdy dy 由此可得特征线方程为 d y x c y x =+=-3 因此作变换 ?? ?+=-=y x y x μξ, 3 从而可得 η ξ???u 2=0 从而有 )()3(),(y x G y x F y x u ++-= 由初始条件可得 )()3(3)()3(' ' 2=+-=+x G x F x x G x F 所以有 C x G x F =-)(3)3(, 从而可得 C x x G C x x F +=-=4 3)(4 9)3(2 2

积分变换的应用

浅谈积分变换的应用 学院:机械与汽车工程学院 专业:机械工程及自动化 年级:12级 姓名:郑伟锋 学号:201230110266 成绩: 2014年1月

目录 1.积分变换的简介 (3) 1.1积分变换的分类 (3) 1.2傅立叶变换 (3) 1.2拉普拉斯变换 (4) 1.3梅林变换和哈尔克变换 (5) 1.3.1梅林变换 (5) 1.3.2汉克尔变换 (6) 2.各类积分变换的应用 (6) 2.1总述 (6) 2.2傅立叶变换的应用 (6) 2.2.1傅立叶变换在图像处理中的应用 (6) 2.2.2傅立叶变换在信号处理中的应用 (7) 2.3拉普拉斯变换的应用 (8) 2.3.1总述 (8) 2.3.2 运用拉普拉斯变换分析高阶动态电路 (8) 参考文献 (9)

1.积分变换的简介 1.1积分变换的分类 通过参变量积分将一个已知函数变为另一个函数。已知?(x),如果 存在(α、b可为无穷),则称F(s)为?(x)以K(s,x)为核的积分变换。 积分变换无论在数学理论或其应用中都是一种非常有用的工具。最重要的积分变换有傅里叶变换、拉普拉斯变换。由于不同应用的需要,还有其他一些积分变换,其中应用较为广泛的有梅林变换和汉克尔变换,它们都可通过傅里叶变换或拉普拉斯变换转化而来。 1.2傅立叶变换 傅里叶变换是一种分析信号的方法,它可分析信号的成分,也可用这些成分合成信号。许多波形可作为信号的成分,比如正弦波、方波、锯齿波等,傅里叶变换用正弦波作为信号的成分。其定义如下 f(t)是t的周期函数,如果t满足狄里赫莱条件:在一个周期内具有有限个间断点,且在这些间断点上,函数是有限值;在一个周期内具有有限个极值点;绝对可积。则有下图①式成立。称为积分运算f(t)的傅里叶变换, ②式的积分运算叫做F(ω)的傅里叶逆变换。F(ω)叫做f(t)的像函数,f(t)叫做 F(ω)的像原函数。F(ω)是f(t)的像。f(t)是F(ω)原像。 ①傅里叶变换 ②傅里叶逆变换

第三章 行波法与积分变换法

第三章行波法与积分变换法 在第二章中,讨论了分离变量法,它是求解有限区域内定解问题的一个常用方法,只要求解的区域很规则(其边界在某种坐标系中的方程能用若干个只含有一个坐标变量的方程表示),对三种典型的方程均可运用。本章介绍另外两个求解定解问题的方法,一是行波法,一是积分变化法。行波法只能用于求解无界域内波动方程的定解问题,积分变换法不受方程类型的限制,主要用于无界域,但对有界域也能应用。 §3.1 一维波动方程的达朗贝尔(D’Alembert) 要求一个常微分方程的特解,惯用的方法是先求出它的通解,然后利用初始条件确定通解中的任意常数得到特解。对于偏微分方程能否采用类似的方法呢?一般来说是不行的,原因之一是在偏微分方程中很难定义通解的概念,原因之二是即使对某些方程能够定义并求出它的通解,但此通解中包含有任意函数,要由定解条件确定出这些任意函数是会遇到很大困难的。但事情不是绝对得,在少数情况下不仅可以求出偏微分方程的通解(指包含有任意函数的解),而且可以由通解求出特解。本节就一维波动方程来建立它的通解公式,然后由它得到初值问题解的表达式。 对于一维波动方程

22 222 u u a t x ??=?? (3.1) 作如下代换: x at x at ξη=+?? =-? (3.2) 利用复合函数微分法则,得 u u u u u x x x ξηξηξη ???????=+=+??????? 222 2 2 22 ()()2u u u u u x x x u u u ξη ξξηηξηξξηη?????????=+++????????????= ++???? (3.3) 同理有 22222 2 2 222 ()()[ 2]u u u u u a a t u u u a ξξηηξηξξηη ???????=---??????????=-+???? (3.4) 将(3.3)及(3.4)代入(3.1)得 20u ξη ?=?? (3.5) 将(3.5)式对η积分得 ()u f ξξ ?=?,(()f ξ是ξ的任意可微函数) 在对此式对ξ积分得 212(,)()()()()u x t f d f f x at f x at ξξη=+=++-? (3.6) 其中1f ,2f 都是任意二次连续可微函数。(3.6)式就是方程(3.1)得通解(包含两个任意函数的解)。 在各个具体问题中,我们并不满足于求通解,还要确定函数1f ,

第三章 行波法(1)

第三章 行波法 §3.1 达朗贝尔公式(P150-152) 1.确定下列初值问题的解 (1)()()20,,00,,01tt xx t u a u u x u x -=== 解:因为 ()()0,1x x ?ψ== 由达朗贝尔公式有: ()()() ()1 ,2 2x at x at x at x at u x t d a ??ψαα+--++= + ? =t (2)()()220,,0sin ,,0tt xx t u a u u x x u x x -=== 解:因为 ()()2 s i n ,x x x x ?ψ== 由达朗贝尔公式有: ()()() ()1 ,2 2x at x at x at x at u x t d a ??ψαα+--++= + ? =2231 sin cos 626x at x at a t a ?? + +? ? =2231sin cos 3 x at x t a t ++ (3)()()230,,0,,0tt xx t u a u u x x u x x -=== 解:因为 ()()3,x x x x ?ψ== 由达朗贝尔公式有: ()()() ()1 ,2 2x at x at x at x at u x t d a ??ψαα+--++= + ? = ()() 1 cos cos 1 2 2x at x at x at x at e d a α+---+++ ? =1cos cos x at e t -+ 2.求解无界弦的自由振动,设弦的初始位移为()x ?,初始速度为()'a x ?-。 解:该问题的数学模型为:

()()()( ) 2 ' ,,0,0,,0t t x x t u a u x t u x x u x a x ???=-∞<<+∞>?? ==-?? 由达朗贝尔公式: ()()() ()' 1,2 2x at x at x at x at u x t a d a ???αα+--++= + -? =()x at ?- 2.求解弦振动方程的古沙问题 ()()()( )()(),,,,tt xx u u u x x x x u x x x x ?ψ=? ? -=-∞<<+∞??=-∞<<+∞ ? 解:该方程的通解为: ()()()12,u x t f x t f x t =++- (1) 令:t x =- ()()()1202x f f x ?=+ 令: t x = ()()()1220x f x f ψ=+ 令2y x =,则有: ()()()()12210202y f y f y f y f ψ??? ? =- ?????? ? ? ?=- ???? ? 所以: ()()1102x t f x t f ψ+??+=- ???,()()2202x t f x t f ?-?? -=- ??? ()()()12,0022x t x t u x t f f ψ?+-???? =+-+?? ? ??????? 又 ()()121(0)(0)002 f f ?ψ+=+???? 所以古沙问题解为: ()()() 00,22 2x t x t u x t ?ψψ?++-????=++ ? ????? 3.求解无限长理想传输线上电压和电流的传播情况。设初始电压分布为 cos A kx cos kx 。

行波法和达朗贝尔公式

第二章数学物理方程的解 §2.1 行波法 达朗贝尔公式 读者已经熟悉常微分方程的常规解法:先不考虑任何附加条件,从方程本身求出通解,通解中含有任意常数(积分常数),然后利用附加条件确定这些常数。偏微分方程能否仿照这种办法求解呢? (一)达朗贝尔公式 试研究均匀弦的横振动方程(7-1-6)、均匀杆的纵振动方程(7-1-9)、理想传输线方程(7-1-14),它们具有同一形式 即 (7-4-1) (1)通解 方程(7-4-1)的形式提示我们作代换 (7-4-2) 因为在这个代换下, 方程(7-4-1)就成为。 但为了以后的书写便利,把代换(7-4-2) 修改为 ,0 2 2 222=??? ????-??u x a t . 0 =??? ????-????? ????+?? u x a t x a t , ),(ηξηξ-=+=t a x , x a t x x t t ??+??=????+????=??ξξξ, ??? ????-??-=????+????=??x a t x x t t ηηη0) /(2=???u ηξ

即 在此代换下,方程(7-4-1)化为 (7-4-3) 就很容易求解了。 先对积分,得 (7-4-4) 其中是任意函数。再对积分,就得到通解 (7-4-5) 其中 和都是任意函数。 式(7-4-5)就是偏微分方程(7-4-1)的通解。不同于常微分方程的情况,式中出现任意函数而不是任意常数。 通解(7-4-5)具有鲜明物理意义。以 而论,改用以速度沿正方 向移动的坐标轴,则新旧坐标和时间之间的关系为 而 与时间T 无关。这就是说,函数的图象在动坐标系中保持不变,亦即是随着动坐标系以速度沿正方向移动的行波。同理,是以速度沿负方向 移动的行波。 这样,偏微分方程(7-4-1)描写以速度向两方向传播的行波。 ?????? ?-=+=),(21),(21ηξηξa t x ?? ?-=+=., at x at x ηξ,0 2=???ηξu η )( ξξf u =??f ξ), ()( )()()()(21212at x f at x f f f f d f u -++=+=+=?ηξηξξ1f 2f )(2at x f -a x X ?? ?=-=, ,t T at x X ),()(22X f at x f =-a x )(1at x f +a x a

第七章 行波法和积分变换法

第七章 行波法和积分变换法 7.2 基础训练 7.2.1 例题分析 一、行波法 例1 求解下列初值问题 () 2 1 000 ,0(,)|cos , (,)|. tt xx t t t u a u x t u x t x u x t e -==?-=-∞<<∞>?? ==?? 解 本题中1()cos ,(),x x x e ?ψ-== 直接应用D ’Alembert 公式(7.2)有 []1 11(,)cos()cos()22cos()cos x at x at u x t x at x at e d a t at x e ξ+--=++-+=+? 例2 试求解有阻尼的波动方程的初值问题 () 220020(1)(,)|(), (,)|() (2) tt xx t t t t v a v v v x v x t x v x t x εε?ψ==-++=-∞<<∞== 解 问题中的泛定方程与无界弦的自由振动相比,多了阻尼项,故不能直接运用D ’Alembert 公式求解。但对于阻力的作用,常常可表示为其解中带一个随时间呈指数衰减的因子,故可令 ()(,)(,) 0t v x t e u x t ββ-=> (3) 其中β为一待定参数,于是 ()2,2, t t t t tt tt t xx xx u v e u v e u u u v e u t ββββββ---???=-=-+= ???? 代入阻尼振动方程(1),得 ()() 222 220tt xx t u a u u u εβεεββ-+-+-+= 显然,若取βε=,则上式变成标准波动方程 20tt xx u a u -= (4) 将变换(3)代入(2),得 0 (,0)(,0)()v x e u x x ε?-?== 0(,0)(,)|()t t t d v x e u x t x dt εψ-?=??= =??

第三章-行波法

第三章 行波法 §3.1 达朗贝尔法(行波法) 考虑无界弦的自由振动问题,有定解问题如下: ???? ?????==??=??)()0,()()0,(22222x x u x x u x u a t u t ψ? ∞+<∞-+∞ <<∞->+∞<<∞-x x t x 0, 对于上面的标准形方程,它有两族特征曲线 1c at x =+,2c at x =- 作变换 at x +=ξ,at x -=η 由上面的方程变为: 02=???η ξu 求上面偏微分方程的解 先对η积分一次得 )(1ξη f u =?? 再对ξ积分一次得: ?+=+=)()()()(2ηξηξξG F f d f u 其中G F ,是具有任意连续可微函数,将原自变量代回得原方程的通解为 )()(),(at x G at x F t x u -++= 下面通过初始条件确定上面的任意函数G F , ∵ )(0x u t ?==,)(0x u t t ψ== ∴ )()()(x x G x F ?=+ (1)

)()()(//x x aG x aF ψ=- (2) 对(2)从0x 到x 积分得: ?-+= -x x x G x F d a x G x F 0)()()(1)()(00ααψ (3) (1)+(3)得 )]()([2 1)(21)(21)(000x G x F d a x x F x x -++=?ααψ? ?---=x x x G x F d a x x G 0)]()([2 1)(21)(21)(00ααψ? ∴ ?+-+++-=at x at x d a at x at x t x u ααψ??)(21)]()([21),( 该公式叫达朗贝尔公式 例:确定初值问题: ?? ???==>∞+<∞??=??-122222)0,( cos )0,(0 e x u x x u ,t x -x u a t u t 解:略。 达朗贝尔方程的物理定义: 先讨论0)(=x ψ (即振动只有初始位移) )]()([2 1),(at x at x t x u ++-= ?? 先看)(at x -?项: 当0=t 时若观察者位于c x =处,此时 )()(c at x ??=- 在x 轴上,若观察者以速度a 沿轴正方向运动,则在t 时刻观察者位于at c x +=处,此时: )()()(c at at c at x ???=-+=- 由于t 是任意的,这说明观察者在运动过程中随时可以看到相同的波形,可见,波形和

相关文档
最新文档