池田写像のソースを表示

[[物理学]]や[[数学]]において、'''池田写像''' (いけだしゃぞう、{{lang-en|Ikeda map}}) は以下の[[複素解析]]写像で与えられる離散時間[[力学系]]である。

:<math> z_{n+1} = A + B z_n e^{i (|z_n|^2 + C)} </math>

オリジナルの池田写像は非線形[[光共振器]] ({{仮リンク|リングレーザー|en|Ring laser}}、[[非線形システム論|非線形]][[誘電体]]を含む) における光の軌跡のモデルとして[[池田研介]]により提案された。 
池田研介、大同寛明、秋元興一により上記の単純化した形に一般化された<ref>K.Ikeda, Multiple-valued Stationary State and its Instability of the Transmitted Light by a Ring Cavity System, Opt. Commun. 30 257-261 (1979);  K. Ikeda, H. Daido and O. Akimoto, Optical Turbulence: Chaotic Behavior of Transmitted Light from a Ring Cavity, Phys. Rev. Lett. 45, 709–712 (1980)</ref>。

<math>z_n</math>は[[共振器]]内の回転のn番目のステップにおける共振器内の電場を表し、<math>A</math>と<math>C</math>はそれぞれ、外部からのレーザー光線と共振器の線形位相を示すパラメータである。

特に、パラメータ<math>B(<=1)</math>は振る舞いを特徴づける散逸パラメータと呼ばれ、<math>B=1</math>の極限で池田写像は[[保存系]]になる。  

非線形誘電体の飽和効果を考慮に入れ、オリジナルの池田写像に修正を施した以下の形で使用される事が多い。

:<math> z_{n+1} = A + B z_n e^{i K/(|z_n|^2 +1)+C} </math>

上記を実2次元平面に描く場合、

:<math> x_{n+1} = 1 + u (x_n \cos t_n - y_n \sin t_n), \, </math>

:<math> y_{n+1} = u (x_n \sin t_n + y_n \cos t_n), \, </math>

となる。ここで、''u''はパラメータであり、 

:<math> t_n = 0.4 - \frac{6}{1+x_n^2+y_n^2}. </math>

である。幾つかの値の''u''では、この系は[[カオス理論|カオス]]な振る舞いを示す。

==アトラクター==
[[media:Ikeda map.ogg|この動画]]はパラメータ<math>u</math>を0.0から1.0まで0.01ずつ変化させた時の挙動を示している。シミュレーションに使用したステップ数は500、ランダムに選択した始点は20000個である。[[アトラクター]]を描くため、それぞれの軌道で20個の点がプロットされている。<math>u</math>の値が増加するにつれ、アトラクターポイントの分岐は大きくなる。 

{|
| [[image:Ikeda0300.png|thumb|left| <math>u=0.3</math>]]
| [[image:Ikeda0500.png|thumb|left| <math>u=0.5</math>]]
|-
| [[image:Ikeda0700.png|thumb|left| <math>u=0.7</math>]]
| [[image:Ikeda0900.png|thumb|left| <math>u=0.9</math>]]
|}

==挙動==
以下の図は、<math>u</math>の値を変化させた時の200個のランダムな点の動きを表したものである。各々の図の左上にある画像は[[アトラクター]]の形成過程を示し、右上の画像は中心部を拡大した画像となっている。

{| 
| [[image:ikeda sim u0.1.png|thumb|250px|left| u = 0.1]]
| [[image:ikeda sim u0.5.png|thumb|250px|left| u = 0.5]]
| [[image:ikeda sim u0.65.png|thumb|250px|left| u = 0.65]]
|-
| [[image:ikeda sim u0.7.png|thumb|250px|left| u = 0.7]]
| [[image:ikeda sim u0.8.png|thumb|250px|left| u = 0.8]]
| [[image:ikeda sim u0.85.png|thumb|250px|left| u = 0.85]]
|-
| [[image:ikeda sim u0.9.png|thumb|250px|left| u = 0.9]]
| [[image:ikeda sim u0.908.png|thumb|250px|left| u = 0.908]]
| [[image:ikeda sim u0.92.png|thumb|250px|right| u = 0.92]]
|-
|}

===Octave/MATLABコード===
このプロットを生成する[[GNU Octave|Octave]]/[[MATLAB]]コードを以下に示す。

<syntaxhighlight lang="matlab">

% u = ikeda parameter
% option = what to plot
%  'trajectory' - plot trajectory of random starting points
%  'limit' - plot the last few iterations of random starting points
function ikeda(u, option)
    P = 200;%how many starting points
    N = 1000;%how many iterations
    Nlimit = 20; %plot these many last points for 'limit' option

    x = randn(1,P)*10;%the random starting points
    y = randn(1,P)*10;

    for n=1:P,
        X = compute_ikeda_trajectory(u, x(n), y(n), N);

        switch option
        case 'trajectory' %plot the trajectories of a bunch of points
            plot_ikeda_trajectory(X);hold on;

        case 'limit'
            plot_limit(X, Nlimit); hold on;

        otherwise
            disp('Not implemented');
        end
    end

    axis tight; axis equal
    text(-25,-15,['u = ' num2str(u)]);
    text(-25,-18,['N = ' num2str(N) ' iterations']);
    end

    % Plot the last n points of the curve - to see end point or limit cycle
    function plot_limit(X,n)
    plot(X(end-n:end,1),X(end-n:end,2),'ko');
end

% Plot the whole trajectory
function plot_ikeda_trajectory(X)
    plot(X(:,1),X(:,2),'k');
    %hold on; plot(X(1,1),X(1,2),'bo','markerfacecolor','g'); hold off
end

%u is the ikeda parameter
%x,y is the starting point
%N is the number of iterations
function [X] = compute_ikeda_trajectory(u, x, y, N)
    X = zeros(N,2);
    X(1,:) = [x y];

    for n = 2:N

        t = 0.4 - 6/(1 + x^2 + y^2);
        x1 = 1 + u*(x*cos(t) - y*sin(t)) ;
        y1 = u*(x*sin(t) + y*cos(t)) ;
        x = x1;
        y = y1;

        X(n,:) = [x y];

    end
end

</syntaxhighlight>

==関連項目==
*[[池田研介]]

==脚注==
{{reflist}}

==外部リンク==
*[http://www.ike-dyn.ritsumei.ac.jp/ 非線形物理学研究室 / 池田研究室 - 立命館大学 理工学部 物理学科 非線形物理学研究室]

{{DEFAULTSORT:いけたしやそう}}
[[Category:数学に関する記事]]
[[Category:カオス力学系]]