数学物理方程复习答案
《数学物理方程》习题参考答案(A)
《数学物理方程》习题参考答案(A)习题一1.判断方程的类型,并将其化成标准形式:0212222=∂∂+∂∂+∂∂y uyu y x u . 解:⎪⎩⎪⎨⎧==><<>-=-≡∆.0,0. ,00,.0,02211212时,抛物型当椭圆型时当时,双曲型当y y y y a a a①当0<y 时,所给方程为双曲型,其特征方程为,0)()(22=+dx y dy 即 ,0])([)(22=--dx y dy就是 0))((=---+dx y dy dx y dy .积分之,得 c y x =-±2,此即两族相异的实特征线.作可逆自变量代换⎪⎩⎪⎨⎧--=-+=,2,2y x y x ηξ则.1 ,1 ,1 ,1yy yy x x -=∂∂--=∂∂=∂∂=∂∂ηξηξ,2 ,2222222ηηξξηξηηξξ∂∂+∂∂∂+∂∂=∂∂∂∂+∂∂=∂∂∂∂+∂∂∂∂=∂∂u u u x u u u y u x u x u ),(1ηξ∂∂+∂∂--=∂∂u u yyu ).1)(2()(121 ]1)1( 1)1([1)()(12122222222222322y u u u u u y y yu yu yuy u y u u y y u -∂∂+∂∂∂-∂∂+∂∂+∂∂---=-∂∂+--∂∂∂++-∂∂∂---∂∂--+∂∂+∂∂--=∂∂ηηξξηξηξηηξξηξ将这些偏导数代入原方程,得附注:若令⎩⎨⎧=-⇒-==0 ,2,ηηξξηξu u y x 碰巧(双曲型的另一标准形),这是巧合.②当0>y 时,所给方程为椭圆型,其特征方程为0)()(22=+dx y dy即 .0))((=-+dx y i dy dx y i dy 其特征线为 )2 ( 2c ix y c y i x =±=±或.作可逆自变量代换 ⎩⎨⎧==,2,y x ηξ则, 1 , 0 , 0 ,1y y y x x =∂∂=∂∂=∂∂=∂∂ηξηξ, 1 , ηξ∂∂=∂∂∂∂=∂∂u y y u u x u . 1121 , 22222222ηηξ∂∂+∂∂-=∂∂∂∂=∂∂u y u y y yu u x u 将这些偏导数代入原方程,得, 021212222=∂∂+∂∂+∂∂-∂∂ηηηξuy u u y u , 0 2222=∂∂+∂∂∴ηξu u 此即(0>y 时)所求之标准形. ③0=y 时,原方程变为 , 02122=∂∂+∂∂y uxu 已是标准形了(不必再化).2.化标准形:. 0222222222222=∂∂∂+∂∂∂+∂∂∂+∂∂∂+∂∂+∂∂t z ut x u z x u y x u zu x u解: u Lu )2222(434131212321δδδδδδδδδδ+++++≡.这是 ⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂∂∂=⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=t z y x4321δδδδδ 的二次型,于是 , u A Lu Tδδ=其中 010*********1111⎪⎪⎪⎪⎪⎭⎫⎝⎛=A 为实对称矩阵.则∃可逆矩阵M ,使 TMAM B = 为对角形. 令 , 'δδT M = 其中 , '4'3'2'1'''''⎪⎪⎪⎪⎪⎭⎫ ⎝⎛=⎪⎪⎪⎪⎪⎪⎪⎪⎪⎭⎫⎝⎛∂∂∂∂∂∂∂∂=δδδδδt z y x 则 u B u MAM Lu T T T '''')()(δδδδ==.M 的找法很多,可配方,可从矩阵入手等.取 ,11000110001100011-=⎪⎪⎪⎪⎪⎭⎫⎝⎛---=N M , 1000110011101111)(1⎪⎪⎪⎪⎪⎭⎫⎝⎛==-TT M N . , 1''''''⎪⎪⎪⎪⎪⎭⎫⎝⎛===⎪⎪⎪⎪⎪⎭⎫ ⎝⎛==-t zy x M MX X N t z y x X N T δδ则.)( )( 2222'2'2'2'2'''tu z uy u x u u B uMAM u A Lu TT T T ∂∂-∂∂+∂∂-∂∂====δδδδδδ这是超双曲型方程的标准形式.习题二1.决定任意函数法:(1).求解第一问题(0))(0) ( ).(),( , 002ψϕψϕ=⎪⎩⎪⎨⎧======-x ux u u a u at x at x xx tt .解:所给方程为双曲型,其特征线为 c at x =±. 令⎩⎨⎧-=+=,,at x at x ηξ 则可将方程化为 0=ξηu .其一般解为)()(),(21at x f at x f t x u -++= (其中21,f f 为二次连续可微函数). 由定解条件有)0()0()0()0( ).()2()0(),()0()2(212121ψϕψϕ==+⇒⎩⎨⎧=+=+f f x x f f x f x f . 则 ⎪⎩⎪⎨⎧-=-=⇒⎩⎨⎧-=-=).0()2()(),0()2()( ),0()()2(),0()()2(12211221f Y Y f f X X f f x x f f x x f ψϕψϕ 故 )()(),(21at x f at x f t x u -++=).0()2()2()]0()0([)2()2(21ϕψϕψϕ--++=+--++=at x at x f f atx at x (2).求解第二问题 ))0()0( ( ).(),( ,101002ϕϕϕϕ=⎪⎩⎪⎨⎧=====x u x u u a u t at x xx tt解:泛定方程的一般解为)()(),(21at x f at x f t x u -++=由定解条件有 (0))(0)(0)( ).()()(),()0()2(021121021ϕϕϕ=+⎩⎨⎧=+=+f f x x f x f x f x f 则 ),0()2()(201f xx f -=ϕ).0()2()()()()(201112f x x x f x x f +-=-=ϕϕϕ故 )()(),(21at x f at x f t x u -++= ).()2()2(100at x atx at x -+--+=ϕϕϕ (3).证明方程22222)1(])1[(tu h x a x u h x x ∂∂-=∂∂-∂∂ 的解可以写成)]()([1),(21at x f at x f xh t x u -++-=. 由此求该方程满足Cauchy 条件 ⎩⎨⎧====)(),(00x u x u t t t ψϕ 的解.解:令 ),,()(),(t x u x h t x v -= 则 ),(t x v 满足方程 xx tt v a v 2=.)()(),( 21at x f at x f t x v -++=∴.故 )]()([1),(21at x f at x f xh t x u -++-=. 因),(t x v 满足 ⎪⎩⎪⎨⎧≡-=≡-====),()()(),()()( ,10002x x x h vx x x h v v a v t t t xx tt ψϕϕϕ由D'Alembert 公式,得⎰+-+-++=atx atx d a at x at x t x v ααψϕϕ)(21)]()([21),( )]())(()())([(2100at x at x h at x at x h ---+++-=ϕϕ+ααϕαd h a atx at x ⎰+--)()(211 故 ),(1),(t x v xh t x u -=[]⎭⎬⎫⎩⎨⎧-+---+++--=⎰+-atx atx d h a at x at x h at x at x h x h ααϕαϕϕ)()(21)())(()())((211100 即为所求之解.2.Poisson 公式及应用:(1).若),,,(t z y x u u =是初值问题 ⎪⎩⎪⎨⎧+=+=>++===)()( , )()(),0( )(002z y uy g x f u t u u u a u t t t zz yy xx tt ψϕ的解,试求解的表达式.解:IIIIIIu u u u ++=(线性叠加原理),其中IIIIII,,u u u 分别满足如下的初值问题:.0 ),(),0( )(:002I ⎪⎩⎪⎨⎧==>++===t t t zz yy xx tt ux f u t u u u a u u).( ),(),0( )(:002II ⎪⎩⎪⎨⎧==>++===y uy g u t u u u a u u t t t zz yy xx tt ϕ).( ,0),0( )(:002III ⎪⎩⎪⎨⎧==>++===z uu t u u u a u u t t t zz yy xx tt ψ由Poisson 公式,可得⎰⎰∂∂=MatS dS f t a t u ])( 41[2I ξπ)].()([21])(21[at x f at x f d f a t atx atx -++=∂∂=⎰+-ξξ.)(21)( 41.)(21)]()([21 ])( 41[)( 412III22II ⎰⎰⎰⎰⎰⎰⎰⎰+-+-==+-++=∂∂+=Mat M atMat S atz at z aty aty S S d a d t a ud aat y g at y g dS g t a t dS t a u ζζψζζψπηηϕηπηϕπ故IIIII I ),,,(u u u t z y x u ++=.)(21)(2a1)]()([21)]()([21 ⎰⎰+-+-++-+++-++=atz at z aty aty d a d at y g at y g at x f at x f ζζψηηϕ(2).求解初值问题 ⎪⎩⎪⎨⎧+==>-+++=== . ,00),(t )(2)(2002yz x u u z y u u u a u t t t zz yy xx tt解: IIIu u u +=,其中I u : ⎪⎩⎪⎨⎧+==>++=== . ,00),(t )(2002yz x u u u u u a u t t t zz yy xx ttII u : ⎪⎩⎪⎨⎧==>-+++===.0 ,00),(t )(2)(002t t t zz yy xx tt uu z y u u u a u由poisson 公式,得32222I 31)()( 41t a t yz x dS t a u Mat S ++=+=⎰⎰ηζξπ. 由Duhamel 原理,得.)( ])(2)( 41[);,,,(2020II)(t z y d dS t a d t z y x w u M t a S tt-=--==⎰⎰⎰⎰-τζητπτττ故 2322)(31)(),,,(t z y t a t yz x t z y x u -+++= 即为所求. 3.降维法:⎪⎩⎪⎨⎧==>++===.0 ,00),(t ),,()(002t t t yy xx tt uu t y x f u u a u 解:把所给初值问题的解),,(t y x u 看作),,,(t z y x 空间中的函数,即与y x ,平面垂直的直线上的函数值都相等:),,(),,,(*t y x u t z y x u =,则 ),,,(*t z y x u 应形式的满足⎪⎩⎪⎨⎧==>+++=== .0 ,00),(t ),,()(0*0****2*t t t zz yy xx tt u u t y x f u u u a u 由推迟势可得dV ra rt f a t z y x u atr ⎰⎰⎰≤-=),,( 41),,,(2*ηξπττηξτπτττηξπττd dS f t a d dS t f a tS tS M t a M t a ]),,([141]),,([ 410202)()(⎰⎰⎰⎰⎰⎰---=-=τηξτηξττηξτπτd y x t a d d t a f t a ty x M t a ])()()( )(),,(2[141222202),()9------∑-=⎰⎰⎰-τηξτηξτηξπτd y x t a d d f a tx M t a ])()()( ),,([ 212222),()(⎰⎰⎰∑-----=-.此即所求初值问题解的积分表达式.习题三1.求解特征值问题 ⎩⎨⎧=+=<<=+ . 0)()( ,0)0(),(0 0)()("''l X l X X l x x X x X λ 解:该特征值问题要有解0≥⇔λ.0>λ时,记2ωλ=,则 x B x A x X ωωsin cos )(+=.x B x A x X ωωωωcos sin )('+-=. 1(*) 由 0)0('=X ,有 0=B .从而 x A x X A ωcos )(,0=≠. 由 0sin cos ,0)()('=-=+l A l A l X l X ωωω有. ωω=l cot . 此即确定 ω(从而确定λ)的超越方程.由图解法,曲线 ωω==y l y cot 和 有无穷个交点,其横坐标<<<<<n ωωω210,从而 ),2,1( 2==n nn ωλ 便是非0特征值,相应的特征函数为2(*) ,2,1 , cos )( ==n x A x X n n n ω.)( , )( 0'A x XB Ax x X =+==时,λ由0)0('=X ,有0=A .由0)()('=+l X l X , 有 0=B .此时只有平凡解 0)(≡x X . 综上,所求特征值问题的解),2,1( , cos )( ==n x A x X n n n ω.其中n ω为超越方程 ωω=l cot 的正根.附注:下证特征函数系{}∞=1cos n n x ω是],0[l 上的正交系:事实上,设x x X n n ωcos )(=和x x X m m ωcos )(=分别是相应于不同特征值2n n ωλ=和2m m ωλ=的特征函数,即)(x X n 和)(x X m 分别满足).()(,0)0(,0)()(:)(''"⎩⎨⎧+==+l X l X X x X x X x X n n nn n n n λ (1) ⎩⎨⎧=+==+.0)()(,0)0(,0)()(:)(''"l X l X X x X x X x X m m m m m m m λ (2) 则[]0 )()2()()1(0=⋅-⋅⎰dx x X x Xln m,即 []⎰-+-=lm n m n n m m n dx x X x X x X x X x X x X"" )()()())()()()((0λλdx x X x X lm n m n ⎰-=0)()()(λλ若,m n λλ≠则 ),2,1,( 0)()(0==⎰m n dx x X x X lm n .即在],0[l 上,不同特征值所对应的特征函数彼此正交. 2.用分离变量法求波动方程混合问题⎪⎩⎪⎨⎧≤≤==>==><<+=== ),0( , ),0( ),( ,),0(),0 ,0( 20022l x x ux u t t t l u t t u t l x g u a u t t t x xx tt的形式解,其中g 为常数.解:(1).边界条件齐次化:令 ),,(),(),(t x Q t x v t x u +=使⎪⎩⎪⎨⎧====,,20t Q t Q l x x x (这不是定解问题),则取 2)(),(t t l x t x Q +-=即可. 这时),(t x v 满足 ⎪⎩⎪⎨⎧≤≤--==>==><<-+===).0( )( , 0),( 0),( ,0),0(),0 ,0( 2200t 2l x l x x vx v t t l v t v t l x g v a v t t x xx tt(2).“拆”——由线性叠加原理:IIIv v v +=,其中⎪⎩⎪⎨⎧+-====><<=== ., ,0),(),0(),0,0( :2002I l x x vx v t l v t v t l x v a v v t t t x xx tt ⎪⎩⎪⎨⎧====><<-+=== .0,0 ,0),(),0(),0,0( 2:002IIt t t x xx tt vv t l v t v t l x g v a v v (3).用分离变量法求得l x n l at n b l at n a t x v n n n 2 )12(cos 2 )12(sin 2 )12(cos ),(1Iπππ-⎥⎦⎤⎢⎣⎡-+-=∑∞=. 其中⎰⎰--=ll n d ln d ln a 022)12(cos2)12(cos 1ξπξξξπξ,ξπξξξξπξπd ln l d l n l a n b lln 2)12(cos )(2)12(cos 2 )12(122-+---=⎰⎰..,2,1 =n (n n b a ,都可算出来).(4).由Duhamel 原理: ττd t x w t x v t⎰=0II),,(),(,其中),,(τt x w 满足 ⎪⎩⎪⎨⎧-====><<=== . 2 , 0 ,0),( ,0),0( ),,(0 2g ww t l w t w t l x w a w t t t x xx tt τττ用分离变量法求得∑∞=---=12 )12(cos 2)( )12(sin),,(n n l xn l t a n c t x w πτπτ.其中 ξπξξπξπd ln g d l n l a n c lln 2)12(cos)2(2)12(cos 2 )12(12----=⎰⎰. ,3,2,1 =n (n c 可算出).综上: ),(),(),(),(),(),(III t x Q t x v t x v t x Q t x v t x u ++=+=.习题四1.用分离变量法求热方程混合问题⎪⎩⎪⎨⎧===><<-== )( ,0),(),0(),0,0( 022x u t l u t u t l x u b u a u t xx t ϕ 的形式解.解:这是齐次方程、齐次边界条件情形,直接分离变量: 令 )()(),(t T x X t x u =,代入泛定方程,得),( )(22'"λ-=+=a bTa T X X 从而 0)()()( , 0)()(2'"=++=+t T b a t T x X x X λλ. 由边界条件,得 ,0)()0(==l X X 于是,特征值问题为⎩⎨⎧==<<=+0.)((0))(0 , 0)()("l X X l x x X x X λ 特征值 2)(l n n πλ=, 特征函数为 x ln x X n πsin )(=,),2,1( =n . 而 )1,2,(n )(])[(22 ==+-t b lan n n eA t T π.取 11])[((*) . sin),(22x ln eA t x u n t b lan n ππ∑∞=+-=利用 ]0[ sinl x ln ,在⎭⎬⎫⎩⎨⎧π上的正交性,可定出 ⎰==ln n d ln l A 0),2,1( sin)(2 ξπξξϕ. 2(*) 1(*),2(*)给出所求混合问题的形式解.附注:若令 ),( ),,(),(2t x v t x v e t x u t b 则-=满足⎪⎩⎪⎨⎧===><<==== ).( ,0),0,0( 002x v v v t l x v a v t l x x xx t ϕ用分离变量法求得lxn eA t x v t lan n n sin),(2)(1ππ-∞=∑=. 而n A 同2(*),这恰与上面结果一致.习题五用Fourier 变换法求初值问题⎩⎨⎧=>++== .0),0( ),(202t xx t u t t x f tu u a u 的形式解.解:方程和初始条件两端关于x 做Fourier 变换(视t 为参数),并记),(~)],([ , ),(~)],([t f t x f F t u t x u F ξξ==.则原问题化为常微分方程的初值问题:⎪⎩⎪⎨⎧=>++-=)( .0)0,(~),0( ),(~~ 2~~22为参数ξξξξu t t f u t u a dtu d 其解为 ττξξτξτξd e f e e e t u a tt a t 2222220),(~),(~⋅⋅⋅=⎰--. 故 )],(~[),(1t uF t x u ξ-= ττξττξττξτξττξτξτξd e f F ee d ef e F e d e f e e e F ta t t a tt t t a t a t t ⎰⎰⎰-----------⋅⋅⋅=⋅⋅⋅=⎥⎦⎤⎢⎣⎡⋅⋅⋅=01)(0101]),(~[]),(~[),(~)(22222222222222ττπτττd et a F x f F F e e tt a x t]])(21[)],([[0)(412222⎰-----⋅⋅=ττπτττd et a x f F F e e tt a x t]])(21*),([[0)(412222⎰-----⋅=τξττξπτξτd d et f e a ett a x t ]1),([20)(4)(2222⎰⎰---∞∞---=即为所求.习题六1.求边值问题⎪⎪⎩⎪⎪⎨⎧≤≤=≤≤==<≤≤<≤=++=== )(0 )( ),0( 0),20 ,0( 01102αθθρπαθρρρραθθθθρρρf u l u u l u u u l 的形式解.解:用分离变量法:令 )()(θρΘ=R u ,代入泛定方程可得)( "'"2λρρ=ΘΘ-=+RR R ,因而 0)()("=Θ+Θθλθ,0)()()('"2=-+ρλρρρρR R R (Euler 方程).由边界条件 00====αθθu u,得 0)()0(=Θ=Θα.于是特征值问题为,0)()0(),0( 0)()("⎩⎨⎧=Θ=Θ<<=Θ+Θααθθλθ 特征值 2)(απλn n =,特征函数为 )1,2,( sin)( ==Θn n n θαπθ.而 Euler 方程 0'"2=-+R R R λρρ 的解 απαπρρρn n D C R -+=)(.为保证有界性应取 0=D ,从而 ),2,1( )( ==n C R n n n απρρ.取 ∑∑∞=∞==Θ=11sin)()(),(n n n n n n n C R u απθρθρθραπ. 1(*)由边界条件 )(θρf ul ==,应有 ∑∞==1sin )(n n n n lC f απθθαπ.由 ⎭⎬⎫⎩⎨⎧απθn sin在 ],0[α上的正交性,可得),2,1( sin)( 2==⎰n d n f l C n n ϕαπϕϕαααπ. 2(*)1(*) ,2(*)给出所求问题的形式解.2.用Green 函数法求解上半平面Dirichlet 问题⎪⎩⎪⎨⎧∞→+=>=+=. ),( ),0( 0220有界时,u y x x f u y u u y yy xx 解:根据二维Poisson 方程Dirichlet 问题⎩⎨⎧=∈-=+∂ ),(D.),( ),,(2y x f u y x y x u u Dyy xx πρ 解的积分表达式P PDDdl n M P G P f dxdy M M G M y x u M u ∂∂-==⎰⎰⎰∂),()(21),()(),()(00000πρ(其中0M 是D 内任一点,P n是边界D ∂上点P 的外法线方向). 其中 满足而 ),( ),,(1ln),(0000M M g M M g r M M G MM -=⎪⎩⎪⎨⎧∂∈=∈=∆).( 1ln ),g(),( 0),(000D P r M P D M M M g PM M),(0M M G 称为Green 函数,找),(0M M G 的问题归结为“特定装置下”找感应电荷所产生的电势),(0M M g -.对上半平面0>y 而言,若在0M 处放置单位正电荷,它在M 处产生的电势为01lnMM r ,则感应电荷应放在0M 关于0=y 的对称点'0M 处,电量为 -1,它于M 处产生的电势为'1lnMM r -,从而Green 函数为'1ln1ln),(0MM MM r r M M G -=20202020)()(ln )()(ln y y x x y y x x ++-+-+--=.故所求解为⎰⎰⎰⎰∞∞-=∞∞-=∞∞-=∞∞-+-=∂∂=-∂∂-=∂∂-=.)()()(21 )()(21)(21),(22000000dx yx x x f y dx yG x f dxy G x f dx n G x f y x u y y y ππππ。
数学物理方程习题解答案
数学物理方程习题解答习题一1,验证下面两个函数:(,)(,)sin x u x y u x y e y ==都是方程0xx yy u u +=的解。
证明:(1)(,)u x y =因为32222222222222223222222222222222222222222211()22()2()()11()22()2()()0()()x xx y yy xx yy x u x x y x y x y x x x y u x y x y yu y x y x y x y y y y x u x y x y x y y x u u x y x y =−⋅⋅=−+++−⋅−=−=++=−⋅⋅=−+++−⋅−=−=++−−+=+=++所以(,)u x y =0xx yy u u +=的解。
(2)(,)sin x u x y e y = 因为sin ,sin cos ,sin x x x xx xxy yy u y e u y e u e y u e y=⋅=⋅=⋅=−⋅所以sin sin 0xxxx yy u u e y e y +=−=(,)sin x u x y e y =是方程0xx yy u u +=的解。
2,证明:()()u f x g y =满足方程0xy x y uu u u −=其中f 和g 都是任意的二次可微函数。
证明:因为()()u f x g y =所以()(),()()()()()()()()()()()()0x y xy xy x y u g y f x u f x g y u f x g y uu u u f x g y f x g y g y f x f x g y ''=⋅=⋅''=⋅''''−=⋅−⋅⋅=得证。
3, 已知解的形式为(,)()u x y f x y λ=+,其中λ是一个待定的常数,求方程 430xx xy yy u u u −+= 的通解。
数学物理方程复习
一、填空题1、物理规律反映同一类物理现象的共同规律,称为___________。
2、在给定条件下求解数学物理方程,叫作____________________。
3、方程20tt xx u a u -=称为_________方程4、方程20t xx u a u -=称为_________方程5、静电场的电场强度E是无旋的,可用数学表示为_____________。
6、方程0j Ñ×=称为_____________的连续性方程。
7、第二类边界条件,就是______________________________________。
8、第一类边界条件,就是______________________________________。
9、00(0,)(0,)x x u x t u x t -=+称为所研究物理量u 的_____________。
10、00(0,)(0,)u x t u x t -=+称为所研究物理量u 的_____________。
11、对于两个自变量的偏微分方程,可分为双曲型、________和椭圆型。
12、对于两个自变量的偏微分方程,可分为双曲型、抛物线型和________。
13、分离变数过程中所引入的常数l 不能为_____________。
14、方程中,特定的数值l 叫作本征值,相应的解叫作_____________。
15、分离变数法的关键是________________________代入微分方程。
16、非齐次振动方程可采用______________和冲量定理法求解。
17、处理非齐次边界条件时,处理非齐次边界条件时,可利用叠加原理,可利用叠加原理,可利用叠加原理,把非齐次边界条件问题转化另一把非齐次边界条件问题转化另一_________的齐次边界条件问题。
18、处理非齐次边界条件时,处理非齐次边界条件时,可利用叠加原理,可利用叠加原理,可利用叠加原理,把非齐次边界条件问题转化另一把非齐次边界条件问题转化另一_________的齐次边界条件问题。
数理方程期末试题及答案
带入微分方程求解得:
k
a2
A 2
则得通解
T1
t
C1
cos
n l
a
t
D1
sin
n l
a
t
a2
A 2
sin t
带入初始条件得: C1
0,
D1
A a2 2
l a
则原定解问题的解为
u x,t
A a2 2
l sin a t cos
a l
l
x
2、 求解下列初值问题:(10 分)
uuttx,0u
xx
数; (3) 将形式解带入泛定方程以及初始条件,求解待定函数 Tn(t).
4、试述行波法的适用范围,并写出无限长弦自由振动的达朗贝尔公式。 答:行波法(特征线法)对双曲型方程是有效的,沿着双曲型方程两条特征线做
自变量替换总可以把双曲型方程化为可积形式,获得通解,由此行波法仅适用于
无界条件的波动方程。
3x x ,t sin x,ut x,0 x
0
解:应用达朗贝尔公式: u 1 (x at) (x at) 1
xat
( )d
其中
2
2a xat
,
x sin x, x x ,带入上式得:
u
1 2
sin
x
at
sin
x
at
1 2a
xat
d
xat
sin x cos at t
数学物理方程期末试题答案
一、 简述题:(每题 7 分,共 28 分) 1、 简述数学物理中的三类典型方程,并写出三类方程在一维情况下的具体形
式。
答:波动方程:
2u t 2
数学物理方程赵忠奎答案
数学物理方程赵忠奎答案
1.赵忠奎试题:物体运动的位移为S(t)=9t^2,求此物体在t=3秒时
的速度?
2.赵忠奎试题:用能量守恒定律求出一物体在没有外力作用下的运动
路径?
答案:答案是匀加速直线运动,即匀速运动。
3.赵忠奎试题:假定一个小球以初始速度V0自由落体,请求出它的
位移量s?
答案:由于小球自由落体,则其加速度为g,位移量s=V0t+1/2gt^2,其中t为小球落体的时间。
4.赵忠奎试题:一个摩擦力为f的物体在水平面上运动,它的动量为p,请求出它的速度?
答案:速度v=p/(m-f/a),其中m为物体的质量,f为摩擦力,a为
物体在水平面上的加速度。
5.赵忠奎试题:物体运动的轨迹方程为x=2t+1,请问该物体的运动速度?
答案:该物体的运动速度为v=2m/s。
6.赵忠奎试题:请利用势能定理求解受外力作用的一个物体的势能?
答案:由势能定理:W=∆K,其中W为受外力作用的一个物体的势能,
K为物体的动能,∆K为物体动能变化的量。
7.赵忠奎试题:一个物体以初始速度v0,加速度a加速运动,请求出这个物体到终点时的速度?
答案:这个物体到终点时的速度为v=v0+at。
数学物理方程答案(全)
SY (ux (x dx,t) ux (x,t)) Sdxutt
utt
Y P
uxx
杆的一端固定,有 u(0,t) 0 ,另一端为自由端有 ux (x,l) 0
由于弦在出事时刻处于静止状态,即初速度为零,故 ut (0,t) 0
在 t 0 时刻,整个杆被纵向拉长 d ,则单位杆长的伸长量为 d ,故 x 点处的伸长 l
(3)特征方程为
4( dy )2 8( dy ) 3 0
dx
dx
解得
dy 3 和 dy 1 dx 2 dx 2
习题 2.2
1.一根半径为 r,密度为 ,比热为 c,热传导系数为 k 的匀质圆杆,如同截面
上的温度相同,其侧面与温度为 1 的介质发生热交换,且热交换的系数为 k1 。
试导出杆上温度 u 满足的方程。 解:
0
x
x+dx
取微元在 (x, x dx) 之间,在时间 t 内
x
从左右两截面流入的热量,有热传导方程可得
dT g dx 对上式进行积分,并且利用在 x 0 处的张力为T x0 gl 可求得 T (x) g(l x)
对于(2)式 sin2 tan2 ux (x dx,t) sin 1 tan 1 ux (x,t)
将上述结果代入(2)式得出
T (x dx,t)ux (x dx,t) T (x)ux (x,t) uttdx
2 x ux)
Y P
x
(x2
u x
)
1 x2
4.一根长为 L、截面面积为 1 的均匀细杆,其 x=0 端固定,以槌水平击其 x=L
端,使之获得冲量 I。试写出定解问题。
解:由 Newton 定律: SYux (x dx,t) YSux (x,t) Sdxutt ,其中,Y 为杨
数学物理方程2-3章课后部分习题答案
数学物理方程 李明奇主编 电子科技大学出版社2-3章部分习题答案习题2.14.一根长为L 、截面面积为1的均匀细杆,其x=0端固定,以槌水平击其x=L 端,使之获得冲量I 。
试写出定解问题。
解:由Newton 定律: tt x x Sdxu t x YSu t dx x SYu ρ=-+),(),(,其中,Y 为杨氏模量,S 为均匀细杆的横截面积,x u 为相对伸长率。
化简之后,可以得到定解问题为:⎪⎪⎩⎪⎪⎨⎧-==========)(|,0|0|,0|)/(0002L x I u u u u u a u Y u t t t L x x x xx xx tt δρρ。
习题2.23.设物体表面的绝对温度为u ,它向外辐射出去的热量,按斯特凡-波尔兹曼定律正比于4u ,即dSdt ku dQ 4=,设物体与周围介质之间,只有热辐射而无热传导,周围介质的绝对温度为已知函数),,,(t z y x ϕ,。
试写出边界条件。
解:由Fourier 热传导实验定律dSdt nuk dQ ∂∂-=1,其中1k 称为热传导系数。
可得dSdt u k dSdt n uk )(441ϕ-=∂∂-,即可得边界条件:)(441ϕ--=∂∂u k k nus。
习题2.34.由静电场Gauss 定理⎰⎰⎰⎰⎰⋅=⋅VsdV dS E ρε01,求证:0ερ=⋅∇E ,并由此导出静电势u 所满足的Poisson 方程。
证明:⎰⎰⎰⎰⎰⎰⎰⎰⋅=⋅=⋅VVsdV dV divE dS E ρε01,所以可以得到:0ερ=divE 。
由E divE ⋅∇=与u E -∇=,可得静电势u 所满足的Poisson 方程:02ερ-=∇u 。
习题2.42.求下列方程的通解:(2):;032=-+yy xy xx u u u (5):;031616=++yy xy xx u u u 解:(2):特征方程:03)(2)(2=--dxdydx dy 解得:1-=dx dy 和3=dxdy。
数学物理方程课后参考答案第二章
第 二 章 热 传 导 方 程§1 热传导方程及其定解问题的提1. 一均匀细杆直径为l ,假设它在同一截面上的温度是相同的,杆的表面和周围介质发生热交换,服从于规律dsdt u u k dQ )(11-= 又假设杆的密度为ρ,比热为c ,热传导系数为k ,试导出此时温度u 满足的方程。
解:引坐标系:以杆的对称轴为x 轴,此时杆为温度),(t x u u =。
记杆的截面面积42l π为S 。
由假设,在任意时刻t 到t t ∆+内流入截面坐标为x 到x x ∆+一小段细杆的热量为t x s xuk t s x u k t s x u k dQ x x x x ∆∆∂∂=∆∂∂-∆∂∂=∆+221 杆表面和周围介质发生热交换,可看作一个“被动”的热源。
由假设,在时刻t 到t t ∆+在截面为x 到x x ∆+一小段中产生的热量为()()t x s u u lkt x l u u k dQ ∆∆--=∆∆--=111124π又在时刻t 到t t ∆+在截面为x 到x x ∆+这一小段内由于温度变化所需的热量为()()[]t x s tuc x s t x u t t x u c dQ t ∆∆∂∂=∆-∆+=ρρ,,3由热量守恒原理得:()t x s u u lk t x s x uk t x s t u c x t ∆∆--∆∆∂∂=∆∆∂∂11224ρ消去t x s ∆∆,再令0→∆x ,0→∆t 得精确的关系:()11224u u l kxu k t u c --∂∂=∂∂ρ或 ()()11222112244u u l c k xu a u u l c k x u c k t u --∂∂=--∂∂=∂∂ρρρ 其中 ρc k a =22. 试直接推导扩散过程所满足的微分方程。
解:在扩散介质中任取一闭曲面s ,其包围的区域 为Ω,则从时刻1t 到2t 流入此闭曲面的溶质,由dsdt nuDdM ∂∂-=,其中D 为扩散系数,得 ⎰⎰⎰∂∂=21t t sdsdt nuDM 浓度由u 变到2u 所需之溶质为()()[]⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰⎰ΩΩΩ∂∂=∂∂=-=2121121,,,,,,t t tt dvdt t uC dtdv t u C dxdydz t z y x u t z y x u C M两者应该相等,由奥、高公式得:⎰⎰⎰⎰⎰⎰⎰⎰ΩΩ∂∂==⎥⎦⎤⎢⎣⎡⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎪⎭⎫ ⎝⎛∂∂∂∂+⎪⎭⎫ ⎝⎛∂∂∂∂=21211t t t t dvdt t uC M dvdt z uD z y u D y x u D x M 其中C 叫做孔积系数=孔隙体积。
数学物理方程考试试题及解答
数学物理方程试题(一)一、填空题(每小题5分, 共20分)1.长为 的两端固定的弦的自由振动, 如果初始位移为 , 初始速度为x 2cos 。
则其定解条件是2.方程.的通解................3.已知边值问题 , 则其固有函数 =4.方程0)(222'"2=-++y n x xy y x α的通解为 二.单项选择题(每小题5分, 共15分)1. 拉普拉斯方程 的一个解是.. )(A )xy e y x u x sin ),(= (B )22),(y x y x u +=(C )221),(y x y x u += (D )22ln),(y x y x u += 2.一细杆中每点都在发散热量, 其热流密度为 ,热传导系数为 , 侧面绝热,体密度为 ,比热为 , 则热传导方程....)(A )ρc t x F x u a t u),(22222+∂∂=∂∂ (B )ρc t x F x u a t u ),(222+∂∂=∂∂ (C ) ρc t x u x F a t F ),(22222+∂∂=∂∂ (D) ρc t x u x F a t F ),(222+∂∂=∂∂ (其中ρc k a =2) 3.理想传输线上电压问题( 其中CL a 12=)的解为( ) (A ))(cos ),(at x A t x u +=ω (B )t a x A t x u ωωcos cos ),(=(C )t a x A t x u ωωsin cos ),(= (D ))(cos ),(t a x A t x u -=ω1. 三.解下列问题2. ( 本题8分) 求问题 ⎪⎩⎪⎨⎧==∂∂+∂∂x ex u yu x u 38)0,(03的解3. ( 本题8分)⎪⎪⎩⎪⎪⎨⎧=-==∂∂∂222),0(,cos 1)0,(6y y u x x u y x y x u...本题8分.求问. 的解1. 四.用适当的方法解下列问题2. ( 本题8分) 解问题 ⎪⎩⎪⎨⎧+-=∂∂=∂∂2222321)0,(x x x u x u a t u 2.( 本题8分) 解问题 ⎪⎪⎩⎪⎪⎨⎧=∂∂+=∂∂+∂∂+∂∂=∂∂==202202222222226,32)(y t uxz y u z u y u x u a t u t t 五. ( 本题10分)解混合问题:六. ( 本题15分)用分离变量法解下列混合问题:⎪⎪⎪⎩⎪⎪⎪⎨⎧=∂∂-===∂∂=∂∂=xt u x x x u t u t u x u a t u t 2sin 3,)(2)0,(0),(),0(022222ππ 一.单项选择题(每小题4分, 共20分)1.(D..2.(B..3.(D..4.(D )二.填空题(每空4分, 共24分)1....2...3.. ,4.)(x X n =cos ,(0,1,2,3,)2n n x B n π= 5.通解为223(,)()()2u x t x y f x g y =++ 三.解下列问..本题7分.1. 求问题 的解解: 设 (2分)代入方程,330,1m m +==- (6分)所以解为 3(,)8x y u x t e -= (7分)2. ( 本题7分) 求问题 ⎪⎪⎩⎪⎪⎨⎧=∂∂=∂∂=∂∂=20222223,2sin )0,(x t ux x u x u a t u t 的解 解: 由达朗贝尔公式, 得211(,)[sin 2()sin 2()]322x at x at u x t x at x at d aξξ+-=++-+⎰(3分) 223cos 2sin 23at x x t a t =++ (7分)四.用适当的方法解下列问题1. .本题7分.解问.解: 设代入方程,令 2066A A a x''=⎧⎨=+⎩ 显然成立 解为 22(,)12366u x t x x a t xt =-+++2.( 本题7分) 解问题 ⎪⎪⎩⎪⎪⎨⎧=∂∂++=∂∂+∂∂+∂∂=∂∂==202202222222226,32)(y t u yz y x u z u y u x u a t u t t 解: 设 (2分)代入方程22326[(212)(12)]A Bt a y At t Bt +=++∆++∆ (4分)令 , 显然成立, 解为322222632),(t a t y t a yz y x t x u +++++=五. ( 本题7分)解混合问题:⎪⎪⎪⎩⎪⎪⎪⎨⎧===∂∂=∂∂x x u t u t u x u a t u πsin 2)0,(0),1(),0(222 解1(,){(,)}u x t L U x s -=222sin a t e x ππ-= 六. ( 本题15分)用分离变量法解下列混合问题:⎪⎪⎪⎩⎪⎪⎪⎨⎧=∂∂-===∂∂=∂∂=xt u x x x u t u t u x u a t u t 2sin 3,)(2)0,(0),(),0(022222ππ 解: 设 代入方程及边界200(0)()0T a T X X X X λλπ''⎧+=⎪''+=⎨⎪==⎩22(),sin n n n n X nx πλπ=== (cos sin )sin n n n u C ant D ant nx =+1(,)(cos sin )sin n n n u x t C ant D ant nx ∞==+∑其中 3028[1(1)]()sin n n C x x nxdx n ππππ--=-=⎰ 00(2)23sin 2sin 3(2)n n D x nxdx n aππ≠⎧⎪==⎨=⎪⎩⎰ 所以解为3138[1(1)](,)sin 2sin 2cos sin n n u x t at x ant nx a n π∞=--=+∑2009-2010学年第一学期数学物理方程试题一、 填空题(每小题4分, 共24分)1.方程.的特征线..........2.长为 的弦做微小的横振动, 、 两端固定, 且在初始时刻处于水平状态, 初始速度为 .则其定解条件.................3.方程 的通解.........4.已知边值问. .. 则其固有函数)(x X n =5.方程 的通解............6...........二. 单项选择题(每小题4分, 共20分)1.微分方程.是..)(A )三阶线性偏微分方程 (B )三阶非线性偏微分方程(C )三阶线性齐次常微分方.....(D )三阶非线性常微分方程2. 拉普拉斯方程 的一个解是.. )(A )xy e y x u x sin ),(= (B )22),(y x y x u +=(C )221),(y x y x u += (D )22ln),(y x y x u += 3.一细杆中每点都在发散热量, 其热流密度为 ,热传导系数为 , 侧面绝热,体密度为 ,比热为 , 则热传导方程....)(A )ρc t x F x u a t u),(22222+∂∂=∂∂ (B )ρc t x F x u a t u ),(222+∂∂=∂∂ (C ) ρc t x u x F a t F ),(22222+∂∂=∂∂ (D) ρc t x u x F a t F ),(222+∂∂=∂∂ (其中ρc k a =2) 4.理想传输线上电压问题(A ))(cos ),(at x A t x u +=ω (B )t a x A t x u ωωcos cos ),(=(C )t a x A t x u ωωsin cos ),(= (D ))(cos ),(t a x A t x u -=ω5.单位半径的圆板的热传导混合问题⎪⎩⎪⎨⎧=<=<∂∂+∂∂=∂∂)()0,(,),(,0),1()1()1(222ρρρρρρρf u M t u t u u u a t u 有形如( )的级数解。
数学物理方程 参考答案
1.求下列波动方程Cauchy 问题的解: (2)2005,tt xx t tt u a u u u x==⎧=⎪⎨==⎪⎩解:根据达朗贝尔公式可得521)55(21),(+=++=⎰+-xt d a t x u at x atx ξξ6.求下列强迫振动的Cauchy 问题的解:(1)⎩⎨⎧==+===2002,5x u u e u a u t t t xxx tt解:令)(),(),(x w t x v t x u +=,代入原方程,得xxx xx tt ew a v a v ++=22令2)(a ex w x-=可得⎪⎩⎪⎨⎧=+====222,5xv ae v v a v t tx t xxtt由达朗贝尔公式可得531)(2121)5()5(21),(3222222++++=+⎥⎦⎤⎢⎣⎡+++=+-+-+-⎰t a tx eead aaea e t x v atx atx atx atx atx at x ξξ所以原问题的解为2232211(,)()523x x atx ate u x t v w eetx a t aa-+=+=++++-7.求解下列定解问题:⎩⎨⎧==>+∞<<-∞=-++==)(),(0,,020022x u x u t x u a u u u t t t xx t tt ψϕεε解:令)0(),,(),(>=-ββt x v et x u t,代入原方程得:)2()(2222=+-+-+-v v v a v t xx tt βεβεβε取εβ=,可得⎩⎨⎧+==>+∞<<-∞=-==)()(),(0,,0002x x v x v t x v a v t t t xx tt εϕψϕ 由达朗贝尔公式得:[][]11(,)()()()()22x at x atv x t x at x at d aϕϕψξεϕξξ+-=++-++⎰所以,原定解问题的解为:[][]11(,)()()()()22x at t tx atu x t x at x at d eaeββϕϕψξεϕξξ+-=++-++⎰习题4.22.求解下列定解问题2000,0,00,0()tt xx t t t x x u a u x t u u u h t ===⎧=<<+∞>⎪==⎨⎪=⎩解:通解为12(,)()()u x t f x at f x at =++-由初始条件1212(,0)()()0(1)(,0)()()0(2)t u x f x f x u x af x af x =+=⎧⎨''=-=⎩对(2)式积分可得121()()f x f x C -=则有1112()2,0()2C f x x C f x ⎧=⎪⎪≥⎨⎪=⎪⎩0x at +≥恒成立,但是x at -可能小于零当0x at -<时1212()()()()()()f at f at h t f f h a ξξξ''+=⎧⎪⎨''+-=⎪⎩令0at ξ=>,积分可得12120()()()(0)(0)f f h d f f aξξξξξ+-=+-⎰令aξη=上式变为12120()()()(0)(0)a f f a h d f f ξξξηη+-=+-⎰21101110()()()()2()2a a a f f a h d C C a h d C C a h d ξξξξξηηηηηη⎡⎤-=-+⎢⎥⎣⎦=--=--⎰⎰⎰所以1210,02()(),02a C f C a h d ξξξηηξ⎧-≥⎪⎪=⎨⎪--<⎪⎩⎰则有1210,2()(),2a C x t a f x at C xa h d t a ξηη⎧-≤⎪⎪-=⎨⎪-->⎪⎩⎰又因为11()2C f x at +=所以00,(,)(),a x t a u x t xa h d t a ξηη⎧≤⎪⎪=⎨⎪->⎪⎩⎰习题4.31.求解下列定解问题200,,,,0,tt t t t u a u x y z t u yz u xz ==⎧=∆-∞<<+∞>⎪⎨==⎪⎩解:对于三维波动方程,其解为1(,,)(,,)(,,,)41(,,)1(,,)44x y z x y z u x y z t dS dS a t atatx y z x y z dS dSa tataatϕψπϕψππ''''''∂⎡⎤=+⎢⎥∂⎣⎦''''''∂⎡⎤=+⎢⎥∂⎣⎦⎰⎰⎰⎰⎰⎰⎰⎰其中2sin cos ,02,0sin sin cos ()sin x x at y y at z z at dS at d d θϕϕπθπθϕθθϕθ'=+≤≤≤≤⎧⎪'=+⎪⎨'=+⎪⎪=⎩在本题目中(,,)x y z yzϕ=,(,,)x y z xz ψ=()()2222222001(,,)41(sin sin )(cos )()sin 41sin sin cos sinsin sin cos sin 412sin 2sin cos 4x y z dS a t aty r z r at d d a tatat yz aty at a t d d a tat yz aty d a t πππππϕπθϕθθϕθπθθθθϕθθϕϕθππθπθθθπ'''∂⎡⎤⎢⎥∂⎣⎦∂++⎡⎤=⎢⎥∂⎣⎦∂⎛⎫=+++ ⎪∂⎝⎭∂=+∂⎰⎰⎰⎰⎰⎰ ()1404at yz a t yzππ⎛⎫ ⎪⎝⎭∂⎛⎫=+⎪∂⎝⎭=⎰1(,,)4x y z dSa atxztψπ'''=⎰⎰则(,,,)u x y z t yz xzt=+3.利用三维泊松公式求解下列问题220,,,,00,tt t t t u a u x y z t u u x yz ==⎧=∆-∞<<+∞>⎪⎨==+⎪⎩ 解:对于三维波动方程,其解为1(,,)(,,)(,,,)41(,,)1(,,)44x y z x y z u x y z t dS dS a tatatx y z x y z dS dSa t ataatϕψπϕψππ''''''∂⎡⎤=+⎢⎥∂⎣⎦''''''∂⎡⎤=+⎢⎥∂⎣⎦⎰⎰⎰⎰⎰⎰⎰⎰在本题目中有(,,,)0x y z t ϕ=,则有()()()()()()()()()222220220021(,,)(,,,)4sin sin sin sin cos 1sin 4sin sin sin sin cos sin 4sin sin sin 4sin 4x y z u x y z t dSaatx at y at z at at d d a at t x at y at z at d d t x at d d yzttx d ππππππψπθϕθϕθθϕθπθϕθϕθθϕθπθϕθϕθπθπ'''=⎛⎫++++ ⎪=⎪⎝⎭=++++=++=⎰⎰⎰⎰⎰⎰⎰⎰ ()()22222322000022223200232sin sin 2sin sin 4sin sin 043d a t d d xat d d yztt x a t d d yzta t x t yztππππππππϕθθϕϕθθϕϕθπθθϕϕπ+++=+++=++⎰⎰⎰⎰⎰⎰⎰⎰习题4.43.导出二维Cauchy 问题解的表达式200(,,),,,00,0tt t t t u a u f x y t x y t u u ==⎧=∆+-∞<<+∞>⎪⎨==⎪⎩ 解:利用齐次化原理求解 如果(,,,)w x y t τ是定解问题20,(,,)tt t tt W a W W W f x y τττ==⎧=∆⎪⎨==⎪⎩的解 则0(,,)(,,,)tu x y t w x y t d ττ=⎰即为定解问题200(,,)0,0tt t t t u a u f x y t u u ==⎧=∆+⎪⎨==⎪⎩的解 对于0(,,)(,,,)tu x y t w x y t d ττ=⎰显然存在如下的关系(,,,)0t uw x y t d ττ===⎰(,,,)t t t u w w w x y t d d tttττττ=∂∂∂=+=∂∂∂⎰⎰此时有00t ut =∂=∂又有222222222000(,,)(,,)t tttu w w w wd f x y t a wd f x y t a d tttx y ττττ=⎛⎫∂∂∂∂∂=+=+∆=++ ⎪∂∂∂∂∂⎝⎭⎰⎰⎰且222222220tuuw wd x y xy τ⎛⎫∂∂∂∂+=+ ⎪∂∂∂∂⎝⎭⎰将上式代入22u t∂∂表达式可得2222222(,,)(,,)u u u f x y t a f x y t a u txy ⎛⎫∂∂∂=++=+∆ ⎪∂∂∂⎝⎭因此齐次化原理得以证明.由齐次方程柯西问题解的泊松公式可得1(,,)(,,,)2Mf w x y t aττπ=⎰⎰所以,原问题的解为()201(cos ,sin ,)(,,)2t a t f x r y r u x y t d aτπθτπ-++=⎰⎰⎰习题5.1 1.若[]()()F g x f ω=,求证[]()2()F f x g πω=-.证明:由傅里叶反变换式1()()2j xg x f ed ωωωπ+∞-∞=⎰,将式中自变量x 换为x -,得1()()2j xg x f ed ωωωπ+∞--∞-=⎰将上式变量x 换为ω,而把ω换为x ,得 1()()2j xg f x e dx ωωπ+∞--∞-=⎰ 即[])(2)(ωπ-=g x f F2.求证 (1)1,0y Fe y ω--⎡⎤>⎣⎦(2)00()()j xF ef x f ωωω⎡⎤=-⎣⎦证明:根据Fourier 变换可得出000++()0()()()()j xj xj xj xF e f x f x e edxf x edxf ωωωωωωω∞--∞∞---∞⎡⎤⎣⎦===-⎰⎰(3)[]()()f aF f at aω=证明:若0>a ,则)(at f 的傅里叶变换为[]+()()j tF f at f at edtω∞--∞=⎰令at x =,则adtdx=代入上式,可得[]+1()()j x adx F f x f x ef aa a ωω-∞-∞⎛⎫==⎪⎝⎭⎰若0<a ,则类似地有[]1()Ff at f aa ω⎛⎫=-⎪⎝⎭综上所述[]()()f aF f at aω=3.求函数的Fourier 变换 (1) ()xf x e -= 证明:2cos sin 22cos 1xxxxj xxF e eedx exdx i exdxexdx ωωωωω+∞+∞+∞------∞-∞-∞+∞-⎡⎤==-⎣⎦==+⎰⎰⎰⎰由于积分区间是关于坐标轴对称,且积分函数是个奇函数故sin 0xexdx ω+∞--∞=⎰因此2022cos 1x xF e e xdx ωω+∞--⎡⎤==⎣⎦+⎰(2) 2()xf x eπ-=证明:直接利用公式[]2222()cos sin 2cos xj xxxxF f x e edxexdx i exdxexdxπωπππωωω+∞---∞+∞+∞---∞-∞+∞-==-=⎰⎰⎰⎰根据公式22240cos xa ba bexd ωωω-+∞-=⎰则[]22441()22Ff x eeωωππ--=⋅=(3)2()cos f x ax = 证明:[]2()cos j xF f x ax edxω+∞--∞=⋅⎰根据cos 2izize ez -+=上式可以变为2222222222222()()2424()42cos 211221122112212j xjaxjaxj xjaxj xjaxj xjax j xjax j xja x jja x jaaa ajja x aaax edxee edxe edx eedxedx edxedx edxeeωωωωωωωωωωωω+∞--∞-+∞--∞+∞+∞----∞-∞+∞+∞----∞-∞--+++∞+∞-∞-∞--⋅+==+=+=+=⎰⎰⎰⎰⎰⎰⎰⎰22()4212jja x aadx eedxωω++∞+∞-∞-∞+⎰⎰令)2x aωξ=-以及)2x aωη=+上式变为222222222222()()4242444401122112211jja x jja x aaaajjj j aajjj j aaeedx eedxeedeede d ed ωωωωωωξηωωξηξη--++∞+∞-∞-∞-+∞+∞-∞-∞-+∞+∞+==+⎰⎰⎰⎰⎰⎰再利用公式2402jj ed πξξ+∞=⎰上式可变为22222222440()()4444()()444422)44jjj j aaj j a a j j a a ed ed ee e e aωωξηωπωπωπωπξηωπ-+∞+∞------+⎤=+⎥⎥⎦⎤+⎥=⎥⎥⎣⎦=-⎰⎰5.求()0axf x ea -=>,,Fourier 正弦与余弦变换.解:由定义,得:2202cos 1cos 11cos cos 1sin 1sin 1sin cos 1cos axaxaxaxaxaxaxaxaxexdxxdeaxe ed xaaexdxa axdea axe xdea a xdeaaωωωωωωωωωωωωωω+∞-+∞-+∞+∞--+∞-+∞-+∞+∞--+∞-=-=-+=-=+⎡⎤=+-⎢⎥⎣⎦=+⎰⎰⎰⎰⎰⎰⎰由此得出222cos axaexdx a ωω+∞-=+⎰即22()c a f a ωω=+同理可得22ˆ()sin axs f exdx a ωωωω+∞-==+⎰习题5.21. 用Fourier 变换法求解定解问题 ⎩⎨⎧==>∈=0)0,(,sin )0,(0,,2x u x x u t R x u a u t xx tt 解:对于初值问题关于x 作Fourier 变换,得:[]2222d (,)(,),,0d (,0)sin ,(,0)0t u t a u t x R t t u F x uωωωωω⎧+∈>⎪⎨⎪==⎩该方程变为带参数ω的常微分方程的初值问题. 解得12(,)ja t ja t ut C e C e ωωω-=+ 于是1212(,0)(sin ),(,0)()0t uF x C C u ja C C ωωω==+=-= 则由[]121sin 2C C F x ==可得[]1(,)sin ()2ja tja tut F x eeωωω-=+作像函数(,)ut ω 的Fourier 逆变换 [][][]11111(,)[(,)]1sin ()21sin (sin )211sin (sin )221[sin()sin()]2sin cos ja t ja t ja t ja t ja t ja t u x t F u t F F x e e F F x e F x e FF x e F F x e x at x at x atωωωωωωω--------=⎡⎤=+⎣⎦⎡⎤=+⎣⎦⎡⎤⎡⎤=+⎣⎦⎣⎦=-++=2.求解下列定解问题2cos ,0,0(,0)0,(,0)0,lim (,)0(0,)0tt xx t x x u a u t x t u x u x u x t u t →+∞⎧=+<<+∞>⎪===⎨⎪=⎩ 解:对自变量t 取Laplace 变换可得⎪⎩⎪⎨⎧=+∞=+=-0),(~,0),0(~1~~22222s u s us s dx ud a u s x求解常微分方程,得)(1~22s s Be Ae u xa sx as+++=-ω于是)1(1,02s s B A +-==所以]1[)1(1~2xas es s u --+=且111222()22211L (1)L L (1)(1)(1)R e s ,R e s ,(1)(1)s x sa x a xs t sta k k kke e s s s s s s e e s s s s s s ------⎡⎤⎡⎤⎡⎤⎢⎥-=-⎢⎥⎢⎥⎢⎥+++⎣⎦⎣⎦⎢⎥⎣⎦⎡⎤⎡⎤⎢⎥=-⎢⎥⎢⎥++⎣⎦⎢⎥⎣⎦∑∑其中k s 是u ~的极点 由于01=s ,js =2,j s -=3都是一级极点,所以222202R e s ,lim lim ()lim ()(1)(1)(1)(1)11()2(1cos )12sin2st st st stk s s j s j kjtjte e e e s s s j s j s s s s s s s s eet t→→→--⎡⎤=⋅+-++⎢⎥++++⎣⎦=-+=-=∑2()2()2sin ,2R e s ,(1)0,x s t a k kat x x t e a a s x s s t a --⎧⎡⎤>⎪⎪⎢⎥=⎨⎢⎥+⎪≤⎢⎥⎣⎦⎪⎩∑所以,最后定解问题为22122sin 2sin ,22[]2sin ,2t at x x t aau L ut x t a--⎧->⎪⎪==⎨⎪≤⎪⎩4.求解定解问题(,),,0(,0)(),(,0)()tt xx t u u f x t x t u x x u x x ϕψ=+-∞<<+∞>⎧⎨==⎩解:首先使用分离变量法,令u VW=+,则可将原定解问题分解为200()(1)()tt xx t t t V a V V x V x ϕψ==⎧=⎪⎪=⎨⎪=⎪⎩200(,)0(2)0tt xx t t t W a W f x t W W ==⎧=+⎪⎪=⎨⎪=⎪⎩对于方程(1).对初值问题关于x 做Fourier 变化,得2222(,)(,)0(,0)(),(,0)()t d v t a v t dt v v ωωωωϕωωψω⎧+=⎪⎨⎪==⎩该方程变为带参数ω的常微分方程的初值问题.解得12(,)j at j at vt C e C e ωωω-=+ 于是1212()(,0)()(,0)()t v C C vj a C C ϕωωψωωω==+==-即有1111(,)()()()()22j at j at vt e ej a j a ωωωϕωψωϕωψωωω-⎡⎤⎡⎤=++-⎢⎥⎢⎥⎣⎦⎣⎦做像函数的Fourier 逆变换[]11111(,)(,)1111()()()()22j at j at j at j at W x t Fut Fe F e F e F e aj a j ωωωωωϕωϕωψωψωωω-------=⎡⎤⎡⎤⎡⎤⎡⎤=++-⎢⎥⎢⎥⎣⎦⎣⎦⎣⎦⎣⎦因为[][]()()()j atj atF x at eF x eωωϕϕϕω±±±== 做逆变换可得 1()()j at Fe x at ωϕωϕ-±⎡⎤=±⎣⎦又因为1()()()x at x j atj at F s ds e F s ds e j ωωψψψωω±±±-∞-∞⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎣⎦⎰⎰ 做逆变换可得11()()x at j at Fe s dsj ωψωψω±-±-∞⎡⎤=⎢⎥⎣⎦⎰因此[][]11(,)()()()()2211()()()22x atx at x at x atV x t x at x at s ds s ds a x at x at s dsaϕϕψψϕϕψ+--∞-∞+-⎡⎤=++-+-⎢⎥⎣⎦=++-+⎰⎰⎰对于方程(2).根据齐次化原理,如果(,,)w x t τ是齐次方程Cauchy 问题的解20(,)tt xx t t t w a ww w f x τττ==⎧=⎪⎪=⎨⎪=⎪⎩则0(,)(,,)tW x t w x t d ττ=⎰是原问题的解.利用变换t t τ'=-则2000(,)t t xx t t t w a w w w f x τ'''=''=⎧=⎪⎪=⎨⎪=⎪⎩ 利用达朗贝尔公式有1(,,)(,)2x at x at w x t f d a τατα'+'-'=⎰ ()()1(,,)(,)2x a t x a t w x t f d aτττατα+---=⎰可求得()0()1(,)(,)2t x a t x a t W x t f d d aττατατ+---=⎰⎰最后,[]()0()111()()()(,)222x at t x a t x atx a t u V W x at x at s ds f d d aaττϕϕψατατ++----=+=++-++⎰⎰⎰习题5.31.求证Laplace 变换的位移定理. 证明: Laplace 变换的位移定理为L ()(),Re()axef x f s a s a σ⎡⎤=-->⎣⎦ 根据Laplace 变换的定义可以求得()00L ()()()(),Re()axaxsxs a xef x ef x edx f x edx f s a s a σ+∞+∞---⎡⎤===-->⎣⎦⎰⎰3.用留数计算1221L (1)()sx ae s s ω--⎡⎤-⎢⎥+⎣⎦解:1122222211L (1)L ()()()sx sa x a e e s s s s s s ωωω----⎡⎤⎡⎤⎢⎥-=-⎢⎥⎢⎥+++⎣⎦⎢⎥⎣⎦根据L 变换的线性性质11122222211L (1)L L ()()()sx sa x ae e s s s s s s ωωω-----⎡⎤⎡⎤⎡⎤⎢⎥-=-⎢⎥⎢⎥⎢⎥+++⎣⎦⎣⎦⎢⎥⎣⎦根据留数定理可得出()12222221L (1)R e s ,R e s ,()()()xs t ssta x a k k kk eee s s s s s s s s ωωω---⎡⎤⎡⎤⎡⎤⎢⎥-=-⎢⎥⎢⎥⎢⎥+++⎣⎦⎣⎦⎢⎥⎣⎦∑∑其中k s 是极点.由于01=s ,ωj s =2,ωj s -=3都是一级极点,所以22222222022222R e s ,lim lim ()lim ()()()()()11()21(1cos )2sin2st st st stk s s j s j kj tj te e e e s s s j s j s s s s s s s s eet tωωωωωωωωωωωωωωωω→→→--⎡⎤=⋅+-++⎢⎥++++⎣⎦=-+=-=∑对于()22R e s ,()xs t a k kes s s ω-⎡⎤⎢⎥⎢⎥+⎢⎥⎣⎦∑需要分情况讨论当xt a >时,()()()()222222220()()22222R e s ,lim lim ()lim ()()()()()1121()1cos 2()sin2x x x x s t s t s t s t a a a a k s s j s j kx xj t j t aae e e e s s s j s j s s s s s s s s e e at x a at x aωωωωωωωωωωωωωωωω----→→→----⎡⎤⎢⎥=⋅+-++⎢⎥++++⎢⎥⎣⎦⎡⎤=-+⎢⎥⎣⎦-⎡⎤=-⎢⎥⎣⎦-=∑当x t a≤时,()22R e s ,0()xs t a k kes s s ω-⎡⎤⎢⎥=⎢⎥+⎢⎥⎣⎦∑综上所述,可以得出2()2222()sin ,2R e s ,()0,x s t a k kat x x t e a a s x s s t a ωωω--⎧⎡⎤>⎪⎪⎢⎥=⎨⎢⎥+⎪≤⎢⎥⎣⎦⎪⎩∑所以,最后结果为22221222222()sin sin ,122L (1)2()sin ,2s x at at x x t a a e t xs s t a ωωωωωωω---⎧->⎪⎡⎤⎪-=⎨⎢⎥+⎣⎦⎪≤⎪⎩7.求下列函数的Laplace 逆变换 (1) 5482+++s ss (2) )0(,)(222>+a a s s解:(1)对原式进行分解,得1)2(61)2(2548222++++++=+++s s s s s s则)sin 6(cos 1)2(61)2(25482212121t t e s L s s L s s s L t+=⎥⎦⎤⎢⎣⎡+++⎥⎦⎤⎢⎣⎡+++=⎥⎦⎤⎢⎣⎡+++----[查表可得](2)对原式进行分解,得22222)(14)(14)(ja s a j ja s a j a s s--+=+由于[]2)(1a s te L at+=-,得:1112222211()4()4()()41sin 2jatjats j jL L L s a a s ja a s ja j t eea t ata ----⎡⎤⎡⎤⎡⎤=-⎢⎥⎢⎥⎢⎥++-⎣⎦⎣⎦⎣⎦=-=[查表可得] 习题5.41.用Laplace 变换法解下列定解问题:(2)2000,0,00,00tt xx t t t x u a u c x t u u u ===⎧=+<<∞>⎪⎪==⎨⎪=⎪⎩解:对时间变量t 坐拉普拉斯变换222200,x x d u c a s u u dx s u =→∞⎡⎤=⎧-=-⎪⎪⎩⎣⎨⎦=求解微分方程上述微分方程. 对应的特征方程为220s a λ⎛⎫-= ⎪⎝⎭特征根为s aλ=±对应齐次方程的通解为s sxxaauAe Be-=+由于00λ=不是特征方程对应的特征根,故非齐次方程的一个特解为*uC =将特解代入原方程可得3c C s=因此原问题的解为*3s sxxa ac uu u Ae Bes-=+=++根据边界条件可得出30c A s B ⎧=-⎪⎨⎪=⎩则33s xac c uess-=-+对其做逆变换可得[]1133(,)sx a c c u x t L u L e ss ---⎡⎤==-+⎢⎥⎣⎦ 根据线性定理可将其变为1133(,)sx ac c u x t L e L s s ---⎡⎤⎡⎤=-+⎢⎥⎢⎥⎣⎦⎣⎦其中,2132c ctL s -⎡⎤=⎢⎥⎣⎦133Re s ,ssx x staa k kc c L e e e s s s ---⎡⎤⎡⎤-=-⋅⎢⎥⎢⎥⎣⎦⎣⎦∑其中,0s =是三阶极点故有()332()32302R e s ,R e s ,1lim (0),(31)!0,(),20,,s xx t s sta a k k kkxt s as c c e e s e s s s d c x s e t ds s a x t a c x x t t a a x t a ---→⎡⎤⎡⎤-⋅=⎢⎥⎢⎥⎣⎦⎣⎦⎧⎡⎤->⎪⎢⎥⎪-⎣⎦=⎨⎪≤⎪⎩⎧->⎪⎪=⎨⎪≤⎪⎩∑∑因此222(),22(,),2ct c x xt t a au x t ct x t a ⎧-->⎪⎪=⎨⎪≤⎪⎩4.用Laplace 变换求解⎪⎩⎪⎨⎧===+∞=>+∞<<=0)0,(,0)0,(0),(),(),0(0,0,2x u x u t u t f t u t x u a u t x xx tt 解:对自变量t 取Laplace 变换22220(0,)(),(,)0x d u s u a dx u s f s us ⎧-=⎪⎨⎪=+∞=⎩微分方程的解为x a sx a sBe Ae u-+=~ 再由(0,)(),(,)0x u s f s us =+∞= 所以()s xa a uef s s-=-由 Laplace 变换的卷积定理,得[][][]()*()()()L g x f x L g x L f x =⋅令xa s es a x g --=)(,对其求逆,得:()()0R e s ,lim (0)0,xxs t s t a a k s kx a t a a a e s s e xs s t a --→⎧->⎪⎡⎤⎡⎤⎪-=--=⎨⎢⎥⎢⎥⎣⎦⎣⎦⎪≤⎪⎩∑,最后定解问题的解是)(*)(x f x g 则最后的解为0()(,)0,xt ax a f d t au x t x t a ττ-⎧->⎪⎪=⎨⎪≤⎪⎩⎰,6.有一根均匀弹性细杆,长为l ,一端固定,另一端受外力sin F A tω=作用.杆的初始位移与速度都为0,求杆的纵向振动规律.解:设Y 与S 分别是细杆杨氏模量与截面积,则定解问题为2,0,0sin (0,)0,(,)(,0)0,(,0)0tt xx xx u a u x l t A t u t u l t SY u x u x ω⎧=<<>⎪⎪==⎨⎪==⎪⎩对自变量t 取Laplace 变换2222220(,)(0,)0,x ld u a s u dx du x s Au s dx SY s ωω=⎧-=⎪⎪⎨⎪==⎪+⎩求解常微分方程可得ss xxaauC eD e-=+代入边界条件可得出2201()l l s s aaC D A C SY s s e e ωω-+=⎧⎪⎪=⎨+⎛⎫⎪+ ⎪⎪⎝⎭⎩所以221()ssx x x x s s a a a a l ls s a a u C e D e e e s s e eωω---⎛⎫=+=- ⎪+⎛⎫⎝⎭+ ⎪⎝⎭对上式取Laplace 逆变换可求得。
数学物理方程全部答案
2.在杆纵向振动时,假设(1)端点固定,(2)端点自由,(3)端点固定在弹性支承上,试分别导出这三种情况下所对应的边界条件。
解:(1)杆的两端被固定在l x x ==,0两点则相应的边界条件为 .0),(,0),0(==t l u t u(2)若l x =为自由端,则杆在l x =的张力xu x E t l T ∂∂=)(),(|l x =等于零,因此相应的边界条件为xu ∂∂|l x ==0同理,若0=x 为自由端,则相应的边界条件为xu∂∂∣00==x(3)若l x =端固定在弹性支承上,而弹性支承固定于某点,且该点离开原来位置的偏移由函数)(t v 给出,则在l x =端支承的伸长为)(),(t v t l u -。
由虎克定律有xuE ∂∂∣)](),([t v t l u k l x --==其中k 为支承的刚度系数。
由此得边界条件)(u x u σ+∂∂∣)(t f l x == 其中Ek =σ 特别地,若支承固定于一定点上,则,0)(=t v 得边界条件)(u xu σ+∂∂∣0==l x 。
同理,若0=x 端固定在弹性支承上,则得边界条件 x uE ∂∂∣)](),0([0t v t u k x -==即 )(u xu σ-∂∂∣).(0t f x -=8.求解波动方程的初值问题⎪⎪⎩⎪⎪⎨⎧=∂∂==∂∂-∂∂==xt u u x t x ut u t t s i n |,0s i n 002222 解:由非齐次方程初值问题解的公式得 τξξτααττd d d t x u t t x t x tx tx ⎰⎰⎰-+--+-+=0)()(s i n 21s i n 21),(=⎰----+---+-td t x t x t x t x 0))](cos())([cos(21)]cos()[cos(21ττττ=⎰-+td t x t x 0)sin(sin sin sin τττ=t t t x t x 0)]sin()cos([sin sin sin τττ-+-+ =x t sin 即 x t t x u sin ),(= 为所求的解。
数学物理方程期末考试题及答案
数学物理方程期末考试题及答案一、选择题(每题2分,共10分)1. 以下哪一项不是数学物理方程的特点?A. 连续性B. 离散性C. 线性D. 非线性答案:B2. 波方程是描述什么的方程?A. 热传导B. 电磁波C. 机械波D. 流体动力学答案:C3. 拉普拉斯方程通常出现在哪种物理现象中?A. 热传导B. 流体流动C. 电磁场D. 弹性力学答案:C4. 以下哪个不是偏微分方程的解的性质?A. 唯一性B. 线性C. 稳定性D. 离散性答案:D5. 波动方程的解通常表示什么?A. 温度分布B. 电荷分布C. 压力分布D. 位移分布答案:D二、填空题(每空2分,共20分)6. 波动方程的基本形式是 _______。
答案:\( \frac{\partial^2 u}{\partial t^2} = c^2 \nabla^2 u \)7. 热传导方程,也称为________方程。
答案:傅里叶8. 拉普拉斯方程 \( \nabla^2 \phi = 0 \) 在静电学中描述的是________。
答案:电势9. 边界条件通常分为________和________。
答案:狄利克雷边界条件;诺伊曼边界条件10. 波动方程的一般解可以表示为________和________的叠加。
答案:基频解;高阶谐波三、简答题(每题10分,共30分)11. 解释什么是边界层的概念,并给出一个实际应用的例子。
答案:边界层是流体力学中的一个概念,指的是流体靠近物体表面处的一层非常薄的流体,其中速度梯度很大。
在边界层内,流体的速度从物体表面的零速度逐渐增加到与外部流体速度相匹配。
一个实际应用的例子是飞机的机翼,边界层的厚度和特性对飞机的升力和阻力有重要影响。
12. 描述什么是格林函数,并解释它在解决偏微分方程中的作用。
答案:格林函数是一种数学工具,用于解决线性偏微分方程。
它是一个特定的函数,当它与方程的算子相乘时,结果是一个狄利克雷问题,其解是原始方程的一个解。
数学物理方程(第三版)第一章练习答案全解
解: k, σ 为正常数
utt − a2uxx + kut = 0, 0 < x < l, t > 0, u|t=0 = φ(x), ut|t=0 = ψ(x), u|x=0 = 0, (ux + σu)|x=l = 0.
齐海涛 (SDU)
数学物理方程
2012-10-3 9 / 67
方程的导出、定解条件
方程的导出、定解条件
.E.xample 1.1
细杆(或弹簧)受某种外界原因而产生纵向振动, 以 u(x, t) 表示静止时在 x 点 处的点在时刻 t 离开原来位置的偏移. 假设振动过程中所发生的张力服从胡克 定律, 试证明 u(x, t) 满足方程
(
)
()
∂ ∂t
ρ(x)
∂u ∂t
=
∂ ∂x
E
∂u ∂x
,
.其中 ρ 为杆的密度, E 为杨氏模量.
齐海涛 (SDU)
数学物理方程
2012-10-3 3 / 67
方程的导出、定解条件
.E.xample 1.1
细杆(或弹簧)受某种外界原因而产生纵向振动, 以 u(x, t) 表示静止时在 x 点 处的点在时刻 t 离开原来位置的偏移. 假设振动过程中所发生的张力服从胡克 定律, 试证明 u(x, t) 满足方程
齐海涛 (SDU)
数学物理方程
2012-10-3 3 / 67
方程的导出、定解条件
u(x + ∆x, t) − ∆x
u(x, t)
=
∆u ∆x
=
∂u ∂x
(x,
t),
∆x → 0.
由 Hooke 定律, B 两端的张力分别为 E(x)ux|x, E(x)ux|x+∆x. B 段的运动方
数学物理方程期末考试试题及答案
数学物理方程期末考试试题及答案一、求解方程(15分)⎧utt -a2uxx=0⎪⎨ux-at=0=ϕ(x)⎪u⎩x+at=0=ψ(x).其中ϕ(0)=ψ(0)。
⎧ξ=x-at解:设⎨则方程变为:η=x+at⎩uξη=0,u=F(x-at)+G(x+at)(8’)由边值条件可得:F(0)+G(2x)=ϕ(x),F(2x)+G(0)=ψ(x)由ϕ(0)=ψ(0)即得:u(x,t)=ϕ(x+at x-at)+ψ()-ϕ(0)。
22二、利用变量分离法求解方程。
(15分)⎧utt -a2uxx=0,(x,t)∈Q,⎪⎨ux=0=ux=l=0,t≥0,⎪u=ϕ(x),ut t=0=ψ(x)⎩t=0其中0≤x≤l。
a>0为常数解:设u=X(x)T(t)代于方程得:X''+λX=0,T''+λa2T=0(8’)X=C1cosλx+C2sinλx,T=C1cosλat+C2sinλat由边值条件得:C 1=0,λ=(∞n π2)ln πx lu =∑(B n cos λat +A n sin λat )sin n =1B n =2l n πx 2l n πx ,ϕ(x )sin dx A =ψ(x )sin dx n ⎰⎰00l l an πl2三.证明方程u t -a u xx -cu =0(c ≥0)具有狄利克雷边界条件的初边值问题解的唯一性与稳定性. (15分)证明:设v =e -ct u 代入方程:⎧v t-a 2v xx =0⎪⎨v t =0=ϕ(x )⎪v (0,t )=g (t ),v (l ,t )=g (t ).12⎩设v 1,v 2都是方程的解设v =v 1-v 2代入方程得:⎧v t-a 2v xx =0⎪⎨v t =0=0⎪v (0,t )=,v (l ,t )=0⎩由极值原理得v =0唯一性得证。
(8’)由v 1-v 2≤v 1-v 2得证。
τ≤ε,稳定性得证由v =e -ct u 知u 的唯一性稳定性四.求解二维调和方程在半平面上的狄利克雷问题(15分).∆u =u xx +u yy +u zz=0,z >0,u z =0=f (x ).解:设p (ξ,η,ζ)是上半平面内一点,在该点放置单位点电荷,其对称点p (ξ,η,-ς)格林函数:G (x ,y ,ξ,η)=-14π14π1(x -ξ)+(y -η)+(z -ς)1(x -ξ)+(y -η)+(z +ς)222222+∂G∂G=-∂n∂z z=0=ς2π[(x-ξ)+(y-η)+ς]2223/2方程的解:u(ξ,η)=ς2πϕ(x,y)⎰[(x-ξ)2+(y-η)2+ς2]3/2dx R2五、证明下列初边值问题解的唯一性.(20分)u utt-a2(uxx+uyy)=f(x,y,t) t=0=ϕ(x,y),=ψ(x,y),ut t=0uΓ=g(x,y,t).其中t>0,(x,y)∈Ω,Γ为Ω的边界.解:设u1,u2都是方程的解设u=u1-u2代入方程得:u tt -a(uxx+uyy)=0u u t t=02 =0=0 t=0uΓ=0.设E(t)=12222[u+a(u+u]dxdy t x y⎰⎰2ΩdE(t)=2⎰⎰[ut utt+a2(uxuxt+uyuyt)]dxdydtΩ=2[ut [utt-a(uxx+uyy)]dxdyΩ⎰⎰2=0(10’)E(t)=E(0)=0,u=C,由边值条件得:u=0。
数学物理方程 习题 答案
数学物理方程习题答案E=mc^2: Exploring the Relationship Between Energy and MassIntroductionThe famous equation E=mc^2, proposed by Albert Einstein in 1905, is one of the most well-known and important equations in the field of physics. It states that energy (E) is equal to mass (m) times the speed of light (c) squared. This equation has had a profound impact on our understanding of the relationship between energy and mass, and has led to numerous advancements in science and technology.Exploring the EquationTo understand the significance of E=mc^2, it is important to break down the components of the equation. The letter E represents energy, which is a fundamental concept in physics. Energy can exist in various forms, such as kinetic energy, potential energy, and thermal energy. The letter m represents mass, which is a measure of the amount of matter in an object. Finally, the letter c represents the speed of light, which is a constant value in the universe, approximately 3.00 x 10^8 meters per second.The equation E=mc^2 suggests that mass and energy are interchangeable, and that a small amount of mass can be converted into a large amount of energy. This concept is central to the field of nuclear physics, where the conversion of mass into energy occurs in processes such as nuclear fission and fusion. The equation also has implications for understanding the behavior of particles athigh speeds, as it demonstrates that as an object approaches the speed of light, its energy and mass become increasingly interrelated.Applications and ImplicationsThe implications of E=mc^2 extend beyond the realm of theoretical physics, and have had a significant impact on technology and society. One of the most well-known applications of this equation is in the development of nuclear power and weapons. The process of nuclear fission, which involves the splitting of atomic nuclei, releases a tremendous amount of energy in accordance with E=mc^2. This has led to the use of nuclear reactors for generating electricity, as well as the development of atomic bombs during World War II.Furthermore, the equation has also played a crucial role in the development of medical imaging technologies such as positron emission tomography (PET) scans. These imaging techniques rely on the conversion of mass into energy, as they detect the annihilation of positrons (antimatter particles) and electrons, which results in the release of gamma rays. By understanding the relationship between mass and energy, scientists have been able to develop non-invasive methods for visualizing internal organs and detecting diseases.ConclusionIn conclusion, the equation E=mc^2 represents a fundamental principle in physics that has far-reaching implications for our understanding of the universe. It demonstrates the interplay between mass and energy, and has led to numerous technological advancements and scientific discoveries. As wecontinue to explore the implications of this equation, we gain a deeper understanding of the fundamental forces that govern the behavior of the universe.。
数学物理方程第三版答案
数学物理方程第三版答案1. 动量守恒定律:$$ P = mv $$其中,$P$为物体的动量,$m$为物体的质量,$v$为物体的速度。
2. 牛顿第二定律:$$ F = ma $$其中,$F$为物体受力,$m$为物体的质量,$a$为物体的加速度。
3. 热力学第二定律:$$ \Delta S \geq \frac{Q}{T} $$其中,$\Delta S$ 是系统的熵变化,$Q$ 是系统吸收的热量,$T$ 是系统的温度。
4. 拉普拉斯方程:$$ \nabla^2\phi=0 $$其中,$\phi$ 是空间的潜在电势,$\nabla^2$ 表示拉普拉斯算子,即 Laplace operator。
5.理想气体状态方程:$$PV=nRT$$其中,$P$是压强,$V$是体积,$n$是物质量的数量,$R$是气体的其中一特定常数,$T$是温度。
6. 位能定律:$$ E_p = \frac{mv^2}{2} $$其中,$E_p$为位能,$m$为物体的质量,$v$为物体的速度。
7. 伽马射线波动方程:$$ \frac{\partial^2 \Psi}{\partial x^2} + \frac{\partial^2 \Psi}{\partial y^2} + \frac{\partial^2\Psi}{\partial z^2} = \frac{1}{c^2} \frac{\partial^2\Psi}{\partial t^2}$$其中,$\Psi$ 是伽马射线的波函数,$c$ 是光速。
8. 洛伦兹力方程:$$ F = \frac{d^2x}{dt^2} = - \frac{kx}{m} $$其中,$F$ 为力,$x$ 为物体位移,$k$ 为弹性力系数,$m$ 为物体的质量,$dt$ 为时间差。
数学物理方程课后参考答案第四章
第四章 二阶线性偏微分方程的分类与总结§1 二阶方程的分类1. 证明两个自变量的二阶线性方程经过可逆变换后它的类型不会改变,也就是说,经可逆变换后2211212a a a -=∆的符号不变。
证:因两个自变量的二阶线性方程一般形式为fcu u b u b u a u a u a y x yy xy xx =+++++212212112经可逆变换 ⎩⎨⎧==),(),(y x y x ηηξξ 0),(),(≠y x D D ηξ化为 f u c u b u a u a u a =++++ηηηξηξξ22212112其中 ⎪⎪⎩⎪⎪⎨⎧++=+++=++=22212211222212111222212211112)(2y y x x y y x y y x x x y y x x a a a a a a a a a a a a ηηηηηξηξηξηξξξξξ所以 y x y x y x y x x y yx a a a aa a aηηξξηηξξηξηξ2211112222122221112222)(+-+=-=∆22221112222222211),(),())(()(⎥⎦⎤⎢⎣⎡∆=--=+-y x D D a a aa a x y y x y x y x ηξηξηξηξξη因0),(),(2>⎥⎦⎤⎢⎣⎡y x D D ηξ,故∆与∆同号,即类型不变。
2. 判定下述方程的类型(1)022=-yy xx u y u x (2)0)(2=++yy xx u y x u (3)0=+yyxx xyuu(4))010001(sgn 0sgn 2sgn ⎪⎩⎪⎨⎧<-=>==++x x x x xuu yu yyxy xx(5) 0424=+++-zz yy xz xy xx u u u u u 解:(1)022=-yy xx u y u x因 022>=∆y x 当0,0≠≠y x 时0,0=>∆x 或0=y 时0=∆。
数学物理方程第三版答案谷超豪
数学物理方程第三版答案谷超豪【篇一:数学物理方程_答案_谷超豪】/p> 1 方程的导出。
定解条件1.细杆(或弹簧)受某种外界原因而产生纵向振动,以u(x,t)表示静止时在x点处的点在时刻t离开原来位置的偏移,假设振动过程发生的张力服从虎克定律,试证明u(x,t)满足方程???u????u????x????e? ?t??t??x??x?其中?为杆的密度,e为杨氏模量。
证:在杆上任取一段,其中两端于静止时的坐标分别为 x与x??x。
现在计算这段杆在时刻t的相对伸长。
在时刻t这段杆两端的坐标分别为:x?u(x,t);x??x?u(x??x,t)其相对伸长等于令?x?[x??x?u(x??x,t)]?[x?u(x,t)]??x?x?ux(x???x,t),取极限得在点x的相对伸长为ux(x,t)。
由虎克定律,张力t(x,t)等于t(x,t)?e(x)ux(x,t)其中e(x)是在点x的杨氏模量。
设杆的横截面面积为s(x),则作用在杆段(x,x??x)两端的力分别为e(x)s(x)ux(x,t);e(x??x)s(x??x)ux(x??x,t).于是得运动方程 ?(x)s(x)??x?utt(x,t)?esu利用微分中值定理,消去?x,再令?x?0得??(x)s(x)u?(esux)?x若s(x)?常量,则得?u?t22x(x??x)|x??x?esux(x)|x?(x)即得所证。
=(e(x)?u?x)2.在杆纵向振动时,假设(1)端点固定,(2)端点自由,(3)端点固定在弹性支承上,试分别导出这三种情况下所对应的边界条件。
解:(1)杆的两端被固定在x?0,x?l两点则相应的边界条件为u(0,t)?0,u(l,t)?0.(2)若x?l为自由端,则杆在x?l的张力t(l,t)?e(x)的边界条件为?u?x?u?x|x?l等于零,因此相应|x?l=0?u同理,若x?0为自由端,则相应的边界条件为?x(3)若x?l端固定在弹性支承上,而弹性支承固定于某点,且该点离开原来位置的∣x?0?0偏移由函数v(t)给出,则在x?l端支承的伸长为u(l,t)?v(t)。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
( 6 cos θ + 1) sin θ cos ϕ
解:原式 = 6 cos θ sin θ cos ϕ + sin θ cos ϕ
3 P21 ( cos θ ) = sin 2θ 2 两项都有: m = 1 ,查表得: 1 P 1 ( cos θ ) = sin θ
1 原式化为: 2 P21 ( cos θ ) cos ϕ + P 1 ( cos θ ) cos ϕ
⎛ m = 0,1, 2," l ⎞ ⎜ ⎟的 ⎝ l = 0,1, 2,3," ⎠
第一项: m = 0 ;第二项: m = 1
3 ⎧ 1 ⎪ P2 ( cos θ ) = 2 sin 2θ ⎪ 1 ⎪ 查表得: ⎨ P20 ( cos θ ) = P2 ( cos θ ) = ( 3cos 2θ + 1) 4 ⎪ 0 ⎪ P0 ( cos θ ) = P0 ( cos θ ) = 1 ⎪ ⎩
∞ n 2π 2 a 2t l2
,
利用初始条件: u t =0 = bx ( l − x ) l 2 ,
⇒ ∑ An sin
n =1
∞
nπ x = bx ( l − x ) l 2 , l
An 是 Fourier 级数的系数。
nπ x 2 l bx ( l − x ) dx sin 2 ∫ l 0 l l nπ x nπ x 2 l blx 2 l bx 2 = ∫ 2 sin dx − ∫ 2 sin dx 0 0 l l l l l l An =
⎧sin mϕ ⎫ 1. 将下列函数展开为球函数 Yl m (θ , ϕ ) = Pl m ( cos θ ) ⎨ ⎬ ⎩cos mϕ ⎭ 形式。 (1) sin θ ( sin θ + cos θ sin ϕ ) 解:
sin θ ( sin θ + cos θ sin ϕ ) = sin 2 θ + sin θ cos θ sin ϕ
0 − 4b ⎡cos ( nπ ) − cos 0 ⎦ ⎤, ⎣ n3π 3 = − , 2 n odd = n = even , 0
{
An =
8b , n = 2 p + 1, n3π 3
∞
p = 0,1, 2,"
求得的解为 u ( x, t ) = ∑
p =0
8b
3 ( 2 p + 1) π 3
1 1 原式 = P 1 ( cos θ ) cos ϕ + P 2 ( cos θ ) cos ϕ
3. 如图所示,长为 l 的弦,两端固定,弦中张力为 T ,在弦的中间点以横向力 F0 把弦拉开,然后突然撤除这力,求解弦的振动。
解:首先求整根弦的初始位移。假设 l 2 点弦的初始位移为 h ,由受力平衡得:
⇒ 本征值 λ =
nπ x n 2π 2 ,本征函数族为 sin , n = 1, 2,3," 2 l l
将所有线性独立的特解叠加,得到方程在边界条件下的通解形式为
∞ nπ at nπ at ⎞ nπ x ⎛ u = ∑ ⎜ An cos + Bn sin , ⎟ sin l l ⎠ l n =1 ⎝
将初始条件代入通解,
ut
t =0
= 0 ⇒ Bn = 0 ,
l ⎧ 2h x, 0 ≤ x ≤ , ⎪ nπ x ⎪ l 2 An sin 将等号右边的函数形式设为 ϕ ( x ) , =⎨ ∑ l l n =1 ⎪ − 2h ( x − l ) , ≤ x ≤ l, ⎪ 2 ⎩ l
∞
将 ϕ ( x ) 展开为 Fourier 级数, An 为 Fourier 级数的系数。
⎧sin mϕ ⎫ 2. 将下列函数展开为球函数 Yl m (θ , ϕ ) = Pl m ( cos θ ) ⎨ ⎬ ⎩cos mϕ ⎭
形式。
⎛ m = 0,1, 2," l ⎞ ⎜ ⎟的 ⎝ l = 0,1, 2,3," ⎠
(1) 3 ( sin 2θ sin ϕ − cos 2θ sin 2ϕ − cos 2θ cos 2ϕ + sin 2ϕ + cos 2ϕ ) − 1 解:含 sin 2ϕ 的项为 ( −3cos 2θ + 3) sin 2ϕ , m = 2, 查表得: P22 ( cos θ ) =
这一项化为 3sin 2θ sin ϕ = 2 P21 ( cos θ ) sin ϕ 常数项 −1 = − P0 ( cos θ ) 原式化为 2 P22 ( cos θ ) sin 2ϕ + 2 P22 ( cos θ ) cos 2ϕ + 2 P21 ( cos θ ) sin ϕ − P0 ( cos θ )
l
du = 2 xdx, 第二项作分部积分: nπ x l nπ x dv = sin dx, v=− cos , l nπ l
u = x2 ,
2b ⎛ l ⎞ 2 nπ x 2b l l nπ x − 3 ⎜− − 3 ∫ 2 x cos dx ⎟ x cos l ⎝ nπ ⎠ l 0 l 0 nπ l
,
n = 2 p + 1, p = 0,1, 2,"
得: u = ∑
p =0
∞
8h
( 2 p + 1) π
2
2
( −1)
p
cos
( 2 p + 1) π at sin ( 2 p + 1) π x
l l
4. 求解细杆导热问题,杆长 l ,两端保持为零度,初始温度分布
u t =0 = bx ( l − x ) l 2
2 p +1 π a t 2 p + 1) π x − ( l) ( sin e
2 2 2 2
l
5. 在球坐标系下将三维波动方程 utt − a 2∇ 2u = 0 分离变量。其中,拉普拉斯算符 在球坐标系下的形式为
∇ 2u = 1 ∂ ⎛ 2 ∂u ⎞ 1 ∂ ⎛ ∂u ⎞ 1 ∂ 2u θ r + sin + ⎜ ⎟ ⎜ ⎟ r 2 ∂r ⎝ ∂r ⎠ r 2 sin θ ∂θ ⎝ ∂θ ⎠ r 2 sin 2 θ ∂φ 2
1 1 1 1 ⎧1 ⎫ − cos 2θ = − ⎨ ⎡ 4 P20 ( cos θ ) − 1⎤ ⎬ ⎣ ⎦ 2 2 2 2 ⎩3 ⎭ 2 2 = P00 ( cos θ ) − P20 ( cos θ ) 3 3 1 sin θ cos θ = P21 ( cos θ ) 3 sin 2 θ =
u = x,
du = dx, nπ x dx, l
l
第一项作分部积分:
dv = sin
v=−
l nπ x cos , nπ l
2b l nπ x 2b l 2 nπ x x cos − 2 + 2 2 2 sin l nπ l 0 l nπ l 0
− 2b cos nπ nπ 0
解:根据题意列出泛定方程和定解条件 k ⎧ 2 a2 = , ⎪ut − a u xx = 0, cρ ⎪ ⎪ ⎨u x =0 = u x =l = 0, ⎪ 2 ⎪u t =0 = bx ( l − x ) l , ⎪ ⎩ 分离变量,设 u ( x, t ) = X ( x ) T ( t )
T '(t ) X '' ( x ) = = −λ , 2 a T (t ) X ( x ) ⎧ ⎪ X '' ( x ) + λ X ( x ) = 0, 构成本征值问题 ⎨ = = X 0 X l 0, ( ) ( ) ⎪ ⎩ ⇒ 本征值 λ =
2 l nπ x dx ϕ ( x ) sin ∫ 0 l l 2 l 2h nπ x 2 l ⎛ 2h ⎞ nπ x = ∫2 x sin dx + ∫ l ⎜ − ⎟ ( x − l ) sin dx l 0 l l l 2⎝ l ⎠ l An =
l l ⎛ ⎞ 2 2h ⎜ l nπ x 2 l2 nπ x 2 ⎟ = −x + cos sin l l ⎜ l 0 n 2π 2 l 0⎟ ⎜ nπ ⎟ ⎝ ⎠ l ⎞ ⎛ l 2 2h ⎜ l nπ x l2 nπ x ⎟ − −x + cos sin l l ⎜ l l n 2π 2 l l⎟ ⎜ nπ ⎟ 2 ⎝ 2 ⎠
根据物理问题列出关于弦的横位移 u 的泛定方程和定解条件得
T ⎧ 2 a2 = , ⎪utt − a u xx = 0, ρ ⎪ ⎪u x =0 = u x =l = 0, ⎪ ⎪ l ⎧ 2h x, 0 ≤ x ≤ , ⎨ ⎪ l 2 ⎪u = ⎪ ⎪ t =0 ⎨ l ⎪ − 2h ( x − l ) , ≤ x ≤ l, ⎪ ⎪ 2 ⎩ l ⎪ ⎪ ⎩ut t =0 = 0
(2) sin θ cos ϕ
1 解: m = 1 ,查表得: P 1 ( cos θ ) = sin θ 1 原式 = P 1 ( cos θ ) cos ϕ
(3) (1 + 3cos θ ) sin θ cos ϕ
3 P21 ( cos θ ) = sin 2θ 2 解: m = 1 ,查表得: 1 P 1 ( cos θ ) = sin θ
2
n 2π 2 nπ x ,本征函数族为 sin , n = 1, 2,3," 2 t ) + λ a T ( t ) = 0 的解为 Ce
,
可得满足方程和边界条件的温度 u ( x, t ) 的通解形式为:
nπ x − u = ∑ An sin e l n =1
用分离变量法,设方程的解具有形式 u ( x, t ) = X ( x ) T ( t ) ,代入泛定方程得:
X '' ( x ) T '' ( t ) = = −λ , −λ 为一个待定的常数 X ( x ) a 2T ( t ) ⎧ ⎪ X '' ( x ) + λ X ( x ) = 0 构成本征值问题 ⎨ X 0 = X l = 0 ( ) ( ) ⎪ ⎩