飞行控制系统大作业《飞行控制系统》课程实验报告 班 级 0314102 学 号 姓 名 孙旭东 成 绩 南京航空航天大学 4 月 (一)飞机纵向飞行控制系统的设计与仿真 1、分析飞机纵向动力学模态,求飞机的长周期与短周期阻尼与自然频率。在 MATLAB 环境下导入数据文件,输入 damp(alon),得出结果: Eigenvalue Damping Freq. (rad/s) -2.29e+000 + 4.10e+000i 4.88e-001 4.69e+000 -2.29e+000 - 4.10e+000i 4.88e-001 4.69e+000 -3.16e-002 1.00e+000 3.16e-002 -7.30e-003 + 3.35e-002i 2.13e-001 3.42e-002 -7.30e-003 - 3.35e-002i 2.13e-001 3.42e-002 长周期的根为 -7.30e-003 + 3.35e-002i 和 -7.30e-003 - 3.35e-002i阻尼为 2.13e-001自然频率为 3.42e-002(rad/s)短 周 期 的 根 为 -2.29e+000 + 4.10e+000i 和 -2.29e+000 - 4.10e+000i阻尼为 4.88e-001自然频率为 4.69e+000(rad/s)2、对升降舵及油门单位阶跃输入下的飞机自然特性进行仿真,画出相应的状态曲线。sys=ss(alon,blon,clon,dlon)[y,t]=step(sys,500)subplot(221)plot(t,y(:,1,1))xlabel('t(s)')ylabel('\Deltau(m/s)')subplot(222)plot(t,y(:,1,2))xlabel('t(s)')ylabel('\Deltau(m/s)')subplot(223)plot(t,y(:,2,1))xlabel('t(s)')ylabel('\Delta\alpha(deg)')subplot(224)plot(t,y(:,2,2))xlabel('t(s)')ylabel('\Delta\alpha(deg)')0200400600-2000200400t(s)u(m/s)0200400600-6-4-20t(s)(deg)0200400600-2000200400t(s)u(m/s)0200400600-2024t(s)(deg)subplot(221)plot(t,y(:,3,1))xlabel('t(s)')ylabel('\Deltaq(deg/s)')subplot(222)plot(t,y(:,3,2))xlabel('t(s)')ylabel('\Deltaq(deg/s)')subplot(223)0200400600-10-505t(s)q(deg/s)0200400600-4-2024t(s)q(deg/s)0200400600-150-100-50050t(s)(deg)0200400600-50050100t(s)(deg)plot(t,y(:,4,1))xlabel('t(s)')ylabel('\Delta\theta(deg)')subplot(224)plot(t,y(:,4,2))xlabel('t(s)')ylabel('\Delta\theta(deg)')subplot(121)plot(t,y(:,5,1))xlabel('t(s)')ylabel('\Deltah(m)')subplot(122)plot(t,y(:,5,2))xlabel('t(s)')ylabel('\Deltah(m)')0200400600-0.500.511.522.5x 104t(s)h(m)0200400600-2.5-2-1.5-1-0.50x 104t(s)h(m) 以上各图为升降舵及油门单位阶跃输入下的飞机自然特性行仿真,左边一列为升降舵的阶...