同济 汽车电子 汽车理论 大作业 动力性能与经济性能计算 Matlab

汽车理论大作业

同济大学汽车学院

xxx xxxxxx 2013.12

%%0. 作业要求

%%1. 发动机外特性曲线

nmin=600; nmax=4000;

n=nmin:0.01:nmax; m=n/1000;

Tq=-19.313+295.27*m-165.44*m.^2+40.874*m.^3-3.8445*m.^4; figure; plot(n,Tq);

title('发动机外特性曲线 Tq-n'); xlabel('n/(r/min)'); ylabel('Tq/Nm');

Pe=(Tq.*n)/9550; figure; plot(n,Pe);

title('发动机外特性曲线 Pe-n'); xlabel('n/(r/min)'); ylabel('Pe/kw');

%%2. 驱动力曲线

nmin=600; nmax=4000;

n=nmin:0.01:nmax; m=n/1000;

Tq=-19.313+295.27*m-165.44*m.^2+40.874*m.^3-3.8445*m.^4; r=0.367; elta=0.85; i0=5.83;

ig_4=[6.09 3.09 1.71 1.00];

ig_5=[5.56 2.769 1.644 1.00 0.793]; ua_4=zeros(4,length(n)); Ft_4=zeros(4,length(n)); ua_5=zeros(5,length(n)); Ft_5=zeros(5,length(n)); for i=1:4

ua_4(i,:)=(0.377*n.*r)/(i0*ig_4(i)); Ft_4(i,:)=Tq*ig_4(i)*i0*elta/(r*1000); end

for j=1:5

ua_5(j,:)=(0.377*n.*r)/(i0*ig_5(j)); Ft_5(j,:)=Tq*ig_5(j)*i0*elta/(r*1000); end

%%figure;

plot(ua_4(1,:),Ft_4(1,:),ua_4(2,:),Ft_4(2,:),ua_4(3,:),Ft_4(3,:),ua_4(4,:),Ft_4(4,:));

title('4档变速箱驱动力曲线 Ft-ua'); xlabel('ua/(km/h)'); ylabel('Ft/kN'); text(15.5,14,'Ft1'); text(22,7.7,'Ft2'); text(35,4.5,'Ft3'); text(75,2.7,'Ft4');

%%figure;

plot(ua_5(1,:),Ft_5(1,:),ua_5(2,:),Ft_5(2,:),ua_5(3,:),Ft_5(3,:),ua_5(4,:),Ft_5(4,:),ua_5(5,:),Ft_5(5,:)); title('5档变速箱驱动力曲线 Ft-ua'); xlabel('ua/(km/h)'); ylabel('Ft/kN'); text(16,13,'Ft1'); text(20,7,'Ft2'); text(35,4.5,'Ft3'); text(75,2.8,'Ft4');

text(105,2.2,'Ft5');

%%3. 计算空载和满载时最高车速

m0=1800; m1=3880; f=0.013; Cd=0.68; A=4.07; g=9.8;

Ff_0=m0*g*f; Ff_1=m1*g*f;

ua4m=0.377*r*nmax/(i0*ig_4(4)); ua5m=0.377*r*nmax/(i0*ig_5(5)); ua4=0:0.01:ua4m; ua5=0:0.01:ua5m;

Fw_4=(Cd*A*(ua4.^2))/21.15; Fw_5=(Cd*A*(ua5.^2))/21.15; Fz4_0=Ff_0+Fw_4; Fz4_1=Ff_1+Fw_4; Fz5_0=Ff_0+Fw_5; Fz5_1=Ff_1+Fw_5; %%figure;

plot(ua_4(1,:),1000*Ft_4(1,:),ua_4(2,:),1000*Ft_4(2,:),ua_4(3,:),1000*Ft_4(3,:),ua_4(4,:),1000*Ft_4(4,:),ua4,Fz4_0,ua4,Fz4_1);

title('4档变速汽车空载&满载驱动力行驶阻力平衡图 Ft Ff+Fw-ua'); xlabel('ua/(km/h)'); ylabel('F/N');

text(15.5,14000,'Ft1'); text(22,7700,'Ft2'); text(35,4500,'Ft3'); text(75,2700,'Ft4'); text(68,1650,'Ff+Fw'); %%figure;

plot(ua_5(1,:),1000*Ft_5(1,:),ua_5(2,:),1000*Ft_5(2,:),ua_5(3,:),1000*Ft_5(3,:),ua_5(4,:),1000*Ft_5(4,:),ua_5(5,:),1000*Ft_5(5,:),ua5,Fz5_0,ua5,Fz5_1);

title('5档变速汽车空载&满载驱动力行驶阻力平衡图 Ft Ff+Fw-ua'); xlabel('ua/(km/h)'); ylabel('Ft/N');

text(16,13000,'Ft1'); text(20,7000,'Ft2'); text(35,4300,'Ft3'); text(75,2800,'Ft4'); text(113,1100,'Ft5'); text(108,2800,'Ff+Fw');

? 由驱动力行驶阻力平衡图得出： 4档空载&满载最高车速为 94.93km/h

5档空载最高车速为 107km/h 满载为99.12km/h

%%4. 计算最高车速时的附着率

hg_0=0.7; hg_1=0.9; a_0=0.3; b_0=0.7; a_1=0.6; b_1=0.4; L=3.2;

ua4_max=94.93;

Cfig_4_0=((m0*g*b_0/L)*f+Cd*A*ua4_max^2/21.15)/(m0*g*a_0/L); Cfig_4_1=((m1*g*b_1/L)*f+Cd*A*ua4_max^2/21.15)/(m1*g*a_1/L); ua5_max0=107.00; ua5_max1=99.12;

Cfig_5_0=((m0*g*b_0/L)*f+Cd*A*ua5_max0^2/21.15)/(m0*g*a_0/L); Cfig_5_1=((m1*g*b_1/L)*f+Cd*A*ua5_max1^2/21.15)/(m1*g*a_1/L);

通过计算得出最高车速时的附着率： 4档空载附着率 0.7434 4档满载附着率0.1741 5档空载附着率0.9363 5档满载附着率0.1890

%%5. 计算空载和满载时最大爬坡度

D4_0=zeros(4,length(n)); D4_1=zeros(4,length(n)); D5_0=zeros(5,length(n)); D5_1=zeros(5,length(n)); for m=1:4

D4_0(m,:)=(1000*Ft_4(m,:)-(Cd*A*(ua_4(m,:).^2))/21.15)/(m0*g); D4_1(m,:)=(1000*Ft_4(m,:)-(Cd*A*(ua_4(m,:).^2))/21.15)/(m1*g); end

for m=1:5

D5_0(m,:)=(1000*Ft_5(m,:)-(Cd*A*(ua_5(m,:).^2))/21.15)/(m0*g); D5_1(m,:)=(1000*Ft_5(m,:)-(Cd*A*(ua_5(m,:).^2))/21.15)/(m1*g); end

i_4_0=zeros(4,length(n)); i_4_1=zeros(4,length(n)); i_5_0=zeros(5,length(n)); i_5_1=zeros(5,length(n)); for k=1:4

ylabel('a/(m/s^2)');

%%figure;

plot((1/3.6)*ua_4(1,:),aa_4_0(1,:),(1/3.6)*ua_4(2,:),aa_4_0(2,:),(1/3.6)*ua_4(3,:),aa_4_0(3,:),(1/3.6)*ua_4(4,:),aa_4_0(4,:));

title('4档汽车空载加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)'); %%figure;

plot((1/3.6)*ua_4(1,:),aa_4_1(1,:),(1/3.6)*ua_4(2,:),aa_4_1(2,:),(1/3.6)*ua_4(3,:),aa_4_1(3,:),(1/3.6)*ua_4(4,:),aa_4_1(4,:));

title('4档汽车满载加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)'); %%figure;

plot((1/3.6)*ua_5(1,:),aa_5_0(1,:),(1/3.6)*ua_5(2,:),aa_5_0(2,:),(1/3.6)*ua_5(3,:),aa_5_0(3,:),(1/3.6)*ua_5(4,:),aa_5_0(4,:),(1/3.6)*ua_5(5,:),aa_5_0(5,:)); axis([0 29.72 0 10]);

title('5档汽车空载加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)'); %%figure;

plot((1/3.6)*ua_5(1,:),aa_5_1(1,:),(1/3.6)*ua_5(2,:),aa_5_1(2,:),(1/3.6)*ua_5(3,:),aa_5_1(3,:),(1/3.6)*ua_5(4,:),aa_5_1(4,:),(1/3.6)*ua_5(5,:),aa_5_1(5,:));

axis([0 27.53 0 15]);

title('5档汽车满载加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)');

uas_4_0_2=find(ua_4(2,:)>=ua_4(1,340001)); uas_4_0_3=find(ua_4(3,:)>=ua_4(2,340001));

uas_4_0_4=find(ua_4(4,:)>=ua_4(3,340001) & ua_4(4,:)<=70); for i=1:340001; if(i

if(i

if(iuas_4_0_4(length(uas_4_0_4))) aa_4_0(4,i)=0; aa_4_1(4,i)=0; end end

t4_0=(1/3.6)*(trapz(ua_4(1,:),aa_4_0(1,:))+trapz(ua_4(2,:),aa_4_0(2,:))+trapz(ua_4(3,:),aa_4_0(3,:))+trapz(ua_4(4,:),aa_4_0(4,:)));

t4_1=(1/3.6)*(trapz(ua_4(1,:),aa_4_1(1,:))+trapz(ua_4(2,:),aa_4_1(2,:))+trapz(ua_4(3,:),aa_4_1(3,:))+trapz(ua_4(4,:),aa_4_1(4,:))); %%figure;

plot(ua_4(1,:),aa_4_0(1,:),ua_4(2,:),aa_4_0(2,:),ua_4(3,:),aa_4_0(3,:),ua_4(4,:),aa_4_0(4,:)); title('4档汽车空载最佳加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)'); %%figure;

plot(ua_4(1,:),aa_4_1(1,:),ua_4(2,:),aa_4_1(2,:),ua_4(3,:),aa_4_1(3,:),ua_4(4,:),aa_4_1(4,:)); title('4档汽车满载最佳加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)');

uas_5_0_2=find(ua_5(2,:)>=ua_5(1,340001)); uas_5_0_3=find(ua_5(3,:)>=ua_5(2,340001));

uas_5_0_4=find(ua_5(4,:)>=ua_5(3,340001) & ua_5(4,:)<=70); for i=1:340001; if(i

if(i

if(iuas_5_0_4(length(uas_5_0_4))) aa_5_0(4,i)=0; aa_5_1(4,i)=0; end end

t5_0=(1/3.6)*(trapz(ua_5(1,:),aa_5_0(1,:))+trapz(ua_5(2,:),aa_5_0(2,:))+trapz(ua_5(3,:),aa_5_0(3,:))+trapz(ua_5(4,:),aa_5_0(4,:)));

t5_1=(1/3.6)*(trapz(ua_5(1,:),aa_5_1(1,:))+trapz(ua_5(2,:),aa_5_1(2,:))+trapz(ua_5(3,:),aa_5_1(3,:))+trapz(ua_5(4,:),aa_5_1(4,:))); %%figure;

plot(ua_5(1,:),aa_5_0(1,:),ua_5(2,:),aa_5_0(2,:),ua_5(3,:),aa_5_0(3,:),ua_5(4,:),aa_5_0(4,:),ua_5(5,:),aa_5_0(5,:));

axis([0 107 0 5]);

title('5档汽车空载最佳加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)'); %%figure;

plot(ua_5(1,:),aa_5_1(1,:),ua_5(2,:),aa_5_1(2,:),ua_5(3,:),aa_5_1(3,:),ua_5(4,:),aa_5_1(4,:),ua_5(5,:),aa_5_1(5,:));

axis([0 99.12 0 10]);

title('5档汽车满载最佳加速度倒数曲线 1/a-ua'); xlabel('ua/(m/s)');

ylabel('1/a /(m/s^2)');

? 通过计算以上最佳换挡时刻图中曲线与坐标轴包围的面积得出：

4档空载加速到70km/h的时间为 11.3492s 4档满载加速到70km/h的时间为 25.7815s 5档空载加速到70km/h的时间为 11.1030s 5档满载加速到70km/h的时间为 25.1317s

%%8. 计算汽车功率平衡图

Pe_4=zeros(4,length(n)); Pe_5=zeros(5,length(n)); for i=1:4

Pe_4(i,:)=(1/3.6)*Ft_4(i,:).*ua_4(i,:)/elta; end

for i=1:5

Pe_5(i,:)=(1/3.6)*Ft_5(i,:).*ua_5(i,:)/elta; end

Pm_4=max(Pe_4(2,:)); x=0:0.01:100; y=Pm_4+x-x;

Pz4_0=Fz4_0.*ua4/(3600*elta); Pz4_1=Fz4_1.*ua4/(3600*elta); %%figure;

plot(ua_4(1,:),Pe_4(1,:),ua_4(2,:),Pe_4(2,:),ua_4(3,:),Pe_4(3,:),ua_4(4,:),Pe_4(4,:),x,y,ua4,Pz4_0,ua4,Pz4_1);

title('4档变速汽车功率平衡图 Pe-ua'); xlabel('ua/(km/h)'); ylabel('Pe/(kW)'); text(6,52,'Ⅰ'); text(17,52,'Ⅱ'); text(33,52,'Ⅲ');

text(61,52,'Ⅳ'); text(50,64,'Pe');

text(72,18,'(Pf+Pw)/elta') Pm_5=max(Pe_5(2,:)); x=0:0.01:120; y=Pm_5+x-x;

Pz5_0=Fz5_0.*ua5/(3600*elta); Pz5_1=Fz5_1.*ua5/(3600*elta); %%figure;

plot(ua_5(1,:),Pe_5(1,:),ua_5(2,:),Pe_5(2,:),ua_5(3,:),Pe_5(3,:),ua_5(4,:),Pe_5(4,:),ua_5(5,:),Pe_5(5,:),x,y,ua5,Pz5_0,ua5,Pz5_1); axis([0,120 0 70])

title('5档变速汽车功率平衡图 Pe-ua'); xlabel('ua/(km/h)'); ylabel('Pe/(kW)'); text(6,52,'Ⅰ'); text(17,52,'Ⅱ'); text(33,52,'Ⅲ'); text(61,52,'Ⅳ'); text(50,64,'Pe');

text(80,52,'Ⅴ');

text(76,20,'(Pf+Pw)/elta');

%%9. 计算最高挡的等速百公里油耗曲线

n0=[815 1207 1614 2012 2603 3006 3403 3804]; b00=[1326.8 -414.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.0002823;...

1129.7 -45.291 0.71113 -0.00075215 -0.000038568]; b_0=(b00(:,1))'; b_1=(b00(:,2))'; b_2=(b00(:,3))'; b_3=(b00(:,4))'; b_4=(b00(:,5))'; bb4=zeros(1,8);

uat_4=(0.377*n0.*r)/(i0*ig_4(4)); up=zeros(1,8); for i=1:8

up(i)=find(ua_4(4,:)>=uat_4(i),1,'first');

bb4(i)=b_0(i)+b_1(i)*Pe_4(4,up(i))+b_2(i)*(Pe_4(4,up(i)))^2+b_3(i)*(Pe_4(4,up(i)))^3+b_4(i)*(Pe_4(4,up(i)))^4; end

b4=myfunb(n,bb4,n0); ua00=[20 30 40 50 60 70]; uap4=zeros(1,6); for i=1:6;

uap4(i)=find(ua_4(4,:)>=ua00(i),1,'first'); end

Qs_4=zeros(1,6); for i=1:6

Qs_4(i)=Pe_4(4,uap4(i))*b4(uap4(i))/(1.02*ua00(i)*8.04); end;

%%figure;

plot(ua00,Qs_4,'x','markersize',10); axis([0 80 0 50]);

title('4档汽车最高档等速燃油消耗量散点图 Qs-ua'); xlabel('ua/(km/h)');

ylabel('Qs/[L (100km)^-1]'); bb5=zeros(1,8);

uat_5=(0.377*n0.*r)/(i0*ig_5(5)); up=zeros(1,8); for i=1:8

up(i)=find(ua_5(5,:)>=uat_5(i),1,'first');

bb5(i)=b_0(i)+b_1(i)*Pe_5(5,up(i))+b_2(i)*(Pe_5(5,up(i)))^2+b_3(i)*(Pe_5(5,up(i)))^3+b_4(i)*(Pe_5(5,up(i)))^4; end

b5=myfunb(n,bb5,n0); ua00=[20 30 40 50 60 70]; uap5=zeros(1,6); for i=1:6;

uap5(i)=find(ua_5(5,:)>=ua00(i),1,'first'); end

Qs_5=zeros(1,6); for i=1:6

Qs_5(i)=Pe_5(5,uap5(i))*b5(uap5(i))/(1.02*ua00(i)*8.04); end;

%%figure;

plot(ua00,Qs_5,'x','markersize',10); axis([0 80 0 50]);

title('5档汽车最高档等速燃油消耗量散点图 Qs-ua'); xlabel('ua/(km/h)');

ylabel('Qs/[L (100km)^-1]');

cftool; %%将两张散点图进行拟合(cubic spline)

%. 计算不同主传动比的经济性-动力性曲线

i__0=[5.17 5.43 5.83 6.17 6.33]; nn=length(n);

ua__41=zeros(4,nn);ua__42=zeros(4,nn);ua__43=zeros(4,nn);ua__44=zeros(4,nn);ua__45=zeros(4,nn); Ft__41=zeros(4,nn);Ft__42=zeros(4,nn);Ft__43=zeros(4,nn);Ft__44=zeros(4,nn);Ft__45=zeros(4,nn); ua__51=zeros(5,nn);ua__52=zeros(5,nn);ua__53=zeros(5,nn);ua__54=zeros(5,nn);ua__55=zeros(5,nn); Ft__51=zeros(5,nn);Ft__52=zeros(5,nn);Ft__53=zeros(5,nn);Ft__54=zeros(5,nn);Ft__55=zeros(5,nn); delta__4=zeros(5,4);delta__5=zeros(5,5);

a__41=zeros(4,nn);a__42=zeros(4,nn);a__43=zeros(4,nn);a__44=zeros(4,nn);a__45=zeros(4,nn); a__51=zeros(5,nn);a__52=zeros(5,nn);a__53=zeros(5,nn);a__54=zeros(5,nn);a__55=zeros(5,nn);

aa__41=zeros(4,nn);aa__42=zeros(4,nn);aa__43=zeros(4,nn);aa__44=zeros(4,nn);aa__45=zeros(4,nn); aa__51=zeros(5,nn);aa__52=zeros(5,nn);aa__53=zeros(5,nn);aa__54=zeros(5,nn);aa__55=zeros(5,nn);

for j=1:5

for i=1:4

delta__4(j,i)=1+Iw/(m0*r*r)+If*ig_4(i)^2*i__0(j)^2*elta/(m0*r*r); end end

for j=1:5

for i=1:5

delta__5(j,i)=1+Iw/(m0*r*r)+If*ig_5(i)^2*i__0(j)^2*elta/(m0*r*r); end end

for i=1:4

ua__41(i,:)=(0.377*n.*r)/(i__0(1)*ig_4(i)); ua__42(i,:)=(0.377*n.*r)/(i__0(2)*ig_4(i)); ua__43(i,:)=(0.377*n.*r)/(i__0(3)*ig_4(i)); ua__44(i,:)=(0.377*n.*r)/(i__0(4)*ig_4(i)); ua__45(i,:)=(0.377*n.*r)/(i__0(5)*ig_4(i)); Ft__41(i,:)=Tq*ig_4(i)*i__0(1)*elta/(r*1000); Ft__42(i,:)=Tq*ig_4(i)*i__0(2)*elta/(r*1000); Ft__43(i,:)=Tq*ig_4(i)*i__0(3)*elta/(r*1000); Ft__44(i,:)=Tq*ig_4(i)*i__0(4)*elta/(r*1000); Ft__45(i,:)=Tq*ig_4(i)*i__0(5)*elta/(r*1000);

a__41(i,:)=(1000*Ft__41(i,:)-Ff_0-Cd*A*(ua__41(i,:).^2)/21.15)/(m0*delta__4(1,i)); a__42(i,:)=(1000*Ft__42(i,:)-Ff_0-Cd*A*(ua__42(i,:).^2)/21.15)/(m0*delta__4(2,i)); a__43(i,:)=(1000*Ft__43(i,:)-Ff_0-Cd*A*(ua__43(i,:).^2)/21.15)/(m0*delta__4(3,i)); a__44(i,:)=(1000*Ft__44(i,:)-Ff_0-Cd*A*(ua__44(i,:).^2)/21.15)/(m0*delta__4(4,i)); a__45(i,:)=(1000*Ft__45(i,:)-Ff_0-Cd*A*(ua__45(i,:).^2)/21.15)/(m0*delta__4(5,i)); aa__41(i,:)=1./(a__41(i,:)); aa__42(i,:)=1./(a__42(i,:)); aa__43(i,:)=1./(a__43(i,:)); aa__44(i,:)=1./(a__44(i,:)); aa__45(i,:)=1./(a__45(i,:)); end

for i=1:5

ua__51(i,:)=(0.377*n.*r)/(i__0(1)*ig_5(i)); ua__52(i,:)=(0.377*n.*r)/(i__0(2)*ig_5(i)); ua__53(i,:)=(0.377*n.*r)/(i__0(3)*ig_5(i)); ua__54(i,:)=(0.377*n.*r)/(i__0(4)*ig_5(i)); ua__55(i,:)=(0.377*n.*r)/(i__0(5)*ig_5(i)); Ft__51(i,:)=Tq*ig_5(i)*i__0(1)*elta/(r*1000); Ft__52(i,:)=Tq*ig_5(i)*i__0(2)*elta/(r*1000); Ft__53(i,:)=Tq*ig_5(i)*i__0(3)*elta/(r*1000); Ft__54(i,:)=Tq*ig_5(i)*i__0(4)*elta/(r*1000); Ft__55(i,:)=Tq*ig_5(i)*i__0(5)*elta/(r*1000);

a__51(i,:)=(1000*Ft__51(i,:)-Ff_0-Cd*A*(ua__51(i,:).^2)/21.15)/(m0*delta__5(1,i)); a__52(i,:)=(1000*Ft__52(i,:)-Ff_0-Cd*A*(ua__52(i,:).^2)/21.15)/(m0*delta__5(2,i)); a__53(i,:)=(1000*Ft__53(i,:)-Ff_0-Cd*A*(ua__53(i,:).^2)/21.15)/(m0*delta__5(3,i)); a__54(i,:)=(1000*Ft__54(i,:)-Ff_0-Cd*A*(ua__54(i,:).^2)/21.15)/(m0*delta__5(4,i)); a__55(i,:)=(1000*Ft__55(i,:)-Ff_0-Cd*A*(ua__55(i,:).^2)/21.15)/(m0*delta__5(5,i)); aa__51(i,:)=1./(a__51(i,:)); aa__52(i,:)=1./(a__52(i,:)); aa__53(i,:)=1./(a__53(i,:)); aa__54(i,:)=1./(a__54(i,:)); aa__55(i,:)=1./(a__55(i,:)); end

%%figure;

plot(ua__41(1,:),aa__41(1,:),ua__41(2,:),aa__41(2,:),ua__41(3,:),aa__41(3,:),ua__41(4,:),aa__41(4,:)); axis([0 110 0 10]);

title('4档 i0=5.17时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__42(1,:),aa__42(1,:),ua__42(2,:),aa__42(2,:),ua__42(3,:),aa__42(3,:),ua__42(4,:),aa__42(4,:)); axis([0 110 0 10]);

title('4档 i0=5.43时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__43(1,:),aa__43(1,:),ua__43(2,:),aa__43(2,:),ua__43(3,:),aa__43(3,:),ua__43(4,:),aa__43(4,:)); axis([0 110 0 10]);

title('4档 i0=5.83时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__44(1,:),aa__44(1,:),ua__44(2,:),aa__44(2,:),ua__44(3,:),aa__44(3,:),ua__44(4,:),aa__44(4,:)); axis([0 110 0 10]);

title('4档 i0=6.17时的加速度倒数曲线');

xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__45(1,:),aa__45(1,:),ua__45(2,:),aa__45(2,:),ua__45(3,:),aa__45(3,:),ua__45(4,:),aa__45(4,:)); axis([0 110 0 10]);

title('4档 i0=6.33时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__51(1,:),aa__51(1,:),ua__51(2,:),aa__51(2,:),ua__51(3,:),aa__51(3,:),ua__51(4,:),aa__51(4,:),ua__51(5,:),aa__51(5,:)); axis([0 110 0 10]);

title('5档 i0=5.17时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__52(1,:),aa__52(1,:),ua__52(2,:),aa__52(2,:),ua__52(3,:),aa__52(3,:),ua__52(4,:),aa__52(4,:),ua__52(5,:),aa__52(5,:)); axis([0 110 0 10]);

title('5档 i0=5.43时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__53(1,:),aa__53(1,:),ua__53(2,:),aa__53(2,:),ua__53(3,:),aa__53(3,:),ua__53(4,:),aa__53(4,:),ua__53(5,:),aa__53(5,:)); axis([0 110 0 10]);

title('5档 i0=5.83时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__54(1,:),aa__54(1,:),ua__54(2,:),aa__54(2,:),ua__54(3,:),aa__54(3,:),ua__54(4,:),aa__54(4,:),ua__54(5,:),aa__54(5,:)); axis([0 110 0 10]);

title('5档 i0=6.17时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)'); %%figure;

plot(ua__55(1,:),aa__55(1,:),ua__55(2,:),aa__55(2,:),ua__55(3,:),aa__55(3,:),ua__55(4,:),aa__55(4,:),ua__55(5,:),aa__55(5,:)); axis([0 110 0 10]);

title('5档 i0=6.33时的加速度倒数曲线'); xlabel('ua/(km/h)'); ylabel('1/a /(m/s^2)');

本文来源：https://www.bwwdw.com/article/34j2.html