




版权说明:本文档由用户提供并上传,收益归属内容提供方,若内容存在侵权,请进行举报或认领
文档简介
弹道计算大作业目录一、初始条件和要求 )三、仿真结果3.1无控飞行弹道仿真根据无控弹道模型,写出s函数,搭建的仿真模块如下图所示:图1无控飞行仿真模块由于初始条件给定,因此模块没有输入;输出有六个,分别为导弹的射程变化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。模块的仿真时间由高度变化决定,当高度降为0(导弹落到地面上)时仿真结束。导出数据后画图如下: 图2无控飞行时各参数变化3.2平衡滑翔弹道仿真平衡滑翔弹道仿真模块如下图所示:图3平衡滑翔模块取仿真时间为150s,无输入,输出分别为:导弹的射程变化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。得到各参量时间变化图如下:图4平衡滑翔飞行时各参数变化3.3最大升阻比滑翔弹道仿真按最大升阻比飞行时弹道仿真模块如下图所示:图5最大升阻比飞行模块取仿真时间为180s,无输入,输出分别为:导弹的射程变化、高度变化、速度变化、弹道倾角变化、攻角变化以及密度变化。得到各参量时间变化图如下:图4最大升阻比飞行时各参数变化附录附表1无控弹道飞行时完整的s函数无控弹道function[sys,x0,str,ts,simStateCompliance]=trace2(t,x,u,flag)switchflag,case0,[sys,x0,str,ts,simStateCompliance]=mdlInitializeSizes;case1,sys=mdlDerivatives(t,x,u);case2,sys=mdlUpdate(t,x,u);case3,sys=mdlOutputs(t,x,u);case4,sys=mdlGetTimeOfNextVarHit(t,x,u);case9,sys=mdlTerminate(t,x,u);otherwiseDAStudio.error('Simulink:blocks:unhandledFlag',num2str(flag));endfunction[sys,x0,str,ts,simStateCompliance]=mdlInitializeSizessizes=simsizes;sizes.NumContStates=4;sizes.NumDiscStates=0;sizes.NumOutputs=5;sizes.NumInputs=0;sizes.DirFeedthrough=0;sizes.NumSampleTimes=1;sys=simsizes(sizes);x0=[0;2000;100;-5/180*pi];str=[];ts=[00];simStateCompliance='UnknownSimState';functionsys=mdlDerivatives(t,x,u)S=1.7;%参考面积,m^2AR=0.86;%展弦比e=0.9;%效率因子;m=115;%质量,kgg0=9.8;%海平面重力加速度,m/s^2Rd=6371000;%地球半径r=Rd/(Rd+x(2));g=g0*r^2;%飞行器所在高度的重力加速度rho0=1.225;%海平面大气密度,kg/m^3T0=288.15;rho=rho0*(1-0.0065*x(2)/T0)^4.2288;%飞行器所在高度的大气密度alpha=-x(4);%无控飞行时CLa=3.141592*AR/(1+sqrt(1+(AR/2)^2));%升力线斜率,/radCDo=0.02;%零升阻力系数epsilon=1/(pi*e*AR);%诱导阻力因子CL=CLa*alpha;%升力系数CD=CDo+epsilon*CL^2;%阻力系数X=CD*1/2*rho*x(3)^2*S;Y=CL*1/2*rho*x(3)^2*S;%以下为飞行器在铅垂平面的运动方程dx=x(3)*cos(x(4));dy=x(3)*sin(x(4));dv=-X/m-g*sin(x(4));dtheta=Y/(m*x(3))-g*cos(x(4))/x(3);sys=[dx;dy;dv;dtheta];functionsys=mdlUpdate(t,x,u)sys=[];functionsys=mdlOutputs(t,x,u)y1=x(1);y2=x(2);y3=x(3);rho0=1.225;T0=288.15;rho=rho0*(1-0.0065*x(2)/T0)^4.2288;sys=[x(1)x(2)x(3)x(4)rho];functionsys=mdlGetTimeOfNextVarHit(t,x,u)sampleTime=1;%Example,setthenexthittobeonesecondlater.sys=t+sampleTime;functionsys=mdlTerminate(t,x,u)sys=[];附表2平衡滑翔飞行部分代码平衡滑翔飞行functionsys=mdlDerivatives(t,x,u)S=1.7;%参考面积,m^2AR=0.86;%展弦比e=0.9;%效率因子;m=115;%质量,kgg0=9.8;%海平面重力加速度,m/s^2Rd=6371000;%地球半径r=Rd/(Rd+x(2));g=g0*r^2;%飞行器所在高度的重力加速度rho0=1.225;%海平面大气密度,kg/m^3T0=288.15;rho=rho0*(1-0.0065*x(2)/T0)^4.2288;%飞行器所在高度的大气密度CLa=3.141592*AR/(1+sqrt(1+(AR/2)^2));%升力线斜率,/radCDo=0.02;%零升阻力系数epsilon=1/(pi*e*AR);%诱导阻力因子alpha=2*m*g*cos(x(4))/(rho*x(3)^2*S*CLa);CL=CLa*alpha;%升力系数CD=CDo+epsilon*CL^2;%阻力系数X=CD*1/2*rho*sqrt(x(3)^2)*S;Y=CL*1/2*rho*sqrt(x(3)^2)*S;dx=x(3)*cos(x(4));dy=x(3)*sin(x(4));dv=-X/m-g*sin(x(4));dtheta=0;sys=[dx;dy;dv;dtheta];functionsys=mdlOutputs(t,x,u)rho0=1.225;T0=288.15;rho=rho0*(1-0.0065*x(2)/T0)^4.2288;S=1.7;%参考面积,m^2AR=0.86;%展弦比m=115;%质量,kgg0=9.8;%海平面重力加速度,m/s^2Rd=6371000;%地球半径r=Rd/(Rd+x(2));g=g0*r^2;%飞行器所在高度的重力加速度CLa=pi*AR/(1+sqrt(1+(AR/2)^2));%升力线斜率,/radalpha=2*m*g*cos(x(4))/(rho*x(3)^2*S*CLa);y(1)=x(1);y(2)=x(2);y(3)=x(3);y(4)=x(4);y(5)=alpha;y(6)=rho;sys=[y];附表3最大升阻比飞行部分代码最大升阻比飞行functionsys=mdlDerivatives(t,x,u)S=1.7;%参考面积,m^2AR=0.86;%展弦比e=0.9;%效率因子;m=115;%质量,kgg0=9.8;%海平面重力加速度,m/s^2Rd=6371000;%地球半径r=Rd/(Rd+x(2)); g=g0*r^2;%飞行器所在高度的重力加速度rho0=1.225;%海平面大气密度,kg/m^3T0=288.15;rho=rho0*(1-0.0065*x(2)/T0)^4.2288;%飞行器所在高度的大气密度CLa=3.141592*AR/(1+sqrt(1+(AR/2)^2));%升力线斜率,/radCDo=0.02;%零升阻力系数epsilon=1/(pi*e*AR);%诱导阻力因子alpha=sqrt(CDo*pi*AR*e)/CLa;CL=CLa*alpha;%升力系数CD=CDo+epsilon*CL^2;%阻力系数X=CD*1/2*rho*x(3)^2*S;Y=CL*1/2*rho*x(3)^2*S;%以下为飞行器在铅垂平面的运动方程dx=x(3)*cos(x(4));dy=x(3)*sin(x(4));dv=-X/m-g*sin(x(4));dtheta=Y/(m*x(3))-g*cos(x(4))/x(3);sys=[dx;dy;dv;dtheta];functionsys=mdlOutputs(t,x,u)rho0=1.225;T0=288.15
温馨提示
- 1. 本站所有资源如无特殊说明,都需要本地电脑安装OFFICE2007和PDF阅读器。图纸软件为CAD,CAXA,PROE,UG,SolidWorks等.压缩文件请下载最新的WinRAR软件解压。
- 2. 本站的文档不包含任何第三方提供的附件图纸等,如果需要附件,请联系上传者。文件的所有权益归上传用户所有。
- 3. 本站RAR压缩包中若带图纸,网页内容里面会有图纸预览,若没有图纸预览就没有图纸。
- 4. 未经权益所有人同意不得将文件中的内容挪作商业或盈利用途。
- 5. 人人文库网仅提供信息存储空间,仅对用户上传内容的表现方式做保护处理,对用户上传分享的文档内容本身不做任何修改或编辑,并不能对任何下载内容负责。
- 6. 下载文件中如有侵权或不适当内容,请与我们联系,我们立即纠正。
- 7. 本站不保证下载资源的准确性、安全性和完整性, 同时也不承担用户因使用这些下载资源对自己和他人造成任何形式的伤害或损失。
最新文档
- 2025年化学纤维加工丝合作协议书
- 企业办公设备采购与使用规范
- 新产品市场推广策略手册
- 知识产权登记的法律法规解读与实务操作
- 2025年海岛自动气象遥测系统项目合作计划书
- 幼儿故事狐狸和乌鸦深度解读
- 砖厂企业文化建设与绿色生产模式
- 西游记成长之道路观后感
- 2025年度电动三轮车租赁与绿色出行推广服务合同
- 电子商务平台内容营销策略探讨
- 学校食堂食品安全知识培训课件
- 合同-劳动主体变更三方协议
- 我国新闻短视频的创新模式及对策研究
- 阿米巴经营-稻盛和夫经营哲学的实学应用
- 八段锦口令标准版合集
- JCT414-2017 硅藻土的标准
- 新版药品管理法培训完整版本课件
- 北师大版高中英语选择性必修四全册课文及翻译(中英文Word)
- 心灵手巧织美好-精美挂件我会编 (教案)-二年级下册劳动浙教版
- IPO项目律师尽职调查查验计划表模版
- 医院信息系统HIS知识培训教学课件-HIS的主要内容
评论
0/150
提交评论