偏微分方程(Partial Differential Equation I)

数学物理方法

  • 偏微分方程的定解问题
    • 定解问题及适定性
    • 一阶线性偏微分方程的通解
    • 二阶线性偏微分方程的分类和标准式
    • 波动方程的行波解
    • 线性叠加原理
  • 分离变量法
    • 齐次边界条件下齐次方程的定解问题
    • 齐次边界条件下非齐次方程的定解问题
    • 非齐次边界条件的齐次化
    • 施图姆-刘维尔本征值问题
  • 附录
    • 热传导方程导出


偏微分方程(Partial Differential Equation I)
偏微分方程(Partial Differential Equation II)
偏微分方程(Partial Differential Equation III)


参考文献:

《数学物理方程》| 季孝达
《数学物理方法》| 吴崇试
《数学物理方法》| 梁昆淼
MOOC北京大学《数学物理方法》| 吴崇试 、高春媛

偏微分方程的定解问题

偏微分方程 (Partial Differential Equations, PDE) ,是指含有多元未知函数及其若干阶偏导数的微分方程。
F ( u , x 1 , x 2 , ⋯   , ∂ u ∂ x 1 , ∂ u ∂ x 2 , ⋯   , ∂ 2 u ∂ x 1 2 , ⋯   ) = 0 (1.1) F(u,x_1,x_2,\cdots,\cfrac{∂u}{∂x_1},\cfrac{∂u}{∂x_2},\cdots,\cfrac{∂^2u}{∂x_1^2},\cdots)=0\tag{1.1} F(u,x1,x2,,x1u,x2u,,x122u,)=0(1.1)
其中最高阶偏导数的阶称为偏微分方程的阶。

在此引入一些重要的微分算子

Hamilton 算子: ∇ = ( ∂ ∂ x , ∂ ∂ y , ∂ ∂ z ) ∇=(\cfrac{∂}{∂x},\cfrac{∂}{∂y},\cfrac{∂}{∂z}) =(x,y,z)
Laplace算子 Δ = ∂ 2 ∂ x 2 + ∂ 2 ∂ y 2 + ∂ 2 ∂ z 2 = ∇ ⋅ ∇ = ∇ 2 Δ=\cfrac{∂^2}{∂x^2}+\cfrac{∂^2}{∂y^2}+\cfrac{∂^2}{∂z^2}=∇\cdot∇=∇^2 Δ=x22+y22+z22==2

定解问题及适定性

数学物理方程 (Equation of Mathematical Physics) ,通常指从物理学及其他各门自然科学、技术科学中产生的偏微分方程。数学物理方程有时还包括常微分方程和积分方程。本课程将着重讨论三类基本的二阶线性偏微分方程。

(1) 弦的横振动方程:未知函数 u ( x , t ) u(x,t) u(x,t) 表示坐标 x x x 处弦的横向位移
u t t − a 2 u x x = f ( x , t ) u_{tt}-a^2u_{xx}=f(x,t) utta2uxx=f(x,t)
其中 a = T ρ , f = F ρ a=\sqrt{\cfrac{T}{\rho}},f=\cfrac{F}{\rho} a=ρT ,f=ρF ρ \rho ρ 为弦的线密度, T T T 为弦的切应力, F F F 表示单位长度受到的横向力。

杆的纵振动方程:未知函数 u ( x , t ) u(x,t) u(x,t) 表示截面相对于平衡位置 x x x 处的纵向位移
u t t − a 2 u x x = f ( x , t ) u_{tt}-a^2u_{xx}=f(x,t) utta2uxx=f(x,t)
其中 a = E ρ , f = F ρ a=\sqrt{\cfrac{E}{\rho}},f=\cfrac{F}{\rho} a=ρE ,f=ρF ρ \rho ρ 为杆的线密度, E E E 为杆的Young模量, F F F 表示单位长度受到的纵向力

更一般的,三维空间的波动方程 (wave equation)是
u t t − a 2 Δ u = f ( r , t ) (1.2) u_{tt}-a^2Δu=f(\mathbf r,t)\tag{1.2} utta2Δu=f(r,t)(1.2)

(2) 热传导方程:(heat equation)未知函数 u ( r , t ) u(\mathbf r,t) u(r,t) 表示温度
u t − a 2 Δ u = f ( r , t ) (1.3) u_t-a^2 Δu=f(\mathbf r,t) \tag{1.3} uta2Δu=f(r,t)(1.3)
其中 κ ρ c , F ρ c \cfrac{\kappa}{\rho c},\cfrac{F}{\rho c} ρcκ,ρcF κ \kappa κ 称为热传导系数, ρ \rho ρ 是介质的密度, c c c 是比热容,与介质的质料有关。 F ( x , y , z , t ) F(x,y,z,t) F(x,y,z,t) 表示介质中有热源,单位时间内在单位体积介质中产生的热量

(3) 如果热传导方程热源不随时间变化,当温度达到稳定时 ( u t = 0 ) (u_t=0) (ut=0) ,温度分布满足 Poisson方程
Δ u = − f / a 2 (1.4) Δu=-f/a^2\tag{1.4} Δu=f/a2(1.4)
如果没有热源 ( f = 0 ) (f=0) (f=0) ,则有 Laplace方程(也称调和方程)
Δ u = 0 (1.5) Δu=0\tag{1.5} Δu=0(1.5)
静电场的电势 u ( r ) u(\mathbf r) u(r)也满足 Poisson方程
Δ u = − ρ / ϵ 0 Δu=-\rho/\epsilon_0 Δu=ρ/ϵ0
其中 ρ \rho ρ 为电荷密度, ϵ 0 \epsilon_0 ϵ0 为真空介电常数。

(4) 如果波动方程 u ( r , t ) u(\mathbf r,t) u(r,t) 随时间周期性变化,频率为 ω \omega ω,则 u ( r , t ) = v ( r ) e − i ω t u(\mathbf r,t)=v(\mathbf r)e^{-i\omega t} u(r,t)=v(r)eiωt,则 v ( r ) v(\mathbf r) v(r) 满足 Helmholtz 方程
Δ v + k 2 v = 0 (1.6) Δv+k^2v=0\tag{1.6} Δv+k2v=0(1.6)
其中 k = ω / a k=\omega/a k=ω/a 称为波数。

通解和特解:如果多元函数 u u u 具有偏微分方程中出现的各阶连续偏导数,并使方程 恒成立,则称此函数为方程的解,也称古典解。和常微分方程类似,称 m m m 阶偏微分方程的含有 m m m 个任意函数的解为通解。通解中的任意函数一旦确定便成为特解。

示例 1:求解偏微分方程
u x y = 0 u_{xy}=0 uxy=0
解:先关于 y y y 积分,得 u x = f ( x ) u_x=f(x) ux=f(x),再关于 x x x 积分,就得到通解
u = ∫ f ( ξ ) + f 2 ( η ) = f 1 ( ξ ) + f 2 ( η ) u =\int f(ξ)+f_2(η)=f_1(ξ)+f_2(η) u=f(ξ)+f2(η)=f1(ξ)+f2(η)
其中 f 1 , f 2 f_1,f_2 f1,f2 是任意函数。

示例 2:二维Laplace方程的通解
u x x + u y y = 0 u_{xx}+u_{yy}=0 uxx+uyy=0
引入变换
{ ξ = x + i y η = x − i y \begin{cases} ξ=x+iy \\ η=x-iy \end{cases} {ξ=x+iyη=xiy
根据复合函数偏导法则,进一步求得
{ u x x = u ξ ξ + 2 u ξ η + u η η u y y = − ( u ξ ξ − 2 u ξ η + u η η ) \begin{cases} u_{xx}=u_{ξξ}+2u_{ξη}+u_{ηη} \\ u_{yy}=-(u_{ξξ}-2u_{ξη}+u_{ηη}) \end{cases} {uxx=uξξ+2uξη+uηηuyy=(uξξ2uξη+uηη)
原方程变为
u ξ η = 0 u_{ξη}=0 uξη=0
于是可求得通解
u ( x , y ) = f ( x + i y ) + g ( x − i y ) u(x,y)=f(x+iy)+g(x-iy) u(x,y)=f(x+iy)+g(xiy)

定解条件:通常把反应系统内部作用导出的偏微分方程称为泛定方程。为了完全描写一个实际物理问题的确定解,我们需要在一定的制约条件,称为定解条件。偏微分方程和相应的定解条件就构成一个定解问题。常见的定解条件有以下几类:

  • 初始条件:应该完全描写初始时刻( t = 0 t=0 t=0)介质内部及边界上任意一点的状况。一般的讲,关于时间 t t t m m m 阶偏微分方程,要给出 m m m 个初始条件才能确定一个特解。
    对于波动方程来说,就是初始时刻的位移和速度
    u ∣ t = 0 = ϕ ( r ) , ∂ u ∂ t ∣ t = 0 = ψ ( r ) u|_{t=0}=\phi(\mathbf r),\quad \cfrac{∂u}{∂t}|_{t=0}=\psi(\mathbf r) ut=0=ϕ(r),tut=0=ψ(r)
    对于热传导方程,只用给出初始时刻物体温度的分布状态
    u ∣ t = 0 = ϕ ( r ) u|_{t=0}=\phi(\mathbf r) ut=0=ϕ(r)
    Laplace 方程和Possion 方程都是描述稳恒状态的,与时间无关,所以不需要初始条件。

  • 边界条件:边界条件形式比较多样化,要由具体问题中描述的具体状态决定。总的原则是:边界条件应该完全描写边界上各点在任意时刻的状况。
    对于弦的横振动问题,约束条件通常有以下三种

  1. 如果弦的两端固定,那么边界条件就是
    u ∣ x = 0 = 0 , u ∣ x = l = 0 u|_{x=0}=0,\quad u|_{x=l}=0 ux=0=0,ux=l=0

  2. 如果一端 ( x = 0 ) (x=0) (x=0) 固定,另一端 ( x = l ) (x=l) (x=l) 受位移方向的外力 F ( t ) F(t) F(t) ,可以推导出
    u ∣ x = 0 = 0 , ∂ u ∂ x ∣ x = l = 1 T F ( t ) u|_{x=0}=0,\quad \cfrac{∂u}{∂x}|_{x=l}=\cfrac{1}{T}F(t) ux=0=0,xux=l=T1F(t)
    如果一端固定,另一端外力为 0,即另一端自由,则
    u ∣ x = 0 = 0 , ∂ u ∂ x ∣ x = l = 0 u|_{x=0}=0,\quad \cfrac{∂u}{∂x}|_{x=l}=0 ux=0=0,xux=l=0

  3. 如果外力是由弹簧提供的弹性力,即 F ( t ) = − k u ( l , t ) F(t)=-ku(l,t) F(t)=ku(l,t),其中 k k k 为弹性系数,则
    ( ∂ u ∂ x + σ u ) ∣ x = l = 0 , σ = k T (\cfrac{∂u}{∂x}+\sigma u)|_{x=l}=0,\quad \sigma=\cfrac{k}{T} (xu+σu)x=l=0,σ=Tk

    对于热传导问题,也有类似的边界条件(以 ∂ Ω ∂\Omega Ω 表示区域 Ω \Omega Ω 的边界)

  4. 如果边界温度分布已知,则
    u ∣ r ∈ ∂ Ω = f ( r , t ) u|_{\mathbf r\in∂\Omega}=f(\mathbf r,t) urΩ=f(r,t)

  5. 如果物体边界和周围介质绝热,则
    ∂ u ∂ n ∣ r ∈ ∂ Ω = 0 \cfrac{∂u}{∂n}|_{\mathbf r\in∂\Omega}=0 nurΩ=0

  6. 设周围介质的温度为 u 1 ( r , t ) u_1(\mathbf r,t) u1(r,t) ,物体通过边界与周围有热量交换,则
    ( ∂ u ∂ n + σ u ) ∣ r ∈ ∂ Ω = σ u 1 (\cfrac{∂u}{∂n}+\sigma u)|_{\mathbf r\in∂\Omega}=\sigma u_1 (nu+σu)rΩ=σu1
    其中 σ = h / k , h \sigma=h/k,h σ=h/k,h表示两种物质间的热交换系数。

    常见的边界条件数学上分为三类(以 ∂ Ω ∂\Omega Ω 表示区域 Ω \Omega Ω 的边界, f f f 为已知函数)

  7. 第一类边界条件 (Dircichlet条件):直接给出边界上各点未知函数 u u u 的值
    u ∣ ∂ Ω = f ( r , t ) u|_{∂\Omega}=f(\mathbf r,t) uΩ=f(r,t)

  8. 第二类边界条件 (Neu-mann条件):给出边界外法线方向上方向导数的数值
    ∂ u ∂ n ∣ ∂ Ω = f ( r , t ) \cfrac{∂u}{∂n}|_{∂\Omega}=f(\mathbf r,t) nuΩ=f(r,t)

  9. 第三类边界条件 (Robin条件):给出边界上各点的函数值与外法线方向上方向导数的某一线性组合值
    ( ∂ u ∂ n + σ u ) ∣ ∂ Ω = f ( r , t ) (\cfrac{∂u}{∂n}+\sigma u)|_{∂\Omega}=f(\mathbf r,t) (nu+σu)Ω=f(r,t)

  • 连接条件:从微分方程的推导知道,方程只在空间区域的内部成立,如果在区域内部出现结构上的跃变,那么在这些跃变点(线,面)上微分方程不再成立,因此需要补充上相应的条件,通常称为连接条件
    例如,由两段不同材质组成的弦在 x 0 x_0 x0 处连接,波动方程需分段讨论,且在连接处位移和张力相等,连接条件为 u 1 ( x 0 , t ) = u 2 ( x 0 , t ) , ∂ u 1 ∂ x ∣ x = x 0 − 0 = ∂ u 2 ∂ x ∣ x = x 0 + 0 u_1(x_0,t)=u_2(x_0,t),\cfrac{∂u_1}{∂x}|_{x=x_0-0}=\cfrac{∂u_2}{∂x}|_{x=x_0+0} u1(x0,t)=u2(x0,t),xu1x=x00=xu2x=x0+0

定解问题的适定性:定解问题来自于实际,它的解也应该回到实际中去。如果一个定解问题的解存在、唯一且稳定,则称该定解问题是适定的 (well-posed) 。稳定性指的是,如果定解条件的数值有细微的改变,解的数值也作细微的改变。

一阶线性偏微分方程的通解

引例:求解右行单波方程的通解
u t + a u x = 0 u_{t}+au_{x}=0 ut+aux=0
引入变换
{ ξ = x − a t η = x \begin{cases} ξ=x-at \\ η=x \end{cases} {ξ=xatη=x
根据复合函数偏导法则,进一步求得
{ u t = − a u ξ u x = u ξ + u η \begin{cases} u_{t}=-au_{ξ}\\ u_{x}=u_{ξ}+u_{η} \end{cases} {ut=auξux=uξ+uη
原方程变为
u η = 0 u_{η}=0 uη=0
于是可求得通解
u ( x , t ) = g ( x − a t ) u(x,t)=g(x-at) u(x,t)=g(xat)

两个自变量一阶线性偏微分方程

先研究两自变量 x , y x,y x,y 的一阶线性偏微分方程
a u x + b u y + c u + f = 0 (2.1) au_x+bu_y+cu+f=0\tag{2.1} aux+buy+cu+f=0(2.1)
其中 a , b , c , f a,b,c,f a,b,c,f 都是自变量 x , y x,y x,y 的函数。
(1) 若 a ≡ 0 , b ≠ 0 a\equiv0,b\neq0 a0,b=0 方程改写为
u y + c b u + f b = 0 u_y+\cfrac{c}{b}u+\cfrac{f}{b}=0 uy+bcu+bf=0
利用一阶线性常微分方程的求解方法可求得通解
u = e − p ( x , y ) [ ∫ e p ( x , y ) f ( x , y ) b ( x , y ) d y + g ( x ) ] u=e^{-p(x,y)}[\int e^{p(x,y)}\cfrac{f(x,y)}{b(x,y)} \mathrm dy+g(x)] u=ep(x,y)[ep(x,y)b(x,y)f(x,y)dy+g(x)]
其中 p ( x , y ) = ∫ c ( x , y ) b ( x , y ) d y , g ( x ) p(x,y)=\displaystyle\int\cfrac{c(x,y)}{b(x,y)} \mathrm dy,g(x) p(x,y)=b(x,y)c(x,y)dy,g(x) 是任意函数。
(2) 若 a ( x , y ) b ( x , y ) ≠ 0 a(x,y)b(x,y)\neq0 a(x,y)b(x,y)=0 不能直接积分求解。我们希望通过适当的自变量变换和未知函数的变换,使方程简化,并在此基础上对方程求解。

作自变量非奇异变换(可逆)
{ ξ = ξ ( x , y ) η = η ( x , y ) \begin{cases} ξ=ξ(x,y) \\ η=η(x,y) \end{cases} {ξ=ξ(x,y)η=η(x,y)
其中雅可比行列式
∣ J ( ξ , η ) ∣ = ∣ ∂ ( ξ , η ) ∂ ( x , y ) ∣ = ∣ ξ x ξ y η x η y ∣ ≠ 0 |J(ξ,η)|=|\cfrac{∂(ξ,η)}{∂(x,y)}|=\begin{vmatrix}ξ_x&ξ_y \\ η_x&η_y\end{vmatrix}\neq0 J(ξ,η)=(x,y)(ξ,η)=ξxηxξyηy=0
以保证新变量 ξ , η ξ,η ξ,η 相互独立,利用链式法则
u x = ∂ u ∂ ξ ∂ ξ ∂ x + ∂ u ∂ η ∂ η ∂ x , u y = ∂ u ∂ ξ ∂ ξ ∂ y + ∂ u ∂ η ∂ η ∂ y u_x=\cfrac{∂u}{∂ξ}\cfrac{∂ξ}{∂x}+\cfrac{∂u}{∂η}\cfrac{∂η}{∂x},\quad u_y=\cfrac{∂u}{∂ξ}\cfrac{∂ξ}{∂y}+\cfrac{∂u}{∂η}\cfrac{∂η}{∂y} ux=ξuxξ+ηuxη,uy=ξuyξ+ηuyη
可得到新的方程
A u ξ + B u η + C u + F = 0 (2.2) Au_ξ+Bu_η+Cu+F=0\tag{2.2} Auξ+Buη+Cu+F=0(2.2)
其中
{ A = a ξ x + b ξ y B = a η x + b η y C = c , F = f (2.3) \begin{cases} A=aξ_x+bξ_y \\ B=aη_x+bη_y \\ C=c,F=f \end{cases} \tag{2.3} A=aξx+bξyB=aηx+bηyC=c,F=f(2.3)
新方程仍然是线性的。
由上式看出,如果取一阶偏微分方程
a z x + b z y = 0 (2.4) az_x+bz_y=0\tag{2.4} azx+bzy=0(2.4)
的特解作为新的自变量 ξ ξ ξ ,此时 A = 0 A=0 A=0,这样方程 (2.2) 得以化简为
B u η + C u + F = 0 Bu_η+Cu+F=0 Buη+Cu+F=0
从而像第一种类型一样积分求解。
上述一阶偏微分方程方程 (2.4) 可变为常微分方程求解,首先修改为
a ( − z x z y ) − b = 0 a(-\cfrac{z_x}{z_y})-b=0 a(zyzx)b=0
如果把 z ( x , y ) = const z(x,y)=\text{const} z(x,y)=const 作为定义隐函数 y = y ( x ) y=y(x) y=y(x) 的方程,则 d y / d x = − z x / z y dy/dx=-z_x/z_y dy/dx=zx/zy,从而
a ( x , y ) d y − b ( x , y ) d x = 0 (2.5) a(x,y)\mathrm dy-b(x,y)\mathrm dx=0\tag{2.5} a(x,y)dyb(x,y)dx=0(2.5)
常微分方程 (2.5) 叫做偏微分方程 (2.1) 的特征方程,特征方程的积分曲线叫做特征线。如果求出了积分曲线簇 z ( x , y ) = const z(x,y)=\text{const} z(x,y)=const ,再任取常数 c c c 使得 ∣ J ( ξ , η ) ∣ ≠ 0 |J(ξ,η)|\neq0 J(ξ,η)=0 ,以此 ξ , η ξ,η ξ,η 作为变量代换,则一阶线性偏微分方程便可求得通解。

二阶线性偏微分方程的分类和标准式

以下研究方程的分类,并把各类方程分别化为标准型,这样以后就只需讨论标准型的解法。

两个自变量方程的分类

先研究两自变量 x , y x,y x,y 的二阶线性偏微分方程
a 11 u x x + a 12 u x y + a 22 u y y + b 1 u x + b 2 u y + c u + f = 0 (3.1) a_{11}u_{xx}+a_{12}u_{xy}+a_{22}u_{yy}+b_1u_x+b_2u_y+cu+f=0\tag{3.1} a11uxx+a12uxy+a22uyy+b1ux+b2uy+cu+f=0(3.1)
其中 a 11 , a 12 , a 22 , b 1 , b 2 , c , f a_{11},a_{12},a_{22},b_1,b_2,c,f a11,a12,a22,b1,b2,c,f 都是自变量 x , y x,y x,y 的函数。我们希望通过适当的自变量变换和未知函数的变换,使方程简化,并在此基础上对方程进行分类。

作自变量非奇异变换(可逆)
{ ξ = ξ ( x , y ) η = η ( x , y ) \begin{cases} ξ=ξ(x,y) \\ η=η(x,y) \end{cases} {ξ=ξ(x,y)η=η(x,y)
其中雅可比行列式
∣ J ∣ = ∣ ∂ ( ξ , η ) ∂ ( x , y ) ∣ = ∣ ξ x ξ y η x η y ∣ ≠ 0 |J|=|\cfrac{∂(ξ,η)}{∂(x,y)}|=\begin{vmatrix}ξ_x&ξ_y \\ η_x&η_y\end{vmatrix}\neq0 J=(x,y)(ξ,η)=ξxηxξyηy=0
通过代换可得到新的方程
A 11 u ξ ξ + 2 A 12 u ξ η + A 22 u η η + B 1 u ξ + B 2 u η + C u + F = 0 (3.2) A_{11}u_{ξξ}+2A_{12}u_{ξη}+A_{22}u_{ηη}+B_1u_ξ+B_2u_η+Cu+F=0\tag{3.2} A11uξξ+2A12uξη+A22uηη+B1uξ+B2uη+Cu+F=0(3.2)
其中
{ A 11 = a 11 ξ x 2 + 2 a 12 ξ x ξ y + a 22 ξ y 2 A 12 = a 11 ξ x η x + a 12 ( ξ x η y + ξ y η x ) + a 22 ξ y η y A 22 = a 11 η x 2 + 2 a 12 η x η y + a 22 η y 2 B 1 = a 11 ξ x x + 2 a 12 ξ x y + a 22 ξ y y + b 1 ξ x + b 2 ξ y B 2 = a 11 η x x + 2 a 12 η x y + a 22 η y y + b 1 η x + b 2 η y C = c , F = f (3.3) \begin{cases} A_{11}=a_{11}ξ_x^2+2a_{12}ξ_xξ_y+a_{22}ξ_y^2 \\ A_{12}=a_{11}ξ_xη_x+a_{12}(ξ_xη_y+ξ_yη_x)+a_{22}ξ_yη_y \\ A_{22}=a_{11}η_x^2+2a_{12}η_xη_y+a_{22}η_y^2 \\ B_1=a_{11}ξ_{xx}+2a_{12}ξ_{xy}+a_{22}ξ_{yy}+b_1ξ_x+b_2ξ_y \\ B_2=a_{11}η_{xx}+2a_{12}η_{xy}+a_{22}η_{yy}+b_1η_x+b_2η_y \\ C=c,F=f \end{cases} \tag{3.3} A11=a11ξx2+2a12ξxξy+a22ξy2A12=a11ξxηx+a12(ξxηy+ξyηx)+a22ξyηyA22=a11ηx2+2a12ηxηy+a22ηy2B1=a11ξxx+2a12ξxy+a22ξyy+b1ξx+b2ξyB2=a11ηxx+2a12ηxy+a22ηyy+b1ηx+b2ηyC=c,F=f(3.3)
新方程仍然是线性的。
系数 A 11 , A 12 , A 22 A_{11},A_{12},A_{22} A11,A12,A22 可用矩阵表示为
D = ( A 11 A 12 A 21 A 22 ) = ( ξ x ξ y η x η y ) ( a 11 a 12 a 21 a 22 ) ( ξ x ξ y η x η y ) T (3.4) \mathbf D= \begin{pmatrix}A_{11}&A_{12} \\ A_{21}&A_{22}\end{pmatrix}= \begin{pmatrix}ξ_x&ξ_y \\ η_x&η_y\end{pmatrix} \begin{pmatrix}a_{11}&a_{12} \\ a_{21}&a_{22}\end{pmatrix} \begin{pmatrix}ξ_x&ξ_y \\ η_x&η_y\end{pmatrix}^{T} \tag{3.4} D=(A11A21A12A22)=(ξxηxξyηy)(a11a21a12a22)(ξxηxξyηy)T(3.4)
其中 A 12 = A 21 , a 12 = a 21 A_{12}=A_{21},a_{12}=a_{21} A12=A21,a12=a21

由 (3.4) 式看出,如果取一阶偏微分方程
a 11 z x 2 + 2 a 12 z x z y + a 22 z y 2 = 0 a_{11}z_x^2+2a_{12}z_xz_y+a_{22}z_y^2=0 a11zx2+2a12zxzy+a22zy2=0
两个线性无关的特解作为新的自变量 ξ , η ξ,η ξ,η ,此时 A 11 = A 22 = 0 A_{11}=A_{22}=0 A11=A22=0,这样方程 (3.2) 得以化简。
上述一阶偏微分方程可变为常微分方程求解,首先修改为
a 11 ( − z x z y ) 2 − 2 a 12 ( − z x z y ) + a 22 = 0 a_{11}(-\cfrac{z_x}{z_y})^2-2a_{12}(-\cfrac{z_x}{z_y})+a_{22}=0 a11(zyzx)22a12(zyzx)+a22=0
如果把 z ( x , y ) = const z(x,y)=\text{const} z(x,y)=const 作为定义隐函数 y = y ( x ) y=y(x) y=y(x) 的方程,则 d y / d x = − z x / z y dy/dx=-z_x/z_y dy/dx=zx/zy,从而
a 11 ( d y d x ) 2 − 2 a 12 ( d y d x ) + a 22 = 0 (3.5) a_{11}(\cfrac{dy}{dx})^2-2a_{12}(\cfrac{dy}{dx})+a_{22}=0\tag{3.5} a11(dxdy)22a12(dxdy)+a22=0(3.5)
常微分方程 (3.5) 叫做偏微分方程 (3.1) 的特征方程,特征方程的积分曲线叫做特征线
特征方程可分为两个方程
d y d x = a 12 ± a 12 2 − a 11 a 22 a 11 (3.6) \cfrac{dy}{dx}=\cfrac{a_{12} \pm \sqrt{a_{12}^2-a_{11}a_{22}}}{a_{11}}\tag{3.6} dxdy=a11a12±a122a11a22 (3.6)
通常根据判别式 Δ = a 12 2 − a 11 a 22 Δ=a_{12}^2-a_{11}a_{22} Δ=a122a11a22 划分方程 (3.1) 的类型
(1) Δ > 0 Δ>0 Δ>0 方程是双曲型;
(2) Δ = 0 Δ=0 Δ=0 方程是抛物型;
(3) Δ > 0 Δ>0 Δ>0 方程是椭圆型

由 (3.4) 矩阵的性质得到
A 12 2 − A 11 A 22 = det ⁡ ( D ) = ∣ J ∣ 2 ( a 12 2 − a 11 a 22 ) A_{12}^2-A_{11}A_{22}=\det(\mathbf D)=|J|^2(a_{12}^2-a_{11}a_{22}) A122A11A22=det(D)=J2(a122a11a22)
这就是说,作自变量变换时,方程的类型不变。

下面就方程中 Δ = a 12 2 − a 11 a 22 Δ=a_{12}^2-a_{11}a_{22} Δ=a122a11a22 符号的不同情况,讨论如何把方程 (3.1) 化为标准型。

  • 双曲型方程:方程 (3.6) 给出两族实的特征线 ξ ( x , y ) = c 1 , η ( x , y ) = c 2 ξ(x,y)=c_1,η(x,y)=c_2 ξ(x,y)=c1,η(x,y)=c2
    ξ = ξ ( x , y ) , η = η ( x , y ) ξ=ξ(x,y),η=η(x,y) ξ=ξ(x,y),η=η(x,y) 作为新的自变量,则 A 11 = A 22 = 0 A_{11}=A_{22}=0 A11=A22=0,从而自变量代换后的新方程 (3.2) 成为

u ξ η = − 1 2 A 12 ( B 1 u ξ + B 2 u η + C u + F ) u_{ξη}=-\cfrac{1}{2A_{12}}(B_1u_ξ+B_2u_η+Cu+F) uξη=2A121(B1uξ+B2uη+Cu+F)

或者,再做自变量变换 ξ = α + β , η = α − β ξ=α+β,η=α-β ξ=α+β,η=αβ ,则方程化为
u α α − u β β = − 1 A 12 [ ( B 1 + B 2 ) u α + ( B 1 − B 2 ) u β + 2 C u + 2 F ] u_{αα}-u_{ββ}=-\cfrac{1}{A_{12}}[(B_1+B_2)u_α+(B_1-B_2)u_β+2Cu+2F] uααuββ=A121[(B1+B2)uα+(B1B2)uβ+2Cu+2F]
上述两方程是双曲型方程的标准形式,弦振动方程即为标准型的双曲方程。

  • 抛物型方程:特征方程只能给出一族特征线 ξ ( x , y ) = c ξ(x,y)=c ξ(x,y)=c
    ξ = ξ ( x , y ) ξ=ξ(x,y) ξ=ξ(x,y) 作为新的自变量,任取一函数 η = η ( x , y ) η=η(x,y) η=η(x,y),与 ξ ξ ξ 线性无关。由于 a 12 2 − a 11 a 22 = 0 a_{12}^2-a_{11}a_{22}=0 a122a11a22=0,带入 (3.3) 可得
    A 11 = a 11 ξ x 2 + 2 a 12 ξ x ξ y + a 22 ξ y 2 = ( a 11 ξ x + a 22 ξ y ) 2 = 0 A 12 = a 11 ξ x η x + a 12 ( ξ x η y + ξ y η x ) + a 22 ξ y η y = ( a 11 ξ x + a 22 ξ y ) ( a 11 η x + a 22 η y ) = 0 A_{11}=a_{11}ξ_x^2+2a_{12}ξ_xξ_y+a_{22}ξ_y^2=(\sqrt{a_{11}}ξ_x+\sqrt{a_{22}}ξ_y)^2=0 \\ A_{12}=a_{11}ξ_xη_x+a_{12}(ξ_xη_y+ξ_yη_x)+a_{22}ξ_yη_y=(\sqrt{a_{11}}ξ_x+\sqrt{a_{22}}ξ_y)(\sqrt{a_{11}}η_x+\sqrt{a_{22}}η_y)=0 A11=a11ξx2+2a12ξxξy+a22ξy2=(a11 ξx+a22 ξy)2=0A12=a11ξxηx+a12(ξxηy+ξyηx)+a22ξyηy=(a11 ξx+a22 ξy)(a11 ηx+a22 ηy)=0
    从而自变量代换后的新方程成为

u η η = − 1 2 A 22 ( B 1 u ξ + B 2 u η + C u + F ) u_{ηη}=-\cfrac{1}{2A_{22}}(B_1u_ξ+B_2u_η+Cu+F) uηη=2A221(B1uξ+B2uη+Cu+F)

这是抛物型方程的标准形式,扩散方程、热传导方程都是标准形式的抛物型方程。

  • 椭圆形方程:特征方程给出两族复数的特征线 ξ ( x , y ) = c 1 , η ( x , y ) = c 2 ξ(x,y)=c_1,η(x,y)=c_2 ξ(x,y)=c1,η(x,y)=c2
    而且 ξ = ξ ( x , y ) , η = η ( x , y ) ξ=ξ(x,y),η=η(x,y) ξ=ξ(x,y),η=η(x,y) 共轭,对 ξ , η ξ,η ξ,η 再做自变量变换 ξ = α + i β , η = α − i β ξ=α+iβ,η=α-iβ ξ=α+iβ,η=αiβ ,则新方程化为

u α α + u β β = − 1 A 12 [ ( B 1 + B 2 ) u α + i ( B 2 − B 1 ) u β + 2 C u + 2 F ] u_{αα}+u_{ββ}=-\cfrac{1}{A_{12}}[(B_1+B_2)u_α+i(B_2-B_1)u_β+2Cu+2F] uαα+uββ=A121[(B1+B2)uα+i(B2B1)uβ+2Cu+2F]

这是椭圆型方程的标准形式,稳定温度分别、静电场方程都是椭圆型的双曲方程。

多自变量方程的分类

考虑多自变量 x 1 , x 2 , ⋯   , x n x_1,x_2,\cdots,x_n x1,x2,,xn 的二阶线性偏微分方程
∑ j = 1 n ∑ i = 1 n a i j u x i x j + ∑ i = 1 n b i u x i + c u + f = 0 (3.7) \displaystyle \sum_{j=1}^{n}\sum_{i=1}^n a_{ij}u_{x_ix_j}+\sum_{i=1}^{n}b_iu_{x_i}+cu+f=0 \tag{3.7} j=1ni=1naijuxixj+i=1nbiuxi+cu+f=0(3.7)
做自变量的非奇异代换
ξ k = ξ k ( x 1 , x 2 , ⋯   , x n ) , k = 1 , 2 , ⋯   , n ξ_k=ξ_k(x_1,x_2,\cdots,x_n),\quad k=1,2,\cdots,n ξk=ξk(x1,x2,,xn),k=1,2,,n
代换的雅克比行列式
∣ J ∣ = ∣ ∂ ( ξ 1 , ξ 2 , ⋯   , ξ n ) ∂ ( x 1 , x 2 , ⋯   , x n ) ∣ ≠ 0 |J|=|\cfrac{∂(ξ_1,ξ_2,\cdots,ξ_n)}{∂(x_1,x_2,\cdots,x_n)}|\neq0 J=(x1,x2,,xn)(ξ1,ξ2,,ξn)=0
通过代换可得到新的方程
∑ j = 1 n ∑ i = 1 n A i j u ξ i ξ j + ∑ i = 1 n B i u ξ i + C u + F = 0 (3.8) \displaystyle \sum_{j=1}^{n}\sum_{i=1}^n A_{ij}u_{ξ_iξ_j}+\sum_{i=1}^{n}B_iu_{ξ_i}+Cu+F=0 \tag{3.8} j=1ni=1nAijuξiξj+i=1nBiuξi+Cu+F=0(3.8)
注意到变换后的方程仍然是线性的,其中系数
{ A i j = ∑ l = 1 n ∑ k = 1 n a k l ∂ ξ i ∂ x k ∂ ξ j ∂ x l B i = ∑ k = 1 n ∂ ξ i ∂ x k + ∑ l = 1 n ∑ k = 1 n ∂ 2 ξ i ∂ x k ∂ x l \begin{cases} A_{ij}=\displaystyle \sum_{l=1}^{n}\sum_{k=1}^na_{kl}\cfrac{∂ξ_i}{∂x_k}\cfrac{∂ξ_j}{∂x_l} \\ \displaystyle B_i=\sum_{k=1}^n\cfrac{∂ξ_i}{∂x_k}+\sum_{l=1}^{n}\sum_{k=1}^n\cfrac{∂^2ξ_i}{∂x_k∂x_l} \end{cases} Aij=l=1nk=1naklxkξixlξjBi=k=1nxkξi+l=1nk=1nxkxl2ξi
系数 A i j A_{ij} Aij 可用矩阵表示为
( A i j ) n × n = J A J T (A_{ij})_{n\times n}=JAJ^T (Aij)n×n=JAJT
其中 A = ( a i j ) n × n A=(a_{ij})_{n\times n} A=(aij)n×n 是由方程 (3.7) 二阶偏导数的系数组成。值得注意的是系数变换公式恰恰是二次齐次式 ∑ j = 1 n ∑ i = 1 n a i j λ i λ j \displaystyle \sum_{j=1}^{n}\sum_{i=1}^n a_{ij}λ_iλ_j j=1ni=1naijλiλj 系数的非奇异线性变换,二次齐次式可以通过适当的变换而对角化,即
( A i j ) n × n = B A B T = d i a g ( A 11 , A 22 , ⋯   , A n n ) (A_{ij})_{n\times n}=BAB^T=\mathrm{diag}(A_{11},A_{22},\cdots,A_{nn}) (Aij)n×n=BABT=diag(A11,A22,,Ann)
其中 A i i ∈ { − 1 , 0 , 1 } A_{ii}\in\{-1,0,1\} Aii{1,0,1} ,根据二次齐次式对角化的惯性定律, A i i A_{ii} Aii − 1 , 0 , 1 -1,0,1 1,0,1 的个数是一定的,设 A i i A_{ii} Aii 中 1 的个数为 p p p ,-1 的个数为 q q q 。据此划分偏微分方程的类型:
(1) 当 p = n , q = 0 p=n,q=0 p=n,q=0,椭圆型(elliptic),标准形式是
∑ i = 1 n ∂ 2 u ∂ ξ i 2 + ∑ i = 1 n B i ∂ u ∂ ξ i + C u + F = 0 \displaystyle\sum_{i=1}^n\cfrac{∂^2u}{∂ξ_i^2}+\sum_{i=1}^{n}B_i\cfrac{∂u}{∂ξ_i}+Cu+F=0 i=1nξi22u+i=1nBiξiu+Cu+F=0
(2) 当 p = n − 1 , q = 0 p=n-1,q=0 p=n1,q=0,抛物型(parabolic),标准形式是
∑ i = 1 n − 1 ∂ 2 u ∂ ξ i 2 + ∑ i = 1 n B i ∂ u ∂ ξ i + C u + F = 0 \displaystyle\sum_{i=1}^{n-1}\cfrac{∂^2u}{∂ξ_i^2}+\sum_{i=1}^{n}B_i\cfrac{∂u}{∂ξ_i}+Cu+F=0 i=1n1ξi22u+i=1nBiξiu+Cu+F=0
(3) 当 p = 1 , q = n − 1 p=1,q=n-1 p=1,q=n1,双曲型(hyperbolic),标准形式是
∂ 2 u ∂ ξ 1 2 − ∑ j = 2 n ∂ 2 u ∂ ξ j 2 + ∑ i = 1 n B i ∂ u ∂ ξ i + C u + F = 0 \displaystyle\cfrac{∂^2u}{∂ξ_1^2}-\sum_{j=2}^n\cfrac{∂^2u}{∂ξ_j^2}+\sum_{i=1}^{n}B_i\cfrac{∂u}{∂ξ_i}+Cu+F=0 ξ122uj=2nξj22u+i=1nBiξiu+Cu+F=0
(4) 当 p > 0 , q > 0 p>0,q>0 p>0,q>0,超双曲型(ultrahyperbolic),标准形式是
∑ i = 1 p ∂ 2 u ∂ ξ i 2 − ∑ j = 1 p + q ∂ 2 u ∂ ξ i 2 + ∑ i = 1 n B i ∂ u ∂ ξ i + C u + F = 0 \displaystyle\sum_{i=1}^p\cfrac{∂^2u}{∂ξ_i^2}-\sum_{j=1}^{p+q}\cfrac{∂^2u}{∂ξ_i^2}+\sum_{i=1}^{n}B_i\cfrac{∂u}{∂ξ_i}+Cu+F=0 i=1pξi22uj=1p+qξi22u+i=1nBiξiu+Cu+F=0
常系数二阶线性偏微分方程:方程 (2.12) 中的系数 a i j , b i , c , f a_{ij},b_i,c,f aij,bi,c,f 都是常数,则方程还可以进一步化简。引入变换
u ( ξ , η ) = V ( ξ , η ) e λ ξ + μ η u(ξ,η)=V(ξ,η)e^{λξ+μη} u(ξ,η)=V(ξ,η)eλξ+μη
其中 λ , μ λ,μ λ,μ 是待定常数,根据方程的类型可选择 λ , μ λ,μ λ,μ,使一阶偏导数项或者函数项的系数为零。

波动方程的行波解

无限长弦振动初值问题:考虑两端为无限长的弦振动方程初值问题
{ u t t − a 2 u x x = 0 u ∣ t = 0 = ϕ ( x ) , ∂ u ∂ t ∣ t = 0 = ψ ( x ) (4.1) \begin{cases} u_{tt}-a^2u_{xx}=0 \\ u|_{t=0}=\phi(x),\cfrac{∂u}{∂t}|_{t=0}=\psi(x) \end{cases}\tag{4.1} utta2uxx=0ut=0=ϕ(x),tut=0=ψ(x)(4.1)
其中 ϕ ( x ) , ψ ( x ) \phi(x),\psi(x) ϕ(x),ψ(x) 分别表示初始位移和初始速度。
(1) 弦振动方程的通解:当 a > 0 a>0 a>0 为常数时,上述方程可分解为
( ∂ ∂ t + a ∂ u ∂ x ) ( ∂ ∂ t − a ∂ ∂ x ) u = 0 (\cfrac{∂}{∂t}+a\cfrac{∂u}{∂x})(\cfrac{∂}{∂t}-a\cfrac{∂}{∂x})u=0 (t+axu)(tax)u=0
等价于两个一阶线性偏微分方程组
{ ∂ u ∂ t + a ∂ u ∂ x = 0 ∂ u ∂ t − a ∂ u ∂ x = v \begin{cases} \cfrac{∂u}{∂t}+a\cfrac{∂u}{∂x}=0 \\ \cfrac{∂u}{∂t}-a\cfrac{∂u}{∂x}=v \end{cases} tu+axu=0tuaxu=v
这两个方程有一组独立的特征线 x + a t = c 1 , x − a t = c 2 x+at=c_1,x-at=c_2 x+at=c1,xat=c2,引入变换 ξ = x + a t , η = x − a t ξ=x+at,η=x-at ξ=x+at,η=xat
利用复合函数的求导法则求得 u t t , u x x u_{tt},u_{xx} utt,uxx 带入弦振动方程,得到新的方程
u ξ η = 0 u_{ξη}=0 uξη=0
上述方程容易求得通解
u = f 1 ( x + a t ) + f 2 ( x − a t ) (4.2) u=f_1(x+at)+f_2(x-at)\tag{4.2} u=f1(x+at)+f2(xat)(4.2)
其中 f 1 , f 2 f_1,f_2 f1,f2 是任意函数。

通解具有明确的物理意义,以 f 2 ( x − a t ) f_2(x-at) f2(xat) 而论,改用以速度 a a a 沿 x x x 轴正方向移动的坐标轴 X X X,则新旧坐标和时间之间的关系为 X = x − a t , T = t X=x-at,T=t X=xat,T=t,而 f 2 ( x − a t ) = f ( X ) f_2(x-at)=f(X) f2(xat)=f(X) 与时间 T T T 无关。这就是说函数图像在动坐标系中保持不变,即是随着动坐标系以速度 a a a 沿 x x x 轴正方向移动的行波。同理, f 1 ( x + a t ) f_1(x+at) f1(x+at) 是以速度 a a a 沿 x x x 轴负方向移动的行波。这样,弦振动方程描述的是以速度 a a a 沿 x x x 轴正负两方向移动的行波。

(2) 函数 f 1 , f 2 f_1,f_2 f1,f2 由初始条件确定
将通解带入初始条件可以得到 { f 1 ( x ) + f 2 ( x ) = ϕ ( x ) a f 1 ′ ( x ) − a f 2 ′ ( x ) = ψ ( x ) \begin{cases} f_1(x)+f_2(x)=\phi(x) \\af'_1(x)-af'_2(x)=\psi(x) \end{cases} {f1(x)+f2(x)=ϕ(x)af1(x)af2(x)=ψ(x) ,由此可解得
f 1 ( x ) = 1 2 ϕ ( x ) + 1 2 a ∫ x 0 x ψ ( ξ ) d ξ + 1 2 [ f 1 ( x 0 ) − f 2 ( x 0 ) ] f 2 ( x ) = 1 2 ϕ ( x ) − 1 2 a ∫ x 0 x ψ ( ξ ) d ξ − 1 2 [ f 1 ( x 0 ) − f 2 ( x 0 ) ] f_1(x)=\frac{1}{2}\phi(x)+\frac{1}{2a}\int_{x_0}^x\psi(ξ)dξ+\frac{1}{2}[f_1(x_0)-f_2(x_0)] \\ f_2(x)=\frac{1}{2}\phi(x)-\frac{1}{2a}\int_{x_0}^x\psi(ξ)dξ-\frac{1}{2}[f_1(x_0)-f_2(x_0)] f1(x)=21ϕ(x)+2a1x0xψ(ξ)dξ+21[f1(x0)f2(x0)]f2(x)=21ϕ(x)2a1x0xψ(ξ)dξ21[f1(x0)f2(x0)]
将此代回 (4.2) 可得,初始问题的解
u ( x , t ) = 1 2 [ ϕ ( x + a t ) + ϕ ( x − a t ) ] + 1 2 a ∫ x − a t x + a t ψ ( ξ ) d ξ (4.3) u(x,t)=\frac{1}{2}[\phi(x+at)+\phi(x-at)]+\frac{1}{2a}\int_{x-at}^{x+at}\psi(ξ)dξ\tag{4.3} u(x,t)=21[ϕ(x+at)+ϕ(xat)]+2a1xatx+atψ(ξ)dξ(4.3)
这个公式称为达朗贝尔(d’Alembert)公式

(3) 物理意义示例:设初速度为零 ψ ( x ) = 0 \psi(x)=0 ψ(x)=0,初始位移 ϕ ( x ) \phi(x) ϕ(x)
如下图 t = 0 t=0 t=0 时刻函数图像实线所示,达朗贝尔公式给出 u ( x , t ) = 1 2 ϕ ( x + a t ) + 1 2 ϕ ( x − a t ) u(x,t)=\frac{1}{2}\phi(x+at)+\frac{1}{2}\phi(x-at) u(x,t)=21ϕ(x+at)+21ϕ(xat) ,即初始位移分为两半(下图虚线),分别以速度 a a a 向左右两端移动,这两个行波的和给出各个时刻的波形(由下向上的实线)。

偏微分方程(Partial Differential Equation I)_第1张图片

半无限弦振动定解问题 ——延拓法
(1) 端点固定弦的自由振动:定解问题为
{ u t t − a 2 u x x = 0 ( x > 0 ) u ∣ t = 0 = ϕ ( x ) , ∂ u ∂ t ∣ t = 0 = ψ ( x ) ( x ⩾ 0 ) u ∣ x = 0 = 0 (4.4) \begin{cases} u_{tt}-a^2u_{xx}=0 \quad (x>0)\\ u|_{t=0}=\phi(x),\cfrac{∂u}{∂t}|_{t=0}=\psi(x) \quad(x\geqslant 0)\\ u|_{x=0}=0 \end{cases}\tag{4.4} utta2uxx=0(x>0)ut=0=ϕ(x),tut=0=ψ(x)(x0)ux=0=0(4.4)
求解上述问题的基本思路是把半无限长当做无限长弦的 x ⩾ 0 x\geqslant0 x0 的部分,然后用达朗贝尔公式求解即可。由定解条件知,振动过程 x = 0 x=0 x=0 处必须保持不动,由微积分知识可知,如果一个连续可微函数 g ( x ) g(x) g(x) ( − ∞ , + ∞ ) (-\infty,+\infty) (,+) 上是奇函数,则必有 g ( 0 ) = 0 g(0)=0 g(0)=0 。因此对函数 ϕ , ψ \phi,\psi ϕ,ψ 进行奇延拓,即
Φ ( x ) = { ϕ ( x ) ( x ⩾ 0 ) − ϕ ( x ) ( x < 0 ) ; Ψ ( x ) = { ψ ( x ) ( x ⩾ 0 ) − ψ ( x ) ( x < 0 ) \Phi(x)=\begin{cases}\phi(x)&(x\geqslant0)\\ -\phi(x)&(x<0)\end{cases};\quad \Psi(x)=\begin{cases}\psi(x)&(x\geqslant0)\\ -\psi(x)&(x<0)\end{cases} Φ(x)={ϕ(x)ϕ(x)(x0)(x<0);Ψ(x)={ψ(x)ψ(x)(x0)(x<0)
将上式带入达朗贝尔公式,可求得
u ( x , t ) = { 1 2 [ ϕ ( x + a t ) + ϕ ( x − a t ) ] + 1 2 a ∫ x − a t x + a t ψ ( ξ ) d ξ ( t ⩽ x a ) 1 2 [ ϕ ( x + a t ) − ϕ ( x − a t ) ] + 1 2 a ∫ a t − x x + a t ψ ( ξ ) d ξ ( t > x a ) (4.5) u(x,t)=\begin{cases}\displaystyle \frac{1}{2}[\phi(x+at)+\phi(x-at)]+\frac{1}{2a}\int_{x-at}^{x+at}\psi(ξ)dξ &(t⩽\cfrac{x}{a})\\ \displaystyle \frac{1}{2}[\phi(x+at)-\phi(x-at)]+\frac{1}{2a}\int_{at-x}^{x+at}\psi(ξ)dξ &(t>\cfrac{x}{a}) \end{cases}\tag{4.5} u(x,t)=21[ϕ(x+at)+ϕ(xat)]+2a1xatx+atψ(ξ)dξ21[ϕ(x+at)ϕ(xat)]+2a1atxx+atψ(ξ)dξ(tax)(t>ax)(4.5)
物理意义示例:设初速度为零 ψ ( x ) = 0 \psi(x)=0 ψ(x)=0,初始位移 ϕ ( x ) \phi(x) ϕ(x)
图中由下向上的实线描述了波形变化,右边的波分别向左右两端移动,左端奇延拓的波也分别向左右两端移动。端点固定的影响表现为反射波,反射波的相位和入射波相反,这就是半波损失
偏微分方程(Partial Differential Equation I)_第2张图片

(2) 端点自由弦的自由振动:定解问题为
{ u t t − a 2 u x x = 0 u ∣ t = 0 = ϕ ( x ) , ∂ u ∂ t ∣ t = 0 = ψ ( x ) ∂ u ∂ x ∣ x = 0 = 0 ( 0 ⩽ x < ∞ ) (4.6) \begin{cases} u_{tt}-a^2u_{xx}=0 \\ u|_{t=0}=\phi(x),\cfrac{∂u}{∂t}|_{t=0}=\psi(x) \\ \cfrac{∂u}{∂x}|_{x=0}=0 \end{cases}\quad(0\leqslant x<\infty)\tag{4.6} utta2uxx=0ut=0=ϕ(x),tut=0=ψ(x)xux=0=0(0x<)(4.6)
同样对函数进行偶延拓
Φ ( x ) = { ϕ ( x ) ( x ⩾ 0 ) ϕ ( − x ) ( x < 0 ) ; Ψ ( x ) = { ψ ( x ) ( x ⩾ 0 ) ψ ( − x ) ( x < 0 ) \Phi(x)=\begin{cases}\phi(x)&(x\geqslant0)\\ \phi(-x)&(x<0)\end{cases};\quad \Psi(x)=\begin{cases}\psi(x)&(x\geqslant0)\\ \psi(-x)&(x<0)\end{cases} Φ(x)={ϕ(x)ϕ(x)(x0)(x<0);Ψ(x)={ψ(x)ψ(x)(x0)(x<0)
将上式带入达朗贝尔公式,可求得
u ( x , t ) = { 1 2 [ ϕ ( x + a t ) + ϕ ( x − a t ) ] + 1 2 a ∫ x − a t x + a t ψ ( ξ ) d ξ ( t ⩽ x a ) 1 2 [ ϕ ( x + a t ) + ϕ ( a t − x ) ] + 1 2 a [ ∫ 0 x + a t ψ ( ξ ) d ξ + ∫ 0 a t − x ψ ( ξ ) d ξ ] ( t > x a ) (4.7) u(x,t)=\begin{cases}\displaystyle \frac{1}{2}[\phi(x+at)+\phi(x-at)]+\frac{1}{2a}\int_{x-at}^{x+at}\psi(ξ)dξ &(t⩽\cfrac{x}{a})\\ \displaystyle \frac{1}{2}[\phi(x+at)+\phi(at-x)]+\frac{1}{2a}[\int_{0}^{x+at}\psi(ξ)dξ+\int_{0}^{at-x}\psi(ξ)dξ] &(t>\cfrac{x}{a}) \end{cases}\tag{4.7} u(x,t)=21[ϕ(x+at)+ϕ(xat)]+2a1xatx+atψ(ξ)dξ21[ϕ(x+at)+ϕ(atx)]+2a1[0x+atψ(ξ)dξ+0atxψ(ξ)dξ](tax)(t>ax)(4.7)
自由端点的影响也是一种反射波,不同的是反射波的相位和入射波相同,没有半波损失。

三维波动方程的球面波:三维波动方程在求坐标系 ( r , θ , ϕ ) (r,θ,ϕ) (r,θ,ϕ) 下可以表示为
1 a 2 ∂ 2 u ∂ t 2 = 1 r 2 ∂ ∂ r ( r 2 ∂ u ∂ r ) + 1 r 2 sin ⁡ θ ∂ ∂ θ ( sin ⁡ θ ∂ u ∂ θ ) + 1 r 2 sin ⁡ 2 θ ∂ 2 u ∂ ϕ 2 (4.8) \cfrac{1}{a^2}\cfrac{∂^2u}{∂t^2}=\cfrac{1}{r^2}\cfrac{∂}{∂r}\left(r^2\cfrac{∂u}{∂r}\right) +\cfrac{1}{r^2\sinθ}\cfrac{∂}{∂θ}\left(\sinθ\cfrac{∂u}{∂θ}\right) +\cfrac{1}{r^2\sin^2θ}\cfrac{∂^2u}{∂ϕ^2}\tag{4.8} a21t22u=r21r(r2ru)+r2sinθ1θ(sinθθu)+r2sin2θ1ϕ22u(4.8)
若函数 u u u 是球对称函数,指 u u u θ , ϕ θ,ϕ θ,ϕ 无关,即 u = u ( r , t ) u=u(r,t) u=u(r,t),方程可化为
∂ 2 u ∂ t 2 = a 2 ( ∂ 2 u ∂ r 2 + 2 r ∂ u ∂ r ) \cfrac{∂^2u}{∂t^2}=a^2(\cfrac{∂^2u}{∂r^2}+\cfrac{2}{r}\cfrac{∂u}{∂r}) t22u=a2(r22u+r2ru)
进一步重写为
∂ 2 ( r u ) ∂ t 2 = a 2 ∂ 2 ( r u ) ∂ r 2 \cfrac{∂^2(ru)}{∂t^2}=a^2\cfrac{∂^2(ru)}{∂r^2} t22(ru)=a2r22(ru)
这时关于 r u ru ru 的一维波动方程,其通解为
r u ( r , t ) = f 1 ( r + a t ) + f 2 ( r − a t ) ru(r,t)=f_1(r+at)+f_2(r-at) ru(r,t)=f1(r+at)+f2(rat)
从而
u = f 1 ( r + a t ) + f 2 ( r − a t ) r (4.9) u=\cfrac{f_1(r+at)+f_2(r-at)}{r}\tag{4.9} u=rf1(r+at)+f2(rat)(4.9)
其中 f 1 , f 2 f_1,f_2 f1,f2 是任意两个二阶连续可微函数,可见,在球对称下,三维波的传播是以球心为中心,沿半径 r r r 传播的球面波。

线性叠加原理

线性叠加原理:考虑 n n n 个自变量(包括时间 t t t x 1 , x 2 , ⋯   , x n x_1,x_2,\cdots,x_n x1,x2,,xn 的二阶线性偏微分方程
∑ j = 1 n ∑ i = 1 n a i j u x i x j + ∑ i = 1 n b i u x i + c u = f \displaystyle\sum_{j=1}^{n}\sum_{i=1}^n a_{ij}u_{x_ix_j}+\sum_{i=1}^{n}b_iu_{x_i}+cu=f j=1ni=1naijuxixj+i=1nbiuxi+cu=f
其中 a i j , b i , c , f a_{ij},b_i,c,f aij,bi,c,f x 1 , x 2 , ⋯   , x n x_1,x_2,\cdots,x_n x1,x2,,xn 的已知函数, f f f 称为方程的非齐次项。

引入微分算子
L = ∑ j = 1 n ∑ i = 1 n a i j ∂ ∂ x i ∂ x j + ∑ i = 1 n b i ∂ ∂ x i + c L=\sum_{j=1}^{n}\sum_{i=1}^n a_{ij}\cfrac{∂}{∂x_i∂x_j} +\sum_{i=1}^{n}b_i\cfrac{∂}{∂x_i}+c L=j=1ni=1naijxixj+i=1nbixi+c
则可简单的表示为
L [ u ] = f L[u]=f L[u]=f
一般的,对任意常数 c 1 , c 2 c_1,c_2 c1,c2 和任意函数 u 1 , u 2 u_1,u_2 u1,u2 ,微分算子满足性质
L [ c 1 u 1 + c 2 u 2 ] = c 1 L [ u 1 ] + c 2 L [ u 2 ] (5.1) L[c_1u_1+c_2u_2]=c_1L[u_1]+c_2L[u_2]\tag{5.1} L[c1u1+c2u2]=c1L[u1]+c2L[u2](5.1)
称为线性微分算子。显然, L L L 为线性微分算子。
对于一般的线性边界条件(包括三类边界条件)也可以写成算子的形式
L 0 [ u ] = ( α ∂ u ∂ n + β u ) ∣ ∂ Ω = ϕ L_0[u]=(\alpha\cfrac{∂u}{∂n}+\beta u)|_{∂\Omega}=\phi L0[u]=(αnu+βu)Ω=ϕ
其中 α , β , ϕ \alpha,\beta,\phi α,β,ϕ 是已知函数,容易证明 L 0 L_0 L0 也是线性微分算子。

叠加原理
(1) 有限叠加原理:若 u i u_i ui 分别满足方程 L [ u i ] = f i ( i = 1 , 2 , ⋯   , m ) L[u_i]=f_i\quad(i=1,2,\cdots,m) L[ui]=fi(i=1,2,,m) ,则他们的线性组合 u = ∑ i = 1 m c i u i u=\displaystyle\sum_{i=1}^mc_iu_i u=i=1mciui 满足方程
L [ u ] = ∑ i = 1 m c i f i L[u]=\displaystyle\sum_{i=1}^mc_if_i L[u]=i=1mcifi
(2) 级数叠加原理:若 u i u_i ui 分别满足方程 L [ u i ] = f i ( i = 1 , 2 , ⋯   ) L[u_i]=f_i\quad(i=1,2,\cdots) L[u

你可能感兴趣的:(数学(Math))