関数 > 求解と最適化 > 例: 1 階常微分方程式の初期値問題を解く
  
例: 1 階常微分方程式の初期値問題を解く
次の常微分方程式を解きます。
クリックしてこの式をコピー
クリックしてこの式をコピー
1. 初期値問題の条件を入力します。
クリックしてこの式をコピー
クリックしてこの式をコピー
クリックしてこの式をコピー
2. 解のパラメータ (解の区間の終点) を入力し、[t0, t1] での解の値の数を入力します。
クリックしてこの式をコピー
クリックしてこの式をコピー
Odesolve
ソルブブロックと odesolve関数を使用して微分方程式を解きます。
1. ソルブブロック内に y の導関数を定義します。
クリックしてこの式をコピー
2. z に対して y をプロットします。
クリックしてこの式をコピー
クリックしてこの式をコピー
3. odesolve 解をソルブブロック内でパラメトリック的に使用します。
クリックしてこの式をコピー
出力 fy(k) は、関数の関数であるため、解を求めるためにはパラメータの値を指定する必要があります。
4. この結果を、独立変数 t がない通常の関数名に割り当てます。
クリックしてこの式をコピー
クリックしてこの式をコピー
5. 2 本の曲線をプロットします。
クリックしてこの式をコピー
クリックしてこの式をコピー
Adams 関数、rkfixed 関数、Rkadapt 関数、Bulstoer 関数、Radau 関数
微分方程式を解く別の方法として、ODE ソルバー Adamsを使用できます。
1. ソルバパラメータ (初期値の解のベクトル) と導関数を定義します。
クリックしてこの式をコピー
クリックしてこの式をコピー
この導関数の 2 つ目の引数は未知関数の値のベクトルでなければなりません。
2. Adams 行列を計算します。
クリックしてこの式をコピー
関数 rkfixedRkadaptBulstoerまたは Radauをここで使用することもできます。
3. 独立変数の値 Y に対して解の関数の値 T をプロットします。
クリックしてこの式をコピー
クリックしてこの式をコピー
クリックしてこの式をコピー
odesolve の結果は、単一行関数ソルバからの結果を補間したものです。ソルブブロックバージョンでは、通常の表記法に従って問題を入力できますが、odesolve によって返される関数は、単一行ソルバによって返されるのと同じ表を補間したものです。