私のコードでは、ソリューション ベクトルのどこでもゼロが返されますが、その理由がわかりません。結合された 2 次 ODE を 4 つの 1 次 ODE に分解しました。
関数を xp.m として定義しています
function zprime = f(t,z)
a = 1;
b = 1;
c = 1.5;
zprime = zeros(4,1);
zprime(1) = z(2);
zprime(2) = -a*z(1) + b*(z(3) - z(1));
zprime(3) = z(4);
zprime(4) = -c*(z(3) - z(1));
end
次のコマンドを使用して、matlab で実行します。
[t,z] = ode45('xp',[1,100],[0 0 0 0]);
私の初期条件はすべて0なので、私の初期条件が0の解または何か他のものを与えるのですか? ic を変更すると、予想どおり、ソリューションが変更されます。
ありがとう