注册 登录
编程论坛 Matlab

导弹追踪问题

剑八 发布于 2008-11-24 12:56, 2420 次点击
clear;clf
h=0.001;%时间步长
k=1;
t(1)=0;x(1)=0;y(1)=0;%初始值
while x<=30
x(k+1)=x(k)+400*h/sqrt(1+((90*t(k)-y(k))/(30-x(k)))^2)+eps;
y(k+1)=y(k)+400*h/sqrt(1+((30-x(k))/(90*t(k)-y(k)))^2)+eps;
t(k+1)=h*k;
k=k+1;
endfor i=1:30000
      plot(x(1:i),y(1:i),'r-',[-1 32],[0 0],'k',[30 30],[0 y(i)],'b.-');
      text(4,7,'导弹追踪'); axis equal; hold on;
      m(i)=getframe;
      hold off;
     end
movie(m);
   导弹追踪问题   运行时
>> clear
>> missile
Warning: Divide by zero.
> In missile at 8
??? Index exceeds matrix dimensions.

Error in ==> missile at 15
      plot(x(1:i),y(1:i),'r-',[-1 32],[0 0],'k',[30 30],[0 y(i)],'b.-');
请问 这是哪里出错了  还有运行时 坐标会动
4 回复
#2
lfphit2008-11-24 13:26
> size(x)

ans =

     1    80
i=30000>80
i超过了x的维数,取i=80就可以了
#3
lfphit2008-11-24 13:38
提示很清楚
#4
剑八2008-11-24 15:16
谢谢  再问一下  
我方战舰发现敌舰位于正东30km,并以90km/h的速度向正北方向行驶,我方迅即向敌舰发射导弹,导弹速度是400km/h,且方向始终对着敌舰。
  这题可以得到解析解吗  
dsolve('(Dx)^2+(Dy)^2=400^2','Dy=Dx*(90*t-y)/(30-x)','x(0)=0','y(0)=0')
  这样不对   哪里出错了
#5
huilichuan2008-11-25 08:15
怎么好像是个圆?另外应该用不到dsolve
1