弦振动偏微分方程的求解

弦振动偏微分方程的求解(郑州航空工业管理学院数理系 田硕 450015)摘要：本文列出了不同情况下的弦振动问题的定解方程及其成立条件，给出了不同情况下偏微分方程的求解方法，对于我们的生活和学习有一定的指导意义。

关键词：数学物理方程；偏微分方程；弦振动；拉普拉斯变换Method for solving partial differential equations of string vibration（Tianshuo Department of mathematics and physics, Zhengzhou Institute ofAeronautics 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 atx 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 ''=''两边除以)()(2t T x X a ，得 )()()()(2x 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<λ，不妨设2k -=λ，这时方程0))(2=-''x X k x X （的通解为：kx kx Be Ae X -+=由边界条件0)()0(==l X X ，得⎩⎨⎧=+=+-00kl kl Be Ae B A 不难求出A=B=0，同样不满足要求。

（3）0>λ，不妨设2k =λ（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 ln 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 由初始条件：∑+∞===1sin),0()(n n lx n C x u x πϕ及∑+∞===1sin ),0()(n n t l x n D l a n x u x ππφ,由正弦展开的系数公式，得： ⎰=l n dx lx n x l C 0sin )(2πϕ ⎰=l n dx lx 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 lx 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),(2002=----==p x f u a u u p u p xx t t t 代入初始条件得：0),(22=----p x f u a p u p xx φϕ该非齐次常微分方程的通解是dp p f p e e a d p p f p e e a Be Aep 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 pe a d p pf pe a d p pf p e a d p p f pe a p x u x x p x a x p x ax p x a x p x ax p x a x p x x p x ax p x x p x ax p )(21)(21),(21),(21)(21)(21)](),()([21)](),()([21)](),()([21)](),()([21),()()()()()()()()()()(对于第一个中括号，运用延迟定理，),(1t 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 xx a x p )(21)(21)(⎰⎰+∞+--⇔ 同理ξξφξφξd ad pe a x at x x a x p )(21)(21)(⎰⎰-∞---⇔ 对第三个中括号，)(ξϕ代替了）（ξφ且多了一个因子p ，则对第一个中括号中原函数中）（ξφ替换行为)(ξϕ并对t 求导即得第三个中括号里的原函数分别为： ξξϕξd p pe a x x p )(21)(⎰∞+--)(21at x a +⇔ϕ ξξϕξd p p e a x x p )(21)(⎰∞---)(21at x a-⇔ϕ 对第二个中括号，运用卷积定理⎰-⇔t d t f f p f p f 02121)()()()(τττ τξτξτξξττξξξτξd d f a d d a x t H f a d p f p e a t t a x xt x x a x p ⎰⎰⎰⎰⎰-+∞+∞+--=---=0)(0)(),(21)(),(21),(21 同理：τξτξξξτξd d f ad p f pe 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 做拉普拉斯变换。