《数学物理方法》第十章作业参考解答 10.1一无限长的弦,弦上x=x0点受到突然的冲击(冲量为D。求此后弦的振 动情况 解:定解问题: 6(x-x0) 根据d’ Alembent公式 l(:)=1 2ajx-at (5-x05=[H(x-x0+a)-H(x-x0-a1) 2ap 10.6长为2l的均匀杆,两端受压而使长度缩为l(2-2E),放手后任其自由振动, 求杆的振动情况。 解一]:设杆的中点为坐标原点,显然杆的振动是关于坐标原点x=0对称 的,因此只需求解0≤x≤l段。定解问题是 0≤x≤1) 0 设u(x,1)=X(x)7(t) (1) 可得,X"+AX=0 把①代入边界条件可得 X(0)=0X(1)=0 (4) 方程(2)和边界条件(4)构成本征值问题,当λ≤0时,方程无非零 解,当λ>0时,方程(2)的解为 X=A, COSv2x+A (5) 由(4)可得A1=0
《数学物理方法》第十章作业参考解答 10.1 一无限长的弦，弦上 0 x = x 点受到突然的冲击（冲量为 I）。求此后弦的振 动情况。 解：定解问题： = = − = = = 0; ( ); 0 0 0 2 x x I u u u a u t t t tt xx δ ρ 根据 d’Alembent 公式， ∫ + − = − = − + − − − x at x at H x x at H x x at a I x d I a u x t [ ( ) ( )] 2 ( ) 2 1 ( , ) 0 0 0 ρ δ ξ ξ ρ 10.6 长为2l 的均匀杆，两端受压而使长度缩为l(2 − 2ε ) ，放手后任其自由振动， 求杆的振动情况。 [解一]：设杆的中点为坐标原点，显然杆的振动是关于坐标原点 x = 0对称 的，因此只需求解0 ≤ x ≤ l 段。定解问题是 ( ) = = = − = = ≤ ≤ = = = = 0, 0 ; 0; 0 0 0 0 2 x x x l t t t tt xx u u u x u u a u x l ε 设u(x,t) = X (x)T(t) （1） = = −λ a T T X X 2 ' ' ' ' 可得， X ' '+λX = 0 ， （2） '' 0 2 T +λa T = ， （3） 把①代入边界条件可得 X (0) = 0 X '(l) = 0 （4） 方程（2）和边界条件（4）构成本征值问题，当λ ≤ 0时，方程无非零 解，当λ > 0 时，方程（2）的解为 X A cos λ x A sin λ x = 1 + 2 （5） 由（4）可得 0 A1 =
X()=√2A2cos√l-=0 所以本征值为2=2「(2n+1)z (n=0,1,2,3…) 相应的本征函数为X=Xn(x)=sin (2n+1)丌 (2n+1) + ]2T=0 T(o=T,(0=Ccos (2n+1)m +1) t+D sin 27 27 n+1)m l(x,1)=∑(C D 27 27 C (2n+1)丌 代入初始条件可得, Dn (2n+1)m.smn(2n+1)z 27 2 系数C,D分别为 (2n+1)mx 8al (2n+1)mx B=0 (x,)=8a/ ∑(-1)-1 (2n+1) (2n+1)mx (2n+1)2 按定解问题,此解仅适用于0≤x≤l。根据对称性上式在-1≤x≤l 也成立。 [解二] 将坐标原点取在杆的左端,则解定解问题为 =0 =0. 其解为 (x,)= ∑ COs cos (2n+1
X '(l) λ A cos λl = 2 =0 ⇒ cos λl =0 所以 本征值为 2 2 (2 1) + = = l n n π λ λ (n=0,1,2,3…) 相应的本征函数为 x l n X X x n 2 (2 1) ( ) sin + π = = (6) ] 0 2 (2 1) '' [ 2 = + + T l n a T π t l n a t D l n a T t Tn t Cn n 2 (2 1) sin 2 (2 1) ( ) ( ) cos π + π + + = = x l n t l n a t D l n a u x t C n n n 2 (2 1) ) sin 2 (2 1) sin 2 (2 1) ( , ) ( cos 0 π π + π ⋅ + + + = ∑ ∞ = 代入初始条件可得， + ⋅ + = + − = ∑ ∑ ∞ = ∞ = x l n l n a D x l n x C n n n n 2 (2 1) sin 2 (2 1) 0 2 (2 1) sin 0 0 π π π ε 系数Cn Dn , 分别为 2 2 1 0 0 2 (2 1) 8 ( 1) 2 (2 1) sin 2 (2 1) sin 1 π π ε ε π + = − + − + = + ∫ ∫ n l dx l n x x dx l n x C n l l n Bn = 0 l n x l n at n l u x t n n 2 (2 1) sin 2 (2 1) cos (2 1) 1 ( 1) 8 ( , ) 2 1 0 2 π π π ε + + + = − + ∞ = ∑ 按定解问题，此解仅适用于0 ≤ x ≤ l 。根据对称性上式在 − l ≤ x ≤ l 也成立。 [解二] 将坐标原点取在杆的左端，则解定解问题为 ( ) = = = − = = ≤ ≤ = = = = 0, 0 ( ); 0; 0 2 0 2 0 0 2 x x x x l t t t tt xx u u u l x u u a u x l ε 其解为 ∑ ∞ = + + + = 0 2 2 2 (2 1) cos 2 (2 1) cos (2 1) 8 1 ( , ) n l n x l n at n l u x t π π π ε
[解三] 将坐标原点取在杆的中点,直接写出杆的解定解问题为 (-l≤x≤1) 0 0 l 设u(x,)=X(x)7(t),可得, T+ha t=o X"+aX=0, X(-1) X()=0 方程(2)和边界条件(3)构成本征值问题,当λ<0时,方程无非零 解,当λ=0时,得X=X0=A(常数),当λ>0时,方程(2)的解 为 X=A cos 由边界条件(3)可得 A1sin√+√A2c (5) √Asin√+√A2 COSTal=0 为使A1,A2不同时为零,必须 sIn Acos√ 0,即,2sin√lcos√l=0,由此得到, sin√ 2 cos√Al 1=1 (n=123.…) 为求出相应的本征函数,将(6)代入(5)中的任意一式,得 由此可见,当n=2m(m=1,2,3,…)时,A2=0,相应本征函数是 X=Xn(x)=X? 当n=2m+1(m=0,12,3…)时,A1=0,相应本征函数是
[解三] 将坐标原点取在杆的中点，直接写出杆的解定解问题为 ( ) = = = − = = − ≤ ≤ =− = = = 0, 0 ; 0; 0 0 2 x x l x x l t t t tt xx u u u x u u a u l x l ε 设u(x,t) = X (x)T(t)， 可得， '' 0 2 T +λa T = ， （1） X ' '+λX = 0 ， （2） X ′(−l) = 0 X '(l) = 0 （3） 方程（2）和边界条件（3）构成本征值问题，当λ < 0时，方程无非零 解，当λ = 0 时，得 X = X0 = A（常数），当λ > 0 时，方程（2）的解 为 X A cos λ x A sin λ x = 1 + 2 （4） 由边界条件（3）可得 − + = + = sin cos 0 sin cos 0 1 2 1 2 A l A l A l A l λ λ λ λ λ λ λ λ （5） 为使 1 2 A , A 不同时为零，必须 0 sin cos sin cos = − l l l l λ λ λ λ ，即，2sin λl cos λl = 0，由此得到， 2 2 = = l n n π λ λ （n = 1,2,3,L） （6） 为求出相应的本征函数，将（6）代入（5）中的任意一式，得 0 2 cos 2 1 sin + 2 = π nπ A n A 由此可见，当n = 2m （m = 1,2,3,L）时， 0 A2 = ，相应本征函数是 x l m X Xn x X m π ( ) cos = = 2 = 当n = 2m +1（m = 0,1,2,3,L）时， 0 A1 = ，相应本征函数是
X=Xn(x)=x2m+I=sin (2m+1) 对应每一个本征值,解方程(1),可得 Do (=4 (1)= t+d sin m7 Tml(t=Emcs- (2m+1)m (2m+1)m t2=λ (2m+1)丌 Fns u(x, t)=Co!+Do+2(Cm cos"/t+Dm sin-1)cos 7 +∑(En0"+ 1) F 2/1)sin(2m+1)z (2m+1)m 代入初始条件可得, a=D0+ (2m+1) cos -x+ 27 0=c0+>Dc3 IIx+>F(2n 1)m(2n+1)x sIn 21 所以系数分别为 arax aX COs dx=o E sin (m+1) I'-Exsin (2m+)7r (2m+1)2z C0=0,Dn=0,Fn=0,因此 (x,D) (2m+1)mt(2m+1)z cos (2m+1) 2 10.9求二维无限深势阱的 Schrodinger方程
( ) x l m X Xn x X m 2 2 1 ( ) 2 1 sin + π = = + = 对应每一个本征值，解方程（1），可得 0 0 0 T (t) = C t + D ( 0) λ = λ0 = t l m a t D l m a T m t Cm m π π 2 ( ) = cos + sin = = 2 2 l m m π λ λ t l m a t F l m a T m t Em m 2 (2 1) sin 2 (2 1) ( ) cos 2 1 π + π + + + = + = + = 2 2 1 2 (2 1) l m m π λ λ x l m t l m a t F l m a E x l m t l m a t D l m a u x t C t D C m m m m m m 2 (2 1) )sin 2 (2 1) sin 2 (2 1) ( cos ( , ) ( cos sin ) cos 0 1 0 0 π π π π π π + + + + + = + + + ∑ ∑ ∞ = ∞ = 代入初始条件可得， + + = + + + − = + + ∑ ∑ ∑ ∑ ∞ = ∞ = ∞ = ∞ = 1 1 0 1 0 0 2 (2 1) sin 2 (2 1) 0 cos 2 (2 1) cos sin m m m m m m m m x l n l n a x F l m l m a C D x l m x E l m x D C π π π π π π ε 所以系数分别为 d 0 d 1 0 = − = ∫ ∫ − − l l l l x x x D ε cos d 0 cos d 1 2 − = = ∫ ∫ − − l l l l m x l m x x x l m x C π ε π 2 2 1 2 (2 1) 8 d ( 1) 2 (2 1) sin d 2 (2 1) sin 1 π π ε ε π + = − + − + = + − − ∫ ∫ m l x l m x x x l m x E m l l l l m C0 = 0 ， Dm = 0 ， Fm = 0，因此， ∑ ∞ = + + + + = − 0 2 1 2 2 (2 1) sin 2 (2 1) cos (2 1) 1 ( 1) 8 ( , ) m m l m x l m at m l u x t π π π ε 10.9 求二维无限深势阱的 Schrodinger 方程
h2 a 2m ax Ev(x,y)(0≤xsa0≤y≤a) =0 0 0 中能级E的最小值和次小值以及相应的归一化波函数 令v(x,y)=(x)Y(y)代入原薛定谔方程 得 (x) r"() 2mE X(x) 可设 X(x) "(y) Y(y) 其中,元+ 2mE h 由边界条件,得 X(0)=0,X(a)=0 F(0)=0,Y(a)=0 (5) 方程(1)和边界条件(4)构成本征值问题,解得 本征值n=(-)2,(n=1,2,3,…) 相应本征函数X nza a 相似地,由方程(2)和边界条件(5)构成的本征值问题,得 本征值 相应本征函数yn(y)=sinm (m=1,2,3,…) 由(3)式得, (n+4n) +m2)丌 相应的本征函数v= Csin -sin
( ) () () ( ) = = = = = ≤ ≤ ≤ ≤ ∂ ∂ + ∂ ∂ − = = = = 0, 0 0, 0 , , , 0 ,0 2 0 0 2 2 2 2 2 y y a x x a x y E x y x a y a y x y m x ψ ψ ψ ψ ψ ψ ψ h 中能级 E 的最小值和次小值以及相应的归一化波函数。 解： 令 ψ (x, y) = X (x)Y( y) 代入原薛定谔方程 得 2 2 ( ) ''( ) ( ) ''( ) h mE Y y Y y X x X x + = − 可设 = −λ ( ) ' '( ) X x X x （1） = −µ ( ) ''( ) Y y Y y （2） 其中， 2 2 h mE λ + µ = ， （3） 由边界条件，得 X (0) = 0, X (a) = 0 （4） Y(0) = 0,Y(a) = 0 （5） 方程（1）和边界条件（4）构成本征值问题，解得 本征值 ( ) , 2 a n n π λ = （n =1,2,3,…） 相应本征函数 a n x X x n π ( ) = sin 相似地，由方程（2）和边界条件（5）构成的本征值问题，得 本征值 2 ( ) a m m π µ = 相应本征函数 a m y Y y m π ( ) = sin （m =1,2,3,…） 由（3）式得， 2 2 2 2 2 2 2 ( ) ( ) 2 ma n m m Enm n m h π h λ µ + = + = 相应的本征函数 a m y a n x nm C π π ψ = sin sin