用C或MATLAB编个程序

来源:百度知道 编辑:UC知道 时间:2024/07/02 08:16:29
已知线性微分方程为 dy(t)/dt + y(t)=u(t),y(0)=1,用C语言或MATLAB语言编写一段程序求y(1000)并输出结果,要求对程序做简要注释

用dsolve(),参数为:方程,初始值和求解的变量:
y=dsolve('Dy+y-u(t)','y(0)=1','t')
例如:u(t)=t,
y=dsolve('Dy+y-t','y(0)=1','t')
y=
t-1+2*exp(-t)
y(1000)~=999 。

D1 = dsolve('Dy + y = u(t)')
matlab中这一句代码足矣,只是这个u(t)是什么具体东西,比如u(t)=t就好办多了。