八度解决微分方程

Question

如何用y（0）= 24解决微分方程y＆＃39; + 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)