missileguidancesimudemo01.m

本程序包为导弹追踪导引法的相对弹道簇图MATLAB仿真计算程序。 关键词： 导弹 弹道 追踪导引 仿真 missile simulation guidance MATLAB

%function [ output_args ] = MissileGuidanceSimuDemo01( input_args )
%%
%%    导弹追踪导引法的相对弹道簇图MATLAB仿真计算程序
%%
%%    COPYRIGHT (C) vxLand , 2002-2008
%%

	clc; clear; close all;
	r0 = 1; p=2;	r0x = r0 /3;	
	q00=0.5*pi; dtor=(180/pi);
for q0=-0.0*pi:10/dtor:1.0*pi; % -2*pi:1:2*pi
	v1 = tan(q0/2);
	c = r0*sin(q0)/ v1^p;
	N = 1;
	for qx=q0:-0.01:0;
		q(N) = qx;
%		r(N) = c*(tan(qx/2))^p/sin(qx);
		r(N) = c*(sin(qx/2))^(p-1) / 2 / (cos(qx/2))^(p+1);
		d1 = pi - qx;
		mx(N) = r(N)*cos(d1);
		my(N) = -r(N)*sin(d1);
		N = N + 1;
	end
	if q0 >= pi
		plot(-mx,my,'m-','linewidth',2);
	else
		plot(mx,my,'m-','linewidth',2);hold on;
		plot(mx,-my,'m-','linewidth',2);
	end;
	title('追踪法导引的相对弹道簇图');
	grid on; set(gca,'YMinorGrid','off') , set(gca,'XMinorGrid','off');	set(gcf,'color','w');
	%xlabel('x / m');ylabel('y / m');%line([tx0 tx1],[0 0],'linewidth',2,'color','k');

	axis([-100 100 -100 100]);%line([-100 100],[0 0],'linewidth',2,'color','k');
	
end

	figure(1);
	axis([-r0 r0 -r0 r0]);line([-r0 r0],[0 0],'linewidth',1,'color','k');
% q(end)*57.3,r(end) %
%r,q*57.3,mx,my
	N1 = 1;
	for gama=0:10/57.3:2*pi
		rx(N1) = r0*cos(gama);
		ry(N1) = r0*sin(gama);
		N1 = N1 + 1;
	end
	plot(rx,ry,'ko-','linewidth',2);

	return