0

ローレンツシステムの最大のLyaponov指数を計算するためにこの手順ode45()を使用して、ダイナミクスのいくつかの作業にmatlabで使用しています。

このプロセスでは、から始まる微分方程式のシステムを解き、x0これをに非常に近い軌道と比較しx0ます。

ode45() 各タイムステップで、タイムステップを進める前に2番目の軌道を再調整する必要があるので、一度だけ呼び出すことができるようにしたいと思います-これは可能ですか？

代替ソリューションの開始の試みがここに示されていますが、機能しません。私が見る限り、以下で得られる行列rとR取得は類似しているはずです。

% Time 
ti = 0; tf = 1; res = 10;
T = linspace(ti, tf, res);

% Solve the system first time
[~,R] = ode45('tbmLorenz',T,x0);

% Alternate trajectory
r = zeros(size(R));

% Temporary cordinates
temp = zeros(3,3);

% Solve for second trajectory
for i = 2:(res-1)
    % Time step
    ts = T((i-1):(i+1));

    % Solve three steps
    [~,temp] = ode45('tbmLorenz',ts,r(i-1,:));

    r_i = temp(2,:)

    % Code to alter r_i goes here

    % Save this
    r(i,:) = r_i;
end

...しかし、そうではありません：

r =
    1.0000    3.0000    4.0000
    9.7011   20.6113    7.4741
   29.9265   16.4290   79.0449
   -5.7096  -15.2075   49.2946
  -12.4917  -13.6448   44.7003
  -13.6131  -13.8826   45.0346
  -13.5061  -13.1897   45.4782
  -13.0538  -13.0119   44.5473
  -13.4463  -13.8155   44.4783
         0         0         0
>> R
R =
    1.0000    3.0000    4.0000
    9.7011   20.6139    7.4701
   29.9663   16.5049   79.1628
   -5.7596  -15.2745   49.3982
  -12.4738  -13.5598   44.7800
  -13.5440  -13.8432   44.9084
  -13.5564  -13.3049   45.4568
  -13.1016  -12.9980   44.6882
  -13.3746  -13.7095   44.4364
  -13.7486  -13.6991   45.4092

の最後の行rがゼロであることは問題ではありません。

何か案は？乾杯！\ T

4

1 に答える 1

0

Ode45は、定義によりステップサイズが変化する適応アルゴリズムです。必要に応じてステップサイズを修正できますが、ode45を使用していない場合は、別のものを使用しています。

この種の問題を解決したい場合は、固定の時間ステップを使用することをお勧めします。このリンクには、必要なもののコードが含まれています。

ode45の使用を主張する場合は、シミュレーションファイル自体を時間補間入力で調整する必要があります。

于 2012-05-18T20:25:14.720 に答える