数学物理方程习题讲义.pdf
- 文档编号:30816748
- 上传时间:2024-01-30
- 格式:PDF
- 页数:49
- 大小:545.25KB
数学物理方程习题讲义.pdf
《数学物理方程习题讲义.pdf》由会员分享,可在线阅读,更多相关《数学物理方程习题讲义.pdf(49页珍藏版)》请在冰豆网上搜索。
数学物理方程习题讲义苏海军山东大学数学学院-2-摘摘摘要要要本讲义所收录习题基于高等教育出版社谷超豪等编著的数学物理方程(第二版)第一章至第四章。
习题解答尚处于完善过程中,作者水平有限,其中难免有疏漏与不足之处,敬请谅解.目录目录第一章波动方程.11.1学习要求.11.2习题选讲.1第二章热传导方程.172.1学习要求.172.2习题选讲.17第三章调和方程.263.1学习要求.263.2习题选讲.26第四章二阶线性偏微分方程的分类与总结.384.1学习要求.384.2习题选讲.38-ii-第一章波动方程第一章波动方程1.1学习要求
(1)理解弦振动方程的物理意义,定解条件的物理意义.
(2)理解波的左右传播,理解依赖区间,决定区域和影响区域的概念,掌握齐次化原理.(3)理解波动方程分离变量法解的物理意义,掌握非齐次边界条件的齐次化方法.(4)理解膜振动方程的物理意义,掌握球平均法和降维法.(5)熟练掌握达郎贝尔公式和分离变量法的推导过程,会应用这两种方法求解定解问题.(6)熟练和非齐次边界条件的齐次化方法.1.2习题选讲1.方程的导出、定解条件1.细杆(或弹簧)受某种外界原因而产生纵向振动,以u(x,t)表示静止时在x处的点在时刻t离开原来位置的偏移.假设振动过程中所发生的张力服从胡克定律,试证明u(x,t)满足方程t?
(x)ut=x?
Eux,其中为杆的密度,E为杨氏模量.证明:
如图建立坐标系,选取杆上一段微元(x,x+x),则微元两端的相对伸长分别为ux(x,t)和ux(x+x,t).假设杆的横截面面积为S,则微元两端所受拉力分别为E(x)ux(x,t)S(x)和E(x+x)ux(x+x,t)S(x+x).因此所受合力为E(x+x)ux(x+x,t)S(x+x)E(x)ux(x,t)S(x),且正向与坐标轴相同.图图图1-1图示设x为微元重心,则重心处加速度为2ut2(x,t),由牛顿第二定律得,(x)S(x)x2ut2(x,t)=E(x+x)S(x+x)ux(x+x,t)E(x)S(x)ux(x,t)=x?
E(x)S(x)ux(x,t)x-1-1.2习题选讲其中x(x,x+x).约去x并令x0,即得t?
(x)S(x)ut=x?
E(x)S(x)ux当S(x)为常数时,即为t?
(x)ut=x?
Eux,2.在杆纵向振动时,假设
(1)端点固定,
(2)端点自由,(3)端点固定在弹性支承上,试分别导出这三种情况下所对应的边界条件.解:
设杆的两个端点坐标分别为0和l.
(1)端点固定:
此时两个端点无位移,即u(0,t)=u(l,t)=0;
(2)端点自由:
此时两个端点无约束,根据上题,拉力E(x)ux(x,t)S=0,即ux(0,t)=ux(l,t)=0;(3)端点固定在弹性支承上:
此时端点所受外力与弹性支承的变形成比例.若支承的弹性系数为k,则支承对杆的左端点x=0处的作用力为E(0)ux(0,t)S,且其正向与x轴方向相反,因此有E(0)ux(0,t)S=ku(0,t),或写为?
ux+uflflflflx=0=0;其中=k/ES.类似的,对x=l端,有?
ux+uflflflflx=l=0.3.试证:
圆锥形枢轴的纵振动方程为Ex?
1xh2ux=1xh22ut2,其中h为圆锥的高.证明:
此时S(x)=S01xh2,其中S0为圆锥枢轴的底面积.根据第1题的推导,即得所证.图图图1-2图示4.绝对柔软而均匀的弦线有一端固定,在它自身重力的作用下,此线处于铅垂的平衡位置,试导出此线的微小横振动方程.-2-第一章波动方程解:
根据弦的微小横振动方程,有2ut2=x?
T(x)ux其中T(x)为弦的内部张力.在本题中,T(x)=g(lx),故有2ut2=gx?
(lx)ux.5.一柔软均匀的细弦,一端固定,另一端是弹性支承.设该弦在阻力与速度成正比的介质中作微小的横振动,试写出弦的位移所满足的定解问题.解:
此时所受外力为阻力F(x)=kut,因而有T2ut22ux2=kut假设固定端为x=0,有u(0,t)=0;对于弹性支承端x=l,有?
ux+uflflflflx=l=0.6.若F(),G()均为其变元的二次连续可导函数,验证F(xat),G(x+at)均满足弦振动方程(1.11).解:
参见第二节.7.验证u(x,y,t)=1pt2x2y2在锥t2x2y20中满足波动方程2ut2=2ux2+2uy2.解:
显然,ut=t(t2x2y2)3/2,2ut2=3t2(t2x2y2)5/2t2x2y23/2类似的,2ux2=3x2(t2x2y2)5/2+t2x2y23/2,2uy2=3y2(t2x2y2)5/2+t2x2y23/2,代入即得所证.2.达朗贝尔公式、波的传播1.证明方程x?
1xh2ux=1a21xh22ut2,的通解可以写成u(x,t)=F(xat)+G(x+at)hx其中h0为常数,F,G为任意的具有二阶连续导数的单变量函数,并由此求解它的初值问题:
t=0:
v=(hx)(x),vt=(hx)(x)解:
(1)令v(x,t)=(hx)u(x,t),则v(x,t)满足方程2vt2=a22vx2-3-1.2习题选讲因此,根据达朗贝尔公式,v(x,t)的通解可写为v(x,t)=F(xat)+G(x+at),从而u(x,t)=F(xat)+G(x+at)hx
(2)根据上述变换,v(x,t)所满足的初始条件为t=0:
v=(hx)(x),vt=(hx)(x)因此v(x,t)=12(hx+at)(xat)+(hxat)(x+at)+12aZx+atxat(h)()d,从而u(x,t)=12(hx)?
(hx+at)(xat)+(hxat)(x+at)+1aZx+atxat(h)()d.2.问初始条件(x)与(x)满足怎样的条件时,齐次波动方程初值问题的解仅由右传播波组成?
解:
由达朗贝尔公式的推导可知,G(x)=(x)2+12aZxx0()d+C而左传播波由G(x+at)构成,要使只有右传播波,则应有G(x)常数,即G0(x)=0,所以所满足的约束条件为a0(x)+(x)03.利用传播波法,求解波动方程的古沙(Goursat)问题2ut2=a22ux2,u|xat=0=(x),u|x+at=0=(x),(0)=(0).解:
方程的通解为u(x,t)=F(xat)+G(x+at),利用定解条件,可得u|xat=0=F(0)+G(2x)=(x)u|x+at=0=F(2x)+G(0)=(x)从而F(x)=x2G(0),G(x)=x2F(0),又因为u(0,0)=(0)=(0),于是F(0)+G(0)=(0)=(0).因此u(x,t)=F(xat)+G(x+at)=?
x+at2+?
xat2(0).4.对非齐次波动方程的初值问题(2.5),(2.6),证明:
当f(x,t)不变时,
(1)如果初始条件在x轴的区间x1,x2上发生变化,那么对应的解在区间x1,x2的影响区域以外不发生变化;
(2)在x轴区间x1,x2上所给的初始条件唯一地确定区间x1,x2的决定区域中解的数值.证明:
(1)根据非齐次问题解的表达式可知,影响区域为(x,t)|t0,x1at6x6x2+at如图所示,影响区域外的任何一点(x,t)必满足xx2+at,其依赖区间与x1,x2完全无关.因此,x1,x2上初始条件的变化不会引起影响区域以外的任何一点的解的变化;-4-第一章波动方程
(2)类似的,只要证明齐次方程若在x1,x2中给出(x)=0,(x)=0,则在x1,x2的决定区域中任一点(x,t)上有u(x,t)=0.事实上,对决定区域中任一点(x,t)成立x1+atxx1以及x+at0,t0,u|t=0=(x),ut|t=0=0,(uxkut)|x=0=0,其中k为正常数.解:
方程的通解为u(x,t)=F(xat)+G(x+at)当xat0时,有u(x,t)=12(xat)+(x+at),xat,当xat0时,根据边界条件,有uxkut|x=0=(1+ak)F0(at)+(1ak)G0(at)=0从而可得,u(x,t)=12(x+at)+1ak2(1+ak)(atx)+ak1+ak(0),0xat6.求解初边值问题uttuxx=0,0t1,u|t=0=0(x),x0,ut|t=0=1(x),x0,ut|t=kx=(x),其中0(0)=(0).解:
方程的通解为u(x,t)=F(xt)+G(x+t)-5-1.2习题选讲当xt0时,根据达朗贝尔公式的推导有,u(x,t)=12(0(x+t)+0(xt)+12Zx+txt1()d当t=x时,它应与达朗贝尔公式一致,即F(0)+G(2x)=12(0(2x)+0(0)+12Z2x01()d当t=kx时,根据边界条件,有F(1k)x)+G(1+k)x)=(x)根据上是可以求得F和G,从而可得u(x,t)=?
xt1k+12?
0(x+t)0?
1+k1k(xt)+12Zx+t1+k1k(xt)1()d,xtkx7.求解下述边值问题uttuxx=0,0tf(x),u|t=x=(x),u|t=f(x)=(x),其中(0)=(0),t=f(x)为由原点出发的,介于特征线x=t与x=t之间的光滑曲线,且对一切x,f0(x)6=1解:
方程的通解为u(x,t)=F(xt)+G(x+t)根据条件有u|t=x=F(0)+G(2x)=(x)解得G(x)=x2F(0)从而u|t=f(x)=F(xf(x)+G(x+f(x)=F(xf(x)+?
x+f(x)2F(0)=(x)令y=xf(x),由于f0(x)6=1,可以解得x=h(y).由此可得,F(y)=(h(y)h(y)y2+F(0)代入通解可得u(x,t)=?
x+t2?
h(xt)xt2+(h(xt)8.求解波动方程的初值问题2ut22ux2=tsinxu|t=0=0,utflflflflt=0=sinx解:
根据非齐次初值问题解的表达式,有a=1,=0,=sinx,f=tsinx,代入可得u(x,t)=12Zx+txtsind+12Zt0Zx+(t)x(t)sindd=tsinx-6-第一章波动方程9.求解波动方程的初值问题utt=a2uxx+tx(1+x2)2,u|t=0=0,ut|t=0=11+x2.解:
根据非齐次初值问题解的表达式,有=0,=11+x2,f=11+x2,代入可得u(x,t)=12aZx+atxat11+2d+12aZt0Zx+a(t)xa(t)(1+2)2dd=18a34atarctan(x)+ln1+(x+t)21+(xt)22x+at2a2arctan(x+t)+2xat2a2arctan(xt)3.初边值问题的分离变量法1.用分离变量法求下列问题的解:
(1)2ut2=a22ux2u|t=0=sin3xl,utflflflflt=0=x(lx),(0xl),u(0,t)=u(l,t)=0
(2)2ut2a22ux2=0,u(0,t)=0,ux(0,t)=0,u(x,0)=hlx,ut(x,0)=0.解:
(1)根据解的表达式有,Ak=2lZl0sin3lsinkld=1,k=3,0,k6=3.Bk=2kaZl0(l)sinkld=8l3(2n1)44a,k=2n10,k=2nn=1,2,3,从而,u(x,t)=cos3atlsin3xl+8l24aXn=11(2n1)4sin(2n1)atlsin(2n1)xl
(2)此时有通解u(x,t)=Xn=1?
Ancos(n+1/2)atl+Bnsin(n+1/2)atlsin(n+1/2)xl-7-1.2习题选讲其中An=2lZl0hlsin(n+1/2)ld=
(1)n2h(n+1/2)22,Bn=0,n=1,2,3,2.设弹簧一端固定,一端在外力作用下作周期振动,此时定解问题归结为2ut2=a22ux2,u(0,t)=0,u(l,t)=Asint,u(x,0)=ut(x,0)=0.求解此问题.解:
令v=uAxlsint,则v满足2vt2a22vx2=Axl2sint,v(0,t)=v(l,t)=0,v(x,0)=0,vt(x,0)=Alx.设v1是对应于上述问题的齐次方程非齐次初始条件的解,v2是相应的非齐次方程齐次初始条件的解,则v=v1+v2.显然,v1=Xn=1Bnsinnatlsinnxl其中Bn=2An2Zl0lsinnxld=
(1)n2Al(n)2a又有v2=Xn=1Zt0Bn()sinna(t)ldsinnxl其中Bn()=2naA2lsinZl0sinnxld=
(1)n+12A2l(n)2a因此Zt0Bn()sinna(t)ld=
(1)n2A2l(2l2n22a2)n22a?
nasintlsinnatl,n6=la
(1)n+1A2ln22a?
1sinttcost,n=la最后得u(x,t)=v1+v2+Axlsint3.求弦振动方程utta2uxx=0,0x0满足以下定解条件的解:
(1)u|x=0=ux|x=l=0,u|t=0=sin32lx,ut|t=0=sin52lx.-8-第一章波动方程
(2)ux|x=0=ux|x=l=0,u|t=0=x,ut|t=0=0.解:
(1)此时有通解u(x,t)=Xn=1?
Ancos(n+1/2)atl+Bnsin(n+1/2)atlsin(n+1/2)xl其中An=2lZl0sin32lsin(n+1/2)ld=(1,n=10,n6=1Bn=l(n+1/2)a2lZl0sin52lsin(n+1/2)ld=2l5a,n=20,n6=2从而有u(x,t)=cos3at2lsin3x2l+2l5asin5at2lsin5x2l
(2)此时有通解u(x,t)=Xn=0?
Ancosnatl+Bnsinnatlcosnxl其中A0=1lZl0()d=1lZl0d=l2An=2lZl0()cosnld=2lZl0cosnld=2l
(1)n1)n22(n6=0)Bn=04.用分离变量法求解初边值问题:
utta2uxx=g,0x0,u|x=0=ux|x=l=0,u|t=0=0,ut|t=0=sinx2l.其中g为常数.解:
令u=vg2a2x(x2l),则v满足方程vtta2vxx=0,0x0,v|x=0=vx|x=l=0,v|t=0=g2a2x(x2l),vt|t=0=sinx2l.方程的通解为v(x,t)=Xn=0?
Ancos2n+12lat+Bnsin2n+12latsin2n+12lx,其中,An=2lRl0()sin2n+12ld,Bn=4(2n+1)aRl0()sin2n+12ld.-9-1.2习题选讲代入可得An=16l2g(2n+1)3a23B0=2la,Bn=0因此,u(x,t)=2lasinat2lsinx2lg2a2x(x2l)16l2ga23Xn=01(2n+1)3cos2n+12latsin2n+12lx.5.用分离变量法求下面问题的解:
2ut2=a22ux2+bsinhx,u|t=0=utflflflflt=0=0,u|x=0=u|x=l=0.解:
根据非齐次方程的齐次化原理,可得u(x,t)=Xn=1Zt0Bn()sinna(t)ldsinnxl其中Bn()=2naZl0bsinhsinnxld=
(1)n+12blsinhla(l2+n22)与无关而Zt0sinna(t)ld=lna?
1cosnatl因此u(x,t)=2bl2sinhla2Xn=1
(1)n+1n(l2+n22)?
1cosnatlsinnxl6.用分离变量法求下面问题的解:
2ut2+2but=a22ux2,(b0),u|x=0=u|x=l=0,u|t=0=hlx,utflflflflt=0=0.解:
设u(x,t)=X(x)T(t),代入方程得XT00+2bXT0=a2X00T有T00a2T+2bT0a2T=X00X=进一步得到特征值问题(X00+X=0,X(0)=X(l)=0,以及T00+2bT0+a2T=0.-10-第一章波动方程可得Xn(x)=Cnsinpnx,n=nl2,以及Tn(t)=ebtAncoshpb2na2t+Bnsinhpb2na2t,nbla根据初始条件,可以求得An=2lZl0hlsinnxld=
(1)n+12hn以及bAn+pb2na2Bn=0,nbla最后可得u(x,t)=Xn=1Tn(t)Xn(x)4.高维波动方程的柯西问题1.利用泊松公式求解波动方程的柯西问题:
(1)(utt=a2(uxx+uyy+uzz),u|t=0=0,ut|t=0=x2+yz;
(2)(utt=a2(uxx+uyy+uzz),u|t=0=x3+y2z,ut|t=0=0.解:
(1)根据泊松公式,利用球面坐标系可得u(x,y,z,t)=14a2tZ20Z0(x+atsincos)2+(y+atsinsin)(z+atcos)a2t2sindd化简得到u(x,y,z,t)=13t3yz+a2t2+3x2
(2)根据泊松公式,利用球面坐标系可得u(x,y,z,t)=t14a2tZ20Z0(x+atsincos)3+(y+atsinsin)2(z+atcos)a2t2sindd化简得到u(x,y,t)=a2t2z+3a2t2x+x3+y2z2.试用降维法导出弦振动方程的达朗贝尔公式.-11-1.2习题选讲解:
将u(x,t),(x),(x)视为三维空间中与y,z无关的函数,由泊松公式可得u(x,t)=t14a2tZ20Z0(x+atsincos)a2t2sindd+14a2tZ20Z0(x+atsincos)a2t2sindd化简即得所证.3.求解平面波动方程的柯西问题utt=a2(uxx+uyy),u|t=0=x2(x+y),ut|t=0=0.解:
根据二维波动方程的泊松公式,可得u(r,t)=12atZat0Z20(x+rcos)2(x+rcos+y+rsin)q(at)2r2rddr化简得到u(x,y,t)=x2(x+y)+a2t2(3x+y)4.求二维波动方程的轴对称解(即二维波动方程的形如u=u(r,t)的解,其中r=px2+y2).解:
二维轴对称波动方程可写为utta21rr?
rur=0,u|t=0=(r),ut|t=0=(r)根据二维波动方程的泊松公式,可得轴对称解为u(r,t)=12atZ20Zat0pr2+22rcosq(at)22dd+12aZ20Zat0pr2+22rcosq(at)22dd5.求解下列柯西问题:
utt=a2(uxx+uyy)+c2u,u|t=0=(x,y),ut|t=0=(x,y).解:
令v(x,y,z,t)=u(x,y,t)ecz/a,则v满足方程vtt=a2v,v|t=0=(x,y)ecz/a,vt|t=0=(x,y)ecz/a.利用泊松公式即可求解.-12-第一章波动方程6.试用齐次化原理导出平面非齐次波动方程utt=a2(uxx+uyy)+f(x,y,t)在齐次初始条件(u|t=0=0,ut|t=0=0下的求解公式.解:
根据齐次化原理,求解方程vtt=a2(vxx+vyy),v|t=0,vt|t=f(x,y,)根据泊松公式可得v(x,y,t;)=12aZa(t)0Z20f(x+rcos,y+rsin,)q(a(t)22dd从而u(x,y,t)=Zt0v(x,y,t;)d.7.用降维法来求解上面的问题.略.8.解非齐次方程的柯西问题:
utt=uxx+uyy+uzz+2(yt),u|t=0=0,ut|t=0=x2+yz.解:
根据三维非齐次波动方程的泊松公式,可以得到u=14ZZSMt2+tdSM+14ZZZr6t2(t+r)rdV,从而u=tx2+yz+yt2.5.波的传播与衰减1.试说明:
对一维波动方程所描述的波的传播过程一般具有后效现象.解:
根据达朗贝尔公式,由于初始速度作用,位移项包括Zx+atxat()d,显然具有后效.2.试说明:
对一维波动方程,即使初始资料具有紧支集,当t时其柯西问题的解没有衰减性.解:
根据达朗贝尔公式,由于初始位移向左右传播,不会产生衰减.-13-1.2习题选讲3.设u为初始资料及具有紧支集的二维波动方程的解.试证明:
对任意固定的(x0,y0)R2,成立limt+u(x0,y0,t)=0.证明:
6.能量不等式、波动方程解的唯一性和稳定性1.对受摩擦力作用且具有固定端点的有界弦振动,满足方程utt=a2uxxcut,其中常数c0,证明其能量是减少的,并由此证明方程utt=a2uxxcut+f,的初边值问题解的唯一性以及关于初始条件及自由项的稳定性.证明:
(1)弦的能量是E(t)=12Zl0u2t+Tu2xdx=2Zl0u2t+a2u2xdx有E0(t)=Zl0ututt+a2uxuxtdx=Zl0ututta2uxxdx+Zl0a2x(uxut)dx因此E0(t)=cZl0u2tdx60所以能量是减少的;
(2)唯一性.假设u1,u2为方程的解,则v=u1u2满足方程vtt=a2vxxcvt,从而能量减少.由于E(0)=0,故E(t)0,因此v0,即u1=u2.稳定性.记E0(t)=Zl0u2dx,则有能量不等式E0(t)etE0(0)+AE(0)et1,进一步可证关于初始条件及自由项的稳定性.2.证明函数f(x,t)在G:
06x6l,06t6T作微小改变时,方程2ut2=x?
k(x)uxq(x)u+f(x,t),(其中k(x)0,q(x)0和f(x,t)都是一些充分光滑的函数)具固定端点边界条件的初边值问题的解在G内的改变也是微小的.证明:
记E(t)=Zl0ut2+kux2+qu2dx,则有能量不等式E(t)etE(0)+etZt0F()ed,-14-第一章波动方程其中F()=Zl0f(x,t)dx.记E0(t)=Zl0u2dx,有能量不等式E0(t)etE0(0)+Zt0F()ed.进一步可证明解的稳定性.3.证明波动方程utt=a2(uxx+uyy)+f(x,y,t),的自由项f在L2(K)意义下作微小改变时,对应的柯西问题的解在L2(K)意义之下改变也是微小的.证明:
类似于课本,在特征锥内积分,考虑到边界条件,有如下能量不等式E1(t)6etE1(0)+etZt0F()ed,E0(t)6etE0(0)+Zt0E1()ed,进一步可证关于自由项的稳定性.4.固定端点有界弦的自由振动可以分解成各种不同固有频率的驻波(谐波)的叠加,试计算各个驻波的动能和位能,并证明弦振动的总能量等于各个驻波能量的叠加.这个物理性质对应的数学事实是什么?
解:
各阶驻波为un(x,t)=?
Ancosnatl+Bnsinnatlsinnxl,其中An=2lZl0()sinnld,Bn=2naZl0()sinnld.因此,动能Un=2Zl0un,t2dx=n22a24l?
Ansin?
natlBncos?
natl2位能Vn=T2Zl0un,x2dx=Tn224l?
Ancos?
natl+Bnsin?
natl2.能量En=Un+Vn=n22a24lAn2+Bn2.在振动过程中,总能量保持不变,故E=E(0)=2Zl0ut2+a2ux2flflt=0dx=2Zl02+a202dx-15-1.2习题选讲根据Fourier级数的性质,有E=Xn=1En,即得所证.5.考虑波动方程的第三类初边值问题utta2(uxx+uyy)=0,t0,(x,y),u|t=0=(x,y),ut|t=0=(x,y),?
un+uflflflfl=0,其中0是常数,是的边界,n为上的单位外法线向量.对于上述定解问题的解,定义能量积分E(t)=ZZut2+a2ux2+uy2dxdy+a2Zu2ds,试证明E(t)常数,并由此证明上述定解问题解的唯一性.证明:
类似于课本,在特征锥内积分,考虑到边界条件,可得dE(t)dt
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数学 物理 方程 习题 讲义