主題
Search

二次對映

下載 Mathematica Notebook下載 Wolfram Notebook

二次對映是如下形式的二次遞推方程

x_(n+1)=a_2x_n^2+a_1x_n+a_0.
(1)

雖然一些二次對映可以以閉合形式求解（例如，logistic 對映的三個可解情況），但大多數都不能。

一個具有閉合形式解的二次對映的簡單例子是

x_n=x_(n-1)^2
(2)

其中 x_0=2，其解為 x_n=2^(2^n)，當 n=0, 1, ... 時，前幾項為 2, 4, 16, 256, 65536, 4294967296, ... (OEIS A001146)。

另一個例子是高度 <=n 的“強”二叉樹的數量，由下式給出

y_n=y_(n-1)^2+1
(3)

其中 y_0=1。 前幾項是 2, 5, 26, 677, 458330, 210066388901, 44127887745906175987802, ... (OEIS A003095) 這個遞推關係具有“解析”解

y_n=|_c^(2^n)_|,
(4)

其中

c=exp[sum_(j=0)^(infty)2^(-j-1)ln(1+y_j^(-2))]
(5)
=1.502836801...
(6)

(OEIS A077496) 且 |_x_| 是向下取整函式（Aho 和 Sloane 1973）。

第三個例子是數字 1 最接近的嚴格下逼近，

S_n=sum_(i=1)^n1/(z_i),
(7)

其中 1<z_1<...<z_n 是整數。 解由遞推關係給出

z_n=z_(n-1)^2-z_(n-1)+1,
(8)

其中 z_1=2。 得到的序列被稱為 Sylvester 序列，前幾項為 2, 3, 7, 43, 1807, 3263443, 10650056950807, ... (OEIS A000058)。 這有一個閉合解為

z_n=|_d^(2^n)+1/2_|
(9)

其中

d=1/2sqrt(6)exp{sum_(j=1)^(infty)2^(-j-1)ln[1+(2z_j-1)^(-2)]}
(10)
=1.2640847353...
(11)

(OEIS A076393；Aho 和 Sloane 1973，Vardi 1991，Graham 等人 1994)。

著名的遞推關係

x_(n+1)=x_n^2+c
(12)

通常被稱為“二次對映”的這個遞推關係通常不能以閉合形式求解。 這是定義 Mandelbrot 集合的復對映的實數版本。 此對映的不動點發生在

x^((1))=[x^((1))]^2+c
(13)
(x^((1)))^2-x^((1))+c=0
(14)
x_+/-^((1))=1/2(1+/-sqrt(1-4c)).
(15)

週期為 2 的不動點發生在

x_(n+2)=x_(n+1)^2+c
(16)
=(x_n^2+c)^2+c
(17)
=x_n^4+2cx_n^2+(c^2+c)
(18)
=x_n.
(19)

去掉下標並因式分解得到

x^4+2x^2-x+(cx^2+c)=(x^2-x+c)(x^2+x+1+c)=0.
(20)

因此解發生在

x_+/-^((2))=1/2(1+/-sqrt(-3-4c)).
(21)

週期為 3 的不動點發生在

x^6+x^5+(3c+1)x^4+(2c+1)x^3+(c^2+3c+1)x^2 +(c+1)^2x+(c^3+2c^2+c+1)=0.
(22)

另一個具有閉合形式解的二次對映的例子是以下情況

a_0=((a_1-4)(a_1+2))/(4a_2).
(23)

這有解

x_n=(r^(2^n)+r^(-2^n)-1/2a_1)/(a_2),
(24)

其中

r=1/4a_1+1/2x_0a_2+1/4sqrt((a_1+4+2x_0a_2)(a_1-4+2x_0a_2)).
(25)

類似地，以下情況

a_0=(a_1(a_1-2))/(4a_2)
(26)

有解

x_n=((2r)^(2^n)-1/2a_1)/(a_2),
(27)

其中

r=1/4a_1+1/2x_0a_2
(28)

(Little)。

在原點具有橢圓不動點的最一般的二階二維對映具有以下形式

x^'=xcosalpha-ysinalpha+a_(20)x^2+a_(11)xy+a_(02)y^2
(29)
y^'=xsinalpha+ycosalpha+b_(20)x^2+b_(11)xy+b_(02)y^2.
(30)

為了保持面積不變，該對映的行列式必須為 1，從而將獨立引數的數量從七個減少到三個。 然後可以透過縮放和旋轉將對映置於標準形式，得到

x^'=xcosalpha-ysinalpha+x^2sinalpha
(31)
y^'=xsinalpha+ycosalpha-x^2cosalpha.
(32)

逆對映是

x=x^'cosalpha+y^'sinalpha
(33)
y=-x^'sinalpha+y^'cosalpha+(x^'cosalpha+y^'sinalpha)^2.
(34)

不動點由下式給出

x_i^2sinalpha+2x_icosalpha-x_(i-1)-x_(i+1)=0
(35)

對於 i=0, ..., n-1。

另請參閱

Bogdanov 對映, Hénon 對映, Julia 集, 線性遞推方程, Logistic 對映, Lozi 對映, Mandelbrot 集, 二次的, 二次方程, 二次公式, 二次遞推方程, 遞推方程, 奇異吸引子, Sylvester 序列

使用 探索

參考文獻

Aho, A. V. 和 Sloane, N. J. A. "一些雙重指數序列。" Fib. Quart. 11, 429-437, 1973。Finch, S. R. "二次遞推常數。" §6.10 in 數學常數。 Cambridge, England: Cambridge University Press, pp. 443-448, 2003。Graham, R. L.; Knuth, D. E.; 和 Patashnik, O. 《具體數學：計算機科學基礎》，第 2 版中的研究問題 4.65。 Reading, MA: Addison-Wesley, 1994。Little, M. "二次對映特殊情況的新穎精確解。" http://homepage.ntlworld.com/little_mm/quadrams.html。Sloane, N. J. A. “整數序列線上百科全書”中的序列 A000058/M0865, A001146/M1297, A003095/M1544, A076393, 和 A077496。Sprott, J. C. "奇異吸引子的自動生成。" Comput. & Graphics 17, 325-332, 1993。 重印於混沌與分形，計算機圖形之旅：十年高階研究彙編 (Ed. C. A. Pickover)。 Amsterdam, Netherlands: Elsevier, pp. 53-60, 1998。Vardi, I. "所有歐幾里得數都是無平方因子數嗎？" 和 "PowerMod 來救援。" 《Mathematica 中的計算娛樂》中的 §5.1 和 5.2。 Reading, MA: Addison-Wesley, pp. 82-89, 1991。

在 中被引用

二次對映

引用為

Weisstein, Eric W. "二次對映。" 來自 —— 資源。 https://mathworld.tw/QuadraticMap.html

主題分類