system: [SYSTEM]

user: 4. 关于性质3)的讨论 性质3)中轨线族的整体性质，也就是轨线族的整体结构，它由轨线族中轨线的演变方式所决定。微分系统的轨线不是孤立的，每一条轨线都是另一条轨线按照每个微分系统各自特有的方式演变而成。在轨线的演变过程中，轨线上每一点的运动轨迹构成了平面微分系统的一条经线。记M为平面微分系统的轨线族，N为经线族。同轨线族M一样，每个光滑平面微分系统的经线族N存在且唯一。 例2 系统 d x d t = y b 2 ， d y d t = − x a 2 ( a b ≠ 0 ) (6) 的轨线族 M = { x 2 a 2 + y 2 b 2 = c ; 0 < c < + ∞ } ，经线族 N = { x = 0 ， y = c x ; − ∞ < c < + ∞ } 。 例3 系统 d x d t = a x + b y ， d y d t = − b x + a y ( a b ≠ 0 ) (7) 的轨线族(用极坐标表示) M = { ρ = c e − a b φ ; 0 < c < + ∞ } ，经线族 N = { ρ = c ; 0 < c < + ∞ } 。 从例2，例3可以看出： 1) 光滑的平面微分系统的经线族N存在且唯一； 2) 轨线上的点沿经线运动时，轨线的局部性质变化最小； 3) 系统的经线反映了系统轨线的演变方式即轨线族的整体性质。系统(5)轨线上的点沿射线 y = c x 运动演变另一条轨线，这种系统我们称为相似系统。系统(7)轨线上的点沿圆周 ρ = c 运动演变成另一条轨线，这种系统我们称为旋转系统。 4) 通过系统的经线，可以分析系统轨线的分布密度，了解微分方程解的稳定性； 5) M ∩ N = ϕ 。 这5条性质，前4条对所有的平面微分系统都成立，唯有性质5)不是对所有微分系统都成立。我们认为：若平面微分系统存在极限环L，则必有 L ∈ M ∩ N 。即存在极限环的微分系统的 M ∩ N 一定非空，而且极限环必定属于 M ∩ N 。这就是存在极限环的平面微分系统轨线族具有特殊整体性质3)。根据这一性质，我们只要求出微分系统的经线族N，就可以从 M ∩ N 中得到系统的所有极限环，从而彻底解决平面微分系统极限环的存在性以及存在个数问题。 如何求平面微分系统的经线呢？根据经线的性质2)，我们考虑： 定义9 设 G = ( g 1 ， g 2 ) 是系统(1)轨线上点 ( x ， y ) 的基因。在以点 ( x ， y ) 为中心， ε 为半径的圆周： x + Δ x = x + ε cos φ ， y + Δ y = y + ε sin φ (8) 上的点 ( x + Δ x ， y + Δ y ) 的基因为 G + Δ G = ( g 1 + Δ g 1 ， g 2 + Δ g 2 ) ，则 | Δ G | = ( Δ g 1 ) 2 + ( Δ g 2 ) 2 在圆周(8)上一定取到最小值。即在圆周(8)上一定存在一点 ( x + Δ x ， y + Δ y ) 是 | Δ G | 的最小值点。当 ε → 0 时，这些最小值点构成一条经过点 ( x ， y ) 的曲线我们称为基因G的遗传曲线。 根据定义9，可求得系统(1)点 ( x ， y ) 基因G的遗传曲线的斜率 d y d x 由方程： ( d y d x ) 2 + H d y d x − 1 = 0 (9) 所确定。其中 H = g 1 x 2 − g 1 y 2 + g 2 x 2 − g 2 y 2 g 1 x g 1 y + g 2 x g 2 y (10) 方程(9)有两个根，其中一个 d y d x 就是基因G变化最小的方向，另一个与 d y d x 正交的 − ( d y d x ) − 1 就是基因G变化最大的方向。 有时轨线及基因的遗传曲线用极坐标方程表示更方便。此时基因G遗传曲线 ρ = ρ ( φ ) 的导数 d ρ d φ 由方程 ( d ρ d φ ) 2 + H d ρ d φ − ρ 2 = 0 (11) 确定。其中 H = − ρ 2 ( g 1 ρ 2 + g 2 ρ 2 ) + ( g 1 φ 2 + g 2 φ 2 ) g 1 ρ g 1 φ + g 2 ρ g 2 φ (12) 定义9基因G的遗传曲线正是例2系统(5)，例3系统(7)的经线。这类以基因G的遗传曲线为经线的系统，我们称为简单系统。而大多数平面微分系统是复杂系统，它们的经线不是基因G的遗传曲线，这是为什么呢？因为G中 g 1 = r 0 sin α 是轨线的拱度， g 2 = r 0 cos α 是轨线的拱偏度，是两种截然不同的几何性质。在轨线演变过程中， g 1 与 g 2 有不同的变化规律，在定义9中把它们象点的坐标的两个分量一样看待，用定义点的距离大小的方法定义轨线性质变化大小是不合适的，在复杂系统中不能真实反映轨线演变过程中点的运动轨迹。于是如何求一般平面微分系统的经线是我们面临的一个问题，本文从以下几个方面进行探讨。 4.1. 曲线基因表示的坐标化 要用定义9定义基因的变化大小就必须使得表示基因的各个量，除了表示距离以外不能有其它几何意义，也就是做到基因的坐标化表示。 如何做到基因表示的坐标化呢？从文 [ 1 ] 的定理1，2，3，4我们知道：足够光滑的平面曲线在它的每一点处有一条圆锥曲线为它的4阶密切曲线，也就是平面曲线与它密切圆锥曲线在切点的基因是相同的，所以只要解决了圆锥曲线基因的坐标化，就解决了所有平面曲线基因的坐标化。又因为圆锥曲线族{ ρ = c 1 − e cos φ ；e是常数， 0 < c < + ∞ }是相似曲线族，我们只要解决其中一条曲线 ρ = c 0 1 − e cos φ 的基因坐标化，就解决了整个曲线族的曲线基因坐标化。为此我们考虑系统： d u d t = v ( u + 1 ) ， d v d t = v 2 − u (13) 与系统 d u d t = v ( u − 1 ) ， d v d t = v 2 − u 。由于两个系统讨论相同，下面以系统(13)为例进行讨论。系统(13)的轨线族 M e = { u 2 + v 2 u + 1 = e ; 0 < e < + ∞ } (14) 用极坐标方程表示 M e = { ρ = e 1 − e co s φ ; 0 < e < + ∞ } 。 M e 中的曲线称为e曲线，记做 L e 。系统(13)中点的坐标用 ( u ， v ) 表示，是因为以后我们经常要用到这一特殊的系统。 系统(13)的轨线族 M e 有以下性质： 1) M e 中的曲线 L e 都以 u = − 1 为准线，原点0为焦点； 2) M e 包含了所有类型的圆锥曲线； 3) M e 中曲线的离心率e两两不相同； 4) 光滑平面曲线任一点的密切圆锥曲线必定和 M e 中唯一一条曲线 L e 相似(有相同的离心率)。 因此我们只要解决了 M e 中曲线 L e 的基因坐标化，就解决了所有平面曲线基因的坐标化。 下面讨论 M e 中曲线基因的坐标表示。 L e 上的点 ( u ， v ) 到焦点 ( 0 ， 0 ) 的距离 u 2 + v 2 ，到准线 u = − 1 的距离为 u + 1 ， u 2 + v 2 u + 1 = e 的点全体就是 L e 。所以 ( u 2 + v 2 ， u + 1 ) 不仅可以确定点 ( u ， v ) 在 L e 上的位置，而且可以确定 L e 。因此有： 定义10 G 1 ( u ， v ) = ( u 2 + v 2 ， u + 1 ) 称为 L e 上点 ( u ， v ) 的线状基因。 另外曲线 L e 也是点 ( u ， v ) 到焦点 ( 0 ， 0 ) 的距离 u 2 + v 2 与点 ( u ， v ) 到另一焦点 ( 2 U ， 0 ) 的距离 ( u − 2 U ) 2 + v 2 之和(或之差)为常数的点的集合(其中 U = e 2 1 − e 2 )，也就是以点 ( u ， v ) ， ( 0 ， 0 ) ， ( 2 U ， 0 ) 为顶点三角形也完全确定了曲线 L e ，我们把这个三角形三个顶点坐标合在一起，并略去其中的常数0，有： 定义11 G 2 ( u ， v ) = ( u ， v ， 2 U ) 称为曲线 L e 上点 ( u ， v ) 的点状基因。 这样我们就得曲线 L e 点基因的两种坐标表示： G 1 ， G 2 。下面我们讨论足够光滑的平面曲线 L : f ( x ， y ) = 0 上点 ( x ， y ) 的基因坐标表示。 定义12 设点 ( x ， y ) 是曲线L上任一点，它的基因是 G ( x ， y ) ，根据曲线族 M e 的性质4，必存在唯一的 L e ∈ M e ， L e 有点 ( u ， v ) 的基因 G ( u ， v ) = G ( x ， y ) ，我们称 G i ( u ， v ) 也为曲线L上点 ( x ， y ) 的 G i 基因。即 G i ( x ， y ) = G i ( u ， v ) ( i = 1 ， 2 ) 。 根据定义12： G i ( x ， y ) = G i ( u ， v ) ⇔ G ( x ， y ) = G ( u ， v ) ⇔ g 1 ( x ， y ) = g 1 ( u ， v ) 且 g 2 ( x ， y ) = g 2 ( u ， v ) ⇔ R 1 ( x ， y ) = R 1 ( u ， v ) 且 R ( x ， y ) R 2 ( x ， y ) = R ( u ， v ) R 2 ( u ， v ) (15) 根据(15)式，有 定理1 足够光滑的平面曲线 L : f ( x ， y ) = 0 在点 ( x ， y ) 的局部性质由： R 1 ( x ， y ) 与 R ( x ， y ) R 2 ( x ， y ) 所决定。 下面求 L e 上点 ( u ， v ) 的 u ， v 与 R 1 ( u ， v ) ， R ( u ， v ) R 2 ( u ， v ) 的关系 从 u 2 + v 2 u + 1 = e ， d v d u = v 2 − u v ( u + 1 ) ，可得 R 1 ( u ， v ) = − 3 v ( v 2 − u ) u + 1 (16) R 1 2 ( u ， v ) = 9 v 2 [ ( 1 + v 2 ) e 2 − v 2 ] (17) R ( u ， v ) R 2 ( u ， v ) = 3 ( 1 + v 2 ) [ e 2 ( 1 + 2 v 2 ) − 2 v 2 ] (18) 从(17)，(18)两式消去 e 2 可得： 9 v 4 + ( 3 R R 2 − 2 R 1 2 ) v 2 − R 1 2 = 0 (19) 从(17)，(18)两式还可得： e 2 = 2 ( 1 + v 2 ) R 1 2 − 3 v 2 R R 2 9 v 2 ( 1 + v 2 ) (20) 解方程(19)，并取定一根 v = v ( R 1 ( u ， v ) ， R ( u ， v ) R 2 ( u ， v ) ) = 1 6 2 ( 2 R 1 2 − 3 R R 2 ) + 2 ( 3 R R 2 − 2 R 1 2 ) 2 + 36 R 1 2 (21) 代入(16)式： u = R 1 + 3 v 3 3 v − R 1 = u ( R 1 ( u ， v ) ， R ( u ， v ) R 2 ( u ， v ) ) (22) 将(21)，(221)式中 R 1 ( u ， v ) ， R ( u ， v ) R 2 ( u ， v ) 分别换成曲线L点 ( x ， y ) 的 R 1 ( x ， y ) ， R ( x ， y ) R 2 ( x ， y ) 得： v = v ( x ， y ) = v ( R 1 ( x ， y ) ， R ( x ， y ) R 2 ( x ， y ) ) (23) u = u ( x ， y ) = u ( R 1 ( x ， y ) ， R ( x ， y ) R 2 ( x ， y ) ) (24) 从(20)式还可得 e = e ( x ， y ) = e ( R 1 ( x ， y ) ， R ( x ， y ) R 2 ( x ， y ) ) (25) 于是得到曲线上点 ( x ， y ) 的线状基因 G 1 ( x ， y ) = ( u 2 ( x ， y ) + v 2 ( x ， y ) ， u ( x ， y ) + 1 ) (26) 点状基因 G 2 ( x ， y ) = ( u ( x ， y ) ， v ( x ， y ) ， 2 U ( x ， y ) ) (27) 其中 U ( x y ) = e 2 ( x ， y ) 1 − e 2 ( x ， y ) 。 定义13 将定义9中点 ( x ， y ) 的基因 G ( x ， y ) 换成点 ( x ， y ) 的基因 G i ( x ， y ) ，得到的曲线称为系统(1)基因 G i 的遗传曲线。 G i 遗传曲线的斜率 d y d x 由方程： ( d y d x ) 2 + H i d y d x − 1 = 0 (28) 确定。 G i 遗传曲线的全体记为 N i ( i = 1 ， 2 ) 。 其中 H 1 = ( u 2 + v 2 ) x 2 − ( u 2 + v 2 ) y 2 + 1 ( u 2 + v 2 ) x ( u 2 + v 2 ) y H 2 = ( u ) x 2 − ( u ) y 2 + ( v ) x 2 − ( v ) y 2 + 4 U x 2 − 4 U y 2 ( u ) x ( u ) y + ( v ) x ( v ) y + 4 U x U y 例4 求系统(13)的 N 1 ， N 2 。 为方便，系统(13)的 G 1 ( u ， v ) = ( u 2 + v 2 ， u + 1 ) 用极坐标表示： G 1 = ( ρ ， 1 + ρ cos φ ) 。由(12)式有： H 1 = − ρ 2 ( 1 + cos 2 φ ) + ρ 2 sin 2 φ − ρ sin φ cos φ = 2 ρ cos φ sin φ 于是 ( d ρ 2 d φ ) + 2 ρ cos φ sin φ − ρ 2 = 0 。 解之得 d ρ d φ = ρ ( cos ρ ± 1 ) sin ρ 。 于是 ρ = c 1 + cos φ ， ρ = c 1 − cos ρ ， 0 < c < + ∞ 。 显然 ρ   = c 1 + cos φ 不是G 1 的遗传曲线，所以 N 1   =   { ρ = c 1 − cos φ ; 0 < c < + ∞ } ， M e   ∩ N 1   = { ρ = 1 1 − cos φ } 从 G 2 ( u ， v ) = ( u ， v ， 2 U ) 得： H 2 = U u 2 − U v 2 U u U v ， ( d v d u ) 2 + U u 2 − U v 2 U u U v d v d u − 1 = 0 。 解之得 d v d u = − U u U v (另一解 d v d u = U v U u )于是 U = e 2 1 − e 2 = 常 数 ，即 N 2 = M e 。 从例4，我们可以分析得到系统(1) M ∩ N 的一些性质： 1) M ∩ N 将M分成若干个独立的子系统； 2) 每个独立的子系统只含有椭圆型轨线，抛物型轨线，双曲型轨线之中的一种轨线，不同型轨线不在同一子系统中。 在分析 G 1 与 G 2 哪个更适合作为基因的坐标表示时，我们发现：系统中基因越接近的点，它们的密切中心 ( X ， Y ) 相隔的距离也越小。根据这一现象，有： 定义14 将定义9点 ( x ， y ) 的基因G换成点 ( x ， y ) 的密切中心 ( X ， Y ) ，得到的曲线称为系统(1)的中心 遗传曲线，它们的全体记为 N 3 。系统(1)中心遗传曲线斜率 d y d x 由方程： ( d y d x ) 2 + H 3 d y d x − 1 = 0 (28) 确定。 H 3 = X x 2 − X y 2 + Y x 2 − Y y 2 X x X y + Y x Y y 。 对于系统(12)而言，点 ( u ， v ) 的密切中心为 ( U ， 0 ) 。所以 H 3 = U u 2 − U v 2 U u U v = H 2 ，于是 N 3 = N 2 = M e 虽然点 ( x ， y ) 的密切中心 ( X ， Y ) 并不表示点 ( x ， y ) 处轨线的性质，但是它在轨线演变过程中，有其自身的变化规律。如能掌握，也能帮助我们找到系统的经线。除了用定义9，我们还可以其它方面考虑系统的经线问题。 4.2. g 1 ，g 2 的梯度与基因突变环 系统(1)轨线在点 ( x ， y ) 处 g r a d g i = ( g i x ， g i y ) 的方向是 g i 变化率最大的方向， g i 沿方向 l → β = ( cos β ， sin β ) 的方向导数为 ∂ g i ∂ l β = g i x cos β + g i y sin β ( i = 1 ， 2 ) 。 定义15 若方向 l → β = ( cos β ， sin β ) 满足： ∂ g 1 ∂ l β / ∂ g 2 ∂ l β = | g r a d g 1 | | g r a d g 2 | 则称方向 l → β 为基因G的(与梯度)同步遗传方向，方程 d y d x = t g β (29) 经过点 ( x ， y ) 的积分曲线称为基因G的同步遗传曲线，它们的全体记为 N 4 。其中 t g β = − g 1 x g 2 x 2 + g 2 y 2 − g 2 x g 1 x 2 + g 1 y 2 g 1 y g 2 x 2 + g 2 y 2 − g 2 y g 1 x 2 + g 1 y 2 显然 M ∩ N 4 = { t g β = q p } 。系统(1)的 N 1 ， N 2 ， N 3 一般比较难求，但 M ∩ N i 都和 M ∩ N 4 一样容易求，因为 M ∩ N i = { H i = p 2 − q 2 p q } ( i = 1 ， 2 ， 3 ) 。 定义16 若系统(1)的轨线 L ∈ M ∩ N i ，则称L系统(1)的 N i 分界环 ( i = 1 ， 2 ， 3 ， 4 ) 。 注意 N i 分界环不一定是极限环。 通过 N 4 的定义可以看出，等拱线： g 1 = 常 数 ，等偏拱线： g 2 = 常 数 ，等径线： r 0 = 常 数 都是系统(1)的重要曲线。 定义16 Λ 1 ( x ， y ) = r 0 x g 1 y − r 0 y g 1 x r 0 x g 1 x + r 0 y g 1 y (30) 称为系统(1)在点 ( x ， y ) 的拱变度。 Λ 2 ( x ， y ) = r 0 x g 2 y − r 0 y g 2 x r 0 x g 2 x + r 0 y g 2 y (31) 称为系统(1)点 ( x ， y ) 的偏拱变度。 Λ ( x ， y ) = g 1 x g 2 y − g 1 y g 2 x g 1 x g 2 x + g 1 y g 2 y (32) 称为系统(1)点 ( x ， y ) 的基因变异度。 如何通过这些变异度去求系统(1)的经线暂不讨论。我们考虑极限环： 定义17 系统(1)中， Λ ( x ， y ) 在等径线 r 0 = 常 数 上的极值点称为系统(1)的基因突变点.由基因突变点组成的曲线称为系统(1)的基因突变线，其全体记为 N * 。 显然 N * 不是系统(1)的经线族，但很重要，因为我们猜测：平面微分系统的任一极限环L都是系统的基因突变线，即任一极限环 L ∈ N * 。 4.3. 平面曲线的全盈度 我们在文 [ 2 ] 定义了平面曲线的盈度，并指出相对于拱度，盈度应用起来不方便，是因为没有合适的量与盈度配对完整地表示曲线的局部性质。下面解决这一问题。由于文 [ 2 ] 定义1.3有一点问题，我们重新定义曲线的盈度。 设L是足够光滑的严格凸平面曲线，O点是曲线L上一点，以过O点的切线为x轴，法线为y轴，如图1建立直角坐标系。设曲线L的方程为： y = y ( x ) 。 图1. 曲线 y = y ( x ) 过O点附近曲线L点 M ( x ， y ) 作x轴的平行线与曲线L交于另一点 M 1 ( x 1 ， y ) ，有 y ( x 1 ) = y ( x ) ，由此确定 x 1 是x的函数并且由文 [ 2 ] 有： x ′ 1 ( 0 ) = − 1 ， x ″ 1 ( 0 ) = − 2 y ‴ ( 0 ) 3 y ″ ( 0 ) ， x ‴ 1 ( 0 ) = − 2 y ‴ 2 ( 0 ) 3 y ″ 2 ( 0 ) 。 曲线L在O点的密切圆： f ( x ) = R 0 − R 0 2 − x 2 ， R 0 = 1 y ″ ( 0 ) 。 在区间 [ x 1 ， x ] 上曲线L的拱高为 y ( x ) ，弦长为 x − x 1 ，在区间 [ − x − x 1 2 ， x − x 1 2 ] 上密切圆的拱高为 f ( x − x 1 2 ) ，弦长也是 x − x 1 。两拱高之差与点O处曲线的曲率 ( R 0 ) − 1 之比 y ( x ) − f ( x − x 1 2 ) ( R 0 ) − 1 = 0 ( y 2 ) 。 因此有： 定义18 ω 1 ( 0 ， 0 ) = lim x → 0 R 0 [ y ( x ) − f ( x − x 1 2 ) ] y 2 ( x ) = 3 y ″ ( 0 ) y ( 4 ) ( 0 ) − 9 y ″ 4 ( 0 ) − 5 y ‴ 2 ( 0 ) 18 y ″ 4 ( 0 ) (32) 称为曲线L在O点(相对于圆弧)的主盈度。 因为 ω 1 ( 0 ， 0 ) = 1 18 ( R 1 2 ( 0 ， 0 ) − 3 R ( 0 ， 0 ) R 2 ( 0 ， 0 ) ) 。 所以称 ω 1 ( x ， y ) = 1 18 ( R 1 2 ( x ， y ) − 3 R ( x ， y ) R 2 ( x ， y ) ) (33) 为曲线L在点 ( x ， y ) 的主盈度。 下面考虑曲线L拱的弦长与密切圆拱的弦长的比较。 设密切圆在 x 0 点的拱高 f ( x 0 ) 与曲线L在x点的拱高相等： f ( x 0 ) = y ( x ) 。由此确定 x 0 与x的函数关系，对 f ( x 0 ) = y ( x ) 两边关于x求导： f ′ ( x 0 ) x ′ 0 = y ′ ( x ) 再求导： f ″ ( x 0 ) x ′ 0 2 + f ( x 0 ) x ″ 0 = y ″ ( x ) 继续求导： f ‴ ( x 0 ) x ′ 0 3 + 3 f ″ ( x 0 ) x 0 x ″ 0 + f ′ ( x 0 ) x ‴ 0 = y ‴ ( x ) 可得 x ′ 0 ( 0 ) = 1 ， x ″ 0 ( 0 ) = y ‴ ( 0 ) 3 y ″ ( 0 ) ， x ‴ 0 ( 0 ) = 3 y ″ ( 0 ) y ( 4 ) ( 0 ) − 9 y ″ 4 ( 0 ) − y ‴ 2 ( 0 ) 12 y ″ 2 ( 0 ) 。 在两拱高相等的情况下，密切圆拱的弦长 2 x 0 ，并且 x 0 2 2 = f ( x 0 ) ( R 0 − f ( x 0 ) 2 ) 。而曲线L拱的弦长为 x − x 1 ，该弦被y轴分成两部分长度的乘积为 | x 1 | x ，因为 | x 1 | x 2 − x 0 2 2 = 0 ( y 2 ) ，所以有 定义19 ω 2 ( 0 ， 0 ) = lim x → 0 | x 1 | x − x 0 2 2 y 2 ( x ) = 1 18 ( 3 R ( 0 ， 0 ) R 2 ( 0 ， 0 ) − 2 R 1 2 ( 0 ， 0 ) ) (34) 称为曲线L在点O的偏盈度。 ω 2 ( x ， y ) = 1 18 ( 3 R ( x ， y ) ) R 2 ( x ， y ) − 2 R 1 2 ( x ， y ) (35) 为曲线L在点 ( x ， y ) 的偏盈度。 Ω ( x ， y ) = ( ω 1 ( x ， y ) ， ω 2 ( x ， y ) ) 为曲线L在点 ( x ， y ) 的全盈度。 显然曲线L点 ( x ， y ) 的全盈度 Ω ( x ， y ) 完全确定了曲线L点 ( x ， y ) 的局部性质，是曲线基因的又一表示式。 讨论了曲线基因这么多不同表示，并定义了平面微分系统的曲线族N 1 ，N 2 ，N 3 ，N 4 。它们中有没有系统的经线族N？还有待确认。即使它们都不是系统的经线族N，也值得我们研究，因为它们从不同方面展示了平面微分系统轨线演变的规律。基于这一点有： 定义19系统(1)中曲线 Λ ( x ， y ) = 常 数 称为系统(1)的基因同变异线，记 N 5 = { Λ ( x ， y ) = c ; − ∞ < c < + ∞ } 显然通过 N 5 ，我们可以知道系统(1)变异相同的点的运动轨迹。 本文根据存在极限环的系统的三个特殊性质1)，2)，3)，讨论了极限环的求解问题，得到的结果是把百年未解决的问题：平面微分系统极限环的存在性与存在个数问题，归结为下面问题之一； 1) 求平面微分系统的经线族N； 2) 求出相似环的所有种类； 3) 证明任一极限环 L ∈ M ∩ N * 。 对这三个问题我们另文讨论。 作为本文的结束，给出本文结果的一个应用。 例5 求方程组 d x d t = y − F ( x ) ， d y d t = − x ，其中 F ( x ) = μ ( x 3 3 − x ) (36) 的 N 3 分界环。 方程组与vanderPol方程等价 [ 3 ]。 从(36)式，得 y ′ = d y d x = − x y − F y ″ = − ( y − F ) 2 + x F ′ ( y − F ) + x 2 ( y − F ) 3 y ‴ = − ( 2 F ′ + x F ″ ) ( y − F ) 3 + ( 3 x + 2 x F ′ 2 ) ( y − F ) 2 + 5 x 2 F ′ ( y − F ) + 3 x 3 ( y − F ) 5 y ( 4 ) = − ( 3 F ″ + x F ‴ ) ( y − F ) 5 + ( 3 + 6 F ′ 2 + 6 x F ′ F ″ ) ( y − F ) 4 + ( 23 x F ′ + 7 x 2 F ″ + 6 x F ′ 3 ) ( y − F ) 3 ( y − F ) 7                 + ( 18 x 2 + 26 x 2 F ′ 2 ) ( y − F ) 2 + 35 x 3 F ′ ( y − F ) + 15 x 4 ( y − F ) 7 其中 F ′ = μ ( x 2 − 1 ) ， F ″ = 2 μ x ， F ‴ = 2 μ 。 代入文 [ 2 ] 的(3)，(4)两式，得到方程(36)中点 ( x ， y ) 的密切中心 ( X ， Y ) 为： X = x − 3 y ″ y ‴ 3 y ″ y ( 4 ) − 5 y ‴ 2 = x − B A Y = y − 3 y ″ ( y ′ y ‴ − 3 y ″ 2 ) 3 y ″ y ( 4 ) − 5 y ‴ 2 = y − C A 其中 A = ∑ i = 0 6 a i ( y − F ) 6 − i B = ∑ i = 0 5 b i ( y − F ) 6 − i C = ∑ i = 0 5 c i ( y − F ) 6 − i a 0 = 9 F ″ + x F ‴ a 1 = 9 − 2 F ′ + 7 x F ′ F ″ + 3 x 2 F ′ F ‴ a 2 = 18 x F ′ − 4 x F ′ 3 − 2 x 2 F ′ 2 F ″ + 3 x 3 F ‴ a i = 0 ( x 2 ) ， i = 3 ， 4 ， 5 . b 0 = 6 F ′ + 3 x F ″ b 1 = 9 x + 12 x F ′ 2 + 3 x 2 F ′ F ″ b i = 0 ( x 2 ) ， i = 2 ， 3 ， 4 ， 5 . c 0 = 9 c 1 = 21 x F ′ − 3 x 2 F ″ c i = 0 ( x 2 ) ， i = 2 ， 3 ， 4 ， 5 因为 p ( x ， y ) = y − F ， q ( x ， y ) = − x 。 设轨线 L ∈ M ∩ N 3 ，则轨线L的斜率 d y d x = q p 和正交曲线斜率 d y d x = − p q 是方程(28)的两个根。 于是 H 3 = X x 2 − X y 2 + Y x 2 − Y y 2 X x X y + Y x Y y = p q − q p 。 即 X x 2 − X y 2 + Y x 2 − Y y 2 X x X y Y x Y y = ( y − F ) 2 − x 2 x ( y − F ) 。 [ ( y − F ) 2 − x ] ( X x X y + Y x Y y ) − x ( y − F ) ( X x 2 − X y 2 + Y x 2 − Y y 2 ) = 0 这是关于 y − F 的23次方程。 令 x = 0 ，得 y 2 = 1 4 ( 1 − 2 μ 2 9 ) 。 当 μ 2 < 9 2 时，方程组(36)有 N 3 分界环L，该环与y轴的两个交点为： y = 1 2 1 − 2 μ 2 9 ， y = − 1 2 1 − 2 μ 2 9 ，这说明该 N 3 分界环L是孤立的闭轨线，即它是方程(36)极限环。