哪位高手有matlab弹道仿真程序
来自：络知识 更新日期：早些时候
求助关于用matlab进行弹道仿真计算的程序~
导弹仿真Matlab源代码(missile simulation)
http://www.matlabsky.com/thread-284-1-1.html
第二例子
%by dynamic
%see also http://www.matlabsky.com
%2009.2.18
%比例导引理想弹道仿真Matlab代码
%给出了比例导引法的差分方程, 采用Matlab语言, 对比例导引法的导弹理想弹道进行了三维数据仿真, 绘制出三维理想弹道。随导弹与目标参数变化情况, 对理想弹道的特性进行了分析。
clear all
clc
clear
tt=0.1;
sm=0.6*tt;
st=0.42*tt;
x(1)=0;y(1)=0;z(1)=0;
pmr(:,1)=[x(1);y(1);z(1)];
ptr(:,1)=[25;5;10];
m=3;
q(1)=0;
o(1)=0;
a(1)=0;
for(k=2:600)
ptr(:,k)=[25-0.42*cos(pi/6)*tt*k;5;10+0.42*sin(pi/6)*k*tt];
r(k-1)=sqrt((ptr(1,k-1)-pmr(1,k-1))^2+(ptr(2,k-1)-pmr(2,k-1))^2+(ptr(3,k-1)-pmr(3,k-1))^2);
c=sqrt((ptr(1,k)-pmr(1,k-1))^2+(ptr(2,k)-pmr(2,k-1))^2+(ptr(3,k)-pmr(3,k-1))^2);
b=acos((r(k-1)^2+st^2-c^2)/(2*r(k-1)*st));
dq=acos((r(k-1)^2-st^2+c^2)/(2*r(k-1)*c));
if abs(imag(b))>0
b=0.0000001;
end
if abs(imag(dq))>0
dq=0.0000001;
end
q(k)=q(k-1)+dq;
o(k)=o(k-1)+m*dq;
a(k)=o(k)-q(k);
c1=r(k-1)*sin(b)/sin(a(k)+b);
c2=r(k-1)*sin(a(k))/sin(a(k)+b);
c3=sqrt((c1-sm)^2+(c2-st)^2+2*(c1-sm)*(c2-st)*cos(a(k)+b));
dq=a(k)-acos(((c1-sm)^2+c3^2-(c2-st)^2)/(2*(c1-sm)*c3));
if abs(imag(dq))>0
dq=0.0000001;
end
q(k)=q(k-1)+dq;
o(k)=o(k-1)+m*dq;
a(k)=o(k)-q(k);
c1=r(k-1)*sin(b)/sin(a(k)+b);
c2=r(k-1)*sin(a(k))/sin(a(k)+b);
c3=sqrt((c1-sm)^2+(c2-st)^2+2*(c1-sm)*(c2-st)*cos(a(k)+b));
dq=a(k)-acos(((c1-sm)^2+c3^2-(c2-st)^2)/(2*(c1-sm)*c3));
if abs(imag(dq))>0
dq=0.0000001;
end
q(k)=q(k-1)+dq;
o(k)=o(k-1)+m*dq;
a(k)=o(k)-q(k);
c1=r(k-1)*sin(b)/sin(a(k)+b);
c2=r(k-1)*sin(a(k))/sin(a(k)+b);
c3=sqrt((c1-sm)^2+(c2-st)^2+2*(c1-sm)*(c2-st)*cos(a(k)+b));
x1(k)=ptr(1,k-1)+c2/st*(ptr(1,k)-ptr(1,k-1));
y1(k)=ptr(2,k-1)+c2/st*(ptr(2,k)-ptr(2,k-1));
z1(k)=ptr(3,k-1)+c2/st*(ptr(3,k)-ptr(3,k-1));
x(k)=pmr(1,k-1)+sm/c1*(x1(k)-pmr(1,k-1));
y(k)=pmr(2,k-1)+sm/c1*(y1(k)-pmr(2,k-1));
z(k)=pmr(3,k-1)+sm/c1*(z1(k)-pmr(3,k-1));
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);
if r(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([0 25 0 5 0 25]);
text(x(80),y(80),z(80),'\leftarrow 比例导引');
grid on
哪位高手有matlab弹道仿真程序视频
相关评论：13420485303：哪位高手有matlab弹道仿真程序
东叙邱导弹仿真Matlab源代码(missile simulation)http:\/\/www.matlabsky.com\/thread-284-1-1.html 第二例子 by dynamic see also http:\/\/www.matlabsky.com 2009.2.18 比例导引理想弹道仿真Matlab代码 给出了比例导引法的差分方程, 采用Matlab语言, 对比例导引法的导弹理想弹道进行了三维数据仿真, 绘制出...
13420485303：matlab弹道仿真问题
东叙邱3、初始条件大多数都没给，包括αβUVW等，r勉强算给了，但含义不明。
13420485303：matlab 内弹道计算问题
东叙邱function ndd%59nian130A=0.87; %枪(炮)膛横断面积A dm^2G=19;%33.4; %弹重 kgW0=2.04; %药室容积 dm^3l_g=25.0; %身管行程 dmP_0 =30000; %起动压力 kpafai1=1.02; %次要功系数K=1.03; %运动阻力系数φ1theta =0.2; %火药热力系数%==...
13420485303：中国兵器工业试验测试研究院待遇
东叙邱1、熟悉各类火工品、作战平台、无人机系统、气象测试系统，具备内弹道建模仿真、内弹道设计、弹药总体设计等能力并能熟练运用相关设计软件（Matlab、UG\/solidworks\/AutoCAD、Ansys等）；2、具备高度工作责任心，较好组织协调和管理能力。职位职责：1、主要负责主持武器装备试验的策划、实施、数据处理与评估等；...
13420485303：导引系数对弹道特性的影响
东叙邱主要研究内容如下:(1)建立无控尾翼火箭弹及固定鸭舵火箭弹物理模型,进行尾翼火箭弹的外流场仿真及气动特性分析,重点研究了固定鸭舵火箭弹的气动特性。结果表明:尾翼火箭弹气动特性符合弹箭气动规律;增大固定鸭舵舵偏角有助于增大火箭弹升力、俯仰力矩,同时会增大全弹空气阻力;当弹体长径比达到某一临界值时鸭舵...
13420485303：飞行器制导与控制及其MATLAB仿真技术目录
东叙邱第6章关注导引律设计，探讨了各种导引法则，以及如何通过MATLAB进行仿真和精度统计分析。最后，第7章详细介绍了飞行器制导与控制弹道的仿真过程，包括目标模型、飞行器模型、制导回路和控制系统仿真，以及控制力和力矩的仿真，为实际工程提供了全面的仿真平台。
13420485303：求解一个matlab弹丸质点外弹道的微分方程组
东叙邱弹丸质点外弹道 计算结果:弹丸落地时间:106.4184 弹丸飞行水平距离:5.5450e+004
13420485303：正侧式合成孔径雷达 (SAR) 成像
东叙邱以实现方位聚焦。6.2 Back-Projection Algorithm时域反投影算法则用于处理最终的成像，通过限制处理区域，提高计算效率。总结与讨论虽然斜视SAR和前视SAR是研究热点，弹载SAR面临更多挑战，如曲线弹道成像、宽域成像、大斜视成像以及目标定位。这些技术要求结合SAR的特性，解决特定环境下的复杂问题。
13420485303：95式狙击步枪垂直向上朝天空开枪,子弹最高能达到多少米?
东叙邱首先更正一下，95式枪族里面的狙击步枪叫85式狙击步枪，85是开始开发的年份，95是定型的年份。第二，世界上还没有把狙击步枪的垂直射高作为数据，所以只能根据子弹初速计算了。85式狙击步枪初速：830米\/秒，重力加速度g=9.80米／秒，在忽略空气阻力的情况下，子弹到速度为零时共用了84.7秒，则飞行...
13420485303：北京航空航天大学宇航学院的考的专业课是什么
东叙邱研究方向名称:不区分研究方向 专业备注:全日制专业学位，学制2.5年 考试科目单元 考试科目代码 考试科目名称 第一门考试科目 101 思想政治理论 第二门考试科目 201 英语一 第三门考试科目 301 数学一 第四门考试科目 931 自动控制原理综合 或 951 力学基础 其他方向见http:\/\/yzb.buaa.edu.cn\/ ...
您好，很高兴为您解答问题
一种是编计算程序，参照故障分析的书；
一种是采用仿真模型做，直接搭建模型，推荐从网上下载PSAT开源元件直接做。
如果对回答满意，请点【采纳答案】，如果还有问题，请点【追问】
希望我的回答对您有所帮助，希望能采纳。
d3.json("uk.json", function(error, uk)
if (error) return console.error(error);
svg.append("path")
.datum(topojson.feature(uk, uk.objects.subunits))
.attr("d", d3.geo.path().projection(d3.geo.mercator()));
});
导弹仿真Matlab源代码(missile simulation)
http://www.matlabsky.com/thread-284-1-1.html
第二例子
%by dynamic
%see also http://www.matlabsky.com
%2009.2.18
%比例导引理想弹道仿真Matlab代码
%给出了比例导引法的差分方程, 采用Matlab语言, 对比例导引法的导弹理想弹道进行了三维数据仿真, 绘制出三维理想弹道。随导弹与目标参数变化情况, 对理想弹道的特性进行了分析。
clear all
clc
clear
tt=0.1;
sm=0.6*tt;
st=0.42*tt;
x(1)=0;y(1)=0;z(1)=0;
pmr(:,1)=[x(1);y(1);z(1)];
ptr(:,1)=[25;5;10];
m=3;
q(1)=0;
o(1)=0;
a(1)=0;
for(k=2:600)
ptr(:,k)=[25-0.42*cos(pi/6)*tt*k;5;10+0.42*sin(pi/6)*k*tt];
r(k-1)=sqrt((ptr(1,k-1)-pmr(1,k-1))^2+(ptr(2,k-1)-pmr(2,k-1))^2+(ptr(3,k-1)-pmr(3,k-1))^2);
c=sqrt((ptr(1,k)-pmr(1,k-1))^2+(ptr(2,k)-pmr(2,k-1))^2+(ptr(3,k)-pmr(3,k-1))^2);
b=acos((r(k-1)^2+st^2-c^2)/(2*r(k-1)*st));
dq=acos((r(k-1)^2-st^2+c^2)/(2*r(k-1)*c));
if abs(imag(b))>0
b=0.0000001;
end
if abs(imag(dq))>0
dq=0.0000001;
end
q(k)=q(k-1)+dq;
o(k)=o(k-1)+m*dq;
a(k)=o(k)-q(k);
c1=r(k-1)*sin(b)/sin(a(k)+b);
c2=r(k-1)*sin(a(k))/sin(a(k)+b);
c3=sqrt((c1-sm)^2+(c2-st)^2+2*(c1-sm)*(c2-st)*cos(a(k)+b));
dq=a(k)-acos(((c1-sm)^2+c3^2-(c2-st)^2)/(2*(c1-sm)*c3));
if abs(imag(dq))>0
dq=0.0000001;
end
q(k)=q(k-1)+dq;
o(k)=o(k-1)+m*dq;
a(k)=o(k)-q(k);
c1=r(k-1)*sin(b)/sin(a(k)+b);
c2=r(k-1)*sin(a(k))/sin(a(k)+b);
c3=sqrt((c1-sm)^2+(c2-st)^2+2*(c1-sm)*(c2-st)*cos(a(k)+b));
dq=a(k)-acos(((c1-sm)^2+c3^2-(c2-st)^2)/(2*(c1-sm)*c3));
if abs(imag(dq))>0
dq=0.0000001;
end
q(k)=q(k-1)+dq;
o(k)=o(k-1)+m*dq;
a(k)=o(k)-q(k);
c1=r(k-1)*sin(b)/sin(a(k)+b);
c2=r(k-1)*sin(a(k))/sin(a(k)+b);
c3=sqrt((c1-sm)^2+(c2-st)^2+2*(c1-sm)*(c2-st)*cos(a(k)+b));
x1(k)=ptr(1,k-1)+c2/st*(ptr(1,k)-ptr(1,k-1));
y1(k)=ptr(2,k-1)+c2/st*(ptr(2,k)-ptr(2,k-1));
z1(k)=ptr(3,k-1)+c2/st*(ptr(3,k)-ptr(3,k-1));
x(k)=pmr(1,k-1)+sm/c1*(x1(k)-pmr(1,k-1));
y(k)=pmr(2,k-1)+sm/c1*(y1(k)-pmr(2,k-1));
z(k)=pmr(3,k-1)+sm/c1*(z1(k)-pmr(3,k-1));
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);
if r(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([0 25 0 5 0 25]);
text(x(80),y(80),z(80),'\leftarrow 比例导引');
grid on
哪位高手有matlab弹道仿真程序视频
相关评论：
东叙邱导弹仿真Matlab源代码(missile simulation)http:\/\/www.matlabsky.com\/thread-284-1-1.html 第二例子 by dynamic see also http:\/\/www.matlabsky.com 2009.2.18 比例导引理想弹道仿真Matlab代码 给出了比例导引法的差分方程, 采用Matlab语言, 对比例导引法的导弹理想弹道进行了三维数据仿真, 绘制出...
东叙邱3、初始条件大多数都没给，包括αβUVW等，r勉强算给了，但含义不明。
东叙邱function ndd%59nian130A=0.87; %枪(炮)膛横断面积A dm^2G=19;%33.4; %弹重 kgW0=2.04; %药室容积 dm^3l_g=25.0; %身管行程 dmP_0 =30000; %起动压力 kpafai1=1.02; %次要功系数K=1.03; %运动阻力系数φ1theta =0.2; %火药热力系数%==...
东叙邱1、熟悉各类火工品、作战平台、无人机系统、气象测试系统，具备内弹道建模仿真、内弹道设计、弹药总体设计等能力并能熟练运用相关设计软件（Matlab、UG\/solidworks\/AutoCAD、Ansys等）；2、具备高度工作责任心，较好组织协调和管理能力。职位职责：1、主要负责主持武器装备试验的策划、实施、数据处理与评估等；...
东叙邱主要研究内容如下:(1)建立无控尾翼火箭弹及固定鸭舵火箭弹物理模型,进行尾翼火箭弹的外流场仿真及气动特性分析,重点研究了固定鸭舵火箭弹的气动特性。结果表明:尾翼火箭弹气动特性符合弹箭气动规律;增大固定鸭舵舵偏角有助于增大火箭弹升力、俯仰力矩,同时会增大全弹空气阻力;当弹体长径比达到某一临界值时鸭舵...
东叙邱第6章关注导引律设计，探讨了各种导引法则，以及如何通过MATLAB进行仿真和精度统计分析。最后，第7章详细介绍了飞行器制导与控制弹道的仿真过程，包括目标模型、飞行器模型、制导回路和控制系统仿真，以及控制力和力矩的仿真，为实际工程提供了全面的仿真平台。
东叙邱弹丸质点外弹道 计算结果:弹丸落地时间:106.4184 弹丸飞行水平距离:5.5450e+004
东叙邱以实现方位聚焦。6.2 Back-Projection Algorithm时域反投影算法则用于处理最终的成像，通过限制处理区域，提高计算效率。总结与讨论虽然斜视SAR和前视SAR是研究热点，弹载SAR面临更多挑战，如曲线弹道成像、宽域成像、大斜视成像以及目标定位。这些技术要求结合SAR的特性，解决特定环境下的复杂问题。
东叙邱首先更正一下，95式枪族里面的狙击步枪叫85式狙击步枪，85是开始开发的年份，95是定型的年份。第二，世界上还没有把狙击步枪的垂直射高作为数据，所以只能根据子弹初速计算了。85式狙击步枪初速：830米\/秒，重力加速度g=9.80米／秒，在忽略空气阻力的情况下，子弹到速度为零时共用了84.7秒，则飞行...
东叙邱研究方向名称:不区分研究方向 专业备注:全日制专业学位，学制2.5年 考试科目单元 考试科目代码 考试科目名称 第一门考试科目 101 思想政治理论 第二门考试科目 201 英语一 第三门考试科目 301 数学一 第四门考试科目 931 自动控制原理综合 或 951 力学基础 其他方向见http:\/\/yzb.buaa.edu.cn\/ ...
