%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第七章例7－1－3程序
% 導彈追蹤運動目標的軌跡
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
%global vt vm
vt=input('vt=(書上取500) ');	%輸入主程序及函數程序共用的參數
vm=input('vm=(書上取1000或800) ');	
z0=input('[x0;y0]=(書上取[3000;4000]) ');				%輸入數值積分函數需要的參數
tspan=input('tspan=[t0,tfinal]=(書上取[0,4.5]) ');		%輸入數值積分函數需要的參數
[t,z] = ode23('exn713f',tspan,z0);		%進行數值積分
plot(z(:,1),z(:,2));grid					%繪圖
% 在慣性坐標中，M點位置的導數是相對速度，而其二次導數則為絕對加速度
dt=diff(t); Ldt=length(dt);				%為了求導數先求t的增量
x=z(:,1);y=z(:,2);							% 把z寫成x,y兩個分量形式
vx=diff(z(:,1))./dt;vy=diff(z(:,2))./dt;	% 對z求導數,注意每導數一次序列長度少一
wx=diff(vx)./dt(1:Ldt-1);wy=diff(vy)./dt(1:Ldt-1); % 對z求二次導數
[t(2:Ldt),x(2:Ldt),y(2:Ldt),wx,wy]					% 顯示數據
set(gcf,'color','w')