你好同学,为你解答如下:
clc;clear
%(a)
odefun = @(t,y) [y(2); -10*y(2)-16*y(1)+16];
your_candidate_number = 202101; % 自己根据自己的学号改
[t,y] = ode45(odefun, [0, your_candidate_number/25000], [0;0]);
%(b)
maxy = max(y(:,1))
maxt = t(y(:,1)==maxy)
maxdy = max(y(:,2))
maxdt = t(y(:,2)==maxdy)
%(c)
figure(1);clf
plot(t,y(:,1)); xlabel('t');ylabel('y');
%(d)
figure(2);clf
plot(t,y(:,2)); xlabel('t');ylabel('dy/dt');
结果:
maxy =
0.999999903803806
maxt =
8.084040000000000
maxdy =
1.258903906957620
maxdt =
0.220406349778054
如有帮助还望题主给个宝贵的采纳支持一下答主答题哟