池田写像

物理学や数学において、池田写像 (いけだしゃぞう、) は以下の複素解析写像で与えられる離散時間力学系である。
オリジナルの池田写像は非線形光共振器 (、非線形誘電体を含む) における光の軌跡のモデルとして池田研介により提案された。 
池田研介、大同寛明、秋元興一により上記の単純化した形に一般化された。
formula_2は共振器内の回転のn番目のステップにおける共振器内の電場を表し、formula_3とformula_4はそれぞれ、外部からのレーザー光線と共振器の線形位相を示すパラメータである。
特に、パラメータformula_5は振る舞いを特徴づける散逸パラメータと呼ばれ、formula_6の極限で池田写像は保存系になる。 
非線形誘電体の飽和効果を考慮に入れ、オリジナルの池田写像に修正を施した以下の形で使用される事が多い。
上記を実2次元平面に描く場合、
となる。ここで、"u"はパラメータであり、 
である。幾つかの値の"u"では、この系はカオスな振る舞いを示す。
アトラクター.
はパラメータformula_11を0.0から1.0まで0.01ずつ変化させた時の挙動を示している。シミュレーションに使用したステップ数は500、ランダムに選択した始点は20000個である。アトラクターを描くため、それぞれの軌道で20個の点がプロットされている。formula_11の値が増加するにつれ、アトラクターポイントの分岐は大きくなる。 
挙動.
以下の図は、formula_11の値を変化させた時の200個のランダムな点の動きを表したものである。各々の図の左上にある画像はアトラクターの形成過程を示し、右上の画像は中心部を拡大した画像となっている。
Octave/MATLABコード.
このプロットを生成するOctave/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