12031180
*****
飞控
1. 攻角 2°干扰 解题思路概括：根据已知量求出传递函数，在 simulink 里搭建模型并仿真。再将得到的 结果导入工作空间，在 matlab 里画出仿真曲线。 传递函数：
V ( s) 0.0057 s 2 0.0049875s 0.3110271 0 ( s) ( s 2 1.464472s 8.945125)( s 2 0.0121326s 0.001546684) ( s) s 3 0.891605s 2 0.0145294s 0.032175324 2 0 ( s) ( s 1.464472s 8.945125)( s 2 0.0121326s 0.001546684)
搭建 simulink 仿真结构图：
12031180
*****
飞控
短周期 10 秒
12031180
*****
飞控
长周期 200s 仿真曲线
12031180
*****
飞控
仿真曲线生成代码： subplot(311) plot(t2,V2,'m'); grid; xlabel('t/s'); ylabel('V/m/s') title('V 受舵偏 1°的变化曲线--10 秒'); subplot(312) plot(t2,alpha2,'b'); grid; xlabel('t/s'); ylabel('alpha/rad') title('alpha 受舵偏 1°的变化曲线--10 秒'); subplot(313) plot(t2,sita2,'r'); grid; xlabel('t/s'); ylabel('sita/rad') title('sita 受舵偏 1°的变化曲线--10 秒');
%%%%%%% subplot(411) plot(t2,gama2,'m');
12031180
*****
飞控
grid; xlabel('t/s'); ylabel('伽马/rad') title('伽马受舵偏 1°的变化曲线--10 秒'); subplot(412) plot(t2,q2,'b'); grid; xlabel('t/s'); ylabel('q/rad/s') title('俯仰角速度 q 受舵偏 1°的变化曲线--10 秒'); subplot(413) plot(t2,H2,'r'); grid; xlabel('t/s'); ylabel('H/m') title('高度 H 受舵偏 1°的变化曲线--10 秒'); subplot(414) plot(t2,delta2,'g'); grid; xlabel('t/s'); ylabel('delta/rad') title('舵偏变化曲线--10 秒');
( s) 8.42892s 0.134992 2 0 ( s) ( s 1.464472s 8.945125)( s 2 0.0121326s 0.001546684)
用 simulink 搭建模型，如下图：
12031180
*****
飞控
10s 短周期仿真曲线：
V ( s ) e ( s ) [( ( s ) e ( s ) ( s ) e ( s )
s s 1)( 1) 1.44 7251.4 s 2 0.0714 s 2 0.493 ) 2( ) s 1][( ) 2( ) s 1] 0.063 0.063 4.27 4.27 s s 2 0.0713 1.68( 1)[( ) 2( ) s 1] 248.5 0.068 0.068 s 2 0.0714 s 2 0.493 [( ) 2( ) s 1][( ) 2( ) s 1] 0.063 0.063 4.27 4.27 s s 4.85( 1)( 1) 0.0098 1.371 s 2 0.0714 s 2 0.493 [( ) 2( ) s 1][( ) 2( ) s 1] 0.063 0.063 4.27 4.27 24.142(
( s) ( s) ( s) e ( s) e ( s) e ( s) q( s) ( s) s e ( s) e ( s)
1 1 H ( s ) tan V ( s ) V0 ( s ) V ( s ) V0 s s
200s 仿真曲线：
仿真曲线生成代码： %攻角 2°干扰 200 仿真 subplot(311)控
plot(t,V,'m'); grid; xlabel('t/s'); ylabel('V/m/s') title('V 受攻角 2°干扰的变化曲线'); subplot(312) plot(t,alpha,'g'); grid; xlabel('t/s'); ylabel('alpha/°') title('alpha 受攻角 2°干扰的变化曲线'); subplot(313) plot(t,sita,'r'); grid; xlabel('t/s'); ylabel('sita/°') title('sita 受攻角 2°干扰的变化曲线'); 2. 1°舵偏 解题思路概括：根据已知量求出传递函数，在 simulink 里搭建模型并仿真。再将得到的 结果导入工作空间，在 matlab 里画出仿真曲线。 传递函数：