CAR_1_car_everything3lc_汽车_汽车动力性能评估_评估_源码

figure(1); n = [1000,1500,1750,2000,3000,4000,4300,4500,5000,6000,6300]; Ttq = [155,270,300,310,310,310,310,307,282,232,217]; plot(n,Ttq); axis([0 7000 100 350]); title('CAR\_1\_5挡手动发动机外特性中的转矩曲线'); xlabel('n/(r·min^{-1})'); ylabel('T_{tq}/(N·m)'); figure(2); r = 0.205*0.55+16*0.0254/2; m = 1900; G = m*9.8; Cd = 0.302; A = 2.19; igi0 = [13.778;7.09;4.322;2.976;2.279]; etaT = 0.9; ua = (0.377*r)./(igi0)*(n); Ft = (igi0*Ttq.*etaT)/r; c = jet(5); for i = 1:5; plot(ua(i,:),Ft(i,:),'color',c(i,:)); hold on; end ub = 0:1:325; f = 0.014*(1+(ub.^2)/19400); Ff = G*f; Fw = Cd*(ub.^2)*A/21.15; plot(ub,Ff+Fw,'m'); title('CAR\_1\_5挡手动驱动力-行驶阻力平衡图'); xlabel('u_a/(km·h^{-1})'); ylabel('F/N'); gtext('F_{t1}'),gtext('F_{t2}'),gtext('F_{t3}'),gtext('F_{t4}'),gtext('F_{t5}'),gtext('F_f+F_w'); figure(3); f = 0.014*(1+(ua.^2)/19400); Ff = G*f; Fw=(ua.^2).*Cd*A/21.15; F = Ft-Ff-Fw; a = F./m; for i = 1:5; plot(ua(i,:),a(i,:),'color',c(i,:)); hold on; end axis([0 250 0 7]); title('CAR\_1\_5挡手动的行驶加速度曲线'); xlabel('u_a/(km·h^{-1})'); ylabel('a/(m·s^{-2})'); gtext('Ⅰ'),gtext('Ⅱ'),gtext('Ⅲ'),gtext('Ⅳ'),gtext('Ⅴ'); figure(4); for i = 1:5; plot(ua(i,:),1./a(i,:),'color',c(i,:)); hold on; end axis([0 150 0 4]); title('CAR\_1\_5挡手动加速度倒数曲线'); xlabel('u_a/(km·h^{-1})'); ylabel('1/a/(s^2·m^{-1})'); gtext('Ⅰ'),gtext('Ⅱ'),gtext('Ⅲ'),gtext('Ⅳ'),gtext('Ⅴ'); figure(5); for i = 1:5; plot(ua(i,:),1./a(i,:),'color',c(i,:)); hold on; end axis([0 150 0 4]); title('CAR\_1\_5挡手动加速度倒数拟合曲线'); xlabel('u_a/(km·h^{-1})'); ylabel('1/a/(s^2·m^{-1})'); gtext('Ⅰ'),gtext('Ⅱ'),gtext('Ⅲ'),gtext('Ⅳ'),gtext('Ⅴ'); p1 = polyfit(ua(1,:),1./a(1,:),4); p2 = polyfit(ua(2,5:10),1./a(2,5:10),2); x1 = 0:0.01:8.65; x2 = 54.46:0.01:100; y1 = 6.8959e-07*x1.^4-9.4592e-05*x1.^3+0.0047.*x1.^2-0.0984.*x1+0.8748; y2 = 9.7120e-05.*x2.^2-0.0116.*x2+0.6645; plot(x1,y1,'--m'); plot(ua(1,:),1./a(1,:),'--m'); plot(x2,y2,'--m'); gtext('拟合曲线'); t1 = trapz(x1./3.6,y1); t2 = trapz(ua(1,:)./3.6,1./a(1,:)); t3 = trapz(x2./3.6,y2); t = t1+t2+t3; figure(6); alpha = asin(F/G); i1 = 100*tan(alpha); for i = 1:5; plot(ua(i,:),i1(i,:),'color',c(i,:)); hold on; end axis([0 250 0 100]); title('CAR\_1\_5挡手动的爬坡度图'); xlabel('u_a/(km·h^{-1})'); ylabel('i(%)'); gtext('Ⅰ'),gtext('Ⅱ'),gtext('Ⅲ'),gtext('Ⅳ'),gtext('Ⅴ');