弦振动偏微分方程的求解
(郑州航空工业管理学院数理系 田硕 450015)
摘要：本文列出了不同情况下的弦振动问题的定解方程及其成立条件，给出了不同情况下偏微分方程的求解方法，对于我们的生活和学习有一定的指导意义。 关键词：数学物理方程；偏微分方程；弦振动；拉普拉斯变换
Method for solving partial differential equations of string vibration （Tianshuo Department of mathematics and physics, Zhengzhou Institute of
Aeronautics Industry Management, henna zhengzhou 450015）
Abstract : This article lists the definite solution of the equation of string vibration problems in different situations and the establishment of conditions, given the method for solving partial differential equations under different circumstances, for our lives and learning have a certain significance. Keywords : mathematical physics equations; partial differential equations; vibrating string; Laplace transform
在数学物理方程中，根据常见物理模型，可以建立求解的偏微分方程。如在很多物理实际问题中要遇到的拉普拉斯方程，泊松方程，波动方程，热传导方程等等。对偏微分方程求解的讨论，有很重要的意义和运用。对不同的偏微分方程，往往有不同的求解方法，这要根据方程本身的特点而定。选取合适的方法不仅可以使问题简化，有时候也能体现出方程背后更深层次的物理意义。理想弦的振动方程就是一个一维波动方程的特例，本文将给出不同情况下的弦振动偏微分方程，并对它们的求解给予一定的讨论。
一、无界弦的自由振动问题
无界弦的自由振动问题既是满足下面条件的偏微分方程[1]
：
??
?+∞<<-∞==>+∞<<-∞=)
(),(),0(),(),0(),
0,(2x x x u x x u t x u a u t xx tt φ? 对于该偏微分方程，我们可以类似常微分方程初始问题的解法，先求出通解，然后把初始条件代入通解，以确定任意常数，从而求得初始问题的解。
做变量代换at x -=ξ，at x +=η，代入偏微分方程，整理可得：
02=???η
ξu
，得方程的通解为：)()()()(at x g at x f g f u ++-=+=ηξ 再代入初始条件，有：
??
?='+'-==+=)
2()
()()(),0()1()()()(),0(x x g a x f a x u x x g x f x u t φ? 对(2)式积分:
)3()(1)()(0c d a
x g x f x
+=
+-?λλφ
将（1）式和（3）式联立，解之则得：
2
)(212)
()(0c
d a x x f x -
-=?λλφ?
2
)(212)
()(0c
d a x x g x +
+=?λλφ? 于是我们便得到了：
?+-+
++-=
++-=at
x at
x d a at x at x at x g at x f x t u λλφ??)(212
)
()()()(),( 这便是一维无界弦的自由振动解的表达式， 称作达朗贝尔公式。由于对u 没有任何限制，只要一维波动方程有解，解必由达朗贝尔公式给出，且解是唯一的。
二、有界弦的自由振动问题。
描述两端固定的有界弦的自由振动的混合问题：
??
?
??====>+∞<<-∞=（初始条件）边界条件），
),(),0(),(),0((0),()0,()0,(2x x u x x u l t u t u t x u a u t xx tt φ? 对于该问题，适合用分离变量方法进行求解。 第一步，分离变量，分析求一族满足泛定方程和边界条件的分离变量形式的非零特解，可以先不估计初始条件。
令：)()(),(t T x X x t u =，把它代入方程，得
)()()()(2t T x X a t T x X ''=''
两边除以)()(2
t T x X a ，得
)
()
()()(2
x X x X t T a t T ''='' 此式左端仅是t 的函数，右端仅是x 的函数，而x 与t 是两个相互独立的变量，所以只有两边都是常数时，等式才能成立，令这个常数为λ-，就得到一个常微分方程：
02=+''T a T λ
及其边值问题（因,0)()0()0,(==t T X t u 所以0)0(=X ;同理,0)()(),(==t T l X l t u 所以0)(=l X ） 故第二个常微分方程是：
??
?==<<=+''0)(0)
0(0))(l X X l x x X x X ）
（（λ
第二步，解固有值问题
怎么找到满足条件的固有值λ，使常微分方程的边值问题有非零解。分三种情况讨论。 （1）0=λ,这时方程为：, 0=''X ,通解为：B Ax X +=,由边界条件
0)()0(==l X X ，得A=0;B=0，0)(≡x X ，不满足要求。
（2）0<λ，不妨设2
k -=λ，这时方程0))(2
=-''x X k x X （的通解为：
kx kx Be Ae X -+=由边界条件0)()0(==l X X ，得
???=+=+-0
kl
kl Be Ae B A 不难求出A=B=0，同样不满足要求。
（3）0>λ，不妨设2
k =λ（0>k ），这时方程0))(2
=+''x X k x X （的通解为：
kx B kx A X sin cos +=
由条件X(0)=0，知，A=0，再由条件0)(=l X ，得
0sin =kl B ，由于B 不能再为零，则必有),2,1( ==n l
n k π
或者：),2,1(2
=??
? ??=n l n πλ
我们把),2,1(2
=??
?
??=n l n n πλ叫做固有值，与固有值对应的非零解为：
l
x n B x X n n πsin
)(=，n B 是任意常数。求固有值和固有函数的边值问题称为固有值
问题。
把固有值2
???
??=l n n πλ代入确定T 的常微分方程：
l
at
n D l at n C t T n n n ππsin
cos
)(+=，n C ,n D 为任意常数。这样得到： ),2,1(sin
sin cos )()(),( =??
? ??
+==n l x n l at n D l at n C t T x X x t u n n n n n πππ 把n B 归入常数n C ,n D
第三步，写出级数形式解
由于方程和边界条件都是线性齐次的，故由叠加原理，级数:
∑∑+∞
=+∞
=??? ?
?
+==11sin
sin cos ),(),(n n n n n l x n l at n D l at n C x t u x t u πππ 仍满足方程和边界条件。
第四步，确定级数解中的系数n C 和n
D
由初始条件：∑+∞
==
=1
sin
),0()(n n
l
x
n C
x u x π?及 ∑
+∞
===1
sin ),0()(n n t l x
n D l a n x u x ππφ,由正弦展开的系数公式，得： ?=l n dx l x
n x l C 0sin )(2π?
?=l n dx l
x n x a n D 0sin )(2πφπ 这样我们得到该问题的定解为：
∑??+∞
=??? ???+?=100sin sin sin )(2cos sin )(2),(n l l
l x n l at n d l n a n l at n d l n l
x t u ππξπξξφππξπξξ?三、无界弦的受迫振动问题 该问题的偏微分方程为：
??
?==>+∞<<-∞+=（初始条件）),
(),0(),(),0()
0,(),(2x x u x x u t x t x f u a u t xx tt φ? 对该问题，用拉普拉斯变换计算比较方便[2]
。
对泛定方程施行拉普拉斯变换
dt e t x u x p u pt ?+∞
-=0
),(),(得：0),(20
02=----==p x f u a u u p u p xx t t
t
代入初始条件得：0),(2
2=----p x f u a p u p xx φ?该非齐次常微分方程的通解是
d p p f p
e e a d p p
f p
e e a Be
Ae
p x u x a
px a px x a px a px a
px a
px )](),()([21)](),()([21),()()(ξ?ξξφξξ?ξξφ+++++-+=??---
考虑到+∞→x 和-∞→x 时u 不应为无穷大，所以A=0，B=0，另为保证积分收敛，第一个积分下限取∞+，第二个积分下限取∞-。所以
?
?
????++??????++?
?????+=+++++=+++++-=??????????∞---∞+--∞---∞+--∞---∞+--∞---∞+--∞---∞--ξξ?ξξ?ξξξξξξφξξφξξ?ξξφξξ?ξξφξξ?ξξφξξ?ξξφξξξξξξξξξξd p p e a d p p e a d p f p e a d p f p e a d p e a d p e a d p p f p e a d p p f p e a d p p f p e a d p p f p
e a p x u x x p x a x p x a
x p x a x p x a x p x a x p x x p x a x p x x p x a
x p )(21)(21),(21),(21)(21)(21)](),()([21)](),()([21)](),()([21)](),()([21),()()()()()()()()()()(对于第一个中括号，运用延迟定理，
),(1
t H P
?则
???+>+<=--?--)
(0
)(1
)()(at x at x a x t H p e a x p ξξξξ
所以ξξφξξφξd a d p e a at
x x x a x p )(21)(21)(??+∞+--? 同理ξξφξφξd a
d p
e a x at x x a x p )(21)(21)(??-∞---? 对第三个中括号，)(ξ?代替了）（ξφ且多了一个因子p ，则对第一个中括号中原函数中）
（ξφ替换行为)(ξ?并对t 求导即得第三个中括号里的原函数分别为： ξξ?ξd p p
e a x x p )(21)(?∞+--)(21
at x a +?? ξξ?ξd p p e a x x p )(21)(?∞---)(21
at x a
-??
对第二个中括号，运用卷积定理
?-?t
d t f f p f p f 0
2121)()()()(τττ
τξτξτξξττξξξτξd d f a d d a
x
t H f a d p f p e a t t a x x
t x x a x p ?????-+∞+∞+--=---=0)(0)(),(21)(),(21),(21
同理：τξτξξξτξd d f a
d p f p
e a t x t a x x a x p ???--∞---=0)
)()(),(21),(21 于是得到该问题解的表达式为：
???+--+--++-++=at x at x t t a x t a x d d f a
d a at x at x t x u τξτξξξφ??ττ0)()(),(21)(21)]()([21),( 四、半无界弦的自由振动问题
该问题即求下面问题的解[3]
：
?
??
??===>+∞<≤=∞→（初始条件）
边界条件）有限，
,),0(,0),0((,0)0,()0,0(2b x u x u u t u t x u a u t x xx tt 对t 做拉普拉斯变换。设dt e t x u x p u pt ?
+∞
-=
),(),(，则有：
)0,()0,(),(2x u x pu p x u p u t tt --→
利用初始条件，上式变为：b p x u p u tt -→),(2
，原方程的像为：
0),(2
22
2
=??--x u a b p x u p ，解为：x a p
x a p Be Ae p b u -++=2。 A 、B 需要用边界条件予以确定。当对t 做拉普拉斯变换时候，像函数所满足的边界条件就是原函数边界条件的像。因此：有限∞→x u 即有限∞
→x u
。故A=0
又有边界条件0)0,(=t u ，可得0)0,(=p u ，代入可得2p
b
B -
=。于是得到 )1()1(),(2p
e p p b e p
b
p x u x a
p x a
p ---=
-=
，对该式做反演。 )(1t H p →，)(a
x
t H p e x a
p -→-（延迟定理） 于是按积分定理，得反演之后的函数表达式为：
dt a x t H b dt a x t H t H b t x u p x u t t
)](1[)]()([),(),(00??--=--=→
当a x t <时，则有00)(00=?=-??dt dt a
x
t H t t
当a x
t >时，则有a x t dt dt a x t H t a
x t -=?=-??1)(0
于是可写作
)()()(0
a
x t H a x t dt a x t H t
--=-? 所以方程的定解可写作
)()(),(a
x
t H a x t b bt t x u ---=
综上所述，我们总结了在不同情况下弦振动偏微分方程的求解，根据不同情况做出了不同的
弦振动偏微分方程求解方法。这对我们学习和巩固偏微分方程在物理学中的应用有很好的应用，便于我们深刻理解物理问题，也对我们的实际生活有一定的指导意义。
参考文献：
[1]严振军.数学物理方法[M].合肥：中国科学技术大学出版社，2005
[2]苗明川.数学物理方法全程导学及习题全解[M].北京：中国时代经济出版社，2011 [3]余恬.数学物理方法解题示例[M].北京：北京邮电大学出版社，2011
（注：可编辑下载，若有不当之处，请指正，谢谢!）