常微分方程式.ppt

範例-3 上述範例畫三度空間之相位圖形 範例11-8：odeLorenz02.m 圖形中只出現一條曲線，此曲線代表以三個狀態變數為座標、以時間為參數的一條三度空間中的曲線 options = odeset(OutputFcn, odephas3); % 使用 odephas3 進行繪圖 ode45(lorenzOde, [0 10], [20 5 -5], options); 提示 要觀看 Lorenz 渾沌方程式隨時間而變的動畫，可在 MATLAB 指令視窗下直接輸入lorenz 指令 範例-4 (I) 假設 OutputFcn 設成“myfunc”： options = odeset(OutputFcn, myfunc) ODE 指令會呼叫 myfunc(tspan, y0, ‘init’) 讓 myfunc 進行各種初始化動作 積分步驟中，ODE 指令會持續呼叫status=myfunc(t, y) 若 status=1，則停止積分 積分結束時，ODE 指令會呼叫 myfunc([ ], [ ], ‘done’)，讓 myfunc 進行收尾動作 OutputSel 可指定要傳送到 OutputFcn 的狀態變數之元素 範例-4 (II) 只要傳送第一及第三個 Lorenz 渾沌方程式的狀態變數至 odeplot 範例11-9：