clc % 清屏
clear all; % 删除workplace变量
close all; % 关掉显示图形窗口
format short
warning off
H = 11; %飞行高度
Ma = 0.8; %飞行马赫数
r = 1.4; % 纯空气 气体绝热指数
% 输入温度给定值
Tind_f=288.15;
Tind_CDFS=428.56862609;
Tind_hcom=473.603961;
% 输入压力给定值
Pind_f=1.01325;
Pind_CDFS=3.5464;
Pind_hcom=4.886;
load(‘D:\搜狗高速下载\MATLAB数学建模\Matlab数学建模经典案例实战源程序\第25章\问题1\程序\fan.mat’); % 加载风扇特性数据——原始数据
% 输入压气机管道总温、总压
tube_TP = tube(H,Ma,r); % function tube_TP=tube(H,Ma,r)
Tin_f = tube_TP(1,1); %风扇输入总温
Pin_f = tube_TP(1,2); %风扇输入总压
%% 风扇
alpha = 0; % 导叶角度
nspeed =0.95; % 风扇和CDFS的物理转速都为0.95
ndata=size(fan);
prc_f=fan(:,1:3:ndata(1,2)-2); % 增压比
prc_f(1,:)=[];
Wc_f=fan(:,2:3:ndata(1,2)-1); % 流量
Wc_f(1,:)=[];
nf=size(prc_f);
figure(1),
for i=1:nf(1,2)
zz(:,i) = (prc_f(:,i) - min(prc_f(:,i)))./(max((prc_f(:,i)))-min(prc_f(:,i)));
plot(zz(:,i),Wc_f(:,i),‘linewidth’,3);
hold on
grid on
end
xlabel(‘压比函数’);ylabel(‘流量’);
grid on;title(‘风扇特性中流量随压比函数值变化图’); axis tight
%%
figure(2),
for i=1:nf(1,2)
Wcf(1,i)=Wc_f(nf(1,1),i);
prcf(1,i)=prc_f(nf(1,1),i);
plot(Wc_f(:,i),prc_f(:,i),‘linewidth’,3);
hold on
grid on
end
plot(Wcf,prcf,‘linewidth’,3)
xlabel(‘换算流量’);ylabel(‘增压比’);
grid on;title(‘典型的压气机流量特性曲线’); axis tight