孫霞 楊慧章

【摘要】對于中學數學競賽和高考中出現的分式線性遞推數列，它其實是復解析動力系統(tǒng)中Mbius變換迭代的一種特殊情形.本文給出了Mbius變換n次迭代的具體表達式，也給出了分式線性遞推數列的通項公式，從而使得求分式線性遞推數列通項公式簡單化.

【關鍵詞】Mbius變換;迭代;分式線性遞推數列

【基金項目】云南開放大學科學研究基金項目（19YNOU01）;云南省教育廳科學研究基金項目（2020J0492，2018JS479）;云南省地方高校聯合專項面上項目（2018FH001-014）;紅河學院第二屆中青年學術骨干（2015GG0207）

分式線性遞推數列是由 a1=a，an+1=can+bdan+e（ce-bd≠0）所確定的數列，它經常出現在中學數學競賽和高考的壓軸題中，我們看到有很多求通項公式的辦法.在本文中，我們站在復解析動力系統(tǒng)的角度，闡述了分式線性遞推數列只是Mbius變換迭代的一種特殊形式，且從動力系統(tǒng)的角度，給出了分式線性遞推數列的通項公式，使得求分式線性遞推數列通項公式簡單化，揭開了該類通項公式求解的神秘面紗.

我們從如下幾個方面來進行介紹.

1 Mbius變換的迭代

在復分析中，形如R（z）=az+bcz+d（ad-bc≠0）的映射稱為分式線性變換或Mbius變換，這樣的映射我們也稱為一次有理函數.

復解析動力系統(tǒng)是復分析中的一個分支，它是研究函數迭代序列或軌道狀態(tài)的.即倘若f是定義域S到其自身的解析映射，對任意初值z0∈S，考慮迭代過程z1=f（z0），z2=f（z1），…，zn=f（zn-1），…稱迭代序列{z0，z1，…，zn，…}為點z0（在f作用下）的軌道，復解析動力系統(tǒng)的研究對象就是上述的函數迭代序列.

復數域C={（a，b），a∈R，b∈R}，而C∪{∞}稱為整個Riemann球面.對z0∈C∪{∞}，我們定義z1=R（z0），z2=R（z1）=R（R（z0））=R2（z0），…，zn=R（zn-1）=R·R·…·R（z0）n=Rn（z0）.

定義 設R（z）為有理函數，點z稱為R（z）的不動點，如果R（z）=z.

引理 若R（z）是度為d的有理函數，則R（z）在C∪{∞}中只有d+1個不動點.

對于Mbius變換R（z）=az+bcz+d（ad-bc≠0）而言，由于它是一次有理函數，所以由上述引理知它在整個Riemann球面就只有2個不動點，接下來，我們從不動點出發(fā)來看Mbius變換的n次迭代表達式.

情形1 若Mbius變換R（z）在C∪{∞}中有且僅有一個二重不動點ζ，如果該不動點為∞，即ζ=∞，則R（z）必定為一次多項式，且R（z）=z+β（β≠0），這時，Rn（z）=z+nβ（n=1，2，3，…）.如果ζ≠∞，則令g（z）=1z-ζ，使得g（ζ）=∞，此時，g-1（z）=1z+ζ，不妨設S（z）=g·R·g-1（z），

那么S（z）也是Mbius變換且它與R（z）共軛，它僅以∞為其不動點，即S（∞）=∞，從而可以定出β ～，使得S（z）=z+β ～.

所以Sn（z）=z+nβ ～（n=1，2，3，…）.

因為Sn（z）=g·Rn·g-1（z），

所以g·Rn·g-1（z）=z+nβ ～，

所以g·Rn（z）=g（z）+nβ ～，

Rn（z）=g-1（g（z）+nβ ～）=1g（z）+nβ ～+ζ

=11z-ζ+nβ ～+ζ.

上述情形中，我們首先處理了較為簡單的情況，即不動點ζ=∞的情況，此時，Rn（z）=z+nβ（n=1，2，3，…）;若ζ≠∞，則做了分式線性變換g（z）=1z-ζ，并讓S（z）與R（z）共軛，從而轉化為以∞為不動點的情形.

情形2 若Mbius變換R（z）在C∪{∞}中有且僅有兩個判別的不動點ζ1和ζ2.先考查一種特殊情況，若ζ1=0，ζ2=∞，這時，R（z）=kz（k≠0，k≠1），于是，Rn（z）=knz（n=1，2，3，…）.對于一般情形的ζ1和ζ2，做分式線性變換h（z）=z-ζ1z-ζ2，使得h（ζ1）=0，h（ζ2）=∞，繼續(xù)令S（z）=h·R·h-1（z），則S（z）是Mbius變換，且0和∞均為其不動點，S（z）不恒等于z，所以存在A∈C，使得S（z）=Az（A≠0，A≠1），于是，Sn（z）=Anz，所以

Rn（z）=h-1·Sn·h（z）.

由于h-1（z）=ζ1-zζ21-z，所以

Rn（z）=h-1·Sn·h（z）=h-1·Snz-ζ1z-ζ2

=h-1Anz-ζ1z-ζ2

=ζ1-Anz-ζ1z-ζ2·ζ21-Anz-ζ1z-ζ2.

同樣地，在該情形中，我們也用了轉化的思想，首先考慮簡單的情形，即以ζ1=0，ζ2=∞為不動點的迭代，接下來，當ζ1≠0，ζ2≠∞時，則做分式線性變換h（z）=z-ζ1z-ζ2，使得h（ζ1）=0，h（ζ2）=∞，從而轉化為我們已經熟知的情形.

2 Mbius變換的迭代與分式線性遞推數列的通項公式

上述情形1中，若R（z）在C∪{∞}中有且僅有一個不動點ζ且ζ=∞，則R（z）=z+β（β≠0），Rn（z）=z+nβ（n∈N），此時，該迭代序列是我們的公差為β的等差數列.情形2中，若R（z）在C∪{∞}中有且僅有兩個判別的不動點ζ1和ζ2，且ζ1=0，ζ2=∞，則R（z）=kz（k≠0，k≠1），此時，Rn（z）=knz（n=1，2，3，…），該迭代序列則是我們的公比為k的等比數列.

接下來，我們通過具體例子來看對于一般的分式線性遞推數列，如何求出它的通項公式.上述討論的Mbius變換是定義在整個Riemann球面的，而我們的數列則只是定義在正整數域上的，為了計算方便，在計算過程中不妨把復變量z換為實變量x.

例1 已知a1=5，an+1=an-4an-3（n≥1），求通項公式an.

解 取R（x）=x-4x-3，令R（x）=x，得不動點為x1=x2=2，所以

g（x）=1x-2，g-1（x）=1x+2，

S（x）=g·R·g-1（x）=g·R1x+2=g1-2x1-x=x-1，

Sn（x）=x-n.

因為Sn（x）=g·Rn·g-1（x）

所以Rn（x）=g-1·Sn·g（x）

=g-1·Sn1x-2=g-11x-2-n

=g-1-nx+（2n+1）x-2

=x-2-nx+（2n+1）+2

=（1-2n）x+4n-nx+（2n+1）.

因為x0=a1=5，所以數列的通項

an=xn-1=Rn-1（x0）=Rn-1（5）

=[1-2（n-1）]·5+4（n-1）-（n-1）·5+[2（n-1）+1]=-6n+11-3n+4=6n-113n-4.

事實上，我們也可以直接用情形1中的公式，因為該函數只有一個二重的不動點.由于ζ=2，β～=-1，所以

Rn-1（x）=11x-2-（n-1）+2=x-2（n-1）（x-2）1-（n-1）（x-2），

an=Rn-1（5）=-6n+11-3n+4.

例2 已知a1=3，an+1=2anan+1，求數列通項an.

解 取R（x）=2xx+1，令R（x）=x，得不動點為x1=0，x2=1，做變換g（x）=xx-1，則g-1（x）=xx-1，從而

S（x）=g·R·g-1（x）=g·Rxx-1=g2x2x-1=2x，

所以Sn（x）=2nx.

由于Sn（x）=g·Rn·g-1（x），

所以Rn（x）=g-1·Sn·g（x）=g-1·Snxx-1

=g-12nxx-1=2nx（2n-1）x+1.

因為x0=a1=3，所以

an=xn-1=Rn-1（x0）=Rn-1（3）=2n-1·3（2n-1-1）·3+1

=3·2n-13·2n-1-2.

我們也可以用情形2中的公式，此時ζ1=0，ζ2=1，A=2，所以

Rn-1（x）=-2n-1xx-11-2n-1xx-1，

an=Rn-1（3）=-2n-1·321-2n-1·32=-3·2n-12-3·2n-1=3·2n-13·2n-1-2.

對于次數為1的有理函數，其迭代后平凡，可以具體寫出迭代的表達式，也就是文中所述的Mbius變換的迭代.而對于次數大于或等于2的有理函數的迭代，其狀況就要復雜得多.從復解析動力系統(tǒng)的角度來看，分式線性遞推數列其實來源于分式線性變換的迭代，它只是一次有理函數迭代中變量取正整數的情形，所以可以寫出具體的表達式.

【參考文獻】

[1]趙嘉璐.用不動點法求分式型遞推數列的通項[J].數學學習與研究，2018（09）：107.

[2]林國夫.利用函數不動點求數列的通項公式[J].數學通報，2008（12）：44-45，47.

[3]史濟懷，劉太順.復變函數[M].合肥：中國科學技術大學出版社，1998.

[4]呂以輦.復解析動力系統(tǒng)[M].北京：科學出版社，1997.

[5]任福堯.復解析動力系統(tǒng)[M].上海：復旦大學出版社，1997.

[6]Alan F.Beardon.Iteration of Rational Functions[M].Berlin：Springer Verlag，1991.