偏微分方程与特征线

偏微分方程与特征线

1函数空间的矢量场

给定一个矢量场i x i v ?=)(x v ，就在空间定义了曲线簇。比如，经过0x 点的积分曲线就可以描述为下列常微分方程的初值问题

)(x i i v x = ，n i ,...,1= 0)0(x x =

这些积分曲线就构成了曲线簇。如果形式地写出这个曲线来就是

x

vt x t v t v vt t x

t x t x x t x )exp(...)!

3!21(...!3!2)(33223

2=++++=++++= 此处x 是0时刻位置，v 是作用于x 的微分算符。

这些曲线，将空间点分成了类，也就是说每条曲线上的点属于一类。曲线集合的维数是n-1维。

矢量场的可积性

那么给定两个矢量场，就会产生两簇曲线，这两簇曲线能否组成面簇呢？我们先 看看从一点出发的曲线是否在一个曲面上的条件：从x 点出发的依此沿两簇直线运动的点若能回到来，就可以认为可以组成面。即

x x vd uc vb ua =)exp()exp()exp()(exp

如果a,b,c,d 都是1级以上的小量，这个表达式有二级以上的精度，就可以找到这样的a,b,c,d,使得方程精确满足。 按照各级展开，有 一级

0a 1111=+=+d b c

二级

v d b u c a vu uv b a )()()(222211+++=-

…

由此，得到条件

v u vu uv v u βα+=-=],[

这就是两个矢量能够构成2维子空间（曲面）的条件，著名的Frobenius 定理。 n 个矢量积分形成n 维积分只空间的条件是，任意两个矢量的对易可以写成这n 个矢量组合。 可以按照下图进行直观理解

给定m 个矢量场，他们线性组合能够形成新的矢量场。组成的矢量场空间一般称为分布。

},{是任意函数i i

i i a v a ∑=?

这个分布中任意两个矢量场对易仍然在这个分布之内，这样满足Frobenius 定理的分布称为闭分布，

????],[

他们积分可以给出m 维积分子流形。

单参数李群

一个矢量场可以构造单参数李群，一个闭分布可以构造李群。

我们先看一下单参数李群的表现，它将1维参数空间（物理上经常是时间），映射为群空间。群元素可以形式地写为算符形式

)exp(vt g t =

在表示空间中也可以写为函数变换

),(t x x g t ?=

这个函数变换是常微分方程的初值问题的解

x

x t x v t x t ==?)0,()

,(),(???

当然这个函数满足如下关系

))

),,((()),(()

()(t s x f s t x f x f g g x f g s t s t ???=+=+

比如平移群)exp(x a a g ?= 表示为 )()()exp()(a x f x f a x f g x a +=?=，

再如 转动群 ))(exp(r r n ???=θθg 表示为)()())(exp()(r r r n r '=???=f f r f g θθ 单参数李群定义了参数空间和实际空间上的变换关系和函数变换关系。

微分形式

一个函数描述为

),...,()(1n x x f x f u ==

可以看做 自变量空间到变量空间的映射

u x R R n →→:

在自变量和因变量联合空间中，可以看做一个超曲面。

如果给自变量微小改变dx x x +?，因变量也有相应的改变

dx f df x ,=

上面下标逗号表示求导。

如果想计算某个方向的导数，仅需要将相应dx 改成相应的矢量分量就可

i x i v f v df i ,=

这就是微分形式。微分形式不再依赖坐标。因此可以认为是客观量。 一般1微分形式可以描述为

i i dx x )(ωω=

不同坐标空间上的微分形式可以通过拉回映射表达出来

)(,:y x y N M n m →→?

那么n

N 空间上的微分形式可以通过映射*?拉回到m

M 空间上的微分形式

j i y i i i dy y x y x y dx y x j

)())(()())((*?==='ωωω?ω

微分形式可以与矢量作用，

i i v v i ωω=

因此可以将微分1形式想象成线元积分场，给定空间某点上一个线元，就给一个值。 当然，给定一条曲线，就可以给一个积分值

一条曲线可以描述为一维空间1

T 向n 维空间n

N 上的映射

)(,:1t x t N T l n →→

??=l

i

i t dx t x l )())((*ωω

微分形式的外积

两个微分形式θω,，相当与两个线元积分场。用这两个线元积分场可以构造一个面元积分场，要求面元大小和方向固定时，这个值是不变的。要求

θωδ

γβ

αθωδγβα∧=

∧++v u u v v u i i i i 因此，j i j i j i i j j i dx dx dx dx ∧=-=

∧θωθωθωθω)(2

1

外微分

观察微分形式ω沿无穷小闭合回路（围出来无穷小面元）的积分值，这可以定义为无穷小面元上的函数（2微分形式）

??

???Ω

Ω

Ω

?∧==j

i i j dx dx d ,ωωω j i i j dx dx d ∧=,ωω

k 形式

对微分形式进行外积或者外微分都可以变成2形式，3形式，。。。 对于m 微空间，可以证明，最高阶是m 形式。

微分形式的可积性

很明显，如果df =ω，那么有0=ωd

一个问题就是如果df ≠ω，那么能否有adf =ω，很明显ωθω∧=d 。也就是说，如果微分形式ω沿无穷小闭合回路（围出来无穷小面元）的面元积分场是由原来的面元积分场合成的，这个线元积分场就可以写成全微分乘以一个因子形式。

另一个问题是给定一些微分形式},...,1,{n =αθα

，能否判定任意一个微分形式的外微分可以表达为这些微分形式的组合形式？ 答案是：01=∧∧=β

βα

θ

θn

d

可以很容易证明这个表达式

将},...,1,{n =αθα

扩充后，形成余切空间*TM 的完备基

},...,1,,,...,1,{n m n -==μσαθμ

α

那么γγαββαασρθτθ∧+∧=d ， 可以肯定01

=∧∧=ββγ

γ

αθσρ

n

，这是关于γ

α

ρ

的线性方程，由于ββγ

θσ

n

1

=∧独立，这个方程

只有0解。 因此ββ

αα

θτ

θ∧=d

我们再看能使},...,1,{n =αθα拉回到0的映射?，

0*=αθ?

能否找到m n m M N →-:?，使得上式成立呢？

这就是Frobenius 定理的另一种描述，当任意α，都有01=∧∧=β

βα

θ

θn

d 时，可以找到

m n m M N →-:?，将},...,1,{n =αθα推回到0.

其实ββ

αα

θτ

θ∧=d 就很能说明问题，几何上讲，绕任意无穷小回路对αθ求和后，都可以

表达为},...,1,{n =αθα的组合形式。因此，使得某点的},...,1,{n =αθα为0的切向场，也可连续延拓到别处。这样的切向量场的积分曲面就是映射形成的曲面。 表达为

},...,1,0:{n i v V v ===αθα

在V 中，可以找到相互对易的m-n 个矢量},..,1,{n m w -=αα，映射可以形式地表示为

x w x M N m n m )exp()(,:ααλλλ?=→→-

很明显

0))(())(()())(()())((*,====ααααλαααλλλθλλλθλλθθ?αd x w x d x x dx x i i i i i i

这些矢量},..,1,{n m w -=αα就构成了方程},...,1,0{n ==αθα

的特征矢量。

微分形式组成的理想

如果给定生成元},...,1,{n =αθα

，我们将}0,{

形式是任意形式，包括αα

ααωθω∑∧=I 成为生成元生成的理想。很明显任意形式ω（包括函数，0形式），只要和理想中的元相乘

（外积），都会变成理想中的元素，即I I ?∧ω。这和常讲的理想意义差不多。 借用理想概念

Frobenius 定理表达为

一个外微分理想I 的具有最大零化子空间的条件是 I dI ?

偏微分方程（组）

表达为0,...),,,(=xx x u u x u F ，可以理解为函数偏导数的约束关系。 Hamilton 力学

L

H H q

H p t q q L p q H t t p q -?=??=-?=-?= ),,(

比如流体（固体）方程

0))2

1

(()21(0

)()(0)(22=?-+??++?=-??+?=??+?v σv v v σvv v v E E t t t ρρρρρρρρ，其中},,{E u u ρ=

再加上本构方程和状态方程才会封闭。 电磁学Maxwell 方程

=?+??=?-??=??=??B E D H B D t t j ρ

，其中}{E B,u =

或者在真空场写为

=?-=?αβμναβ

μνμνμε

F j F

μννμμνA A F ?-?=，其中},{A A u φ==

加上电磁学本构和电流方程才会封闭。 量子力学的薛定谔方程

ψψ))(2(2

2r V m

i t +?=? ，其中ψ=u

相对论电子运动的狄拉克方程

0)(=-

?ψγμμ

mc

i 刘维尔方程

00=??-??+?=?+?+?ρρρρρρp q q p t p q t H H p q

],[ρρH i t =?

相对论电磁学

?e A ev c v c m L -?+??

?

??-=2

2

01/

224

202

20)(1/eA p c c m e e c v c m H -++=+??

? ??-=??

2

24

202

2

24202)()()()

(eA p c c m A

eA p e c e p

eA p c c m eA p c v -+??--?-=-+-=

?

切触空间

为了从几何上描述偏微分方程的意义我们定义切触空间。

我们定义：

自变量x 的空间称为域空间D

自变量x ，因变量u 构成的空间称为图空间G

自变量x ，因变量u ，和因变量对自变量的导数p ，构成的空间},,{p u x 叫做切触空间K 。切触空间是对图空间的拓展。带来一些自然结构，即切触形式

i i dx p du C ααα-=

任何函数),(,:u x →→λφG M ，扩充为到切触空间的映射

),,(,:p u x →→ΦλK M

都会满足切触关系0=ΦC

这样，一阶偏微分方程组描述为

}),,,({βββββαdx p du p u x F PE i i i -=

如果一个映射满足0*=ΦPE ，这个映射就是0=PE 的解。

同样地可以定义高阶切触空间

...

j ij i i i

i dx p dp C dx p du C α

ααααα-=-=

高阶偏微分方程表示

,...},,...),,,,({αααααβi ij i i C C p p u x F PE =

方程解是满足 0*=ΦPE

的映射。

一阶偏微分方程（组）的特征线

一阶偏微分方程

}),,,({pdx du p u x F PE -=

为了寻找它的解法，我们寻找合适的微分形式，对函数微分，得到

},,{,,,1pdx du dp F du F dx F F PE p u x -++=

很自然地想到微分形式组合的特性矢量，就是}0:{1=PE i w w 的矢量。这里有一个问题需要解决，1PE 封闭吗？也就是说是否满足)(11PE I dPE ?？

很明显0=ddF ，但是dp dx dC ∧==ω并不在理想中，因此}0:{1=PE i w w 的矢量，有可能不能够积分出一个子空间来，因此不是偏微分方程的解。

},,,{,,,2dp dx pdx du dp F du F dx F F PE p u x ∧-++=是封闭的，定义 ),,,()(,,,dp dx pdx du dp F du F dx F F I PE I p u x ∧-++=

)}()(:{PE I PE I i w w ?的矢量就是偏微分方程的特征矢量，它们的积分组成偏微分方程的

解。

我们考虑只有一个因变量情况0),,(=i i p u x H 的偏微分方程的

),,,()(,,,dp dx pdx du dp H du H dx H H I PE I p u x ∧-++=

设i i x i

u p i c b a w ?+?+?=是上述理想的特征矢量 因此，i i i i i i w dp c dx a dp dx i +-=∧)(，可以描述为

??

?-++pdx

du dp

H du H dx H p u x ,,, 线性组合的形式。消去du 的项，得到，

i p i i u i x dp H dx p H dx H i i ,,,++，

因此可以选i

i p i i u x i H c p H H a ,,,)(=+-=，再有0)(=-=-i i i i w c p b dx p du i ，得到i p i H p b ,=

因此u p i x p p i u x i i i i i H p H p H H w ?+?+?+-=

,,,)(

因此特征线可以采用常微分方程积分

i i i p i p i i u x i H p u

H x p H H p ==+-= ,,,)(

如果u 当成作用量，H 当成Hamilton 量，x 粒子位置，p 粒子动量，这正是经典粒子运动方程。

李导数

一个由矢量i x i v v ?=形成的单参数李群可以表述为m m

M M →的映射，也是m M R →的

映射

x vt x x t t )exp(),(:=→φ。

对于一个函数)(x f ，这个映射可以将其拉回到)(*t x f φ，即))(exp()()exp()(*x vt f x f vt x f t ==φ，后者可以方便地用算符运算验证

))(exp(1))exp()(exp(1)exp()()exp()()exp(x vt f vt x vt f vt x f vt x f vt =-=-=

我们计算函数的无穷小变换即李导数

df i x vf t

x f x f vt f L v v ==-=

)()

()()exp(

也就是说，函数的李导数就是其方向导数 对于微分形式i i dx ωω=，

))(exp())(exp()(*x vt d x vt dx x i t t i ωωωφ==

ω

ωωωωωωωωωωω)()()()()()()()())(exp())(exp(v v v v i i i i v i i i i v i i i i v di d i di d i v d dx d i dv dx d i t

dx x vxt x d vxt x t

dx

x x vt d x vt L +=+=+∧=+=-++=-=

这个公式也可包含前一个，只要认为0=f i v 就可以了，因此

v v v di d i L +=

可以证明，李导数满足莱布尼兹法则

θωθωθωv v v L L L ∧+∧=∧)(

接着我们推导矢量场的李导数

i x i u u ?=

在t x 处的矢量要与x 处的矢量比较，首先必须通过映射将其映射到x 点上，映射是

x vt x x t )exp(=→，这个映射诱导的映射只能将x 处的矢量推向t x 。因此采用这个映射的

逆映射将t x 处的矢量推向x 。t t x vt x x )exp(:-=→φ

]

,[)()()()(*)()(v u v u u v t

x u x vxt x u t

x u x u u L i x j i x j x i x j t x i x

i t x t i v j j i

j i i

i

=?-?=?-??+=

?-?=φ

微分方程的对称性

对于在切触空间中微分形式理想I 表达的微分方程，如果存在坐标变换（映射）使得I 不变，这个变换就是微分方程的对称性。我们考虑一种连续变换，它形成单参数李群，写为

))exp(,)exp(,)(exp(),,(,:p Vt u Vt x Vt p u x t →→K K ?

其中ααα

αi

i p i u x i V V V V ?+?+?=表示切触空间中的矢量。

因此我们需要研究这个的无穷小变换将微分形式如何改变。

实际上如果I 的李导数仍然属于I 就可以保证李变换是微分方程的对称变换。

β

αβααααC

A C L dx p du C V i

i =-=

i

i i i i i i i

i

i i i i i V v V V dV

p dV dx V V p V d dp V dx V V p V d dx dp i C di d i C L ααααααα

ααααα-+-=-++-=-+∧-=+=)()()()(

i i dx p A du A C A βαββαββαβ-=

对比βj dp ，βdu ,i

dx 的系数，可以得到

0=?-?i j i j V p V βααβ αβααβA V p V j j =-?)(

βαβαααi j j i i p A V p V V -=-?+-)(

约去系数A,

)()()(j j i i j j j j i i V p V D p V p V V p V V ααβααβααα-=-?+-?= 0=?-?i j i j V p V βααβ

由上式兼容性，可以得到αV V i ,与β

j P 的依赖关系 写为

j i i j V V p V αβααβδ-=-?)(

兼容性条件是

k j i i k j j k i i j k V V p V V V p V βαγααγβγαβααβγδδ?-=-??=?-=-??)()(

因变量数目大于1时，选βγα≠=，有k j V β-?=0，再代入前公式0=?αβV j 。也就是说在因变量数目大于1时i V V ,α仅是),(u x 的任意函数，而)(j j i i V p V D V ααα-=。 对于只有一个因变量情况k j j k V V 11?=?，这表明j i i j V V p V -=-?)(111

令i i V V p V C i 11-==ξ

ξ

ξξξj j p j p j D V p V V j j =?-=-?=11 对照特征矢量，可以看出，原来这个对称性仅仅是将1个自变量的微分方程作了规范变换。

这就是满足切触条件必须满足的对称性约束。 如果知道方程对称性，比如V

就可以根据方程的解构造新的解了，比如m n n m s C R ++→:φ是方程的一只解，那么

)exp(V t φ也是方程的一只解，且m n n m s C R Vt +++→1:)exp( φ，这样可以从低维解子空

间构造高维解子空间了。

对于常微分方程，可以减少未知数或降阶。

常微分方程的对称性及其精确解

给定常微分方程组，确定其微分形式理想。如果确定一个对称性，就可以将其降阶。具体地， 如果确定了其对称矢量v ，并解析求出其李变换场

)),,(),,,(),,,((),,(p u x p p u x u p u x x p u x t t t →，将t 消去，可以得到n-1个表达式，加上t,

以它们为变量，方程不再包含t （由于在对称性变换时，有a t t +→，又由于系统具有对称性，因此，新变量不能显含t 。），变量数目就减少了1个。 比如其次方程

)/(x y f y =',有1个因变量和一个自变量，方程对称性为 )exp(),exp(00t y y t x x ==,消

去t,得到表达式x y w /=，)/ln(0x x t =将这个表达式带入方程，方程变得0),(=t w w F ，不再含有未知数t,因此减少了一个未知数（可能是自变量，也可能是自变量）。对于目前方程，可以解析求解了

)(w f w w t =+

)/ln()(0x x t w f w dw

==-?