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

《汽车理论大作业动力性能与经济性能计算Matlab》由会员分享，可在线阅读，更多相关《汽车理论大作业动力性能与经济性能计算Matlab（46页珍藏版）》请在装配图网上搜索。

1、汽车理论大作业 同济大学汽车学院xxxxxxxxx %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、%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);fori=1:4ua_4(i,:)=

3、(0.377*n.*r)/(i0*ig_4(i);Ft_4(i,:)=Tq*ig_4(i)*i0*elta/(r*1000);endfor j=1:5ua_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);t

4、ext(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)

5、; 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

6、=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,16

7、50,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(

8、75,2800,Ft4);text(113,1100,Ft5);text(108,2800,Ff+Fw); 由驱动力行驶阻力平衡图得出：4档空载&满载最高车速为 94.93km/h5档空载最高车速为 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_max2/21.15)/(m0*g*a_0/L);Cfig_4_1=(m1*g*b_1/L)*f+Cd*A*ua4_max2

9、/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_max02/21.15)/(m0*g*a_0/L);Cfig_5_1=(m1*g*b_1/L)*f+Cd*A*ua5_max12/21.15)/(m1*g*a_1/L);通过计算得出最高车速时的附着率：4档空载附着率 0.74344档满载附着率0.17415档空载附着率0.93635档满载附着率0.1890%5. 计算空载和满载时最大爬坡度D4_0=zeros(4,length(n);D4_1=zeros(4,length(n

10、);D5_0=zeros(5,length(n);D5_1=zeros(5,length(n);for m=1:4D4_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);endfor m=1:5D5_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

11、)/(m1*g);endi_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:4i_4_0(k,:)=100*tan(asin(D4_0(k,:)-f*sqrt(1-D4_0(k,:).2+f*f)/(1+f*f);i_4_1(k,:)=100*tan(asin(D4_1(k,:)-f*sqrt(1-D4_1(k,:).2+f*f)/(1+f*f);endfor k=1:5i_5_0(k,:)=100*tan(asin(D5_0(k,:)-f

12、*sqrt(1-D5_0(k,:).2+f*f)/(1+f*f);i_5_1(k,:)=100*tan(asin(D5_1(k,:)-f*sqrt(1-D5_1(k,:).2+f*f)/(1+f*f);endimax1=max(i_4_0(1,:);uai1=ua_4(1,(find(i_4_0(1,:)=imax1);%figure;plot(ua_4(1,:),i_4_0(1,:),ua_4(2,:),i_4_0(2,:),ua_4(3,:),i_4_0(3,:),ua_4(4,:),i_4_0(4,:);title(4档变速汽车空载爬坡度图i-ua);xlabel(ua/(km/h);y

13、label(i/(%);text(15,110,);text(20,47,);text(35,25,);text(75,12,);imax2=max(i_4_1(1,:);uai2=ua_4(1,(find(i_4_1(1,:)=imax2);%figure;plot(ua_4(1,:),i_4_1(1,:),ua_4(2,:),i_4_1(2,:),ua_4(3,:),i_4_1(3,:),ua_4(4,:),i_4_1(4,:);title(4档变速汽车满载爬坡度图i-ua);xlabel(ua/(km/h);ylabel(i/(%);text(17,33,);text(20,19,);t

14、ext(40,10,);text(80,4,);imax3=max(i_5_0(1,:);uai3=ua_5(1,(find(i_5_0(1,:)=imax3);%figure;plot(ua_5(1,:),i_5_0(1,:),ua_5(2,:),i_5_0(2,:),ua_5(3,:),i_5_0(3,:),ua_5(4,:),i_5_0(4,:),ua_5(5,:),i_5_0(5,:);title(5档变速汽车空载爬坡度图i-ua);xlabel(ua/(km/h);ylabel(i/(%);text(17,90,);text(20,42,);text(35,25,);text(75,

15、12,);text(110,1,);imax4=max(i_5_1(1,:);uai4=ua_5(1,(find(i_5_1(1,:)=imax4);%figure;plot(ua_5(1,:),i_5_1(1,:),ua_5(2,:),i_5_1(2,:),ua_5(3,:),i_5_1(3,:),ua_5(4,:),i_5_1(4,:),ua_5(5,:),i_5_1(5,:);title(5档变速汽车满载爬坡度图i-ua);xlabel(ua/(km/h);ylabel(i/(%);text(19,30,);text(20,18,);text(40,10,);text(80,4,);te

16、xt(110,1,); 由各种情况的爬坡度图得出：4档空载最大爬坡度：1.374档满载最大爬坡度为：0.3935档空载最大爬坡度：1.085档满载最大爬坡度为：0.353%6. 计算最大爬坡度时的附着率Cfid_4_0=(imax1/100)/(a_0/L+hg_0*imax1/(100*L);Cfid_4_1=(imax2/100)/(a_1/L+hg_1*imax2/(100*L);Cfid_5_0=(imax3/100)/(a_0/L+hg_0*imax3/(100*L);Cfid_5_1=(imax4/100)/(a_1/L+hg_1*imax4/(100*L); 通过计算得出最大爬坡

17、度时的附着率：4档空载附着率 3.48154档满载附着率1.31955档空载附着率3.27695档满载附着率1.2311%7. 计算空载和满载时用1挡起步加速到70km/h的最短加速时间Iw=2*0.899+4*0.8995;If=0.218;delta4_0=zeros(1,4);delta4_1=zeros(1,4);delta5_0=zeros(1,5);delta5_1=zeros(1,5);for q=1:4;delta4_0(q)=1+Iw/(m0*r*r)+If*ig_4(q)2*i02*elta/(m0*r*r);delta4_1(q)=1+Iw/(m1*r*r)+If*ig_

18、4(q)2*i02*elta/(m1*r*r);endfor q=1:5delta5_0(q)=1+Iw/(m0*r*r)+If*ig_5(q)2*i02*elta/(m0*r*r);delta5_1(q)=1+Iw/(m1*r*r)+If*ig_5(q)2*i02*elta/(m1*r*r);enda_4_0=zeros(4,length(n);a_4_1=zeros(4,length(n);a_5_0=zeros(5,length(n);a_5_1=zeros(5,length(n);aa_4_0=zeros(4,length(n);aa_4_1=zeros(4,length(n);aa_

19、5_0=zeros(5,length(n);aa_5_1=zeros(5,length(n);for j=1:4a_4_0(j,:)=(1000*Ft_4(j,:)-Ff_0-Cd*A*(ua_4(j,:).2)/21.15)/(m0*delta4_0(j);a_4_1(j,:)=(1000*Ft_4(j,:)-Ff_1-Cd*A*(ua_4(j,:).2)/21.15)/(m1*delta4_1(j);aa_4_0(j,:)=1./(a_4_0(j,:);aa_4_1(j,:)=1./(a_4_1(j,:);endfor j=1:5a_5_0(j,:)=(1000*Ft_5(j,:)-Ff_

20、0-Cd*A*(ua_5(j,:).2)/21.15)/(m0*delta5_0(j);a_5_1(j,:)=(1000*Ft_5(j,:)-Ff_1-Cd*A*(ua_5(j,:).2)/21.15)/(m1*delta5_1(j);aa_5_0(j,:)=1./(a_5_0(j,:);aa_5_1(j,:)=1./(a_5_1(j,:);end%figure;plot(ua_4(1,:),a_4_0(1,:),ua_4(2,:),a_4_0(2,:),ua_4(3,:),a_4_0(3,:),ua_4(4,:),a_4_0(4,:);title(4档汽车空载加速度曲线 a-ua);xlab

21、el(ua/(km/h);ylabel(a/(m/s2);plot(ua_4(1,:),a_4_1(1,:),ua_4(2,:),a_4_1(2,:),ua_4(3,:),a_4_1(3,:),ua_4(4,:),a_4_1(4,:);title(4档汽车满载加速度曲线 a-ua);xlabel(ua/(km/h);ylabel(a/(m/s2);plot(ua_5(1,:),a_5_0(1,:),ua_5(2,:),a_5_0(2,:),ua_5(3,:),a_5_0(3,:),ua_5(4,:),a_5_0(4,:),ua_5(5,:),a_5_0(5,:);axis(0 107 0 4)

22、;title(5档汽车空载加速度曲线 a-ua);xlabel(ua/(km/h);ylabel(a/(m/s2);plot(ua_5(1,:),a_5_1(1,:),ua_5(2,:),a_5_1(2,:),ua_5(3,:),a_5_1(3,:),ua_5(4,:),a_5_1(4,:),ua_5(5,:),a_5_1(5,:);axis(0 99.12 0 2.5);title(5档汽车满载加速度曲线 a-ua);xlabel(ua/(km/h);ylabel(a/(m/s2);%figure;plot(1/3.6)*ua_4(1,:),aa_4_0(1,:),(1/3.6)*ua_4(

23、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/s2);%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);

24、xlabel(ua/(m/s);ylabel(1/a /(m/s2);%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/s2);%figure;plot(1/3.6)*u

25、a_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/s2);uas_4_0_2=find(ua_4(2,:)=ua_4(1,340001);uas_4_0_3=find(ua_4(3,:)=ua_4(2,34000

26、1);uas_4_0_4=find(ua_4(4,:)=ua_4(3,340001) & ua_4(4,:)=70);fori=1:340001;if(iuas_4_0_2(1) aa_4_0(2,i)=0; aa_4_1(2,i)=0;endif(iuas_4_0_3(1) aa_4_0(3,i)=0; aa_4_1(3,i)=0;endif(iuas_4_0_4(length(uas_4_0_4) aa_4_0(4,i)=0; aa_4_1(4,i)=0;endendt4_0=(1/3.6)*(trapz(ua_4(1,:),aa_4_0(1,:)+trapz(ua_4(2,:),aa_4

27、_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档汽

28、车空载最佳加速度倒数曲线 1/a-ua);xlabel(ua/(m/s);ylabel(1/a /(m/s2);%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/s2);uas_5_0_2=find(ua_5(2,:)=ua_5(1,340001);uas_5_0_3=find(ua_5(3,:)=ua_5(2,340001

29、);uas_5_0_4=find(ua_5(4,:)=ua_5(3,340001) & ua_5(4,:)=70);fori=1:340001;if(iuas_5_0_2(1) aa_5_0(2,i)=0; aa_5_1(2,i)=0;endif(iuas_5_0_3(1) aa_5_0(3,i)=0; aa_5_1(3,i)=0;endif(iuas_5_0_4(length(uas_5_0_4) aa_5_0(4,i)=0; aa_5_1(4,i)=0;endendt5_0=(1/3.6)*(trapz(ua_5(1,:),aa_5_0(1,:)+trapz(ua_5(2,:),aa_5_

30、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,:),

31、aa_5_0(5,:);axis(0 107 0 5);title(5档汽车空载最佳加速度倒数曲线 1/a-ua);xlabel(ua/(m/s);ylabel(1/a /(m/s2);%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/s2);

32、 通过计算以上最佳换挡时刻图中曲线与坐标轴包围的面积得出：4档空载加速到70km/h的时间为 11.3492s4档满载加速到70km/h的时间为 25.7815s5档空载加速到70km/h的时间为 11.1030s5档满载加速到70km/h的时间为 25.1317s%8. 计算汽车功率平衡图Pe_4=zeros(4,length(n);Pe_5=zeros(5,length(n);fori=1:4 Pe_4(i,:)=(1/3.6)*Ft_4(i,:).*ua_4(i,:)/elta;endfori=1:5 Pe_5(i,:)=(1/3.6)*Ft_5(i,:).*ua_5(i,:)/elta

33、;endPm_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,5

34、2,);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

35、,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

36、.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

37、=(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);fori=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;endb4=myfunb(n,bb4,n0);

38、ua00=20 30 40 50 60 70;uap4=zeros(1,6);fori=1:6;uap4(i)=find(ua_4(4,:)=ua00(i),1,first);endQs_4=zeros(1,6);fori=1:6Qs_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);

39、bb5=zeros(1,8);uat_5=(0.377*n0.*r)/(i0*ig_5(5);up=zeros(1,8);fori=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;endb5=myfunb(n,bb5,n0);ua00=20 30 40 50 60 70;uap5=zeros(1,6);fori=1:6;uap5(i)=find(ua_5(5,

40、:)=ua00(i),1,first);endQs_5=zeros(1,6);fori=1:6Qs_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)%10. 计算不同主传动比的经济性-动力性曲线i_0=5.17 5.4

41、3 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);

42、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

43、_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:5fori=1:4delta_4(j,i)=1+Iw/(m0*r*r)+If*ig_4(i)2*i_0(j)2*elta/(m0*r*r);endendfor j=1:5fori=1:5delta_5(j,i)=1+Iw/(m0

44、*r*r)+If*ig_5(i)2*i_0(j)2*elta/(m0*r*r);endendfori=1:4ua_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_

45、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.

46、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,:);a

47、a_43(i,:)=1./(a_43(i,:);aa_44(i,:)=1./(a_44(i,:);aa_45(i,:)=1./(a_45(i,:);endfori=1:5ua_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*i

48、g_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,

49、:)-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(

50、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/s2);%figure;plo

51、t(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/s2);%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=

52、5.83时的加速度倒数曲线);xlabel(ua/(km/h);ylabel(1/a /(m/s2);%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/s2);%figure;plot(ua_45(1,:),aa_45(1,:),ua_45(2,:),aa_45(2,:),ua_45(3,:),

53、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/s2);%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);yla

54、bel(1/a /(m/s2);%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/s2);%figure;plot(ua_53(1,:),aa_53(1,:),ua_53(2,:),aa_53(2,:),ua_53(3,:),aa_53(3,:),ua

55、_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/s2);%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/(

56、km/h);ylabel(1/a /(m/s2);%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/s2);% 画出所有曲线之后发现70km/h时均在4档处！uar_412=find(ua_41(2,:)=ua_41(1,nn);uar_413=find(ua_41(3,:)=ua_41(2,nn);uar_414=find(ua_41(4,:)=ua_41(3,nn) & ua_41(4,:)=ua_42(1,nn);uar_423=find(ua_42(3,:)=ua_42(2,nn);uar_424=find(ua_42(4,:)=ua_42(3,nn) & ua_42(4,:)=70);uar_432=find(