分析ソリューション:
- 方程式を解く:
[x1、x2、x3、...] = resolve( 'eq1'、 'eq2'、 'eq3'、......、 'x1'、 'x2'、 'x3'、 '...。 ')
S = resolve( 'eq1'、 'eq2'、 'eq3'、......、 'x1'、 'x2'、 'x3'、 '....')
- 最初の方法では、ソルブの括弧内にx1、x2、およびx3の順序が必要であり(前の[]と一致)、2番目の方法には要件がありません。eqは方程式であり、式ではありません。
- 引用符を付けるか、引用符を付けないのが最善ですが、事前にsym、syms、その他の設定および事前定義された記号を使用する必要があります。
- 微分方程式を解きます。
dsolve( 'eq1'、 'eq2'、....)
同上
数値解法
通常、ode45(@ fun、tspan、y0)を使用します
1.サブ関数を記述し、行列の形式を使用して微分方程式を表します。元の問題はn個の1次微分方程式に変換する必要があります。
2.y0は通常(n * 1)ベクトルであり、tspanは通常バイナリベクトルです
2.関数ハンドル(無名関数)を使用することもできます
例:
y '= 2tの場合、微分方程式を解きます
tspan = [0,5]; y0 = 0; [t、y] = ode45(@(t、y)2 * t、tspan、y0);%関数ハンドルの使用方法