查看“︁特征线法”︁的源代码

[[数学]]中的'''特征线法'''是求解[[偏微分方程]]的一种方法，适用于准线性偏微分方程的求解。只要初始值不是沿着特征线给定，即可通过特征线法获得偏微分方程的精确解。 其基本思想是通过把双曲线型的准线性偏微分方程转化为两组常微分方程，再对常微分方程进行求解。两组常微分方程中的一组用于定义特征线，另一组用以描述解沿给定特征线变化。

== 基本方法 ==
设所需求解的准线性偏微分方程为
{{NumBlk|:|<math>a_1 \left( \boldsymbol{x}, u \right) \frac{\partial u}{\partial x_1} + a_2 \left( \boldsymbol{x}, u \right) \frac{\partial u}{\partial x_2} + \cdots + a_N \left( \boldsymbol{x}, u \right) \frac{\partial u}{\partial x_N} = b \left( \boldsymbol{x}, u \right)</math>|{{EquationRef|1}}}}

其中 <math>u \left( \boldsymbol{x} \right) = u \left( x_{1,} x_2, \ldots, x_N \right)</math>。

取某变量 <math>s</math>，令 <math>u(\boldsymbol{x})</math> 对 <math>s</math> 求导數，可得
{{NumBlk|:|<math>\frac{\text{d}u}{\text{d}s} = \left( \frac{\partial x_1}{\partial s} \right) \frac{\partial u}{\partial x_1} + \left( \frac{\partial x_2}{\partial s} \right) \frac{\partial u}{\partial x_2} + \ldots + \left( \frac{\partial x_N}{\partial s} \right) \frac{\partial u}{\partial x_N}
</math>|{{EquationRef|2}}}}

若定义 <math>\frac{\partial x_k}{\partial s} = a_k \left( \boldsymbol{x}, u \right)</math>，可知
{{NumBlk|:|<math>\frac{\text{d}u}{\text{d}s} = a_1 \left( \boldsymbol{x}, u \right) \frac{\partial u}{\partial x_1} + a_2 \left( \boldsymbol{x}, u \right) \frac{\partial u}{\partial x_2} + \ldots + a_N \left( \boldsymbol{x}, u \right) \frac{\partial u}{\partial x_N} = b \left( \boldsymbol{x}, u \right)</math>|{{EquationRef|3}}}}

即，求解的偏微分方程({{EquationNote|1}})的过程变作对联立的常微分方程组作积分
{{NumBlk|:|<math>\left\{ 
   \begin{array}{rcl}
     \dfrac{\partial x_k}{\partial s} &=& a_k \left( \boldsymbol{x}, u \right)\\[1em]
     \dfrac{\text{d}u}{\text{d}s} &=& b \left( \boldsymbol{x}, u \right)
   \end{array} 
\right.</math>|{{EquationRef|4}}}}

积分过程需要给定初始条件。一般初始条件给定的形式为 <math>x</math> 空间中的流形
{{NumBlk|:|<math>g \left( \boldsymbol{x}, u \right) = 0</math>|{{EquationRef|5}}}}

将此曲面对应为 <math>s=0</math>。

设想 <math>\boldsymbol{x}</math> 和 <math>u</math> 依赖于变量 <math>\{s,t_1,t_2,...,t_{N-1}\}</math>，则 <math>\{t_1,t_2,...,t_{N-1}\}</math> 可作方程({{EquationNote|5}})中的初始值，即
{{NumBlk|:|<math>
\begin{array}{rcl}
x_1 \left( s = 0 \right) &=& h_1 \left( t_1, t_2, \ldots, t_{N - 1} \right)\\
x_2 \left( s = 0 \right) &=& h_2 \left( t_1, t_2, \ldots, t_{N - 1} \right)\\
\vdots \\
u \left( s = 0 \right) &=& v \left( t_1, t_2, \ldots, t_{N - 1} \right)
\end{array}
</math>|{{EquationRef|6}}}}

从方程组({{EquationNote|4}})和初始条件({{EquationNote|6}})确定 <math>\boldsymbol{x}</math> 和 <math>u</math> 后，可以得到解的隐式形式。如果可以解析消掉 <math>s,t_1,t_2,...,t_{N-1}</math>，则可获得显式形式的解。

==一阶偏微分方程的特征线法==
沿着一阶偏微分方程的特征线, 偏微分方程简化为一个[[常微分方程]]. 沿着特征线求出对应常微分方程的解就可以得到偏微分方程的解.

为了更好地解释这一方法, 考虑具有以下形式的偏微分方程

{{NumBlk|:|<math>a(x,y,u) \frac{\partial u}{\partial x}+b(x,y,u) \frac{\partial u}{\partial y}=c(x,y,u).</math>|{{EquationRef|1}}}}

假设解 ''u'' 已知, 考虑'''R'''<sup>3</sup>中的曲面 ''z''&nbsp;=&nbsp;''u''(''x'',''y''). 曲面的[[法向量]]为

:<math>(u_x(x,y),u_y(x,y),-1).\,</math>

那么,<ref>{{harvnb|John|1991}}</ref> 方程 ({{EquationNote|1}}) 表示向量场

:<math>(a(x,y,z),b(x,y,z),c(x,y,z))\,</math>

与曲面  ''z''&nbsp;=&nbsp;''u''(''x'',''y'') 在任意点处相切. 换句话说, 解函数的图像必定是该向量场的[[积分曲线]]的并. 这些积分曲线被称作偏微分方程的特征线.

==注记==
<div class="references-small" >
<references />
</div>

== 参考资料 ==
* {{citation|first1=Richard|last1=Courant|authorlink1=Richard Courant|first2=David|last2=Hilbert|authorlink2=David Hilbert|title=Methods of Mathematical Physics, Volume II|publisher=Wiley-Interscience|year=1962}}
*{{citation|url=http://www.jstor.org/pss/2133111|last=Delgado|first=Manuel|title=The Lagrange-Charpit Method|journal=SIAM Review|volume=39|number=2|year=1997|pages=298–304|doi=10.1137/S0036144595293534}}
* {{citation|first=Lawrence C.|last=Evans|title=Partial Differential Equations|publisher=American Mathematical Society|publication-place=Providence|year=1998|isbn=0-8218-0772-2}}
* {{citation|first=Fritz|last=John|authorlink=Fritz John|title=Partial differential equations|publisher=Springer|edition=4th |year=1991|isbn=978-0387906096}}
* {{citation|first1=A. D.|last1=Polyanin|first2=V. F.|last2=Zaitsev|first3=A.|last3=Moussiaux|title=Handbook of First Order Partial Differential Equations|publisher=Taylor & Francis|publication-place=London|year=2002|isbn=0-415-27267-X}}
* {{citation|first=A. D.|last=Polyanin|title=Handbook of Linear Partial Differential Equations for Engineers and Scientists|publisher=Chapman & Hall/CRC Press|publication-place=Boca Raton|year=2002|isbn=1-58488-299-9}}
* {{citation|last=Sarra|first=Scott|title=The Method of Characteristics with applications to Conservation Laws|journal=Journal of Online Mathematics and its Applications|year=2003|url=http://www.scottsarra.org/shock/shock.html|accessdate=2010-04-06|archive-date=2020-02-22|archive-url=https://web.archive.org/web/20200222133849/http://www.scottsarra.org/shock/shock.html|dead-url=no}}.
*{{citation|last1=Streeter|first1=VL|last2=Wylie|first2=EB|title=Fluid mechanics|publisher=McGraw-Hill Higher Education|edition=International <math>9^{th} </math> Revised|year=1998}}

== 外部链接 ==
* [http://www.scottsarra.org/shock/shock.html Prof. Scott Sarra tutorial on Method of Characteristics] {{Wayback|url=http://www.scottsarra.org/shock/shock.html |date=20200222133849 }}
* [http://www-solar.mcs.st-and.ac.uk/~alan/MT2003/PDE/node5.html Prof. Alan Hood tutorial on Method of Characteristics] {{Wayback|url=http://www-solar.mcs.st-and.ac.uk/~alan/MT2003/PDE/node5.html |date=20200221213917 }}

[[Category:偏微分方程]]