1、数学物理方程答案第一章1 方程的导出。定解条件1细杆（或弹簧）受某种外界原因而产生纵向振动，以 u(x,t)表示静止时在 x 点处的点在时刻 t 离开原来位置的偏移，假设振动过程发生的张力服从虎克定律，试证明 满足方程),(tuxEtuxt其中 为杆的密度， 为杨氏模量。E证：在杆上任取一段，其中两端于静止时的坐标分别为 与 。现在计算这段杆在时刻 的xt相对伸长。在时刻 这段杆两端的坐标分别为：t ),();,(txutxu其相对伸长等于 ),(, txux令 ，取极限得在点 的相对伸长为 。由虎克定律，张力 等于0xxx),(tT)(,(uEtT其中 是在点 的杨氏模量。)(xE设杆的横截
2、面面积为 则作用在杆段 两端的力分别为),(xS),(xxuS)( xut);,.,t于是得运动方程 txs)(ES,( xxESu|)(|)(利用微分中值定理，消去 ，再令 得0tu)(x()若 常量，则得)(xs=2)(tx)(xE即得所证。2在杆纵向振动时，假设(1)端点固定，(2)端点自由， (3)端点固定在弹性支承上，试分别导出这三种情况下所对应的边界条件。解：(1)杆的两端被固定在 两点则相应的边界条件为lx,0.)(,)0(tlut数学物理方程答案(2)若 为自由端，则杆在 的张力 | 等于零，因此相应的边界lxlxxuEtlT)(,l条件为 | =0 ulx同理，若 为自由端，
3、则相应的边界条件为 0x0(3)若 端固定在弹性支承上，而弹性支承固定于某点，且该点离开原来位置的偏移由函数l给出，则在 端支承的伸长为 。由虎克定律有)(tvx)(,tvluEklx其中 为支承的刚度系数。由此得边界条件k 其中)(u)(tflxEk特别地，若支承固定于一定点上，则 得边界条件,0tv 。)(xlx同理，若 端固定在弹性支承上，则得边界条件0xuE)(,0tvkx即 )(.f3. 试证：圆锥形枢轴的纵振动方程为 22)1()1(tuhxuhxE其中 为圆锥的高(如图 1)h证：如图，不妨设枢轴底面的半径为 1，则点处截面的半径 为：lhx1所以截面积 。利用第 1 题，得2)
4、()xs)(22xuhExtuhx若 为常量，则得Ex)( 22)1()1(tuhxuhx4. 绝对柔软逐条而均匀的弦线有一端固定，在它本身重力作用下，此线处于铅垂平衡位置，试导出此线的微小横振动方程。解：如图 2，设弦长为 ，弦的线密度为 ，则 点处的张力 为lx)(xT数学物理方程答案)()xlgT且 的方向总是沿着弦在 点处的切线方向。仍以 表示弦上各点在时刻 沿垂直于 轴方向)(xT ),(txutx的位移，取弦段 则弦段两端张力在 轴方向的投影分别为),(x )(sin)();(sinxxlglg 其中 表示 方向与 轴的夹角)(x)(Tx又 .siut于是得运动方程 xltx)(2
5、 xulgg利用微分中值定理，消去 ，再令 得0。)(2xulgtu5. 验证 在锥 0 中都满足波动方程221),(ytyx2yxt证：函数 在锥 0 内对变量22yuxt 221),(yxtx2yxt有yx,二阶连续偏导数。且 txtu232)(252 )()( tytytt )2()(3xtxtytu2)(252232 xyxtxt 5tyt数学物理方程答案同理 22522 yxtyxtyu所以 .22 tuttx 即得所证。6. 在单性杆纵振动时,若考虑摩阻的影响,并设摩阻力密度涵数(即单位质量所受的摩阻力)与杆件在该点的速度大小成正比(比例系数设为 b), 但方向相反,试导出这时位移
6、函数所满足的微分方程 .解: 利用第 1 题的推导,由题意知此时尚须考虑杆段 上所受的摩阻力.由题设,单位质量所受x摩阻力为 ,故 上所受摩阻力为tubxtusp运动方程为: tuxsbxuEStEStxsx2利用微分中值定理，消去 ,再令 得0 .2 tuxsbxuStuxs 若 常数，则得)(xstuxbExtu2若 则 得 方 程令也 是 常 量是 常 量 ,., 2ax .22xutbtu2 达朗贝尔公式、 波的传抪1. 证明方程 常 数011222 htuxaxuhx的通解可以写成xtGtFu数学物理方程答案其中 F,G 为任意的单变量可微函数,并由此求解它的初值问题: .,:0xt
7、ut解：令 则vuxhxvuhxhxvu2,)()()()()( 222 xvuhxhxuxvuxh 又 2t代入原方程，得 221tvxhaxvh即 22t由波动方程通解表达式得 atxGtFtxv,所以 hau为原方程的通解。由初始条件得 )1(1xGFxha/所以 )2(10 cdhxx由 两式解出)2(,122121cdhaxhxFxo数学物理方程答案22121cdhaxhxGxo所以 )()()()(),( atxttttu+ atxh()(21.d即为初值问题的解散。问初始条件 与 满足怎样的条件时，齐次波动方程初值问题的解仅由右传播波组成？)(x解：波动方程的通解为u=F(x-a
8、t)+G(x+at)其中 F， G 由初始条件 与 决定。初值问题的解仅由右传播组成，必须且只须对)(x于任何 有 G(x+at) 常数.tx,即对任何 x, G(x) C 0又 G（x）= xaCda02)(21)(所以 应满足)(,（常数）xx01)(或 (x)+ =01a3.利用传播波法，求解波动方程的特征问题（又称古尔沙问题）).(022xutatx)0(解：u(x,t)=F(x-at)+G(x+at)令 x-at=0 得 =F（0）+G（2x）)(令 x+at=0 得 =F（2x）+G(0)x所以 F(x)= -G(0).)2(数学物理方程答案G（x）= -F(0).)2(且 F（
9、0）+G(0)= ).0(所以 u(x,t)= + -()attx.即为古尔沙问题的解。4对非齐次波动方程的初值问题 )(),(,0,022 xxtut tfxau证明：（1） 如果初始条件在 x 轴的区间x ,x 上发生变化，那末对应的解在区间 ，12 1x的影响区域以外不发生变化；2x（2） 在 x 轴区间 上所给的初始条件唯一地确定区间 的决定区2,1 21,x域中解的数值。证：（1） 非齐次方程初值问题的解为u(x,t)= atxtxat2)()(2d)(+ ttxadf0)(.,1当初始条件发生变化时，仅仅引起以上表达式的前两项发生变化，即仅仅影晌到相应齐次方程初值的解。当 在 上发
10、生变化，若对任何 t0,有 x+atx ,则区间x-at,x+at 整个落在),(x2,1x12区间 之外，由解的表达式知 u(x,t)不发生变化，即对 t0,当 xx +at,也就是（x,t）落在2,1区间 的影响域x)0(2taxat之外，解 u(x,t)不发生变化。 （1）得证。(2). 区间 的决定区域为 2, atxtxt21,在其中任给（x,t）,则21atxtx故区间x-at,x+at完全落在区间 中。因此 上所给的初绐1,21,x数学物理方程答案条件 代入达朗贝尔公式唯一地确定出 u(x,t)的数值。)(,x5. 若电报方程 GRuLCuttx 具体形如为 常 数GRLC,at
11、xftx的解（称为阻碍尼波） ，问此时 之间应成立什么关系？GRLC,解 txftxu,aatxfttxftt atxftut 22代入方程，得 012 atxftGRttLGCRtL tfLCaaxfa由于 是任意函数，故 的系数必需恒为零。即ff,0201tGRtLCRtLa于是得 2atu所以tLGCRaec20代入以上方程组中最后一个方程，得 0242 GRLa又 CRCLa21,得即 02G数学物理方程答案最后得到 RGLC6利用波的反射法求解一端固定并伸长到无穷远处的弦振动问题 0, 0,02tuxxatt 解：满足方程及初始条件的解，由达朗贝尔公式给出：。atxtdtxattx
12、2121,由题意知 仅在 上给出，为利用达朗贝尔解，必须将 开拓到,0 x,上，为此利用边值条件，得0x。attdat21因此对任何 必须有t tt0attd即 必须接奇函数开拓到 上，记开拓后的函数为 ；x, xx,0,0, xx所以 atxtdtattxu 2121,。 0,2121 ,xatdaxtatx ttttxtatx7求方程 形如 的解（称为球面波）其中222zuyxtutrf,。2zyxr解： trfu,数学物理方程答案 xrurxu 32221r3222yuryu)1(3222rzzz代入原方程，得)(32222 rzyxruatu即 )(22t令 ，则vru222, rvu
13、rvurt ，代入方程，得 v 满足22rvat故得通解 )()(),(atrGtFv所以 1ru8求解波动方程的初值问题xtuttt sin|,02解：由非齐次方程初值问题解的公式得ddtxutxttxt0)()(si21sin21),(=数学物理方程答案t dtxtxxtx0 )(cos)(cos21)cos()cs(21 =tdt0)in(ini= ttxt 0)si(cosisi =t即 为所求的解。xtuin),(9求解波动方程的初值问题。 20021|,| )(xutattxt解： tatxatxt ddt 0)()(221),( atxt trcgrcgd)(12 t taxta
14、xt dd0)(20)()(21(1= txtax0 22)()(2= xat atxdudut )1()1(22= xat atxxatz222= 222 )(1ln4)()( txrcgtrcg+ atxtatxtatt = )()(21)()(21 trcgrcgx+ 2)(ln4atxat数学物理方程答案所以 )(1ln2)( )2()()41, 23 atxatrcgxxarctg tttttu 3 混合问题的分离变量法1. 用分离变量法求下列问题的解：(1) 0),(,0 )0()13sin22tlut lxxxuatott解：边界条件齐次的且是第一类的，令 )(),(tTxXt得
15、固有函数 ，且lnnsi，tlaBtaAtTnnco)()2,1(n于是 1 siics(,nnxlntltltxu 今由始值确定常数 及 ，由始值得1sin3sinxlAlx1i)(nnlBlal所以 当,3A,03ln xdlnxa0si)(2 xlnllll cossico22数学物理方程答案 )1(4cos2sin2 303 nlanxlnlxll 因此所求解为1443 sii)(si3co),( nnxltlalxltlatxu(2) 0),(,)0,(22xtuxlhutlt解：边界条件齐次的，令)(),(tTxXt得： (1) 0)(,0)(lX及 。)2(2aT求问题(1)的非
16、平凡解，分以下三种情形讨论。时，方程的通解为10xxeCX21)(由 得)(21c由 得0l 02lleeC解以上方程组，得 ， ，故 时得不到非零解。10时，方程的通解为2xcX21)(由边值 得 ，再由 得 ，仍得不到非零解。0)(X1cl0时，方程的通解为3xcxcxsinos)(21由 得 ，再由 得 0)(X1c0lX数学物理方程答案0cos2l为了使 ，必须 ，于是02c21ln)2,10(n且相应地得到 xlxXnsi)( ,将 代入方程(2)，解得talnBtalnAtTn 21si21co)( )2,10(n于是 0 si)is,nnxltltltxu 再由始值得 021si
17、n21nnxlBalAxlh容易验证 构成区间 上的正交函数系：xlsi ),( ,0l nmlxdlnlml 当当212siin0利用 正交性，得xl21sixdlnxlhAln21si0lxnnllll 022 1si)1(cos)1( nnh)()(820B所以 022 21sin1cos)(18),(nn xltalhtxu 数学物理方程答案2。设弹簧一端固定，一端在外力作用下作周期振动，此时定解问题归结为求解此问题。0),()0,( sin),(22xtutAtlat 解：边值条件是非齐次的，首先将边值条件齐次化，取 ，则 满足txltUsin),(),(txU，0),(tUtAtl
18、sin),(令 代入原定解问题，则 满足),(xvtxu),(txv)1()0,()0,(sin222xlAxtvxvlt ta满足第一类齐次边界条件，其相应固有函数为 ， ),(t xlnXnsi)()2,10(故设 2sin)(),(1nxltTtxv将方程中非齐次项 及初始条件中 按 展成级数，得tlAi2 lAxlnsi12sin)(sinxltftxl 其中 ln dltAtf02i)(lxnlxlntl 022 sicossin lAtAni)1(xlnsi1其中 nln AdlA)1(2i20 数学物理方程答案将(2)代入问题(1)，得 满足)(tTn nnnATttlat )1
19、(2)0(,)0( si2解方程，得通解 212)(si)sicos)( lattlaBtlaAt nnnn由始值，得 0n2231)(1)(2)1( lanAlanAAaBn 所以 12si)(),(n tlltxvxlntnlaAi1)(2xlnttlalnl sisi)(212 因此所求解为12)(si),(nlalAtxltuxlttlasii3用分离变量法求下面问题的解0|022lxxttubshxu解：边界条件是齐次的，相应的固有函数为),21(sin)(xlxXn设 1i),nltTtu数学物理方程答案将非次项 按 展开级数，得bshxinxl1sin)(nxltf其中 shlb
20、nldlhlbtfn 2)1(i2)( 20将 代入原定解问题，得 满足1si)(),(nxltTtxu)(tTn 0)(,)0( 2)1(2nn nshlbtlat 方程的通解为 shllnbaltlnBtlaAtT nn 122)()(sicos)( 由 ，得：0n shlbn 122)(由 ，得)(所以 )cos1()()1(22 tlanllnatTnn 所求解为1212 si)()(),(nn xltllshlbtxu4用分离变量法求下面问题的解：0|,| )(022ttlxxuhubatbt解：方程和边界条件都是齐次的。令)(),(tTxXt代入方程及边界条件，得 ab“2“0)(
21、lX数学物理方程答案由此得边值问题 0)(“lX因此得固有值 ，相应的固有函数为2ln,1,si)(xlxn又 满足方程)(tT02“Tab将 代入，相应的 记作 ，得 满足n)(t)(tn)(tn022“ lbTn一般言之， 很小，即阻尼很小，故通常有,21,2la故得通解 )sinco()( tBtAetTnbtn 其中 2la所以 xlntBtAetxunnbt si)icos(),(1再由始值，得 xlnBbAlxlhnnsi)(0i1所以102)(2sinlnhxdhA1)(nnnbB所求解为数学物理方程答案 .sin)i(cos)12),( xltbtnehtxunbt 4 高维波
22、动方程的柯西问题1 利用泊松公式求解波动方程)(2zyxt uau的柯西问题 023t解：泊松公式dsradsratuStMStM4141现 zyx23,0且 0|sin),(atrs drdrMat其中 )cos,sin,coi),( rzyx()isn()sin23r3223 csicsci rxrxzyxosis22yzrcsnincoi23ry计算 0si),(dr)(4 )cos(2)(in2302 03zyxr zyxr 02 0202cossin3sincosi drxdr 02 02032iii3xr数学物理方程答案 2003 sin41cos1xr driin43023320
23、04cossi xdr 020202 siniini dyzrryz zrzrd320033 23022 4sin41cos1 iiisi 00202 isincodrydy 2002302 sincosiniii dyr 020340223sincosiniidrr所以 3144)( 22232ztatxzyat rrdsStMatu(x,y,z)= SatMrt1ztaxtzyxt222331即为所求的解。2 试用降维法导出振动方程的达朗贝尔公式。数学物理方程答案解：三维波动方程的柯西问题),(),(002 zyxuzyxuatt zt 当 u 不依赖于 x,y,即 u=u(z),即得弦振
24、动方程的柯西问题： )(),(002zztt利用泊松公式求解SatMSatMdsrdsrtu4141因只与 z 有关，故SatMdttzdsr 202sin)(cos( dattzdsin)co(20令 ,= atcos+z tin-得 SatMatztddr)(2所以 atzt atztddtzu)(21)(21),( atzttzatz)()()(即为达郎贝尔公式。3. 求解平面波动方程的柯西问题： 0|20ttyxuua解： 由二维波动方程柯西问题的泊松公式得：mat dyxtyx 222,21,数学物理方程答案mat dyx222, 02sin,co1rtartt又 sicosin,c
25、oyxryrx22cosrryxcsin2sicrxno23因为 20020 cos,sin,cosddd.0sin,i 2020320 所以 at rdtyrx022sin,coa atrdyxrt0 02322 3又 at atrd0 02|at att rdrt0 022223 |3023|trta于是 yxayxttyxu 321,2a322即为所求的解。4. 求二维波动方程的轴对称解（即二维波动方程的形如 的解，tru,数学物理方程答案.)2yxr解： 解法一：利用二维波动方程柯西问题的积分表达式,21, 222 222matmat yxdtxu由于 u 是轴对称的 故其始值 , 只
26、是 r 的函数， ,trurut0|，记圆上任一点 的矢径为.,| 220 tayxrmatt 为 圆又 p圆心 其矢径为 记 则由余弦定2),(yM2r22yxs理知， ，其中 为 与 的夹角。选极坐标 。cosrsroMp),(scos22rsr，,于是以上公式可写成sdatrsrtatyxut202co1, stt202由上式右端容易看出，积分结果和 有关，因此所得的解为轴对称解，即),(trat sdtrsrtu022)(co1, + )(cos022datrrat解法二：作变换 , .波动方程化为xiny)1(22ruratu用分离变量法，令 u(r,t)=R(r)T(t).代入方程
27、得数学物理方程答案02“2RrrtaT解得：)()sincos0rJrtaBtAt 令 叠加得duJtBtAtru )(sin)(cos)(),( 00 5.求解下列柯西问题),(),(0022yxrvycvattxt 提示：在三维波动方程中，令 ,tvezuacz解：令 ),(),(tyxvetzyxuacz则 yacztact veveuzz2代入原问题，得 ),(),(002 yxeuyxeuaacztcztt dsdsatzyx MatcMtac SrSr),(4141),( , 222)()(: zySMat 记 为上半球， 为下半球， 为 在 平面上的投影。t atatto,则dy
28、xtads222)()( Mat Mat MatSSSccc dsersersre ),(1),(1), 数学物理方程答案Mat dyxetazc ),()()(222)()( 22yxteMat tazc ),()(222)()( 22 dyxtcheMatcz ),()()(222 rdrrtacchetaz )sin,co(202所以 xrtahetzyxutacz ()21),(02sin,cofrr rdyrxrtacchetaz )sin,co()2102于是 xrtatyxvt ()21),(02xrtacchrdyr t ()21)sin,co02i,即为所求的解。6试用 第七
29、段中的方法导出平面齐次波动方程4),()(2tyxfuauyxt 在齐次初始条件0,0tt数学物理方程答案下的求解公式。解：首先证明齐次化原理：若 是定解问题),(tyxw),(,02fatott的解，则 即为定解问题tdyxwyxu0),(),(0,),()02tt yxtutxfa的解。显然， 0tudtwttyxwttt00),(（ ）.所以0t tu又 dttt022dywudxuyftt 2202202,),(因为 w 满足齐次方程，故 u 满足)(),(222 yxatyft u齐次化原理得证。由齐次方程柯西问题解的泊松公式知 dyxtaftyxwMta)( 222)()(),1)
30、,(所以 rdtyrxftyxut0)(222)(),sin,co1),(即为所求的解。所以 ta rtayrxftyx0)( 22)(),si,c2),( 数学物理方程答案7用降维法来解决上面的问题解：推迟势dvratfatzyxutr),(41),(2其中积分是在以 为中心， 为半径的球体中进行。它是柯西问题,0, ),()(02tt zyxutzyxf的解。对于二维问题 ， 皆与 无关，故fdsrratfatyxutSMr02),(41),(其中 为以 为中心 r 为半径的球面，即 Mrs)0,( 222)(: ryxSr drds222)()( dsratfsratfdsratf Mr
31、MrMr SSS ),(), dyxrtfM222)()(,其中 分别表示 的上半球面与下半球面， 表示 在 平面上的投影。rs,rsMrrso所以 atrM dyxatftyxu0 2222 )()(,1),( drrartfatr 022),sin,co1 数学物理方程答案在最外一层积分中，作变量置换，令 ，即 ，当 时 ，art ),(taadr0rt当 时， ，得atr0t tayxftyxu0)(222)(),sin,co1),( 即为所求，与 6 题结果一致。8 非齐次方程的柯西问题 yzxuttt 200,)(解：由解的公式得 )1(),(4141),( 2 adVrtfadsr
32、atzyxuMtStr计算 MtS zydsr 022 )cos)(sin()cosin( 02 22ici(i rxrztr trdzy sin)osinisnco202,4sid020ci0223,3coin 02osind0,sid02 .0ici所以 MtStyzxtdr324)(4计算 tr tr drrrtdVtf sin)sin(, 2dry02si)sin(数学物理方程答案t dry0sin)(4.343)(21820tytrtt所以 1, ttzxtzyxu)(2t即为所求的解。5 能量不等式，波动方程解的唯一和稳定性1 设受摩擦力作用的固定端点的有界弦振动，满足方程txtc
33、uau2证明其能量是减少的，并由此证明方程 ftxt的混合问题解的唯一性以及关于初始条件及自由项的稳定性。证： 首先证明能量是减少。1能量 lxtduatE02)()l xttdt )( |200 dxuadulttxtl ltxxtl 20 |)(因弦的两端固定， 所以,|,|lxxu0|0ltt于是 dxuadtEtl)(2)(20(xclt )0c数学物理方程答案因此，随着 的增加， 是减少的。t)(tE证明混合问题解的唯一性.2混合问题: )(|),(|002xuxfcattltt 设 是以上问题的解。令 则 满足21,u,210|,|0ttlxxttuca能量 daExlt)()20当 利用初始条件有 由 得,0t ,|tu,0|t|0tx所以 )(E又 是减少的，故当 又由 的表达式知)(tE,)(,tt )(tE,0)(t所以 0)(t由此得 及 于是得到,0tu,x常量再由初始条件 得 因此 即混合问题解的唯一的。,|0t,u,21u3 证明解关于初始条件的稳定性，即对任何 可以找到 只要初始条件之差. 0.,0满足2121,222 |,|,| 111 LLxL则始值 所对应的解 及 所对应的解 之差 满足 ),(1u)(u2|1L
