微分方程式の解決において、オクターブ/マットラボでイベントを実装する方法を理解するのにいくつかの問題があります。
たとえば、微分方程式 y' = -y を解く次の単純なコードを考えてみましょう。
function dy = odefun(t,y)
dy = -y;
endfunction
options = odeset('RelTol',1e-4,'AbsTol',[1e-4]);
[T,Y] = ode45(@odefun,[0 12],[1],options);
今回はランダムイベントを紹介したいと思います。たとえば、一定の時間ステップで、値をランダムに変更yし、微分方程式に従って進化を続けたいと思います。これどうやってするの？