作业帮 > 综合 > 作业

写了一个MATLAB程序,可是除了B=0以外,其余取值都是无限死循环T

来源:学生作业帮 编辑:作业帮 分类:综合作业 时间:2024/07/19 06:06:22
写了一个MATLAB程序,可是除了B=0以外,其余取值都是无限死循环T
function [xp,yp,xm,ym,T]=euler6(h,B)
Vp=135;Vm=450;
xp(1)=0;
yp(1)=120;
xm(1)=0;
ym(1)=0;
d(1)=120;
T(1)=0;
A(1)=pi/2;
C(1)=A(1)-B;
k=1;
while d>0.1
k=k+1;
xp(k)=xp(k-1)+Vp*h*cos(C(k-1));
yp(k)=yp(k-1)+Vp*h*sin(C(k-1));
xm(k)=xm(k-1)+Vm*h*cos(A(k-1));
ym(k)=ym(k-1)+Vm*h*sin(A(k-1));
d(k)=sqrt((xp(k)-xm(k))^2+(yp(k)-ym(k))^2);
A(k)=asin((yp(k)-ym(k))/d(k));
C(k)=A(k)-B;
T(k)=T(k-1)+h;
end;
xp,yp,xm,ym,T
这个就是那个程序,当h=0.001,B=0时能得出正确答案,可是代入B的别的取值后就是无限死循环T T不知道到底是哪里的问题T T
粗略看了下,应该是while后面条件写错了,我估计你的意思应该是写成d(k)>0.1这样应该不会死循环了
再问: 不是这个的问题的说。。。。。因为。。。h=0.0001,B=0的时候算出来的答案是对的。。。顺便可不可以问下。。。怎么改才能求出B在区间【0,pi】上那个值的时候,T最大?