2.6 固有值和固有函数
数学物理方程练习题第七版(学生用)

= u(0, t) 0= , ux (2,t) 1,
u(x= ,0)
cos π x + x3 − 3x2 − x.
2
3.求定解问题的解:
u
x= x + u yy
sinπ x,
0 < x < 1, 0 < y < 1,
= u(0, y) 1,= u(1, y) 2,
u(x,0) =1+ x,
7
u
rr
+
1 u
r
r
+
1 r2
uθθ
= 0,
u= (1,θ ) A cosθ (−π < θ ≤ π ).
4. 设 A, B 为常数,用试探法求如下定解问题的解:
u rr
1 +rur
+
1 r2
u
θθ
=
0,
r < a,
u r= =a A cosθ + B sinθ (−π < θ ≤ π ).
练习十五
练习六
1.求解如下定解问题:
ut = uxx + cosπ x, (0 < x < 1, t > 0), u= x (0,t) u= x (1,t) 0, u(x,0) = 0.
3
2.求解如下定解问题:
= u tt
a2u
xx
+
t
sin
π l
x
,
u= (0,t) u= (l,t) 0, t ≥ 0,
X= ′(0)
X= (l)
0.
3. 求如下定解问题的解:
= ut uxx , 0 < x < 2, t > 0, ux= (0, t) u= (2, t) 0,
25-26具有非齐次边界条件的问题以及固有值和固有函数资料

6
utt a 2uxx f (x, t) (0 x l, t 0), (79) 若边界条件不全是第一类,也可采用类似方法 把非齐次边界条件化成齐次的。 我们就下列几种
非齐次边界条件的情况,分别给出相应辅助函数
w(x,t) 的表达式:
(1) u(0,t) u1(t),
于是可得
w(t, x)
x l
[u2
(t
)
u1
(t
)]
u1
(t
).
因此,令
u( x, t )
v( x, t )
x l
[u2 (t)
u1(t)]
u1 (t ).
则问题(79)-(81)可化成v(x,t) 的定解问题
(79) (80) (81) (82)
(85)
4
utt a 2uxx f (x, t) (0 x l, t 0),
4
vt (x,0) sin l x.
为了将 v(x,t) 的边界条件也化成齐次,则 w(x)满足
w(0) 3,
w(l) 6.
13
utt a 2uxx u(0,t) 3,
sin 2 x cos 2
l
l
u(l,t) 6,
x
(0
x l,
t
0),
(91)
u ( x,0)
l
x,
为了将此方程化成齐次的,自然选取w(x) 满足
a2w sin 2 x cos 2 x 0.
l
l
12
例2 求解下列问题:
utt a 2uxx u(0,t) 3,
固有值问题和特殊函数Bessel方程

(1) y j (2) yi 得 d dyi d dy j [x ]y j [x ] yi (ki 2 k j 2 ) xyi y j 0,即 dx dx dx dx d (ki 2 k j 2 ) xyi y j [ x( y j yi yi y j )]. dx 再从0到1积分得, (ki k j ) xyi y j dx d [ x( y j yi yi y j )]
因为a1 0,故由上式可得 a2 k 1 0,k 0,1, 2, a2 k a2 k 2 a 2 2 k 2 , k 1, 2, 2k (2 2k ) 2 k ( k ) (5)
反复应用(5),有 (1) 2 a2 k 22 a2 k 22 2 k (k 1)( k )( k 1) (1)3 a2 k 23 23 2 k (k 1)(k 2)( k )( k 1)( k 2) (1) k a0 2k 2 k !( k )( k 1) ( 1) (1) k 2 ( 1)a0 2 k . 2 k !( k 1)
四、Bessel函数的正交性及模
1.含参数的Bessel方程 x 2 y xy ( 2 x 2 n 2 ) y 0 (*) y (0) 且满足端点条件 的固有值问题。 y (a) 0 在(*)中令 x,方程变为
2 2 y y ( n ) y 0. 2
(ii)s1 s2 2 2n(n 1,2),即s1 n, s2 n.
(1)k x2k n (1)k x2k n 此时J n ( x) 2k n , J n ( x) 2k n . k !(n k 1) k !(n k 1) k 0 2 k 0 2 1 而当 n k 1为零和负整数时(n k 1) ,即 0, (n k 1)
数理方程常规例题I

数学物理方程常规例题I(1-20题)一、数学模型例题例1. 密度为ρ均匀柔软的细弦线x =0端固定,垂直悬挂,在重力作用下,横向拉它一下,使之作微小的横振动。
试导出振动方程。
解:考虑垂直悬挂的细弦线上一段微元ds ,该微元在坐标轴上投影为区间[x ,x+d x ],在微元的上端点处有张力:)(1x L g T -=ρ,在下端点处有张力:)(2dx x L g T --=ρ考虑张力在位移方向的分解,应用牛顿第三定律,有tt u m T T =-1122sin sin αα 由于细弦作微小振动,所以有近似)(tan sin 22dx x u x +=≈αα )(tan sin 11x u x =≈αα代入牛顿第三定律的表达式,有tt x x u ds t x u x L g t dx x u dx x L g ρρρ≈--+--),()(),()(上式两端同除以ds ρ,得tt x x u dsx u x L dx x u dx x L g≈--++-)()()())((由于dx ds ≈,而x x x x x u x L dxx u x L dx x u dx x L )]()[()()()())((-≈--++-所以,细弦振动的方程为tt x x u u x L g =-])[(例2. 长为L 密度为ρ底半径为R 的均匀圆锥杆(轴线水平)作纵振动,锥顶点固定在x =0处。
导出此杆的振动方程。
(需要包括假设在内的具体推导) 解:设均匀圆锥杆作纵振动时位移函数为u (x ,t )则在点x 处,弹力与相对伸长量成正比,即),(),(t x Yu t x P x = 其中,Y 为杨氏模量。
在截面上张力为T (x , t ) = S (x ) P (x , t )这里,S (x )为x 处圆锥截面积。
考虑圆锥杆上对应于区间[x ,x+dx ]处的微元(如右图所示)。
应用牛顿第二定律,得),()]()()[(31),(),(t x u x xS dx x S dx x t x T t dx x T tt -++=-+ρ 由于圆锥截面积2)()(x LR x S π= 微元(圆台)体积)33()(31)]()()[(313222dx xdx dx x LRx xS dx x S dx x ++=-++ρπρ 所以),()33()(31)],(),()[()(3222222t x u dx xdx dx x L Rt x u x t dx x u dx x L R Y tt x x ++=-++ρππ两端除dx ,并取极限,得),()],([22t x u x t x u x Y tt x x ρ=记ρ/2Y a =,则有方程)2(2x xx tt u xu a u += 二、二阶偏微分方程化简与求通解只考虑未知函数是两个自变量情形,即),(y x u 。
数学物理方程 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 Iu 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 nuk )(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 方程:2ερ-=∇u 。
习题2.42.求下列方程的通解:(2):;032=-+yy xy xx u u u (5):;031616=++yy xy xx u u u解:(2):特征方程:03)(2)(2=--dx dy dx dy解得:1-=dx dy 和3=dxdy。
2.6 固有值和固有函数

x 2 y ′′ + xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0 1 固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。
{ (1)首先求出固有函数系 (1)首先求出固有函数系 y n (x)}的具体表达式 t t = ln x 作变换 x = e 则有
y1 ( x), y 2 ( x), ⋯, y n ( x), ⋯.
(2) 如果把对应于固有值 λn 的固有函数记为 y n (x), 那么所有y n (x)组成一个带权函数 ρ (x) 的正交函数 系,即 b (m ≠ n). (96) ∫a ρ ( x)y m ( x) y n ( x)dx = 0 例如
解
1 y x = yt ⋅ , x
1 1 1 1 1 y xx = ( y tt ⋅ ) ⋅ + y t ⋅ (− 2 ) = 2 y tt − 2 y t , x x x x x
代入原方程有
y tt − y t + y t + λy = 0
y tt + λy = 0
10
15. 试证问题
齐次欧拉方程
x 2 y ′′ + xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0 1 固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。
解
y tt + λy = 0
y (0) = y (1) = 0.
λ n = ( nπ ) 2 ,
f ( x) = ∑ c n y n ( x),
n =1 ∞
(97)
数学物理方程与特殊函数课后答案

29.0(,)11cos ,sin (,)(cos ,sin ),cos sin ;sin cos .sin cos ;s xx yy rr r r x y x y x r y laplace u u r u u u r rx r y r u x y u r r u u u u r u r u u u u ru θθθθθθθθθθθθθθθ+=++==⎧⎨=⎩∴==+⎧⎪⎨=−+⎪⎩=−⇒=∵ 证明方程在极坐标下为 证明: sin cos ;cos cos in .sin .sin ()cos ()sin sin cos cos r xx x r r u u r y r r u u u x x r r x u u r r r r θθθθθθθθθθθθθθθθθθ⎧∂∂∂⎛⎞⎧=−⎜⎟⎪⎪∂∂∂⎝⎠⎪⎪⇒⎨⎨∂∂∂⎛⎞⎪⎪+=+⎜⎟⎪⎪⎩∂∂∂⎝⎠⎩∂∂∂∂∂⎛⎞==−⎜⎟∂∂∂∂∂⎝⎠∂∂∂∂⎛⎞⎛=−−⎜⎟⎜∂∂∂∂⎝⎠⎝ 从而2222222222222sin cos sin cos sin cos sin cos sin cos sin .cos ()sin ()sin yy u u u u r r r r r r u u ur r r r u u u y y r r y θθθθθθθθθθθθθθθθθθθ⎞⎟⎠∂∂∂∂=+−+∂∂∂∂∂∂∂∂−++∂∂∂∂∂∂∂∂∂⎛⎞==+⎜⎟∂∂∂∂∂⎝⎠= 2222222222222cos cos sin sin cos sin cos cos sin sin cos sin cos cos .1u u r r r r u u u u r r r r r r u u ur r r r u u u u θθθθθθθθθθθθθθθθθθθθθθ∂∂∂∂⎛⎞⎛⎞++⎜⎟⎜⎟∂∂∂∂⎝⎠⎝⎠∂∂∂∂=−++∂∂∂∂∂∂∂∂+−+∂∂∂∂+=+ 所以 10.u +=21.(01,0),(0,)(1,)0,1,0.(2)2(,0)11,1,2(,0)(1);tt xx tu a u x t u t u t x x u x x x u x x x ⎧=<<>⎪==⎪⎪⎧⎪<≤⎪⎨⎪=⎨⎪⎪⎪−<<⎪⎩⎪⎪=−⎩求下列问题的解22(,)()().()()0,()()0.(0)(1)0.()()0,(0)(1)0.(),()si n n n u x t X x T t T t a T t X x X x X X X x X x X X n X x B λλλλπ=′′+=′′+===′′+=⎧⎨==⎩==解:应用分离变量法,令 代入方程分离变量,得 由边界条件分离变量,得 求解固有值问题得, 111212202n (1,2,).()cos sin (1,2,).(,)(cos sin )sin .42sin (1)sin sin .2n n n n n n n n x n T t C an t D an t n u x t a an t b an t n x n a x n xdx x n xdx n ππππππππππ∞===+==+⎡⎤=+−=⎢⎥⎣⎦∑∫∫ 代入另一常微分方程,得则其中 ()()14402244124(1)sin 11.44(,)(sin cos 11sin )sin .2nn nn b x x n xdx an n a n u x t an t an t n x n n a πππππππππ∞=⎡⎤=−=−−⎣⎦⎡⎤=+−−⎣⎦∫∑ 因此,所求定解问题的解为2(0,0),(0,)(,)0,(3)35(,0)3sin6sin ,22(,0)0.tt xx x t u a u x l t u t u l t x xu x l l u x ππ⎧=<<>⎪==⎪⎪⎨=+⎪⎪=⎪⎩ ()22(,)()().()()0,()()0.(0)()0.()()0,(0)()0.21(),(2n n u x t X x T t T t a T t X x X x X X l X x X x X X l n X l λλλπλ=′′+=′′+=′==′′+=⎧⎨′==⎩+=解:应用分离变量法,令 代入方程分离变量,得 由边界条件分离变量,得 求解固有值问题得, ()()()()()()121)sin (0,1,2,).22121()cossin (0,1,2,).22212121(,)(cossin )sin .222235(3sin6sin 22n n n n n n n n n x B x n la n a n T t C t D t n l la n a n n u x t a tb t x l l l x x a l l ππππππππ∞=+==++=+=+++=+=+∑ 代入另一常微分方程,得则 其中 ()03,1;21)sin 6,2;20,12.0.3355(,)3cos sin 6cos sin .2222l n n n xdx n l l n b a a u x t t x t x l l l lπππππ=⎧+⎪==⎨⎪≠⎩==+∫、 因此,所求定解问题的解为3.4(0,0),(2)(0,)0,(,)0,(,0)().t xx x x u u x l t u t u l t u x x l x =<<>⎧⎪==⎨⎪=−⎩求下列定解问题的解:2(,)()().()4()0,()()0.(0)()0.()()0,(0)()0.(),()n n u x t X x T t T t T t X x X x X X l X x X x X X l n X x A lλλλπλ=′+=′′+=′′==′′+=⎧⎨′′==⎩==解:应用分离变量法,令 代入方程分离变量,得 由边界条件分离变量,得 求解固有值问题得, 222()2()012000cos (0,1,2,).()(0,1,2,).1(,)cos .222().62()cos n n t ln n n t ln n l l n n x n l T t D e n n u x t a a e x l l a x l x dx l n a x l x xd l l πππππ−∞−=====+=−==−∑∫∫ 代入另一常微分方程,得则 其中 2222222()2212[1(1)].2[1(1)](,)cos .6n n n t ln l x n l l n u x t e x n lππππ∞−=−−+−=−−+−=+∑ 因此,所求定解问题的解为2110(01),,0,(1,)0,.,.rr r u u u r r r A u A θθθαθαθπα⎧++=<<⎪⎪⎨⎧≤≤⎪⎪=⎨⎪<≤⎪⎩⎩其中为已知常数22(,)()().()()()0,()()0.()()0,()(2).(),()cos sin n n n n u r R r r R r rR r R r n X x A n B n θθλθλθθλθθθπλθθ=Φ′′′+−=′′Φ+Φ=′′Φ+Φ=⎧⎨Φ=Φ+⎩==+解:应用分离变量法,令 代入方程分离变量,得求解固有值问题得,()2010(0,1,2,).()()()0,(0).()(0,1,2,).1(,)cos sin .212n n n n n n n n n r R r rR r R r R R r C r n u r a a n b n r Aa Ad a ααλθθθαθππ∞=−=′′′⎧+−=⎨<+∞⎩===++==∑∫ 代入另一常微分方程的定解问题得, 则 其中 112cos sin ,1sin 0.2(,)sin cos .n nn AA n d n n b A n d A A u x t r n n n ααααθθαππθθπααθππ−−∞======+∫∫∑ 因此,所求定解问题的解为0(0,0),(0,)0,(,)0(0),(,0)(1),lim (,)0(0),.xx yy y u u x l y u y u l y y x u x A u x y x l l A →∞⎧+=<<<<∞⎪⎪==≤<∞⎨⎪⎪=−=<<⎩其中为已知常数 2(,)()().()()0,()()0.(0)()0.()()0,(0)()0.(),()sin n n n u x y X x Y y X x X x Y y Y y X X l X x X x X X l n X x B lλλλπλ=′′+=′′−===′′+=⎧⎨==⎩==解:应用分离变量法,令 代入方程分离变量,得 由边界条件分离变量,得 求解固有值问题得, 10(1,2,).()(1,2,).(,)sin.22()sin .lim (,)0n n y y lln n n n n y y l ln n n l n n y n x n l Y y C e D e n n u x y a e b e x l x n A a b A l xdx l l l n u x y a ππππππππ−∞−=→∞==+=⎛⎞=+⎜⎟⎝⎠+=−==⇒∑∫ 代入另一常微分方程,得则 其中 10.2(,)sin .n n y l n A n u x t e x n l πππ∞−===∑因此,所求定解问题的解为()22228.-10.cos ,sin ,111(0),0.{cos sin }.,()xx yy x y a rr r r an a u u u x r y r u u u r a r r u A n B n u r a r θθθθθθθ+==+====⎧++=−<<⎪⎨⎪=⎩+= 在以原点为心,为半径的圆内,试求泊松方程 的解,使它满足边界条件解:令作极坐标变换,得由固有函数法,相应的固有函数系为 因此,设方程的解为[]()()()()()()()0002222cos ()sin .11,110,0210,323()0()n n n n n n n n n nn n nn n n n b r n a a r n a a a n r r nb b b r r a r A r B r n b r C r D θθ∞=−+⎧′′′+=−⎪⎪⎪′′′+−=≠⎨⎪⎪′′′+−=⎪⎩=+≠=+∑ 代入方程,得方程,的通解:, ()()2000(0),()0;(0),()0.()00()0.11()ln ,4(0),()n n n n n n n n r a a a b b a a r n b r a r A r B r a a a −<+∞=<+∞==≠==+−<+∞=. 由有界性条件及边界条件,得 , 方程的通解: 由有界性条件及边界条件,()()()()()220222220.1().41,.41,.a r a r u r a r u x y a x y θ=−=−⎡⎤=−+ 得 则定解问题的解为 化成直角坐标,则得21210.sin ,(2)(0,)0,(,)0(0),(,0)0,(,0)0(0);{sin }.(,)()sin .tt xx tn n n u a u t x l u t u l t t u x u x x l n x ln u x t u t x l n a u u l ππππ∞=⎧=+⎪⎪==≥⎨⎪==≤≤⎪⎩=⎛⎞′′+⎜⎟⎝⎠∑求下列问题的解:解:由固有函数法,相应的固有函数系为 设方程的解为 代入原方程,得()2111020(1),.(0)(0)0(1,2,),1()0;1()sin sin .n n n n t n a u u t l u u n n u t l an u t t d al l l a t t a a l ππτττππππ=≠⎛⎞′′+=⎜⎟⎝⎠′===≠===−⎛⎞⎛⎞=−⎜⎟⎜⎟⎝⎠⎝⎠∫"" 由初始条件,得当时, 当时, 2(,)sin sin l l a u x t t t x a a l l ππππ⎛⎞⎛⎞=−⎜⎟⎜⎟⎝⎠⎝⎠ 故所求的解为2110(0,0),(3)(0,)0,(,)0,(,0)0.,{sin}.(,)()sin .sin 22sin [1(t xx n n n n l n u a u A x l t u t u l t u x n x ln u x t u t x l n A A A x l n A A A xdx l l n πππππ∞=∞=⎧=+<<>⎪==⎨⎪=⎩====−∑∑∫ 解:由固有函数法相应的固有函数系为 设方程的解为 并将展为: ,其中 222()023321)].2[1(1)],(0)0.2()[1(1)]2[1(1)][1].(,n n n n n n a t tn l n n a t n ln a A u u l n u Au t e d n Al e n au x πτπππτππ⎛⎞−−⎜⎟⎝⎠⎛⎞−⎜⎟⎝⎠−⎧⎛⎞′+=−−⎪⎜⎟⎨⎝⎠⎪=⎩=−−=−−−∫ 代入原方程可得得: 故所求的解为2233212)[1(1)][1]sin .n a tnl n Al n t e x n alπππ⎛⎞∞−⎜⎟⎝⎠==−−−∑()2211.224sin cos ,(2)(0,)0,(,)(0),(,0),(,0)()(0).(,)(,)().224sin cos ,(0,)(0ttxx t ttxx u a u x x l lu t u l t B t Bu x x u x x l x x l l u x t v x t w x v a v w x x l lv t w ππππ⎧=+⎪⎪==≥⎨⎪⎪==−≤≤⎩=+′′=+++求下列问题的解解:设问题的解为 将其代入上面的定解问题,得22222)0,(,)(),(,0)(),(,0)().224sincos 0,(0)0().4()sin.8(0,)0,(,)0,(,0)t tt xx v l t w l B Bv x w x x v x x l x l a w x x l lw w l B B l w t x x l a l v a v v t v l t v x ππππ⎧⎪⎪=+=⎨⎪⎪+==−⎩⎧′′+=⎪⎨⎪==⎩=+==== 化成下面两个问题:(1) , 解得: (2) 12222022340(),(,0)().(,)cos sin sin .0,4;24sin sin 8, 4.824()sin t n n n l n l n Bx w x v x x l x l n a n a n v x t a t b t x l l l n l n a x xdx l l a l l n an l b x l x xdx n a l n ππππππππππ∞=⎧⎪⎪⎨⎪⎪−=−⎩⎛⎞=+⎜⎟⎝⎠≠⎧⎪=−⋅=⎨−=⎪⎩=−⋅=∑∫∫ 解得: 其中, ()()43222441222[11].4[11]44(,)cos sin sin sin .844(,)(,)()1cossin 8nn n al l a n a n v x t t x t x a l l n a l l B l a u x t v x t w x x t x l a l l πππππππππ∞=−−−−=−+⎛⎞=+=+−⎜⎟⎝⎠∑ 则 因此,原问题的解为14..0,(2)(-)(),(-)().0().:0X X X X X X X x Be Ae Be A B λππππλ′′+=⎧⎨′′==⎩<=++=+−=−==⇒求下列问题的固有值与固有函数解:当时,方程的通解为 由边界条件,有, ; 得0()0.0().-0.:().0().sin ,X x X x Ax B A B A B A X x C X x A B A B A Bλππλ===++=+⇒==>=+−=++=− 当时,方程的通解为 由边界条件,有 得当时,方程的通解为 由边界条件,有22sin ;()0sin 0(1,2,);()cos sin .(0,1,2,),()cos sin .n n n n n n n n X x n n X x A nx B nx n n X x A nx B nx λλ+====+===+"""" 要不恒等于,则,得故,固有值 固有函数222()()0,(3)(1)()0.ln ,()0.0()00:x y x xy x y y y e x e x d y y d y x Be Bx A B Be τλτλττλ′′′⎧++=⎨==⎩==+=<=+=++=+=解:方程通过自变量代换 或 得: 当时,方程的通解为 由边界条件,有 , ; 得))0()0.0()ln .0,0.:()0.0()cos ln sin ln .0,A B y x y x A B A x B B A y x y x A B A x B x A λτλ==⇒===+=+===>=+=+= 当时,方程的通解为 由边界条件,有 得当时,方程的通解为 由边界条件,有()()2220;()00(1,2,);()sin ln .(1,2,),()sin ln .n n n n n n B y x n n y x B n x n n y x B n x λππλπ========"""" 要不恒等于,则,得 故,固有值 固有函数。
数学物理方程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、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
3
类似于傅里叶级数, 固有函数系展开有下 (3) 类似于傅里叶级数,按固有函数系展开有下 面的收敛性 收敛性: 面的收敛性: 若函数 f (x)在 (a, b) 内有一阶连续导数及分段 连续的二阶导数,并且满足所给的边界条件, 连续的二阶导数,并且满足所给的边界条件, 内可以按固有函数展开为绝对且 则 f (x)在 (a, b) 内可以按固有函数展开为绝对且 一致收敛的级数 的级数: 一致收敛的级数:
这样的问题称为固有值问题。 这样的问题称为固有值问题。 固有值问题 也属于施图姆也属于施图姆-刘维尔问题 施图姆
1
施图姆施图姆-刘维尔方程的一般形式
d dy p ( x) − q ( x) y + λρ ( x) y = 0 dx dx
(95)
其中 1. p( x), p ′( x) ∈ C[a, b], p( x) > 0 (a < x < b); 2. q( x) ∈ C[a, b], 或者 q( x) ∈ C (a, b), 而在 区间端点处至多有一阶极点, 区间端点处至多有一阶极点,且 q( x) ≥ 0; 3. ρ ( x) ∈ C[a, b], ρ ( x) > 0. 方程(95)加上边界条件就称为施图姆方程(95)加上边界条件就称为施图姆-刘维尔问题 (95)加上边界条件就称为施图姆 那些使施 刘问题存在 存在非 那些使施-刘问题存在非0解的 λ 值,称为该问题 固有值,而相应于给定的固有值的非 的固有值,而相应于给定的固有值的非0解,称为 固有函数。 固有函数。
∫
e
1
1
x y n ( x) y m ( x)dx
t 作变换 x = e
= ∫ e 2t y n (t ) y m (t )dt
0
= ∫ sin nπt sin mπtdt = 1
0
1
0, m ≠ n,
2 ,
m = n.
15
固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。
{ (1)首先求出固有函数系 (1)首先求出固有函数系 y n (x)}的具体表达式 t t = ln x 作变换 x = e 则有
解
1 y x = yt ⋅ , x
1 1 1 1 1 y xx = ( y tt ⋅ ) ⋅ + y t ⋅ (− 2 ) = 2 y tt − 2 y t , x x x x x
f ( x) = ∑ c n y n ( x),
n =1 ∞
(97)
内是分段连续函数, 若函数 f (x), f ′(x) 在 (a, b) 内是分段连续函数, 则级数(97) (97)在 则级数(97)在 f (x) 的间断点 x 0 处收敛于
1 [ f ( x0 + 0) + f ( x0 − 0)], 2
齐次欧拉方程
x 2 y ′′ + xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0 1 固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。
解
{ (2)现在验证固有函数系 (2)现在验证固有函数系 y n (x)}的函数正交性
x 2 y ′′ + xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0 1 固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。
解
y tt + λy = 0
y (0) = y (1) = 0.
λ n = ( nπ ) 2 ,
7
X ′′( x ) + λX ( x ) = 0, X ′(0) = αX (0), X ′ | = 0. x =l
(1)边界条件中的 (1)边界条件中的 α = 0 时固有值问题简化为
X ′′( x) + λX ( x) = 0, X ′ | x =0 = 0, X ′ | x =l = 0.
∫
e
1
1 y n ( x) y m ( x)dx x
t 作变换 x = e
= ∫ y n (t ) y m (t )dt
0
1
= ∫ sin nπt sin mπtdt = 1
0
1
0, m ≠ n,
2 ,
m = n.
12
思考 试证问题
x 2 y ′′ + 3xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0
2.6
固有值与固有函数
在本章的前三节我们应用分离变量法求解弦振 动方程、 动方程、一维热传导方程和二维拉普拉斯方程的 有关定解问题时, 有关定解问题时,都需要解决一个含参变量 λ 的 常微分方程的边值问题, 常微分方程的边值问题,
X ′′( x) + λX ( x) = 0, X (0) = X ′(l ) = 0.
x 2 y ′′ + xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0 1 固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。
{ (1)首先求出固有函数系 (1)首先求出固有函数系 y n (x)}的具体表达式 t t = ln x 作变换 x = e 则有
此时对应的固有值和固有函数为
nπ 2 nπx λn = ( ) , X n ( x) = cos (n = 0, 1, 2, ⋯). l l
8
X ′′( x ) + λX ( x ) = 0, X ′(0) = αX (0), X ′ | = 0. x =l
(2)边界条件中的 (2)边界条件中的α → ∞ 时固有值问题简化为
x
(n = 1, 2, ⋯)
14
思考 试证问题
x 2 y ′′ + 3xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0
固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。 解
{ (2)现在验证固有函数系 (2)现在验证固有函数系 y n (x)}的函数正交性
2
关于固有值和固有函数的几点结论: 关于固有值和固有函数的几点结论: 存在无穷多个实的固有值: (1) 存在无穷多个实的固有值: λ1 ≤ λ 2 ≤ ⋯ ≤ λ n ≤ ⋯ ,
λ 当 q( x) ≥ 0 时, n ≥ 0 (n = 1, 2, 3, ⋯); 对应于这些固有值 有无穷多个固有函数: 有无穷多个固有函数:
6
下列定解问题: 13. 下列定解问题: u t = a 2 u xx (0 < x < l , t > 0), (u x − αu ) | x =0 = 0, u x | x =l = 0, u ( x,0) = ϕ ( x) 的固有值问题为 X ′′( x ) + λX ( x ) = 0, X ′′( x) + λX ( x) = 0, ( X ′ − αX ) | x =0 = 0, X ′(0) = αX (0), X ′ | = 0. X ′ | = 0. x =l x =l
代入原方程有
y tt − y t + 3 y t + λy = 0
y tt + 2 y t + λy = 0
13
思考 试证问题
x 2 y ′′ + 3xy ′ + λy = 0, (1 < x < e) y (1) = y (e) = 0
固有函数系 {y n (x)} 在 [1, e] 上带权函数 x 正交。 正交。 y tt + 2 y t + λy = 0 解
y1 ( x), y 2 ( x), ⋯, y n 的固有函数记为 y n (x), 那么所有y n (x)组成一个带权函数 ρ (x) 的正交函数 系,即 b (m ≠ n). (96) ∫a ρ ( x)y m ( x) y n ( x)dx = 0 例如
解
1 y x = yt ⋅ , x
1 1 1 1 1 y xx = ( y tt ⋅ ) ⋅ + y t ⋅ (− 2 ) = 2 y tt − 2 y t , x x x x x
代入原方程有
y tt − y t + y t + λy = 0
y tt + λy = 0
10
15. 试证问题
齐次欧拉方程
y (0) = y (1) = 0.
λ n = (nπ ) 2 + 1,
y n (t ) = Bn e − t sin nπt (n = 1, 2, ⋯).
1 将 t = ln x 代入即得 y n ( x) = Bn x sin(nπ ln x), (n = 1, 2, ⋯)
则原问题的固有函数系 1 为 {y n ( x)} = sin(nπ ln x)
f ( x) = ∑ c n y n ( x),
n =1 ∞
(97)
(n = 1, 2, 3, ⋯);
4
其中
cn
∫ =
b
a
ρ ( x) f ( x) y n ( x)dx
∫
b
a
ρ ( x) y ( x)dx
2 n
类似于傅里叶级数, 固有函数系展开有下 (3) 类似于傅里叶级数,按固有函数系展开有下 面的收敛性 收敛性: 面的收敛性: