用matlab对微分方程求解实验报告(2)

利用二分法更改tf tf=5时

6

tf=2.5时

7

tf=3.15时：

8

所以在t=3.15时刻恰好追上！ W=5 M文件代码

function dy=eq4(t,y) dy=zeros(2,1);

dy(1)=5*(10+20*cos(t)-y(1))/sqrt((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2);

dy(2)=5*(20+15*sin(t)-y(2))/sqrt((10+20*cos(t)-y(1))^2+(20+15*sin(t)-y(2))^2); 命令：t0=0;tf=10;

9

[t,y]=ode45('eq4',[t0 tf],[0 0]); T=0:0.1:2*pi; X=10+20*cos(T); Y=20+15*sin(T); plot(X,Y,'-') hold on

plot(y(:,1),y(:,2),'*') 运行结果

更改tf=20 运行结果

10