比例导引+弹道成型导引源程序

参考文献：

《几种增大空地导弹落角的制导方式比较》

《Strategic and tactical missile guidance》

clear
clc%-----------导弹参数---------
V_m=260;%飞行速度
X_m=0;
Y_m=1000; %初始飞行高度
theta_m=0*pi/180; %弹道倾角 %----------目标参数---------
V_t=0; %静止目标
X_t=5000;
Y_t=0;
theta_t=0;
dtheta_t=0;
n_t=0;
A_t=0;R=sqrt((X_m-X_t)^2+(Y_m-Y_t)^2);
q=atan((Y_t-Y_m)/(X_t-X_m));
dR=((X_m-X_t)*(V_m*cos(theta_m)-V_t*cos(theta_t))+(Y_m-Y_t)*(V_m*sin(theta_m)-V_t*sin(theta_t)))/sqrt((X_m-X_t)^2+(Y_m-Y_t)^2);
dq=((X_t-X_m)*(V_t*sin(theta_t)-V_m*sin(theta_m))-(Y_t-Y_m)*(V_t*cos(theta_t)-V_m*cos(theta_m)))/((X_m-X_t)^2+(Y_m-Y_t)^2);n_m=-q+theta_m;c=3;
qf=-90*pi/180;
g=9.8;n=1;
t=0;
dt=0.01;while (dR<0)if R>500t_go=R/abs(dR);Am=4*dq*abs(dR)+2*abs(dR)*(q-qf)/t_go+g*cos(theta_m); %考虑重力补偿因素的弹道成型制导律elseAm=4*dq*abs(dR); %比例导引enddtheta_m=Am/V_m; %纵向通道：弹道倾角变化函数theta_m=theta_m+dtheta_m*dt;%----------------------------计算坐标----------------------------X_m=X_m+V_m*cos(theta_m)*dt;Y_m=Y_m+V_m*sin(theta_m)*dt;alpha=Am/g/(0.3*g);n_m=-q+theta_m; R=sqrt((X_m-X_t)^2+(Y_m-Y_t)^2);q=atan((Y_t-Y_m)/(X_t-X_m));dR=((X_m-X_t)*(V_m*cos(theta_m)-V_t*cos(theta_t))+(Y_m-Y_t)*(V_m*sin(theta_m)-V_t*sin(theta_t)))/sqrt((X_m-X_t)^2+(Y_m-Y_t)^2);dq=((X_t-X_m)*(V_t*sin(theta_t)-V_m*sin(theta_m))-(Y_t-Y_m)*(V_t*cos(theta_t)-V_m*cos(theta_m)))/((X_m-X_t)^2+(Y_m-Y_t)^2);theta_m_store(n)=theta_m;  %保存弹道倾角 Am_store(n)=Am;   %保存纵向过载alpha_store(n)=alpha;%保存攻角P_m_store(:,n)=[X_m;Y_m]; %保存拦截弹坐标n=n+1;t=t+dt;
end
disp('脱靶量为（m）：')
R
disp('飞行时间为（s）：')
tfigure(1)
plot(P_m_store(1,:),P_m_store(2,:),X_t,Y_t,'r+')
hold on
xlabel('X/m')
ylabel('Y/m')figure(2)
plot((1:n-1)*dt,Am_store/g)
hold on
xlabel('time/s')
ylabel('Acceleration/g')
title('加速度')figure(3)
plot((1:n-1)*dt,theta_m_store*180/pi)
hold on
xlabel('time/s')
ylabel('\theta_m/°')
title('弹道倾角')figure(4)
plot((1:n-1)*dt,alpha_store)
hold on
xlabel('time/s')
ylabel('\alpha/°')
title('攻角')

比例导引+弹道成型导引源程序相关推荐

MATLAB 仿真 n年后的比例,基于Matlab的比例导引弹道仿真分析.PDF
基于Matlab的比例导引弹道仿真分析.PDF ·60· 战术导弹技术Ta(·ti(·alMissile TechnologyMay,2009,(3):60-64 [文章编号]1009．1300(20 ...
鱼雷导引仿真matlab
鱼雷导引弹道仿真与分析 1. 实验目的 2. 实验内容 3. 实验要求 4. 附加题代码下载 1. 实验目的通过鱼雷追踪目标的导引弹道仿真过程,初步掌握系统数学仿真方法. 2. 实验内容已知 ...
一种高超声速飞行器弹道的仿真方法
一种高超声速飞行器弹道的仿真方法目录一.相关研究二.基本思路三.巡航段弹道计算四.俯冲段弹道计算五.仿真结果参考文献一.相关研究临近空间高超声速飞行器(Near Space Hype ...
计算机模拟电子云密度,Intro2ComtSimu.ppt
Intro2ComtSimu.ppt 求职应注意的礼仪求职时最礼貌的修饰是淡妆面试时最关键的神情是郑重无论站还是坐,不能摇动和抖动对话时目光不能游弋不定要控制小动作不要为掩饰紧张情绪而散淡 ...
计算机网络笔记——概述、物理层、链路层（方老师408课程）
文章目录前言互联网概述互联网发展的三个阶段互联网标准化机构互联网的组成边缘部分的通信方式核心部分的交换方式我国计算机网络的发展计算机网络的类别计算机网络的性能速率.带宽.吞吐量 ...
定位技术课程设计-微信小程序校园导游系统
定位技术课程设计课程设计教学目的课程设计要求课程设计题目原题目拓展内容需求分析原理分析微信小程序API定位原理WIFI指纹定位原理路径规划算法调研详细设计总述主页面介绍学校简介页面介绍导引页面概述导引 ...
布点算法的原理和实现
在数据可视化的过程中,绘制网络拓扑图是很重要的,它能清晰呈现一个复杂网络的结构,节点的重要性和关系.比如下面几张图: 下面这张图是我的软件绘制的: 这些都有一个共同的问题,就是如何让图绘制的更加美观? ...
【计网】计算机网络期末总复习-谢希仁（个人总结）理论概念
文章目录计算机网络期末复习(个人总结) 第一章概述第二章.物理层第三章数据链路层第四章网络层第五章运输层第六章应用层其他补充算法和协议题目计算机网络期末复习(个人总结) ...
计算机网络「二」—— 物理层（多图详解）
本文自学计算机网络时所写笔记,网课为B站湖科大教书匠的计算机网络微课堂.(强烈安利这个课程,讲课思路条理清晰,PPT美轮美奂.通俗易懂) 本文为第二章笔记,陆续会更新余下内容计算机网络「一」计算机 ...

比例导引+弹道成型导引源程序

比例导引+弹道成型导引源程序相关推荐

最新文章

热门文章