Processing math: 30%

数学物理学报, 2019, 39(6): 1352-1364 doi:

论文

多项式型迭代方程的多项式解

余志恒1, 龚小兵,2

Polynomial Solutions of the Polynomial-Like Iterative Equation

Yu Zhiheng1, Gong Xiaobing,2

通讯作者: 龚小兵, E-mail: xbgong@163.com

收稿日期: 2018-09-20  

基金资助: 国家自然科学基金.  11701476
中央高校基本科研业务费专项资金.  2682018CX63
四川省教育厅自然科学基金.  17ZA0217

Received: 2018-09-20  

Fund supported: the NSFC.  11701476
the Fundamental Research Funds for the Central Universities.  2682018CX63
the Key Project of Sichuan Provincial Department of Education.  17ZA0217

摘要

关于多项式型迭代方程的绝大多数结果都是在已知函数为单调函数情形下给出的.该文研究了多项式型迭代方程在已知函数为特殊的非单调函数-多项式函数情形下的多项式解.首先,在一维情形下,利用计算机代数系统SINGULAR分解代数簇的方法分别给出了二次和三次多项式型迭代方程有二次多项式解的充分必要条件,以及解的具体形式,并进一步给出计算多项式解的算法.最后利用一维情形的思想研究了二维情形下多项式型迭代方程的多项式解,在已知函数为平面二次齐次多项式映射时,得到了二次多项式型迭代方程有平面二次保次多项式解的几个充要条件.

关键词: 多项式型迭代方程 ; 首项系数问题 ; 多项式解 ; 代数簇 ; 最小不可约分解

Abstract

Most of known results for the polynomial-like iterative equation were given for monotone functions. In this paper, we discuss this equation for a polynomial function, which is non-monotonic. In one-dimensional case, we apply the method of computer algebra system SINGULAR decomposing algebraic varieties to find a sufficient and necessary condition for the polynomial-like iterative equation of orders 2 and 3 having quadratic polynomial solutions respectively and give quadratic polynomial solutions of both two equations. Then we give a procedure for computing polynomial solutions of the polynomial-like equation. In two-dimensional case, applying the idea of one-dimensional case we obtain several sufficient and necessary conditions for second order polynomial-like iterative equation having quadratic degree-preserving polynomial solutions when the given function is a two-dimensional homogeneous polynomial mapping of degree 2.

Keywords: Polynomial-like iterative equation ; Leading Coefficient Problem ; Polynomial solutions ; Algebraic variety ; Minimal irreducible decomposition

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

本文引用格式

余志恒, 龚小兵. 多项式型迭代方程的多项式解. 数学物理学报[J], 2019, 39(6): 1352-1364 doi:

Yu Zhiheng, Gong Xiaobing. Polynomial Solutions of the Polynomial-Like Iterative Equation. Acta Mathematica Scientia[J], 2019, 39(6): 1352-1364 doi:

1 引言

S是复数域C的一个非空子集.自映射f:SSi次迭代通过递推公式定义为fi(x)=f(fi1(x)), f0(x)=x,xS,i=1,2,.多项式型迭代方程

λ1f(x)+λ2f2(x)++λnfn(x)=F(x), xS
(1.1)

是一类重要的迭代方程.很多动力学问题,如迭代根问题[18, 35],不变曲线[18, 28],费根鲍姆方程[7, 26],微分同胚的横截同宿相交[10],动力系统中的正规形[2]和二次映射的动力学[16]均可转化为此方程.关于该方程的研究已有很多结果:当S为区间, F为线性函数时的结果见文献[8-9, 17, 23-25, 27, 34, 40, 49], F为非线性函数时的结果见文献[14, 29, 32, 38, 44-45, 47],在高维空间[13, 15, 19, 33]和多值函数[21, 39]情形下也有许多结果.

随着对迭代方程的深入研究,我们希望考虑当|λn|较大时方程(1.1)解的存在性和惟一性.这个问题被称为首项系数问题[43, 48].若λn0,则方程(1.1)可转化为如下形式

fn(x)+κn1fn1(x)++κ1f(x)=F(x),  xS.
(1.2)

2004年张伟年[46]利用Schröder变换f(x)=ϕ(cϕ1(x))将方程(1.2)转化为一个不含未知函数迭代的辅助方程,并利用Schauder不动点定理证明了在已知函数F扩张和非双曲条件下,该方程在不动点附近解的存在性.随后,在更多扩张和压缩条件下局部C1解被研究[6]. 2007年徐冰[37]利用逐段定义法构造出了无限多个递增连续解. 2013年,利用局部线性解逼近的方法全局解的存在性被证明[50].同样在2013年递减解也被构造出来[20].

关于首项系数问题的所有结果都是在已知函数为单调函数的情形下给出的. 2017年,在已知函数为PM函数(逐段单调函数)情形下,一类特殊的非单调解被讨论了并利用延拓的方法给出了多项式型迭代方程解的一般构造[22].

多项式函数作为一类特殊的非单调函数,在迭代根问题中已进行了讨论[41-42].受到文献[41-42]中思想的启发,在本文中我们在给定函数F为多项式函数情形下,在一维空间和二维空间中分别研究了方程(1.2).在一维空间中,利用计算机代数系统SINGULAR对代数簇进行分解得到了方程(1.2)在n=2n=3时有二次多项式解的充要条件以及两个方程的多项式解的具体形式,并进一步给出了计算方程(1.2)多项式解的可执行算法.在二维空间中,利用一维空间的思想得到了方程(1.2)在n=2和已知函数为平面二次齐次多项式时有平面二次保次多项式解的几个充要条件.

2 一维情形

设方程(1.2)中Ff均为多项式.并设degF=m, degf=l,则m=ln.为了方便给出符号计算的步骤,我们研究方程(1.2)的多项式解,其中F是一个次数为2n的多项式,即

F(x)=a2nx2n+a2n1x2n1++a2x2+a1x+a0, a2n0, aiC.
(2.1)

2.1 F为4次多项式时的多项式解

首先考虑如下方程

f2(x)+κ1f(x)=F(x), κ1C
(2.2)

的多项式解,其中

F(x)=a4x4+a3x3+a2x2+a1x+a0,a40.
(2.3)

定理2.1  设F形如(2.3), βa4的立方根.则方程(2.2)有多项式解的充要条件为

a0=16a24κ21β232a24κ1β2+16a22a248a2a23a4+32a2a24β+a4312a23a4β16a3a24κ116a3a2464a34;a1=a3(4a2a4a23)8a24.

并且其任一多项式解f形如

f(x)=βx2+a32β2x4κ1β44a2β3+2a3β2+a328β5.
(2.4)

  令f(x)=b2x2+b1x+b0.f(x)F(x)代入方程(2.2)比较两边的系数,得到多项式集F={f1,f2,f3,f4,f5},其中

f1:=b32a4=0;f2:=2b22b1a3=0;f3:=2b22b0+b1b2+κ1b2+b2b21a2=0;f4:=b21+κ1b1+2b2b1b0a1=0;f5:=b1b0+b2b20+κ1b0+b0a0=0.

APS1表示这些多项式方程构成的代数系统.计算APS1的约化Gröbner基[1, 3],然后选出不含b0b1的多项式(事实上,我们需要计算消元理想APS1C[ais,κ1,b2]的Gröbner基),即

g1:=b32a4;g2:=8b62a1+4b32a2a3a33;g3:=16b82κ2116b62κ1a332b82κ164b92a0+16b62a228b32a2a23+32b72a2+a43       12b42a2316b62a3.

最后,在计算机代数系统SINGULAR中利用minAssGTZ程序包(此方法是基于Gianni, Trager和Zacharias的工作[12])得到理想g1,g2,g3的极小相关素理想,我们可知,该理想即g1,g2,g3本身.计算过程表明结果无冗余.因此,由上述素理想定义出了消元理想代数簇的极小不可约分解如定理2.1所示(即方程(2.2)存在多项式解的条件).将条件代入代数系统APS1

b2=β,b1=a32β2,b0=4κ1β44a2β3+2a3β2+a328β5,

其中βa4的一个立方根.这样我们得到二次多项式解的表达式(2.4).

进一步,将得到的二次多项式f(x)代入方程(2.2)的左边得

f2(x)+κ1f(x)=β3x4+a3x3+a2x2+a3(4β3a2a23)8β6+16β8κ2132β8κ1+16β6a228β3a2a23+32β7a2+a4312β4a2316β6a3κ116β6a364β9.

注意到a4=β3.因此当

f(x)=βx2+a32β2x4κ1β44a2β3+2a3β2+a328β5

时方程(2.2)成立.证毕.

注2.1  在定理2.1中如果a0, a1,,a4, κ1R,则βa4的立方根,即β=a1/34.因此此定理在实数集R上也成立.

2.2 寻找多项式解的算法

为了给出计算方程(1.2)次数更高时的多项式解的算法,我们引入代数簇和不可约分解的基本定义.

K为一个代数闭域, AnK上的一个仿射n维空间.在环K[x1,,xn]上的多项式f可看成An上的一个k值函数在An中的点对f求值,即对每个xi选取K中的值.对K[x1,,xn]中的每个多项式集S,定义零点轨迹Z(S)An中点集,在其上S中的函数同时为零,即

Z(S)={xAnf(x)=0, fS}.

An中的子集V称为仿射代数集[1, 5, 11, 30]如果对某些SV=Z(S).非空仿射代数集V称为不可约的如果它不能表成两个适当的代数子集的并.不可约仿射代数集也称为仿射簇[5, 30].一般地,考虑环K[x1,,xn]上的多项式集F:={f1,,fm}.求解由集合F定义的多项式代数系统f1==fm=0可约化为寻找代数簇V(F).事实上,有三种行之有效的算法可用以寻找代数簇.其中一种重要的经典方法是建立在特征列方法上[30, 36].另外一种重要方法是对包含有n个变元的n次齐次多项式代数系统的结式法[5, 11],它是一种确定多项式是否有相同的非平凡解的工具.第三种方法是通过利用Gröbner基的方法来求解多项式方程[1, 3, 31],尤其是0维的情形.现在,随着Gröbner基的广泛应用,几乎所有的常用计算机代数系统(Mathematica, Maple, REDUCE, SINGULAR, Macaulay等)都能用来计算Gröbner基.下面,我们如定理2.1证明过程那样给出利用Gröbner基的方法来求解方程(1.2)的算法.

算法

第一步  计算fi, i=2,,n并比较方程(1.2)两边的系数得到m个关于变元ai, biκi (其中m=2n+1)的多项式方程(记为F:={f1,..,fm}).PS{f1=0,,fm=0}.

第二步  计算理想f1,,fm的约化Gröbner基,选出不含b0b1的多项式,记为PS.在计算机代数系统SINGULAR中利用程序包minAssGTZ求出PS对应的代数簇的极小不可约分解,并由此得到aiκi的最简代数关系,从而给出f(x)为方程(1.2)的解的条件.

第三步  将aiκi的关系代入PS,再返解PS,则bi可通过aiκi表示出来,进而最终得出方程解f(x)的表达式.

2.3 F为8次多项式时的多项式解

本节将利用上节给出算法求解三次多项式型迭代方程,即如下多项式型迭代方程

f3(x)+κ2f2(x)+κ1f(x)=F(x), κiC, i=1,2,
(2.5)

并给出方程(2.5)有二次多项式解的充分必要条件,其中

F(x)=a8x8+a7x7+a6x6+a5x5+a4x4+a3x3+a2x2+a1x+a0,a80.
(2.6)

定理2.2  设F是形如(2.6)的多项式, γa8的7次方根.则方程(2.5)有多项式解的充要条件为

a0=11048576a78(131072a2a7a58γ665536a4a6a7a48γ6+16384a4a37a38γ6+16384a36a7a38γ66144a26a37a28γ6512a6a57a8γ6+224a77γ6262144a4a6a58γ5+114688a4a27a48γ5+98304a36a48γ567584a26a27a38γ5+7168a6a47a28γ5+1568a67a8γ5+131072a4a7a58γ449152a26a7a48γ4+12288a6a37a38γ4+1792a57a28γ4+65536a26a58γ357344a6a27a48γ3+12544a47a38γ365536a6a7a58γ2+28672a37a48γ2262144a2a6a68+98304a2a27a58+65536a4a26a5832768a4a6a27a48+3072a4a47a3812288a46a48+6144a36a27a38+128a26a47a28160a6a67a8262144a6a68γ39a87+114688a27a58γ+131072a7a68);a1=a7(4096a2a58256a4a27a38+48a6a47a817a67)16384a68;a2=4096a58κ1γ+2048a4a6a48512a4a27a38512a36a38+192a26a27a28+16a6a47a87a674096a58;a3=a7(128a4a3820a6a27a8+7a47)256a48;a4=256a6a28γ5112a27a8γ5+512a38κ2γ3+192a26a2848a6a27a87a47512a38;a5=a7(24a6a87a27)32a28.
(2.7)

并且方程(2.5)的任一多项式解f形如

f(x)=γx2+a74γ2x+8a6γ74a7γ63a2732γ13.
(2.8)

  由上节给出的算法第一步令f(x)=b2x2+b1x+b0.将函数f(x)F(x)代入方程(2.5)并比较方程(2.5)两边系数得到多项式集F={f1,f2,,f9},其中

f1:=b72a8;f2:=4b62b1a7;f3:=6b52b21+2b52b1+4b62b0a6;f4:=4b42b31+6b42b21+12b52b1b0a5;f5:=b32b41+2b42b0+b21b32+κ2b32+6b42b1b0+6b32b31+b1b32+12b42b21b0+6b52b20a4;f6:=2b22b41+4b32b1b0+12b42b1b20+12b32b21b0+4b32b31b0+2κ2b22b1+2b31b22+2b22b21a3;f7:=4b22b21b0+6b32b20b1+κ1b2+2κ2b22b0+6b22b31b0+κ2b1b2+4b32b20+b2b21+4b42b30       +4b1b22b0+6b32b21b20+κ2b2b21+b2b41+b2b31a2;f8:=4b2b21b0+4b32b1b30+6b22b21b20+κ2b21+b31+2κ2b2b1b0+2b2b31b0+4b22b1b20       +κ1b1a1;f9:=b21b0+κ2b0+κ1b0+κ2b2b20+b1b0+κ2b1b0+b0+b2b20+2b22b30b1+b2b21b20       +b32b40+2b22b30+3b2b1b20a0.

APS2为这些多项式所构成的代数系统,由算法的第二步得到不含b0b1的新多项式代数系统APS2,即APS2:={g1=0,,g7=0},其中

g1:=a8b72;g2:=32b142a5+24b72a6a77a37;g3:=512b242κ2512b212a4+192b142a2648b72a6a27+256b192a67a47112b122a27;g4:=256b282a3+128b212a4a720b72a6a37+7a57;g5:=4096b362κ14096b352a2+2048b282a4a6512b212a4a27512b212a36+192b142a26a27       +16b72a6a477a67;g6:=16384b422a1+4096b352a2a7256b212a4a37+48b72a6a5717a77;g7:=1048576b492a0+262144b422a2a698304b352a2a27131072b412a2a765536b352a4a26       +32768b282a4a6a27+65536b342a4a6a7+262144b402a4a63072b212a4a4716384b272a4a37       114688b332a4a27131072b392a4a7+12288b282a466144b212a36a2716384b272a36a7       98304b332a36128b124a26a47+6144b202a26a37+67584b262a26a27+49152b322a26a7       65536b382a26+160b72a6a67+512b132a6a577168b192a6a4712288b252a6a37       +57344b312a6a27+65536b372a6a7+262144b432a6+39a87224b62a771568b122a67       1792b182a5712544b224a4728672b302a37114688b362a27131072b422a7.

然后利用计算机代数系统SINGULAR中的程序包minAssGTZ计算出APS2对应的代数簇的极小不可分解,我们可知该分解即V(g1,g2,,g7)本身.进一步,通过上述极小不可约分解我们得到aiκi的最简代数关系,从而得到f(x)为方程(2.5)解的条件,即定理2.2中所述.再由算法中的第三步将条件代入代数系统APS2得到

b2=γ,b1=a74γ2,b0=8a6γ74a7γ63a2732γ13,

其中γa8的7次方根.这样就得到多项式型迭代方程(2.5)的二次多项式解的表达式(2.8).

类似定理2.1的证明,将上面得到的f(x)的表达式代入方程(2.5)两边,注意到a8=γ7,则当

f(x)=γx2+a74γ2x+8a6γ74a7γ63a2732γ13

时方程(2.5)成立.证毕.

同样,定理2.2在实数集R上也成立.事实上,上面给出的算法也可以用来计算三次以上的多项式型迭代方程.

3 二维情形

本节我们将利用一维情形的思想研究多项式型迭代方程(1.2)在二维空间中存在多项式解的条件.由于平面多项式的迭代计算比较复杂,迭代后的次数无规律,因此在二维空间中要判断多项式型迭代方程(1.2)有多项式解要比一维情形复杂得多.为了简洁地给出符号计算的过程,考虑一类特殊的平面二次多项式映射,称为保次映射[4],这类多项式经过迭代后次数不超过2.我们考虑如下迭代方程

ϕ2ρ(x)+κ1ϕρ(x)=Fλ(x), κ1C, κ10,
(3.1)

其中

Fλ:[xy][a3x2+a4xy+a5y2b3y2+b4yx+b5x2],
(3.2)

λ=(λ1,λ2), λ1=(a3,a4,a5), λ2=(b3,b4,b5), ai, biC, i=3,4,5.我们研究方程(3.1)有如下二次多项式解

ϕρ:[xy][c0+c1x+c2y+c3x2+c4xy+c5y2d0+d1y+d2x+d3y2+d4yx+d5x2]
(3.3)

的条件,其中ρ=(ρ1,ρ2), ρ1=(c0,c1,c2,c3,c4,c5), ρ2=(d0,d1,d2,d3,d4,d5), ci, diC, i=0,,5.显然, ϕρ必须是一个保次映射,由文献[4]可知ϕ2ρ的2次迭代被分为11种类型

A11:={ρ:c3=c4=c5=d3=d4=d5=0},A12:={ρ:c1=c2=c3=c4=c5=0,d3=d5=0,d40},A13:={ρ:d1=d2=d3=d4=d5=0,c3=c5=0,c40},A21:={ρ:c3=c4=c5=d3=d4=0,d50},A22:={ρ:c1=c2=c3=c4=c5=0,d3=0,d4d50},A23:={ρ:c5=d3=0,c1=d2d4/d5,c2=d1d4/d5,c3=d4,c4=d24/d5,d4d50},A31:={ρ:c3=c4=d3=d4=d5=0,c50},A32:={ρ:d1=d2=d3=d4=d5=0,c3=0,c4c50},A33:={ρ:c3=d5=0,d1=c2c4/c5,d2=c1c4/c5,d3=c4,d4=c24/c5,c4c50},A41:={ρ:c1=c3d2/d5,c4=(c5d25+c33)/(c3d5),d1=c2d5/c3,d3=c5d5/c3,d4=(c5d25+c33)/c23,c5c33/d25,c3c5d50},A42:={ρ:c4=2c23/d5,c5=c33/d25,d3=c23/d5,d4=2c3,c5d50}.

因此,我们将根据上述11种类型给出方程(3.1)存在二次多项式解的条件.为了方便起见,如果\rho\in A_{ij},记\phi_{\rho}A_{ij}类型,其中(i, j)\in \{(m, n):m, n=1, 2, 3\}\cup \{(4, 1), (4, 2)\}.

定理3.1   (ⅰ)方程(3.1)没有A_{11}类型的多项式解.

(ⅱ)方程(3.1)有A_{12}类型的二次多项式解的充要条件为a_3=a_{4}=a_{5}=b_3=b_5=0b_4\neq0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}\phi_{\rho}(x, y)=(0, \frac{b_4}{\kappa_1}xy).\end{eqnarray*}

(ⅲ)方程(3.1)有A_{13}类型的二次多项式解的充要条件为a_3=a_{5}=b_3=b_4=b_5=0a_4\neq0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}\phi_{\rho}(x, y)=(\frac{a_4}{\kappa_1}xy, 0).\end{eqnarray*}

(ⅳ)方程(3.1)有A_{21}类型二次多项式解的充要条件为a_{4}=a_{5}=4\, b_{3}b_{5}+2\, a_{3}b_{4}-b_{4}^{2}=0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}&&\phi_{\rho}(x, y)=(0, \frac{b_5}{\kappa_1}x^2)~{\rm {\mbox {或}} }~\\&& \phi_{\rho}(x, y)=(-\kappa_1x, -\kappa_1y+\frac{b_5}{\kappa_1^2}x^2)~{\rm {\mbox {或}} }~\\&& \phi_{\rho}(x, y)=(-\kappa_1x, \frac{b_5}{\kappa_1^2+\kappa_1}x^2)~ {\rm {\mbox {或}} }~\\&&\phi_{\rho}(x, y)=(\frac{b_4}{2a_3}x+\frac{b_3}{a_3}y, -\frac{2a_3\kappa_1+b_4}{2a_3}y-\frac{b_4(2a_3\kappa_1+b_4)}{4a_3b_3}x+\frac{a_3^2}{b_3}x^2).\end{eqnarray*}

(ⅴ)方程(3.1)有A_{22}类型二次多项式解的充要条件为b_4b_5\neq 0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}&&\phi_{\rho}(x, y)=(0, \frac{b_4}{\kappa_1}yx+\frac{b_5}{\kappa_1}x^2)~{\rm {\mbox {或}} }~\\&&\phi_{\rho}(x, y)=(\frac{b_5(b_5+k_1)}{k_1^2b_4}, -\frac{b_5^2(b_5+k_1)}{k_1^2b_4^2}+\frac{k_1b_4}{b_5}yx+k_1x^2), ~\kappa_1=-1, ~\forall k_1\in{\Bbb C}\backslash\{0\}.\end{eqnarray*}

(ⅵ)方程(3.1)有A_{23}类型二次多项式解的充要条件为a_{3}+b_{4}=a_3^2+a_4b_5=0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}&&\phi_{\rho}(x, y)=(\frac{a_3}{\kappa_1}x^2-\frac{a_3^2}{\kappa_1b_5}xy, -\frac{a_3}{\kappa_1}yx+\frac{b_5}{\kappa_1}x^2)~{\rm {\mbox {或}} }~\\&&\phi_{\rho}(x, y)=(\frac{k_2a_3}{b_5}-a_3x^2+\frac{a_3^2}{b_5}xy, k_2+a_3yx-b_5x^2), \kappa_1=-1, ~\forall k_2\in{\Bbb C}~{\rm {\mbox {或}} }~\\&&\phi_{\rho}(x, y)=(-\frac{k_3a_4}{a_3}x^2-\frac{k_3a_4^2}{a_3^2}xy, \frac{b_5(b_5+k_3)}{k_3^2a_4}+\frac{k_3a_4}{a_3}yx+k_3x^2), \kappa_1=-1, ~\forall k_3\in{\Bbb C}\backslash\{0\}.\end{eqnarray*}

(ⅶ)方程(3.1)有A_{31}类型的二次多项式解的充要条件为b_{4}=b_{5}=4\, a_{5}a_{3}+2\, a_{4}b_{3}-a_{4}^{2}=0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}&&\phi_{\rho}(x, y)=(\frac{a_5}{\kappa_1}y^2, 0)~{\rm {\mbox {或}} }~\\&& \phi_{\rho}(x, y)=(-\kappa_1x+\frac{a_5}{\kappa_1^2}y^2, -\kappa_1y)~{\rm {\mbox {或}} }~\\&& \phi_{\rho}(x, y)=(\frac{a_5}{\kappa_1^2+\kappa_1}y^2, -\kappa_1y)~{\rm {\mbox {或}} }~\\&&\phi_{\rho}(x, y)=(-\frac{2b_3\kappa_1+a_4}{2b_3}x-\frac{a_4(2b_3\kappa_1+a_4)}{4a_3b_3}y+\frac{b_3^2}{a_3}y^2, \frac{a_4}{2b_3}y+\frac{a_3}{b_3}x).\end{eqnarray*}

(ⅷ)方程(3.1)有A_{32}类型二次多项式解的充要条件为a_4a_5\neq0,并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}&&\phi_{\rho}(x, y)=(\frac{a_4}{\kappa_1}xy+\frac{a_5}{\kappa_1}y^2, 0)~{\rm {\mbox {或}} }~ \\&&\phi_{\rho}(x, y)=(-\frac{a_5^2(a_5+k_4)}{k_4^2a_4^2}+\frac{k_4a_4}{a_5}xy+k_4y^2, \frac{a_5(a_5+k_4)}{k_4^2a_4}), \kappa_1=-1, ~\forall k_4\in{\Bbb C}\backslash\{0\}.\end{eqnarray*}

(ⅸ)方程(3.1)有A_{33}类型的二次多项式解的充要条件为a_{4}+b_{3}=a_4^2+a_5b_4=0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}&&\phi_{\rho}(x, y)=(\frac{a_4}{\kappa_1}xy+\frac{a_5}{\kappa_1}y^2, -\frac{a_4}{\kappa_1}x^2-\frac{a_4^2}{\kappa_1a_5}xy)~{\rm {\mbox {或}} }~\\&&\phi_{\rho}(x, y)=(k_4-a_4xy-a_5y^2, \frac{k_4b_4}{a_4}+a_4y^2+\frac{a_4^2}{a_5}yx), \kappa_1\!=\!-1, ~\forall k_4\in{\Bbb C}~{\rm {\mbox {或}} }~ \\ && \phi_{\rho}(x, y)=(\frac{k_5a_4}{a_5}xy+k_5y^2-\frac{a_5(a_5+k_5)}{k_5^2b_4}, -\frac{k_5a_4}{a_5}y^2-\frac{k_5a_4^2}{a_5^2}yx), \kappa_1=-1, ~\forall k_5\in{\Bbb C}\backslash\{0\}. \end{eqnarray*}

(ⅹ)方程(3.1)有A_{41}类型的二次多项式解的充要条件为a_4b_5-a_3b_4=a_5b_4-a_4b_3=a_5b_5-a_3b_3=b_3^2+b_3a_4+a_5a_3=b_4b_3+b_4a_4+a_4a_3=b_5b_3+b_4a_3+a_3^2=0a_3a_5b_5\ne0.并且方程(3.1)任一二次多项式解\phi_{\rho}形如

\begin{eqnarray*}&&\phi_{\rho}(x, y)=(\frac{a_3}{\kappa_1}x^2-\frac{a_3^3+a_5b_5^2}{\kappa_1a_3b_5}xy+\frac{a_5}{\kappa_1}y^2, \frac{a_5b_5}{\kappa_1a_3}y^2-\frac{a_3^3+a_5b_5^2}{\kappa_1a_3^2}yx+\frac{b_5}{\kappa_1}x^2)~{\rm {\mbox {或}} }~\\ && \phi_{\rho}(x, y)=( -a_3x^2+\frac{a_3^3+a_5b_5^2}{a_3b_5}xy-a_5y^2+k_6, -\frac{a_5b_5}{a_3}y^2+\frac{a_3^3+a_5b_5^2}{a_3^2}yx-b_5x^2+\frac{k_6b_3}{a_5}), \\&&~~~~~~~~~~~~\quad \kappa_1=-1, ~\forall k_6\in{\Bbb C}\backslash\{0\}. \end{eqnarray*}

(ⅹⅰ)方程(3.1)有A_{42}类型的二次多项式解的充要条件为a_4b_4-4a_5b_5=b_4^2-2a_3b_4+2a_4b_5-4b_3b_5=a_4^2-2a_4b_3+2a_5b_4-4a_3a_5=0.

  令F_{\lambda}\phi_\rho分别为形如(3.2)和(3.3)式的平面二次多项式.计算2次迭代\phi_{\rho}^2,将\phi_\rho^{2}(x), \kappa_{1}\phi_\rho(x), F_{\lambda}(x)代入方程(3.1)并比较\phi_{\rho}^2+\kappa_{1}\phi_\rhoF_{\lambda}对应的系数,得到多项式集合PS= \{g_1, g_2, \cdots, g_{30}\},其中

\begin{eqnarray*}g_1&:=&c_{{4}}c_{{3}}d_{{5}}+c_{{5}}d_{{5}}^{2}+c_{{3}}^{3}=0, \\g_2&:=&2c_{{3}}^{2}c_{{4}}+c_{{4}}^{2}d_{{5}}+c_{{4}}c_{{3}}d_{{4}}+2c_{{5}}d_{{4}}d_{{5}}=0, \\g_3&:=&2c_{{3}}^{2}c_{{5}}+c_{{3}}c_{{4}}^{2}+c_{{4}}^{2}d_{{4}}+c_{{5}}d_{{4}}^{2}+c_{{4}}c_{{3}}d_{{3}}+c_{{4}}c_{{5}}d_{{5}}+2c_{{5}}d_{{3}}d_{{5}}=0, \\&&\vdots, \\g_{10}&:=&2c_0c_3^2+c_0c_4d_5+c_1^2c_3+c_1c_4d_2+c_3c_4d_0+2c_5d_0d_5+c_5d_2^2+c_1c_3+c_2d_5\nonumber\\& & {}+c_3\kappa_1-a_3=0, \\g_{11}&:=&2c_0c_3c_4+c_0c_4d_4+2c_1c_2c_3+c_1c_4d_1+c_2c_4d_2+c_4^2d_0+2c_5d_0d_4+2c_5d_1d_2\nonumber\\& & {}+c_1c_4+c_2d_4+c_4\kappa_1-a_4=0, \\g_{12}&:=&2c_0c_3c_5+c_0c_4d_3+c_2^2c_3+c_2c_4d_1+c_4c_5d_0+2c_5d_0d_3+c_5d_1^2+c_1c_5+c_2d_3\nonumber\\& & {}+c_5\kappa_1-a_5=0, \\&&\vdots, \\g_{25}&:=&2c_0c_3d_5+c_0d_4d_5+c_1^2d_5+c_1d_2d_4+c_3d_0d_4+2d_0d_3d_5+d_2^2d_3+c_3d_2+d_1d_5\nonumber\\& & {}+d_5\kappa_1-b_5=0, \\g_{26}&:=&2c_0c_4d_5+c_0d_4^2+2c_1c_2d_5+c_1d_1d_4+c_2d_2d_4+c_4d_0d_4+2d_0d_3d_4+2d_1d_2d_3\nonumber\\& & {}+c_4d_2+d_1d_4+d_4\kappa_1-b_4=0, \\g_{27}&:=&2c_0c_5d_5+c_0d_3d_4+c_2^2d_5+c_2d_1d_4+c_5d_0d_4+2d_0d_3^2+d_1^2d_3+c_5d_2+d_1d_3\nonumber\\& & {}+d_3\kappa_1-b_3=0, \\g_{28}&:=&2c_0c_1d_5+c_0d_2d_4+c_1d_0d_4+2d_0d_2d_3+c_1d_2+d_1d_2+d_2\kappa_1=0, \\g_{29}&:=&2c_0c_2d_5+c_0d_1d_4+c_2d_0d_4+2d_0d_1d_3+c_2d_2+d_1^2+d_1\kappa_1=0, \\g_{30}&:=&c_0^2d_5+c_0d_0d_4+d_0^2d_3+c_0d_2+d_0d_1+d_0\kappa_1+d_0=0.\end{eqnarray*}

PS中的每一个多项式在有理数域{\bf Q}中都不可约.显然,由于二次多项式不可能由线性函数迭代给出,因此方程(3.1)无A_{11}类型的二次多项式解.注意到A_{12}A_{13}, A_{21}A_{31}, A_{22}A_{32}, A_{23}A_{33}关于变量c_id_i具有相同的对称性,因此只考虑A_{13}, A_{31}, A_{32}, A_{33}类型.下面我们详细讨论A_{31}类型,其它类型可类似讨论.

将保次条件A_{31}代入PS得半代数系统

PS_1:=\{{\mathfrak g}_1=0, {\mathfrak g}_2=0, \cdots , {\mathfrak g}_{10}=0, ~b_4=0, ~b_5=0, ~c_5\ne0\},

其中

\begin{eqnarray*} {\mathfrak g}_1&:=&c_5d_2^2-a_3, \\ {\mathfrak g}_2&:=&2c_5d_1d_2-a_4, \\ {\mathfrak g}_3&:=&c_5d_1^2+c_1c_5+c_5\kappa_1-a_5, \\ {\mathfrak g}_{4}&:=&2c_5d_0d_2+c_1^2+c_1\kappa_1+c_2d_2, \\ {\mathfrak g}_{5}&:=&2c_5d_0d_1+c_1c_2+c_2d_1+c_2\kappa_1, \\ {\mathfrak g}_{6}&:=&c_5d_0^2+c_0c_1+c_0\kappa_1+c_2d_0+c_0, \\ {\mathfrak g}_{7}&:=&c_5d_2-b_3, \\ {\mathfrak g}_{8}&:=&c_1d_2+d_1d_2+d_2\kappa_1, \\ {\mathfrak g}_{9}&:=&c_2d_2+d_1^2+d_1\kappa_1, \\ {\mathfrak g}_{10}&:=&c_0d_2+d_0d_1+d_0\kappa_1+d_0. \end{eqnarray*}

由于半代数系统的消元十分困难,注意到c_5\ne 0当且仅当存在w\in{\Bbb C}使得1-wc_5=0,因此,我们考虑环{\Bbb C}[a_i\, {\rm s}, b_i\, {\rm s}, c_i\, {\rm s}, d_i\, {\rm s}, \kappa_1, w]中的代数系统

\widehat{PS_1}:=\{{\mathfrak g}_1=0, {\mathfrak g}_2=0, \cdots , {\mathfrak g}_{10}=0, ~b_4=0, ~b_5=0, ~1-wc_5=0\}.

这样问题就转化为计算代数簇V({\bf I}),其中{\bf I}:=\langle {\mathfrak g}_1, {\mathfrak g}_2, \cdots, {\mathfrak g}_{10}, ~b_4, ~b_5, ~1-wc_5\rangle.在计算机代数系统SINGULAR中,我们采用分次字典序作为项序且变元序为

c_i\succ d_i\succ w \succ \kappa_1 \succ a_i\, {\rm s}\succ b_i\, {\rm s},

运用命令eliminate计算出代数系统\widehat{PS_1}的消元理想{\bf I}\cap {\Bbb C}[a_i\, s, b_i\, s, \kappa_1]的Gröbner基,记为{\bf I}'.并进一步消去变量c_i (i=0, 1, 2, 5), d_j (j=0, 1, 2)w.最后利用计算机代数系统SINGULAR中的程序包minAssGTZ找出消元理想的极小相关素理想,得到

\Lambda:=\{\lambda:~b_{4}=b_{5}=0, ~4\, a_{5}a_{3}+2\, a_{4}b_{3}-a_{4}^{2}=0\},

即方程(3.1)的A_{31}类型多项式解存在的条件.反之,我们将条件\lambda\in \Lambda代入PS_1

\begin{eqnarray*} &&{\rm (i)}\quad c_0=c_1=c_2=c_3=c_4=d_0=d_1=d_2=d_3=d_4=d_5=0, ~c_5=\frac{a_5}{\kappa_1}, \\ &&{\rm (ii)}\quad c_0=c_2=c_3=c_4=d_0=d_2=d_3=d_4=d_5=0, ~c_1=d_1=-\kappa_1, ~c_5=\frac{a_5}{\kappa^2}, \\ &&{\rm (iii)}\quad c_0=c_1=c_2=c_3=c_4=d_0=d_2=d_3=d_4=d_5=0, ~c_5=\frac{a_5}{\kappa_1^2+\kappa_1}, ~d_1=-\kappa_1, \\ &&{\rm (iv)}\quad c_0=c_3=c_4=d_0=d_3=d_4=d_5=0, ~c_1=-\frac{2b_3\kappa_1+a_4}{2b_3}, ~c_2=-\frac{a_4(2b_3\kappa_1+a_4)}{4a_3b_3}, \\ &&~~~~~\quad c_5=\frac{b_3^2}{a_3}, ~d_1=\frac{a_4}{2b_3}, ~d_2=\frac{a_3}{b_3}, \end{eqnarray*}

这样就给出方程(3.1)的A_{31}类型的四个解的表达式.进一步类似于定理2.1和2.2的证明,将上述得到的四个解分别代入方程(3.1)的左边结合条件\Lambda,得到当\phi_{\rho}为如下形式

\begin{eqnarray*} &&\phi_{\rho}(x, y)=(\frac{a_5}{\kappa_1}y^2, 0) ~{\rm {\mbox{或}} }\\&& \phi_{\rho}(x, y)=(-\kappa_1x+\frac{a_5}{\kappa_1^2}y^2, -\kappa_1y)~{\rm {\mbox{或}} }~\\&& \phi_{\rho}(x, y)=(\frac{a_5}{\kappa_1^2+\kappa_1}y^2, -\kappa_1y)~{\rm {\mbox{或}} }~\\&&\phi_{\rho}(x, y)=(-\frac{2b_3\kappa_1+a_4}{2b_3}x-\frac{a_4(2b_3\kappa_1+a_4)}{4a_3b_3}y+\frac{b_3^2}{a_3}y^2, \frac{a_4}{2b_3}y+\frac{a_3}{b_3}x)\end{eqnarray*}

时方程(3.1)成立.

利用A_{21}A_{31}的对称性得到集合\Lambda_3的表达式.进一步得到方程(3.1)的A_{21}类型的四个解.方程(3.1)的A_{41}A_{42}类型与A_{31}类型的多项式解存在的条件类似可得.证毕.

注3.1  在定理3.1中,虽然给出了方程(3.1)有A_{42}类型的二次多项式解的充要条件,但并没有给出这种情形下解的表达式,其原因是计算解的表达式需要计算机有较大的内存,而我们的民用计算机的内存不足以计算出解的具体表达式.

参考文献

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

[本文引用: 3]

Bessis D , Marmi S , Turchetti G .

On the singularities of divergent majorant series arising from normal form theory

Rend Math Ser, 1989, 9: 645- 659

[本文引用: 1]

Buchberger B. Ein Algorithmus zum Auffinden der Basiselemente des Restklassenringes nach einem Nulldimensional Polynomideal[D]. Austria: Universität Innsbruck, 1965

[本文引用: 2]

Chen X , Shi Y , Zhang W .

Planar quadratic degree-preserving maps and their iteration

Results Math, 2009, 55: 39- 63

DOI:10.1007/s00025-009-0389-6      [本文引用: 2]

Chen X , Zhang W .

Decomposition of algebraic sets and applications to weak centers of cubic systems

J Comput Appl Math, 2009, 232: 565- 581

DOI:10.1016/j.cam.2009.06.029      [本文引用: 3]

Chen J , Zhang W .

Leading coefficient problem for polynomial-like iterative equations

J Math Anal Appl, 2009, 349: 413- 419

DOI:10.1016/j.jmaa.2008.09.015      [本文引用: 1]

Collet P , Eckmann J P , Lanford O E .

Universal properties of maps on an interval

Comm Math Phys, 1980, 76: 211- 254

DOI:10.1007/BF02193555      [本文引用: 1]

Draga S .

A note on the polynomial-like iterative equations order

Comment Math, 2016, 56: 243- 249

URL     [本文引用: 1]

Draga S , Morawiec J .

Reducing the polynomial-like iterative equations order and a generalized Zoltán Boros' problem

Aequationes Math, 2016, 90: 935- 950

DOI:10.1007/s00010-016-0420-4      [本文引用: 1]

Fontich E .

Transversal homoclinic points of a class of conservative diffeomorphisms

J Diffe Equ, 1990, 87: 1- 27

DOI:10.1016/0022-0396(90)90012-E      [本文引用: 1]

Gelfand I M, Kapranov M M, Zelevinsky A V. Discriminants, Resultantsand Multidimensional Determinants. Boston: Birkhäuser, 1994

[本文引用: 2]

Gianni P , Trager B , Zacharias G .

Gröbner basis and primary decomposition of polynomials

J Symbolic Comput, 1988, 6: 146- 167

[本文引用: 1]

Gong X .

Convex solutions of the polynomial-like iterative equation on open set

Bull Korean Math Soc, 2014, 51: 641- 651

DOI:10.4134/BKMS.2014.51.3.641      [本文引用: 1]

Gong X , Zhang P .

Convex solutions of the polynomial-like iterative equation with variable coefficents

J Appl Anal Comput, 2016, 6: 306- 321

URL     [本文引用: 1]

Gong X , Zhang W .

Convex solutions of the polynomial-like iterative equation in Banach spaces

Publ Math Debrecen, 2013, 82: 341- 358

DOI:10.5486/PMD.2013.5305      [本文引用: 1]

Greenfield S J , Nussbaum R D .

Dynamics of a quadratic map in two complex variables

J Diffe Equ, 2001, 169: 57- 141

DOI:10.1006/jdeq.2000.3895      [本文引用: 1]

Jarczyk W .

On an equation of linear iteration

Aequationes Math, 1996, 51: 303- 310

DOI:10.1007/BF01833285      [本文引用: 1]

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

[本文引用: 2]

Kulczycki M , Tabor J .

Iterative functional equations in the class of Lipschitz functions

Aequationes Math, 2002, 64: 24- 33

DOI:10.1007/s00010-002-8028-2      [本文引用: 1]

Li L , Zhang W .

Continuously decreasing solutions for polynomial-like iterative equations

Sci China Ser, 2013, 56A: 1051- 1058

URL     [本文引用: 1]

Liu J, Duan H, Yu Z. On the blended soltuions of polynomial-like iterative equation with multivalued functions. Discrete Dyn Nat Soc, 2018, Art ID: 9862794

[本文引用: 1]

Liu L , Gong X .

The polynomial-like iterative equation for PM functions

Sci China Math, 2017, 60: 1503- 1514

DOI:10.1007/s11425-016-0319-x      [本文引用: 1]

Matkowski J , Zhang W .

Method of characteristic for functional equations in polynomial form

Acta Math Sinica, 1997, 13: 421- 432

DOI:10.1007/BF02560023      [本文引用: 1]

Matkowski J , Zhang W .

Characteristic analysis for a polynomial-like iterative equation

Chinese Sci Bull, 1998, 43: 192- 196

DOI:10.1007/BF02898908     

Matkowski J , Zhang W .

On linear dependence of iterates

J Appl Anal, 2000, 6: 149- 157

URL     [本文引用: 1]

McCarthy P J .

The general exact bijective continuous solution of Feigenbaum's functional equation

Comm Math Phys, 1983, 91: 431- 443

DOI:10.1007/BF01208784      [本文引用: 1]

Nabeya S .

On the functional equation f(p+qx+rf(x))=a+bx+cf(x)

Aequationes Math, 1974, 11: 199- 211

DOI:10.1007/BF01834919      [本文引用: 1]

Ng C T , Zhang W .

Invariant curves for planar mappings

J Differ Equ Appl, 1997, 3: 147- 168

DOI:10.1080/10236199708808092      [本文引用: 1]

Ng C T , Zhao H .

Periodic and continuous soltuions of a polynomial-like iterative equation

Aequationes Math, 2017, 91: 185- 200

DOI:10.1007/s00010-016-0456-5      [本文引用: 1]

Ritt J F. Differential Algebra. New York: Amer Math Soc, 1950

[本文引用: 3]

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

[本文引用: 1]

司建国.

迭代方程\sum\limits_{i = 1}^n {{\lambda _i}f^i(x) = F(x)} 局部解析解的存在性

数学学报, 1994, 37: 590- 600

[本文引用: 1]

Si J G .

The existence of local analytic solutions of the iterated equations \sum\limits_{i = 1}^n {{\lambda _i}f^i(x) = F(x)}

Acta Math Sin (Chinese Series), 1994, 37: 590- 600

[本文引用: 1]

Tabor J , Żołdak M .

Iterative equations in Banach spaces

J Math Anal Appl, 2004, 299: 651- 662

DOI:10.1016/j.jmaa.2004.06.011      [本文引用: 1]

Tabor J , Tabor J .

On a linear iterative equation

Results Math, 1995, 27: 412- 421

DOI:10.1007/BF03322847      [本文引用: 1]

Targonski G. Topics in Iteration Theory. Göttingen: Vandenhoeck and Ruprecht, 1981

[本文引用: 1]

Wu W T .

On zeros of algebraic equations-An appliation of Ritt's principle

Chin Sci Bul, 1986, 31: 1- 5

DOI:10.1360/csb1986-31-1-1      [本文引用: 1]

Xu B , Zhang W .

Construction of continuous solutions and stability for the polynomial-like iterative equation

J Math Anal Appl, 2007, 325: 483- 497

URL     [本文引用: 1]

Xu B , Zhang W .

Decreasing solutions and convex solutions of the polynomial-like iterative equation

J Math Anal Appl, 2007, 329: 1160- 1170

URL     [本文引用: 1]

Xu B , Nikodem K , Zhang W .

On a multivalued iterative equation of order n

J Convex Anal, 2011, 18: 673- 686

[本文引用: 1]

Yang D , Zhang W .

Characteristic solutions of polynomial-like iterative equations

Aequationes Math, 2004, 67: 80- 105

DOI:10.1007/s00010-003-2708-4      [本文引用: 1]

Yang L , Yang L , Yu Z , Zhang W .

Real polynomial iterative roots in the case of nonmonotonicity height \ge 2

Sci China Math, 2012, 55: 2433- 2446

DOI:10.1007/s11425-012-4537-1      [本文引用: 2]

Yu Z , Yang L , Zhang W .

Discussion on polynomials having polynomial iterative roots

J Symb Comput, 2012, 47: 1154- 1162

DOI:10.1016/j.jsc.2011.12.038      [本文引用: 2]

Zhang J , Yang L , Zhang W .

Some advances on functional equations

Adv Math (China), 1995, 24: 385- 405

URL     [本文引用: 1]

Zhang W .

Discussion on the iterated equation \sum\limits^{n}_{i=1}\lambda_{i}f^{i}(x)=F(x)

Chinese Sci Bull, 1987, 32: 1444- 1451

DOI:10.1360/csb1987-32-19-1444      [本文引用: 1]

Zhang W .

Discussion on the differentiable solutions of the iterated equation \sum\limits_{i=1}^{n}\lambda_{i}f^{i}(x)=F(x)

Nonlinear Anal, 1990, 15: 387- 398

DOI:10.1016/0362-546X(90)90147-9      [本文引用: 1]

Zhang W .

On existence for polynomial-like iterative equations

Results Math, 2004, 45: 185- 194

DOI:10.1007/BF03323006      [本文引用: 1]

Zhang W .

Solutions of equivariance for a polynomial-like iterative equation

Proc Roy Soc Edinburgh Sect A, 2000, 130: 1153- 1163

DOI:10.1017/S0308210500000615      [本文引用: 1]

Zhang W , Baker J A .

Continuous solutions of a polynomial-like iterative equation with variable coefficients

Ann Polon Math, 2000, 73: 29- 36

DOI:10.4064/ap-73-1-29-36      [本文引用: 1]

Zhang W , Zhang W .

On continuous solutions of n-th order polynomial-like iterative equations

Publ Math Debrecen, 2010, 76 (1/2): 117- 134

[本文引用: 1]

Zhang W , Xu B , Zhang W .

Global solutions for leading coefficient problem of polynomial-like iterative equations

Results Math, 2013, 63: 79- 93

DOI:10.1007/s00025-011-0162-5      [本文引用: 1]

/