blob: 79be1d056df50572990d6e44ef63e62999d4e704 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
|
function F = Trim_Cessna_Body(x)
m = 1043.26;
S_ref = 16.1651;
d = 1.225;
V = 60;
g = 9.8;
C_L_0 = 0.25;
C_L_alpha = 4.47;
C_L_delta_e = 0.3476;
C_D_0 = 0.036;
C_D_alpha = 0.13;
C_D_delta_e = 0.06;
C_M_0 = -0.02;
C_M_alpha = -1.8;
C_M_delta_e = -1.28;
C_L = C_L_0 + C_L_alpha*x(1) + C_L_delta_e*x(2);
% C_D = C_D_0 + C_D_alpha*abs(x(1)) + C_D_delta_e*abs(x(2));
C_D = C_D_0 + 0.3*C_L^2;
C_M = C_M_0 + C_M_alpha*x(1) + C_M_delta_e*x(2);
F = [C_M;
0.5*d*V*V*S_ref*(-C_D*cos(x(1)) + C_L*sin(x(1))) - m*g*sin(x(1)) + x(3);
0.5*d*V*V*S_ref*(-C_D*sin(x(1)) - C_L*cos(x(1))) + m*g*cos(x(1))];
|
