clc;
清晰;
关闭所有;
%定义变量
V = 202;%的速度(米/秒)
R = 1000;%曲率半径(米)
g = 9.8;%重力加速度(m / s ^ 2)
ρ= 1.058;%空气密度(公斤/米^ 3)
S = 24;%机翼面积(m ^ 2)
m = 7800;%飞机重量(公斤)
CDzero = 0.02;%阻力系数(无量纲)
K = 0.05;%拖动调整因素(无量纲)
Iy = 57000;%惯性矩(公斤* m ^ 2)
厘米= 0.1;%力矩系数(无量纲)
Sc = 62;
α= (0 5 10 15 20);%入射角度(度)
γ= [0 10 30 60 90 120 150 180 210 240 270 300 330 360];%翼偏转角度(度)
%计算值CL, CD, T为每个伽马和α的价值
为我= 2:长度(γ)
为j = 1:长度(α)
%计算CL
CL (i) = 2 * (m * g * cosd(γ(i)) + (m * V ^ 2 / (R * g)) * g) /(ρ* V ^ 2 * S);
CL_vector (i) = CL (i);
%计算CD
CD = CDzero + K * CL(我)^ 2;
CD_vector (i) = CD;
%计算T
T = m * g *信德(γ(i)) +(ρ* V ^ 2 * * CD) / 2;
T_vector (i) = T;
%计算方程的最终值
diff_alpha_gamma =(α(j) -γ(我))* 180(π/;%将γ的区别和α转换为弧度
eqn = ((Iy * (V ^ 2 / R) * (diff_alpha_gamma) ^ 2) /(γ(i)γ(张))*π/ 180)-(0.5 *ρ* V ^ 2 * Sc *厘米)
托尔= 1 e-6;
如果(eqn) <托尔%使用相对宽容
流(“对γ= %。1f and alpha = %d, the equation is satisfied\n"γ(i),α(j));
其他的
流(γ= %的方程是不满意。1 f和α= % d \ n 'γ(i),α(j));
结束
结束
结束
%计算俯仰力矩(电梯偏转)
CL_alpha_alpha = 2 *π;%的导数CL入射角,无量纲
CL_zero = 0;%发病率升力系数为零,无量纲
CL_deltae = 0.3;%的导数CL电梯偏转角,无量纲
deltae = 1 / CL_deltae * (CL-CL_alpha_alpha-CL_zero)
流(“电梯时刻γ= %。1 f和α= % d埃斯特%。4 f m \ n”,αγdeltae);
图
情节(γ,deltae)
几个点,关于方程和代码结构
- 似乎也写评估不正确,使用“-”来重新排列方程和RHS assigni值变量eqn在韩
eqn = ((Iy * (V ^ 2 / R) * (diff_alpha_gamma) ^ 2) /(γ(i)γ(张))*π/ 180)-(0.5 *ρ* V ^ 2 * Sc *厘米)
%的差异使用“-”
- 使用一个相对宽容eqn <托尔代替abs宽容,因为α可以大范围的值,CL, CD,和CM发生显著的变化。
- 在for循环中,您可以使用因为我= 2:长度(γ)作为δ射线方程涉及术语。
- 最后,您需要用((γ(i)γ(张))* 180(π/术语π/ 180,因为所有角变量需要在相同的单位。在一个地方你在度和其他地方的弧度。
