%/////////////////////////////////////////////////////////////////////
%升力系数
%CLBasic
Alphadgtable=[-5:1:25];
Flapstable=[0,1,5,10,20,25,30];
CLBasicgg=[-999,0,1,5,10,20,25,30;
-5,-0.48,-0.44,-0.56, -0.40,-0.26,-0.04,0.20 ;
-4,-0.40,-0.36,-0.46,-0.30,-0.13,0.10,0.36;
-3,-0.30,-0.25,-0.34,-0.19,-0.01,0.22,0.50;
-2,-0.21,-0.16,-0.24,-0.08,0.12,0.35,0.62;
-1,-0.12,-0.06,-0.14,0.03,0.24,0.47,0.76;
0,-0.04,0.02,-0.04,0.13,0.35,0.58,0.87;
1,0.06,0.12,0.08,0.23,0.47,0.70,0.99;
2,0.14,0.21,0.18,0.34,0.59,0.80,1.09;
3,0.22,0.31,0.28,0.44,0.70,0.90,1.19;
4,0.32,0.40,0.39,0.55,0.81,1.01,1.30;
5,0.41,0.50,0.50,0.65,0.92,1.11,1.40;
6,0.49,0.60,0.60,0.75,1.04,1.22,1.50;
7,0.59,0.68,0.70,0.85,1.14,1.32,1.61;
8,0.68,0.77,0.81,0.95,1.24,1.41,1.72;
9,0.77,0.86,0.91,1.05,1.35,1.52,1.82;
10,0.86,0.96,1.01,1.16,1.45,1.62,1.93;
11,0.94,1.04,1.12,1.26,1.55,1.72,2.02;
12,1.02,1.12,1.22,1.36,1.65,1.83,2.13;
13,1.07,1.20,1.33,1.47,1.74,1.93,2.21;
14,1.11,1.26,1.44,1.57,1.83,2.02,2.28;
15,1.14,1.32,1.54,1.66,1.90,2.10,2.33;
16,1.18,1.37,1.63,1.75,1.96,2.16,2.38;
17,1.20,1.41,1.70,1.81,2.01,2.21,2.42;
18,1.22,1.44,1.74,1.84,2.06,2.25,2.45;
19,1.24,1.46,1.78,1.86,2.08,2.26,2.46;
20,1.26,1.48,1.80,1.88,2.08,2.25,2.45;
21,1.28,1.50,1.80,1.89,2.07,2.22,2.41;
22,1.30,1.52,1.81,1.90,2.06,2.18,2.35;
23,1.32,1.54,1.82,1.90,2.04,2.15,2.29;
24,1.34,1.55,1.82,1.91,2.03,2.12,2.23;
25,1.36,1.57,1.83,1.92,2.02,2.09,2.19;];
CLBasictable=CLBasicgg(2:32,2:8)';
%CLAlphadot
Mtable=[0:0.025:0.995];
CLAlphadotgg=[
0, -6.70;0.025, -6.70;0.050, -6.70;0.075, -6.70;0.100, -6.70;
0.125, -6.70;0.150, -6.70;0.175, -6.70;0.200, -6.70;0.225, -6.70;
0.250, -6.70;0.275, -6.70;0.300, -6.70;0.325, -6.70;0.350, -6.67;
0.375, -6.63;0.400, -6.62;0.425, -6.61;0.450, -6.60;0.475, -6.60;
0.500, -6.60;0.525, -6.58;0.550, -6.55;0.575, -6.53;0.600, -6.50;
0.625, -6.48;0.650, -6.46;0.675, -6.43;0.700, -6.38;0.725, -6.30;
0.750, -6.23;0.775, -6.15;0.800, -6.04;0.825, -5.92;0.850, -5.80;
0.875, -5.63;0.900, -5.48;0.925, -5.67;0.950, -6.03;0.970, -6.17];
CLAlphadottable=CLAlphadotgg(:,2)';
%CLqdot
Hinchtable=[0,20000,40000];
CLqdotgg=[-999, 0, 20000, 40000;
0, 6.78, 6.78, 6.78;
0.025, 6.60, 6.60, 6.73;
0.050, 6.45, 6.53, 6.68;
0.075, 6.30, 6.40, 6.62;
0.100, 6.20, 6.30, 6.60;
0.125, 6.00, 6.20, 6.55;
0.150, 5.90, 6.12, 6.50;
0.175, 5.75, 6.00, 6.48;
0.200, 5.65, 5.93, 6.40;
0.225, 5.50, 5.85, 6.38;
0.250, 5.40, 5.75, 6.33;
0.275, 5.30, 5.70, 6.30;
0.300, 5.20, 5.60, 6.24;
0.325, 5.10, 5.50, 6.20;
0.350, 5.00, 5.43, 6.15;
0.375, 4.90, 5.37, 6.10;
0.400, 4.80, 5.30, 6.07;
0.425, 4.73, 5.28, 6.04;
0.450, 4.65, 5.25, 6.00;
0.475, 4.57, 5.20, 5.98;
0.500, 4.50, 5.18, 5.95;
0.525, 4.40, 5.17, 5.93;
0.550, 4.30, 5.15, 5.90;
0.575, 4.20, 5.13, 5.88;
0.600, 4.15, 5.12, 5.85;
0.625, 4.05, 5.10, 5.85;
0.650, 3.95, 5.10, 5.83;
0.675, 3.85, 5.08, 5.82;
0.700, 3.75, 5.06, 5.83;
0.725, 3.65, 5.05, 5.85;
0.750, 3.55, 5.02, 5.87;
0.775, 3.45, 5.00, 5.90;
0.800, 3.40, 4.98, 5.97;
0.825, 3.33, 5.00, 6.10;
0.850, 3.40, 5.10, 6.30;
0.875, 3.50, 5.30, 6.65;
0.900, 3.60, 5.50, 6.90;
0.925, 3.15, 5.18, 6.65;
0.950, 2.75, 4.70, 6.35;
0.970, 2.30, 4.35, 6.10;];
CLqtable=CLqdotgg(2:41,2:4)';
%CLdEi
Kalphatablegg=[-999,0,1,5,10,20,25,30;
-5,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
-4,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
-3,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
-2,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
-1,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
0,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
1,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
2,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
3,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
4,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
5,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
6,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
7,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
8,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
9,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
10,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
11,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
12,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
13,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
14,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
15,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000,1.0000;
16,1.0000,0.9968,1.0000,0.9900,0.9913,0.9882,0.9924;
17,1.0000,0.9903,1.0000,0.9750,0.9704,0.9592,0.9714;
18,0.9980,0.9772,1.0002,0.9500,0.9402,0.8750,0.9359;
19,0.9900,0.9576,0.9963,0.9100,0.8996,0.7816,0.8860;
20,0.9767,0.9315,0.9897,0.8700,0.8604,0.7211,0.8296;
21,0.9647,0.9027,0.9819,0.8300,0.8263,0.7039,0.7705;
22,0.9600,0.8675,0.9635,0.8000,0.8027,0.6961,0.7206;
23,0.9540,0.8296,0.9385,0.7750,0.7975,0.6961,0.7088;
24,0.9450,0.7969,0.9070,0.7500,0.7962,0.6961,0.7338;
25,0.9250,0.7799,0.8703,0.7350,0.7949,0.6961,0.7692];
Kalphatable=Kalphatablegg(2:32,2:8)';
Hinchtable2=[0, 10000, 20000, 30000, 40000];
CLdEigg=[-999, 0, 10000, 20000, 30000, 40000;
0, 0.00325, 0.00325, 0.00325,0.00325,0.00325;
0.025, 0.00322, 0.00324, 0.00324,0.00326,0.00327;
0.050, 0.00320, 0.00323, 0.00325,0.00327,0.00328;
0.075, 0.00318, 0.00322, 0.00325,0.00328,0.00329;
0.100, 0.00315, 0.00321, 0.00325,0.00329,0.00330;
0.125, 0.00313, 0.00320, 0.00325,0.00330,0.00331;
0.150, 0.00310, 0.00317, 0.00325,0.00331,0.00332;
0.175, 0.00305, 0.00315, 0.00325,0.00332,0.00334;
0.200, 0.00300, 0.00312, 0.00325,0.00333,0.00336;
0.225, 0.00295, 0.00310, 0.00325,0.00334,0.00338;
0.250, 0.00290, 0.00307, 0.00324,0.00335,0.00340;
0.275, 0.00285, 0.00304, 0.00323,0.00335,0.00342;
0.300, 0.00280, 0.00302, 0.00322,0.00336,0.00343;
0.325, 0.00273, 0.00298, 0.00320,0.00336,0.00345;
0.350, 0.00268, 0.00295, 0.00319,0.00337,0.00346;
0.375, 0.00260, 0.00292, 0.00318,0.00337,0.00348;
0.400, 0.00255, 0.00290, 0.00317,0.00337,0.00350;
0.425, 0.00245, 0.00285, 0.00316,0.00337,0.00351;
0.450, 0.00240, 0.00282, 0.00315,0.00337,0.00352;
0.475, 0.00235, 0.00278, 0.00313,0.00337,0.00354;
0.500, 0.00225, 0.00275, 0.00311,0.00337,0.00355;
0.525, 0.00220, 0.00272, 0.00310,0.00336,0.00355;
0.550, 0.00213, 0.00268, 0.00308,0.00335,0.00356;
0.575, 0.00205, 0.00264, 0.00306,0.00334,0.00356;
0.600, 0.00200, 0.00260, 0.00304,0.00334,0.00356;
0.625, 0.00193, 0.00255, 0.00302,0.00332,0.00355;
0.650, 0.00187, 0.00251, 0.00300,0.00331,0.00354;
0.675, 0.00180, 0.00248, 0.00298,0.00331,0.00354;
0.700, 0.00175, 0.00245, 0.00295,0.00330,0.00353;
0.725, 0.00170, 0.00240, 0.00294,0.00329,0.00351;
0.750, 0.00170, 0.00238, 0.00290,0.00326,0.00350;
0.775, 0.00170, 0.00237, 0.00289,0.00325,0.00348;
0.800, 0.00167, 0.00236, 0.00286,0.00324,0.00345;
0.825, 0.00165, 0.00235, 0.00284,0.00319,0.00342;
0.850, 0.00163, 0.00232, 0.00278,0.00315,0.00336;
0.875, 0.00162, 0.00223, 0.00268,0.00307,0.00330;
0.900, 0.00160, 0.00210, 0.00253,0.00292,0.00313;
0.925, 0.00158, 0.00193, 0.00225,0.00257,0.00284;
0.950, 0.00157, 0.00167, 0.00190,0.00226,0.00254;
0.970, 0.00156, 0.00150, 0.00175,0.00206,0.00235];
CLdEitable=CLdEigg(2:41,2:6)';
%CLdEotable
CLdEogg=[-999, 0, 10000, 20000, 30000, 40000;
0, 0.00380,0.00380,0.00380,0.00380,0.00380;
0.025,0.00372,0.00375,0.00375,0.00378,0.00380;
0.050,0.00364,0.00370,0.00375,0.00378,0.00380;
0.075,0.00358,0.00365,0.00370,0.00375,0.00380;
0.100,0.00350,0.00360,0.00366,0.00375,0.00380;
0.125,0.00344,0.00352,0.00365,0.00375,0.00380;
0.150,0.00335,0.00348,0.00360,0.00375,0.00386;
0.175,0.00326,0.00342,0.00360,0.00375,0.00385;
0.200,0.00317,0.00335,0.00355,0.00375,0.00385;
0.225,0.00307,0.00330,0.00355,0.00375,0.00385;
0.250,0.00299,0.00325,0.00352,0.00375,0.00385;
0.275,0.00389,0.00320,0.00350,0.00375,0.00388;
0.300,0.00280,0.00312,0.00345,0.00375,0.00385;
0.325,0.00270,0.00305,0.00342,0.00375,0.00385;
0.350,0.00260,0.00300,0.00340,0.00370,0.00390;
0.375,0.00250,0.00292,0.00335,0.00370,0.00390;
0.400,0.00240,0.00285,0.00330,0.00370,0.00390;
0.425,0.00228,0.00280,0.00325,0.00365,0.00390;
0.450,0.00318,0.00270,0.00320,0.00365,0.00390;
0.475,0.00207,0.00265,0.00318,0.00360,0.00390;
0.500,0.00195,0.00255,0.00310,0.00360,0.00385;
0.525,0.00185,0.00250,0.00305,0.00355,0.00385;
0.550,0.00175,0.00240,0.00300,0.00350,0.00385;
0.575,0.00165,0.00230,0.00295,0.00345,0.00380;
0.600,0.00155,0.00222,0.00285,0.00340,0.00375;
0.625,0.00144,0.00215,0.00280,0.00335,0.00370;
0.650,0.00132,0.00205,0.00270,0.00325,0.00365;
0.675,0.00122,0.00195,0.00260,0.00315,0.00360;
0.700,0.00111,0.00182,0.00250,0.00305,0.00350;
0.725,0.00100,0.00170,0.00235,0.00290,0.00335;
0.750,0.00089,0.00160,0.00215,0.00280,0.00325;
0.775,0.00075,0.00140,0.00200,0.00260,0.00310;
0.800,0.00062,0.00125,0.00180,0.00245,0.00295;
0.825,0.00043,0.00105,0.00160,0.00225,0.00275;
0.850,0.00026,0.00080,0.00140,0.00205,0.00250;
0.875,0.00009,0.00055,0.00115,0.00180,0.00230;
0.900,-0.00004,0.00030,0.00095,0.00165,0.00215;
0.925,-0.00025,0.00010,0.00076,0.00150,0.00196;
0.950,-0.00043,-0.00005,0.00075,0.00125,0.00175;
0.970,-0.00055,-0.00020,0.00045,0.00105,0.00158];
CLdEotable=CLdEogg(2:41,2:6)';
%Clspoilers
Clsp58gg=[-999,0,1,5,10,20,25,30;
-5, -0.030,-0.030,-0.028,-0.032,-0.042,-0.071,-0.081;
-4, -0.031,-0.031,-0.029,-0.034,-0.045,-0.075,-0.087;
-3, -0.032,-0.032,-0.030,-0.036,-0.049,-0.079,-0.091;
-2, -0.033,-0.033,-0.031,-0.038,-0.052,-0.082,-0.093;
-1, -0.034,-0.034,-0.033,-0.040,-0.056,-0.085,-0.098;
0, -0.035,-0.035,-0.034,-0.040,-0.058,-0.088,-0.100;
1, -0.036,-0.036,-0.035,-0.042,-0.061,-0.091,-0.101;
2, -0.036,-0.036,-0.035,-0.043,-0.063,-0.093,-0.102;
3, -0.037,-0.037,-0.036,-0.044,-0.063,-0.095,-0.103;
4, -0.035,-0.035,-0.034,-0.044,-0.064,-0.096,-0.102;
5, -0.037,-0.037,-0.036,-0.044,-0.062,-0.096,-0.101;
6, -0.038,-0.038,-0.036,-0.044,-0.060,-0.093,-0.100;
7, -0.038,-0.038,-0.035,-0.043,-0.058,-0.091,-0.098;
8, -0.038,-0.038,-0.034,-0.043,-0.058,-0.091,-0.097;
9, -0.038,-0.038,-0.034,-0.042,-0.059,-0.089,-0.098;
10, -0.038,-0.038,-0.034,-0.042,-0.061,-0.090,-0.099;
11, -0.039,-0.039,-0.035,-0.042,-0.062,-0.091,-0.099;
12, -0.038,-0.038,-0.036,-0.043,-0.064,-0.093,-0.101;
13, -0.036,-0.036,-0.037,-0.044,-0.065,-0.096,-0.104;
14, -0.034,-0.034,-0.038,-0.044,-0.066,-0.097,-0.107;
15, -0.030,-0.030,-0.039,-0.044,-0.066,-0.101,-0.110;
16, -0.026,-0.026,-0.040,-0.043,-0.065,-0.105,-0.113;
17, -0.020,-0.020,-0.040,-0.041,-0.063,-0.109,-0.117;
18, -0.015,-0.015,-0.039,-0.038,-0.060,-0.112,-0.121;
19, -0.010,-0.010,-0.038,-0.036,-0.056,-0.115,-0.124;
20, -0.008,-0.008,-0.036,-0.031,-0.048,-0.112,-0.118;
21, -0.007,-0.007,-0.034,-0.025,-0.036,-0.094,-0.095;
22, -0.006,-0.006,-0.030,-0.017,-0.026,-0.067,-0.067;
23, -0.006,-0.006,-0.025,-0.011,-0.020,-0.051,-0.051;
24, -0.006,-0.006,-0.019,-0.006,-0.014,-0.041,-0.041;
25, -0.006,-0.006,-0.013,-0.003,-0.009,-0.033,-0.033];
Clsp58table=Clsp58gg(2:32,2:8)';
Clsp91011234gg=[-999, 0, 1, 5, 10, 20, 25, 30;
-5, -0.029, -0.029,-0.040,-0.054,-0.071,-0.071,-0.076;
-4, -0.032, -0.032,-0.040,-0.053,-0.071,-0.079,-0.083;
-3, -0.035, -0.035,-0.040,-0.052,-0.071,-0.086,-0.091;
-2, -0.037, -0.037,-0.040,-0.052,-0.072,-0.092,-0.098;
-1, -0.039, -0.039,-0.040,-0.052,-0.072,-0.097,-0.103;
0, -0.041, -0.041,-0.041,-0.052,-0.073,-0.100,-0.108;
1, -0.044, -0.044,-0.042,-0.053,-0.074,-0.103,-0.111;
2, -0.046, -0.046,-0.042,-0.054,-0.075,-0.104,-0.113;
3, -0.047, -0.047,-0.042,-0.055,-0.077,-0.106,-0.116;
4, -0.049, -0.049,-0.043,-0.056,-0.080,-0.109,-0.120;
5, -0.050, -0.050,-0.043,-0.058,-0.084,-0.114,-0.125;
6, -0.052, -0.052,-0.043,-0.060,-0.089,-0.121,-0.132;
7, -0.053, -0.053,-0.044,-0.062,-0.094,-0.129,-0.140;
8, -0.054, -0.054,-0.045,-0.064,-0.099,-0.136,-0.146;
9, -0.054, -0.054,-0.046,-0.067,-0.103,-0.141,-0.151;
10, -0.054, -0.054,-0.046,-0.070,-0.108,-0.145,-0.153;
11, -0.054, -0.054,-0.048,-0.074,-0.112,-0.147,-0.154;
12, -0.053, -0.053,-0.051,-0.077,-0.116,-0.148,-0.154;
13, -0.052, -0.054,-0.051,-0.080,-0.120,-0.148,-0.153;
14, -0.049, -0.050,-0.054,-0.083,-0.122,-0.146,-0.150;
15, -0.038, -0.047,-0.058,-0.086,-0.124,-0.142,-0.146;
16, -0.024, -0.038,-0.061,-0.088,-0.125,-0.137,-0.140;
17, -0.015, -0.029,-0.063,-0.090,-0.126,-0.130,-0.133;
18, -0.011, -0.021,-0.064,-0.090,-0.126,-0.121,-0.124;
19, -0.009, -0.017,-0.064,-0.089,-0.123,-0.110,-0.112;
20, -0.007, -0.014,-0.063,-0.087,-0.120,-0.098,-0.101;
21, -0.006, -0.012,-0.062,-0.084,-0.113,-0.084,-0.088;
22, -0.005, -0.010,-0.058,-0.080,-0.105,-0.070,-0.073;
23, -0.004, -0.008,-0.055,-0.074,-0.097,-0.055,-0.057;
24, -0.004, -0.007,-0.052,-0.070,-0.092,-0.045,-0.047;
25, -0.004, -0.006,-0.050,-0.068,-0.089,-0.040,-0.042;];
Clsp91011234table=Clsp91011234gg(2:32,2:8)';
Clsp67gg=[-999,0,1,5,10,20,25,30;
-5,-0.052,-0.052,-0.063,-0.085,-0.0131,-0.0180,-0.0195;
-4,-0.054,-0.054,-0.065,-0.086,-0.0132,-0.0180,-0.0195;
-3,-0.054,-0.054,-0.067,-0.087,-0.0133,-0.0180,-0.0194;
-2,-0.056,-0.056,-0.068,-0.088,-0.0134,-0.0180,-0.0194;
-1,-0.058,-0.058,-0.070,-0.089,-0.0134,-0.0180,-0.0194;
0,-0.060,-0.060,-0.074,-0.090,-0.0135,-0.0178,-0.0190;
1,-0.061,-0.061,-0.073,-0.092,-0.0135,-0.0178,-0.0188;
2,-0.062,-0.062,-0.074,-0.094,-0.036,-0.0178,-0.0185;
3,-0.063,-0.063,-0.076,-0.092,-0.0136,-0.0178,-0.0182;
4,-0.064,-0.064,-0.076,-0.02,-0.0136,-0.0177,-0.0180;
5,-0.064,-0.064,-0.078,-0.093,-0.0137,-0.0177,-0.0177;
6,-0.065,-0.065,-0.080,-0.093,-0.0137,-0.0176,-0.0174;
7,-0.065,-0.065,-0.080,-0.094,-0.0136,-0.0175,-0.0170;
8,-0.065,-0.065,-0.080,-0.094,-0.0138,-0.0174,-0.0168;
9,-0.064,-0.064,-0.080,-0.094,-0.0136,-0.0172,-0.0164;
10,-0.064,-0.064,-0.080,-0.093,-0.0135,-0.0171,-0.0160;
11,-0.063,-0.063,-0.080,-0.092,-0.0134,-0.0170,-0.0156;
12,-0.062,-0.062,-0.080,-0.092,-0.0132,-0.0168,-0.0152;
13,-0.060,-0.060,-0.079,-0.091,-0.0129,-0.0165,-0.0148;
14,-0.058,-0.058,-0.078,-0.088,-0.0126,-0.0162,-0.0142;
15,-0.054,-0.054,-0.074,-0.086,-0.0120,-0.0157,-0.0137;
16,-0.048,-0.048,-0.071,-0.082,-0.0114,-0.0150,-0.0130;
17,-0.037,-0.037,-0.064,-0.076,-0.0106,-0.0143,-0.0123;
18,-0.025,-0.025,-0.056,-0.067,-0.098,-0.0132,-0.0115;
19,-0.010,-0.010,-0.046,-0.058,-0.089,-0.0121,0.0105;
20,0.004,0.004,-0.034,-0.048,-0.078,-0.0110,-0.096;
21,0.012,0.012,-0.020,-0.037,-0.066,-0.096,-0.086;
22,0.020,0.020,-0.009,-0.023,-0.053,-0.086,-0.074;
23,0.027,0.027,0.004,-0.012,-0.040,-0.073,-0.063;
24,0.032,0.032,0.012,0.000,-0.027,-0.060,-0.052;
25,0.037,0.037,0.021,0.010,-0.014,-0.049,-0.039;];
Clsp67table=Clsp67gg(2:32,2:8)';
Ksp20=[0.29,0.30,0.40,0.44,0.50,0.53,0.58];
KGEBggtable=[-999,0,1,5,10,20,25,30;
0.0,0.178,0.178,0.163,0.140,0.060,0.000,-0.045;
0.5,0.175,0.175,0.161,0.135,0.054,-0.006,-0.052;
1.0,0.173,0.173,0.156,0.130,0.052,-0.010,-0.059;
1.5,0.170,0.170,0.151,0.120,0.049,-0.015,-0.069;
2.0,0.166,0.166,0.137,0.118,0.047,-0.020,-0.076;
2.5,0.161,0.161,0.132,0.110,0.042,-0.027,-0.083;
3.0,0.156,0.156,0.127,0.105,0.040,-0.035,-0.093;
3.5,0.149,0.149,0.120,0.100,0.035,-0.039,-0.100;
4.0,0.144,0.144,0.115,0.090,0.033,-0.044,-0.110;
4.5,0.137,0.137,0.108,0.082,0.028,-0.051,-0.119;
5.0,0.132,0.132,0.101,0.078,0.023,-0.059,-0.126;
5.5,0.128,0.128,0.096,0.070,0.018,-0.066,-0.136;
6.0,0.120,0.120,0.089,0.063,0.013,-0.071,-0.146;
6.5,0.116,0.116,0.082,0.059,0.006,-0.080,-0.158;
7.0,0.111,0.111,0.077,0.052,0.001,-0.088,-0.165;
7.5,0.104,0.104,0.073,0.044,-0.010,-0.095,-0.174;
8.0,0.096,0.096,0.065,0.037,-0.016,-0.105,-0.186;
8.5,0.089,0.089,0.056,0.030,-0.020,-0.114,-0.198;
9.0,0.085,0.085,0.051,0.021,-0.028,-0.121,-0.210;
9.5,0.080,0.080,0.042,0.014,-0.037,-0.131,-0.220;
10.0,0.073,0.073,0.034,0.004,-0.047,-0.143,-0.232;
10.5,0.065,0.065,0.025,0.000,-0.057,-0.155,-0.244;
11.0,0.058,0.058,0.018,-0.010,-0.066,-0.165,-0.258;
11.5,0.051,0.051,0.008,-0.020,-0.078,-0.177,-0.270;
12.0,0.046,0.046,0.001,-0.029,-0.090,-0.191,-0.285;
12.5,0.039,0.039,-0.010,-0.036,-0.107,-0.203,-0.299;
13.0,0.032,0.032,-0.020,-0.048,-0.119,-0.223,-0.316;
13.5,0.025,0.025,-0.030,-0.057,-0.136,-0.237,-0.330;
14.0,0.000,0.000,-0.037,-0.067,-0.151,-0.252,-0.347;];
KGEBtable=KGEBggtable(2:30,2:8)';
Alphawdptable=0:0.5:14;
%CLgear
CLgeargg=[-999,0,1,5,10,20,25,30;
-5,0.0465,0.0465,0.0465,0.0370,0.0420,-0.0060,-0.0060;
-4,0.0450,0.0450,0.0450,0.0355,0.0390,-0.0100,-0.0100;
-3,0.0440,0.0440,0.0435,0.0335,0.0360,-0.0150,-0.0150;
-2,0.0430,0.0430,0.0420,0.0320,0.0335,-0.0190,-0.0190;
-1,0.0425,0.0425,0.0410,0.0305,0.0310,-0.0215,-0.0215;
0,0.0415,0.0415,0.0400,0.0290,0.0290,-0.0250,-0.0250;
1,0.0400,0.0400,0.0380,0.0270,0.0260,-0.0270,-0.0270;
2,0.0390,0.0390,0.0370,0.0255,0.0240,-0.0285,-0.0285;
3,0.0380,0.0380,0.0350,0.0240,0.0220,-0.0290,-0.0290;
4,0.0370,0.0370,0.0340,0.0225,0.0200,-0.0300,-0.0300;
5,0.0360,0.0360,0.0320,0.0210,0.0180,-0.0300,-0.0300;
6,0.0350,0.0350,0.0310,0.0190,0.0160,-0.0300,-0.0300;
7,0.0335,0.0335,0.0290,0.0175,0.0145,-0.0300,-0.0300;
8,0.0325,0.0325,0.0280,0.0160,0.0125,-0.0295,-0.0295;
9,0.0310,0.0310,0.0260,0.0145,0.0105,-0.0290,-0.0290;
10,0.0300,0.0300,0.0245,0.0130,0.0090,-0.0285,-0.0285;
11,0.0285,0.0285,0.0230,0.0115,0.0070,-0.0280,-0.0280;
12,0.0275,0.0275,0.0215,0.0100,0.0060,-0.0270,-0.0270;
13,0.0260,0.0260,0.0195,0.0085,0.0045,-0.0265,-0.0265;
14,0.0245,0.0245,0.0185,0.0070,0.0030,-0.0260,-0.0260;
15,0.0235,0.0235,0.0160,0.0060,0.0015,-0.0250,-0.0250;
16,0.0210,0.0210,0.0160,0.0045,0.0000,-0.0245,-0.0245;
17,0.0190,0.0190,0.0125,0.0035,-0.0010,-0.0240,-0.0240;
18,0.0175,0.0175,0.0105,0.0020,-0.0020,-0.0240,-0.0240;
19,0.0155,0.0155,0.0080,0.0001,-0.0040,-0.0240,-0.0240;
20,0.0135,0.0135,0.0060,-0.0005,-0.0050,-0.0245,-0.0245;
21,0.0120,0.0120,0.0040,-0.0015,-0.0060,-0.0250,-0.0250;
22,0.0090,0.0090,0.0020,-0.0030,-0.0075,-0.0255,-0.0255;
23,0.0070,0.0070,-0.0005,-0.0040,-0.0090,-0.0265,-0.0265;
24,0.0045,0.0045,-0.0030,-0.0050,-0.0100,-0.0270,-0.0270;
25,0.0020,0.0020,-0.0050,-0.0060,-0.0110,-0.0280,-0.0280];
CLgeartable=CLgeargg(2:32,2:8)';
KCLgearM=[1,1,1,1,1,1,1,1,1,1,1,1,1,1.05,1.12,1.18,1.25,1.3,1.4,1.45,1.55,1.63,1.72,1.8,1.9,2,2.13,2.25,2.4,2.6,2.85,3.15,3.5,3.85,3.85,3.85,3.85,3.85,3.85,3.85];
%CLgroundeffect
Kgatable=[1.0000,0.9900,0.9700,0.9400,0.8900,0.8300,0.7500,0.6600,0.5700,0.4800,0.3700];
Hgatable=[0,5,10,15,20,25,30,35,40,45,50];
CLgagg=[-999,0,1,5,10,20,25,30;
0, -0.075,-0.060,-0.045, 0.040,0.075,0.245,0.225;
0.5,-0.070,-0.050,-0.037, 0.050,0.105,0.255,0.237;
1.0,-0.067,-0.045,-0.030, 0.057,0.118,0.257,0.233;
1.5,-0.053,-0.037,-0.015, 0.065,0.120,0.259,0.235;
2.0,-0.045,-0.025,-0.005, 0.075,0.130,0.265,0.235;
2.5,-0.035,-0.015, 0.005, 0.085,0.143,0.270,0.235;
3.0,-0.025,-0.007, 0.015, 0.100,0.157,0.275,0.235;
3.5,-0.015, 0.007, 0.027, 0.110,0.160,0.276,0.230;
4.0,-0.005, 0.018, 0.035, 0.120,0.180,0.280,0.228;
4.5, 0.010, 0.030, 0.050, 0.130,0.190,0.280,0.225;
5.0, 0.023, 0.045, 0.065, 0.145,0.200,0.280,0.220;
5.5, 0.035, 0.060, 0.078, 0.158,0.210,0.277,0.215;
6.0, 0.050, 0.075, 0.095, 0.170,0.215,0.275,0.208;
6.5, 0.065, 0.090, 0.108, 0.178,0.220,0.272,0.200;
7.0, 0.080, 0.109, 0.123, 0.187,0.223,0.270,0.187;
7.5, 0.095, 0.123, 0.137, 0.200,0.225,0.265,0.180;
8.0, 0.110, 0.135, 0.150, 0.207,0.225,0.260,0.168;
8.5, 0.120, 0.147, 0.163, 0.215,0.225,0.250,0.157;
9.0, 0.130, 0.159, 0.173, 0.220,0.225,0.243,0.145;
9.5, 0.135, 0.168, 0.180, 0.225,0.225,0.235,0.130;
10.0,0.148, 0.173, 0.190, 0.230,0.225,0.225,0.115;
10.5,0.155, 0.177, 0.196, 0.235,0.224,0.215,0.100;
11.0,0.160, 0.183, 0.200, 0.235,0.222,0.200,0.085;
11.5,0.165, 0.185, 0.203, 0.235,0.220,0.185,0.070;
12.0,0.168, 0.189, 0.205, 0.235,0.215,0.175,0.053;
12.5,0.070, 0.190, 0.205, 0.235,0.213,0.157,0.033;
13.0,0.172, 0.190, 0.205, 0.234,0.205,0.143,0.018;
13.5,0.175, 0.192, 0.203, 0.233,0.200,0.125,0.000;
14.0,0.175, 0.192, 0.200, 0.230,0.192,0.102,-0.020];
CLgatable=CLgagg(2:30,2:8)'
%////////////////////////////////////////////////////////////////
%俯仰力矩系数
%CMBasictable
CMBasicgg=[-999,0,1,5,10,20,25,30;
-5,0.2780,0.235,0.266,0.280,0.196,0.149,0.117;
-4,0.2508,0.213,0.243,0.259,0.172,0.129,0.095;
-3,0.2272,0.188,0.215,0.229,0.147,0.105,0.068;
-2,0.2017,0.162,0.188,0.202,0.0121,0.083,0.040;
-1,0.1742,0.139,0.160,0.174,0.091,0.060,0.012;
0,0.1487,0.113,0.131,0.145,0.064,0.034,-0.015;
1,0.1232,0.090,0.106,0.119,0.040,0.012,-0.043;
2,0.0977,0.066,0.078,0.092,0.012,-0.013,-0.075;
3,0.0722,0.039,0.051,0.064,-0.015,-0.039,-0.104;
4,0.0467,0.013,0.021,0.035,-0.041,-0.063,-0.134;
5,0.0212,-0.010,-0.006,0.007,-0.069,-0.088,-0.167;
6,-0.0024,-0.034,-0.032,-0.020,-0.094,-0.114,-0.197;
7,-0.0298,-0.057,-0.061,-0.047,-0.122,-0.140,-0.231;
8,-0.0553,-0.081,-0.089,-0.077,-0.150,-0.167,-0.266;
9,-0.0828,-0.102,-0.118,-0.104,-0.177,-0.195,-0.298;
10,-0.1083,-0.122,-0.146,-0.134,-0.203,-0.223,-0.329;
11,-0.1299,-0.138,-0.171,-0.161,-0.231,-0.252,-0.361;
12,-0.1515,-0.151,-0.195,-0.193,-0.254,-0.278,-0.387;
13,-0.1731,-0.163,-0.218,-0.226,-0.280,-0.304,-0.408;
14,-0.1927,-0.175,-0.242,-0.261,-0.306,-0.329,-0.428;
15,-0.2123,-0.191,-0.265,-0.281,-0.327,-0.351,-0.442;
16,-0.2319,-0.208,-0.281,-0.287,-0.347,-0.375,-0.456;
17,-0.2574,-0.234,-0.295,-0.293,-0.365,-0.389,-0.458;
18,-0.3006,-0.259,-0.303,-0.297,-0.359,-0.351,-0.434;
19,-0.3693,-0.289,-0.316,-0.303,-0.331,-0.315,-0.390;
20,-0.4340,-0.320,-0.365,-0.312,-0.308,-0.306,-0.361;
21,-0.4949,-0.356,-0.442,-0.332,-0.304,-0.335,-0.375;
22,-0.5576,-0.393,-0.532,-0.391,-0.349,-0.400,-0.475;
23,-0.6204,-0.432,-0.636,-0.475,-0.428,-0.500,-0.625;
24,-0.6871,-0.479,-0.750,-0.565,-0.550,-0.625,-0.805;
25,-0.7500,-0.525,-0.850,-0.654,-0.659,-0.790,-1.100;];
CMBasictable=CMBasicgg(2:32,2:8)';
%CMAlphadottable
CMAlphadotgg=[-999,0,20000,40000;
0, -3.40,-3.40, -3.40;
0.025,-3.45,-3.45, -3.44;
0.050,-3.42,-3.45, -3.46;
0.075,-3.40,-3.45, -3.49;
0.100,-3.37,-3.45, -3.52;
0.125,-3.32,-3.45, -3.54;
0.150,-3.32,-3.45, -3.57;
0.175,-3.32,-3.45, -3.60;
0.200,-3.29,-3.45, -3.62;
0.225,-3.26,-3.45, -3.67;
0.250,-3.24,-3.48, -3.70;
0.275,-3.21,-3.48, -3.73;
0.300,-3.18,-3.51, -3.75;
0.325,-3.16,-3.53, -3.80;
0.350,-3.13,-3.53, -3.80;
0.375,-3.11,-3.56, -3.83;
0.400,-3.08,-3.56, -3.86;
0.425,-3.08,-3.58, -3.91;
0.450,-3.08,-3.58, -3.94;
0.475,-3.08,-3.61, -3.99;
0.500,-3.08,-3.61, -4.04;
0.525,-3.08,-3.61, -4.09;
0.550,-3.08,-3.64, -4.15;
0.575,-3.08,-3.69, -4.22;
0.600,-3.11,-3.77, -4.33;
0.625,-3.16,-3.87, -4.43;
0.650,-3.21,-3.98, -4.57;
0.675,-3.26,-4.13, -4.78;
0.700,-3.34,-4.34, -5.01;
0.725,-3.42,-4.58, -5.33;
0.750,-3.58,-4.87, -5.69;
0.775,-3.79,-5.21, -6.17;
0.800,-3.97,-5.49, -6.69;
0.825,-4.11,-5.73, -7.06;
0.850,-4.16,-6.02, -7.61;
0.875,-4.39,-6.54, -8.53;
0.900,-4.78,-6.93, -8.87;
0.925,-4.42,-5.52, -7.43;
0.950,-3.08,-4.66, -6.06;
0.970,-3.15,-4.70, -6.10;];
CMAlphadottable=CMAlphadotgg(2:41,2:4)';
%CMqdottable
CMqdotgg=[
-999,0,20000,40000;
0, -22.60,-22.60,-22.60;
0.025,-22.57,-22.58,-22.62;
0.050,-22.49,-22.55,-22.60;
0.075,-22.44,-22.53,-22.62;
0.100,-22.31,-22.53,-22.62;
0.125,-22.15,-22.47,-22.60;
0.150,-21.94,-22.42,-22.60;
0.175,-21.68,-22.37,-22.60;
0.200,-21.42,-22.29,-22.57;
0.225,-21.14,-22.21,-22.57;
0.250,-20.88,-22.09,-22.60;
0.275,-20.59,-21.98,-22.57;
0.300,-20.33,-21.85,-22.57;
0.325,-20.07,-21.72,-22.57;
0.350,-19.83,-21.59,-22.57;
0.375,-19.57,-21.44,-22.57;
0.400,-19.37,-21.31,-22.60;
0.425,-19.18,-21.13,-22.60;
0.450,-18.95,-20.97,-22.60;
0.475,-18.77,-20.81,-22.57;
0.500,-18.53,-20.66,-22.57;
0.525,-18.35,-20.50,-22.62;
0.550,-18.17,-20.43,-22.62;
0.575,-17.96,-20.32,-22.62;
0.600,-17.80,-20.30,-22.65;
0.625,-17.60,-20.30,-22.65;
0.650,-17.39,-20.32,-22.70;
0.675,-17.18,-20.30,-22.78;
0.700,-16.97,-20.35,-22.99;
0.725,-16.79,-20.37,-23.18;
0.750,-16.58,-20.40,-23.41;
0.775,-16.37,-20.43,-23.73;
0.800,-16.14,-20.50,-24.12;
0.825,-15.88,-20.58,-24.62;
0.850,-15.57,-20.63,-25.14;
0.875,-15.18,-20.56,-25.37;
0.900,-14.60,-20.40,-25.25;
0.925,-13.82,-19.50,-24.40;
0.950,-13.02,-18.37,-23.00;
0.970,-12.20,-17.30,-21.40;];
CMqtable=CMqdotgg(2:41,2:4)';
%CMdEi
CMdEigg=[-999,0,10000,20000,30000,40000;
0, -0.0130,-0.0130,-0.0130,-0.0130,-0.0130;
0.025,-0.0129,-0.0130,-0.0130,-0.0130,-0.0130;
0.050,-0.0128,-0.0129,-0.0130,-0.0130,-0.0130;
0.075,-0.0127,-0.0129,-0.0130,-0.0130,-0.0131;
0.100,-0.0126,-0.0128,-0.0130,-0.0131,-0.0131;
0.125,-0.0125,-0.0127,-0.0130,-0.0131,-0.0131;
0.150,-0.0123,-0.0127,-0.0129,-0.0130,-0.0131;
0.175,-0.0121,-0.0125,-0.0128,-0.0131,-0.0132;
0.200,-0.0120,-0.0124,-0.0128,-0.0131,-0.0132;
0.225,-0.0118,-0.0123,-0.0128,-0.0132,-0.0133;
0.250,-0.0116,-0.0123,-0.0128,-0.0132,-0.0134;
0.275,-0.0114,-0.0121,-0.0128,-0.0132,-0.0135;
0.300,-0.0112,-0.0120,-0.0127,-0.0132,-0.0136;
0.325,-0.0109,-0.0119,-0.0127,-0.0132,-0.0137;
0.350,-0.0106,-0.0117,-0.0127,-0.0133,-0.0137;
0.375,-0.0104,-0.0117,-0.0126,-0.0133,-0.0138;
0.400,-0.0101,-0.0116,-0.0126,-0.0133,-0.0139;
0.425,-0.0098,-0.0114,-0.0125,-0.0133,-0.0139;
0.450,-0.0095,-0.0113,-0.0124,-0.0133,-0.0140;
0.475,-0.0093,-0.0111,-0.0124,-0.0133,-0.0141;
0.500,-0.0090,-0.0110,-0.0123,-0.0133,-0.0141;
0.525,-0.0087,-0.0108,-0.0123,-0.0133,-0.0142;
0.550,-0.0085,-0.0107,-0.0122,-0.0133,-0.0142;
0.575,-0.0082,-0.0105,-0.0122,-0.0132,-0.0142;
0.600,-0.0079,-0.0103,-0.0120,-0.0131,-0.0142;
0.625,-0.0077,-0.0101,-0.0120,-0.0131,-0.0142;
0.650,-0.0075,-0.0099,-0.0119,-0.0131,-0.0141;
0.675,-0.0072,-0.0098,-0.0117,-0.0130,-0.0141;
0.700,-0.0070,-0.0097,-0.0117,-0.0130,-0.0140;
0.725,-0.0069,-0.0096,-0.0116,-0.0129,-0.0140;
0.750,-0.0068,-0.0095,-0.0115,-0.0128,-0.0140;
0.775,-0.0067,-0.0094,-0.0114,-0.0128,-0.0139;
0.800,-0.0066,-0.0094,-0.0113,-0.0127,-0.0138;
0.825,-0.0065,-0.0093,-0.0113,-0.0126,-0.0137;
0.850,-0.0064,-0.0092,-0.0111,-0.0124,-0.0135;
0.875,-0.0063,-0.0090,-0.0107,-0.0122,-0.0132;
0.900,-0.0062,-0.0084,-0.0100,-0.0115,-0.0125;
0.925,-0.0061,-0.0075,-0.0089,-0.0103,-0.0114;
0.950,-0.0060,-0.0067,-0.0078,-0.0090,-0.0101;
0.970,-0.0060,-0.0060,-0.0070,-0.0080,-0.0088];
CMdEitable=CMdEigg(2:41,2:6)';
%CMgear
CMgeargg=[-999,0,1,5,10,20,25,30;
-5,0.010,0.010,0.005,0.002,-0.002,-0.015,0.015;
-4,0.010,0.010,0.005,0.002,-0.003,-0.015,0.015;
-3,0.010,0.010,0.005,0.001,-0.004,-0.015,0.015;
-2,0.010,0.010,0.005,0.000,-0.005,-0.015,0.015;
-1,0.010,0.010,0.005,-0.001,-0.005,-0.014,-0.014;
0,0.009,0.009,0.004,-0.001,-0.006,-0.013,-0.013;
1,0.009,0.009,0.004,-0.002,-0.006,-0.013,-0.013;
2,0.008,0.008,0.004,-0.002,-0.006,-0.013,-0.013;
3,0.008,0.008,0.003,-0.003,-0.007,-0.013,-0.013;
4,0.008,0.008,0.003,-0.003,-0.007,-0.013,-0.013;
5,0.007,0.007,0.002,-0.003,-0.007,-0.013,-0.013;
6,0.007,0.007,0.002,-0.004,-0.008,-0.013,-0.013;
7,0.006,0.006,0.001,-0.004,-0.008,-0.013,-0.013;
8,0.006,0.006,0.001,-0.005,-0.009,-0.013,-0.013;
9,0.005,0.005,0.000,-0.006,-0.009,-0.013,-0.013;
10,0.005,0.005,-0.001,-0.007,-0.010,-0.014,-0.014;
11,0.004,0.004,-0.002,-0.008,-0.011,-0.014,-0.014;
12,0.004,0.004,-0.003,-0.009,-0.011,-0.015,-0.015;
13,0.003,0.003,-0.004,-0.010,-0.012,-0.016,-0.016;
14,0.002,0.002,-0.005,-0.011,-0.012,-0.017,-0.017;
15,0.001,0.001,-0.006,-0.011,-0.013,-0.019,-0.019;
16,0.000,0.000,-0.007,-0.011,-0.013,-0.021,-0.021;
17,-0.001,-0.001,-0.008,-0.012,-0.013,-0.024,-0.024;
18,-0.002,-0.002,-0.008,-0.012,-0.014,-0.028,-0.028;
19,-0.002,-0.002,-0.007,-0.011,-0.014,-0.033,-0.033;
20,-0.002,-0.002,-0.005,-0.011,-0.014,-0.040,-0.040;
21,0.000,0.000,-0.003,-0.011,-0.014,-0.049,-0.049;
22,0.005,0.005,0.002,-0.010,-0.014,-0.060,-0.060;
23,0.011,0.011,0.007,-0.009,-0.014,-0.062,-0.062;
24,0.017,0.017,0.012,-0.008,-0.014,-0.044,-0.044;
25,0.022,0.022,0.017,-0.075,-0.013,-0.030,-0.030];
CMgeartable=CMgeargg(2:32,2:8)';
KCMgearM=[1,1,1,1,1,1,1,1,1,1,1,1,1,1.013,1.039,1.052,1.078,1.091,1.1,1.15,1.208,1.26,1.299,1.364,1.442,1.519,1.61,1.7,1.818,1.9,2.052,2.195,2.364,2.5,2.5,2.5,2.5,2.5,2.5,2.5;];
%CMdEo
CMdEogg=[-999,0,10000,20000,30000,40000;
0, -0.0150,-0.0150,-0.0150,-0.0150,-0.0150;
0.025,-0.0148,-0.0149,-0.0149,-0.0151,-0.0151;
0.050,-0.0145,-0.0147,-0.0148,-0.0150,-0.0151;
0.075,-0.0142,-0.0144,-0.0146,-0.0150,-0.0151;
0.100,-0.0139,-0.0142,-0.0145,-0.0149,-0.0151;
0.125,-0.0136,-0.0140,-0.0144,-0.0149,-0.0151;
0.150,-0.0133,-0.0138,-0.0143,-0.0149,-0.0151;
0.175,-0.0129,-0.0135,-0.0142,-0.0149,-0.0151;
0.200,-0.0126,-0.0133,-0.0141,-0.0149,-0.0151;
0.225,-0.0122,-0.0131,-0.0140,-0.0149,-0.0152;
0.250,-0.0118,-0.0129,-0.0139,-0.0149,-0.0153;
0.275,-0.0115,-0.0126,-0.0138,-0.0149,-0.0153;
0.300,-0.0111,-0.0123,-0.0137,-0.0148,-0.0154;
0.325,-0.0107,-0.0121,-0.0136,-0.0148,-0.0154;
0.350,-0.0103,-0.0118,-0.0134,-0.0148,-0.0154;
0.375,-0.0099,-0.0116,-0.0133,-0.0147,-0.0154;
0.400,-0.0095,-0.0113,-0.0132,-0.0146,-0.0154;
0.425,-0.0090,-0.0110,-0.0130,-0.0145,-0.0154;
0.450,-0.0086,-0.0107,-0.0128,-0.0144,-0.0154;
0.475,-0.0082,-0.0105,-0.0126,-0.0143,-0.0154;
0.500,-0.0077,-0.0102,-0.0123,-0.0142,-0.0153;
0.525,-0.0073,-0.0098,-0.0120,-0.0140,-0.0152;
0.550,-0.0069,-0.0095,-0.0118,-0.0138,-0.0151;
0.575,-0.0065,-0.0092,-0.0115,-0.0136,-0.0150;
0.600,-0.0061,-0.0088,-0.0112,-0.0134,-0.0148;
0.625,-0.0057,-0.0084,-0.0109,-0.0131,-0.0146;
0.650,-0.0053,-0.0081,-0.0105,-0.0128,-0.0144;
0.675,-0.0049,-0.0077,-0.0101,-0.0125,-0.0141;
0.700,-0.0045,-0.0072,-0.0097,-0.0120,-0.0137;
0.725,-0.0041,-0.0067,-0.0092,-0.0116,-0.0133;
0.750,-0.0035,-0.0062,-0.0086,-0.0111,-0.0129;
0.775,-0.0031,-0.0056,-0.0080,-0.0105,-0.0124;
0.800,-0.0025,-0.0048,-0.0071,-0.0099,-0.0117;
0.825,-0.0018,-0.0041,-0.0063,-0.0091,-0.0109;
0.850,-0.0014,-0.0031,-0.0054,-0.0083,-0.0100;
0.875,-0.0006,-0.0022,-0.0045,-0.0074,-0.0092;
0.900,-0.0002,-0.0013,-0.0037,-0.0065,-0.0085;
0.925,0.0004,-0.0005,-0.0031,-0.0058,-0.0078;
0.950,0.0010,-0.0002,-0.0025,-0.0051,-0.0071;
0.970,0.0016,0.0008,-0.0018,-0.0044,-0.0062];
CMdEotable=CMdEogg(2:41,2:6)';
%CMgroundeffect
CMgatablegg=[-999,0,1,5,10,20,25,30;
0, 0.038,0.038,0.005,-0.015,-0.045,-0.081,-0.093;
0.5,0.035,0.035,0.003,-0.016,-0.051,-0.084,-0.097;
1.0,0.034,0.034,0.001,-0.019,-0.052,-0.088,-0.103;
1.5,0.031,0.031,-0.001,-0.020,-0.055,-0.092,-0.109;
2.0,0.028,0.028,-0.003,-0.022,-0.057,-0.097,-0.114;
2.5,0.026,0.026,-0.007,-0.024,-0.059,-0.101,-0.119;
3.0,0.026,0.026,-0.011,-0.028,-0.063,-0.106,-0.125
3.5,0.018,0.018,-0.015,-0.032,-0.066,-0.111,-0.130;
4.0,0.014,0.014,-0.020,-0.035,-0.070,-0.117,-0.134;
4.5,0.007,0.007,-0.026,-0.040,-0.076,-0.122,-0.138;
5.0,0.001,0.001,-0.031,-0.045,-0.081,-0.126,-0.140;
5.5,-0.006,-0.006,-0.037,-0.052,-0.086,-0.131,-0.142;
6.0,-0.012,-0.012,-0.044,-0.057,-0.092,-0.136,-0.142;
6.5,-0.019,-0.019,-0.051,-0.064,-0.098,-0.140,-0.142;
7.0,-0.024,-0.024,-0.057,-0.070,-0.103,-0.146,-0.143;
7.5,-0.031,-0.031,-0.063,-0.077,-0.109,-0.150,-0.143;
8.0,-0.037,-0.037,-0.069,-0.084,-0.115,-0.154,-0.144;
8.5,-0.043,-0.043,-0.074,-0.092,-0.121,-0.158,-0.146;
9.0,-0.048,-0.048,-0.081,-0.098,-0.127,-0.163,-0.147;
9.5,-0.053,-0.053,-0.088,-0.105,-0.134,-0.167,-0.150;
10.0,-0.060,-0.060,-0.094,-0.113,-0.146,-0.171,-0.154;
10.5,-0.066,-0.066,-0.099,-0.121,-0.147,-0.175,-0.158;
11.0,-0.074,-0.074,-0.106,-0.128,-0.152,-0.181,-0.163;
11.5,-0.082,-0.082,-0.114,-0.136,-0.159,-0.185,-0.165;
12.0,-0.092,-0.092,-0.125,-0.146,-0.165,-0.188,-0.165;
12.5,-0.101,-0.101,-0.136,-0.156,-0.171,-0.189,-0.163;
13.0,-0.113,-0.113,-0.147,-0.167,-0.176,-0.188,-0.159;
13.5,-0.125,-0.125,-0.159,-0.180,-0.179,-0.187,-0.156;
14.0,-0.136,-0.136,-0.171,-0.192,-0.183,-0.181,-0.147];
CMgatable=CMgatablegg(2:30,2:8)';
%///////////////////////////////////////////////////////////////////
%阻力系数
%CDBasic
CDBasicgg=[-999,0,1,5,10,20,25,30;
-5,0.090,0.090,0.120,0.100,0.100,0.090,0.135;
-4,0.076,0.076,0.105,0.090,0.090,0.085,0.126;
-3,0.063,0.063,0.090,0.076,0.081,0.078,0.120;
-2,0.052,0.052,0.076,0.066,0.070,0.072,0.116;
-1,0.043,0.043,0.065,0.055,0.065,0.070,0.115;
0, 0.035,0.035,0.055,0.050,0.060,0.070,0.117;
1, 0.030,0.030,0.048,0.045,0.058,0.071,0.121;
2, 0.028,0.028,0.042,0.042,0.058,0.075,0.128;
3, 0.027,0.027,0.038,0.041,0.061,0.080,0.138;
4, 0.029,0.029,0.039,0.044,0.065,0.088,0.150;
5, 0.031,0.031,0.040,0.048,0.073,0.098,0.162;
6, 0.035,0.035,0.045,0.053,0.082,0.108,0.177;
7, 0.040,0.040,0.050,0.060,0.092,0.120,0.196;
8, 0.045,0.045,0.058,0.070,0.104,0.135,0.212;
9, 0.054,0.054,0.068,0.080,0.116,0.150,0.232;
10,0.061,0.061,0.078,0.093,0.130,0.166,0.250;
11,0.071,0.071,0.090,0.108,0.147,0.185,0.270;
12,0.084,0.084,0.104,0.120,0.160,0.205,0.293;
13,0.098,0.098,0.110,0.138,0.178,0.230,0.314;
14,0.113,0.113,0.137,0.155,0.197,0.255,0.338;
15,0.135,0.135,0.155,0.172,0.215,0.280,0.360;
16,0.153,0.153,0.170,0.192,0.237,0.310,0.385;
17,0.190,0.190,0.190,0.213,0.260,0.340,0.413;
18,0.225,0.225,0.210,0.235,0.280,0.375,0.450;
19,0.265,0.265,0.230,0.260,0.300,0.410,0.490;
20,0.295,0.295,0.258,0.285,0.325,0.440,0.530;
21,0.325,0.325,0.290,0.320,0.360,0.490,0.570;
22,0.355,0.355,0.335,0.350,0.395,0.525,0.600;
23,0.390,0.390,0.375,0.410,0.440,0.560,0.633;
24,0.415,0.415,0.410,0.430,0.478,0.595,0.660;
25,0.445,0.445,0.440,0.458,0.510,0.624,0.688];
CDBasictable=CDBasicgg(2:32,2:8)';
%CDbata
Batadgtable=[-15:1:15];
CDBatadggg=[-999,0,1,5,10,20,25,30;
-15,-0.0225,-0.0225,-0.0225,-0.0225,-0.0225,-0.0280,-0.0310;
-14,-0.0205,-0.0205,-0.0205,-0.0205,-0.0205,-0.0255,-0.0280;
-13,-0.0180,-0.0180,-0.0180,-0.0180,-0.0180,-0.0240,-0.0255;
-12,-0.0160,-0.0160,-0.0160,-0.0160,-0.0160,-0.0200,-0.0225;
-11,-0.0140,-0.0140,-0.0140,-0.0140,-0.0140,-0.0175,-0.0200;
-10,-0.0120,-0.0120,-0.0120,-0.0120,-0.0120,-0.0150,-0.0154;
-9,-0.0100,-0.0100,-0.0100,-0.0100,-0.0100,-0.0125,-0.0145;
-8,-0.0080,-0.0080,-0.0080,-0.0080,-0.0080,-0.0100,-0.0120;
-7,-0.0065,-0.0065,-0.0065,-0.0065,-0.0065,-0.0085,-0.0095;
-6,-0.0050,-0.0050,-0.0050,-0.0050,-0.0050,-0.0060,-0.0060;
-5,-0.0035,-0.0035,-0.0035,-0.0035,-0.0035,-0.0040,-0.0050;
-4,-0.0025,-0.0025,-0.0025,-0.0025,-0.0025,-0.0030,-0.0035;
-3,-0.0015,-0.0015,-0.0015,-0.0015,-0.0015,-0.0015,-0.0015;
-2,-0.0010,-0.0010,-0.0010,-0.0010,-0.0010,-0.0010,-0.0010;
-1,-0.0005,-0.0005,-0.0005,-0.0005,-0.0005,-0.0005,-0.0005;
0,0,0,0,0,0,0,0;
1,-0.0005,-0.0005,-0.0005,-0.0005,-0.0005,-0.0005,-0.0005;
2,-0.0010,-0.0010,-0.0010,-0.0010,-0.0010,-0.0010,-0.0010;
3,-0.0015,-0.0015,-0.0015,-0.0015,-0.0015,-0.0015,-0.0015;
4,-0.0025,-0.0025,-0.0025,-0.0025,-0.0025,-0.0030,-0.0035;
5,-0.0035,-0.0035,-0.0035,-0.0035,-0.0035,-0.0040,-0.0050;
6,-0.0050,-0.0050,-0.0050,-0.0050,-0.0050,-0.0060,-0.0060;
7,-0.0065,-0.0065,-0.0065,-0.0065,-0.0065,-0.0085,-0.0095;
8,-0.0080,-0.0080,-0.0080,-0.0080,-0.0080,-0.0100,-0.0120;
9,-0.0100,-0.0100,-0.0100,-0.0100,-0.0100,-0.0125,-0.0145;
10,-0.0120,-0.0120,-0.0120,-0.0120,-0.0120,-0.0150,-0.0154;
11,-0.0140,-0.0140,-0.0140,-0.0140,-0.0140,-0.0175,-0.0200;
12,-0.0160,-0.0160,-0.0160,-0.0160,-0.0160,-0.0200,-0.0225;
13,-0.0180,-0.0180,-0.0180,-0.0180,-0.0180,-0.0240,-0.0255;
14,-0.0205,-0.0205,-0.0205,-0.0205,-0.0205,-0.0255,-0.0280;
15,-0.0225,-0.0225,-0.0225,-0.0225,-0.0225,-0.0280,-0.0310];
CDBatadgtable=-CDBatadggg(2:32,2:8)';
%CDspwdp
CDspwdp121=[0.0045,0.0045,0.0047,0.0047,0.0043,0.0031,0.0000];
CDspwdp91011234=[0.0060,0.0060,0.0064,0.0064,0.0076,0.0082,0.0085];
CDspwdp58=[0.0008,0.0008,0.0010,0.0010,-0.0002,-0.0006,-0.0030];
CDspwdp67=[0.0008,0.0008,0.0020,0.0020,-0.0008,-0.0036,-0.0054];
CDsppad121=[0.00027,0.00027,0.00044,0.00056,0.00064,0.00052,0.00022];
CDsppad58=[0.00004,0.00004,0.00013,0.00021,0.00026,0.00021,0.00015];
CDsppad91011234=[0.00070,0.00070,0.00071,0.00076,0.00115,0.00164,0.00224];
CDsppad67=[-0.00013,-0.00013,-0.00012,-0.00016,-0.00032,-0.00056,-0.00106];
%CDgear
CDgeargg=[-999,0,1,5,10,20,25,30;
-5,0.0275,0.0275,0.0260,0.0255,0.0245,0.0185,0.0155;
-4,0.0275,0.0275,0.0260,0.0255,0.0245,0.0185,0.0150;
-3,0.0275,0.0275,0.0260,0.0255,0.0245,0.0180,0.0140;
-2,0.0275,0.0275,0.0260,0.0255,0.0245,0.0180,0.0130;
-1,0.0280,0.0280,0.0260,0.0255,0.0245,0.0175,0.0125;
0, 0.0280,0.0280,0.0260,0.0255,0.0245,0.0170,0.0115;
1, 0.0280,0.0280,0.0260,0.0255,0.0245,0.0170,0.0110;
2, 0.0280,0.0280,0.0260,0.0255,0.0240,0.0165,0.0100;
3, 0.0280,0.0280,0.0260,0.0255,0.0235,0.0155,0.0095;
4, 0.0280,0.0280,0.0260,0.0255,0.0235,0.0150,0.0090;
5, 0.0280,0.0280,0.0260,0.0255,0.0235,0.0145,0.0085;
6, 0.0280,0.0280,0.0260,0.0245,0.0225,0.0135,0.0080;
7, 0.0280,0.0280,0.0260,0.0245,0.0220,0.0125,0.0075;
8, 0.0280,0.0280,0.0260,0.0240,0.0210,0.0115,0.0075;
9, 0.0275,0.0275,0.0255,0.0235,0.0205,0.0105,0.0070;
10,0.0270,0.0270,0.0250,0.0230,0.0200,0.0100,0.0070;
11,0.0265,0.0265,0.0250,0.0225,0.0190,0.0095,0.0070;
12,0.0260,0.0260,0.0245,0.0220,0.0190,0.0095,0.0070;
13,0.0255,0.0255,0.0240,0.0215,0.0185,0.0095,0.0070;
14,0.0250,0.0250,0.0235,0.0210,0.0180,0.0090,0.0075;
15,0.0245,0.0245,0.0230,0.0205,0.0180,0.0100,0.0080;
16,0.0240,0.0240,0.0225,0.0200,0.0175,0.0105,0.0085;
17,0.0235,0.0235,0.0220,0.0195,0.0175,0.0115,0.0090;
18,0.0235,0.0235,0.0215,0.0190,0.0170,0.0125,0.0100;
19,0.0230,0.0230,0.0210,0.0190,0.0165,0.0130,0.0115;
20,0.0225,0.0225,0.0205,0.0185,0.0160,0.0140,0.0125;
21,0.0220,0.0220,0.0200,0.0180,0.0155,0.0140,0.0125;
22,0.0215,0.0215,0.0195,0.0175,0.0155,0.0130,0.0115;
23,0.0210,0.0210,0.0190,0.0170,0.0150,0.0120,0.0100;
24,0.0205,0.0205,0.0185,0.0170,0.0145,0.0100,0.0085;
25,0.0200,0.0200,0.0180,0.0165,0.0140,0.0080,0.0065];
CDgeartable=CDgeargg(2:32,2:8)';
%CDgroundeffect
Kgatable2=[1.0000,0.9900,0.9700,0.9400,0.8900,0.8300,0.7500,0.6600,0.5700,0.4800,0.3700];
CDgagg=[-999,0,1,5,10,20,25,30;
0, -0.021,-0.017,-0.020,-0.017,-0.027,-0.030,-0.033;
0.5, -0.022,-0.017,-0.020,-0.018,-0.028,-0.031,-0.034;
1.0, -0.023,-0.017,-0.021,-0.018,-0.028,-0.033,-0.035;
1.5, -0.023,-0.017,-0.021,-0.018,-0.028,-0.033,-0.037;
2.0, -0.023,-0.017,-0.022,-0.018,-0.029,-0.035,-0.038;
2.5, -0.023,-0.017,-0.022,-0.019,-0.030,-0.037,-0.040;
3.0, -0.024,-0.017,-0.023,-0.020,-0.031,-0.038,-0.042;
3.5, -0.024,-0.017,-0.023,-0.020,-0.033,-0.040,-0.043;
4.0, -0.025,-0.018,-0.024,-0.021,-0.033,-0.042,-0.046;
4.5, -0.026,-0.018,-0.024,-0.022,-0.035,-0.044,-0.048;
5.0, -0.027,-0.018,-0.025,-0.023,-0.036,-0.046,-0.050;
5.5, -0.027,-0.019,-0.025,-0.023,-0.037,-0.048,-0.052;
6.0, -0.028,-0.020,-0.026,-0.024,-0.038,-0.050,-0.055;
6.5, -0.029,-0.020,-0.027,-0.025,-0.040,-0.053,-0.057;
7.0, -0.029,-0.021,-0.028,-0.025,-0.042,-0.055,-0.060;
7.5, -0.030,-0.021,-0.028,-0.026,-0.043,-0.058,-0.063;
8.0, -0.031,-0.022,-0.029,-0.027,-0.044,-0.060,-0.065;
8.5, -0.032,-0.023,-0.030,-0.028,-0.045,-0.063,-0.068;
9.0, -0.033,-0.023,-0.030,-0.029,-0.047,-0.066,-0.072;
9.5, -0.034,-0.024,-0.031,-0.030,-0.049,-0.069,-0.075;
10.0,-0.035,-0.024,-0.032,-0.031,-0.052,-0.072,-0.078;
10.5,-0.035,-0.025,-0.033,-0.032,-0.053,-0.075,-0.082;
11.0,-0.036,-0.026,-0.034,-0.033,-0.055,-0.078,-0.083;
11.5,-0.037,-0.026,-0.034,-0.034,-0.057,-0.080,-0.086;
12.0,-0.038,-0.027,-0.035,-0.035,-0.058,-0.082,-0.088;
12.5,-0.039,-0.028,-0.036,-0.036,-0.060,-0.084,-0.089;
13.0,-0.040,-0.028,-0.037,-0.037,-0.062,-0.085,-0.090;
13.5,-0.041,-0.029,-0.038,-0.038,-0.063,-0.085,-0.091;
14.0,-0.043,-0.030,-0.038,-0.038,-0.064,-0.086,-0.092];
CDgatable=CDgagg(2:30,2:8)'
%--------------------------------------------------------------------------
%q,alpha,Alphadot,PJedg,H,Vtrue,(M)
dE=-17; %Degree
AlphaDG=12; %Degree
BataDG=0; %Degree
Flaps=10; %Degree
Alphadot=0; %rad/s
q=0 %rad/s
ALpha=AlphaDG/57.2958; %Rad
Bata=AlphaDG/57.2958; %Rad
Vtrue=80; %m/s
H=0 %m
Ksp=0 %扰流板开1,关0
%cbar=27.37; %ft
%求解中间量
%ALphaDG=Alpha*57.2958;
M=0.3
Alphadotdl=Alphadot*27.37/(2*Vtrue*3.28); %(rad/s*ft)/(ft/s)
qdl=q*27.37/(2*Vtrue*3.28); %(rad/s*ft)/(ft/s)
%--------------------------------------------------------------------------
%求值
CLBasic=interp2(Alphadgtable,Flapstable,CLBasictable,(AlphaDG),Flaps,'linear');
CLAlphadot=interp1(Mtable,CLAlphadottable,M,'linner')*Alphadotdl;
CLq=interp2(Mtable,Hinchtable,CLqtable,M,H*3.28,'linner')*qdl;
CLdEi=interp2(Mtable,Hinchtable2,CLdEitable,M,H*3.28,'linner')*interp2(Alphadgtable,Flapstable,Kalphatable,(AlphaDG),Flaps,'linear')*dE;
CLdEo=interp2(Mtable,Hinchtable2,CLdEotable,M,H*3.28,'linner')*interp2(Alphadgtable,Flapstable,Kalphatable,(AlphaDG),Flaps,'linear')*dE;
Clsp45=Ksp*(interp2(Alphadgtable,Flapstable,Clsp58table,(AlphaDG),Flaps,'linear')+interp2(Alphadgtable,Flapstable,Clsp58table,(AlphaDG),Flaps,'linear')*interp1(Flapstable,Ksp20,Flaps,'linear')*2+interp2(Alphadgtable,Flapstable,Clsp67table,(AlphaDG),Flaps,'linear')*2*interp1(Flapstable,Ksp20,Flaps,'linear')+interp2(Alphadgtable,Flapstable,Clsp91011234table,(AlphaDG),Flaps,'linear')*5.4)*(1+interp2(Alphawdptable,Flapstable,KGEBtable,(AlphaDG),Flaps));
CLgear=interp2(Alphadgtable,Flapstable,CLgeartable,(AlphaDG),Flaps,'linear')*interp1(Mtable,KCLgearM,M,'linear');
if AlphaDG<=14
CLga=interp1(Hgatable,Kgatable,H*3.28,'linear')*interp2(Alphawdptable,Flapstable,CLgatable,(AlphaDG),Flaps,'linear');
elseif AlphaDG>14
CLga=interp1(Hgatable,Kgatable,H*3.28,'linear')*interp2(Alphawdptable,Flapstable,CLgatable,14,Flaps,'linear');
end
CL=CLBasic+CLAlphadot+CLq+CLdEi+CLdEo+Clsp45+CLgear+CLga
CMBasic=interp2(Alphadgtable,Flapstable,CMBasictable,(AlphaDG),Flaps,'linear');
CMAlphadot=interp2(Mtable,Hinchtable,CMAlphadottable,M,H*3.28,'linner')*Alphadotdl;
CMq=interp2(Mtable,Hinchtable,CMqtable,M,H*3.28,'linner')*qdl;
CMdEi=interp2(Mtable,Hinchtable2,CMdEitable,M,H*3.28,'linner')*interp2(Alphadgtable,Flapstable,Kalphatable,(AlphaDG),Flaps,'linear')*dE;
CMdEo=interp2(Mtable,Hinchtable2,CMdEotable,M,H*3.28,'linner')*interp2(Alphadgtable,Flapstable,Kalphatable,(AlphaDG),Flaps,'linear')*dE;
CMgear=interp2(Alphadgtable,Flapstable,CMgeartable,(AlphaDG),Flaps,'linear')*interp1(Mtable,KCMgearM,M,'linear');
if AlphaDG<=14
CMga=interp1(Hgatable,Kgatable,H*3028,'linear')*interp2(Alphawdptable,Flapstable,CMgatable,(AlphaDG),Flaps,'linear');
elseif AlphaDG>14
CMga=interp1(Hgatable,Kgatable,H*3.28,'linear')*interp2(Alphawdptable,Flapstable,CMgatable,14,Flaps,'linear');
end
CM=CMBasic+CMAlphadot+CMq+CMdEi+CMdEo+CMgear+CMga
CDBasic=interp2(Alphadgtable,Flapstable,CDBasictable,(AlphaDG),Flaps,'linear')
CDBata=interp2(Batadgtable,Flapstable,CDBatadgtable,BataDG,Flaps,'linear')
CDspwdp=Ksp*(interp1(Flapstable,CDspwdp121,Flaps,'linear')+interp1(Flapstable,CDspwdp91011234,Flaps,'linear')*5.3+interp1(Flapstable,CDspwdp58,Flaps,'linear')*2+interp1(Flapstable,CDspwdp67,Flaps,'linear')*2+[interp1(Flapstable,CDsppad121,Flaps,'linear')+interp1(Flapstable,CDsppad91011234,Flaps,'linear')*5.3+interp1(Flapstable,CDsppad67,Flaps,'linear')*2+interp1(Flapstable,CDsppad58,Flaps,'linear')*2]*(AlphaDG-4));
CDgear=interp2(Alphadgtable,Flapstable,CDgeartable,(AlphaDG),Flaps,'linear')
if AlphaDG<=14
CDga=interp1(Hgatable,Kgatable2,H*3028,'linear')*interp2(Alphawdptable,Flapstable,CMgatable,(AlphaDG),Flaps,'linear')
elseif AlphaDG>14
CDga=interp1(Hgatable,Kgatable2,H*3.28,'linear')*interp2(Alphawdptable,Flapstable,CMgatable,14,Flaps,'linear')
end
CD=CDBasic+CDBata+CDspwdp+CDgear+CDga |