\[\newcommand{\half}[1]{\dfrac{#1}{2}} \newcommand{\PD}[2]{\dfrac{\partial #1}{\partial #2}} \newcommand{\pd}[2]{\frac{\partial #1}{\partial #2}} \newcommand{\eps}{\varepsilon} \newcommand{\norm}[1]{\left\Vert#1\right\Vert} \newcommand{\set}[1]{\left\{#1\right\}} \newcommand{\BX}{\mathbf{B}(X)} \newcommand{\Real}{\mathbb R} \newcommand{\A}{\mathcal{A}} \newcommand{\To}{\longrightarrow} \newcommand{\vecn}{\vec{\mathbf{n}}} \newcommand{\Vector}[1]{\left[\begin{array}{c} #1 \end{array}\right]} \newcommand{\Reci}[1]{\dfrac{1}{#1}} \newcommand{\reci}[1]{\frac{1}{#1}} \newcommand{\dx}{\Delta x} \newcommand{\dy}{\Delta y} \newcommand{\DS}{\displaystyle} \newcommand{\limz}[1]{\lim_{#1\to0}} \newcommand{\dlimz}[1]{\displaystyle\lim_{#1\to0}} \newcommand{\dint}{\displaystyle\int} \newcommand{\dlim}{\displaystyle\lim} \newcommand{\casefunc}[1]{\left\{\begin{aligned} #1 \end{aligned}\right.} \newcommand{\det}[1]{\left|\begin{aligned} #1 \end{aligned}\right|} \]

分离变量法

3. 固有值问题的Sturm-Liouville理论

张瑞
中国科学技术大学数学科学学院

固有值问题的Sturm-Liouville理论

一般的二阶齐次线性方程的分离变量

\[L_tu+c(t)L_xu=0, x\in[a,b], t\in I \] 这里$I$是有限或无限区间,$L_t$$L_x$是二阶线性偏微分算子

\[L_t=a_0(t)\dfrac{\partial^2}{\partial t^2}+a_1(t)\dfrac{\partial }{\partial t}+a_2(t) \]

\[L_x=b_0(x)\dfrac{\partial^2}{\partial x^2}+b_1(x)\dfrac{\partial }{\partial x}+b_2(x) \] 其中$a_i(t), b_i(x), i=0,1,2$$c(t)$是满足一定条件的已知函数。

分离变量,设$u=X(x)T(t)$,代入方程,可得

\[\dfrac{L_t T}{cT}+\dfrac{L_x X}{X}=0 \] 这样,有2个常微分方程

\[L_x X(x)+\lambda X(x)=0 \]\[L_t T(t)-\lambda c(t) T(t)=0 \]

$X(x)$的常微分方程和边界条件,构成了固有值问题。

$y(x)=X(x)$,则固有值问题的方程可以写成

\[b_0(x)y''(x)+b_x(x)y'(x)+b_2(x)y(x)+\lambda y(x)=0 \]$\rho(x)$,满足

\[(\rho(x)b_0(x))'=\rho(x)b_1(x) \] 则有

\[\rho(x)=\dfrac1{b_0(x)}e^{\int\frac{b_1(x)}{b_0(x)}dx} \]

$k(x)=\rho(x)b_0(x)$$-q(x)=\rho(x)b_2(x)$,则方程可以写成 Sturm-Liouville型方程

\[\dfrac{d}{dx}\left(k(x)\dfrac{dy}{dx}\right)-q(x)y+\lambda \rho(x)y=0 \]

对系数的要求

对 Sturm-Liouville型方程

\[\dfrac{d}{dx}\left(k(x)\dfrac{dy}{dx}\right)-q(x)y+\lambda \rho(x)y=0 \] 的系数有如下假定:

  • $[a,b]$上,$k(x), k'(x), \rho(x)$连续;当$x\in(a,b)$时,$k(x)>0, \rho(x)>0, q(x)\geq0$,而$a,b$至多是$k(x), \rho(x)$的1级零点
  • $q(x)$$(a,b)$上连续,而在端点处至多有$1$级极点

常见的5种边界条件

  • $k(a)>0$,且$q(x)$$a$点连续,对$u$可以给出第一、二、三类边界条件,即 \[(\alpha_1 u-\beta_1 u)|_{x=a}=0 , \alpha_1, \beta_1\geq0, \alpha_1^2+\beta_1^2\neq 0 \] 此时,得到的固有值问题的边界条件为 \[\alpha_1 y'(a)-\beta_1y(a)=0 \]

  • $k(b)>0$,且$q(x)$$b$点连续,对$u$可以给出第一、二、三类边界条件,即 \[(\alpha_2 u+\beta_2 u)|_{x=b}=0 , \alpha_2, \beta_2\geq0, \alpha_2^2+\beta_2^2\neq 0 \] 此时,得到的固有值问题的边界条件为 \[\alpha_2 y'(b)+\beta_2y(b)=0 \]

  • $k(a)=k(b)>0$时,可以加上$u$周期的。可以得到$y(x)$的条件为 \[y(a)=y(b), y'(a)=y'(b) \]$k(x), \rho(x)$$[a,b]$上为正时,称方程为正则的;否则,称为奇异的

  • $k(a)=0$时,要附加有界性条件。此时$y(a)$有界,即 \[|y(a)|<+\infty \] 称这种边界条件为自然边界条件

  • $k(b)=0$时,要附加$b$处的自然边界条件,即 \[|y(b)|<+\infty \]

例 1. $k(x)=1$, $q(x)=0$, $\rho(x)=1$,则方程为

\[y''+\lambda y=0 \]

$a=0$, $b=l$,则

$k(a)=1>0$,则$a$处可以取第一、二、三类边界条件,如

\[y(0)=0 \]

$k(l)=1>0$,则$b$处处可以取第一、二、三类边界条件,如

\[y(l)=0 \]

例 2. $k(x)=1$, $q(x)=0$, $\rho(x)=1$,则方程为

\[y''+\lambda y=0 \]

$a=0$, $b=l$,则

$k(a)=k(b)=1>0$,则可以加上周期边界条件:

\[y(x)=y(x+l) \]

例 3. $k(x)=1-x^2$, $q(x)=0$, $\rho(x)=1$,则方程为

\[((1-x^2)y')'+\lambda y=0 \]

$a=0$, $b=1$,则

$k(a)=k(b)=0$,则可以加上自然边界条件:

\[|y(\pm 1)|<\infty \]

例 4. $k(x)=x$, $q(x)=\frac{m^2}x$, $\rho(x)=x$,则方程为

\[(x y')'-\frac{m^2}x y+\lambda xy=0 \]

$a=0$, $b=1$,则

$k(a)=0$,则可以加上自然边界条件:

\[|y(0)|<\infty \]

$k(b)>0$,则可以加上第一、二、三类边界条件,如

\[y(1)=0 \]

Sturm-Liouville定理

$k(x),\rho(x),q(x)$满足前面的条件,则固有值问题 \[\casefunc{ & \dfrac{d}{dx}\left(k(x)\dfrac{dy}{dx}\right)-q(x)y+\lambda \rho(x)y=0 \\ & \mbox{5种边界条件中的任一种} } \] 的固有值和固有函数有如下性质:

  1. 可数性。存在可数无穷个固有值$\lambda_1<\lambda_2<\cdots$,且 \[\lim_{n\to+\infty}\lambda_n=+\infty \] 与第一个固有值对应的线性无关的固有函数有且只有一个(对周期边界条件不成立)
  2. 非负性。$\lambda_n\geq0$;且固有值$\lambda=0$的充要条件为:$q(x)\equiv0$,有固有值问题中$a,b$两端都不取第一、三类边界条件,这时,相应 的固有函数为常数;

非负性的证明

\[\int_a^b(\dfrac{d}{dx}\left(k(x)\dfrac{dy}{dx}\right)-q(x)y+\lambda \rho(x)y)ydx=0 \] 则有

\[\lambda\int_a^b\rho(x)y^2dx=\int_a^bq(x)y^2dx-\int_a^by(k(x)y')'dx \\ =\int_a^b q(x)y^2dx-k(x)y'y|_a^b+\int_a^b k(x)(y')^2dx \\ =\int_a^b q(x)y^2dx+\int_a^b k(x)(y')^2dx \\ +k(a)y(a)y'(a)-k(b)y(b)y'(b) \]

若能证明

\[k(a)y(a)y'(a)-k(b)y(b)y'(b)\geq0 \]$\lambda>0$。上式跟边界条件有关。


(1) $k(a)>0$$a$端有第一、二、三类边界条件

\[\alpha_1y'(a)-\beta_1y(a)=0, \alpha_1,\beta_1\geq0,\alpha_1^2+\beta_1^2\neq 0 \] 不防设$\alpha_1\neq0$,则

\[y'(a)=\frac{\beta_1}{\alpha_1}y(a), \frac{\beta_1}{\alpha_1}\geq0 \] 即,$y'(a)$$y(a)$同号。这样有

\[k(a)y(a)y'(a)\geq0 \]

(2) $k(b)>0$$b$端有第一、二、三类边界条件

\[\alpha_1y'(a)+\beta_1y(a)=0, \alpha_1,\beta_1\geq0,\alpha_1^2+\beta_1^2\neq 0 \] 不防设$\alpha_1\neq0$,则

\[y'(b)=-\frac{\beta_1}{\alpha_1}y(a), \frac{\beta_1}{\alpha_1}\geq0 \] 即,$y'(a)$$y(a)$异号。这样有

\[-k(b)y(b)y'(b)\geq0 \]

(3) $a$端有自然边界,即$|y(a)|<+\infty$,且$k(a)=0$,则

\[k(a)y(a)y'(a)=0 \]

(4) $b$端有自然边界,即$|y(b)|<+\infty$,且$k(b)=0$,则

\[k(b)y(b)y'(b)=0 \]

(5) 周期边界条件。此时$y(a)=y(b)$, $y'(a)=y'(b)$, $k(a)=k(b)$,即

\[k(a)y(a)y'(a)-k(b)y(b)y'(b)=0 \]


下证$\lambda=0$时的充要条件。

  1. 正交性。$\lambda_m, \lambda_n$为任意两个不同的固有值,则相应的固有函数满足 \[\int_a^b\rho(x)y_m(x)y_n(x)dx=0 \]
  2. 完备性。固有函数系$\{y_n(x)\}$是完备的。即,对任意有一阶连续导数及分段二阶连续导数的函数$f(x)$,只要它满足固有值问题中的边界条件,则它可以按固有函数系$\{y_n(x)\}$展开成绝对且一致收敛的广义Fourier级数 \[f(x)=\sum_{n=1}^{+\infty}f_ny_n(x) \]

正交性的证明:

\[\dfrac{d}{dx}\left(k(x)\dfrac{dy_m}{dx}\right)-q(x)y_m+\lambda_m \rho(x)y_m=0 \\ \dfrac{d}{dx}\left(k(x)\dfrac{dy_n}{dx}\right)-q(x)y_n+\lambda_n \rho(x)y_n=0 \]

第一式乘$y_n$与第二式乘$y_m$相减后,得到

\[(k(x)y'_m)'y'_n-(k(x)y_n)y'_m+(\lambda_m-\lambda_n)\rho(x)y_my_n=0 \]$[a,b]$上积分,有

\[(\lambda_m-\lambda_n)\int_a^b\rho(x)y_m(x)y_n(x)dx \\ =\int_a^b(-(k(x)y'_m)'y'_n+(k(x)y_n)y'_m)dx \]

右端使用分部积分,有

\[=k(a)(y_n(a)y'_m(a)-y_m(a)y'_n(a))\\ -k(b)(y_n(b)y'_m(b)-y_m(b)y'_n(b)) \]

就5种边界条件,讨论可得,值为$0$。即有

\[\int_a^b\rho(x)y_m(x)y_n(x)dx=0 \]

其中 \[f_k=\dfrac{1}{\|y_k(x)\|^2}\int_a^b\rho(x)f(x)y_k(x)dx \]

这里,

\[\|y_k(x)\|^2=\int_a^b\rho(x)y^2_k(x)dx \]

例 5. 解固有值问题

\[\casefunc{ & y''+\lambda y=0 , x\in(-l,l) \\ & y'(-l)=y'(l)=0 } \]

例 6. 解固有值问题

\[\casefunc{ & x^2y''+ xy'(x)+\lambda y=0 , x\in(1,e) \\ & y(1)=y(e)=0 } \]

例 7. 长为$l$的导热细杆,侧面绝热,且杆的左端绝热,右端与周围保持常温为零的介质的热交换,杆身没有热源,杆的初始温度为$\phi(x)$,求杆内的温度分布。

vertical slide1

vertical slide 2

例 8. 求定解问题

\[\casefunc{ & \triangle_3u=0 , x\in(0,a),y\in(0,b),z\in(0,c) \\ & u(0,y,z)=u(a,y,z)=u(x,0,z)=u(x,b,z)=0 \\ & u(x,y,0)=0, u(x,y,c)=\phi(x,y) } \]

9.

Step 1

分离变量

第二步

解固有值问题

二阶方程的固有值问题的解

对于二阶方程的固有值问题

\[\casefunc{ & u''+\lambda u=0 , x\in(0,l) \\ & u \mbox{在} x=0,l \mbox{处的边界条件} } \]

总结如下的固有值和固有函数:

(1) $u(0)=u(l)=0$时, \[\lambda_n=(\frac{n\pi}l)^2 , u_n(x)=\sin\frac{n\pi x}l \]

(2) $u'(0)=u'(l)=0$时,$0$也是特征值 \[\lambda_n=(\frac{n\pi}{l})^2 , u_n(x)=\cos\frac{n\pi x}{l} \]

(3) $u(0)=u'(l)=0$时, \[\lambda_n=(\frac{(2n-1)\pi}{2l})^2 , u_n(x)=\sin\frac{(2n-1)\pi x}{2l} \]

(4) $u'(0)=u(l)=0$时, \[\lambda_n=(\frac{(2n-1)\pi}{2l})^2 , u_n(x)=\cos\frac{(2n-1)\pi x}{2l} \]

(5) $u(0)=u(l)$, $u'(0)=u'(l)$时(周期边界), $0$也是固有值 \[\lambda_n=(\frac{2n\pi}{l})^2 , u_n(x)=\{\sin(\frac{2n\pi x}l),\cos\frac{2n\pi x}{l}\} \]

(6) $u(0)=u'(l)+\sigma u(l)=0$时, \[\lambda_n=(\frac{\gamma_n}{l})^2 , u_n(x)=\sin\frac{\gamma_n x}l \] 其中$\gamma_n$是方程$\tan\gamma=-\frac{\gamma}{\sigma l}$的第$n$个正根,$\sigma>0$

(7) $u'(0)=u'(l)+\sigma u(l)=0$时, \[\lambda_n=(\frac{\gamma_n}{l})^2 , u_n(x)=\cos\frac{\gamma_n x}l \] 其中$\gamma_n$是方程$\cot\gamma=\frac{\gamma}{\sigma l}$的第$n$个正根,$\sigma>0$

(8) $u'(0)-\sigma u(0)=u(l)=0$时, \[\lambda_n=(\frac{\gamma_n}{l})^2 , u_n(x)=\frac{\gamma_n}{\sigma l}\cos\frac{\gamma x}l+\sin\frac{\gamma_n x}l \] 其中$\gamma_n$是方程$\tan\gamma=-\frac{\gamma}{\sigma l}$的第$n$个正根,$\sigma>0$

(9) $u'(0)-\sigma u(0)=u(l)=0$时, \[\lambda_n=(\frac{\gamma_n}{l})^2 , u_n(x)=\frac{\gamma_n}{\sigma l}\cos\frac{\gamma x}l+\sin\frac{\gamma_n x}l \] 其中$\gamma_n$是方程$\cot\gamma=\frac{\gamma}{\sigma l}$的第$n$个正根,$\sigma>0$

(10) $u'(0)-\sigma_1 u(0)=u'(l)+\sigma_2 u(l)=0$时, \[\lambda_n=(\frac{\gamma_n}{l})^2 , u_n(x)=? \]

目录

谢谢

例 9. 本节读完

9.

注解

注解. 什么是分离变量