%% 汽车动力性计算(自己编的动力性计算程序，供大家计算动力性时参考，具体参数大家

根据所给程序对应输入，并对坐标轴数值按需要进行修改)

clc;

clear;

close all;

%%根据所给发动机数据拟合外特性曲线 (发动机数据按照你所得到的数据进行输入)

n_test=[1000 1200 1400 1800 2200 2400 2600 2800 3000 3200 3400 3600 3800 4000 4200

4400 4600 5000 ];

T_test=[201 206 216 2 18 221 221 227 228 235 246 228 227 223 215 210 208 205 109];

figure(1)

plot(n_test,T_test, 'g');

hold on

grid on

p=polyfit(n_test,T_test,7);

n=[1000:1:5000];

Ttq=polyval(p,n);

plot(n,Ttq, 'k');

xlabel('发动机转速n(r/min)');

ylabel('发动机转矩Ttq(N*m)');

title('发动机转矩曲线');

legend('测试曲线','拟合曲线');

axis([1000 5000 100 300]);

%%所给车型动力总成相关参数

ig=[4.148 2.37 1.556 1.115 0.859 0.686];

i0=3.394;

eta=0.9;

r=0.307;

M=1700;

g=9.8;

c=1.2;

f0=0.009;

f1=0.00018;

f4=0.0003;

CD=0.3;

A=2.54;

Iw=3.0;

If=0.04;

%% 发动机外特性曲线图

figure(2)

hold on

grid on

for i=length(n);

Pe=Ttq.*n/9550;

end

[AX,H1,H2]=plotyy(n,Ttq,n,Pe);

xlabel('发动机转速n(r/min)');

ylabel('发动机转矩Ttq(N*m)');

ylabel(AX(2),'发动机功率Pe(Kw)');

title('发动机外特性曲线');

%% 各挡位速度曲线

%计算各挡位车速

for i=1:length(ig);

ua(i,:)=0.377*r*n/ig(i)/i0;

end

%计算各档位最高车速

uamax=ua(:,length(ua(1,:)));

figure(3)

hold on

for i=1:length(ig);

plot(n,ua(i,:),'k');

end

hold on

grid on

xlabel('转速n(r/min)');

ylabel('各挡位车速(km/h)');

title('各挡位车速-转速表');

legend('1挡车速','2挡车速','3挡车速','4挡车速','5挡车速');

axis([1000 5000 0 270]);

%% 驱动力和行驶阻力平衡图

%计算滚动阻力系数

for i=1:length(ig);

f(i,:)=f0+f1*(ua(i,:)/100)+f4*(ua(i,:)/100).^4;

end

%计算滚动阻力

for i=1:length(ig);

Ff(i,:)=c*M*g*f(i,:);

end

%计算空气阻力

for i=1:length(ig);

Fw(i,:)=CD*A*(ua(i,:).^2)/21.15;

end

%计算行驶阻力

for i=1:length(ig);

F(i,:)=Ff(i,:)+Fw(i,:);

end

%计算汽车驱动力

for i=1:length(ig);

Ft(i,:)=Ttq*ig(i)*i0*eta/r;

end

figure(4)

hold on

for i=1:length(ig);

plot(ua(i,:), Ft(i,:), 'k');

plot(ua(i,:), F(i,:), 'r');

plot(ua(i,:), Ff(i,:), 'b ');

end

hold on

grid on

xlabel('车速(km/h)');

ylabel('驱动力、行驶阻力(N)');

legend('驱动力Ft','行驶阻力Ff+Fw','滚动阻力Ff');

title('驱动力-行驶阻力