關(guān)于Halley法Julia集的對稱性

劉 剛，陳少林，李瀏蘭

（衡陽師范學(xué)院 數(shù)學(xué)與計(jì)算科學(xué)系，湖南 衡陽 421002）

0 引言及主要結(jié)果

給定多項(xiàng)式f（z），則如下定義的公式

稱為關(guān)于多項(xiàng)式f的Halley法。該法是一種重要的迭代找根算法，其被關(guān)注的程度僅次牛頓法。計(jì)算數(shù)學(xué)工作者對該法進(jìn)行了大量的研究并取得非常豐富的成果，而從復(fù)動力系統(tǒng)的角度進(jìn)行的研究以及取得的成果相對較少。K.Kneis［1］研究了Hf不動點(diǎn)的性質(zhì)并給出了一些特殊多項(xiàng)式的Halley法的Julia集的分形圖集。文獻(xiàn)［2］對Hf的Julia集的拓?fù)浣Y(jié)構(gòu)進(jìn)行了研究。

Julia集是復(fù)動力系統(tǒng)的基本研究對象，我們關(guān)心Julia集的拓?fù)?，分形，幾何等性質(zhì)，鑒于多項(xiàng)式在迭代找根算法下一般都為有理函數(shù)，本文只在有理函數(shù)范圍內(nèi)進(jìn)行研究。有理函數(shù)Julia集的對稱性所研究的是有理函數(shù)的Julia集的自相似性，涉及的是Julia集的幾何性質(zhì)。A.F.Beardon［3］研究了臨界有限的有理函數(shù)的Julia集的對稱性。隨后G.Levin［4］研究了臨界有限的有理函數(shù)的Julia集的對稱性。D.Boyd［5］分析了有理函數(shù)的Julia集具有平移不變性的情形。文獻(xiàn)［6］對三類迭代找根算法Julia集對稱性進(jìn)行了研究。

給定有理函數(shù)R（z），則R的Julia集的對稱群∑（R）為保持J（R）不變的共形歐氏變換σ構(gòu)成的群，即∑（R）＝｛σ：σ（z）＝e2πθiz＋b（θ∈［0，1］）且σ（J（R））＝J（R）｝。給定非線性多項(xiàng)式f（z）＝anzn＋an－1zn－1＋…＋a1z＋a0（an≠0），我們稱為f（z）的形心。如果an＝1且an－1＝0，則這樣的多項(xiàng)式稱為標(biāo)準(zhǔn)多項(xiàng)式。若f（z）為非標(biāo)準(zhǔn)的多項(xiàng)式，令，則T?f?T－1為標(biāo)準(zhǔn)多項(xiàng)式。結(jié)合關(guān)于多項(xiàng)式的Halley法的動力學(xué)性質(zhì)和關(guān)于有理函數(shù)Julia集對稱群的知識，我們獲得了如下結(jié)果：

定理1 設(shè)f（z）為度大于1的標(biāo)準(zhǔn)多項(xiàng)式，則∑（f）?∑（Hf）。

定理2 設(shè)f（z）為度大于1的多項(xiàng)式，則平移變換σ＝z＋1∈∑（Hf）當(dāng)且僅當(dāng)J（Hf）為一水平直線。

定理3 設(shè)f（z）為度大于1的多項(xiàng)式，則J（Hf）為一條直線當(dāng)且僅當(dāng)f（z）＝c（z＝－z1）n（z－z2）n（c∈?＼｛0｝，z1≠z2）。此時(shí)J（Hf）為連接z1與z2線段垂直平分線。

推論1 設(shè)f（z）為度大于1的多項(xiàng)式，則σ＝z＋1∈∑（Hf）當(dāng)且僅當(dāng)f（z）＝c（z－z1）n（z－z2）n（c∈?＼｛0｝，z1≠z2，Re（z1）＝Re（z2））。

1 相關(guān)定義及引理

設(shè)R為復(fù)球面?－上的有理函數(shù)，，其中P和Q為互素的多項(xiàng)式，degR＝max（degP，degQ）稱為R的度。記Rn為R的第n次迭代。若存在m，k＞0使得Rk（Rm（z0））＝Rm（z0），則稱z0為R的預(yù)周期點(diǎn)。特別地，當(dāng)m＝0時(shí)，z0稱為周期點(diǎn)，相應(yīng)地最小的的k稱為周期點(diǎn)z0的周期。若｜（Rk）′（z0）｜＜1，｜（Rk）′（z0）｜＝1或者｜（Rk）′（z0）｜＞1，則z0稱為（超）吸性，中性或者斥性周期點(diǎn)。R的Julia集J（R）是定義為其斥性周期點(diǎn)的閉包。J（R）補(bǔ)集，即?－＼J（R），稱為R的Fatou集F（R）。J（R）是完全不變的，即R（J（R））＝J（R）＝R－1（J（R）。另外，J（R）沒有完全不變的真閉子集。關(guān)于有理函數(shù)動力系統(tǒng)更多的定義和結(jié)果參見文獻(xiàn)［7，8］。

為證明文中結(jié)果，我們還需以下引理。

引理1［3］設(shè)f（z）為度大于1的多項(xiàng)式，則對稱群∑（f）是由以f（z）的形心為中心的旋轉(zhuǎn)變換構(gòu)成的。如果f（z）是標(biāo)準(zhǔn)的，而且∑（f）是有限的，則∑（f）的階為使得f（z）可表示成f（z）＝zrQ（zm）的最大整數(shù)m，其中Q（z）為一多項(xiàng)式。如果對稱群∑（f）是無限的，那么J（f）是一個(gè)圓周，并且f（z）共軛于zn，其中n＝deg（f）。

引理2［5］設(shè)R（z）為度大于1的有理函數(shù)，且J（R）＋1＝J（R）。如果無窮遠(yuǎn)點(diǎn)∞是預(yù)周期的，那么J（R）是一水平直線或全平面?－。

引理3［1］設(shè)f（z）為度大于1的多項(xiàng)式，則多項(xiàng)式f的每一個(gè)根都是Hf的（超）吸性不動點(diǎn)。

引理4［2］設(shè)f為多項(xiàng)式，g（z）＝az＋b（a≠0）為仿射變換，則g?Hcfog?g－1＝Hf，（c∈?＼｛0｝）。

2 定理的證明

定理1的證明 如果∑（f）是有限的，由引理1知，對任意的σ∈∑（f），存在正整數(shù)k使得m為使得f表示成f（z）＝zrQ（zm）的最大整數(shù)。從而

易見σ?Hf?σ－1（z）＝Hf（z），則σJ（Hf）＝J（Hf），從而有σ∈∑（Hf），即證

如果∑（Hf）是無限的，由引理1有f（z）＝zn，其中n＝deg（f）。從而，即知J（Hf）＝∞，顯然有∑（Hf）＝｛σ∶σ（z）＝e2πθiz＋b（θ∈［0，1］）｝。注意到J（f）為單位圓周，從而∑（Hf）＝｛σ∶σ（z）＝e2πθiz（θ∈［0，1］）｝?∑（Hf）。

定理2的證明 當(dāng)J（Hf）是一水平直線時(shí)，顯然平移變換σ＝z＋1∈∑（Hf）。

如果平移變換σ＝z＋1∈∑（Hf），則Hf（z）滿足引理2的條件，故J（Hf）是一水平直線或全平面，然而多項(xiàng)式f（z）至少有一個(gè)零點(diǎn)，由引理3，Hf至少有一個(gè)吸性不動點(diǎn)，故J（Hf）≠?－，即知J（Hf）為一水平直線。

以下證明當(dāng)n1≠n2時(shí)，J（Hf＊）不能為一直線。注意到

因此｜Hf＊（z）－i｜＜｜z－i｜。因?yàn)閕為Hf＊的（超）吸性不動點(diǎn)，故｛z∶Im（z）≥3｝包含在i的所在的Fatou分支內(nèi)。再注意到任意的水平線｛z∶Im（z）＝t｝∪｛∞｝在Hf＊下不是前向不變的，此時(shí)J（Hf＊）不可能為一條直線。

接下來證明當(dāng)n1＝n2?n時(shí)，J（Hf＊）為帶無窮點(diǎn)的實(shí)軸。此時(shí)

易見Hf＊（?∪｛∞｝）＝?∪｛∞｝。

下證Hf－1＊（?∪｛∞｝）＝?∪｛∞｝。因Hf＊（∞）＝∞且，故∞為Hf＊的斥性不動點(diǎn)，從而∞∈J（Hf＊），又由所以∈J（Hf＊）。對于任意的x，因

故Hf＊（x）在三個(gè)區(qū)間以及都是從－∞單調(diào)遞增到0；而Hf＊（x）在三個(gè)區(qū)間以及從0單調(diào)遞增到＋∞。從而對任意的x∈?，H－1f＊（x）∈?，可知J（Hf＊）＝?∪｛∞｝.由仿射變換的性質(zhì)，可知J（Hf）＝g（J（Hf＊））為連接z1與z2線段的垂直平分線。

推論1的證明 結(jié)合定理2和定理3的結(jié)論即證。

［1］Kneisl K.Julia sets for the super Newton method，Cauchy's method，and Halley's method［J］.Chaos，2001，11：359－370.

［2］Wang X Y，Yu X J.Julia sets for the standard Newton's method，Halley's method，and Schr?der's method for multiple root［J］.Appl.Math.Comput.，2007，189（2）：335－338.

［3］Beardo A F.Symmetries of Julia sets［J］.Bull.London Math.Soc.，1990，22：575－582.

［4］Levin.G.Symmetries on a Julia set［J］.Adv.in Sov.Math.，1991，3：131－141.

［5］Boyd D.Translation invariant Julia sets［J］.Proc.Amer.Math.Soc.，2000，128：803－812.

［6］Yang W F.Symmetries of the Julia sets of Newton's method for multiple root［J］.Appl.Math.Comput.，2010，217（6）：2490－2494.

［7］Beardon A F.Iteration of Rational Functions［M］.Berlin：Springer，1991.

［8］喬建永.重整化變換的復(fù)動力學(xué)［M］.北京：科學(xué)出版社，2010.