ロジスティック写像
変数とパラメータ
ロジスティック写像を再掲する。
(1-1)
ここで、n はステップ数、n = 0, 1, 2... である。離散時間の発展を意味する。記号 t とも表記する場合もある。n = 1 を初期状態として n = 1, 2, 3... と数える場合もあるが、本記事中では前者で統一する。
xₙ は第 n 項目の x の値を示す。式(1-1)は、xₙ の値とその次の値である xₙ₊₁ の関係を示すものである。すなわち、xₙ から xₙ₊₁ が一意的に決定されることを示している。一方で、xₙ の値からその前の値である xₙ−1 の値は一意的に決定できない。二次関数の解が2つあることからも明らかなように、xₙ の値に辿り着く xₙ₋₁ の値については2つの値が有り得る。このような現在の値から過去の値への一意性が無いことは、1変数差分方程式でカオスを生み出す必要条件でもある。2変数以上になるとこの条件は必要ではなくなる。
x₀ を初期値として、ロジスティック写像では一般的に 0 ≤ x₀ ≤ 1 の範囲で初期値を取る。文献によっては 0 < x₀ ≤ 1 の範囲であるが、本記事中では前者で統一する。
a は定数で、パラメータとも呼ばれ、非線形性の強さを表すとされる。xₙ をある環境下における生物個体数とし、ロジスティック写像を個体群の成長モデルとみなすような文脈では、a は最大の生物個体数増加率 (biotic potential) とも呼ばれる。記号には r や μ なども割り当てることもある。パラメータは一般的に 0 ≤ a ≤ 4 の範囲で選択される。文献によっては 0 < a ≤ 4 の範囲であるが、本記事中では前者で統一する。
写像のグラフ
式(1-1)について、縦軸を xₙ₊₁、横軸を xₙ としてグラフを描くと、単純な二次関数の放物線が得られる。0 ≤ xₙ ≤ 1 の範囲でこのグラフを描くと、両端の xₙ = 0, 1 で、xₙ₊₁ = 0 となる。真ん中の xₙ = 0.5 で、xₙ₊₁ = a / 4 となる極大値かつ最大値を取る。このような xₙ から xₙ₊₁ への写像をリターン・マップと呼ぶ。
ロジスティック写像のような山形の曲線を描く写像は単峰写像と呼ばれる写像の一種で、その中でもロジスティック写像はシュワルツ導関数が負であるという条件を満たす。これらの条件を満たす写像では、後述するような分岐のパターンが発生する。他の例としては、正弦波を用いた xₙ₊₁ = a sin xₙ という写像でもこの条件を満たし、周期倍分岐、カオス、窓など、ロジスティック写像同様な振る舞いを見せる。このような曲線の写像ではカオスを生み出す基本機構である引き伸ばしと折り畳みという操作を備えているのが特徴である。a = 4 のときの xₙ から xₙ₊₁ への変換を考えると、まず区間 [0, 1] の xₙ が全体として2倍に引き伸ばされ、伸びた区間が半分に折り畳まれ、区間 [0, 1] の xₙ₊₁ として与えられる、というような操作を行っていることになる。このような引き伸ばしと折り畳みを繰り返す操作は、料理のパイをこねる操作に似ていることからパイこね変換と呼ばれ、低次元カオスの典型的な発生機構である。
定義域と値域
一般的には、ロジスティック写像は、式(1-1)による単位区間 [0, 1] から [0, 1] への写像とされる。この区間の写像であるために必要な初期値 x₀ とパラメータ a の範囲を確認する。xₙ が区間 [0,1] 内の値であれば、xₙ₊₁ が取り得る最小値は xₙ = 0, 1 のとき xₙ₊₁ = 0、最大値は xₙ = 0.5 のとき xₙ₊₁ = a/4 である。したがって、初期値が 0 ≤ x₀ ≤ 1 を満たし、なおかつパラメータが 0 ≤ a ≤ 4 を満たせば、次の値は 0 ≤ x₂ ≤ 1 となる。よって、その次の値も 0 ≤ x₃ ≤ 1 となり、数学的帰納法からすべての n で常に 0 ≤ xₙ ≤ 1 が成立する。
a が 4 を超えると、後述のように初期値に関わらず負の無限大へ発散する。a が負のときは、xₙ が負の値を取るようになるが、a = −2 までは有限な範囲で変動して発散はしない。ただし、多くの文献では前述の 0 ≤ a ≤ 4 の範囲を扱い、ロジスティック写像が生態学上のモデルとして研究された経緯もあってパラメータ a が負の場合を論じることは少ない。
0 ≤ a ≤ 1 のとき
パラメータ a が大きくなるに従った、変数 xₙ の振る舞いの変化を順に説明する。まず、a が 0 < a ≤ 1 の範囲では、初期値 x₀ が 0 から 1 までのどの値を取ったとしても、xₙ は n の増加と共に単調減少し、最終的に 0 に収束するこの収束する点を安定な固定点、安定な不動点、あるいは吸引的不動点、沈点などと呼ぶ。一般に、十分な時間経過後にその力学系が漸近する集合はアトラクタや吸引集合と呼ばれる。この場合は、x = 0 がアトラクタに相当する。
1次元写像の軌道を視覚的に表すのには、クモの巣図法と呼ばれる手法が利用される。これはグラフ上に原点を通る傾き1の直線を引き、この対角線と写像の曲線を交互に行き来することで1次元写像の軌道を2次元的に視覚化する方法である。a = 0.9 におけるクモの巣図と、n の増加に伴ったxₙの時系列データを示す。
上記のようにある1点に落ち着いて、その値から動かなくなったとき、写像は常に xₙ₊₁ = xₙ を満たす状態にある。よって、このときの値を xf とすると、
(2-1)
という関係式を満たす。単に不動点あるいは固定点と言えば、式(2-1)を満たす値 xf を指す。式(2-1)を解くと、xf は、
(2-2)
となる。0 < a ≤ 1 のとき、不動点 1 − 1/a は負または0の値となるので、0 ≤ a ≤ 1 のときは不動点1 − 1/aは現れず、不動点は0のみである。
式(2-2)から不動点の値は求められるが、xₙ ≠ xf のときに、n の増加と共に xₙ = xf に収束するかは保証されない。xf の近傍にある xₙ が最終的に xf に収束するならば、この不動点 xf は安定であるいう。逆に n の増加と共に xf から離れていくならば、この不動点 xf は不安定であるという。また、離れも近づきもしない場合は中立安定であるという。0 < a ≤ 1 のとき、0 は上記通り安定不動点に該当する。一方の不動点 1 − 1/a は、0 ≤ xₙ ≤ 1 の範囲では現れない不安定不動点である。
1 < a ≤ 2 のとき
1 < a ≤ 2 のときは、初期値が 0 または 1 である場合を除き、xₙ が単調増加あるいは単調減少しながら 1 − 1/a に収束する。以下にそれぞれの例を示す。
a < 1 から a を増加させていくと、安定不動点 0 と不安定不動点 1 − 1/a が a の増加と共に近づき、a = 1 で衝突し、1 − 1/a が安定不動点に、0 が不安定不動点に変わる安定性の交替が発生する。ただし、a = 1 ではまだ 0 へ収束し、a = 1 を超えたところで漸近的な振る舞いは変化する。このようなパラメータの変化により力学系の解の性質が変わることを分岐と呼ぶ。分岐が発生するパラメータの値を、この場合は a = 1 を分岐点と呼ぶ。a = 1 における分岐はトランスクリティカル分岐と呼ばれるタイプの分岐に相当する。このようなパラメータの変化に伴った分岐の発生の様子を示すために、横軸に写像のパラメータを取り、縦軸に長期間後の変数の値を示した、分岐図と呼ばれる図が用いられる。以下にパラメータ a が 0 から 3 までの分岐図を示す。
0も不動点として残っているが、不安定なので0に収束するのは初期値0、1のときだけである。これは以下の1 < a < 4 の範囲でも同様である。
2 < a ≤ 3 のとき
2 < a ≤ 3 のときは、初期値0、1を除いて、1 < a ≤ 2 のときと同様に 1 − 1/a に収束する。ただし 2 < a ≤ 3 のときは単調増加または単調減少ではなく、1 − 1/a の値を一端通り過ぎ、周囲で変動しながら 1 − 1/a に収束していく振る舞いとなる。a = 2.8 における例を以下に示す。
3 < a ≤ 3.4494897... のとき
3 < a ≤ 3.4494897... の範囲では、ほぼ全ての初期値において、xₙ は1つの値には収束せず、ある2つの値を交互に取る振る舞いを起こすようになる。交互に取られる2つの値は、周期点や平衡点と呼ばれる。今のような2周期のときは、周期点を2周期点、振る舞いを2周期軌道などと呼ぶ。不動点と同様の考え方で、周期点にも安定なもの、不安定なものが存在する。今の範囲における2周期点は安定周期点である。a = 3.2 のときの例を以下に示す。
a = 3 に達すると不動点 1 − 1/a が安定不動点から不安定不動点に移行する。a = 3 のときはまだ 1 − 1/a (= 2/3) に収束するが、2 < a < 3 のときよりも収束速度が遅い。a ≤ 3 における不動点への収束を安定な1周期軌道と見なせば、a = 3 を超えた点で、1周期から2周期へ周期が倍増する。この振る舞いの変化も分岐の一種で、周期倍分岐と呼ばれる。
xₙ₊₂ と xₙ の関係、すなわち式(1-1)の2回の反復合成写像は、次のように示される。
(2-3)
この2周期点をそれぞれ xf1、xf2 とする。xf1、xf2 は、その定義より、
(2-4)
を満たす値である。これを解けば xf1、xf2 の値は、
(2-5)
で示され、パラメータ a から求めることができる。
ほぼ全ての初期値から2周期点に到達するが、0 と 1 − 1/a は不安定不動点として残っている。そのため、初期値が、0、1 − 1/a、さらに、n を増加させる内に 0 または 1 − 1/a にたどり着く初期値であれば、2周期点ではなくこれらの不動点に到達する。ただし、これらの値の範囲は、初期値 0 < x₀ < 1 の実数集合に対して無視できるほど小さい。
3.4494897... < a ≤ 3.5699456... のとき
2周期軌道の状態から a を増加させていくと、a = 1 + √6 = 3.4494897... で安定な2周期軌道から安定な4周期軌道の状態に変化する。これは2周期から4周期への周期倍分岐である。この4周期点の値も2周期点と同様に a に依存して決まるが、2周期点のように陽関数の形で得ることはできない。さらに a の値を大きくしていくと4周期から8周期、8周期から16周期、16周期から32周期... と周期倍分岐が続けて発生していく。
この周期倍分岐は無限に続くが、一方で、周期倍分岐が発生する分岐点の間隔は等比数列的に減少する。1周期から2周期への分岐を1番目の周期分岐として数えれば、k 番目の分岐点で 2ᵏ 周期が発生する。k 番目の分岐点パラメータ a を aₖ と表すと、実際に16周期までの分岐点は a₁ = 3, a₂ = 3.44949..., a₃ = 3.54409..., a₄ = 3.56441... と、分岐点の間隔は減少していく。このため、k → ∞ の極限における分岐点を ac とすると、ac = 3.5699456... という有限な値に収束する。ac は集積点やファイゲンバウム点と呼ばれる。aₖ の減少の割合の極限は次式で示すような定数値となる。
(2-6)
この δ の値は、ミッチェル・ファイゲンバウムにより発見されたことからファイゲンバウム定数と呼ばれる。このような周期倍分岐を無限回繰り返した末にカオスへと遷移する現象はカオスへ至る道筋の一つで、ファイゲンバウムのシナリオや周期倍分岐ルートと呼ばれる。
a = ac のときのロジスティック写像は、無限個の周期点が存在する閉じることの無い周期軌道を取る。この値における軌道はファイゲンバウム・アトラクタや臨界2∞アトラクタと呼ばれる。ファイゲンバウム・アトラクタの構造はフラクタルになっており、容量次元は約0.54である。一方で、カオスの要件の一つである初期値鋭敏性は持たない。他のカオスではフラクタル構造と初期値鋭敏性を同時に備えるのに対し、ファイゲンバウム・アトラクタではフラクタル構造は持つが初期値鋭敏性は持たないのが特徴である。
式(2-4)と式(2-5)は k 周期のときにそのまま拡張できるので、k 周期点の j 番目の値を xfj とすれば、以下の関係を満たす値がそれぞれの周期点（k = 1 のときは不動点）となる。
(2-7)
3.5699456... < a < 4 のとき
ac = 3.5699456... を超えるとカオスが発生し、収束も周期性も無い不規則で複雑な挙動を示すようになる。初期値鋭敏性の指標であるリアプノフ指数 λ を計算すると、a < ac の範囲では λ は負または 0 の値の範囲に留まっていたが、a > ac の範囲から λ が正の値も取るようになる。ロジスティック写像のような1次元の区間力学系が生み出すカオスは、特に1次元カオスと呼ばれる。
ac < a < 4 の範囲で分岐図を見ると、例えば a = 3.8 などでは1つの区間内に軌道が収まっているが、a = 3.65 などでは2つの区間に軌道が分かれている。このような区間はバンドと呼ばれる。a がバンドを複数持つ領域にあるとき、軌道は規則的に順番に各バンドを巡り、しかし各バンド内での取る値は不規則である。バンドが2つのときを例にすると、n が偶数のときは下側のバンド内に存在し、奇数のときは上側のバンド内に存在するが、偶数、奇数それぞれのときの値はバンド内で不規則に決まる。このような状態のカオス軌道を、バンドカオスや周期的カオスと呼ぶ。a = 3.65 からさらに値を小さくしていくと、周期倍加分岐のときと同じように、4, 8, 16, ...2ᵏ というふうにバンドの数が2倍に増えていく。バンドが増えるときのパラメータ a の値を eₖ で表すとする。eₖ の値の間隔も、周期倍加分岐と同じく、バンドの増加と共に急激に小さくなっていき、k → ∞ で e∞ = ac となり、周期倍加分岐の集積点と一致する。バンドとバンドが一緒になることをバンド融合、分かれることをバンド分裂と呼び、これらもカオスの典型的な分岐の1つである。
以上のように、ac < a < 4 の範囲でバンドカオスの融合や分裂が発生しているが、軌道の最大値と最小値についてはバンドの数に関わらず決定される。これは後述の窓領域でも同様である。このような xₙ の最大最小範囲は、パラメータ a により以下のように得ることができる。
(2-8)
ac < a < 4 の範囲でカオス的振る舞いが発生するようになるが、常にカオス軌道を示す訳ではなく、a の値の領域によって、カオス軌道になったり周期軌道になったりする。このようなカオス軌道に至った後に周期軌道に変わるaの領域を窓や周期窓、窓領域と呼ぶ。カオスの非周期的領域と窓の周期的領域は交互に出現する。カオスから窓への分岐は、接線分岐と呼ばれるタイプの分岐により発生する。これはサドルノード分岐とも呼ばれ、1次元写像におけるサドルノード分岐を特に接線分岐と呼ぶ。
ac < a < 4 の範囲で発生する窓の詳細について見ていく。まず、窓の個数は、この範囲に無限個存在している。窓の周期については、ac < a < 4 の範囲に、3周期以上の全ての自然数 k の周期の窓が存在する。各 k 周期の軌道がそれぞれ1回ずつ発生する訳ではなく、kの値が大きいほど多く発生する。k 周期の窓の個数を Nₚ としたとき、k を素数に限定すると、k の値から Nₚ の値が以下の式より得られる。
(2-9)
式(2-9)は、k を素数に限定した式だが、実際の Nₚ を調べて比較すると、素数でない k に対しても非常に近い値の Nₚ を式(2-9)で得ることができる。
窓の幅（窓が始まる a と窓が終わる a の差）は、3周期の窓が最も広く、周期が大きいほど幅は狭まっていくと推定される。よって、窓は無限個存在するが、窓の周期的領域の範囲はカオス軌道が支配する領域に比べて小さい。概算によると、区間 [ac, 4] の約90%はカオス軌道が支配する範囲で、残りが窓領域の範囲である。
窓領域で現れる周期軌道は、窓領域ではカオスが背後に存在するが不安定のため観測されず周期軌道のみが可観測となって現れるもので、この点において集積点前で現れる周期軌道と異なる。窓領域中では、初期値から安定周期軌道に落ち着くまでの間のみに、この潜在的なカオスが過渡的に現れる。このようなカオスを過渡カオスと呼ぶ。
最も広い3周期窓を例にして、窓領域内の振る舞いについて見ていく。まずカオス発生領域から3周期窓発生領域へ移る過程を観察する。3周期窓の発生分岐点は 1 + 2√2 = 3.82843... で与えられる。この3周期窓分岐点を a₃ とする。a₃ よりもわずかに小さい a = 3.8282 のときの時系列を見てみると、xₙ が不規則的振る舞いを起こしている時間域と3周期の周期的振る舞いを起こしている時間域が交互に発生する。この振る舞いの周期的な部分をラミナーと呼び、不規則部分をバーストと呼ぶ。この振る舞いはカオスに相当するが、間欠的にしかカオス挙動が発生しないことから間欠カオスや間欠性カオスと呼ばれる。また、バーストとラミナーの時間域の長さに規則性は無く、不規則に変化する。次に、3.8282より大きいが a₃ には満たない a = 3.828327 のときの振る舞いを観察すると、a = 3.8282 のときよりラミナーの平均的な時間域長さが長くなり、バーストの平均的な時間域長さが短くなる。さらに a を大きくしていくとラミナーの長さがどんどん大きくなっていき、a₃に至ったところで完全な3周期に変わる。このような振る舞いの変化を、a を大きい方から小さい方へ変化させる方向で見てみると、周期軌道からカオス軌道への移行していることになる。これも前述の周期倍分岐ルートと同じくカオスに至る道筋の一種で、このような接線分岐による間欠カオスの発生を特徴とした道筋を、ポモウ・マンヴィルのシナリオやインターミッテンシールート(intermittency route)と呼ぶ。
次に、窓内部の構造を観察する。完全に3周期窓に入った後は、既に述べたように振る舞いは3周期軌道になる。分岐図で見ていくと、a₃ を過ぎたところからカオス領域から3本の曲線に変わり伸びていくが、a が大きくなるとある値でそれぞれの曲線が2つに分かれる。これは 3 < a ≤ 3.5699456... で見た周期倍加分岐と同じ現象が発生しており、全体で見ると、3周期、6周期、12周期...という風に 3 × 2ᵏ 周期が発生していく。最終的には、2ᵏ 周期分岐のときと同じくある値の集積点で無限周期に至り、カオスが発生する。集積点通過後も a 増加に伴う振る舞いの変化は、2ᵏ 周期分岐終了後と同じく、小領域を巡る周期的カオスから始まり、小領域の結合、窓の発生（大きな窓の中に窓がある構造で子窓などと呼ぶ）が起こる。1つの窓の中に、このような全体と同じような振る舞いの変化が存在することにより、ロジスティック写像の分岐図は、分岐図の中に同じ形状の分岐図が存在して、それが無限に続くというフラクタル的な自己相似形状になっている。ただしこの自己相似形状は、無限の入れ子構造となる厳密な自己相似形となってはいない。
最後に、3周期窓領域の終了部分について見ていく。全体と同じような振る舞いの変化を辿り、3周期窓領域の最終部では、3つの領域を巡る3バンドカオスとなる。この状態からさらに a を大きくしていくと、ある値で窓領域が終了し、3バンドカオスから1バンドカオスに戻る。この分岐点の値は a = 3.8568... で与えられる。このような窓の終わりで軌道が一気に大きな区間へ移る分岐はクライシスと呼ばれる。アトラクタが不安定周期点に接触して崩壊し、崩壊したアトラクタの外側のアトラクタに捕まることで、区間の変化が発生する。後述の境界クライシスと呼び分けるために内部クライシスと呼ばれる。あるいはバンド融合クライシスとも呼ぶ。
a = 4 のとき
a = 4 のときは、区間[0, 1]全体で変動する。このとき、ロジスティック写像は全ての自然数k周期軌道が存在するが、それら全ては不安定周期軌道となっており、周期軌道に漸近することなく非周期で区間 [0, 1] を巡り続ける軌道となる。この a = 4 におけるカオスは特別にピュアカオスとも呼んだりもする。
a = 4 のとき、初期値鋭敏性の指標であるリアプノフ指数の値は ln2 で、写像の複雑性を示す位相的エントロピーの値は log2 である。また、このときの区間 [0, 1] での xₙ の分布関数あるいは不変測度 μ(x) は次式で示される。
(2-10)
ここで、xₙ の値が微小区間 dx に訪れる確率は μ(x)dx で得られる。よって、xₙ の発生頻度分布は、区間 [0, 1] で一様分布ではなく、0 と 1 の両端が発生確率が高いU字型の分布となっている。分岐図のプロットの濃淡からも明らかなように、a < 4 のカオス発生範囲でも発生分布は一様ではなく偏りが存在している。
a = 4 のときの軌道が取り得る複雑さは、確率的なランダムに匹敵すると説明される。具体的には以下のようになっている。n = 0, 1, 2... に対して時系列xₙ = x₀, x₁, x₂... が得られる。次に、この時系列の値それぞれを、x = 0.5 以下のものは0、x = 0.5 を超えるものは1に変換する。変換された 0 か 1 の値を、時系列順に並べて数列sLを作る。例えば、初期値 x₀ = 0.2 とすれば、 x₁ = 0.64, x₂ = 0.9216, x₃ = 0.28901, ... となるので、sL = 0110... が得られることになる。一方、1/2 の確率でそれぞれ表裏が出るコイントスを想定する。表が出たときを 0、裏が出たときを 1 として、試行順に数列 sC を作る。例えば1回目表、2回目表、3回目裏、4回目裏...であれば、sC=0011...が得られる。試行は無限回行うとして、sC の 0 と 1 が無限に並ぶ場合も含める。この試行では確率論的なコイントスで記号を決めていったので、数列 sC は完全にランダムであらゆるパターンの数列を得ることができる。これに対して、ロジスティック写像から得られた数列 sL は決定論的な法則に従って決めていったものである。しかしながら、コイントスで得られる可能性のある全てパターンの数列 sC を、初期値を適切に選びさえすればロジスティック写像による sL によって実現できる。すなわち、無限列 sC に対して、sC = sL となる初期値x₀ が区間 [0, 1] 中にただ1点存在することが証明されている。
以上のような性質を持つ解は「真にカオス的」「ピュアカオス」とも形容される。一方で、a = 4 では、例外的なケースとして、カオスでありながらも n について明示的な厳密解を得ることができる。式(5-2)にそれらの解を示す。
4 < a のとき
a が 4 を超えるとほとんど全ての点 x₀ から始まる挙動はマイナス無限大へ発散する。この分岐は、窓領域の終わりで発生していた内部クライシスと同じ現象だが、アトラクタが不安定化しても外側の大きなアトラクタに移行することなく軌道が無限遠に発散することが特徴である。内部クライシスと呼び分けるために境界クライシスと呼ぶ。あるいはアトラクタ破壊のクライシスとも呼ぶ。
4 < a のときでも、マイナス無限大へ発散しない挙動を取る初期値 x₀ の集合が存在する。この集合を A とすると、単位区間 I から差し引いた I − A 内にある点が発散の一途をたどることになる。一方、発散しない集合 A は、長さが 0 でありながら非可算無限の点が存在するカントール集合の一種を形成する。
a < 0 のとき
以上まではパラメータ a が正の場合について説明してきた。ロジスティック写像は生態学上のモデルとして研究された経緯もあり、一般的に a が負の場合については論じられることは少ない。a が負の範囲のときの写像の分岐を確認すると、正の場合と似たような分岐を経て発散に至る。
a の値を0から減少させていくと、−1 までは xf = 0 の安定不動点に漸近するが、−1 を超えたところから2周期点に分岐し、正のときと同じく周期倍分岐を経てカオスへ至る。a = −2 のときはカオスだが、a = 4 のときと同じように時間発展の厳密解を得ることができ、式(5-3)で表される。最終的には、a が −2 を下回るとプラス無限大へ発散するようになる。以下に −2 から 4 までのロジスティック写像分岐図の全体像を示す。
個体群成長モデル
ロジスティック写像は、生物の個体数が世代ごとにどう変化するか予測する計算モデルとしてロバート・メイにより提案された。後述の#研究史に記す通り、それ以前からも式(1-1)については既に研究されていたが、この1970年代前半のメイの研究によりロジスティック写像が広く知られるようになった。
数理生物学的には、ロジスティック写像は、生物の個体群における個体数あるいは個体群密度の変動モデルを意味する。細かくは、離散世代型の単一種個体群モデルに相当する。ロジスティック写像の数理生物学的な導出の仕方は様々だが、まず最も単純な指数成長型の個体群モデルを導入して、それの拡張として生物学的条件を与えて導入する場合について説明する。
生物の個体数を N として、繁殖過程に関して世代が重ならないような生物種とする。このようなとき、n 世代における N の値に比例して増えていく単純なモデルが考えられる。
(6-1)
ここで c は、n によらない定数とし、正味の増殖率を意味する。個体当たりの結実率・出産率と繁殖可能になるまでの生存率との積とも考えられる。また、c が n によらない定数とは、環境条件が時間的に変化しない場合である。
この単純なモデルでは、個体数 N は指数関数的に無制限に増加することになる。一方で、密度効果を考えると、個体数が多くなると資源の奪い合いが激しくなり、繁殖率が低下すると仮定できる。これも単純なモデルの1つとして、次のような増殖率 c と個体数 N の関係を考えることができる。
(6-2)
ここで、a が最大繁殖率、b が繁殖率に対する密度効果の強さを表している。
式(6-1)と式(6-2)を合成すると、
(6-3)
となり、xₙ = b Nₙ/a と変数変換すれば、ロジスティック写像と同形の次式が得られる。
(6-4)
N の値が負のとき、個体数としての意味が成さないので、この制約から、a ≤ 4 という制約が生まれる。
1941年、昆虫学者の内田俊郎は離散世代型の生物種であるマメゾウムシの密閉容器中の繁殖実験で、個体数変動が振動しながら収束していくデータが得ている。これはロジスティック方程式では再現できないものだったが、ロジスティック写像の 2 < a ≤ 3 の場合で再現できる。また、ロジスティック写像の 3 < a ≤ 3.44949... に相当するような実験データとしては、内田俊郎がヨツモンマメゾウムシによる実験で、昆虫学者のニコルソンがヒツジキンバエによる実験で2周期振動型の個体数変化のデータを得ている。
ただし、ロジスティック写像のみが離散世代型単一種個体群モデルに適用できるモデルというわけではない。1953年に内田俊郎らも自身の実験データを別のモデルで再現している。また、ロジスティック写像ではパラメータの生物学的解釈が難しいとする意見もある。ロジスティック写像型のモデル以外にも、様々な種類のモデルが数理生物学で研究されており。リッカーモデルなどもパラメータ変化によるカオス的振る舞いも見せる。また、実験室の厳密に管理された条件では実際の生物個体群がカオス的変動を示した例はあるが、野外環境での個体群が本当にカオス的変動を起こすのかどうかは、様々な環境変動の影響があることからまだ論争が続いている。
研究史
式(1-1)がロジスティック写像と呼ばれる前の時代としては、1947年にスタニスワフ・ウラムとジョン・フォン・ノイマンが、
(7-1)
を使用すれば疑似乱数を得ることができることを指摘している。
式(7-1)は a = 4 のロジスティック写像に相当する。この頃はカオスという言葉もまだ存在していない。また、ウラムとノイマンはこの式とテント写像が位相同型の関係にあること、この式の不変測度が一様分布ではないことも明らかにしていた。
1950年にプルキン(S. P. Pulkin) が連続写像や区分的に連続な写像から振動の繰り返し系列が生じることを示したが、この論文中にもロジスティック写像と同じ式が現れる。その後、1950年代前半から1960年代後半にかけて、他の研究者からも任意のパラメータ a を備えた形式についても研究されている。ミュルバーク(Pekka Juhana Myrberg) は1958年から1963年にかけて、xₙ₊₁ = xₙ² − λ という形式の1次元2次写像の分岐現象について調べ、1963年にはこの写像による周期倍加分岐も研究している。カオス現象自体の発見（可視化）は、1961年の上田睆亮や1963年のエドワード・ローレンツなどで達成されるが、まだこの時点でもカオスという共通語は存在していない。
その後、数理生物学者のロバート・メイが、生態学の問題に取り組む過程で式(1-1)を用いた研究を1970年初頭から開始する。1973年に、メイはロジスティック方程式を離散時間化により差分方程式に変換させて、式(1-1)を求めた。メイはこの式を基に数値実験を行い、パラメータaによる変数の振る舞いを変化を研究した。1974年と1976年にメイはこれらの研究成果を発表し、これらの論文によって、カオスと関係の深い写像としてロジスティック写像が広く知られるようになる。さらに、ロジスティック写像型の関数が持つ非周期軌道と周期軌道の存在について独立に研究していたジェームス・ヨーク(James A. Yorke) らが、メイの研究に刺激を受ける格好で、1975年に "Period Three Implies Chaos"（3周期はカオスを意味する）の題名でそれらの研究を発表し、カオスという共通語が広まっていく。
その後、メイは先行研究があったことを尊重した上で自身の功績について次のように述べている。
このようにしてわれわれの研究は、カオスを科学の本流に押し上げる第2の主役となった。ジム・ヨークが言ったように、われわれは2次写像の奇怪な数学的挙動を独立に最初に発見したのではなく、科学におけるその広範な意味づけを最後に強調した研究者たちなのだ!
—Robert M. May、(エイブラハム・ウエダ 2002, p. 160)、稲垣耕作・赤松則男訳
また、ロバート・デバニー(Robert L. Devaney) は、著作でロジスティック写像の解説に入る前に次のように述べている。
これにより、単に2次関数 fλ(x) = λx(1 − x)（これもまたロジスティック写像と呼ばれる）を反復すれば、最初の個体数 x₀ の運命が予測できるようになるというわけだ。簡単な話のように聞こえるだろうが、あえて言い添えておくならば、この単純な2次関数の反復が完全に理解できるようになったのは、何百人もの数学者の努力の末、やっと1990年代の終わり頃になってからのことなのである。
—Robert L. Devaney、(Hirsch et al. 2007, p. 347)、三波篤朗訳