汽车理论课程设计,百公里燃油消耗率计算程序

m=3880;G=38024;r=0.367;n=0.85;f=0.013;CdA=2.77;If=0.218;
Iw1=1.798;Iw2=3.598;ig=1;pg=7;Qi=0.299;miles=1075;i0=6.50;
syms nrnr
PPee=(-3.8445/1000^4)*nrnr^5+(40.874/1000^3)*nrnr^4-(165.44/1000^2)*nrnr^3+(295.27/1000)*nrnr^2-19.313*nrnr;%%依据外特性计算的发动机能发出的最大功率
glds=diff(PPee,nrnr,1);%%功率对转速的导函数
Nmax=solve(glds);%%功率对转速的导函数等于0时的转速
Nmax=double(Nmax);
Nmax=Nmax(2)%%取出发动机功率最大值时对应的转速
UA=[25,40,50];%速度矩阵，可变
S=[50,250,250];%距离矩阵 ，可变
N=[815,1207,1614,2012,2603,3006,3403,3804];%转速矩阵
B=[1326.8 -416.46 72.379 -5.8629 0.17768;
1354.7 -303.98 36.657 -2.0553 0.043072;
1284.4 -189.75 14.524 -0.51184 0.0068164;
1122.9 -121.59 7.0035 -0.18517 0.0018555;
1141.0 -98.893 4.4763 -0.091077 0.00068906;
1051.2 -73.714 2.8593 -0.05138 0.00035032;
1233.9 -84.478 2.9788 -0.047449 0.00028230;
1129.7 -45.291 0.71113 -0.00075215 -0.000038568];%拟合系数矩阵
I=[5.0 5.17 5.43 5.83 6.17 6.33 6.50];%主减速器传动比矩阵,可变
Qq=[];%空矩阵，存放各个主减速比时的百公里油耗，即最终结果

for jsxx=1:7 %计数变量，因为有7个主减速比，可变
Qs=0;
smqa=0;
Qd=0;
i0=I(jsxx);
%%%匀速行驶阶段燃油消耗计算
x=1:3; %匀速行驶阶段数目，可变
for k=x
ua=UA(1,k);
Pe=(G*f*ua/3600+CdA*ua*ua*ua/76140)/n;%匀速行驶阶段阻力功率
PE=[1;Pe;Pe*Pe;Pe*Pe*Pe;Pe*Pe*Pe*Pe];
BE=B*PE;%各给定转速下燃油消耗率矩阵
nr=ig*i0*ua/(0.377*r);
Tq1=-19.313+295.27*nr/1000-165.44*nr^2/10000+40.874*nr^3/1000^3-3.8445*nr^4/1000^4;%%%
Pmax1=9550*nr;%%%
if nr<=Nmax%最大功率对应的转速，由功率拟合曲线方程解出，可变%%%
if Pmax1ig=1.71;%三挡传动比，可变
nr=ig*i0*ua/(0.377*r);
else
end
else
end
j=1;
be=0;
while j<=8 && be==0
x1=nr-N(1,j);
x2=nr-N(1,j+1);
if x1*x2<0
x3=BE(j+1,1)-BE(j,1);
x4=nr-N(1,j);
x5=N(1,j+1)-N(1,j);
x6=BE(j,1);
be=x3*x4/x5+x6;
else j=j+1;
end
end
Qs=Pe*be*S(1,k)/(102*ua*pg)+Qs;
end
%%%匀速行驶阶段燃油消耗计算

%%%加速阶段燃油消耗率计算
UUA=[25 40];%各加速阶段初始速度矩阵
UB=[40 50 25];%各加速阶段最终速度矩阵及最后的速度
A=[0.25 0.2 0.36];%各加速阶段加速度矩阵
for mm=1:2 %mm为计数变量，根据加速阶段选取
ua=UUA(mm);ub=UB(mm);a=A(mm);
t=1;s=0;
deltaT=0.5/(3.6*a);
delta=1+(Iw1+Iw2)/(m*r*r)+If*ig*ig*i0*i0*n/(m*r*r);%系数δ
while ua<=ub
Pe=(G*f*ua/3600+CdA*ua*ua*ua/76140+delta*m*ua*a/3600)/n;%加速阶段阻力功率
PE=[1;Pe;Pe*Pe;Pe*Pe*Pe;Pe*Pe*Pe*Pe];
BE=B*PE;%各给定转速下燃油消耗率矩阵
nr=ig*i0*ua/(0.377*r);%转速
Tq2=-19.313+295.27*nr/1000-165.44*nr^2/10000+40.874*nr^3/1000^3-3.8445*nr^4/1000^4;%%%
Pmax2=9550*nr;%%%%%%%
if nr<=Nmax%最大功率对应的转速，由功率拟合曲线方程解出，可变%%%
if

Pmax2ig=1.71;%三挡传动比，可变
nr=ig*i0*ua/(0.377*r);
else
end
else
end
j=1;
be=0;
while j<=8 && be==0
x1=nr-N(1,j);
x2=nr-N(1,j+1);
if x1*x2<0 %判断转速介于哪两个转速之间
x3=BE(j+1,1)-BE(j,1);
x4=nr-N(1,j);
x5=N(1,j+1)-N(1,j);
x6=BE(j,1);
be=x3*x4/x5+x6;
else j=j+1;
end
end
Q(t)=Pe*be/(367.1*pg);%各个速度点对应的燃油消耗
t=t+1;
ua=ua+0.5;
end
for l=1:(t-1)
s=s+Q(l)*deltaT;
end
Qa=s-0.5*[Q(1)+Q(t-1)]*deltaT;
smqa=smqa+Qa;%加速阶段总的燃油消耗
end
%%%加速阶段燃油消耗率计算


%怠速阶段燃油消耗
ts=(UB(mm)-UB(mm+1))/(3.6*A(mm+1));
Qd=Qi*ts;%怠速阶段油耗
Qq(jsxx)=(smqa+Qs+Qd)*100/miles
end