生態系のバランスや、バタフライエフェクトの問題をODE45で解く：無償ソフトで技術計算しよう【シミュレーション応用編】（1）（2/2 ページ）

「微分方程式を解く」といえば、無味乾燥なイメージ。そこで今回は、生態系のバランスや、気象など私たちの生活に関わっているような数式や問題を取り上げる。

[伊藤孝宏，MONOist] PC用表示関連情報

LINE

Hatena

前のページへ | 　　　　　　

ローレンツ方程式

　大気変動を研究していた気象学者のエドワード・ローレンツは、初期値のわずかな違いが気象の予測結果を大きく変えることに気付きました。これが、バタフライ効果と呼ばれているもので、下記のローレンツ方程式が、カオス研究の先駆けとなった式です。

　では、ローレンツ方程式をFreeMatで解いてみます。下記のlorentz.mを保存し、コマンドウィンドウでlorentz（x初期値、y初期値、z初期値）と入力すると、図3に示すような結果が得られます。図3は蝶が羽根を広げたような形状で、ローレンツ・アトラクタと呼ばれています。

　初期値を変えて計算してみると分かりますが、わずかな違いでアトラクタは形状を複雑に変化させます。参考までに、初期値を原点周りで半径0.1の円周上で変化させた場合のアトラクタの形状変化を動画にしたものを図4に添付しましたので、ご覧ください。

function lorentz(x0,y0,z0)
    xyz=@(t,y) [-10*y(1)+10*y(2);-y(1)*y(3)+28*y(1)-y(2);y(1)*y(2)-8/3*y(3)];
    span=[0,50];
    t0=[x0;y0;z0];
    [t,y]=ode45(xyz,span,[x0;y0;z0]);
    plot3(y(:,1),y(:,2),y(:,3));
    view(3);grid('on');
    title(['X0=',num2str(x0),' Y0=',num2str(y0),' Z0=',num2str(z0)]);
    xlabel('x');ylabel('y');zlabel('z');
    axis([-30,30,-30,30,0,50]);

lorentz.m

＞＞「lorentz.m」ダウンロード