末制导大作业答案V3.0

1、0=丄Vma一gcos0=vyVm-1Vcos0Lasiny一acosyVvzvyaVcos0mx=Vcos0cos屮mmVy=Vsin0mmz=-Vcos0sin屮mmVV为导弹速度，0弹道倾角，屮弹道偏角，y速度滚转角,mVV弹位移，X为气动力阻力，具体表达式为x,y,z为导mmmX=CqSxrefa和a分别为导弹在弹道坐标系内的法向加速度和侧向加速度，vyvz弹对过载指令的响应加速度。该加速度为导2三维空间弹目相对运动方程弹目相对运动方程如下：末制导大作业答案1运动学和动力学方程为了验证制导律，首先简化设计思路，暂不考虑导弹绕质心的转动，不考虑姿态控制回路，暂不考虑加速度指令的跟踪时间常

2、数，同时假设攻角和侧滑角为小量，同时假设速度滚转角为0，得运动学和动力学方程如下：V=(PcosacosP-X一Gsin0)/m=(P一X)/m一gsin0macosy+asiny一gcos0vyVvzVe_VV-vxxtxme=V一VvyytymeVVvzzyzm弹目相对距离为e=Je2+e2+e2r*xyz弹目相对速度为V=pe2+e2+e2c耳vxvyvze俯仰视线角九二arcsin4Der偏航视线角九二Te俯仰视线角九De(e2+e2)一e(ee+ee)vyxz.yxvxz_vze2-e2+e2rxz偏航视线角九Te-e-e-exvzzvxe2+e2xz3制导律根据比例导引理论，设计比

3、例导引律nN-V-q+N-cos九-gcypycegyDn一N-V-qczpzcb其中，n、n为加速度指令，当忽略自动驾驶仪动态特性和指令饱和时，有cycza=n和a=n；q、q为导引头测量的俯仰、偏航视线角速度，当忽略导引vycyvzczeb头动态测量特性时，有q=九和q=九，N为俯仰比例导引系数，N为重力eDbTpygy补偿系数，N为偏航比例导引系数。pz4仿真条件仿真初始条件m1.06kg005,1001)屮0。v0a卩。000V20m/s0 x=y=z=0TOC o 1-5 h z2)m0m0m0=V=V=0mx0my0mz0y=z=0t0t0 x=1500m,2000mt0=V=V=

4、0 xt0yt0zt0阻力系数马赫数0.30.60.9Cx0.5018430.5060340.8068645）其他参数弹体长度L0.485m参考面积Sref0.0012566m2绕Z轴转动惯量0.02kg-m2全弹质量m01.06kg发动机推力P140N秒耗量ms0.06kg/s无量纲化后弹体俯仰阻尼系数mQz-1.5由舵面引起的升力系数C5zy0.006注：发动机工作时间为2s。2s后推力与秒耗量均为05仿真结果5.1状态10=5，x=1500m。采样时间间隔0.01s。0t0脱靶量：1.0317m。飞行时间7.95s。天邙平面弹交会20040060080010D01200144)01&00纵向位置(m)头Z平面弾；交会-51806040h-o0.10-20040060080010001200144)01&00纵向位置(m)才弹連度变化曲线弹门相对距离变化曲线E西.ENWE注5.2状态20二10,x二2000m,z二500m。采样时间间隔0.01s。0t0t0脱靶量：1.0072m。飞行时间12.2s。至】W起一上闽天邙平面弹交会10001500纵向悝置(m)X-Z平面弾；交至200-02&005004003M)2001000C/500W起一上龜10001500纵向位置(m)20