Question

这货编的一个MATLAB程序一直死循环T T急求帮助

function [xp,yp,xm,ym]=euler4(h)
Vp=135;Vm=450;
xp(1)=0;
yp(1)=120;
xm(1)=0;
ym(1)=0;
d(1)=120;
T(1)=0;
k=1;
while d>0.1
k=k+1;
xp(k)=xp(k-1)+Vp*(yp(k-1)-ym(k-1))*h/d(k-1);
yp(k)=yp(k-1)-Vp*(xp(k-1)-xm(k-1))*h/d(k-1);
xm(k)=xm(k-1)+Vm*(xp(k-1)-xm(k-1))*h/d(k-1);
ym(k)=ym(k-1)+Vm*(yp(k-1)-ym(k-1))*h/d(k-1);
d(k)=sqrt((xp(k)-xm(k))^2+(yp(k)-ym(k))^2);
T(k)=T(k-1)+h;
end;
end

Answer 1

有可能是你的h取得太大，导致，上一步未相遇，下一步已经走远
未能满足小于一定距离

或许可以让判断条件定为距离差是否在变小

如果在变小说明还在靠近
一旦距离差变大，说明这一步已经让两者又重新远离
这时候循环就该结束了
那样就不会死循还

Answer 2

也就是d始终大于0.1，可以修改一下阈值判定 比如0.1改为0.01什么的

追问

d如果始终大于0.1，那就肯定始终大于0.01啊。。。。

追答

看错了。。。改成0.2 什么的，反正是阈值的问题

追问

Q Q我试了一下，只有d>18才能算出东西来。。。可是。。。相遇问题= =d>18要怎么破？！

追答

不知道你死循环持续了多久，建议
d(k)=sqrt((xp(k)-xm(k))^2+(yp(k)-ym(k))^2);

这句末尾的分号去掉，看看d每次循环的数据