Loading [MathJax]/jax/element/mml/optable/BasicLatin.js

数学物理学报, 2025, 45(2): 554-566

平面三次多项式的分类

李林1, 刘玲伶2, 余志恒,3,*

1嘉兴大学数学系 浙江嘉兴 314001

2西南石油大学理学院 成都 610500

3西南交通大学数学学院 成都 611756

On Classification of Planar Cubic Polynomials

Li Lin1, Liu Lingling2, Yu Zhiheng,3,*

1Faculty of Mathematics, Jiaxing University, Zhejiang Jiaxing 314001

2School of Sciences, Southwest Petroleum University, Chengdu 610500

3School of Mathematics, Southwest Jiaotong University, Chengdu 611756

通讯作者: * 余志恒,E-mail:yuzhiheng9@163.com.

收稿日期: 2023-10-31   修回日期: 2024-10-15  

基金资助: 国家自然科学基金(11871041)
国家自然科学基金(12171337)
四川省自然科学基金(2022NSFSC1834)
四川省自然科学基金(2023NSFSC0064)
中央高校基本科研业务费(2682021ZTPY041)
四川省中央引导地方科技发展专项(2024ZYD0059)

Received: 2023-10-31   Revised: 2024-10-15  

Fund supported: NSFC(11871041)
NSFC(12171337)
National Science Foundation of Sichuan Province(2022NSFSC1834)
National Science Foundation of Sichuan Province(2023NSFSC0064)
Fundamental Research Funds for the Central Universities(2682021ZTPY041)
Central Government Guided Local Science and Technology Development Projects(2024ZYD0059)

摘要

该文将讨论平面多项式的正规形理论. 作者利用多项式代数理论求出相应变量的最小不可约分解, 并通过全局共轭得到一类平面三次多项式的光滑分类. 作者的定理还应用于讨论平面迭代根与嵌入流问题.

关键词: 平面三次多项式; 共轭; 光滑正规形; 最小不可约分解; 迭代根

Abstract

In this paper, we consider cubic polynomial normal forms on R2. By employing the theory of polynomial algebra to find the minimal irreducible decomposition of the corresponding varieties, we obtain smooth classifications of a kind of planar cubic polynomials via global conjugations. Our results are also applied to study the problems of iterative roots and embedding flows on R2.

Keywords: planar cubic polynomial; conjugacy; smooth normal form; minimal irreducible decomposition; iterative root

PDF (549KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

李林, 刘玲伶, 余志恒. 平面三次多项式的分类[J]. 数学物理学报, 2025, 45(2): 554-566

Li Lin, Liu Lingling, Yu Zhiheng. On Classification of Planar Cubic Polynomials[J]. Acta Mathematica Scientia, 2025, 45(2): 554-566

1 引言

在动力系统的研究中, 最重要的问题之一是在保留所有必要的动力学性质的情况下, 将给定系统化为简单的系统, 例如著名的 Hartman-Grobman 定理[3],[7],[27]. 为了简化系统, 一种经典的方法是寻找光滑的等价关系以避免动力学性质的丢失[30],[31],[34],[36]. 在通常情况下, 寻找等价关系下的微分同胚类比向量场类更加复杂[22],[31]. 尽管对于向量场的正规形理论已有丰富的结果[2],[9],[10],[12],[14],[15],[26],[28],[29],[39],[40],[20], 微分同胚的相应结果却少之又少[11],[17],[18],[31],[32]. 两者正规形理论的关系可参考文献 [4],[17],[38].

对于任意的整数 k1, 令 Hk 为所有系数为实数的 k 阶齐次多项式构成的 n 维向量空间. 考虑一个 Cr 映射F:RnRn, 其中 r1. 假设原点是 F 的一个不动点, 那么映射 F 在原点的一个邻域 ΩnRn 内可以写成

F(x)=Ax+F2(x)++Fr(x)+o(xr),
(1.1)

其中 AGL(n,R) 并且 FkHk 对所有 2kr 成立. 根据著名的 Poincaré-Dulac 定理, 任何具有固定线性部分的形式映射在坐标变换下都可以化为共振形式[3]. 进一步, 如果存在有限共振关系, 则微分同胚的胚芽在形式上是有限确定的. 如果 F 有无限多个共振关系, Ichikawa 定理指出 F 有无限个共振关系是有限确定的当且仅当 F 是 1-共振且它的非线性部分不属于某个无限余维集[21] (亦可见文献 [6]). 事实上, Takens[37] 已得到低维 (RR2) 有限确定同胚的完整分类. 对于 Rn, n3, 任等人在最近的工作中考虑所谓的强 1-共振微分同胚的分类[33]. 关于非强 1-共振和 0-共振情况的更多工作见参考文献 [18],[32].

值得一提的是, 上述关于微分同胚的结果是基于 Poincaré-Dulac 理论, 特别是共振条件的使用, 那么经典正规形理论的主要思想是通过变量替换简化系统 (1.1). 在文献 [11],[17] 中, 作者从代数角度得到经典正规形的进一步简化. 也就是说, 考虑一个近恒同坐标变换

x=φ(y)=y+φk(y),
(1.2)

其中 φkHk, 2kr. 那么转换后的映射 (1.1) 在原点的一个小邻域内变成

G(x)=Ay+F2(y)++Fk1(y)+Gk(y)+o(y2k1),
(1.3)

这里 Gk(y)=Fk(y)[φk(Ay)Aφk(y)]. 进一步, 我们定义线性算子 LkA:HkHk

LkAφk(x)=φk(Ax)Aφk(x),φkHk.
(1.4)

RkLkAHk 中的值域, 而 CkRk 的任意互补子空间, 则 Hk=RkCk.

为了求解 (1.4) 式, 作者[11]提出了一种计算 Rk 的互补子空间基的方法. 另外, 基于经典的 Poincaré-Dulac 定理, 在一维情况下有以下两个有趣的定理[11].

定理1.1 F(x)=x+arx2r+1+jr+1ajx2j+1 满足 ajR, ar0, r1xΩ1. 那么 F 简化至 N2r+1 阶的形式为 G(x)=x+arx2r+1+bx4r+1, 其中 bRxΩ1.

定理1.2F(x)=x+jrajxj 满足 ajR, ar0, r2xΩ1. 那么 F 简化至 Nr 阶的形式为 G(x)=x+arxr+bx2r1, 其中 bRxΩ1.

显然, 定理 1.1-1.2 给出一维非双曲不动点正规形问题的完整答案, 同时也是非 Poincaré 类型[31]. 然而, 由于混合类型变量的复杂计算, 高维甚至 2 维的情况仍然未知. 在本文中, 我们继续考虑维度为 2 的情形. 具体地说, 我们将研究以下平面映射 F:R2R2 满足 F(0)=0, 具有形式

F(z)=Az+F3(z)+o(z3),
(1.5)

其中 A=E (或者 E).

注意到上述关于系统 (1.5) 的正规形结果都只考虑不动点附近的小邻域,即 Ωn, n1 (同样可见定理 1.1-1.2). 其主要原因是保证 (1.2) 式中变换的可逆性, 这也意味着此等价关系是局部的. 在本文中, 我们将在第 2 节给出整体共轭的描述. 为了克服计算的复杂性, 我们采用多项式代数理论[5],[8],[16],[35]计算相应代数系统的 Gröbner 基以得到其极小不可约分解, 从而求出映射 F 系数之间的代数关系, 并给出其正规形的精确表达式. 作为应用, 我们的定理也将用于研究迭代根与嵌入流问题.

2 共轭

在给出本节的主要结果之前, 我们首先从等式 (1.3) 中得到以下事实

此变换不影响 F 中阶小于等于 k1 的项.

A=E, 则 Gk(y)=Fk(y). 若 A=E, 则 Gk(y)=Fk(y)[(1)i+j+1]φk(y), 这意味着 Gk(y)=Fk(y)i+j 是奇数, 且 Gk(y)=Fk(y)2φk(y)i+j 是偶数.

k=2, 那么 x=y+φ2(y)G(x)=Ay+G2(y)+o(y3), 其中 G2(y)=F2(y)[φ2(Ay)Aφ2(y)].

利用多项式代数理论, 我们得到平面系统 (1.5) 经典正规形理论的简化结果.

定理2.1 对于平面三次映射 (1.5) 式, 存在一个线性可逆映射 g:R2R2 可将其转化为一个 3 阶的正规形

φ1(xy)=A(xy)+(x3y3)+o((x,y)3).

如果

(I)A=E, 当且仅当以下一个条件满足

(i) h300;

(ii) h30=0,g30h21;

(iii) h30=0,g30=h21,g21h12;

(iv) h30=h21=h12=0,g2213g30g12=0,g2123g21g033g12h03=0,g21g129g30g033g21h03=0.

(II)A=E, 当且仅当以下一个条件满足

(i) g30=g221g123g21g03h21+g212h213g212,h30=g21h213g12,h12=3g03h21g12,h03=3g21g12g03g3129g203h213g212;

(ii) g21=g12=g03=0,h30=h12(9g30h03h212)27h203,h21=h2123h03;

(iii) g21=g12=g03=h12=h03=0,h21=3g30,g300;

(iv) g12=h21=0,g21=h12,g30=h12h033g03,h30=h2129g03, 其中 gkl=1k!l!k+1kxlyF13(x,y), hkl=1k!l!k+1kxlyF23(x,y)F3=(F13,F23)T, k,l=0,1,2,3.

令线性可逆映射

g(xy)=(a1x+a2yb1x+b2y),

其中 a1,a2,b1,b2R. 由于 a1b2a2b10, 我们有

g1:[xy][b2a1b2a2b1x+(a2a1b2a2b1)yb1a1b2a2b1x+a1a1b2a2b1y].

A=E, 通过计算 g1Fg 并且比较等式两边的系数, 我们可得

g1Fg=[x+x3y+y3]+O((x,y)3),

则该问题转化为求解下述半代数系统 SPS

P1:=a31a2h30a31b2g30+a21a2b1h21a21b1b2g21+a1a2b21h12a1b21b2g12+a2b31h03b31b2g03+a1b2a2b1=0,P2:=3a21a22h303a21a2b2g30+a21a2b2h21a21b22g21+2a1a22b1h212a1a2b1b2g21+2a1a2b1b2h122a1b1b22g12+a22b21h12a2b21b2g12+3a2b21b2h033b21b22g03=0,P3:=3a1a32h303a1a22b2g30+2a1a22b2h212a1a2b22g21+a1a2b22h12a1b32g12+a32b1h21a22b1b2g21+2a22b1b2h122a2b1b22g12+3a2b1b22h033b1b32g03=0,P4:=a42h30a32b2g30+a32b2h21a22b22g21+a22b22h12a2b32g12+a2b32h03b42g03=0,P5:=a41h30a31b1g30+a31b1h21a21b21g21+a21b21h12a1b31g12+a1b31h03b41g03=0,P6:=3a31a2h30+a31b2h213a21a2b1g30+2a21a2b1h21a21b1b2g21+2a21b1b2h122a1a2b21g21+a1a2b21h122a1b21b2g12+3a1b21b2h03a2b31g123b31b2g03=0,P7:=3a21a22h30+2a21a2b2h21+a21b22h123a1a22b1g30+a1a22b1h212a1a2b1b2g21+2a1a2b1b2h12a1b1b22g12+3a1b1b22h03a22b21g212a2b21b2g123b21b22g03=0,P8:=a1a32h30+a1a22b2h21+a1a2b22h12+a1b32h03a32b1g30a22b1b2g21a2b1b22g12b1b32g03a1b2+a2b1=0,P9:=a1b2a2b10.

为了解决上述含参数 uisvis 的半代数系统 SPS, 我们注意到 P5a1 的系数为 h30, 那么我们有以下两种情形

(Ia)h300,(Ib)h30=0.

在情况 (Ia) 中, 根据 Pi(a1), i=1,2,3,4,6,7,8P5(a1) (见文献 [16]) 可知, I:=P1,P2,,P8 的理想可由 I:=P4,P5 的理想生成. 因此, 通过求解半代数系统 SPS

P4:=a42h30a32b2g30+a32b2h21a22b22g21+a22b22h12a2b32g12+a2b32h03b42g03=0,P5:=a41h30a31b1g30+a31b1h21a21b21g21+a21b21h12a1b31g12+a1b31h03b41g03=0,P9:=a1b2a2b10,

我们得到定理中的相应结果.

在情况 (Ib) 中, 我们把 h30=0 代入半代数系统 SPS. 根据 P9 的表达式, 系数 b1b2 不会同时等于 0, 这意味着我们需要考虑以下三种可能性

(Ib1)b2=0,b10,g1(a1)=0,(Ib2)b1=0,b20,g2(a1)=0,(Ib3)b10,b20,g1(a1)=0,g2(a1)=0.

事实上, 情况 (Ib3) 的讨论与情况 (I) 类似, 那么我们同样可以得到定理中的条件 (iii)-(iv). 接下来我们考虑情况 (Ib1) 并将 b2=0 代入 SPS, 于是 h21=h12=0, 我们也得到一个更为简单的半代数系统 {\mathbf{SPS}"}

\begin{eqnarray*} Q_1&:=&b_1^2h_{03} - 1=0, \\ Q_2&:=&a_1^3g_{30} + a_1^2b_1g_{21} + a_1b_1^2g_{12} - a_1b_1^2h_{03} + b_1^3g_{03}=0, \\ Q_3&:=&3a_1^2g_{30} + 2a_1b_1g_{21} + b_1^2g_{12}=0, \\ Q_{4}&:=&3a_1g_{30} + b_1g_{21}=0, \\ Q_{5}&:=&a_2^2g_{30} - 1=0, \\ Q_{6}&:=&a_1b_2-a_2b_1\ne0. \end{eqnarray*}

利用文献 [35],[41] 中的方法解决半代数系统 {\mathbf{SPS}"}, 我们只需考虑以下代数系统

\widetilde{{\mathbf{SPS}}}:=\{Q_1=0,Q_2=0,Q_3=0,Q_4=0,Q_5=0,1-\kappa Q_{6}=0\},

其中 \kappa\in\mathbb{R}. 基于对 \mathbb{C}[a_i\,{\rm s},b_i\,{\rm s},u_i\,{\rm s},v_i\,{\rm s},\kappa] 上任意项顺序理想 \mathbf{J}:=\langle Q_1,\cdots,Q_{5},1-\kappa Q_{6}\rangle 的简化 Gröbner 基 G 的计算[5],[8], 我们得出 G 不同于 \{1\} 的结论. 进而, 利用消去定理[35]消除 G\kappa, a_1, a_2, b_1, 得到消去理想 \mathbf{J}_1. 然后, 计算 \mathbf{J}_1 对应变量的最小不可约分解, 再由 g 得到定理所示的 f\varphi 之间共轭的条件 (ii). 另外, 将条件 (ii) 代入 {\mathbf{SPS}"} 并求解 a_1, a_2, b_1, 我们有

a_1=-\frac{b_1g_{21}}{3g_{30}}, a_2=\pm\frac{1}{\sqrt{g_{30}}}, b_1=\pm\frac{1}{\sqrt{h_{03}}}.

因此, 情况 (Ib2) 的讨论与 (Ib1) 类似, 我们同样可得条件 (ii).

情形 A\!=\!-E 的证明是类似的, 我们省略细节部分. 各种情况的变换可见本文最后的附录.

推论2.1 若 g_{03} h_{30}>0, 平面三次映射 (1.5) 在 A=E 的条件下 C^\infty 共轭于定理 2.1 中的映射 \varphi_1: \mathbb{R}^2\to \mathbb{R}^2.

推论2.2 定义映射 f: \mathbb{R}^2\to \mathbb{R}^2

\begin{eqnarray*} f\! \left( \begin{array}{cc} x \\ y \end{array} \right) \!=\! \left(\begin{array}{cc} -x+(\frac{g_{30}^2}{3}-g_{30}g_{21}g_{12}+\frac{g_{12}}{3})x^3+g_{30}x^2y+xy^2+g_{21}y^3 \\ -y+\frac{g_{30}g_{12}}{3}x^3+g_{12}x^2y+3g_{21}g_{12}xy^2+(\frac{1}{3}-g_{30}g_{21}+3 g_{21}^2 g_{12} )y^3 \end{array}\right)\!+\!o(\|(x,y)\|^3), \end{eqnarray*}

其中 g_{30},g_{21},g_{12}\in \mathbb{R}, 那么 fC^\infty 共轭于定理 2.1 中的映射 \varphi_2: \mathbb{R}^2\to \mathbb{R}^2.

进一步, 我们可知

定理2.2 定义映射 f: \mathbb{R}^2\to \mathbb{R}^2

\begin{eqnarray*} f \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} x+x^3+\Sigma_{i,j=0}^{i+j=4}u_{ij}x^iy^j \\ y+y^3+\Sigma_{i,j=0}^{i+j=4}v_{ij}x^iy^j \end{array}\right)+o(\|x,y\|^4), \end{eqnarray*}

其中 u_{ij},v_{ij}\in \mathbb{R}, i,j=0,1,2,3,4. 那么存在一个可逆映射 h: \mathbb{R}^2\to \mathbb{R}^2,

\begin{eqnarray*} h \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} x+a_3x^2+a_4xy+a_5y^2 \\ y+b_3x^2+b_4xy+b_5y^2 \end{array}\right), \end{eqnarray*}

a_3,a_4,a_5, b_3,b_4,b_5\in \mathbb{R} 使得

\begin{eqnarray*} h^{-1}\circ f\circ h=\varphi_3 \end{eqnarray*}

当且仅当以下一个条件成立

\begin{eqnarray*} &&{\rm (i)} 2v_{31}+v_{13}=3v_{40}+2v_{22}=u_{13}-2v_{04}=2u_{22}+3u_{04}=u_{31}+4v_{04} \\ && =4u_{40}+v_{13}=3v_{13}^2-16v_{22}v_{04}=u_{04}v_{13}-8v_{04}^2=2u_{04}v_{22}-3v_{13}v_{04}=0, \\ &&{\rm (ii)} u_{04}=-\frac{2u_{22}}{3}\ne0, u_{40}=u_{31}=u_{13}=v_{40}=v_{31}=v_{22}=v_{13}=v_{04}=0, \\ &&{\rm (iii)} v_{40}=-\frac{2v_{22}}{3}\ne0, u_{40}=u_{31}=u_{22}=u_{13}=u_{04}=v_{31}=v_{13}=v_{04}=0, \end{eqnarray*}

其中映射 \varphi_3: \mathbb{R}^2\to \mathbb{R}^2 定义为

\begin{eqnarray*} \varphi_3 \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} x+x^3 \\ y+y^3 \end{array}\right)+o(\|(x,y)\|^4). \end{eqnarray*}

这里我们可知, 在条件 (i) 中,

a_3=-u_{40}, a_4=u_{13}, a_5=-\frac{u_{13}^2}{4u_{40}}, b_3=-\frac{2u_{40}^2}{u_{13}}, b_4=2u_{40}, b_5=-\frac{u_{13}}{2};

在条件 (ii) 中,

a_5=\frac{u_{04}}{2}, a_3=a_4=b_3=b_4=b_5=0;

在条件 (iii) 中,

b_3=\frac{v_{40}}{2}, a_3=a_4=a_5=b_4=b_5=0.

进一步, 从定理 2.2 条件 (i) 的结论可得, 若 u_{40},v_{04}\neq 0 我们有

u_{04}=-\frac{2v_{04}^2}{u_{40}}, v_{40}=-\frac{2u_{40}^2}{v_{04}}, v_{31}=2u_{40}, v_{13}=-4u_{40},
v_{22}=\frac{3u_{40}^2}{v_{04}}, u_{22}=\frac{3v_{04}^2}{u_{40}}, u_{13}=2v_{04}, u_{31}=-4v_{04}.

因此, 我们有下面的推论.

推论2.3 定义映射 f: \mathbb{R}^2\to \mathbb{R}^2

\begin{eqnarray*} f \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} x+x^3-\frac{2v^2}{u}y^4+2vxy^3+\frac{3v^2}{u}x^2y^2-4vx^3y+ux^4 \\ y+y^3+vy^4-4uxy^3+\frac{3u^2}{v}x^2y^2+2ux^3y-\frac{2u^2}{v}x^4 \end{array}\right)+o(\|(x,y)\|^4), \end{eqnarray*}

其中 u,v\in \mathbb{R} \backslash \{0\}, 那么 fC^\infty 共轭于定理 2.2 中的映射 \varphi_3: \mathbb{R}^2\to \mathbb{R}^2.

注2.1 定理 2.2 中通过可逆映射 h 共轭于 \varphi_3 的映射 f 不是唯一的. 事实上, 我们可以考虑映射 f_1: \mathbb{R}^2\to \mathbb{R}^2,

\begin{eqnarray*} f_1 \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} x+\Sigma_{i,j=0}^{i+j=4}u_{ij}x^iy^j \\ y+\Sigma_{i,j=0}^{i+j=4}v_{ij}x^iy^j \end{array}\right)+o(\|(x,y)\|^4) \end{eqnarray*}

\varphi_3 通过 h 建立共轭关系, 即 h^{-1}\circ f_1\circ h=\varphi_3. 通过比较两边的系数, 特别是 x^3 (或者 y^3) 的系数, 我们得到 1=0, 这是一个矛盾. 换句话说, f_1 必须包含一个阶数小于4的非线性项, 以保证其共轭关系.

注2.2 根据定理 2.2 可知映射 f_2: \mathbb{R}^2\to \mathbb{R}^2,

\begin{eqnarray*} f_2 \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} -x+x^3+\Sigma_{i,j=0}^{i+j=4}u_{ij}x^iy^j \\ -y+y^3+\Sigma_{i,j=0}^{i+j=4}v_{ij}x^iy^j \end{array}\right)+o(\|(x,y)\|^4) \end{eqnarray*}

与映射 \varphi_4: \mathbb{R}^2\to \mathbb{R}^2,

\begin{eqnarray*} \varphi_4 \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} -x+x^3 \\ -y+y^3 \end{array}\right)+o(\|(x,y)\|^4) \end{eqnarray*}

无法通过 h 共轭, 其中 u_{ij} s 和 v_{ij} s 不全为 0. 事实上, 比较等式 h^{-1}\circ f\circ h=\varphi_4 两边对应的系数, 我们可知 a_3=a_4=a_5=b_3=b_4=b_5=0, 这意味着所有 u_{ij} s 和 v_{ij} s 都等于 0.

3 一些应用

在本文的最后部分, 我们将应用共轭结果研究迭代根和嵌入流问题.

对于一个给定的映射 F:\mathbb{R}^2\to \mathbb{R}^2, 其 n 次 (n\in \mathbb{N}) 迭代根是指寻找一个映射 f:\mathbb{R}^2\to \mathbb{R}^2 满足下述函数方程

f^n=F,
(3.1)

其中 f^n 表示 fn 次迭代, 即 f^n(x,y)=f\circ f^{n-1}(x,y)f^0(x,y)\equiv (x,y) 对任意(x,y)\in \mathbb{R}^2 成立.

关于一维映射的迭代根, 特别是连续单调映射的迭代根, 已有大量的结果[23], [24]. 显然, 第二部分和第三部分中定义的映射 \hat{\varphi}: \mathbb{R}^2\to \mathbb{R}^2, \hat{\varphi}(x,y)=(x+x^3,y+y^3), x,y\in \mathbb{R}, 是关于变量 xy 严格递增的. 那么, 根据定理 15.7[23] (亦可见定理 11.2.2[24]) 中关于单变量函数的结果和变量之间的独立性, 映射 \hat{\varphi} 在整个平面上存在任意次的连续迭代根. 因此, 基于第三部分的全局共轭关系, 我们直接得到下列结果.

定理3.1 令映射 f: \mathbb{R}^2\to \mathbb{R}^2

\begin{eqnarray*} f \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} x+g_{30}x^3+g_{21}x^2y+g_{12}xy^2+g_{03}y^3 \\ y+h_{30}x^3+g_{30}x^2y+g_{21}xy^2+g_{12}y^3 \end{array}\right), \end{eqnarray*}

其中 g_{30},g_{21},g_{12},g_{03},h_{30}\in \mathbb{R}g_{03} h_{30}>0. 那么, f 存在任意次连续的迭代根.

定理3.2 令映射 f: \mathbb{R}^2\to \mathbb{R}^2

\begin{eqnarray*} f \left( \begin{array}{cc} x \\ y \end{array} \right) = \left(\begin{array}{cc} x+x^3-\frac{2v^2}{u}y^4+2vxy^3+\frac{3v^2}{u}x^2y^2-4vx^3y+ux^4 \\ y+y^3+vy^4-4uxy^3+\frac{3u^2}{v}x^2y^2+2ux^3y-\frac{2u^2}{v}x^4 \end{array}\right), \end{eqnarray*}

其中 u,v\in \mathbb{R}\backslash\{ 0\}. 那么, f 存在任意次的连续迭代根.

\hat{\varphi} 不同, 第二部分和第三部分中的映射 \tilde{\varphi}: \mathbb{R}^2\to \mathbb{R}^2, \tilde{\varphi}(x,y)=(-x+x^3,-y+y^3), x,y\in \mathbb{R}, 关于每个变量 xy 都是非单调的. 因此, 单调迭代根的相关理论就不再适用[23], [24]. 事实上, 对于一维三次多项式映射, 由于阶数是素数他们不存在任意次的多项式型迭代根. 那么有个问题就很自然产生: 平面三次多项式映射 \tilde{\varphi} 是否存在多项式型迭代根? 这个答案是肯定的, 因为下列映射\omega_1,\omega_2,\omega_3,\omega_4: \mathbb{R}^2\to \mathbb{R}^2,

\begin{aligned} \omega_{1}\binom{x}{y} & =\binom{y-y^{3}}{-x}, \omega_{2}\binom{x}{y}=\binom{-y+y^{3}}{x}, \\ \omega_{3}\binom{x}{y} & =\binom{y}{-x+x^{3}}, \omega_{4}\binom{x}{y}=\binom{-y}{x-x^{3}} \end{aligned}

都是 \tilde{\varphi} 的两次连续迭代根, 即 \omega_i^2=\tilde{\varphi} 对所有 i=1,2,3,4 成立. 尽管每个映射 \omega_i 都是 3 次, 但由于 \omega_i^3 的阶数大于 3, 他们在迭代下都无法保次, 这使得高维多项式迭代根问题变成更加有趣[13].

另一方面, 我们也对嵌入流问题, 即迭代根问题的推广形式感兴趣. 对于一个给定的原点附近邻域 \Omega^2 上的完整向量场 \nu, 我们说自治微分方程 \dot{x}=\nu(x) 确定一个流 \Phi_{\nu}, 即对于每一个 t\in \mathbb{R} 映射 \Phi_{\nu}(t,.) 都是 \Omega^2 上的一个微分同胚. 然后, 它的逆问题是: 给定一个微分同胚 f\in C^r(\Omega^2,\mathbb{R}^2), r\geq 3, 是否可以找到一个完全向量场 \nu 满足

\Phi_{\nu}(1,x)=f(x), x\in \Omega^2.
(3.2)

著名的 Sternberg 线性化定理表明, 如果 \nu 满足某些非共振条件, 等式 (3.2) 必然成立[22]. 更多关于微分同胚可嵌入流的条件请参考文献 [1],[19],[25].

{\cal A}=\{f\in C^{r}(\Omega^2,\mathbb{R}^2)|f\left( \begin{array}{cc} x \\ y \end{array} \right)=\left(\begin{array}{cc} x+x^3 \\ y+y^3 \end{array}\right)+o(\|(x,y)\|^3), r\geq 3\},
{\cal B}=\{f\in C^{r}(\Omega^2,\mathbb{R}^2)|f\left( \begin{array}{cc} x \\ y \end{array} \right)=\left(\begin{array}{cc} -x+x^3 \\ -y+y^3 \end{array}\right)+o(\|(x,y)\|^3), r\geq 3\}.

作为引理 2.11[27]的推广, 我们得到本节的最后一个定理.

定理3.3f 是以下系统

\dot{X}=\left(\begin{array}{cc}\dot{x} \\\dot{y}\end{array}\right)=(1+a_1x^2+a_2y^2+a_3xy)^{-1}\left(\begin{array}{cc}x^3 \\y^3\end{array}\right)
(3.3)

的时间-1 映射, 其中 X=(x,y)^T\in \Omega^2 并且 a_i\in \mathbb{R}, i=1,2,3. 那么 f\in {\cal A}.

f^tX 为系统 (3.3) 的流, 即

f^tX=\left(\begin{array}{cc}\sum_{i+j=1}^{\infty} a_{ij}(t)x^iy^j\\\sum_{i+j=1}^{\infty} b_{ij}(t)x^iy^j\end{array}\right),
(3.4)

其中

a_{10}(0)=1, a_{ij}(0)=0, (i,j)\not\in \{(1,0)\},
(3.5)
b_{01}(0)=1, b_{ij}(0)=0, (i,j)\not\in \{(0,1)\}
(3.6)

为系统的初始条件.

将 (3.4) 式带入 (3.3) 式, 我们可得

\begin{matrix}\label{eqasum} \sum_{i+j=1}^\infty a_{ij}'(t)x^iy^j&&=[1+a_1\big(\sum_{i+j=1}^\infty a_{ij}(t)x^iy^j\big)^2 +a_2(\sum_{i+j=1}^\infty b_{ij}(t)x^iy^j)^2\nonumber\\ && +a_3\sum_{i+j=1}^\infty a_{ij}(t)x^iy^j\sum_{i+j=1}^\infty b_{ij}(t)x^iy^j]^{-1}(\sum_{i+j=1}^\infty a_{ij}(t)x^iy^j)^3 \end{matrix}
(3.7)

\begin{matrix}\label{equbsum} \sum_{i+j=1}^\infty b_{ij}'(t)x^iy^j&&=[1+a_1(\sum_{i+j=1}^\infty a_{ij}(t)x^iy^j)^2 +a_2(\sum_{i+j=1}^\infty b_{ij}(t)x^iy^j)^2\nonumber\\ && +a_3\sum_{i+j=1}^\infty a_{ij}(t)x^iy^j\sum_{i+j=1}^\infty b_{ij}(t)x^iy^j]^{-1}(\sum_{i+j=1}^\infty b_{ij}(t)x^iy^j)^3. \end{matrix}
(3.8)

将上述两个方程 (3.7)-(3.8) 的两边分别对 x,y,xy,x^2y^2 项进行比较, 我们得到

\begin{eqnarray*} && a'_{10}(t)=a'_{01}(t)=a'_{11}(t)=a'_{20}(t)=a'_{02}(t) \\ &&=b'_{10}(t)=b'_{01}(t)=b'_{11}(t)=b'_{20}(t)=b'_{02}(t)=0. \end{eqnarray*}

再利用初始条件 (3.5)-(3.6), 我们有

a_{01}(t)=a_{11}(t)=a_{20}(t)=a_{02}(t)=b_{10}(t)=b_{11}(t)=b_{20}(t)=b_{02}(t)=0,
a_{10}(t)=b_{01}(t)=1.

进一步, 将上述两个方程 (3.7)-(3.8) 的两边对 xy^2,x^2y, x^3y^3 项再进行比较, 我们得到

\begin{eqnarray*} a'_{12}(t)=a'_{21}(t)=a'_{03}(t)=b'_{12}(t)=b'_{21}(t)=b'_{30}(t)=0, a'_{30}(t)=b'_{03}(t)=1. \end{eqnarray*}

根据初始条件 (3.5)-(3.6) 可知

a_{12}(t)=a_{21}(t)=a_{03}(t)=b_{12}(t)=b_{21}(t)=b_{30}(t)=0, a_{30}(t)=b_{03}(t)=1

因此,

fX=f^1X=\left(\begin{array}{cc}x+x^3 \\y+y^3\end{array}\right)+o(\|(x,y)\|^3)

并且 f\in {\cal A}.

类似的, 我们同样可得

推论3.1 令 f 为以下系统

\dot{X}=\left(\begin{array}{cc}\dot{x} \\\dot{y}\end{array}\right)=(-1+a_1x^2+a_2y^2+a_3xy)^{-1}\left(\begin{array}{cc}x^3 \\y^3\end{array}\right)

的时间-1 映射, 其中 X=(x,y)^T\in \Omega^2a_i\in \mathbb{R}, i=1,2,3. 那么 f\in {\cal B}.

附录

以下是定理 2.1 中第二种情形下的具体变换.

在情况 (Ii) 中,

a_1=\alpha b_1, a_2=\beta b_2,

其中 \alpha\beta 是方程

\begin{eqnarray*} h_{30}x^4+(h_{21}-g_{30})x^3+(h_{12}-g_{21})x^2+(h_{03}-g_{12})x-g_{03}=0 \end{eqnarray*}

两个不同的实根.

在情况 (Iii) 中,

b_2=0, a_1=-\frac{b_1g_{21}}{3g_{30}}, a_2=\pm\frac{1}{\sqrt{g_{30}}}, b_1=\pm\frac{1}{\sqrt{h_{03}}}

b_1=0, a_1=\pm\frac{1}{\sqrt{g_{30}}}, a_2=-\frac{b_2g_{21}}{3g_{30}}, b_2=\pm\frac{1}{\sqrt{h_{03}}}.

在情况 (Iiii) 中,

a_1=0, b_1\ne 0, a_2=\gamma_{1,2}b_2, a_1b_2-a_2b_1\ne0, g_{03}=0,

其中

\gamma_{1,2}=-{\frac {{g_{21}}-{h_{12}}\pm\sqrt {-4\,{g_{12}}\,{g_{30}}+4\,{h_{03}}\,{g_{30}}+g_{21}^{2}-2\,{g_{21}}\,{h_{12}}+4\,{g_{12}}\,{h_{21}}-4 \,{h_{03}}\,{h_{21}}+h_{12}^{2}}}{2\,{g_{30}}-2\,{h_{21}}}}

-4\,{g_{12}}\,{g_{30}}+4\,{h_{03}}\,{g_{30}}+g_{21}^{2}-2\,{g_{21}}\,{h_{12}}+4\,{g_{12}}\,{h_{21}}-4 \,{h_{03}}\,{h_{21}}+h_{12}^{2}\geq 0

a_2=0, b_2\ne 0, a_1=\gamma_{1,2}b_1, a_1b_2-a_2b_1\ne0, g_{03}=0

a_1=\gamma_{1,2}b_1, a_2=\gamma_{1,2}b_2, a_1b_2-a_2b_1\ne0, g_{03}=0

b_1=\rho_{1}a_1, b_2=\rho_{2}a_2, g_{03}\ne0, \left( -{h_{12}}+{g_{21}} \right) {g_{03}}-\frac{g_{12}^{2}}{3}+\frac{{2h_{03}}\,{g_{12}}}{3}-\frac{h_{03}^{2}}{3}=0,

其中 \rho_1\rho_2 是方程

\begin{eqnarray*} 3g_{03}^2x^3 + (3g_{12}g_{03} - 3g_{03}h_{03})x^2 + (g_{12}^2 - 2g_{12}h_{03} + h_{03}^2)x + 3g_{30}g_{03} - 3g_{03}h_{21}=0 \end{eqnarray*}

两个不同的实根, 或

b_1=\varrho_{1}a_1, b_2=\varrho_{2}a_2, g_{03}\ne0, \left( -{h_{12}}+{g_{21}} \right) {g_{03}}-\frac{g_{12}^{2}}{3}+\frac{{2h_{03}}\,{g_{12}}}{3}-\frac{h_{03}^{2}}{3}\ne0,

其中 \varrho_1\varrho_2 是方程

\begin{eqnarray*} g_{03}x^3 + (g_{12}-h_{03})x^2 + (g_{21}-h_{12})x + g_{30} - h_{21}=0 \end{eqnarray*}

两个不同的实根.

在情况 (Iiv) 中,

a_1=0, b_1\ne 0, a_2=\frac{h_{03}-g_{12}}{g_{21}-h_{12}}b_2, a_1b_2-a_2b_1\ne0, g_{03}=0

a_1=\frac{h_{03}-g_{12}}{g_{21}-h_{12}}b_1, b_2\ne 0, a_2=0, a_1b_2-a_2b_1\ne0, g_{03}=0

b_1=\xi_1a_1, b_2=\xi_2a_2, a_1b_2-a_2b_1\ne0, g_{03}\ne0,

其中 \xi_1\xi_2 是方程

\begin{eqnarray*} g_{03}x^2 + (g_{12} - h_{03})x+ g_{21} - h_{12} =0 \end{eqnarray*}

两个不同的实根.

在情况 (IIi) 中,

b_1=-\frac{a_1g_{12}g_{21} - 3a_1g_{03}h_{21} \pm \sqrt{a_1^2g_{21}^2g_{12}^2 - 6a_1^2g_{21}g_{12}g_03h_{21} + 4a_1^2g_{12}^3h_{21} + 9a_1^2g_{03}^2h_{21}^2}}{2g_{12}^2},
b_2=-\frac{a_2g_{21}g_{12} - 3a_2g_{03}h_{21} \pm \sqrt{a_2^2g_{21}^2g_{12}^2 - 6a_2^2g_{21}g_{12}g_{03}h_{21} + 4a_2^2g_{12}^3h_{21} + 9a_2^2g_{03}^2h_{21}^2}}{2g_{12}^2}

g_{21}^2g_{12}^2 - 6g_{21}g_{12}g_{03}h_{21} + 4g_{12}^3h_{21} + 9g_{03}^2h_{21}^2\leq 0

b_1=-\frac{a_1g_{12} \pm \sqrt{a_1^2g_{12}^2-3a_1^2g_{21}g_{03}}}{3g_{03}}, b_2=-\frac{a_2g_{12} \pm \sqrt{a_2^2g_{12}^2-3a_2^2g_{21}g_{03}}}{3g_{03}}

g_{12}^2-3g_{21}g_{03}\leq 0

b_1=-\frac{a_1g_{21}}{2g_{12}}, b_2=-\frac{a_2g_{21}}{2g_{12}}

g_{03}=0.

在情况 (IIii) 中,

a_1=0, a_2=\pm\frac{1}{\sqrt{g_{30}}}, b_1=\pm\frac{1}{\sqrt{h_{03}}}, b_2=\mp\frac{h_{12}}{3h_{03}\sqrt{g_{30}}}

a_2=0, a_1=\pm\frac{1}{\sqrt{g_{30}}}, b_2=\pm\frac{1}{\sqrt{h_{03}}}, b_1=\mp\frac{h_{12}}{3h_{03}\sqrt{g_{30}}}.

在情况 (IIiii) 中,

b_1= -\frac{a_1h_{30}}{2g_{30}}, b_2= -\frac{a_2h_{30}}{2g_{30}}.

在情况 (IIiv) 中,

b_1=\pm\frac{\sqrt{-3g_{03}h_{12}}a_1}{3g_{03}}, b_2=\frac{(3h_{03} \pm \sqrt{12g_{03}h_{12} + 9h_{03}^2})a_2}{6g_{03}}

g_{03}h_{12}\leq 0, 4g_{03}h_{12} + 3h_{03}^2\leq 0

b_1=\frac{(3h_{03} \pm \sqrt{12g_{03}h_{12} + 9h_{03}^2})a_1}{6g_{03}}, b_2=\frac{(3h_{03} \pm \sqrt{12g_{03}h_{12} + 9h_{03}^2})a_2}{6g_{03}}

4g_{03}h_{12} + 3h_{03}^2\leq 0.

参考文献

Arango J, Gómez A.

Diffeomorphisms as time one maps

Aequationes Math, 2002, 64: 304-314

[本文引用: 1]

Arcet B, Giné J, Romanovski V.

Linearizability of planar polynomial Hamiltonian systems

Nonlinear Anal-Real, 2022, 63: Article 103422

[本文引用: 1]

Arnold V I, Yu S I. Ordinary Differential Equations. Berlin: Springer-Verlag, 1988

[本文引用: 2]

Ashkenazi M, Chow S.

Normal forms near critical points for differential equations and maps

IEEE Trans Circuits Systems, 1988, 35: 850-862

[本文引用: 1]

Becker T, Weispfenning V. Gröbner Bases:A Computational Approach to Commutative Algebra. New York: Springer, 1993

[本文引用: 2]

Belitskii G R.

Equivalence and normal forms of germs of smooth mappings

Russ Math Surv, 1978, 33: 107-177

[本文引用: 1]

Bronstein I U, Kopanskii A Y. Smooth Invariant Manifolds and Normal Forms. Singapore: World Scientific, 1994

[本文引用: 1]

Buchberger B. Ein Algorithmus Zum Auffinden Der Basiselemente Des Restklassenringes Nach Einem Nulldimensional Polynomideal. Austria: Universität Innsbruck, 1965

[本文引用: 2]

Chen G.

Further reduction of normal forms for vector fields

Numerical Algorithms, 2001, 27: 1-33

[本文引用: 1]

Chen G, Dora J.

Further reductions of normal forms for dynamical systems

J Differ Equ, 2000, 166: 79-106

[本文引用: 1]

Chen G, Dora J.

Normal forms for differentiable maps near a fixed point

Numerical Algorithms, 1999, 22: 213-230

[本文引用: 4]

Chen G, Dora J.

An algorithm for computing a new normal form for dynamical systems

J Symbolic Computation, 2000, 29: 393-418

[本文引用: 1]

Chen X, Shi Y, Zhang W.

Planar quadratic degree-preserving maps and their iteration

Result Math, 2009, 55: 39-63

[本文引用: 1]

Chen G, Wang D, Yang J.

Unique normal forms for Hopf-zero vector fields

C R Acad Sci Paris, 2003, 336: 345-348

[本文引用: 1]

Cuong L, Doan T, Siegmund S.

A sternberg theorem for nonautonomous differential equations

J Dyn Diff Equ, 2019, 31: 1279-1299

[本文引用: 1]

Gelfand M, Kapranov M M, Zelevinsky A V. Discriminants, Resultants and Multidimensional Determinant. Boston: Birkhäuser, 1994

[本文引用: 2]

Gramchev T, Walcher S.

Normal forms of maps: formal and algebraic aspects

Acta Appl Math, 2005, 87: 123-146

[本文引用: 3]

Guo L, Ren Z.

Polynomial normal forms for some germs of nonstrongly 1-resonant diffeomorphisms

Int J Bifurcation and Chaos, 2015, 25: Article 1550174

[本文引用: 2]

Han M.

Conditions for a diffeomorphism to be embedded in a C^r flow

Acta Math Sinica, 1988, 4: 111-123

[本文引用: 1]

Han Z, Zhang Z, Zheng Y.

New normal for nonlinear system (II)-controllable form

Acta Math Sci, 1993, 13B(2): 139-146

[本文引用: 1]

Ichikawa F.

Finitely determined singularities of formal vector fields

Invent Math, 1982, 66: 199-214

[本文引用: 1]

Katok A, Hasselblatt B. Introduction to the Modern Theory of Dynamical Systems. New York: Cambridge University Press, 1995

[本文引用: 2]

Kuczma M. Functional Equations in a Single Variable. Warszawa: Polish Scientific Publishers, 1968

[本文引用: 3]

Kuczma M, Choczewski B, Ger R. Iterative Functional Equations. Cambridge: Cambridge University Press, 1990

[本文引用: 3]

Lam P F.

Embedding a differentiable homeomorphism in a flow subject to a regularity condition on the derivatives of the positive transition homeomorphisms

J Differ Equ, 1978, 30: 31-40

[本文引用: 1]

Li J, Kou L, Wang D, Zhang W.

Unique normal form and the associated coefficients for a class of three-dimensional nilpotent vector fields

Int J Bifurcation and Chaos, 2017, 27(14): 1750224

[本文引用: 1]

李伟固. 正规形理论及其应用. 北京: 科学出版社, 2000

[本文引用: 2]

Li W G. Normal Form Theory and Its Application. Beijing: Science Press, 2000

[本文引用: 2]

Li W, Llibre J, Wu H.

Polynomial and linearized normal forms for almost periodic differential systems

Discrete and Continuous Dynamical Systems, 2015, 36(1): 345-360

[本文引用: 1]

Liu L, Tang Y, Zhang W.

Versal unfolding of homogeneous cubic degenerate centers in strong monodromic family

J Differ Equ, 2021, 283: 136-162

[本文引用: 1]

Rayskin V.

Hölder linearization

J Differ Equ, 1998, 147: 271-284

[本文引用: 1]

Ren Z, Xia L, Yang J.

On classification of the Poincaré type maps on \mathbb{R}^3

Bull Sci Math, 2015, 139(5): 582-598

[本文引用: 4]

Ren Z, Yang J, Yu F.

On conjugating equivalence of 0-resonant diffeomorphisms on \mathbb{R}^3

Int J Bifurcation and Chaos, 2013, 23(6): Article 1350100

[本文引用: 2]

Ren Z, Yan S, Yang J.

Finite determinacy and polynomial normal forms for diffeomorphisms near a strongly 1-resonant fixed point

J Math Anal Appl, 2014, 413(1): 476-481

[本文引用: 1]

Rodrigues H M, Sola-Morales J.

Linearization of class C^1 for contractions on Banach spaces

J Differ Equ, 2004, 201: 351-382

[本文引用: 1]

Romanovski V, Shafer D. The Center and Cyclicity Problems:A Computational Algebra Approach. Boston: Birkhäuser, 2009

[本文引用: 3]

Strien S V.

Smooth linearization of hyperbolic fixed points without resonance conditions

J Differ Equ, 1990, 85(1): 66-90

[本文引用: 1]

Takens F.

Singularities of vector fields

Publ Math de lÍnstit des Hautes Etudes Sci, 1974, 43: 47-100

[本文引用: 1]

Takens F.

Normal forms for certain singularities of vector fields

Ann Inst Fourier, 1973, 23(2): 163-195

[本文引用: 1]

Turqui A, Dali D.

Normal forms of planar polynomial differential systems

Qual Theory Dyn Syst, 2019, 18: 11-33

[本文引用: 1]

Yang J.

Polynomial normal forms for vector fields on \mathbb{R}^3

Duke Math J, 2001, 106(1): 1-18

[本文引用: 1]

Yu Z, Yang L, Zhang W.

Discussion on polynomials having polynomial iterative roots

J Symbolic Comput, 2012, 47(10): 1154-1162

[本文引用: 1]

/