)+c2/s_t*(ptr(3,k)-ptr(3,k-1)); %计算出角b所对边和目标运动轨迹交点:(x1,y1,z1)%%%参见公式三%%%x(k)=pmr(1,k-1)+s_m/c1*(x1(k)-pmr(1,k-1)); y(k)=pmr(2,k-1)+s_m/c1*(y1(k)-pmr(2,k-1)); z(k)=pmr(3,k-1)+s_m/c1*(z1(k)-pmr(3,k-1)); %计算出导弹k时刻所运动到位置:(x,y,z)%%%参见公式三%%%pmr(:,k)=[x(k);y(k);z(k)]; r(k)=sqrt((ptr(1,k)-pmr(1,k))^2+(ptr(2,k)-pmr(2,k))^2+(ptr(3,k)-pmr(3,k))^2); ifr(k)<0.06; break; end; end sprintf('遭遇时间:%3.1f',0.1*k); figure(1); plot3(pmr(1,1:k),pmr(2,1:k),pmr(3,1:k),'k',ptr(1,:),ptr(2,:),ptr(3,:)); axis([02505025]); text(x(180),y(180),z(180),'\rightarrow百分比导引律制导下导弹运动轨迹'); text(ptr(1,280),ptr(2,280),ptr(3,280),'\rightarrow目标运动轨迹'); gridon 以后,鉴于程序中很多地方不结合模型图也极难了解,将其中关键图例和公式提取以下: 最终,程序运行过程分析完成,具体细节详见注释,运行结果以下图所表示: 总而言之,本工作对百分比导引法求解三维制导问题仿真程序进行了具体分析和注释,程序运行正常,期望对大家了解百分比导引法有所帮助。4月17日于试验室Createdbylylogn