樊龍

(山西大同大學(xué)煤炭工程學(xué)院，山西大同037000)

一維波動(dòng)方程混合問題的通解

樊龍

(山西大同大學(xué)煤炭工程學(xué)院，山西大同037000)

一維情形下波動(dòng)方程的混合問題(初邊值問題)是一類重要的物理模型，常用求解方法是波的反射原理，計(jì)算特征線在邊界上的反射次數(shù)得出問題的解，但是弊端在于計(jì)算量大，且沒有通用的求解公式，并不能反映出波的反射實(shí)質(zhì)，另一種方法是Fourier級(jí)數(shù)法，利用分離變量將原方程化為常微分方程組，再利用常微分方程特征理論得出級(jí)數(shù)解，同樣不易計(jì)算。為了簡(jiǎn)化計(jì)算過程，先對(duì)初值條件φ(x),ψ(x)根據(jù)邊值條件進(jìn)行相應(yīng)的奇偶性延拓，可將原問題化簡(jiǎn)為初值問題，由D’Alambert公式給出問題在R上的解，再將問題的全局解限定在有限區(qū)間[0,l]上得出通解公式，結(jié)果具有一般性。

波動(dòng)方程；D’Alambert公式；延拓

引言

本文主要討論以下兩類混合問題，

(1)

邊值條件為

ux(0,t)=0，ux(l,t)=0，t≥0

(2)

或者邊值條件為

u(0,t)=0，ux(l,t)=0，t≥0

(3)

對(duì)于以上問題求解，常用的方法有波的反射原理[1]以及Fourier級(jí)數(shù)方法(分離變量法)[2-3]。文獻(xiàn)[4]利用特征線法，給出方程在可解區(qū)域內(nèi)的通解公式，提供了另一種計(jì)算方法，文獻(xiàn)[5]中提及了通過將函數(shù)延拓，將求解區(qū)域進(jìn)行劃分，然后逐個(gè)進(jìn)行討論，利用波的反射原理，給出各個(gè)可解區(qū)域內(nèi)解的顯示表達(dá)，最近，在文獻(xiàn)[6]中，利用分離變量方法求解了幾類特殊的波動(dòng)方程，在文獻(xiàn)[7]中利用Laplace變換對(duì)有限一維空間彈性動(dòng)力學(xué)邊值問題給出解的嚴(yán)格推導(dǎo)。

以往的工作都集中在具體的求解過程及方法，并未給出一個(gè)確切的解的表達(dá)式，本文的工作完善了波動(dòng)方程在通解公式方面的內(nèi)容，并且采用不同于以往的方法，直接利用D’Alambert公式計(jì)算，從而避開各區(qū)域討論的繁瑣過程，直接給出方程在各個(gè)區(qū)域的解的顯示表達(dá)，且結(jié)果具有一般性。

1 主要結(jié)論

引理1對(duì)于問題(1)，若初始函數(shù)φ(x),ψ(x)都是關(guān)于x0的奇函數(shù)，則方程的解在任何時(shí)間都有u(x0,t)=0。

證明由達(dá)朗貝爾公式，方程的解為

由于φ(x),ψ(x)都是關(guān)于x0的奇函數(shù)，有

φ(x0+C)=-φ(x0-C)，ψ(x0+C)=-ψ(x0-C)

令τ=2x0-ξ

證畢。

引理2對(duì)于問題(1)，若初始函數(shù)φ(x),ψ(x)都是關(guān)于x0的偶函數(shù)，則方程的解在任何時(shí)間都有ux(x0,t)=0。

證明

φ(x)為關(guān)于x0偶函數(shù)，則φ′(x)為關(guān)于x0的奇函數(shù)，所以u(píng)x(x0,t)=0，證畢。

由引理1和引理2，問題(1)、(2)以及(1)、(3)可對(duì)初始函數(shù)φ(x),ψ(x)進(jìn)行相應(yīng)的延拓，再利用達(dá)朗貝爾公式得出通解。

定理1對(duì)于問題(1)、(2)，對(duì)函數(shù)φ(x),ψ(x)對(duì)于點(diǎn)x=0,x=l均作偶延拓，即

(-1)m-1(m+δm)l]+φ[(-1)n(x-at)+

(-1)n-1(n+δn)l]}+

定理2對(duì)于問題(1)、(3)，對(duì)函數(shù)φ(x),ψ(x)關(guān)于x=0作奇延拓，對(duì)函數(shù)φ(x),ψ(x)關(guān)于x=l作偶延拓即

u(x,t)=

其中δk定義同定理1。

2 結(jié)論證明

定理1的證明：首先問題(1)、(2)的解可用達(dá)朗貝爾公式表示

(4)

由引理1、引理2以及Φ(x)，Ψ(x)均為關(guān)于x=1,x=l的偶函數(shù)，可知ux(0,t)=0，ux(l,t)=0，滿足邊值條件。

由Φ(x),Ψ(x)的定義，容易得到

(5)

(6)

分析積分部分，將問題分情況討論：

(1)n為偶數(shù)，m為偶數(shù)

由于Ψ(x)是關(guān)于x=0,l的偶函數(shù)，所以對(duì)于k∈Ζ有

進(jìn)而

綜合

(2)n為奇數(shù)，m為偶數(shù)

同理，

(3)n為偶數(shù)，m為奇數(shù)

(4)n為奇數(shù)，m為奇數(shù)

所以

(7)

將式(5)、式(7)代入式(4)，定理1證畢。

定理2的證明：方法類似于定理1，由引理1、引理2以及Φ(x)，Ψ(x)的定義，可知u(0,t)=0,ux(l,t)=0；易知

(8)

(9)

(1)n為偶數(shù)，m為偶數(shù)

由(9)式可知

因?yàn)棣?x)是關(guān)于2kl(k∈Z)的奇函數(shù)，且m,n均為偶數(shù)，所以

同理

所以可得

(2)n為奇數(shù)，m為偶數(shù)

同理，

(3)n為偶數(shù)，m為奇數(shù)

(4)n為奇數(shù)，m為奇數(shù)

綜合可得

(10)

將式(8)、式(10)代入式(4)，定理2證畢。

推論1當(dāng)邊值條件為u(0,t)=0，u(l,t)=0，t≥0，解為

u(x,t)=

(-1)nφ[(-1)n(x-at)+(-1)n-1(n+δn)l]}+

推論2當(dāng)邊值條件為ux(0,t)=0，u(l,t)=0，t≥0，解為

u(x,t)=

證明過程類似定理1，定理2，在此省略。

3 實(shí)例

求解如下混合問題

其他區(qū)域同理，這也體現(xiàn)了在不同區(qū)域，由于沿各邊界反射次數(shù)不同，所以得到不同形式的解，反映出波的反射原理的實(shí)質(zhì)。

[1]朱長(zhǎng)江,鄧引斌.偏微分方程教程.北京:科學(xué)出版社,2005.

[2]谷超豪,李大潛,陳恕行,等.數(shù)學(xué)物理方程.2版.北京:高等教育出版社,1996.

[3]陳恕行,秦鐵虎,周憶.數(shù)學(xué)物理方程.上海:復(fù)旦大學(xué)出版社,2003.

[4]姜玲玉.關(guān)于波動(dòng)方程混合問題的特征線方法.數(shù)學(xué)雜志,2004,24(5):577-580.

[5]趙天玉,毛占軍.求解波動(dòng)方程混合問題的通解函數(shù)延拓法.長(zhǎng)江大學(xué)學(xué)報(bào):自然科學(xué)版,2005,2(1):4-9.

[6]王心平.分離變量法在求解波動(dòng)方程中的應(yīng)用.科技視界,2014,35:60-61.

[7]趙福垚,宋二祥.有限一維彈性波動(dòng)問題的公理化求解及其在自由場(chǎng)中的應(yīng)用.工程力學(xué),2015(4):47-53.

[8]杜心華.一類非線性波動(dòng)方程混合問題整體解的存在唯一性.四川師范大學(xué)學(xué)報(bào):自然科學(xué)版,1994(4):10-14.

[9]楊萌.可控一維波動(dòng)方程的邊值混合問題.湖北理工學(xué)院學(xué)報(bào),2008(6):55-58.

[10]梁志輝,李之杰.一維波動(dòng)方程混合問題的差分解.內(nèi)蒙古民族大學(xué)學(xué)報(bào),2009(2):5-6.

[11]趙天玉.求解波動(dòng)方程混合問題的曲線積分法.長(zhǎng)沙大學(xué)學(xué)報(bào),2005(2):1-3.

[12]林瓊桂.高維波動(dòng)方程與熱傳導(dǎo)方程非齊次邊界條件的一般處理.大學(xué)物理,2015(5):1-4.

[13]WALTER A S.Partial Differential Equations-An Introduction.New York:John Wiley&Sons Inc,1992.

[14]RAUCH J.Partial Differential Equations.New York:Springer-Verlag,1991.

[15]洪潔.一類非線性波動(dòng)方程混合問題解的爆破.蘭州理工大學(xué)學(xué)報(bào),2004(3):124-126.

The General Solution to Mixed Problem of a 1D Wave Equation

FANLong

(School of Coal Engineering, Shanxi Datong University, Datong 037000, China)

The mixed problem of 1D wave equation is an important model in physics. The most commonly used method is the principle of reflection wave, which calculates the number of reflection on the boundary to obtain the solution of the problem, but the disadvantage of this method is that the scale of calculation is large and general formula is not given, at the same time, the essence of reflection is not reflected. Another method is the Fourier series which is also hard to calculate. With method of Fourier series, ordinary differential equations can be obtained by separation of variables, then the problem is solved by using characteristic method. In the purpose of simplifying the calculation, taking corresponding extension of initial valueφ(x),ψ(x)accordingtotheboundaryvalue,theproblemcanbesimplifiedintotheformofinitialproblem,thenthesolutioninRcanbegotbyusingD’Alambert’sformula,andthegeneralsolutionistherestrictionofglobalsolutionon[0,l],meanwhiletheresultisgeneral.

wave equation; D’Alambert’s formula; extensions

2016-11-27

山西大同大學(xué)教學(xué)改革項(xiàng)目(XJY2013211)

樊 龍(1989-)，男，山西忻州人，助教，碩士，主要從事偏微分方程方面的研究，(E-mail)fanlongmath@163.com

1673-1549(2017)01-0092-05

10.11863/j.suse.2017.01.18

O175.27

A