汽车理论matlab编程,(百度)汽车理论matlab编程.ppt

画驱动力与行驶阻力平衡图的MATLAB源程序： ig=[6.09,3.09,1.71,1]; for i=1:4 n=600:4000; t=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; ua=0.377*0.367*n./ig(i)/5.83; u=ua.*1000/3600; pe=t.*ig(i)*5.83.*u./0.367; str=['r','m','b','c']; plot(ua,pe,str(i)); grid on hold on; end ua=0:120; pf=((3880*9.8*0.013*ua./3.6)+(2.77*ua.*ua.*ua./76.14))./0.85; plot(ua,pf,'g-') xlabel('速度(Km/h)'),ylabel('功率(Kw)') title('汽车功率平衡图') legend('一档发动机功率','二档发动机功率','三档发动机功率','四档发动机功率','汽车行驶阻力功率',2) 画汽车坡度曲线图的MATLAB源程序： ig=[6.09,3.09,1.71,1]; for i=1:4 n=600:1:4000; t=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; ft=t*ig(i)*5.83*0.85/0.367; ua=0.377*0.367*n/ig(i)/5.83; ff=3880*9.8*0.013+2.77/21.15*ua.*ua; m=ft-ff; I=tan(asin(m/(3880*9.8))); str=['r','m','b','c']; plot(ua,I,str(i)); grid on; hold on; end xlabel('速度(Km/h)'),ylabel('坡度') title('汽车爬坡度图') legend('一档爬坡度','二档爬坡度','三档爬坡度','四档爬坡度') 画汽车行驶加速度倒数曲线图的MATLAB源程序： ig=[6.09,3.09,1.71,1]; for i=1:4 n=600:1:4000; t=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; ft=t*ig(i)*5.83*0.85/0.367; ua=0.377*0.367*n/ig(i)/5.83; ff=3880*9.8*0.013+2.77/21.15*ua.*ua; mh=3880+5.4/0.367/0.367+0.218*5.83*5.83*ig(i)^2*0.85/0.367/0.367; a=(ft-ff)/mh; m=1./a; str=['r','m','b','c']; plot(ua,m,str(i)); grid on; hold on; end xlabel('速度(Km/h)'),ylabel('加速度倒数(s^2/m)') title('汽车行驶加速度倒数曲线') legend('一档加速度倒数','二档加速度倒数','三档加速度倒数','四档加速度倒数',2) 画汽车功率平衡图的MATLAB源程序： ig=[6.09,3.09,1.71,1]; for i=1:4 n=600:4000; t=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; ua=0.377*0.367*n./ig(i)/5.83;u=ua.*1000/3600; pe=t.*ig(i)*5.83.*u./0.367; str=['r','m','b','c'];plot(ua,pe,str(i)); grid on,hold on; end ua=0:120;pf=((3880*9.8*0.013*ua./3.6)+(2.77*ua.*ua.*ua./76.14))./0.85; plot(ua,pf,'g-'),xlabel('速度(Km/h)'),ylabel('功率(Kw)') title('汽车功率平衡图') legend('一档发动机功率','二档发动机功率','三档发动机功率','四档发动机功率','汽车行驶阻力功率',2)