一.无界问题的特征线法求解求解1.一维无界弦振动方程的达朗贝尔公式(特征线法在弦振动方程的应用)求解法 1.1齐次方程两端无界弦振动方程的求解 齐次弦振动方程及初始条件:⎩⎨⎧==+∞<<-∞>=-)()0,(),()0,(,0,02x x u x x u x t u a u t xx tt ψϕ其方程为+∞<<-∞>=-x t u a u xx tt ,0,02,其特征方程为022=-⎪⎭⎫⎝⎛a dt dx ,2,1c at x =±所以at x +=ξ,at x -=ηηξu u u x +=,ηξu a u a u t ⨯-⨯=,ηηξηξξu u u u xx ++=2,ηηξηξξu a u a u a u tt 2222+-=)()()()(),(0042at x G at x F G F t x u u u u a u xx tt -++=+=⇒=⇒=-=-ηξξηξη由初始条件)()(')(')0,(),()()()0,(x x aG x aF x u x x G x F x u t ψϕ=-==+=来确定⎰=---xx dbb x G x G a x F x F a 0)()]0()([)]0()([ψ)0()0()(1)()(0x G x F db b a x G x F xx -+=-⎰ψ)()()(x x G x F ϕ=+)(2)0()0()(21)(0x x G x F db b a x F xx ϕψ+-+=⎰)(212)0()0()(21)(0at x x G x F db b aat x F at x x ++-+=+⎰+ϕψ)(2)0()0()(21)(0x x G x F db b a x G xx ϕψ+---=⎰ )(2)()()(21)(0at x at x G at x F db b a at x G atx x -+-----=-⎰-ϕψ)()(),(at x G at x F t x u -++=⎰+-+-++=atx atx db b a at x at x t x u )(21)]()([21),(ψϕϕ(1)此公式为达朗贝尔公式 1.2单侧无界弦振动齐次方程的求解⎪⎩⎪⎨⎧>=>==>>=-0,0),0(),()0,(),()0,(0,0,02t t u t t x x u x x u x t u a u t xx tt ψϕ先求出对应双侧无界弦振动方程⎩⎨⎧ψ=Φ=+∞<<-∞>=-)()0,(),()0,(,0,02x x u x x u x t u a u t xx tt 其中要求)(),(x x ψΦ为奇函数又已知其右侧函数表达式可以求出求出左侧表达式⎩⎨⎧<--≥=Φ0),(0),()(x x x x x ϕϕ,⎩⎨⎧<--≥=ψ0),(0),()(x x x x x ψψ 将其带入达朗贝尔公式可求出对应双侧无界弦振动方程的解⎰+-ψ+-Φ++Φ=atx atx db b a at x at x t x u )(21)]()([21),( 只要令0)(21)]()([210),(,0=Φ+Φ-Φ⇒==⎰-db b a at at t x u x atat又令0>x ,⎪⎪⎩⎪⎪⎨⎧<+---+>+-++=⎰⎰+--+-atx at x atx at x at x db b a at x at a a at x db b a at x at x t x u )(,)(21))](()([21,)(21)]()([21),(ϕϕϕϕϕϕ 此),(t x u 即为单侧无界弦振动齐次方程的解 1.3零初始条件的非齐次弦振动方程的求解⎩⎨⎧==>=-0)0,(,0)0,(0),,(2x u x u t t x f u a u t xx tt 设);,(τt x w 为下面齐次方程的解⎩⎨⎧==>=-),(),(,0),(,02ττττx f x u x u t u a u t xx tt 则⎰=td t x w t x u 0);,(),(ττ为零初始条件的非齐次弦振动方程的解(将),(t x f 作用延时效果累积为将齐次化思想)转换计时器的初始时刻将齐次方程初始时刻换为0需要τ-=t t '可得0','>⇒>=t t dt dt τ 齐次方程可以化简为⎩⎨⎧===>=-0'),,()0,(,0)0,(0',0'2''t x f x w x w t w a w t xx t t τ 使用达朗贝尔公式可以求得⎰+-+-++='')(21)]'()'([21)',(at x at x db b a at x at x t x w ψϕϕ其中),()(,0)(τψϕx f x x ==则⎰-+--=)()(),(21),(τττt a x t a x db b f a t x w ⎰⎰⎰++--==t t a x t a x td db b f a d t x w t x u 0)()(0),(21),(),(τττττ 1.4有初始条件的非齐次无界弦波动方程的求解⎩⎨⎧==+∞<<-∞>=-)()0,(),()0,(,0),,(2x x u x x u x t t x f u a u t xx tt ψϕ 此方程要使用叠加原理进行求解设),(),(),(t x z t x v t x u +=其中分别满足以下方程⎩⎨⎧==+∞<<-∞>=-)()0,(),()0,(,0,02x x v x x v x t v a v t xx tt ψϕ(1)和⎩⎨⎧==+∞<<-∞>=-0)0,(,0)0,(,0),,(2x y x y x t t x f y a y t xx tt (2) 对于方程(1),使用达朗贝尔公式可以求得:其特征方程为022=+⎪⎭⎫⎝⎛a dt dx ,2,1c at x =±所以at x +=ξ,at x -=ηηξv v v x +=,ηξv a v a v t ⨯-⨯=,ηηξηξξv v v v xx ++=2,ηηξηξξv a v a v a v tt 2222+-=)()()()(),(0042at x G at x F G F t x v v v v a v xx tt -++=+=⇒=⇒=-=-ηξξηξη由初始条件)()(')(')0,(),()()()0,(x x aG x aF x v x x G x F x v t ψϕ=-==+=来确定⎰=---xx dbb x G x G a x F x F a 0)()]0()([)]0()([ψ)0()0()(1)()(0x G x F db b a x G x F xx -+=-⎰ψ)()()(x x G x F ϕ=+)(2)0()0()(21)(0x x G x F db b a x F xx ϕψ+-+=⎰)(212)0()0()(21)(0at x x G x F db b aat x F at x x ++-+=+⎰+ϕψ)(2)0()0()(21)(0x x G x F db b a x G xx ϕψ+---=⎰)(2)()()(21)(0at x at x G at x F db b a at x G atx x -+-----=-⎰-ϕψ)()(),(at x G at x F t x v -++=⎰+-+-++=atx atx db b a at x at x t x v )(21)]()([21),(ψϕϕ对于方程2,使用齐次化原理可以求得⎩⎨⎧==>=-0)0,(,0)0,(0),,(2x y x y t t x f y a y t xx tt 设);,(τt x w 为下面齐次方程的解⎩⎨⎧==>=-),(),(,0),(,02ττττx f x y x y t y a y t xx tt 则⎰=td t x w t x y 0);,(),(ττ为零初始条件的非齐次弦振动方程的解(将),(t x f 作用延时效果累积为将齐次化思想)转换计时器的初始时刻将齐次方程初始时刻换为0需要τ-=t t '可得0','>⇒>=t t dt dt τ 齐次方程可以化简为⎩⎨⎧===>=-0'),,()0,(,0)0,(0',0'2''t x f x w x w t w a w t xx t t τ 使用达朗贝尔公式可以求得⎰+-+-++='')(21)]'()'([21)',(at x at x db b a at x at x t x w ψϕϕ其中),()(,0)(τψϕx f x x ==则⎰-+--=)()(),(21),(τττt a x t a x db b f a t x w ⎰⎰⎰++--==t t a x t a x td db b f a d t x w t x y 0)()(0),(21),(),(τττττ最后,根据叠加原理求得⎰⎰⎰++--+-++-++=+=t t a x t a x at x at x d db b f a db b a at x at x t x y t x v t x u 0)()(),(21)(21)]()([21),(),(),(ττψϕϕττ1.5.无界弦振动方程的决定区域与影响区域 决定区域:对于特定u(x,t)依赖的(x,t)的取值范围对于(x,t )的取值能影响u(x,t)的取值范围为影响区域2.只含二阶导的2阶偏微分方程的特征线法求解 2.1只含二阶导的二阶偏微分方程的初步化简⎩⎨⎧===++)(),0(),(),0(0y y u y y u Cu Bu Au x yy xy xx ψϕ其特征方程为00,0222=+-⎪⎭⎫ ⎝⎛⇒-=⇒=+==++C dx dy B dx dy A dx dy dy dx d C B A y x y x y y x x ϕϕϕϕϕϕϕϕϕ根据特征方程解的三种不同情况将其进行进一步的化简 2.2特征方程存在两个不同实根时的化简 先用公式法求出特征方程两个不同的实根A ACB B dx dy 242-±=,g A AC B B dx dy =-+=⎪⎭⎫ ⎝⎛2421,e A AC B B dx dy =--=⎪⎭⎫⎝⎛24221c gx y +=2c ex y +=可以用换元法对此偏微分方程进行化简x A AC B B y 242-+-=ξxAACB B y 242---=η将其带入=++yy xy xx Cu Bu Au=ξηu例1.化简下列方程并求解⎩⎨⎧===-+σφ)0,(,)0,(032t u t u u u u x xx tx tt3/2)/(032032222=-+⇒=-+⇒=-+x t x t x x t t xx tx tt u u u ϕϕϕϕϕϕϕϕdtdx dx dt d x t x t //0-=⇒=+=ϕϕϕϕϕ03/2)/(03)/(2)/(22=--⇒=--+dt dx dt dx dt dx dt dx,0,0,3,10,0,0,1,13)2(,)2(22121242===-=======-=+-=+=--=+±=⇒±=+±=tt xt xx t x tt tx xx t x tx t t x t x t t x c t t x dt dx ηηηηηξξξξξηξηηξηξξηξηηηξξηξξηηξηξξηξηηηξξηξξηηξηξξηξηηξηξηξξηξηξηξηξηξηξηξηξηξηξηξu u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u u xt xt x x tx xx xx x x xx tt tt tt tt x x x t t t 32)3()3(2)()(96)3(3)3(1,3--=++-+-=++=+++++=+-=++---=+=+=-=+=)()(),(00)369()646()321(32ηξξηηηξηξξg f t x u u u u u u u u xx tx tt +==⇒=--+---+-+=-+2.3当特征方程存在2个相等实根A B dx dy 2)(2,1=12c x AB y =-),0(,2≠=-=B y x A By ηξ 0,0·,0,00====⇒=xx yy u C u A B 或如例1化简下列方程44=++xx tx tt u u u4/4)/(044044222=++⇒=++⇒=++x t x t x x t t xx tx tt u u u ϕϕϕϕϕϕϕϕdtdx dx dt d x t x t //0-=⇒=+=ϕϕϕϕϕ2/,04/4)/(04)/(4)/(22==+-⇒=+-+dt dx dt dx dt dx dt dx dt dx,0,10,2,1,,2========-===-=xt xx tt t x tt xt xx t x x t x ηηηηηξξξξξηξηηξηξξηξηηξηξξηηξηξξηξηηξηξξξξηξηηξηξξηξηηηξηξξξηξηηξξηξηηξξu u u u u u u u u u u u u u u u u u u u u u u u u tx tx x t t x x t x t tx xx xx x x x x xx tt tt t t t t tt 222)(22422222---=+++++=++=++++==++++=0)480()880()4244(=⇒=+-++-+⨯-+ηηηηξηξξu u u u)2()2()()()(t x g t x xf g f u f u -+-=+=⇒=ξξηξη2.4当特征方程存在一对共轭复根时二.积分变换法求解无界一维波动方程、1维热传导方程和二维Laplace 方程 1.傅立叶变换的定义与性质 1.1傅立叶变换的定义)()())((w F dx e x f x f F iwx ==⎰+∞∞-1.2傅立叶变换的位移性质)()()()]([)(c x d ee c xf dx e c x f c x f F iwcRRc x iw iwx --=-=-----⎰⎰)()]([)()()]([)(w F e x f F e c x d e c x f e c x f F iwc Riwc c x iw iwc -----==--=-⎰1.3.傅立叶变换的相似性质dcx e cx f c dcx c ecx f dx ecx f cx f F Rcx c wi Rcx cw i Riwx⎰⎰⎰---===)(11)()()]([)(1)(1)]([1c wF c du e u f c cx f F u c wR ==-⎰1.3傅立叶变换的微分性质⎰⎰⎰-+∞∞-----===RiwxRiwx iwx Riwx dex f e x f x df e dx e x f x f F )(|)()()('))('( )())(()())((0))('(w iwF x f iwF dx e x f iw dx e iw x f x f F Riwx iwx R===--=⎰⎰--⎰⎰⎰-+∞∞-----===Riwx iwx Riwx Riwx dex f e x f x df e dx e x f x f F )('|)(')(')(''))(''( )()())(()())('())(''(22w F iw x f F iw x f iwF x f F ===dx e x f iw e x f x df e dx e x f x f F iwx Rn iwx n n Riwx Riwx n n -------⎰⎰⎰+===)()()()())(()1()1()1()()()()())(()())(())((1)(w F iw x f F iw x f iwF x f F n n n n ===-1.3.傅立叶变换的乘多项式性质⎰⎰⎰---=-==R Riwx iwx iwx Rdx e x f dw di dx e x f dw d i dx e x xf x xf F ))(())((1)())(( ))(())((())(())((w F dwdi x f F dw d i dx e x f dw d ix xf F R iwx ===⎰- ⎰⎰⎰---===R Riwx iwx Riwxdx e x f dw d i dx e x xf dw d i dx ex xxf x f x F ))(())(()())((2222)())(())(())((2222222222w F dw d i dx e x f dw d i dx e x f dw d i x f x F R iwx iwx R===⎰⎰-- dx e x f x dwd idx e x f xx dx e x f x x f x F iwx n RRiwx n Riwx n n ))(()()())((11-----⎰⎰⎰=== ⎰⎰====--Rn nn n n n R iwx n n n iwx n n nnw F dw d i x f F dw d i dx e x f dw d i dx e x f dw d i x f x F ))(()))((())(())(())((1.4傅立叶变换积分性质由傅立叶变换的微分性质)())((x f dt t f dx dx=⎰∞- ⎰∞-=xdt t f iw x f F )())(()(1))((1))((w F iwx f F iw dt t f F x==⎰∞- 1.5傅立叶变换的卷积性质卷积定义式⎰-=*Rdt t x g t f x g f )()()(卷积公式1)()()(w G w F g f F =*先做卷积再变换系数不变 证明:⎰⎰⎰⎰-----=-=*R iwt t x iw Riwx R Rdx e e dt t x g t f dx dte t x g t f x g f F )()()()()())((⎰⎰⎰⎰-----=-=*RRiwu iwt Rt x iw Riwt dt du e u g e t f dt dx e t x g e t f x g f F )()()()())(()()()())(())(())(()()(w G w F t f F u g F dt u g F e t f g f F Riwt ===*⎰-卷积公式2))()((2)()(x g x f F w G w F π=*先傅立叶变换再做卷积系数要乘系数2π 1.6 主要函数的傅立叶变换)(0,00,)(指数信号⎩⎨⎧<>=-x x e x f x β iw e iw dx e dx eex f F iw x iw x iwxx +=+-===∞++-+∞+-+∞--⎰⎰βββββ1|1))((0)(0)(02)(x ex f -=2.傅立叶变换法求解一维波动方程 2.1无界齐次波动方程的求解⎪⎩⎪⎨⎧==>∈=-)3)(()0,()2)(()0,()1(0,,02x x u x x u t R x u a u txx tt ψϕ 分别对(1)、(2)、(3)式进行傅立叶变换)4(0),()()),((0),()()),((22=+⇒=-t w F aw t w u F t w F iaw t w u F tt tt)5))((())0,((x F w u F ϕ=)6))((())0,((x F w u F t ψ=)7()()()),((21iawt iawt e w C e w C t w u F -+=将(5)、(6)代入(7)式⎩⎨⎧-=+=--iawtawt t iawtiawt e awiC e w awiC t w u F e w C e w C t w u F 2121)()),(()()()),(( ⎩⎨⎧=-=+))(()()())(()()(2121x F w awiC w awiC x F w C w C ψϕ ⎪⎩⎪⎨⎧-=+=)))((1))(((21)()))((1))(((21)(21x F iaw x F w C x F iaw x F w C ψϕψϕ iawt iawt e x F iawx F e x F iaw x F t w u F --++=)))((1))(((21)))((1))(((21)),((ψϕψϕ又由傅立叶变换的位移性质))(()())((x f F e dx e c x f c x f F iwc Riwx --=-=-⎰左边的项的位移系数可以求出at c iwat iwc -=⇒=-)8))(((21))((21at x F e x F iawt +=ϕϕ iwaw F w G at x G e w G e w G F e x F iwaiawt iawt iawt 2))(()()()())(())((21ψψ=+===用傅立叶变换的积分性质进一步化简))((1))(()())((x f F iw dy y f F x f dy x f dx d xx =⇒=⎰⎰∞-∞- ))((21))((1212))(()()(⎰+∞-===+=atx dy y F a w F iw a iwa w F at x G w G ψψψ右边第一项的系数也可以用位移性质求出at c iwat iwc =⇒-=-))((21))((21at x F e x F iwt -=-ϕϕ iwaw F w H at x H e w H e x F iwaiwat iwat 2))(()()()())((21ψψ=-==--继续用傅立叶变换积分性质来化简))((1))(()())((x f F iw dy y f F x f dy x f dx d xx =⇒=⎰⎰∞-∞-))((21))((1212))(()()(⎰-∞-===-=atx dy y F a w F iw a iwa w F at x H w H ψψψ 四项全部求和 )))((21))(((21)))((21))(((21)),((⎰⎰-∞-+∞---+++=atx at x dy y F a at x F dy y F a at x F t w u F ψϕψϕ ))((21))(()(((21)),((⎰+-+-++=atx atx dy y F a at x F at x F t w u F ψϕϕ 对此式施加傅立叶逆变换 ⎰+-+-++=at a at x dy y a at x at x t x u )(21))()((21),(ψϕϕ 2.2非齐次方程的无界波动方程(不用齐次化原理)2.3半无界波动方程的求解3.傅立叶变换法求解一维热传导方程4.傅立叶变换法求解2维Laplace 方程place 变换的定义与性质place 变换求解一维波动方程place 变换求解一维热传导方程place 变换求解2维Laplace 方程二.有限边界的分离变量法求解(正弦初始条件以及二次初始条件)1.第一类边界条件和第二类边界条件第三类边界条件的特征值问题2.齐次化方程(可以用傅里叶级数展开或用齐次化原理)3.齐次化边界条件4.齐次方程,齐次边界条件第一类边界条件和第二类边界条件的波动方程和热传导方程推导与例子5.齐次方程,非齐次边界条件第一类边界条件和第二类边界条件的波动方程和热传导方程推导与例子6.非齐次方程,非齐次边界条件第一类边界条件和第二类边界条件的波动方程和热传导方程推导与例子7.非齐次方程,非齐次边界条件第一类边界条件和第二类边界条件的波动方程和热传导方程推导与例子8.圆域LAPLACE 问题求解9.矩形域Laplace 方程。