性能.m

用MATLAB编写关于飞机设计的matlab程序

clc;
clear;
%提示输入控制电压向量和前飞速度向量
%调用Pavailable函数，创建某一电压下的可用功率曲线，并在同一张图上绘出曲线
fprintf('********************  最小输入电压不得低于2.4伏  *************************\n');
fprintf('********************  最大输入电压不得高于9.6伏  *************************\n');
Umin=input('********************  请输入最小电压 （1.2倍数） *************************\n');
Umax=input('********************  请输入最大电压 （1.2倍数） *************************\n');
%********************  创建控制电压向量  ********************
U=Umin:1.2:Umax;
I=(Umax-Umin)/1.2+1;
%****************************************
Vmin=input('********************  请输入最小飞行速度（米/秒）  ********************\n');
Vmax=input('********************  请输入最大飞行速度（米/秒）  ********************\n');
%********************  创建飞行速度向量  ********************
V=Vmin:0.1:Vmax;
J=(Vmax-Vmin)/0.1+1;
%****************************************
Hf1=figure(3);
Ha1=axes('Parent',Hf1);
Hf2=figure(2);
Ha2=axes('Parent',Hf2);

set(Ha1,'nextplot','add');
set(Ha2,'nextplot','add');
for i=1:I
    Uin=U(i);
    for j=1:J
        Vf=V(j);
        [P,T]=Pavailable3(Uin,Vf);
        Pky(j)=P;
        Tky(j)=T;
    end
    line('xdata',V,'Ydata',Pky,'Parent',Ha1,'color','blue');
    line('xdata',V,'Ydata',Tky,'Parent',Ha2,'color','blue');
end

%全机极曲线及平飞需用推力曲线和平飞需用功率曲线
AR_w=8;Wto=2.4;g=9.81;den=1.225;
S_w=0.575;S_fus=0.071*0.075;S_tl=0.0812+0.0235;S_bld=0.0094;S_lg=0.0029;S_wh=0.0013*3;
Cd_w=0.0091;Cd_fus=0.29;Cd_tl=0.021;Cd_bld=0.7;Cd_lg=1.10;Cd_wh=1.460;
Cd_w=1.1*(Cd_w+Cd_fus*S_fus/S_w+Cd_tl*S_tl/S_w+Cd_bld*S_bld/S_w+Cd_lg*S_lg/S_w+Cd_wh*S_wh/S_w);
Cl=2*Wto*g./(den*V.^2*S_w);
Cd_t=Cd_w+Cl.^2/(pi*AR_w);
figure(1)
plot(Cd_t,Cl);
xlabel('Cd');
ylabel('Cl');
title('全机极曲线');
Kmax=max(Cl./Cd_t);
fprintf('全机最大升阻比为 %.3f \n',Kmax);
K=Cl./Cd_t;
Tpx=Wto*g./K;
f=polyfit(V,Tpx,2);
t=polyval(f,V);
figure(2)
plot(V,Tpx);
xlabel('V(m/s)');
ylabel('Tpx(N)               Tky(N)');
title('平飞可用－需用推力曲线');
Ppx=Tpx.*V;
f=polyfit(V,Ppx,2);
t=polyval(f,V);
figure(3)
plot(V,Ppx);
xlabel('V(m/s)');
ylabel('Ppx(w)                Pky(w)');
title('平飞可用－需用功率曲线');
%起飞滑跑距离
Clmax=1.13;WL=4.17;f=0.035;
P=1*g;
Tz=10*P/(Wto*g);
Ppj=0.95*Tz;
Lqh=0.908*g*WL/(Clmax*(Ppj-f));
fprintf('起飞滑跑距离为 %.3f m\n',Lqh);