张瑞
中国科学技术大学数学科学学院
rui@ustc.edu.cn |
求无界弦的自由振动问题的解 \[\casefunc{ u_{tt}=a^2 u_{xx} , x\in\Real, t>0 \\ u(0,x)=\phi(x), \\ u_t(0,x)=\psi(x) , x\in\Real } \] 前面已经求得了通解 \[u=f(x-at)+g(x+at) \] $f$,$g$是任意的二次可微函数。把给定的初始条件代入,有 \[u(0,x)=f(x)+g(x)=\phi(x) \\ u_t(0,x)=-af'(x)+ag'(x)=\psi(x) \] 对第2个式子积分,有 \[-f(x)+g(x)=\dfrac1{a}\int_0^x\psi(\xi)d\xi+c \]
与第1个式子联立,可得 \[f(x)=\dfrac{\phi(x)}2-\dfrac1{2a}\int_0^x\psi(\xi)d\xi-\dfrac{c}2 \\ f(x)=\dfrac{\phi(x)}2+\dfrac1{2a}\int_0^x\psi(\xi)d\xi + \dfrac{c}2 \] 这样,解为 \[u(t,x)=f(x-at)+g(x+at) \\ =\dfrac{\phi(x-at)+\phi(x+at)}2+\dfrac1{2a}\int_{x-at}^{x+at}\psi(\xi)d\xi \] 这个公式,称为d'Alembert公式
求解一个定解问题,一般分为3步:
解释:
一维波动方程的一般解
\[u=f(x-at)+g(x+at) \] 表示左行波与右行波
d'Alembert公式中,
若$\phi(x), \psi(x)$为奇函数时。则 \[u(t,0)=\dfrac12[\phi(-at)+\phi(-at)]+\dfrac1{2a}\int_{-at}^{at}\psi(\xi)d\xi=0 \]
若$\phi(x), \psi(x)$为偶函数时。则 \[u_x(t,0)=\dfrac12[\phi'(-at)+\phi'(-at)]+\dfrac1{2a}(\psi(at)-\psi(-at)) \\ =0 \]
例 1. 求半无界弦的自由振动问题的解 \[\casefunc{ u_{tt}=a^2 u_{xx} , x>0, t>0 \\ u(0,x)=\phi(x), u_t(0,x)=\psi(x) , x>0 \\ u(t,0)=0 } \]
1 不能直接用d'Alembert公式
注意到奇函数满足边界条件,所以用奇函数延拓
\[\Phi(x)=\casefunc{ \phi(x), x\geq 0 \\ -\phi(x), x<0 } \]
\[\Psi(x)=\casefunc{ \psi(x), x\geq 0 \\ -\psi(x), x<0 } \] 得到解为
\[u(t,x)=\dfrac12[\Phi(x+at)+\Phi(x-at)]+\dfrac1{2a}\int_{x-at}^{x+at}\Psi(\xi)d\xi \] 所以,
例 2. 球对称三维波动方程
\[\casefunc{ u_{tt}=a^2\triangle_3 u , t>0 \\ u(0,r)=\phi(r) , r=\sqrt{x^2+y^2+z^2} \\ u_t(0,r)=\psi(r) } \]
例2. 解:
若初始函数 $\phi(x), \psi(x)$达不到光滑性要求,可以用Weierstrass逼近定理得到极限解。
例 3. 解方程
\[\casefunc{ & u_{xx}-u_{yy}+\cos x=0 , y>0, x\in\Real \\ & u(x,0)=0 \\ & u_y(x,0)=4x } \]
3.