作图：建立程序tuxing41.m如下： ezplot('t + 1/exp(t)',[0,1]) title('t + 1/exp(t)') （2）先建立M文件fun3.m，如下： function dy=fun3(t,y) dy=zeros(1,1); dy(1)=-y(1)+t+1; 再取t0=0,tf=1,建立程序tuxing42.m如下： t0=0;tf=1; [T,Y]=ode45('fun3',[0 1],[1]); plot(T,Y,'ro'); title('比较图'); t=0:0.1:1; y=t+1./exp(t); hold on plot(t,y,'b'); 四. 上机结果 题1结果为： ans = C4*exp(2*x) + C2*exp(x*(5^(1/2)/2 - 1/2)) + C3/exp(x*(5^(1/2)/2 + 1/2)) 题2结果为： x = 4*cos(t) - 2/exp(2*t) + 3*sin(t) - (2*sin(t))/exp(t) y = sin(t) - 2*cos(t) + (2*cos(t))/exp(t) 题3结果为： 题4结果为： 解析解为： y = t + 1/exp(t)
微分方程(组)的解析解和数值解
作图：建立程序tuxing41.m如下：ezplot('t+1/exp(t)',[0,1])title('t+1/exp(t)')（2）先建立M文件fun3.m，如下：functiondy=fun3(t,y)dy=zeros(1,1);dy(1)=-y(1)+t+1;再取t0=0,tf=1,建立程序tuxing42.m如下：t0=0;tf=1;[T,Y
推荐度：
点击下载文档文档为doc格式
