Octave 求解微分方程

Question

我如何用 y(0)=24 求解微分方程 y'+y=t？

我需要用 .m 格式的文件定義微分方程嗎？

Answer 1

要求解常微分方程，您可以使用 lsode 函數（運行 lsode 尋求幫助）。

f = @(y,t) t-y;
t = linspace(0,5,50)';
y=lsode(f, 24, t);
plot(t,y);
print -djpg figure-lsnode.jpg

Answer 2

dy/dt + y = t  

在兩側乘以積分因子e^t 。 然后，

d/dt [(e^t)y] = te^t
y = t-1+ce^(-t)  

因為， y(0)=24那么c=25

y = t-1+25e^(-t)