xhxn.m

来自「飞机续航性能计算」· M 代码 · 共 133 行

% No.02
% 飞行性能计算－续航性能
% 确定飞机的续航性能，包括航程、航时
clc;
clear;
close;


G = 8150.0*9.81;
S = 23.0;

TCx0 = [
    0.013478
    0.013493
    0.013507
    0.013524
    0.014611
    0.025896
    0.034436
    0.034441
    0.031589
    0.030538
    0.029696
    0.029255
    0.028396
    0.027771
    0.027273
    0.026845];
TAi = [
    0.206530
    0.206289
    0.206047
    0.205806
    0.205560
    0.209018
    0.223322
    0.244748
    0.267197
    0.290234
    0.313367
    0.336425
    0.359560
    0.382798
    0.406161
    0.429703];

TQh_max = 1.0e+003 * [
    3.0445       NaN       NaN       NaN       NaN       NaN       NaN       NaN
    3.1757    2.9130    2.4096       NaN       NaN       NaN       NaN       NaN
    3.3423    3.0610    2.5006    2.0944       NaN       NaN       NaN       NaN
    3.5478    3.2274    2.6504    2.1879       NaN    1.1817    1.0181       NaN
    3.7629    3.4310    2.7975    2.2867    1.5275    1.2338    1.0652       NaN
    3.9986    3.6676    2.9929    2.4398    1.6443    1.3171    1.1179       NaN
    4.2705    3.9022    3.2076    2.5965    1.8255    1.3967    1.2259    0.8178
    4.5420    4.1495    3.4555    2.7899    1.9746    1.5253    1.3261    0.8840
    4.8058    4.4469    3.7097    3.0245    2.1790    1.6739    1.4514    0.9877
    5.0670    4.7402    3.9768    3.2722    2.3658    1.8408    1.5975    1.0916
    5.3712    5.0628    4.2602    3.5131    2.5606    2.0183    1.7658    1.2001
       NaN       NaN    4.5773    3.7777    2.7942    2.2331    1.9516    1.3147
       NaN       NaN       NaN    4.0691    3.0417    2.4477    2.1611    1.4521
       NaN       NaN       NaN    4.3942    3.2763    2.6784    2.3727    1.5895
       NaN       NaN       NaN       NaN    3.5489    2.9094    2.6072    1.7499
       NaN       NaN       NaN       NaN    3.8330    3.1439    2.8417    1.9274
       NaN       NaN       NaN       NaN       NaN    3.3784    3.0857    2.0866
       NaN       NaN       NaN       NaN       NaN    3.6371    3.2934    2.2323
       NaN       NaN       NaN       NaN       NaN       NaN    3.4695    2.3439
       NaN       NaN       NaN       NaN       NaN       NaN    3.5987    2.4481]';
TH = [0, 1, 3, 5, 8, 10, 11];

TQh_cru = TQh_max*0.9;

Pxy = [];
AA = [];
M = 0.1:0.1:2.0;

TVy = [
   64.9072
   59.5824
   50.1006
   40.3454
   25.2266
   13.7551
    8.4080];

tt = [];
LL = [];
for i = 1:7
    H = TH(i)*1000;
    [T,P,D,A] = GetAir(H);
	
	AA = [AA, A];
    V = M * A;
    Vy = TVy(i);
	
	Cy = 2.0*G./(D*S*V.^2);
	
	Cx0 = interp1(0.5:0.1:2.0, TCx0, M, 'cube');
	Ai = interp1(0.5:0.1:2.0, TAi, M, 'cube');
	
	Cx = Cx0 + Ai.*Cy.^2;
	
	K = Cy./Cx;
	
	% 求有利飞行M数
    [Kmax, k] = max(K);
    Mo = k*0.1;
    Vo = Mo*A*3.6;
    
    Qh_cru = TQh_cru(i, :);
    Qh = interp1(0.1:0.1:2.0, Qh_cru, Mo, 'cube')
    
    % 上升到巡航高度所需时间
    th =H/Vy/3600.0
    
    % 取空中飞行可用燃油量为1900.0kg，假设上升段和下滑段耗时、前进距离相同
    t = (1900.0-th*Qh/0.9)./Qh;
    tt = [tt ,t];
    L = Vo*t;
    LL = [LL, L];
end

plot(TH, tt);
xlabel('飞行高度/km');
ylabel('航时/h');
title('不同飞行高度上的续航时间');
grid on;

figure;
plot(TH, LL);
xlabel('飞行高度/km');
ylabel('航程/km');
title('不同飞行高度上的航程');
grid on;