我的代碼給了我零到處爲我的解決方案向量,但我不知道爲什麼。我已經將一個耦合的二階ODE分解爲四階一階ODE。MATLAB ODE解算器給0的地方
我有我的函數定義爲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解決方案還是別的什麼?當我改變集成電路時,解決方案會改變,正如預期的那樣。
感謝