pfxn.m

来自「飞机平飞性能计算」· M 代码 · 共 143 行

% No.02
% 飞行性能计算－平飞性能
% 确定飞机的平飞性能，包括平飞最大速度、平飞最小速度
clc;
clear;
close;

G = 8150.0;
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];

TPky = 1.0e+003 * [
%     3.0018       NaN       NaN       NaN       NaN       NaN       NaN       NaN
%     2.9834    2.6958    2.2556    1.8459       NaN       NaN       NaN       NaN
%     3.0028    2.7438    2.2770    1.8663       NaN       NaN       NaN       NaN
%     3.0734    2.8048    2.3136    1.8978       NaN    1.0376    0.9077       NaN
%     3.1595    2.8952    2.3781    1.9582    1.4069    1.0496    0.9085    0.5996
%     3.2652    2.9864    2.4786    2.0206    1.4376    1.0818    0.9394    0.6302
%     3.3858    3.0974    2.5804    2.1168    1.4891    1.1367    0.9873    0.6710
%     3.5073    3.2184    2.7050    2.2374    1.5753    1.2122    1.0714    0.7229
%     3.6507    3.3497    2.8463    2.3672    1.7220    1.3176    1.1485    0.7993
%     3.7713    3.4930    3.0275    2.5310    1.8560    1.4487    1.2692    0.8664
%     3.8956    3.6268    3.2375    2.7315    2.0466    1.6191    1.4003    0.9616
%        NaN    3.7866    3.4661    2.9697    2.2376    1.7794    1.5532    1.0714
%        NaN       NaN    3.7057    3.1722    2.4185    1.9324    1.7333    1.1672
%        NaN       NaN       NaN    3.3063    2.5722    2.0853    1.8772    1.2629
%        NaN       NaN       NaN    3.3503    2.7230    2.2533    2.0273    1.3632
%        NaN       NaN       NaN       NaN    2.8362    2.3590    2.1399    1.4443
%        NaN       NaN       NaN       NaN    2.9239    2.4692    2.2545    1.5247
%        NaN       NaN       NaN       NaN    2.9707    2.5564    2.3514    1.5657
%        NaN       NaN       NaN       NaN       NaN    2.6373    2.4285    1.6041
%        NaN       NaN       NaN       NaN       NaN    2.6899    2.4714    1.6058]';
    3.0018       NaN       NaN       NaN       NaN       NaN
    2.9834    2.2556    1.8459       NaN       NaN       NaN
    3.0028    2.2770    1.8663       NaN       NaN       NaN
    3.0734    2.3136    1.8978       NaN    1.0376       NaN
    3.1595    2.3781    1.9582    1.4069    1.0496    0.5996
    3.2652    2.4786    2.0206    1.4376    1.0818    0.6302
    3.3858    2.5804    2.1168    1.4891    1.1367    0.6710
    3.5073    2.7050    2.2374    1.5753    1.2122    0.7229
    3.6507    2.8463    2.3672    1.7220    1.3176    0.7993
    3.7713    3.0275    2.5310    1.8560    1.4487    0.8664
    3.8956    3.2375    2.7315    2.0466    1.6191    0.9616
       NaN    3.4661    2.9697    2.2376    1.7794    1.0714
       NaN    3.7057    3.1722    2.4185    1.9324    1.1672
       NaN       NaN    3.3063    2.5722    2.0853    1.2629
       NaN       NaN    3.3503    2.7230    2.2533    1.3632
       NaN       NaN       NaN    2.8362    2.3590    1.4443
       NaN       NaN       NaN    2.9239    2.4692    1.5247
       NaN       NaN       NaN    2.9707    2.5564    1.5657
       NaN       NaN       NaN       NaN    2.6373    1.6041
       NaN       NaN       NaN       NaN    2.6899    1.6058]';
TH = [0, 3, 5, 8, 10, 13.5];

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

% 求有利飞行M数
disp('不同高度上的有利M数和有利速度');
for i = 1:6
    
    H = TH(i)*1000
    [T,P,D,A] = GetAir(H);
	
	AA = [AA, A];
    V = M * A;
	
	Cy = 2.0*9.8*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;
	
	[Kmax, k] = max(K);
    Mo = k*0.1
    Vo = Mo*A
    P = G./K;
    Pxy = [Pxy; P];
end

%%彩色图线%%
% plot(   M, Pxy(1,:), 'r', M, TPky(1,:), 'r',...
% 		M, Pxy(2,:), '-+g', M, TPky(2,:), '-+g',...
% 		M, Pxy(3,:), '-ob', M, TPky(3,:), '-ob',...
% 		M, Pxy(4,:), '-*m', M, TPky(4,:), '-*m',...
% 		M, Pxy(5,:), '-xc', M, TPky(5,:), '-xc',...
% 		M, Pxy(6,:), '-sk', M, TPky(6,:), '-sk',...
% 		M, Pxy(7,:), '-dy', M, TPky(7,:), '-dy',...
% 		M, Pxy(8,:), '-^r', M, TPky(8,:), '-^r');

%%黑白图线%%
plot(   M, Pxy(1,:), '-k', M, TPky(1,:), '-k',...
		M, Pxy(2,:), '-+k', M, TPky(2,:), '-+k',...
		M, Pxy(3,:), '-ok', M, TPky(3,:), '-ok',...
		M, Pxy(4,:), '-*k', M, TPky(4,:), '-*k',...
		M, Pxy(5,:), '-xk', M, TPky(5,:), '-xk',...
		M, Pxy(6,:), '-sk', M, TPky(6,:), '-sk',...
        'MarkerSize', 5);
    
axis([0.0, 2.0, 0.0, 5000.0]);
xlabel('Ma');
ylabel('P/kgf');
%title('根据此图求平飞最大M数和平飞最小M数');
grid on;