北航飞行力学理论与应用课程大作业第组

1、1飞行力学理论及应用大作业 院（系）自动化科学与电气工程学院 专 业控制科学与工程 学生姓名 朱付涛 BY1403158李琛 SY1403520左京兴 B1月4日飞行力学理论大作业大作业分工情况：李琛：飞机运动方程建模左京兴：小扰动线性化朱付涛：求解飞机时域响应和附录第一部分 飞机运动建模在平面地球假设下，飞机飞行在平面大地上空，不考虑地球曲率和地球自转。首先列出气动力和气动力矩方程：其中，为气动力在飞机机体系的分量，坐标系选为美式坐标系。为飞机气动力矩在机体上的分量。为动压，为飞机的机翼面积，为翼展，为平均气动弦长。1.1 飞机质心动力学方程由式(5.1.7)和式(

2、5.4.9)出发，推导飞机的质心动力学方程。本文推导了两种情况下的质心动力学方程：风轴系中的质心动力学方程和体轴系中的质心动力学方程。(1) 风轴系中的质心动力学方程式（）在风轴系中写为式（）是一般坐标系中惯性加速度的表示，将其写为地轴系中的表达，此时就是质心相对地球的速度，用来表示。假定地轴固定于惯性空间，则。原点的加速度就是与地球转动有关的向心加速度，它与g相比可以略去。其中的向心加速度项通常也可略去。另外令,哥氏加速度为，则有由于相对于的角速度为，利用将式转换到风轴系得到其中最后一项可以化简为：于是得到下式：其中利用矩阵表示式可得到和的表示为：将式和式代入到式，并假设大气静止，即W=0，

3、则可得： 将式代入到式中可得到标量方程为： 大气飞行的力矢量由两部分组成：气动反作用力A和重力mg，即 在风轴系中，A的分量可以写成下面的形式：上式中，D是阻力，C是侧力，L是升力。重力的表达如下式所示：因此利用风轴系的欧拉角，可得 将式、式和式。代入到式中可得如下方程组：当忽略地球转动时，和项为零，则有由此得到了风轴系下飞机的质心动力学方程。(2) 体轴系中的质心动力学方程式（）在体轴系中可写为设相对于的角速度为，则根据式及式可得：其中最后一项根据教材中式（）可化简为:则式可改为：式矢量的标量分量为根据式可得将式和式代入到式中，并考虑大气静止的情况，即W=0，得到如下方程组：类似的可得：考虑

4、到体轴系中的气动力表达式为：并利用式可得：即将式和式代入到式中得到体轴系下的质心动力学方程为：当忽略地球转动时，为零，则有：1.2飞机转动动力学方程由式（）出发，推导飞机的转动动力学方程（考虑发动机转子产生的角动量）。仅在体轴系下考虑：，式中 将式转换成矩阵符号为：根据式则可得h在中的分量为根据教材中式（），有，则上式可以化简为：利用可将上式中第二项改写为则式可重新写为：在以上分析的基础上，设两个体轴系和，相对于有一角速度，并考虑转子产生的角动量，则对这两个参考系的角动量为：令，则可将其转换成下列形式： 当为零，即，则得到了平均轴系下的关系式：下面再回到(5.4.12)式，它在参考系中的形式为

5、：那么转换到体轴系中为：其中为：当采用平均轴系时，将式代入到式中得到： 对于刚体，则上式中可以忽略，并进行标量展开为：将上式写成方程组的形式为：通常是飞行器的对称面，此时,则上式可以简化为列出整个方程的标量展开式没有什么好处，对于忽略的且没有转子项的局限情况，即对于刚体，式（1.46）标量方程为 1.3 飞机质心运动学方程由速度矢量三角形出发，推导飞机的质心运动学方程。飞机的质心运动学方程在风轴系和体轴系展开。(1) 风轴系的质心运动学方程由可写出风轴系中的方程,进而可以得到牵连垂直坐标系中的分量为：式中 将式展开成标量形式为： 当忽略地球曲率时（），即可认为平行于，并且假设大气相对于地球静止

6、，即W=0，而质心的位置坐标可写成下面的形式：上式便是风轴系中飞机质心运动学方程。(2) 体轴系中运动学方程由可写出体轴系中的方程,进而得到牵连垂直坐标系中的分量：式中当忽略地球曲率时（），即可认为平行于，并且假设大气相对于地球静止，即W=0，而质心的位置坐标可写成下面的形式：上式便是体轴系中的飞机质心运动学方程。1.4飞机转动运动学方程由(5.2.5)式出发，推导飞机的转动运动学方程。(1) 体轴系转动运动学方程由式(5.2.5)可得所以将上式变换得到欧拉角的变化率的表示为：上式就是体轴系中的飞机转动运动学方程。(2) 风轴系转动运动学方程与体轴系中的推导过程完全相同，最终得到的转动运动学方

7、程为：1.5六自由度全量方程组当无风时( W = 0 )，对于具有对称面的刚体飞机，建立飞机的6自由度全量运动方程组飞机的质心动力学方程：对于具有对称面的刚体飞机的转动动力学方程为：飞机的质心运动学方程为：飞机的转动运动学方程为：第二部分 小扰动线性化设基准运动为对称定常直线水平飞行，对飞机全量运动方程组进行小扰动线性化处理。2.1飞机全量运动方程组飞机的全量运动方程组如下所示，在此基础上进行小扰动线性化。由所给的F-16气动参数表，可以得到气动力和气动力矩的表达式为：气动力系数需要按如下方法通过插值求出：以上为飞机的全量运动方程组。2.2小扰动线性化方程组设基准运动为对称定常直线水平飞行，对

8、刚体飞机全量运动方程组进行小扰动线性化处理。对如下纵向非线性方程进行小扰动线性化，有其中上述四个方程的基准运动方程为：用式减去式并化简得到：然后将所给的力和力矩的导数形式代入到式中得到：然后以为状态变量，以为控制变量，可将式线化成的形式，即以上是纵向方程组的线化形式，同理可得到横侧向方程组（其中状态变量，控制变量）的线化形式为：2.3求气动导数在第二部分中已经建立了分立的小扰动线性化方程，式子中涉及到了无因次量的导数，需要求出这些导数的计算式，计算过程如下，此过程中。要建立小扰动线性化方程，在纵向，需要求出气动力和气动力矩关于的导数，在横侧向，需要求出气动力和气动力矩关于的导数。(1) 轴向力

9、的导数导数是在基准状态下求得的，所以有轴向力的导数有、：其中，另外有：(2) 法向力的导数轴向力的导数有、：(3) 力矩的导数力矩的导数有、： (4) 侧向力的导数侧向力的导数有、： (5) 力矩的导数力矩的导数有、：(6) 力矩的导数力矩的导数有、：第三部分 飞机时域响应求解已知某飞机数据，建立描述飞机运动特性的全量运动方程。飞机的惯性数据如下表所示：表1 飞机的惯性数据重量(lbs)转动惯量(slug-ft2)WIxIyIzIxz2500094965581463100982重量(kg)转动惯量(kg-m2)WIxIyIzIxz11339.8093飞机机翼几何尺寸如下表所示：表2 飞机的几何

10、尺寸翼展b(ft)机翼面积S(ft2)平均气动弦(ft)3030011.32翼展b(m)机翼面积S(m2)平均气动弦(m)由于不同的大气密度计算公式算出的大气密度差别很大，所以采用常用的大气模型，1976年美国标准大气模型USSA76，该模型以地理纬度45°3233地区海平面为基准，全年实际大气参数的统计平均值为标准大气参数。大气基准值为：温度，密度，声速。USSA适用高度范围为，并分成低层和高层部分。低层高度为，在这里只介绍低层大气。低层大气空气分子量为常数，大气是完全气体，处在静平衡状态，其位势高度为，重力加速度为，为海拔，与对应的地球半径为。于是，。 大气参数分段计算公式：(1

11、). (2). 的大气参数在此不列出。声速。3.1飞机全量运动方程组其中，气动力和气动力矩为：气动力系数需要按如下方法通过插值求出：下面以定直平飞状态作为基准运动状态，完成纵向力和力矩的配平，得到配平迎角、配平升降舵偏角和配平油门的大小。3.2纵向力和力矩的配平地面的重力加速度g=9.8066 m/s2、大气密度。，定直平飞状态下，。纵向方向配平需要用到的力和力矩方程为：纵向气动力和气动力矩的表达式为又有、和，可将上式进一步化简为由气动参数表查表可知，Cz=A*(1-(beta/57.3)2)-0.19*(elevator/25.0)，在定直平飞状态中，beta为零，于是配平方程进一步写为：在

12、给定飞行状态时，上式为关于的非线性方程组。用matlab的fsolve函数可以先解出，再解出。每一种飞行状态对应的速度和高度如下表所示：表3 飞机的飞行状态飞行状态速度(m/s)高度(m)飞行状态速度(m/s)高度(m)1500131705000260014180100003700151901000048001620010000590017210100006100018220110007110100019230110008120100020240110009130100021250120001014010002226013000111505000232701400012160500024280

13、15000配平攻角、升降舵偏角及油门的值为：图1 对应所有飞行状态的配平结果我们这一组（第四组）所取飞行状态的速度为120m/s，高度为1000m，其配平攻角、舵偏角和油门开度为：6.3296°、-0.5714°和0.169。3.3小扰动线性化方程组(1) 纵向小扰动方程在第二部分中已经建立了分立的小扰动线性化方程和气动导数的求解公式，根据3.2节的配平结果，可以求得每一个气动导数，进一步求得纵向和横侧向的小扰动线性化方程。其中，纵向的小扰动线性化方程为：其中，纵向系统矩阵的特征值和对应的特征向量为：和(2) 横侧向小扰动方程横侧向的小扰动线性化方程为： 横侧向系统矩阵的特

14、征值和对应的特征向量为：和至此，求出了小扰动线性化方程的具体形式。3.4计算时域响应基于全量运动方程和小扰动线化方程，计算飞机各状态变量的时域响应。飞机的状态变量为：体轴系速度；角速度；姿态角；位移。输入为舵偏角；发动机油门开度。配平攻角、升降舵和油门为6.3296°、-0.5714和0.169。 (1). 状态变量初值为零，外界输入为零对于全量方程：初值为零：=120*cos(6.3296/57.3) 0 120*sin(6.3296/57.3) 0 0 0 0 6.3296/57.3 0 0 0 0。外界输入为零：=0 -0.5714°0 0.169。对于小扰动方程：初

15、值为零：和。 外界输入为零：。(2). 纵向模态激励下，外界输入为零纵向特征值为和。对应的特征向量的实部分别为和。对于全量方程：对于长周期模态，模态激励为=120*cos(6.3296/57.3) 0 120*sin(6.3296/57.3) 0 0 0 0 6.3296/57.3 0 0 0 0+，外界输入为=0 -0.5714°0 0.169；对于短周期模态，模态激励为=120*cos(6.3296/57.3) 0 120*sin(6.3296/57.3) 0 0 0 0 6.3296/57.3 0 0 0 0+ ，外界输入为=0 -0.5714° 0 0.169。对于

16、小扰动方程：对于长周期模态，模态激励为，外界输入为；对于短周期模态，模态激励为，外界输入为。(3). 侧向模态激励下，外界输入为零侧向特征值为、和，现在取一个稳定模态和一个不稳定模态。对应的特征向量分别为和。对于全量方程：对于长周期模态，模态激励为=120*cos(6.3296/57.3) 0 120*sin(6.3296/57.3) 0 0 0 0 6.3296/57.3 0 0 0 0+ ，外界输入为=0 -0.5714°0 0.169；对于短周期模态，模态激励为=120*cos(6.3296/57.3) 0 120*sin(6.3296/57.3) 0 0 0 0 6.3296

17、/57.3 0 0 0 0+ ，外界输入为=0 -0.5714° 0 0.169。对于小扰动方程：对于长周期模态，模态激励为，外界输入为；对于短周期模态，模态激励为，外界输入为。(4). 升降舵单位阶跃输入（1º0.01745 rad）对于全量方程：初值为零：=120*cos(6.3296/57.3) 0 120*sin(6.3296/57.3) 0 0 0 0 6.3296/57.3 0 0 0 0。外界输入为=0 -0.5714° 0 0.169+ 0 step 0 0。对于小扰动方程：初值为零：和。 外界输入为。(5). 副翼单位阶跃输入（1º0.

18、01745 rad）对于全量方程：初值为零：=120*cos(6.3296/57.3) 0 120*sin(6.3296/57.3) 0 0 0 0 6.3296/57.3 0 0 0 0。外界输入为=0 -0.5714° 0 0.169+ step 0 0 0。对于小扰动方程：初值为零：和。 外界输入为。初值为零输入为零(1) 全量方程模型全量方程模型在MATLAB中的实现如下图所示：图2 零状态零输入的全量方程模型全量方程在初值为零无输入情况下的响应结果如下图所示：图3 初值为零输入为零全量方程的三轴速度图4 初值为零输入为零全量方程的三轴角速度图5 初值为零输入为零全量方程的三

19、轴姿态角(2) 小扰动线性化方程模型小扰动线性化模型在MATLAB中的实现如下图所示：图6 零状态零输入的小扰动线性化模型小扰动线性化模型在初值为零无输入的情况下的响应如下图所示：图7 初值为零输入为零小扰动方程的三轴姿速度图8 初值为零输入为零小扰动方程的三轴姿角速度图9 初值为零输入为零小扰动方程的姿态角结果的分析：从图3-5和图7-9可以看出，零状态零输入情况下全量方程和小扰动方程的响应结果是一致的。小扰动方程的响应一直为零，这与小扰动方程的特性是一致的。全量方程的响应与小扰动方程的响应结果一致，在开始的地方有振荡调整过程，其纵向状态的最终值有微小的振荡。纵向模态激励下全量方程在纵向模态

20、激励下的响应，分为2种情况：以对应最大模特征值的特征向量的实部为状态初始值；以对应最小模特征值的特征向量的实部为状态初始值。 (1) 全量方程模型全量方程模型的MATLAB实现如下图所示：图10纵向模态激励全量方程模型以对应最大模特征值的特征向量的实部为状态初始值，纵向模态激励全量模型的响应如下图所示：图11 初值为零大特征值纵向模态激励全量方程的三轴速度图12 初值为零大特征值纵向模态激励全量方程的三轴角速度图13 初值为零大特征值纵向模态激励全量方程的三轴姿态角以对应最小模特征值的特征向量的实部为状态初始值：图14 初值为零小特征值纵向模态激励全量方程的三轴速度图15 初值为零小特征值纵向

21、模态激励全量方程的三轴角速度图16 初值为零小特征值纵向模态激励全量方程的三轴姿态角(2) 小扰动线性化方程模型小扰动线性化模型的MATLAB实现如下图所示：图17纵向模态激励小扰动方程模型小扰动方程的响应，大特征值对应的模态激励下的响应： 图18 初值为零大特征值纵向模态激励小扰动方程的三轴速度图19 初值为零大特征值纵向模态激励小扰动方程的三轴角速度图20 初值为零大特征值纵向模态激励小扰动方程的姿态角小扰动方程的响应，小特征值对应的模态激励下的响应：图21 初值为零小特征值纵向模态激励小扰动方程的三轴速度图22 初值为零小特征值纵向模态激励小扰动方程的三轴角速度图23 初值为零小特征值纵

22、向模态激励小扰动方程的姿态角分析：对比图11-13和图14-16可知，全量方程模型在两种纵向激励作用下的响应差别不大。对比图18-20和图21-23可知，小扰动方程模型在长周期模态激励和短周期激励作用下的响应差别明显，短周期响应在4秒时趋于稳定，而长周期响应需要近1000秒时间才能趋于稳定。从图11-23可以看出，小扰动线性化方程模型便于分析飞机的运动特性，能够准确分析飞机的长周期模态和短周期模态。而全量方程模型虽然能够体现出一定的短周期特性，但整体上分析不出飞机的运动特性。横侧向模态激励下飞机横侧向特征值为、和，现在取一个稳定模态和一个不稳定模态来求取响应。对应的特征向量分别为和。(1) 全

23、量方程模型全量方程的MATLAB实现如下图所示：图24横侧向模态激励全量方程模型以大特征值对应的特征向量为初始状态求得响应如下图所示：图25 初值为零大特征值横侧向模态激励全量方程的速度图26 初值为零大特征值横侧向模态激励全量方程的角速度图27 初值为零大特征值横侧向模态激励全量方程的姿态角以小特征值对应的特征向量为初始状态求得响应如下图所示：图28 初值为零小特征值横侧向模态激励全量方程的速度图29 初值为零小特征值横侧向模态激励全量方程的角速度图30 初值为零小特征值横侧向模态激励全量方程的姿态角(2) 小扰动线性化方程模型小扰动线性化模型在MATLAB中的实现如下图所示：图31横侧向模

24、态激励小扰动方程模型小扰动线性化模型在横侧向模态激励下的响应，大实根的特征向量为状态初始值：图32 初值为零大特征值横侧向模态激励小扰动方程的速度图33 初值为零大特征值横侧向模态激励小扰动方程的角速度图34 初值为零大特征值横侧向模态激励小扰动方程的姿态角小实根的特征向量为状态初始值：图35 初值为零小特征值横侧向模态激励小扰动方程的速度图36 初值为零小特征值横侧向模态激励小扰动方程的角速度图37 初值为零小特征值横侧向模态激励小扰动方程的姿态角分析：横侧向长周期模态对应的特征值为稳定特征点，而短周期模态对应的特征值是不稳定的，这从图25-30和图32-37也能看出，无论是全量方程模型还是

25、小扰动线性化模型，短周期模态激励下的状态响应都是不稳定的。对于长周期模态激励下的响应，小扰动线性化模型的结果能够很快的收敛；而全量方程模型的结果收敛速度较慢。 升降舵单位阶跃输入全量方程模型在MATLAB中的实现如下图所示。其中，配平升降舵偏角为-0.5714，配平油门开度为0.169，升降舵阶跃输入的幅值为0.05。(1) 全量方程模型图38 零状态升降舵阶跃输入全量方程模型全量模型在升降舵阶跃输入下的响应结果如下图所示：图39 零状态升降舵阶跃输入全量方程的速度图40 零状态升降舵阶跃输入全量方程的角速度图41 零状态升降舵阶跃输入全量方程的姿态角(2) 小扰动线性化方程模型小扰动线性化方

26、程模型在MATLAB中的实现如下图所示：图42 零状态升降舵阶跃输入小扰动方程模型小扰动线性化模型在升降舵单位阶跃输入下的响应结果如下图所示：图43 零状态升降舵阶跃输入小扰动方程的速度图44 零状态升降舵阶跃输入小扰动方程的角速度图45 零状态升降舵阶跃输入小扰动方程的姿态角分析：从图39-45可以看出，飞机状态量最终趋于稳定，说明飞机纵向动态稳定；对比图39-41和图43-45可以看出全量方程模型和小扰动模型的动态调整时间都较长。飞机的纵向动态稳定和2对复根的实部都小于0相对应。副翼单位阶跃输入副翼阶跃输入下全量方程在MATLAB中的实现如下图所示，副翼阶跃输入幅值为0.01。(1) 全量

27、方程模型图46 零状态副翼阶跃输入全量方程模型全量模型在副翼阶跃输入下的响应如下图所示：图47 零状态副翼阶跃输入全量方程的速度图48 零状态副翼阶跃输入全量方程的角速度图49 零状态副翼阶跃输入全量方程的姿态角(2) 小扰动线性化方程模型小扰动线性化方程在MATLAB中的实现如下图所示：图50 零状态副翼阶跃输入小扰动方程模型小扰动线性化模型在副翼单位阶跃输入下的响应如下图所示：图51 零状态副翼阶跃输入小扰动方程的速度图52 零状态副翼阶跃输入全量方程的角速度图53 零状态副翼阶跃输入全量方程的姿态角分析：从图47-49和图51-53可以看出，全量方程模型和小扰动线性化模型的副翼阶跃响应都

28、是不稳定的，这与飞机横侧向有正实根是相对应的。附录5.1 配平程序5.1.1 配平主函数clear all;close all;clc%-所给数据及矩阵将需要插值求出的气动数据设为全局变量-%global given_Alpha given_Elevator given_Cz given_Cx given_Cm given_Cl given_Cn given_Cl_delta_a given_Cl_delta_rglobal given_Cn_delta_a given_Cn_delta_r given_throtte_000 given_throtte_077 given_throtte_1

29、00global given_Cxq given_Cyr given_Cyp given_Czq given_Clr given_Clp given_Cmq given_Cnr given_Cnp given_delta_beta given_beta given_height given_mahglobal m g c_ S b r2d H rou Ix Iy Iz Ixz V0 Mah0global Result_alpha Result_elevator Result_Throtte T resultalfaelevator R g0%-地球半径，地面的重力加速度-%R = 635676

30、6;g0 = 9.80665;%-飞机的常数（国际单位制）-%m = 9307.8093;c_= 3.4503;S = 27.8709;b = 9.144;r2d=57.29577951;%-飞机的惯性数据（国际单位制）-%Ix = 12874.8446; Iy = 75673.6077; Iz = 85552.0953; Ixz = 1331.4130; Isum=Ix*Iz-Ixz2;%-F-16的气动参数数据-%given_Alpha = -10 -5 0 5 10 15 20 25 30 35 40 45;given_Elevator = -24 -12 012 24;given_Cz

31、 =0.7700 0.2410 -0.1000 -0.4160 -0.7310 -1.0530 -1.3660 -1.6460 -1.9170 -2.1200 -2.2480 -2.2290;given_Cx = -0.099-0.048-0.022-0.04-0.083 -0.081-0.038-0.02-0.038-0.073 -0.081-0.04-0.021-0.039-0.076 -0.063-0.021-0.004-0.025-0.072 -0.025 0.016 0.032 0.006-0.046 0.044 0.083 0.094 0.062 0.012 0.097 0.127

32、 0.128 0.087 0.024 0.113 0.137 0.13 0.085 0.025 0.145 0.162 0.154 0.1 0.043 0.167 0.177 0.161 0.11 0.053 0.174 0.179 0.155 0.104 0.047 0.166 0.167 0.138 0.091 0.04;given_Cm = 0.205 0.081-0.046-0.174-0.259 0.168 0.077-0.02-0.145-0.202 0.186 0.107-0.009-0.121-0.184 0.196 0.11-0.005-0.127-0.193 0.213 0

33、.11-0.006-0.129-0.199 0.251 0.141 0.01-0.102-0.15 0.245 0.127 0.006-0.097-0.16 0.238 0.119-0.001-0.113-0.167 0.252 0.133 0.014-0.087-0.104 0.231 0.108 0 -0.084-0.076 0.198 0.081-0.013-0.069-0.041 0.192 0.093 0.032-0.006-0.005;given_Cl = 0-0.001-0.003-0.001 0 0.07 0.009 0-0.004-0.009-0.01-0.01-0.01-0

34、.011 0-0.008-0.017-0.02-0.022-0.023-0.023 0-0.012-0.024-0.03-0.034-0.034-0.037 0-0.016-0.03-0.039-0.047-0.049-0.05 0-0.019-0.034-0.044-0.046-0.046-0.047 0-0.02-0.04-0.05-0.059-0.068-0.074 0-0.02-0.037-0.049-0.061-0.071-0.079 0-0.015-0.016-0.023-0.033-0.06-0.091 0-0.008-0.002-0.006-0.036-0.058-0.076

35、0-0.013-0.1-0.014-0.035-0.062-0.077 0-0.015-0.19-0.027-0.035-0.059-0.076;given_Cn = 00.0180.0380.0560.0640.0740.079 00.0190.0420.0570.0770.0860.09 00.0180.0420.0590.0760.0930.106 00.0190.0420.0580.0740.0890.106 00.0190.0430.0580.0730.080.096 00.0180.0390.0530.0570.0620.08 00.0130.030.0320.0290.0490.

36、068 00.0070.0170.0120.0070.0220.03 00.0040.0040.0020.0120.0280.064 0 -0.014 -0.035 -0.046 -0.034 -0.0120.015 0 -0.017 -0.047 -0.071 -0.065 -0.0020.011 0 -0.033 -0.057 -0.073 -0.041 -0.013 -0.001;given_Cl_delta_a = -0.041-0.041-0.042-0.04-0.043-0.044-0.043 -0.052-0.053-0.053-0.052-0.049-0.048-0.049 -

37、0.053-0.053-0.052-0.051-0.048-0.048-0.047 -0.056-0.053-0.051-0.052-0.049-0.047-0.045 -0.05-0.05-0.049-0.048-0.043-0.042-0.042 -0.056-0.051-0.049-0.048-0.042-0.041-0.037 -0.082-0.066-0.043-0.042-0.042-0.02-0.003 -0.059-0.043-0.035-0.037-0.036-0.028-0.013 -0.042-0.038-0.026-0.031-0.025-0.013-0.01 -0.0

38、38-0.027-0.016-0.026-0.021-0.014-0.003 -0.027-0.023-0.018-0.017-0.016-0.011-0.007 -0.017-0.016-0.014-0.012-0.011-0.01-0.008;given_Cl_delta_r = 0.0050.0070.0130.0180.0150.0210.023 0.0170.0160.0130.0150.0140.0110.01 0.0140.0140.0110.0150.0130.010.011 0.01 0.0140.0120.0140.0130.0110.011 -0.0050.0130.01

39、10.0140.0120.010.011 0.0090.0090.0090.0140.0110.0090.01 0.0190.0120.0080.0140.0110.0080.008 0.0050.0050.0050.0150.010.010.01 0 0 -0.0020.0130.0080.0060.006 -0.0050.0040.0050.0110.0080.0050.014 -0.0110.0090.0030.0060.0070 0.02 0.0080.0070.0050.0010.0030.0010 ;given_Cn_delta_a = 0.001 0.002-0.006-0.01

40、1-0.015-0.024-0.022 -0.027-0.014-0.008-0.011-0.015-0.01 0.002 -0.017-0.016-0.006-0.01-0.014-0.004-0.003 -0.013-0.016-0.006-0.009-0.012-0.002-0.005 -0.012-0.014-0.005-0.008-0.011-0.001-0.003 -0.016-0.019-0.008-0.006-0.008 0.003-0.001 0.001-0.021-0.005 0 -0.002 0.014-0.009 0.017 0.002 0.007 0.004 0.00

41、2 0.006-0.009 0.011 0.012 0.004 0.007 0.006-0.001-0.001 0.017 0.015 0.007 0.1 0.012 0.004 0.003 0.008 0.015 0.006 0.004 0.011 0.004-0.002 0.016 0.011 0.006 0.1 0.011 0.006 0.001;given_Cn_delta_r = -0.018-0.028-0.037-0.048-0.043-0.052-0.062 -0.052-0.051-0.041-0.045-0.044-0.034-0.034 -0.052-0.043-0.03

42、8-0.045-0.041-0.036-0.027 -0.052-0.046-0.04-0.045-0.041-0.036-0.028 -0.054-0.045-0.04-0.044-0.04-0.035-0.027 -0.049-0.049-0.038-0.045-0.038-0.028-0.027 -0.059-0.057-0.037-0.047-0.034-0.024-0.023 -0.051-0.052-0.03-0.048-0.035-0.023-0.023 -0.03-0.03-0.027-0.049-0.035-0.02-0.019 -0.037-0.033-0.024-0.04

43、5-0.029-0.016-0.009 -0.026-0.03-0.019-0.033-0.022-0.01-0.025 -0.013-0.008-0.013-0.016-0.009-0.014-0.01;given_throtte_000=1060 670 880 114015001860 635 425 690 101013301700 60 25 345 75511301525 -1020 -710-300 350910 1360 -2700 -1900-1300-247600 1100 -3600 -1400-595-342 -200 700;given_throtte_077=126

44、8091506200395024501400 1268091506313404024701400 1261093126610429026001560 1264098397090466028401660 12390101767750532032501930 1168098488050610038002310;given_throtte_100=200001500010800700040002500 214201570011225732344352600 227001686012250815450002835 242401891013760928557003215 2607021075159751

45、111568603950 2888623319183001348486425057;given_Cxq = -0.267 -0.11 0.308 1.34 2.08 2.91 2.76 2.05 1.5 1.49 1.83 1.21;given_Cyr = 0.882 0.852 0.876 0.958 0.962 0.974 0.819 0.483 0.59 1.21 -0.493 -1.04;given_Cyp = -0.108 -0.108 -0.188 0.11 0.258 0.226 -0.344 0.362 0.611 0.529 0.298 -2.27;given_Czq = -

46、8.8 -25.8 -28.9 -31.4 -31.2 -30.7 -27.7 -28.2 -29 -29.8 -38.3 -35.3;given_Clr = -0.126 -0.126 0.063 0.113 0.208 0.23 0.319 0.437 0.68 0.1 0.447 -0.33;given_Clp = -0.36 -0.359 -0.443 -0.42 -0.383 -0.375 -0.329 -0.294 -0.23 -0.21 -0.12 -0.1;given_Cmq = -7.21 -0.54 -5.23 -5.26 -6.11 -6.64 -5.69 -6 -6.2 -6.4 -6.6 -6;given_Cnr = -0.38 -0.363 -0.378 -0.386 -0.37 -0.453 -0.55 -0.582 -0.595 -0.637 -1.02 -0.804;given_Cnp = 0.061 0.052 0.052 -0.012 -0.013 -0.024 0.05 0.15 0.13 0.158 0