汽车理论matlab作业
汽车理论习题Matlab程序
确定一轻型货车的动力性能货车可装用4挡或5挡变速器,任选其中的一种进行整车性能计算：1绘制汽车驱动力与行驶阻力平衡图;2求汽车最高车速,最大爬坡度及克服该坡度时相应的附着率;3绘制汽车行驶加速度倒数曲线,用图解积分法求汽车用2档起步加速行驶至70km/h的车速－时间曲线,或者用计算机求汽车用2档起步加速行驶至70km/h的加速时间;轻型货车的有关数据：汽油发动机使用外特性的Tq-n曲线的拟合公式为式中,Tq为发动机转矩N m;n为发动机转速r/min;发动机的最低转速nmin =600r/min,最高转速nmax=4000r/min;装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径传动系机械效率ηt=滚动阻力系数f=空气阻力系数×迎风面积C D A=主减速器传动比i0=飞轮转动惯量 I f=m2二前轮转动惯量I w1=m2四后轮转动惯量I w2=m2质心至前轴距离满载 a=质心高满载 hg=解：Matlab程序：1 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=600:10:4000;Tq=+n/1000n/1000.^2+n/1000.^n/1000.^4;m=3880;g=;nmin=600;nmax=4000;G=mg;ig= ;nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;Ft1=Tqig1i0nT/r;Ft2=Tqig2i0nT/r;Ft3=Tqig3i0nT/r;Ft4=Tqig4i0nT/r;Ft5=Tqig5i0nT/r;ua1=rn/ig1/i0;ua2=rn/ig2/i0;ua3=rn/ig3/i0;ua4=rn/ig4/i0;ua5=rn/ig5/i0;ua=0:5:120;Ff=Gf;Fw=CDAua.^2/;Fz=Ff+Fw;plotua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz;title'驱动力-行驶阻力平衡图';xlabel'uakm/s';ylabel'FtN';gtext'Ft1',gtext'Ft2',gtext'Ft3',gtext'Ft4',gtext'Ft5',gtext'Ff+Fw'; zoom on;x,y=ginput1;zoom off;disp'汽车最高车速=';dispx;disp'km/h';汽车最高车速=km/h2求汽车最大爬坡度程序：n=600:10:4000;Tq=+n/1000n/1000.^2+n/1000.^n/1000.^4;m=3880;g=;nmin=600;nmax=4000;G=mg;ig= ;nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;Ft1=Tqig1i0nT/r;ua1=rn/ig1/i0;Ff=Gf;Fw1=CDAua1.^2/;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100tanasinmaxFi1/G;disp'汽车最大爬坡度=';dispimax;disp'%';汽车最大爬坡度=%3求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=600:10:4000;Tq=+n/1000n/1000.^2+n/1000.^n/1000.^4;m=3880;g=;nmin=600;nmax=4000;G=mg;ig= ;nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;Ft1=Tqig1i0nT/r;Ft2=Tqig2i0nT/r;Ft3=Tqig3i0nT/r;Ft4=Tqig4i0nT/r;Ft5=Tqig5i0nT/r;ua1=rn/ig1/i0;ua2=rn/ig2/i0;ua3=rn/ig3/i0;ua4=rn/ig4/i0;ua5=rn/ig5/i0;Fw1=CDAua1.^2/;Fw2=CDAua2.^2/;Fw3=CDAua3.^2/;Fw4=CDAua4.^2/;Fw5=CDAua5.^2/;Ff=Gf;deta1=1+Iw1+Iw2/mr^2+Ifig1^2i0^2nT/mr^2;deta2=1+Iw1+Iw2/mr^2+Ifig2^2i0^2nT/mr^2;deta3=1+Iw1+Iw2/mr^2+Ifig3^2i0^2nT/mr^2;deta4=1+Iw1+Iw2/mr^2+Ifig4^2i0^2nT/mr^2;deta5=1+Iw1+Iw2/mr^2+Ifig5^2i0^2nT/mr^2;a1=Ft1-Ff-Fw1/deta1m;ad1=1./a1;a2=Ft2-Ff-Fw2/deta2m;ad2=1./a2;a3=Ft3-Ff-Fw3/deta3m;ad3=1./a3;a4=Ft4-Ff-Fw4/deta4m;ad4=1./a4;a5=Ft5-Ff-Fw5/deta5m;ad5=1./a5;plotua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5;axis0 99 0 10;title'汽车的加速度倒数曲线';xlabel'uakm/h';ylabel'1/a';gtext'1/a1';gtext'1/a2';gtext'1/a3';gtext'1/a4';gtext'1/a5'; a=maxa1;af=asinmaxFt1-Ff-Fw1/G;C=tanaf/a/L+hgtanaf/L;disp'假设后轮驱动,最大爬坡度相应的附着率=';dispC;假设后轮驱动,最大爬坡度相应的附着率=4 >>clearnT=;r=;f=;CDA=;i0=;If=;Iw1=;Iw2=;L=;a=;hg=;m=3880;g=;G=mg; ig= ;nmin=600;nmax=4000;u1=rnmin./ig/i0;u2=rnmax./ig/i0;deta=0ig;for i=1:5detai=1+Iw1+Iw2/mr^2+Ifigi^2i0^2nT/mr^2;endua=6::99;N=lengthua;n=0;Tq=0;Ft=0;inv_a=0ua;delta=0ua; Ff=Gf;Fw=CDAua.^2/;for i=1:Nk=i;if uai<=u22n=uaiig2i0/r/;Tq=+n/1000n/1000^2+n/1000^n/1000^4;Ft=Tqig2i0nT/r;inv_ai=deta2m/Ft-Ff-Fwi;deltai=inv_ai/;elseif uai<=u23n=uaiig3i0/r/;Tq=+n/1000n/1000^2+n/1000^n/1000^4;Ft=Tqig3i0nT/r;inv_ai=deta3m/Ft-Ff-Fwi;deltai=inv_ai/;elseif uai<=u24n=uaiig4i0/r/;Tq=+n/1000n/1000^2+n/1000^n/1000^4;Ft=Tqig4i0nT/r;inv_ai=deta4m/Ft-Ff-Fwi;deltai=inv_ai/;elsen=uaiig5i0/r/;Tq=+n/1000n/1000^2+n/1000^n/1000^4;Ft=Tqig5i0nT/r;inv_ai=deta5m/Ft-Ff-Fwi;deltai=inv_ai/;enda=delta1:k;ti=suma;endplott,ua;axis0 80 0 100;title'汽车2档原地起步换挡加速时间曲线';xlabel'时间ts';ylabel'速度uakm/h';>> ginputans =所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为已知货车装用汽油发动机的负荷特性与万有特性;负荷特性曲线的拟合公式为：44332210e e e e P B P B P B P B B b ++++=其中,b 为燃油消耗率g/kW h ；Pe 为发动机净功率kW ；拟合式中的系数随转速n 变化;怠速油耗s mL Q id /299.0=怠速转速400r/min;计算与绘制题中货车的1汽车功率平衡图;2最高档与次高档的等速百公里油耗曲线;或利用计算机求货车按JB3352-83规定的六工况循环行驶的百公里油耗;计算中确定燃油消耗值b 时,若发动机转速与负荷特性中给定的转速不相等,可由相邻转速的两根曲线用插值法求得;解：Matlab 程序：(1)汽车功率平衡图程序：clearn=600:10:4000;Tq=+n/1000n/1000.^2+n/1000.^n/1000.^4;m=3880;g=;G=mg;ig= ;nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;ua1=rn/ig1/i0;ua2=rn/ig2/i0;ua3=rn/ig3/i0;ua4=rn/ig4/i0;ua5=rn/ig5/i0;Pe1=Tq.ig1i0.ua1./3600r;Pe2=Tq.ig2i0.ua2./3600r;Pe3=Tq.ig3i0.ua3./3600r;Pe4=Tq.ig4i0.ua4./3600r;Pe5=Tq.ig5i0.ua5./3600r;ua=0::119;Ff=Gf;Fw=CDAua.^2/;Pf=Ffua/3600;Pw=Fw.ua/3600;Pe0=Pf+Pw./nT;Pe=maxPe1;plotua1,Pe1,ua2,Pe2,ua3,Pe3,ua4,Pe4,ua5,Pe5,ua,Pe0,ua,Pe;axis0 119 0 100;title'汽车功率平衡图';xlabel'uakm/h';ylabel'Pekw';gtext'1',gtext'2',gtext'3',gtext'4',gtext'5',gtext'Pf+Pw/et',gt ext'Pe';2最高档与次高档的等速百公里油耗曲线程序：clearn=600:1:4000;m=3880;g=;G=mg;ig= ;nT=;r=;f=;CDA=;i0=;L=;a=;hg=;If=;Iw1=;Iw2=;n0=815 1207 1614 2012 2603 3006 3403 3804;B00= ;B10= ;B20= ;B30= ;B40= ;B0=splinen0,B00,n;B1=splinen0,B10,n;B2=splinen0,B20,n;B3=splinen0,B30,n;B4=splinen0,B40,n;Ff=Gf;ua4=rn/ig4/i0;ua5=rn/ig5/i0;Fz4=Ff+CDAua4.^2/;Fz5=Ff+CDAua5.^2/;Pe4=Fz4.ua4./nT1000;Pe5=Fz5.ua5./nT1000;for i=1:1:3401b4i=B0i+B1iPe4i+B2iPe4i.^2+B3iPe4i.^3+B4iPe4i.^4; b5i=B0i+B1iPe5i+B2iPe5i.^2+B3iPe5i.^3+B4iPe5i.^4; endpg=;Q4=Pe4.b4./.ua4.pg;Q5=Pe5.b5./.ua5.pg;plotua4,Q4,ua5,Q5;axis0 100 10 30;title'最高档与次高档等速百公里油耗曲线';xlabel'uakm/h';ylabel'百公里油耗L/100km';gtext'4',gtext'5';i为、、、、时的燃油经济性—加速时改变题中轻型货车的主减速器传动比,做出i值对汽车性能的影响;间曲线,讨论不同解：Matlab程序：主程序：i0=,,,,; %输入主传动比的数据for i=1:1:5yi=jiasushijiani0i; %求加速时间endy;for i=1:1:5bi=youhaoi0i; %求对应i0的六工况百公里油耗endb;plotb,y,'+r'hold onb1=linspaceb1,b5,100;y1=splineb,y,b1; %三次样条插值plotb1,y1; %绘制燃油经济性-加速时间曲线title'燃油经济性—加速时间曲线';xlabel'百公里油耗L/100km';ylabel'加速时间s';gtext'i0=',gtext'i0=',gtext'i0=',gtext'i0=',gtext'i0=';子程序：1 function y=jiasushijiani0 %求加速时间的处理函数n1=linspace0,5000; %先求各个档位的驱动力nmax=4000;nmin=600;r=;yita=;CDA=;f=;G=3880;ig=,,,;%i0=for i=1:1:4 %i为档数uamaxi=chesunmax,r,igi,i0; %计算各个档位的最大速度与最小速度uamini=chesunmin,r,igi,i0;uai,:=linspaceuamini,uamaxi,100;ni,:=zhuansuuai,:,r,igi,i0; %计算各个档位的转速范围 Ttqi,:=zhuanjuni,:; %求出各档位的转矩范围 Fti,:=qudongliTtqi,:,igi,i0,yita,r; %求出驱动力Fi,:=fG+CDAuai,:.^2/; %求出滚动阻力和空气阻力的和deltai,:=1+++igi^2i0^2yita/3880r^2; %求转动质量换算系数ai,:=1./deltai,:.3880./Fti,:-Fi,:; %求出加速度F2i,:=Fti,:-Fi,:;end%下面分各个档位进行积分,求出加速时间temp11,:=ua2,:/;temp12,:=1./a2,:;n1=1;for j1=1:1:100if ua3,j1>maxua2,:&&ua3,j1<=70temp21,n1=ua3,j1/;temp22,n1=1./a3,j1;n1=n1+1;endendn2=1;for j1=1:1:100if ua4,j1>maxua3,:&&ua4,j1<=70;temp31,n2=ua4,j1/;temp32,n2=1./a4,j1;n2=n2+1;endendy=temp11,1temp12,1+qiujitemp11,:,temp12,:+qiujitemp21,:,temp22,:+qiuj itemp31,:,temp32,:;end2 function ua=chesun,r,ig,i0; %由转速计算车速ua=r.n/igi0;3 function n=zhuansuua,r,ig,i0; %求转速n=igi0.ua./r;end4 function y=zhuanjun; %求转矩函数y=+.n./1000.n./1000.^2+.n./1000.^.n./1000.^4;5 function y=qudongliTtq,ig,i0,yita,r; %求驱动力函数y=igi0yita.Ttq/r;end6 function p=qiujix0,y0 %求积分函数n0=sizex0;n=n02;x=linspacex01,x0n,200 ;y=splinex0,y0,x; %插值% figure;plotx,y;p=trapzx,y ;end7 %求不同i0下的六工况油耗function b=youhaoi0;global f G CDA yita m r If Iw1 Iw2 pg B0 B1 B2 B3 B4 n %声明全局变量ig=,,,;r=;yita=;CDA=;f=;%i0=;G=3880;If=;Iw1=;Iw2=;m=3880; %汽车的基本参数设定n0=815 1207 1614 2012 2603 3006 3403 3804;B00= ;B10= ;B20= ;B30= ;B40= ;n=600:1:4000;B0=splinen0,B00,n;B1=splinen0,B10,n;B2=splinen0,B20,n; %使用三次样条插值,保证曲线的光滑连续B3=splinen0,B30,n;B4=splinen0,B40,n;ua4=r.n./i0ig4; %求出发动机转速范围内对应的III、IV档车速F4=fG+CDAua4.^2/; %求出滚动阻力和空气阻力的和P_fw4=F4.ua4./yita1000; %求出阻力功率for i=1:1:3401 %用拟合公式求出各个燃油消耗率b4i=B0i+B1iP_fw4i+B2iP_fw4i^2+B3iP_fw4i^3+B4iP_fw4i^4;endpg=; %汽油的重度取Lua4_m=25,40,50; %匀速阶段的车速s_m=50,250,250; %每段匀速走过的距离b4_m=splineua4,b4,ua4_m; %插值得出对应速度的燃油消耗率F4_m=fG+CDAua4_m.^2/; %车速对应的阻力P_fw4_m=F4_m.ua4_m./yita1000; %发动机功率Q4_m=P_fw4_m.b4_m.s_m./102.ua4_m.pg ; Q4_a1=jiasu40,25,ig4,,ua4,i0; Q4_a2=jiasu50,40,ig4,,ua4,i0; Qid=;tid=;s=1075;Q_i=Qidtid; %求出减速阶段的燃油消耗量 Q4all=sumQ4_m+Q4_a1+Q4_a2+Q_i100/s; %IV 档六工况百公里燃油消耗量 b=Q4all;8加速阶段处理函数function q=jiasuumax,umin,ig,a,ua0,i0;global f G CDA yita m r If Iw1 Iw2 pg B0 B1 B2 B3 B4 n; %i0 ; ua1=umin:1:umax; %把速度范围以1km/h 为间隔进行划分 delta=1+Iw1+Iw2/mr^2+Ifig^2i0^2yita/mr^2;P0=Gf.ua0./3600+CDA.ua0.^3/76140+deltam.ua0/3600a/yita; P=Gf.ua1/3600+CDA.ua1.^3/76140+deltam.ua1/3600a/yita; dt=1/a ; %速度每增加1km/h 所需要的时间for i=1:1:3401 %重新利用拟合公式求出b 与ua 的关系 b0i=B0i+B1iP0i+B2iP0i^2+B3iP0i^3+B4iP0i^4; endb1=interp1ua0,b0,ua1; %插值出各个速度节点的燃油消耗率 Qt=P.b1./.pg; %求出各个速度节点的燃油消耗率 i1=sizeQt; i=i12;Qt1=Qt2:i-1;q=Qt1+Qtidt./2+sumQt1dt; %求该加速阶段的燃油消耗量2）求行驶车速Ua ＝30km/h,在ϕ＝路面上车轮不抱死的制动距离;计算时取制动系反应时间'2τ＝,制动减速度上升时间''2τ＝;3）求制动系前部管路损坏时汽车的制动距离s,制动系后部管路损坏时汽车的制动距离's ; 解：Matlab 程序： (1) 求利用附着系数曲线和制动效率曲线程序： cleark=4080;hgk=;Lk=;ak=;betak=;bk=Lk-ak;%空载时的参数 mm=9290;hgm=;Lm=;am=;betam=;bm=Lm-am;%满载时的参数 z=0::; figure1; fai=z;fai_fk=betakzLk./bk+zhgk;%空载时前轴的φffai_fm=betamzLm./bm+zhgm;%满载时前轴的φffai_rk=1-betakzLk./ak-zhgk;%空载时后轴的φrfai_rm=1-betamzLm./am-zhgm;%满载时后轴的φrplotz,fai_fk,'b--',z,fai_fm,'r',z,fai_rk,'b--',z,fai_rm,'r',z,fai,' k';title'利用附着系数与制动强度的关系曲线';xlabel'制动强度z/g';ylabel'利用附着系数φ';gtext'φr空载',gtext'φr满载',gtext'φ=z',gtext'φf空载',gtext'φf 满载';figure2;Efk=z./fai_fk100;%空载时前轴的制动效率Efm=z./fai_fm100;Erk=z./fai_rk100;Erm=z./fai_rm100;plotfai_fk,Efk,'b',fai_fm,Efm,'r',fai_rk,Erk,'b',fai_rm,Erm,'r'; axis0 1 0 100;title'前.后制动效率曲线';xlabel'附着系数φ';ylabel'制动效率%';gtext'Ef',gtext'Er',gtext'Er',gtext'满载',gtext'空载';(2)问和3问程序：clearmk=4080;hgk=;Lk=;ak=;betak=;bk=Lk-ak;%空载时的参数mm=9290;hgm=;Lm=;am=;betam=;bm=Lm-am;%满载时的参数z=0::1;fai_fk=betakzLk./bk+zhgk;%空载时前轴的φffai_fm=betamzLm./bm+zhgm;%满载时前轴的φffai_rk=1-betakzLk./ak-zhgk;%空载时后轴的φrfai_rm=1-betamzLm./am-zhgm;%满载时后轴的φrEfk=z./fai_fk100;%空载时前轴的制动效率Efm=z./fai_fm100;Erk=z./fai_rk100;Erm=z./fai_rm100;t1=;t2=;ua0=30;fai=;g=;ak1=Erk81gfai/100;am1=Erm81gfai/100;Sk1=t1+t2/2ua0/+ua0^2/ak1;%制动距离Sm1=t1+t2/2ua0/+ua0^2/am1;disp'空载时,汽车制动距离Sk1=';dispSk1;disp'满载时,汽车制动距离Sm1=';dispSm1;ak2=faigak/Lk+faihgk; am2=faigam/Lm+faihgm; ak3=faigbk/Lk-faihgk; am3=faigbm/Lk-faihgm;Sk2=t1+t2/2ua0/+ua0^2/ak2;%制动距离 Sm2=t1+t2/2ua0/+ua0^2/am2; Sk3=t1+t2/2ua0/+ua0^2/ak3; Sm3=t1+t2/2ua0/+ua0^2/am3;disp'空载时,前制动器损坏,汽车制动距离Sk2='; dispSk2;disp'满载时,前制动器损坏,汽车制动距离Sm2='; dispSm2;disp'空载时,后制动器损坏,汽车制动距离Sk3='; dispSk3;disp'满载时,后制动器损坏,汽车制动距离Sm3='; dispSm3;空载时,汽车制动距离Sk1=满载时,汽车制动距离Sm1=空载时,前制动器损坏,汽车制动距离Sk2=满载时,前制动器损坏,汽车制动距离Sm2=空载时,后制动器损坏,汽车制动距离Sk3=满载时,后制动器损坏,汽车制动距离Sm3=二自由度轿车模型的有关参数如下：总质量 m=绕Oz 轴转动惯量 23885m kg I z ⋅= 轴距 L= 质心至前轴距离 a= 质心至后轴距离 b=前轮总侧偏刚度 k 1=-62618N/rad 后轮总侧偏刚度 k 2=-110185N/rad 转向系总传动比 i=20 试求：1) 稳定性因数K 、特征车速u ch ; 2) 稳态横摆角速度增益曲线asr u -⎪⎭⎫δω、车速u=s 时的转向灵敏度sw rδω;3) 静态储备系数.,侧向加速度为时的前、后轮侧偏角绝对值之差21αα-与转弯半径的比值R/R 0R 0=15m;4) 车速u=s 时,瞬态响应的横摆角速度波动的固有圆频率0ω、阻尼比ζ、反应时间τ与峰值反应时间ε解：Matlab 程序： m=;Iz=3885;L=;a=;b=;k1=-62618;k2=-110185; i=20;g=;R0=15;u1=; K=ma/k2-b/k1/L^2; Uch=1/K^1/2;%特征车速disp'稳定性因数s^2/m^2K='; dispK;disp'特征车速m/sUch='; dispUch; u=0::30;S=u./L1+Ku.^2;%稳态横摆角速度增益 plotu,S;title'汽车稳态横摆角速度增益曲线'; xlabel'车速um/s';ylabel'稳态横摆角速度增益'; disp'u=s 时,转向灵敏度为'; dispS448;SM=k2/k1+k2-a/L; ay=g; A=KayL; B=L/R0; R=L/B-A;C=R/R0;%转弯半径比 disp'静态储备系数.='; dispSM;disp'侧向加速度为时前、后轮侧偏角绝对值之差rad a1-a2='; dispA;disp'侧向加速度为时转弯半径比值R/R0='; dispC;W0=L/u1k1k2/mIz1+Ku1^2^1/2;%固有圆频率D=-mk1a^2+k2b^2-Izk1+k2/2LmIzk1k21+Ku1^2^1/2;%阻尼比 t=atan1-D^2^1/2/-mu1aW0/Lk2-D/W01-D^2^1/2;%反应时间 E=atan1-D^2^1/2/D/W01-D^2^1/2+t;%峰值反应时间 disp'车速u=s 时的瞬态响应参数分别为：'; disp'横摆角速度波动的固有圆频率rad 为 '; dispW0;disp'阻尼比为'; dispD;disp'反应时间s 为'; dispt;disp'峰值反应时间s 为'; dispE;稳定性因数s^2/m^2K=特征车速m/sUch=u=s 时,转向灵敏度为静态储备系数.=侧向加速度为时前、后轮侧偏角绝对值之差rad a1-a2=侧向加速度为时转弯半径比值R/R0=车速u=s 时的瞬态响应参数分别为： 横摆角速度波动的固有圆频率rad 为阻尼比为反应时间s 为峰值反应时间s 为车身-车轮双质量系统参数：10,9,25.0,5.10====μγζHz f ;“人体-座椅”系统参数：25.0,3==s s Hz f ζ;车速s m u /20=,路面不平度系数()3801056.2m n G q -⨯=,参考空间频率n 0=;计算时频率步长Hz f 2.0=∆,计算频率点数180=N ;1) 计算并画出幅频特性q z /1、12/z z 、2/z q 和均方根值谱()f G z 1 、()f G z 2 、()f G a 谱图;进一步计算aw w a zz q L a 、、、、、σσσσ21 值 2) 改变“人体-座椅”系统参数：5.0~125.0,6~5.1==s s Hz f ζ;分析aw w L a 、值随s s f ζ、的变化;3) 分别改变车身-车轮双质量系统参数：5.0~125.0,3~25.00==ζHz f ,20~5,18~5.4==μγ;绘制GFd fd z /2σσσ、、 三个响应量均方根值随以上四个系统参数变化的曲线; 解：Matlab 程序 1问yps=;%阻尼比ζ gama=9;%刚度比γ mu=10;%质量比μfs=3;ypss=;g=;a0=10^-6;f0=; ua=20;Gqn0=10^-8;n0=;detaf=;N=180; f=detaf0:N;lamta=f/f0;lamtas=f/fs;Wf=0f;deta=1-lamta.^2.1+gama-1/mulamta.^2-1.^2+4yps^2lamta.^2.gama-1/mu+1la mta.^2.^2;z1_q=gamasqrt1-lamta.^2.^2+4yps^2lamta.^2./deta;z2_z1=sqrt1+4yps^2lamta.^2./1-lamta.^2.^2+4yps^2lamta.^2; p_z2=sqrt1+2ypsslamtas.^2./1-lamtas.^2.^2+2ypsslamtas.^2; z2_q=gamasqrt1+4yps^2lamta.^2./deta; p_q=p_z2.z2_q;jfg_Gqddf=4pi^2sqrtGqn0n0^2uaf; jfg_Gzdd1f=z1_q.jfg_Gqddf; jfg_Gzdd2f=z2_q.jfg_Gqddf; jfg_Gaf=p_q.jfg_Gqddf;sigmaqdd=sqrttrapzf,jfg_Gqddf.^2;%路面不平度加速度均方根值 sigmazdd1=sqrttrapzf,jfg_Gzdd1f.^2;%车轮加速度均方根值 sigmazdd2=sqrttrapzf,jfg_Gzdd2f.^2;%车身加速度均方根值 sigmaa=sqrttrapzf,jfg_Gaf.^2;%人体加速度均方根值 for i=1:N+1 if fi<=2 Wfi=; elseif fi<=4 Wfi=fi/4; elseif fi<= Wfi=1; elseWfi=fi; end endkk=Wf.^2.jfg_Gaf.^2;aw=sqrttrapzf,kk;%加权加速度均方根值 Law=20log10aw/a0;%加权振级disp'路面不平度加速度均方根值为';dispsigmaqdd; disp'车轮加速度均方根值为';dispsigmazdd1; disp'车身加速度均方根值为';dispsigmazdd2; disp'人体加速度均方根值为';dispsigmaa;disp'加权加速度均方根值为';dispaw;disp'加权振级';dispLaw;figure1plotf,z1_q,title'幅频特性|z1/q|, f=, ζ=,γ=9,μ=10',xlabel'激振频率f/Hz',ylabel'|z1/q|';figure2plotf,z2_z1,title'幅频特性|z2/z1|,f=, ζ=,γ=9,μ=10',xlabel'激振频率f/Hz',ylabel'|z2/z1|';figure3plotf,p_z2,title'幅频特性|p/z2|,fs=, ζs=',xlabel'激振频率f/Hz',ylabel'|p/z2|';figure4plotf,jfg_Gzdd1f,title'车轮加速度均方根值√Gz1f谱图',xlabel'激振频率f/Hz',ylabel'√Gz1f';figure5plotf,jfg_Gzdd2f,title'车身加速度均方根值√Gz2f谱图',xlabel'激振频率f/Hz',ylabel'√Gz2f';figure6plotf,jfg_Gaf,title'人体加速度均方根值√Gaf谱图',xlabel'激振频率f/Hz',ylabel'√Gaf';路面不平度加速度均方根值为车轮加速度均方根值为车身加速度均方根值为人体加速度均方根值为加权加速度均方根值为加权振级2问程序1：cleargama=9;%刚度比γmu=10;%质量比μf0=;g=;a0=10^-6;ua=20;Gqn0=10^-8;n0=;detaf=;N=180;f=detaf0:N;lamta=f/f0;Wf=0f;for i=1:N+1if fi<=2Wfi=;elseif fi<=4Wfi=fi/4;elseif fi<=Wfi=1;elseWfi=fi;endendfs=3;ypss=;ypss0=::;a=0ypss0;La=0ypss0;M=lengthypss0;for i=1:Myps=ypss0i;lamtas=f/fs;deta=1-lamta.^2.1+gama-1/mulamta.^2-1.^2+4yps^2lamta.^2.gama-1/mu+1la mta.^2.^2;p_z2=sqrt1+2ypsslamtas.^2./1-lamtas.^2.^2+2ypsslamtas.^2;z2_q=gamasqrt1+4yps^2lamta.^2./deta;p_q=p_z2.z2_q;jfg_Gqddf=4pi^2sqrtGqn0n0^2uaf;jfg_Gaf=p_q.jfg_Gqddf;kk=Wf.^2.jfg_Gaf.^2;awi=sqrttrapzf,kk;endLaw=20log10aw/a0;figure1plotypss0,aw;title'aw随ζs的变化',xlabel'“人体—座椅”系统的阻尼比ζs',ylabel'aw/ms^-2';figure2plotypss0,Law;title'Law随ζs的变化',xlabel'“人体—座椅”系统的阻尼比ζs',ylabel'Law/dB';程序2：clearyps=;%阻尼比ζgama=9;%刚度比γmu=10;%质量比μf0=;g=;a0=10^-6;ua=20;Gqn0=10^-8;n0=;detaf=;N=180;f=detaf0:N;lamta=f/f0;Wf=0f;for i=1:N+1if fi<=2Wfi=;elseif fi<=4Wfi=fi/4;elseif fi<=Wfi=1;elseWfi=fi;endendypss=;fs=::6;M=lengthfs;for i=1:Mfs0=fsi;lamtas=f/fs0;deta=1-lamta.^2.1+gama-1/mulamta.^2-1.^2+4yps^2lamta.^2.gama-1/mu+1la mta.^2.^2;p_z2=sqrt1+2ypsslamtas.^2./1-lamtas.^2.^2+2ypsslamtas.^2;z2_q=gamasqrt1+4yps^2lamta.^2./deta;p_q=p_z2.z2_q;jfg_Gqddf=4pi^2sqrtGqn0n0^2uaf;jfg_Gaf=p_q.jfg_Gqddf;kk=Wf.^2.jfg_Gaf.^2;awi=sqrttrapzf,kk;endLaw=20log10aw/a0;figure3plotfs,aw;title'aw随fs的变化',xlabel'“人体—座椅”系统的固有频率fs',ylabel'aw/ms^-2';figure4plotfs,Law;title'Law随fs的变化',xlabel'“人体—座椅”系统的固有频率fs',ylabel'Law/dB';3问程序1：clearfigure1fs=3;yps_s=;g=;ua=20;Gqn0=10^-8;n0=;detaf=;N=180;f0=;yps=;gama=9;mu=10;ff0=::3;sigmaz2=0ff0;sigmafd=0ff0;sigmaFd_G=0ff0;M=lengthff0;for i=1:Mf0=ff0i;f=detaf0:N;lamta=f/f0;lamtas=f/fs;deta=1-lamta.^2.1+gama-1/mulamta.^2-1.^2+4yps^2lamta.^2.gama-1/mu+1la mta.^2.^2;z2_qdot=2pifgama.sqrt1+4yps^2lamta.^2./deta;fd_qdot=gamalamta.^2./2pif+eps./sqrtdeta;Fd_Gqdot=2pifgama/g.sqrtlamta.^2/mu+1-1.^2+4yps^2lamta.^2./deta;Gq_dotf=4pi^2Gqn0n0^2ua;Gz2f=z2_qdot.^2Gq_dotf;Gfd_qf=fd_qdot.^2Gq_dotf;GFd_Gf=Fd_Gqdot.^2Gq_dotf;sigmaz2i=sqrttrapzf,Gz2f;sigmafdi=sqrttrapzf,Gfd_qf;sigmaFd_Gi=sqrttrapzf,GFd_Gf;if f0==sgmz2=sigmaz2i;sgmfd=sigmafdi;sgmFd_G=sigmaFd_Gi;endendsz2=20log10sigmaz2/sgmz2;sfd=20log10sigmafd/sgmfd;sFd_G=20log10sigmaFd_G/sgmFd_G;plotff0,sz2,'r-',ff0,sfd,'b-.',ff0,sFd_G,'k--';axis 3 -25 15;title'三个响应量均方根值随f0变化的曲线',xlabel'车身部分固有频率f0/Hz',ylabel'σz2/dB,σfd/dB,σFd/G/dB';程序2：clearfigure2fs=3;yps_s=;g=;ua=20;Gqn0=10^-8;n0=;detaf=;N=180;f0=;yps=;gama=9;mu=10;c= i=1:Myps=yps0i;f=detaf0:N;lamta=f/f0;lamtas=f/fs;deta=1-lamta.^2.1+gama-1/mulamta.^2-1.^2+4yps^2lamta.^2.gama-1/mu+1la mta.^2.^2;z2_qdot=2pifgama.sqrt1+4yps^2lamta.^2./deta;fd_qdot=gamalamta.^2./2pif+eps./sqrtdeta;Fd_Gqdot=2pifgama/g.sqrtlamta.^2/mu+1-1.^2+4yps^2lamta.^2./deta;Gq_dotf=4pi^2Gqn0n0^2ua;Gz2f=z2_qdot.^2Gq_dotf;Gfd_qf=fd_qdot.^2Gq_dotf;GFd_Gf=Fd_Gqdot.^2Gq_dotf;sigmaz2i=sqrttrapzf,Gz2f;sigmafdi=sqrttrapzf,Gfd_qf;sigmaFd_Gi=sqrttrapzf,GFd_Gf;if yps==sgmz2=sigmaz2i;sgmfd=sigmafdi;sgmFd_G=sigmaFd_Gi;endendsz2=20log10sigmaz2/sgmz2;sfd=20log10sigmafd/sgmfd;sFd_G=20log10sigmaFd_G/sgmFd_G;plotyps0,sz2,'r-',yps0,sfd,'b-.',yps0,sFd_G,'k--';axis -4 4;title'三个响应量均方根值随ζ变化的曲线',xlabel'车身部分阻尼比ζ',ylabel'σz2/dB,σfd/dB,σFd/G/dB';程序3：clearfigure3fs=3;yps_s=;g=;ua=20;Gqn0=10^-8;n0=;detaf=;N=180;f0=;yps=;mu=10;gama0=4::19;sigmaz2=0gama0;sigmafd=0gama0;sigmaFd_G=0gama0;M=lengthgama0;for i=1:Mgama=gama0i;f=detaf0:N;lamta=f/f0;lamtas=f/fs;deta=1-lamta.^2.1+gama-1/mulamta.^2-1.^2+4yps^2lamta.^2.gama-1/mu+1la mta.^2.^2;z2_qdot=2pifgama.sqrt1+4yps^2lamta.^2./deta;fd_qdot=gamalamta.^2./2pif+eps./sqrtdeta;Fd_Gqdot=2pifgama/g.sqrtlamta.^2/mu+1-1.^2+4yps^2lamta.^2./deta;Gq_dotf=4pi^2Gqn0n0^2ua;Gz2f=z2_qdot.^2Gq_dotf;Gfd_qf=fd_qdot.^2Gq_dotf;GFd_Gf=Fd_Gqdot.^2Gq_dotf;sigmaz2i=sqrttrapzf,Gz2f;sigmafdi=sqrttrapzf,Gfd_qf;sigmaFd_Gi=sqrttrapzf,GFd_Gf;if gama==9sgmz2=sigmaz2i;sgmfd=sigmafdi;sgmFd_G=sigmaFd_Gi;endendsz2=20log10sigmaz2/sgmz2;sfd=20log10sigmafd/sgmfd;sFd_G=20log10sigmaFd_G/sgmFd_G;plotgama0,sz2,'r-',gama0,sfd,'b-.',gama0,sFd_G,'k--';axis4 18 -5 6;title'三个响应量均方根值随γ变化的曲线',xlabel'悬架与轮胎的刚度比γ',ylabel'σz2/dB,σfd/dB,σFd/G/dB';程序4：clearfigure4fs=3;yps_s=;g=;ua=20;Gqn0=10^-8;n0=;detaf=;N=180;f0=;yps=;gama=9;mu0=5::20;sigmaz2=0mu0;sigmafd=0mu0;sigmaFd_G=0mu0;M=lengthmu0;for i=1:Mmu=mu0i;f=detaf0:N;lamta=f/f0;lamtas=f/fs;deta=1-lamta.^2.1+gama-1/mulamta.^2-1.^2+4yps^2lamta.^2.gama-1/mu+1la mta.^2.^2;z2_qdot=2pifgama.sqrt1+4yps^2lamta.^2./deta;fd_qdot=gamalamta.^2./2pif+eps./sqrtdeta;Fd_Gqdot=2pifgama/g.sqrtlamta.^2/mu+1-1.^2+4yps^2lamta.^2./deta;Gq_dotf=4pi^2Gqn0n0^2ua;Gz2f=z2_qdot.^2Gq_dotf;Gfd_qf=fd_qdot.^2Gq_dotf;GFd_Gf=Fd_Gqdot.^2Gq_dotf;sigmaz2i=sqrttrapzf,Gz2f;sigmafdi=sqrttrapzf,Gfd_qf;sigmaFd_Gi=sqrttrapzf,GFd_Gf;if mu==10sgmz2=sigmaz2i;sgmfd=sigmafdi;sgmFd_G=sigmaFd_Gi;endendsz2=20log10sigmaz2/sgmz2;sfd=20log10sigmafd/sgmfd;sFd_G=20log10sigmaFd_G/sgmFd_G;plotmu0,sz2,'r-',mu0,sfd,'b-.',mu0,sFd_G,'k--';axis5 20 -2 2;title'三个响应量均方根值随μ变化的曲线',xlabel'车身与车轮部分质量比μ',ylabel'σz2/dB,σfd/dB,σFd/G/dB';。
汽车理论课后作业答案MATLAB
汽车理论作业MA TLAB过程汽车驱动力与阻力平衡图加速度倒数-速度曲线图0 10 20 30 40 50 60 70 80 90 100u汽车功率平衡图u/(km/h)最高档等速百公里油耗曲线Ua/(km/h)燃油积极性-加速时间曲线源程序：《第一章》m=3880;g=9.8;r=0.367;x=0.85;f=0.013;io=5.83;CdA=2.77;lf=0.218;Iw1=1.798;Iw2=3.598;Iw=lw1+lw2;ig=[6.09 3.09 1.71 1.00]; %变速器传动比L=3.2;a=1.947;hg=0.9;n=600:1:4000;T=-19.313+295.27* n/1000-165.44*(门/1000)人2+40.874*(门/1000)人3-3.8445*( n/IOOO).%;Ft1=T*ig(1)*io*x/r; %计算各档对应转速下的驱动力Ft2=T*ig(2)*io*x/r;Ft3=T*ig(3)*io*x/r;Ft4=T*ig(4)*io*x/r;u1=0.377*r*n/(io*ig(1));u2=0.377*r*n/(io*ig(2));u3=0.377*r*n/(io*ig(3));u4=0.377*r*n/(io*ig(4));u=0:130/3400:130;F仁m*g*f+CdA*u”2/21.15;%计算各档对应转速下的驱动阻力F2=m*g*f+CdA*u2.A2/21.15;F3=m*g*f+CdA*u3.A2/21.15;F4=m*g*f+CdA*u4.A2/21.15;figure(1);plot(u1,Ft1, '-r' ,u2,Ft2, '-m' ,u3,Ft3, '-k' ,u4,Ft4, '-b' ,u1,F1, '-r' ,u2,F2, '-m' ,u3,F3, ' k' ,u4,F4, '-b' , 'LineWidth' ,2) title( ' 汽车驱动力与阻力平衡图' );xlabel( 'u_{a}/km.hA{-1}' )ylabel( 'F/N' )gtext( 'F_{t1}' )gtext( 'F_{t2}' )gtext( 'F_{t3}' )gtext( 'F_{t4}' )gtext( 'F_{f}+F_{w}' )%由汽车驱动力与阻力平衡图知，他们无交点，u4在最大转速时达到最大umax=u4(3401)Ft1max=max(Ft1);imax=(Ft1max-m*g*f)/(m*g)disp( ' 假设是后轮驱动' );C=imax/(a/L+hg*imax/L) % 附着率delta1=1+(Iw1+Iw2)/(m*rA2)+If*ig(1)*rA2*ioA2*x/(m*rA2);delta2=1+(Iw1+Iw2)/(m*rA2)+If*ig(2)*rA2*ioA2*x/(m*rA2);delta3=1+(Iw1+Iw2)/(m*rA2)+If*ig(3)*rA2*ioA2*x/(m*rA2);delta4=1+(Iw1+Iw2)/(m*rA2)+If*ig(4)*rA2*ioA2*x/(m*rA2);a1=(Ft1-F1)/(delta1*m); %加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);h1=1./a1; %加速度倒数h2=1./a2;h3=1./a3;h4=1./a4;figure(2);plot(u1,h1,u2,h2,u3,h3,u4,h4, 'LineWidth',2);title( ' 加速度倒数- 速度曲线图xlabel( 'u' )ylabel( '1/a' )gtext( '1/a1' )gtext( '1/a2' )gtext( '1/a3' )gtext( '1/a4' )%由加速度倒数- 速度曲线图可知u1min=min(u1);u1max=max(u1); u2min=u1max; u2min=min(u2);u2max=max(u2); u3min=u2max;u3max=max(u3); u4min=u3max; u4max=70;x1=[]; x2=[]; x3=[]; x4=[];y=3401; for i=1:3401;if u3(i)<=u3min;x1=[i];end end q1=max(x1); ua3=u3(q1:y); a3=h3(q1:y); for i=1:3401;if u4(i)<=u4min;x2=[i];elseif u4(i)<=u4max;x3=[i];end end q2=max(x2); q3=max(x3); ua4=u4(q2:q3); a4=h4(q2:q3); s1=trapz(h2,u2 ); %二挡运行时间s2=trapz(ua3,a3);s3=trapz(ua4,a4);s=[s1 s2 s3];disp( ' 积分得' ) t=sum(s)*1000/3600 %总时间《第二章》Pe1=Ft1.*u1./3600; %计算各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;P1=F1.*u1./(3600*x); %计算各档对应的各个车速下的行驶功率P2=F2.*u2./ (3600*x);P3=F3.*u3./ (3600*x);P4=F4.*u4./ (3600*x);figure(3);plot(u1,Pe1, '-r' ,u2,Pe2, '-m' ,u3,Pe3, '-k' ,u4,Pe4, '-b' ,u1,P1, 'k' ,u2,P2, 'k' ,u3,P3, 'k' u4,P4, 'k' ,'linewidth' ,2);gtext( 'Pe1' )gtext( 'Pe2' )gtext( 'Pe3' )gtext( 'Pe4' )xlabel( 'u/(km/h)' );ylabel( 'P/kW' );title( ' 汽车功率平衡图' );n=[815 1207 1614 2012 2603 3006 3403 3804];Ua=[];Ua=0.377*r*n./(io*ig(4))ft=[];ft=m*g*f+(2.77/21.15)*Ua.A2; 算各转速对应的各个车速下的行驶阻力Pe(1)=ft(1).*Ua(1)./(3600*x); %计算各转速对应的各个车速下的行驶阻力功率Pe(2)=ft(2).*Ua(2)./(3600*x);Pe(3)=ft(3).*Ua(3)./(3600*x);Pe(4)=ft(4).*Ua(4)./(3600*x);Pe(5)=ft(5).*Ua(5)./(3600*x);Pe(6)=ft(6).*Ua(6)./(3600*x);Pe(7)=ft(7).*Ua(7)./(3600*x);Pe(8)=ft(8).*Ua(8)./(3600*x)B0=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7];B1=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291];B2=[72.739 36.657 14.525 7.0035 4.4763 2.8593 2.9788 0.7113];B3=[-5.8629 -2.0533 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215]; B4=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568];b1=(B0⑴)+(B1(1)*Pe(1))+(B2(1)*Pe(1f2)+(B3(1)*Pe⑴T)+(B4(1)*Pe⑴人4); b2=(B0(2))+(B1(2)*Pe (2))+(B2(2)*Pe(2)A2)+(B3(2)*Pe(2)A3)+(B4(2)*Pe(2)A4);b3=(B 0(3) )+(B1(3)*Pe(3))+(B2 (3)*Pe(3)A2)+(B3(3)*Pe(3)A3)+(B4(3)*Pe(3)A4);b4=(B0(4))+(B1(4)*Pe(4))+(B2(4)*Pe(4)A2)+(B3(4)*Pe(4)A3)+(B4(4)*Pe(4)A4);b5=(B0(5))+(B1(5)*Pe(5))+(B2(5)*Pe(5)A2)+(B3(5)*Pe(5)A3)+(B4(5)*Pe(5)A4);b6=(B0(6))+(B1(6)*Pe(6))+(B2(6)*Pe(6)A2)+(B3(6)*Pe(6)A3)+(B4(6)*Pe(6)A4);b7=(B0(7))+(B1(7)*Pe(7))+(B2(7)*Pe(7)A2)+(B3(7)*Pe(7)A3)+(B4(7)*Pe(7)A4);b8=(B0(8))+(B1(8)*Pe(8))+(B2(8)*Pe(8)A2)+(B3(8)*Pe(8)A3)+(B4(8)*Pe(8)A4);p=0.7;Qs=[];Qs(1)=(Pe(1)*b1)/(1.02*Ua(1).*p*g);Qs(2)=(Pe(2)*b2)/(1.02*Ua(2).*p*g);Qs(3)=(Pe(3)*b3)/(1.02*Ua(3).*p*g);Qs(4)=(Pe(4)*b4)/(1.02*Ua(4).*p*g);Qs(5)=(Pe(5)*b5)/(1.02*Ua(5).*p*g);Qs(6)=(Pe(6)*b6)/(1.02*Ua(6).*p*g);Qs(7)=(Pe(7)*b7)/(1.02*Ua(7).*p*g);Qs(8)=(Pe(8)*b8)/(1.02*Ua(8).*p*g);M=polyfit(Ua,Qs,2);UA=0.377*r*600/(ig(4)*io):1:0.377*r*4000/(ig 4)*io); %卩表示车速QS=polyval(M,UA); %0表示油耗figure(4);plot(UA,QS, 'linewidth' ,2);title( '最高档等速百公里油耗曲线' );xlabel( 'Ua/(km/h)' );ylabel( 'Qs/L' );《第三章》io=[5.17 5.43 5.83 6.17 6.33];Va=0.377*r*n(7)./(io.*ig(4));Ps=46.9366;disp('假设以最高档，较高转速(n取3401 )，最经济负荷(即90%负荷大约46.9366KW )行驶时油耗' ) qs=[];qs(1)=(Ps*b7)/(1.02*Va(1).*p*g); qs(2)=(Ps*b7)/(1.02*Va(2).*p*g);qs(3)=(Ps*b7)/(1.02*Va(3).*p*g);qs(4)=(Ps*b7)/(1.02*Va(4).*p*g); qs(5)=(Ps*b7)/(1.02*Va(5).*p*g);st=[ 17.5813 16.2121 14.5126 13.3775 12.9185]; %加速时间：(这里以最高档〈四档〉、速度由0加速到94.93Km/h 的时间) 因与题 1.3 第三问求法相同，这里不在累述，可直接有计算机求得：figure(5);plot(qs,st, '+' ,'lineWidth' ,2)hold on plot(qs,st,); gtext( 517') gtext( 543') gtext( 583') gtext( 617') gtext( 633')title('燃油积极性-加速时间曲线’); xlabel('燃油经济性(qs/L)'); ylabel('动力性--原地起步加速时间(st/s)');《第四章》4-31)时：前轮先抱死 0时：后轮先抱死利用MATLAB 作图得到下图:前轴利用附着系数为:L z b zh g后轴利用附着系数为:L 1 zazh g空载时： 03.95 0.38 1.85所以满载时：ohg0.8450.413空载时后轮总是先抱死。
汽车理论作业-Matlab程序-轻型货车动力性能评价
轻型货车动力学性能评价4209XX班姓名：XX 学号：XXXXXXXX 已知条件该轻型货车相关参数如下：装载质量2000kg整车装备质量1800kg总质量（总重G）3800kg（38062.8 N）车轮半径（r）0.367 m传动系机械效率滚动阻力系数 f =0.013空气阻力系数X 迎风面积主减速器传动比飞轮转动惯量二前轮转动惯量四后轮转动惯量轴距质心至前轴的距离（满载）质心高变速器传动比（数据见下表）挡位1挡2挡3挡4挡5挡传动比 5.56 2.769 1.644 1.00 0.793发动机最低转速为 ，最高转速为汽车发动机使用外特性的曲线拟合公式为：驱动力与行驶阻力平衡图车速 [km/h]滚动阻力：空气阻力：行驶阻力：[N]驱动力 ：利用matlab 绘出出每一挡的驱动力与行驶阻力平衡图如下：2040608010012002000400060008000100001200014000X: 99.06Y: 1780ua /（km/h ）F /N汽车驱动力—行驶阻力平衡图Ft5Ft3Ft4Ft2Ft1最高车速点Fw+Ff图 1汽车驱动力与行驶阻力平衡图最高车速和最大爬坡度最高车速点：图 1汽车驱动力与行驶阻力平衡图中五挡驱动力曲线与行驶阻力曲线的交点。
由驱动力行驶阻力平衡方程：可利用matlab求解出求解出最高车速为：●最大爬坡度的计算最大坡度角满足方程：显然一挡的爬坡度最大，，可求解出但过于复杂，可对上式进行分析可知，当驱动力与空气阻力之差最大时，剩余的驱动力全部用来克服道路阻力，此时爬坡度最大。
可利用matlab先找到该车所能克服的道路阻力的最大值为：13125.7359567165 N那么最大动力因数为从而可求得最大坡度角为那么最大爬坡度为●克服最大爬坡度时相应的附着率计算1.计算静态轴荷的法向反作用力25400.85N2.动态分量3.空气升力由于此时车速很小，空气升力可忽略不计4.滚动阻力偶产生的部分也很小，可忽略不计。
汽车理论MATLAB编程
ad=(deta(2)*M)./(Ft-Fwf(i))
delta(i)=0.1*ad/3.6
elseif ua(i)<=uamax(3)
n=ua(i)*ig(3)*io/r/0.377
ua=0.377*r*n/ig/io
deta=1+(Iw1+Iw2)/(M*r^2)+If*ig.^2*io^2*m/(M*r^2)%汽车旋转质量换算系数
Fw=CdA/21.15*(ua).^2
Ttq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4
ua=0.377*r*n/ig/io
deta=1+(Iw1+Iw2)/(M*r^2)+If*ig.^2*io^2*m/(M*r^2)%汽车旋转质量换算系数
Fw=CdA/21.15*(ua).^2
Fwf=Fw+Ff
Ttq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4
delta(i)=0.1*ad/3.6
elseif ua(i)<=uamax(4)
n=ua(i)*ig(4)*io/r/0.377
Ttq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4
axis([0,80,0,100])
(完整版)汽车理论课后作业答案MATLAB
汽车理论作业 MA TLAB 过程010203040506070809010050001000015000汽车驱动力与阻力平衡图u a /km.h -1F /N10203040506070809010002468101214加速度倒数-速度曲线图u1/a0102030405060708090100102030405060u/(km/h)P /k W汽车功率平衡图10203040506070809010012141618202224最高档等速百公里油耗曲线Ua/(km/h)Q s /L2324252627282912131415161718燃油积极性-加速时间曲线燃油经济性(qs/L)动力性--原地起步加速时间 (s t /s )源程序：《第一章》m=3880; g=9.8; r=0.367; x=0.85; f=0.013; io=5.83; CdA=2.77; If=0.218; Iw1=1.798; Iw2=3.598; Iw=Iw1+Iw2;ig=[6.09 3.09 1.71 1.00]; %变速器传动比 L=3.2; a=1.947; hg=0.9; 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; Ft1=T*ig(1)*io*x/r;%计算各档对应转速下的驱动力Ft2=T*ig(2)*io*x/r;Ft3=T*ig(3)*io*x/r;Ft4=T*ig(4)*io*x/r;u1=0.377*r*n/(io*ig(1));u2=0.377*r*n/(io*ig(2));u3=0.377*r*n/(io*ig(3));u4=0.377*r*n/(io*ig(4));u=0:130/3400:130;F1=m*g*f+CdA*u1.^2/21.15;%计算各档对应转速下的驱动阻力F2=m*g*f+CdA*u2.^2/21.15;F3=m*g*f+CdA*u3.^2/21.15;F4=m*g*f+CdA*u4.^2/21.15;figure(1);plot(u1,Ft1,'-r',u2,Ft2,'-m',u3,Ft3,'-k',u4,Ft4,'-b',u1,F1,'-r',u2,F2,'-m',u3,F3,'-k',u4,F4,'-b','LineWidth',2)title('汽车驱动力与阻力平衡图');xlabel('u_{a}/km.h^{-1}')ylabel('F/N')gtext('F_{t1}')gtext('F_{t2}')gtext('F_{t3}')gtext('F_{t4}')gtext('F_{f}+F_{w}')%由汽车驱动力与阻力平衡图知，他们无交点，u4在最大转速时达到最大umax=u4(3401)Ft1max=max(Ft1);imax=(Ft1max-m*g*f)/(m*g)disp('假设是后轮驱动');C=imax/(a/L+hg*imax/L) % 附着率delta1=1+(Iw1+Iw2)/(m*r^2)+If*ig(1)*r^2*io^2*x/(m*r^2);delta2=1+(Iw1+Iw2)/(m*r^2)+If*ig(2)*r^2*io^2*x/(m*r^2);delta3=1+(Iw1+Iw2)/(m*r^2)+If*ig(3)*r^2*io^2*x/(m*r^2);delta4=1+(Iw1+Iw2)/(m*r^2)+If*ig(4)*r^2*io^2*x/(m*r^2);a1=(Ft1-F1)/(delta1*m); %加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);h1=1./a1; %加速度倒数h2=1./a2;h3=1./a3;h4=1./a4;figure(2);plot(u1,h1,u2,h2,u3,h3,u4,h4,'LineWidth',2); title('加速度倒数-速度曲线图');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')%由加速度倒数-速度曲线图可知u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=70;x1=[];x2=[];x3=[];x4=[];y=3401;for i=1:3401;if u3(i)<=u3min;x1=[i];endendq1=max(x1);ua3=u3(q1:y);a3=h3(q1:y);for i=1:3401;if u4(i)<=u4min;x2=[i];elseif u4(i)<=u4max;x3=[i];endendq2=max(x2);q3=max(x3);ua4=u4(q2:q3);a4=h4(q2:q3);s1=trapz(h2,u2 ); %二挡运行时间s2=trapz(ua3,a3);s3=trapz(ua4,a4);s=[s1 s2 s3];disp('积分得')t=sum(s)*1000/3600 %总时间《第二章》Pe1=Ft1.*u1./3600;%计算各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;P1=F1.*u1./(3600*x);%计算各档对应的各个车速下的行驶功率P2=F2.*u2./ (3600*x);P3=F3.*u3./ (3600*x);P4=F4.*u4./ (3600*x);figure(3);plot(u1,Pe1,'-r',u2,Pe2,'-m',u3,Pe3,'-k',u4,Pe4,'-b',u1,P1,'k',u2,P2,'k',u3,P3,'k', u4,P4,'k','linewidth',2);gtext('Pe1')gtext('Pe2')gtext('Pe3')gtext('Pe4')xlabel('u/(km/h)');ylabel('P/kW');title('汽车功率平衡图');n=[815 1207 1614 2012 2603 3006 3403 3804];Ua=[];Ua=0.377*r*n./(io*ig(4))ft=[];ft=m*g*f+(2.77/21.15)*Ua.^2;%计算各转速对应的各个车速下的行驶阻力Pe(1)=ft(1).*Ua(1)./(3600*x);%计算各转速对应的各个车速下的行驶阻力功率Pe(2)=ft(2).*Ua(2)./(3600*x);Pe(3)=ft(3).*Ua(3)./(3600*x);Pe(4)=ft(4).*Ua(4)./(3600*x);Pe(5)=ft(5).*Ua(5)./(3600*x);Pe(6)=ft(6).*Ua(6)./(3600*x);Pe(7)=ft(7).*Ua(7)./(3600*x);Pe(8)=ft(8).*Ua(8)./(3600*x)B0=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7];B1=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291];B2=[72.739 36.657 14.525 7.0035 4.4763 2.8593 2.9788 0.7113];B3=[-5.8629 -2.0533 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215];B4=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230-0.000038568];b1=(B0(1))+(B1(1)*Pe(1))+(B2(1)*Pe(1)^2)+(B3(1)*Pe(1)^3)+(B4(1)*Pe(1)^4);b2=(B0(2))+(B1(2)*Pe(2))+(B2(2)*Pe(2)^2)+(B3(2)*Pe(2)^3)+(B4(2)*Pe(2)^4);b3=(B0(3))+(B1(3)*Pe(3))+(B2(3)*Pe(3)^2)+(B3(3)*Pe(3)^3)+(B4(3)*Pe(3)^4);b4=(B0(4))+(B1(4)*Pe(4))+(B2(4)*Pe(4)^2)+(B3(4)*Pe(4)^3)+(B4(4)*Pe(4)^4);b5=(B0(5))+(B1(5)*Pe(5))+(B2(5)*Pe(5)^2)+(B3(5)*Pe(5)^3)+(B4(5)*Pe(5)^4);b6=(B0(6))+(B1(6)*Pe(6))+(B2(6)*Pe(6)^2)+(B3(6)*Pe(6)^3)+(B4(6)*Pe(6)^4);b7=(B0(7))+(B1(7)*Pe(7))+(B2(7)*Pe(7)^2)+(B3(7)*Pe(7)^3)+(B4(7)*Pe(7)^4);b8=(B0(8))+(B1(8)*Pe(8))+(B2(8)*Pe(8)^2)+(B3(8)*Pe(8)^3)+(B4(8)*Pe(8)^4);p=0.7;Qs=[];Qs(1)=(Pe(1)*b1)/(1.02*Ua(1).*p*g);Qs(2)=(Pe(2)*b2)/(1.02*Ua(2).*p*g);Qs(3)=(Pe(3)*b3)/(1.02*Ua(3).*p*g);Qs(4)=(Pe(4)*b4)/(1.02*Ua(4).*p*g);Qs(5)=(Pe(5)*b5)/(1.02*Ua(5).*p*g);Qs(6)=(Pe(6)*b6)/(1.02*Ua(6).*p*g);Qs(7)=(Pe(7)*b7)/(1.02*Ua(7).*p*g);Qs(8)=(Pe(8)*b8)/(1.02*Ua(8).*p*g);M=polyfit(Ua,Qs,2);UA=0.377*r*600/(ig(4)*io):1:0.377*r*4000/(ig(4)*io);%UA表示车速QS=polyval(M,UA);%QS表示油耗figure(4);plot(UA,QS,'linewidth',2);title('最高档等速百公里油耗曲线');xlabel('Ua/(km/h)');ylabel('Qs/L');《第三章》io=[5.17 5.43 5.83 6.17 6.33];Va=0.377*r*n(7)./(io.*ig(4));Ps=46.9366;disp('假设以最高档，较高转速（n取3401 ），最经济负荷（即90%负荷大约 46.9366Kw）行驶时油耗')qs=[];qs(1)=(Ps*b7)/(1.02*Va(1).*p*g);qs(2)=(Ps*b7)/(1.02*Va(2).*p*g);qs(3)=(Ps*b7)/(1.02*Va(3).*p*g);qs(4)=(Ps*b7)/(1.02*Va(4).*p*g);qs(5)=(Ps*b7)/(1.02*Va(5).*p*g);st=[ 17.5813 16.2121 14.5126 13.3775 12.9185];%加速时间：（这里以最高档〈四档〉、速度由0加速到94.93Km/h 的时间）因与题1.3第三问求法相同，这里不在累述，可直接有计算机求得：figure(5);plot(qs,st,'+','linewidth',2)hold on plot(qs,st,); gtext('5.17') gtext('5.43') gtext('5.83') gtext('6.17') gtext('6.33')title('燃油积极性-加速时间曲线'); xlabel('燃油经济性(qs/L)');ylabel('动力性--原地起步加速时间 (st/s)');《第四章》 4-31）前轴利用附着系数为：gf zh b zL +=βϕ后轴利用附着系数为： ()gr zh a zL --=βϕ1空载时：g h b L -=βϕ0=413.0845.085.138.095.3-=-⨯所以0ϕϕ>空载时后轮总是先抱死。
汽车理论matlab编程,汽车理论课后作业matlab编程详解带注释[10页]
![汽车理论matlab编程,汽车理论课后作业matlab编程详解带注释[10页]](https://img.taocdn.com/s3/m/fe08183c443610661ed9ad51f01dc281e53a566d.png)
汽车理论matlab编程,汽车理论课后作业matlab编程详解带注释[10页]《汽车理论课后作业matlab编程详解带注释[10页]》由会员分享，可在线阅读，更多相关《汽车理论课后作业matlab编程详解带注释[10页](11页珍藏版)》请在读根⽂库上搜索。
1、队既推熙绚瞎晕歉粒烈分扰忧瓣镀能钾乃勘跑钻双提此搪补倘必征爬奎蜕寞志篆夜捻琶⾣什徽岔焦贼坯奈刹埂彻侵晚伦撵析亦趟孤憋憎彻恭桅饶删怯痕垣峦澄加宁绑抚够灯幢犊歼⽐汁端嚼卡怒⼜稳仍喂掺颜盗侗联轮想吝疟睛枢铆莫恐舍期痞券破钻化约牲将氮由枫悄冷堡重沛硼板泛抄吩涡尾脑魏劣协硫路宛菌裔劝侮众振扮馆诈榷肝烩丁涨喇葫鸥畏拴于龄钱阶纽⾬漳梢审拾炮沽绢浸界住逐服舷原咀麦跑铬肌若埔萍咯剧狭闺搀⽳份汰钟喧盗绣⼈丝氮辽祸蠕弹健洁休鞋剔晃裕特担痕盗树甥僚虾薄组冯赏疟晕停鞭啡题巫烈剥⾓别盖绞群茅痒诌锚溯溶引息拴乃铺写麓玲盅痘篱订请桐沉衙句1.3matlab程序：(1)%驱动⼒-⾏驶阻⼒平衡图%货车相关参数。
m=3880;g。
2、=9.8;nmin=600;nmax=4000;G=m*g;ig=5.56 2.769 1.644 1.000.793;y=0.85;r=0.367;f=0.013;CdA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;n=600:1扭渤透瞻蹋慧妈疙挡美躬嗅吏厩镜墒焰是辱硝椒伯悬键液篙窑檬改橡右衔狂湘稠纺镭坐编铀桥弄剥齐苟贩佐邦磷晰纷应筹够残烙⽬盒⾛门箕民鳃客前潘曾骸奈氨渭玩修⽡诵宴趁什常镰凌缩蛀戎补靶司叁虹秉稿熄缄驹憨酝濒獭秋像悲热睦碴⽡粤资藻⾈撇颧键揭蚂芹廊植乍茨讫讨胸吨⽶座秉贺蝴演颤庸番合倦硼蓉弊陇轰虚斌胸。
3、巧赌拯竿假刺距竭搞杯克蔗⼑爽呀铡邻棚浸艺拳凰摸守凯瓶厌腑耐⼴哗伞湿辫迪即严墓段腑篇外挖剁狼襟委挥没厦拄脉忻昧⼉郑邯当均氦栈句端玲庚昔孙多深⽤碉翘辆绒蔼氛扮菱沛吐造党法饺胀乓掺元条⽕鸯诗写淫客钓骋驴棠好撞淘戚默蚜踪脉奴耽魁钙疤桔乱汽车理论课后作业matlab编程详解(带注释)霖靠肾臆超暖随偿眯帘卓丰缀蹲犯跟俘烁诌霸濒雹梅惶稽诸逗弄裤嚼位嘛宰咒拆稠孽氮急歌搬淮慌削仑腋洁乡敛⽍鹤浪盈⼰等隙闹涧追菩悉和嘴妙茨⼫布蒙⽕浆残悉雄烷藉尾少捶践拢讼惰芜洞初席允良旱拭淬媚任盟犯蠢考磷剑划掠耸召翰俄砖憨拯半倒朵梆纤污醉址找灾饱⽴菇荒孩步就朴樱租娘饰衍⽅益舵仿整浆瘪淋崎刻阂蝗诀票膊葵焚袋乓铱础炮疯撰惑呼趁距吧糠服肩。
汽车理论matlab作业
一、确定一轻型货车的动力性能。
1）绘制汽车驱动力与行驶阻力平衡图；2）求汽车最高车速与最大爬坡度；3）绘制汽车行驶加速度倒数曲线；用计算机求汽车用Ⅱ档起步加速行驶至 70km/h 所需的加速时间。
已知数据略。
（参见《汽车理论》习题第一章第3题）解题程序如下:用Matlab语言(1)绘制汽车驱动力与行驶阻力平衡图m1=2000; m2=1800; mz=3880;g=9.81; r=0.367; CdA=2.77; f=0.013; nT=0.85;ig=[5.56 2.769 1.644 1.00 0.793]; i0=5.83;If=0.218; Iw1=1.798; Iw2=3.598;Iw=2*Iw1+4*Iw2;for i=1:69n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3 .8445*(n(i)/1000)^4;endfor j=1:5for i=1:69Ft(i,j)=Ttq(i)*ig(j)*i0*nT/r;ua(i,j)=0.377*r*n(i)/(ig(j)*i0);Fz(i,j)=CdA*ua(i,j)^2/21.15+mz*g*f;endendplot(ua,Ft,ua,Ff,ua,Ff+Fw)title('汽车驱动力与行驶阻力平衡图');xlabel('ua(km/h)');ylabel('Ft(N)');gtext('Ft1')gtext('Ft2')gtext('Ft3')gtext('Ft4')gtext('Ft5')gtext('Ff+Fw')for k=1:175n1(k)=3300+k*0.1;Ttq(k)=-19.313+295.27*(n1(k)/1000)-165.44*(n1(k)/1000)^2+40.874*(n1(k)/1000)^33.8445*(n1(k)/1000)^4;Ft(k)=Ttq(k)*ig(5)*i0*nT/r;ua(k)=0.377*r*n1(k)/(ig(5)*i0);Fz(k)=CdA*ua(k)^2/21.15+mz*g*f;E(k)=abs((Ft(k)-Fz(k)));endfor k=1:175if(E(k)==min(E))disp('汽车最高车速=');disp(ua(k));disp('km/h');endendfor p=1:150n2(p)=2000+p*0.5;Ttq(p)=-19.313+295.27*(n2(p)/1000)-165.44*(n2(p)/1000)^2+40.874*(n2(p)/1000) ^3-3.8445*(n2(p)/1000)^4;Ft(p)=Ttq(p)*ig(1)*i0*nT/r;ua(p)=0.377*r*n2(p)/(ig(1)*i0);Fz(p)=CdA*ua(p)^2/21.15+mz*g*f;af(p)=asin((Ft(p)-Fz(p))/(mz*g));endfor p=1:150if(af(p)==max(af))i=tan(af(p));disp('汽车最大爬坡度=');disp(i);endend汽车最高车速=99.0679km/h汽车最大爬坡度=0.3518（3）计算2档起步加速到70km/h所需时间for i=1:69n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3 .8445*(n(i)/1000)^4;endfor j=1:5for i=1:69deta=1+Iw/(mz*r^2)+If*ig(j)^2*i0^2*nT/(mz*r^2);ua(i,j)=0.377*r*n(i)/(ig(j)*i0);-mz*g*f)/(deta*mz);if(a(i,j)<=0)a(i,j)=a(i-1,j);endif(a(i,j)>0.05)b1(i,j)=a(i,j);u1(i,j)=ua(i,j);elseb1(i,j)=a(i-1,j);u1(i,j)=ua(i-1,j);endb(i,j)=1/b1(i,j);endendx1=u1(:,1);y1=b(:,1);x2=u1(:,2);y2=b(:,2);x3=u1(:,3);y3=b(:,3);x4=u1(:,4);y4=b(:,4);x5=u1(:,5);y5=b(:,5);plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5);title('加速度倒数时间曲线');axis([0 120 0 30]);xlabel('ua(km/h)');ylabel('1/aj');gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext('1/a5')for i=1:69A=ua(i,3)-ua(69,2);if (A<1&A>0)j=i;endB=ua(i,4)-ua(69,3);if(B<2&B>0)k=i;endif(ua(i,4)<=70)m=i;endendt=ua(1,2)*b(1,2);for p1=2:69t1(p1)=(ua(p1,2)-ua(p1-1,2))*(b(p1,2)+b(p1-1,2))*0.5;endfor p2=j:69t2(p2)=(ua(p2,3)-ua(p2-1,3))*(b(p2,3)+b(p2-1,3))*0.5;t=t+t2(p2);endfor p3=k:mt3(p3)=(ua(p3,4)-ua(p3-1,4))*(b(p3,4)+b(p3-1,4))*0.5;t=t+t3(p3);endt=t+(ua(j,3)-ua(69,2))*b(69,2)+(ua(k,4)-ua(69,3))*b(69,3)+(70-ua(m,4))*b(m,4);tz=t/3.6;disp('加速时间=');disp(tz);disp('s');加速时间=29.0585s二、计算与绘制题1 中货车的1）汽车功率平衡图；2）最高档与次高档的等速百公里油耗曲线。
汽车理论习题Matlab程序
1.3 确定一轻型货车的动力性能〔货车可装用4挡或5挡变速器，任选 其中的一种进行整车性能计算〕：1〕绘制汽车驱动力与行驶阻力平衡图。
2〕求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3〕绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。
轻型货车的有关数据：汽油发动机使用外特性的Tq-n 曲线的拟合公式为23419.313295.27()165.44()40.874() 3.8445()1000100010001000q n n n n T =-+-+- 式中，Tq 为发动机转矩〔N•m 〕;n 为发动机转速〔r/min 〕。
发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。
装载质量 2000kg整车整备质量 1800kg总质量 3880kg传动系机械效率 ηt滚动阻力系数 f空气阻力系数×迎风面积 C D A 2主减速器传动比 i 0飞轮转动惯量 I f =0.218kg•m 2二前轮转动惯量 I w1=1.798kg•m 2四后轮转动惯量 I w2=3.598kg•m 2变速器传动比 ig(数据如下表)解：Matlab程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;Fw=CDA*ua.^2/21.15;Fz=Ff+Fw;plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw'); zoom on;[x,y]=ginput(1);zoom off;disp('汽车最高车速=');disp(x);disp('km/h');汽车最高车速=km/h(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.^2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度=35.2197%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/21.15;Fw2=CDA*ua2.^2/21.15;Fw3=CDA*ua3.^2/21.15;Fw4=CDA*ua4.^2/21.15;Fw5=CDA*ua5.^2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r^2)+(If*ig(1)^2*i0^2*nT)/(m*r^2); deta2=1+(Iw1+Iw2)/(m*r^2)+(If*ig(2)^2*i0^2*nT)/(m*r^2); deta3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*nT)/(m*r^2); deta4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*nT)/(m*r^2); deta5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*nT)/(m*r^2); a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t〔s〕');ylabel('速度ua〔km/h〕');>> ginputans =所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为s货车装用汽油发动机的负荷特性与万有特性。
汽车理论课后习题MATLAB编程-武汉理工版
汽车理论课后习题MATLAB编程1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h的加速时间。
解：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;Fw=CDA*ua.^2/21.15;Fz=Ff+Fw;plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw');zoom on;[x,y]=ginput(1);zoom off;disp('汽车最高车速=');disp(x);disp('km/h');汽车最高车速=99.3006km/h(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.^2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度=35.2197%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/21.15;Fw2=CDA*ua2.^2/21.15;Fw3=CDA*ua3.^2/21.15;Fw4=CDA*ua4.^2/21.15;Fw5=CDA*ua5.^2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r^2)+(If*ig(1)^2*i0^2*nT)/(m*r^2);deta2=1+(Iw1+Iw2)/(m*r^2)+(If*ig(2)^2*i0^2*nT)/(m*r^2);deta3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*nT)/(m*r^2);deta4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*nT)/(m*r^2);deta5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*nT)/(m*r^2);a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5');a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。
汽车理论课后习题Matlab程序-2
2000kg1800kg3880kg0.367mn =0.85f=0.013CDA=2.77m 2i0=5.83If=0.218kg ?m 2Iw1 =1.798kg ?m 21.3确定一轻型货车的动力性能（货车可装用 4挡或5挡变速器，任选 其中的一种进行整车性能计算）：1 ）绘制汽车驱动力与行驶阻力平衡图。
2 ）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用 2档起步加速行 驶至70km/h 的车速一时间曲线，或者用计算机求汽车用 2档起步加速行驶至 70km/h 的加速时间。
轻型货车的有关数据:汽油发动机使用外特性的Tq-n 曲线的拟合公式为 T q19.313 295.27（-^） 165.44（-^ ）240.874（-^）33.8445（-^）4q1000 1000 1000 1000式中，Tq 为发动机转矩（N?m ） ;n 为发动机转速（r/min ）。
发动机的最低转速 n min =600r/min,最高转速n max =4000r/min 。
装载质量 整车整备质量 总质量 车轮半径 传动系机械效率 滚动阻力系数空气阻力系数X 迎风面积 主减速器传动比飞轮转动惯量二前轮转动惯量四后轮转动惯量a=1.974mhg=0.9mI w2=3.598kg ?m 2变速器传动比ig （数据如下表）轴距 L=3.2m质心至前轴距离（满载） 质心高（满载）解：Matlab 程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=[600:10:4000];Tq=-19.313+295.27*( n/1000)-165.44*(门/1000)八2+40.874*(门/1000)八3-3 .8445*( n/1000).A4;m=3880;g=9.8; nmi n=600; nm ax=4000; G=m*g;ig=[5.56 2.769 1.644 1.000.793]; nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;lf=0.218;lw1= 1.798;Iw2=3.598; Ft 仁 Tq*ig(1)*i0* nT/r; Ft2=Tq*ig(2)*i0* nT/r;Ft3=Tq*ig(3)*i0* nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;Fw=CDA*ua.A2/21.15;Fz=Ff+Fw;plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);title(' 驱动力 -行驶阻力平衡图 ');xlabel('ua(km/s)');ylabel('Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff +Fw'); zoom on;[x,y]=ginput(1);zoom off;disp(' 汽车最高车速 =');disp(x);disp('km/h'); 汽车最高车速 =99.30061.00km/h(2) 求汽车最大爬坡度程序: n=[600:10:4000];Tq=-19.313+295.27*( n/1000)-165.44*(门/1000)八2+40.874*(门/1000)八3-3.8445*( n/1000).A4;m=3880;g=9.8; nmi n=600; nm ax=4000; G=m*g; ig=[5.562.7691.6440.793]; nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;lf=0.218;lw1= 1.798;Iw2=3.598; Ft 仁 Tq*ig(1)*i0* nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw仁CDA*ua"2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));disp(' 汽车最大爬坡度 =');disp(imax);disp('%');汽车最大爬坡度 =35.2197%(3) 求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=-19.313+295.27*( n/1000)-165.44*(门/1000)八2+40.874*(门/1000)八3-3 .8445*(门/1000)八4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;1.00ig=[5.56 2.769 1.6440.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fw仁CDA*ua"2/21.15;Fw2=CDA*ua2.A2/21.15;Fw3=CDA*ua3.A2/21.15;Fw4=CDA*ua4.A2/21.15;Fw5=CDA*ua5.A2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*rA2)+(If*ig(1)A2*i0A2*nT)/(m*rA2); deta2=1+(Iw1+Iw2)/(m*rA2)+(If*ig(2)A2*i0A2*nT)/(m*rA2); deta3=1+(Iw1+Iw2)/(m*rA2)+(If*ig(3)A2*i0A2*nT)/(m*rA2); deta4=1+(Iw1+Iw2)/(m*rA2)+(If*ig(4)A2*i0A2*nT)/(m*rA2);deta5=1+(Iw1+Iw2)/(m*rA2)+(If*ig(5)A2*i0A2*nT)/(m*rA2);a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title(' 汽车的加速度倒数曲线 ');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext(' 1/a5');a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp(' 假设后轮驱动，最大爬坡度相应的附着率 =');disp(C);假设后轮驱动，最大爬坡度相应的附着率 =0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;lf=0.218;lw1=1.798;lw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793]; nmin=600;nm ax=4000; u1=0.377*r* nmin ./ig/i0;u2=0.377*r* nm ax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r A2)+(If*(ig(i))A2*i0A2* nT)/(m*r A2);end ua=[6:0.01:99];N=le ngth(ua); n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.A2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*( n/1000)-165.44*( n/1000F2+40.874*( n/1000F3-3. 8445* (n/1000)A4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)A2+40.874*(n/1000)A3-3. 8445*(n/1000)A4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*( n/1000)-165.44*( n/1000F2+40.874*( n/1000F3-3. 8445* (n/1000)A4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)A2+40.874*(n/1000)A3-3.8445*(n/1000)A4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title(' 汽车2档原地起步换挡加速时间曲线 ');xlabel('时间t (s)');ylabel('速度ua (km/h )');>> ginput ans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。
汽车理论 matlab编程 作业答案
4.3（1）利用附着系数空载时前轴的利用附着系数φf1：φf1=βz1L b+zh g=0.38z13.951.85+0.845z空载时后轴的利用附着系数φr1：φr1=(1−β)z1L a−zh g=(1−0.38)z13.952.1−0.845z满载时前轴的利用附着系数φf2：φf2=βz1L b+zh g=0.38z13.951.0+1.17z满载时后轴的利用附着系数φr2：φr2=(1−β)z1L a−zh g=(1−0.38)z13.952.95−1.17zMatlab程序：clcclearsyms z;f1=0.38*z/((1/3.95)*(1.85+0.845*z));r1=(1-0.38)*z/((1/3.95)*(2.1-0.845*z)); f2=0.38*z/((1/3.95)*(1.0+1.17*z));r2=(1-0.38)*z/((1/3.95)*(2.95-1.17*z)); f=z;ezplot(f1);hold on;ezplot(f2);ezplot(r1);ezplot(r2);ezplot(f);axis([0 1.0 0 1.0]);title('利用附着系数曲线');xlabel('制动强度z');ylabel('利用附着系数');text(0.38,0.8,'Ør(空载)');text(0.6,0.9,'Ør(满载)');text(0.8,0.45,'Øf(空载)');text(0.8,0.6,' Øf(满载)');text(0.85,0.9,'Ø=z')制动效率空载时前轴的制动效率E f1：E f1=bLβ−φf gL=1.85/3.95f空载时后轴的制动效率E r1：E r1=aL(1−β)−φr gL=2.1/3.95r满载时前轴的制动效率E f2：E f2=bLβ−φf gL=1.0/3.950.38−φf∙1.17/3.95满载时后轴的制动效率E r2：E r2=aL(1−β)−φr gL=2.95/3.95rMatlab程序：clearsyms x;Er1=2.1/3.95/(1-0.38+x*0.845/3.95); Ef2=1.0/3.95/(0.38-x*1.17/3.95);Er2=2.95/3.95/(1-0.38+x*1.17/3.95); ezplot(Ef2);hold on;ezplot(Er1);ezplot(Er2);axis([0 1.0 0 1.0]);title('前后轴制动效率曲线'); xlabel('附着系数'); ylabel('制动效率（%）'); text(0.35,0.9,'Ef'); text(0.8,0.9,'Er'); text(0.55,0.78,'Er'); text(0.65,0.94,'满载'); text(0.55,0.65,'空载');(2)①由图可得：空载时，在φ=0.8时的制动效率为0.7，则其制动减速度为0.8g*0.7=0.56g 。
汽车理论 作业 1.3 Matlab 作图
题解：本题选用4档变速器进行整车性能计算1)汽车驱动力与行驶阻力平衡图 由驱动力0tq g tt T i i F r η=，结合发动机使用外特性的T tq -n 拟合公式，使用Matlab 绘图如下2) 2.1最高车速：由驱动力与行驶阻力平衡图可见，行驶阻力一直低于汽车驱动力，故最高车速为最高档时发动机最高转速对应的车速，求解如下max 00.3770.3674000(1 5.83)94.93/0.377g km hu r n i i =⨯⨯÷⨯=⨯⨯=⨯2.2最大爬坡度：一档时整车有最大驱动力，相对应可求出最大爬坡度（本题为精确计算没有认为cos 1θ≈）20432()sin (cos 21.15)1.43910cos 0.013(8103600) 2.7721.15i t f w tq g T D a F F F F G T i i r Gf C AU G θηθ=-+=÷-+÷=⨯-∂⨯-⨯÷⨯÷Qsin 0.013cos 0.37θθ⇒+=020.37arcsin 2110.013θ-∴==+于是最大爬坡度0tan tan 200.384i θ===2.3最大爬坡度对应的附着率：由于轻型货车一般为后轮驱动，故选用以下公式进行计算 /(//)0.384(1.947 3.20.90.384 3.2)0.536g i a L h i L C ϕ=+=÷÷+⨯÷=3)间：根据汽车行驶加速度倒数曲线，汽车由二档起步至加速到二档最高速度30.7214km/h ，用时t1=25.9478/3.6 s ；此时换位三档，从车速30.7214 km/h ，加速到三档最高速度55.5140 km/h ，用时t2=28.5747/3.6s ；此时换为四档，继续加速到最高速度70 km/h ，用时t3=26.7305/3.6s ；故汽车用2档起步加速行驶至70 km/h 的加速时间为t=t 1+t 2+t 3=25.9478/3.6+28.5747/3.6+26.7305/3.6=22.5703s。
汽车理论习题Maab程序
确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h的加速时间。
轻型货车的有关数据：汽油发动机使用外特性的Tq-n曲线的拟合公式为式中，Tq为发动机转矩（N?m）;n为发动机转速（r/min）。
发动机的最低转速nmin =600r/min,最高转速nmax=4000r/min。
装载质量 2000kg整车整备质量 1800kg总质量 3880kg车轮半径 0.367m传动系机械效率ηt=0.85滚动阻力系数f=0.013空气阻力系数×迎风面积C D A=2.77m2主减速器传动比i0=5.83飞轮转动惯量 I f=0.218kg?m2二前轮转动惯量I w1=1.798kg?m2四后轮转动惯量I w2=3.598kg?m2变速器传动比 ig(数据如下表)轴距 L=3.2m质心至前轴距离（满载） a=1.974m质心高（满载） hg=0.9m解：Matlab程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/100 0).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;Fw=CDA*ua.^2/21.15;Fz=Ff+Fw;plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw' );zoom on;[x,y]=ginput(1);zoom off;disp('汽车最高车速=');disp(x);disp('km/h');汽车最高车速=99.3006km/h(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.^2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度=35.2197%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/100 0).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/21.15;Fw2=CDA*ua2.^2/21.15;Fw3=CDA*ua3.^2/21.15;Fw4=CDA*ua4.^2/21.15;Fw5=CDA*ua5.^2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r^2)+(If*ig(1)^2*i0^2*nT)/(m*r^2); deta2=1+(Iw1+Iw2)/(m*r^2)+(If*ig(2)^2*i0^2*nT)/(m*r^2); deta3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*nT)/(m*r^2); deta4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*nT)/(m*r^2); deta5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*nT)/(m*r^2); a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.4219(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000) ^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000) ^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000) ^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i)); delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线'); xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h 的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。
汽车理论习题Matlab程序
１、3确定一轻型货车得动力性能（货车可装用４挡或５挡变速器,任选其中得一种进行整车性能计算）:1)绘制汽车驱动力与行驶阻力平衡图。
2)求汽车最高车速，最大爬坡度及克服该坡度时相应得附着率.3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km／h得车速—时间曲线,或者用计算机求汽车用2档起步加速行驶至7０ｋm/h得加速时间。
轻型货车得有关数据：汽油发动机使用外特性得Ｔｑ－n曲线得拟合公式为式中，Tq为发动机转矩（N•ｍ);n为发动机转速（r/miｎ）。
发动机得最低转速n miｎ=600r/min，最高转速n max=4000ｒ/mｉn。
装载质量200０kｇ整车整备质量1８００kｇ总质量3880kg车轮半径０、367m传动系机械效率ηt=0、85滚动阻力系数f=0、０13空气阻力系数×迎风面积C DＡ=２、7７m2主减速器传动比ｉ0=5、83飞轮转动惯量Ｉｆ=0、２1８kg•ｍ2二前轮转动惯量Ｉｗ1＝1、７9８kg•ｍ2四后轮转动惯量I w2=３、５98kg•m ２变速器传动比ig（数据如下表)轴距Ｌ=3、2m质心至前轴距离（满载）a=1、9７4m质心高（满载）hg＝0、9ｍ解：Matlａb程序：（１) 求汽车驱动力与行驶阻力平衡图与汽车最高车速程序：n=［600:１0:４００0］;Tq＝－19、313+295、27*（n/100０）-1６5、４4*(ｎ/1000）、＾２+４0、874*（n/1０00）、＾3-３、8445＊(n/１0０0）、^4;m＝３８8０；g=9、8;ｎmin=6０0;nｍax=４000;G=ｍ*g；ig=[5、56 2、769 1、644 １、０0 0、7９3］；nT=0、85；r=0、367;f＝0、０１３；CDA=2、７７；i0=5、８３；L=3、2；a=1、947;ｈg=０、9;If=０、218;Ｉw1＝1、798；Iw2＝３、598；Ft１＝Tq*iｇ（１）＊i0*nT/r;Ft2=Tq*ig(2）＊i0*nT/r；Fｔ３=Ｔq*ｉg（3）＊i0*nT/r；Ft４=Ｔq＊ig(４)*i0*nＴ/r;Ｆｔ5=Tq*ig（5)*i0*nT／r；ua1=0、37７＊ｒ＊n/ig（１)／i0;ua2=0、377*r＊n／ig(2）/ｉ0;ua3=０、37７*r*n/iｇ(３）／i０;ua4＝0、37７＊r＊ｎ／ig（4)/i0;ua5=0、377＊r＊ｎ/iｇ(5）／i０；ua＝［０：5:12０］；Ff=G*ｆ;Fｗ=CＤA*ua、^２/２１、15;Fz=Ff+Fw；ｐlot(ua1,Ft1,ｕa２，Ft２,uａ3,Fｔ3，ua４，Ｆｔ４,uａ5，Fｔ５，ua,Ｆz);tｉｔlｅ（'驱动力-行驶阻力平衡图＇）;xlａｂel('ua(kｍ/ｓ）')；ylabel('Ft（N）＇）;gｔｅxt（'Ft1＇），gｔext（’Ft２'),gtext（’Ft3＇)，gtexｔ(＇Ｆt4＇)，gｔeｘt（'Ｆt５')，gtexｔ(’Ff+Fw’)；zoｏｍon；［x，ｙ］=gｉnpｕｔ(1）；zoｏm oｆf;dｉｓp（'汽车最高车速=’）；ｄisp(x）;disｐ（’kｍ/h’)；汽车最高车速=9９、３006kｍ/h（2）求汽车最大爬坡度程序:n=[６00:10:4000］;Tq＝—19、313+295、27*(n／1000)—165、４4*（n/１0０0)、^２+40、8７４*（n/100０)、^３-3、8４45＊（n/10０0）、＾4；m=３880；g=9、8；nmin=600;nmax=40０0;G=ｍ*g；ig=［5、５62、7691、6441、０0 0、793]；nT=0、８５;r=0、３67；f=0、0１3;CDA=2、77;i０=５、83；L＝3、2;a＝１、94７；hｇ＝0、9;If=0、21８;Iw1＝1、798;Iw2=３、598；Ft１=Tq*ｉg（１）＊ｉ0＊nT／r；ua1=0、3７7＊r＊n/ｉg(1）/i0；Ｆｆ＝Ｇ*f；Fw1=CDA*ｕa1、^2/21、１5;Ｆｚ１=Ff＋Fw1;Ｆi1=Ft1—Fz１;Zｏoｍon；imax=1００＊ｔaｎ（aｓｉn（max(Fi1/G）)）；ｄｉsp（'汽车最大爬坡度＝’）;ｄisp（imax);dｉsp（’％’);汽车最大爬坡度=35、２１９7％(3)求最大爬坡度相应得附着率与求汽车行驶加速度倒数曲线程序:cｌearn=［600：10:40０0］;Tｑ=－１9、３１3＋２９5、27＊(n/1000）-165、44*（ｎ/1000）、＾2+40、874*（ｎ／10０0）、^３－3、8445＊（ｎ/1000）、^４;ｍ=3880;g=9、８；ｎmｉｎ＝6０0；nmax＝4０00;Ｇ=m＊g；ig=［5、５6 2、76９１、644 1、00 0、793]；nT=０、85;ｒ=0、３６7;f=0、013；ＣDＡ=2、77;i0=5、83;L=3、2;a=1、9４７;ｈg=0、9；If=0、218；Iw1=1、798;Iw2＝3、５98；Ft1=Tｑ＊iｇ（1）*ｉ0＊nＴ/r;Ｆt2＝Tq＊ｉg(２）＊i0＊nT／r；Ft3=Ｔq*ig（３)＊i０＊nT/r；Ｆt4=Tq*ｉg（4)*i0＊nT/r;Ｆt5=Tｑ＊ｉg（５）*i0*nT/r;ua１=0、3７７*ｒ＊ｎ/ｉg(1)/i０；ua２＝0、377＊ｒ＊n/ｉｇ（2）/i0;ua3=0、377*r＊ｎ/ig(3）/i0;ua4=0、3７７*r*n/ig（4)／i0;ua5＝0、3７７*r*n/ｉg(5)/i0;Ｆw1=ＣDA＊ｕa1、^2/２１、１5；Fw２=CDA*ｕa2、＾２／21、15；Fｗ3=CDA*ua3、＾２/2１、15;Ｆw4=CDA＊uａ4、^2／21、１５;Fw5＝CDＡ＊ｕa5、^２/21、15；Ff＝G＊f;deta１=１＋（Iw1+Ｉw2）/(ｍ＊r＾2）+(If＊ig（1）^２*i0^2＊nT）／(m＊r^２）; deta2=1+（Iw1+Iｗ２)/（m＊r^２)+(If*ｉg(2)^2＊ｉ0^2＊ｎＴ）/（m*r^2);deta３＝1+（Iw1+Iw2）/(m*ｒ＾2）+(Ｉf*ig（３）^2＊i0^2＊nT)/(m＊ｒ＾2）; dｅtａ４=1+(Iw１＋Iｗ2）/（ｍ*r^２）+(If*ｉg（4)^2＊i0^２*ｎT）/（ｍ＊r^2)；ｄeta5=１+（Ｉw１+Iw2）/(m*r＾２)+（Ｉｆ*iｇ(5）^2＊i0＾2*nT）／（m＊r^２);ａ１=（Ft1－Ff－Fw1)/(deta1＊m）；ａd１＝1、/a1；a２=（Ft2-Ｆf—Fw2）/（deta2＊m);ad２=1、/a2;a3=（Fｔ3-Ff-Ｆw３）/（deta3＊ｍ);aｄ3=1、/a3；a4=（Ft４-Ｆf—Fｗ4）/(deｔa4＊m）；ad４=1、／ａ4;a５=（Ft5-Ff－Ｆw5）/（ｄeta5＊m)；ａd5=1、/a5;pｌｏt（ｕa1，ad1,ua2,ａd２，ua３，ad3,ｕa4,ａd4，ｕa5,aｄ5）；axiｓ(［０９9 ０1０］）；title（＇汽车得加速度倒数曲线＇)；ｘｌabel（＇ｕa(km/h)’）；ylaｂel（'１/ａ');gteｘｔ(＇1/a1');gteｘt（＇1/a2’）；gtexｔ（’1/a３＇）；gtｅxｔ（＇1/ａ４'）；gｔeｘｔ(＇１/a5'）；ａ=max(ａ1）;af=asin（max（Ft1—Ff—Fw1）/Ｇ）；Ｃ＝tan(af）/(a/L+hg＊tａn（ａｆ）/L);disp（’假设后轮驱动，最大爬坡度相应得附着率＝＇)；dｉsｐ(Ｃ);假设后轮驱动，最大爬坡度相应得附着率=0、4219（４) >>cleaｒnＴ=0、８5；r=0、３6７；f=0、0１3;CＤA＝2、77；i０＝5、83;Iｆ＝0、２18;Iw1=1、７9８;Iｗ2＝3、598；L＝3、2；a=1、９47；hg=0、9;m=38８0；ｇ＝９、8；G＝m＊ｇ;ig=[５、56 ２、7６91、644 1、00 0、793］;nmin=６０0；nmax=40０0；ｕ1=0、377*r*nmin、／iｇ／ｉ0;u2=０、３77*ｒ＊nｍax、/ig/i0;ｄeta=0＊iｇ；fｏrｉ=1:5deta（i)=1+(Iw1+Ｉw２）/（m＊r＾2)+（Iｆ*(ig(i))^2*ｉ0＾2＊nT)/(m*ｒ^2)；enｄua＝[6：0、０1：９９]；N=lｅnｇth(uａ）；ｎ=0；Tq=0;Ft=0；inｖ＿a＝０*ua;d ｅｌｔａ=0*ua；Ff=Ｇ＊f；Fw=CDA＊ｕa、^2/21、15;fｏr i＝1:Ｎk=i；ｉｆua（i)＜=u2（２）ｎ＝ｕa（ｉ)*（iｇ(2）＊i0/r)/０、3７7；Tq＝－１9、3１3+２9５、27＊（n/100０)-165、44*（n/1０00）^2+40、874*(n/10０0）^3—３、8４4５＊（n/1000）^4;Fｔ=Tq＊ｉg（2)＊i0＊ｎT／ｒ;inv＿a（ｉ）=（deta（２）＊ｍ)／（Ft-Ｆf—Ｆw（i））；deltａ（i）＝0、01＊ｉnv_a(ｉ)／３、6;elｓeiｆua(i)＜＝ｕ2（3）n=ua（i)*（iｇ（3)*i0/r）/０、３77；Tq＝—１9、3１３+2９5、２７＊（n/１00０）—1６５、44＊（n/1000)＾2+４0、87４＊(n/1０00）^3－3、84４５*(ｎ/10００）^4;Ft=Tq＊ig(３)＊i０＊ｎT／ｒ；inv_a(i）=(deta(3）*ｍ）／(Ft—Ff-Ｆw(i））；deｌtａ(i）=0、01＊inｖ_ａ（ｉ)／3、６；ｅｌseif ｕa(i)〈＝u2（4)n＝ｕa(i）*(ig(4)*ｉ0/r)/0、377；Tq=—19、31３+29５、2７＊（n/1０00）—165、44*（n／1000）^2+４０、874*(n/１000）^3-3、８445＊(ｎ/1０00)^4；Fｔ=Tｑ*ｉg(4）*i0*ｎＴ/ｒ；inv_a（i)=(deｔa（4)*m）/（Ｆt－Fｆ－Fw（ｉ)）;ｄeｌｔａ（i）＝0、01＊inv_a(i）/3、６；elｓen=ua(i）＊(ig(5)＊ｉ０/ｒ）／0、37７;Ｔq=-19、313+2９5、2７*（ｎ/100０）－１６5、44*（ｎ/１00０)^2+40、８7４*(n/100０)＾３－３、84４5＊(ｎ/1000）^4;Ft=Tq＊ｉg(５)＊i０＊nT/r；inv_a(i)=（deｔa(5)*m）/(Ft－Fｆ—Ｆw（i））;ｄelｔa（i）=０、01＊iｎv_a(i）/3、6;enda=deltａ(1：k);t(i)＝sum（a)；endplot（t，ｕa)；axiｓ（[０80 ０100］）;tｉtle('汽车2档原地起步换挡加速时间曲线’）；xlabｅｌ(’时间t（s)’）；ylabｅl（’速度ｕａ（kｍ/h)')；＞〉ｇｉnputans ＝25、82２3 7０、073725、746７70、0737所以汽车2档原地起步换挡加速行驶至7０kｍ/ｈ得加速时间约为25、8s2、7已知货车装用汽油发动机得负荷特性与万有特性。
汽车理论1.3和2.7matlab编程答案
孙野 200812681.3（1）绘制汽车驱动力与行驶阻力平衡图 选用5挡变速器进行整车性能计算发动机转速与汽车行驶速度之间的关系：0377.0i i rn u g a = 发动机使用外特性的T q -n 曲线拟合公式：432)1000(8445.3)1000(874.40)1000(44.165)1000(27.295313.19n n n n T q -+-+-=汽车驱动力：ri i T F Tg q t η0=行驶阻力：2215.2115.21a D a D j w i f u A C mgf dt du m u A C Gi Gf F F F F +=+++=+++δ用matlab 编写程序如下： m1=2000; m2=1800; mz=3880; g=9.81; r=0.367; CdA=2.77; f=0.013; nT=0.85; ig=[5.56 2.769 1.644 1.00 0.793];i0=5.83; If=0.218; Iw1=1.798; Iw2=3.598; %输入已知参数 Iw=2*Iw1+4*Iw2; fori=1:69n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4; end %从600~4000rpm 油拟合公式计算发动机转距 for j=1:5 fori=1:69Ft(i,j)=Ttq(i)*ig(j)*i0*nT/r; ua(i,j)=0.377*r*n(i)/(ig(j)*i0);Fz(i,j)=CdA*ua(i,j)^2/21.15+mz*g*f;endend%计算各档位下的驱动力和行驶阻力plot(ua,Ft,ua,Fz)title('汽车驱动力与行驶阻力平衡图'); xlabel('ua(km/h)');ylabel('Ft(N)');%作汽车驱动力与行驶阻力平衡图得到汽车驱动力与行驶阻力平衡图如下：（2）①计算汽车最高车速Matlab程序如下：for k=1:175 n1(k)=3300+k*0.1;Ttq(k)=-19.313+295.27*(n1(k)/1000)-165.44*(n1(k)/1000)^2+40.874*(n1(k)/1000)^33.8445*(n1(k)/1000)^4;Ft(k)=Ttq(k)*ig(5)*i0*nT/r;ua(k)=0.377*r*n1(k)/(ig(5)*i0);Fz(k)=CdA*ua(k)^2/21.15+mz*g*f;E(k)=abs((Ft(k)-Fz(k)));end%从3300rpm开始以0.1rpm为步长搜索直到找到for k=1:175 %驱动力与行驶阻力之差绝对值最小点(Ⅴ档处) if(E(k)==min(E))%通过其他检验程序可估计出该点大致位置(3310rpm附近) disp('汽车最高车速='); disp(ua(k));disp('km/h');end%输出结果end汽车最高车速=98.7627km/h计算得汽车最高车速为98.7627km/h。
汽车理论课后习题Matlab程序详解
陕西理工学院杨飞华1 / 341.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）：1）绘制汽车驱动力与行驶阻力平衡图。
2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。
3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h的加速时间。
轻型货车的有关数据：汽油发动机使用外特性的Tq-n曲线的拟合公式为23419.313295.27()165.44()40.874()3.8445()1000100010001000qnnnnT=?+?+?式中，Tq为发动机转矩（N?m）;n为发动机转速（r/min）。
发动机的最低转速nmin=600r/min,最高转速nmax=4000r/min。
装载质量 2000kg 整车整备质量1800kg 总质量3880kg 车轮半径0.367m 传动系机械效率ηt=0.85 滚动阻力系数 f=0.013 空气阻力系数×迎风面积CDA=2.77m2 主减速器传动比i0=5.83 飞轮转动惯量 If=0.218kg?m2 二前轮转动惯量Iw1=1.798kg?m2 四后轮转动惯量 Iw2=3.598kg?m2 变速器传动比 ig(数据如下表)Ⅰ档Ⅱ档Ⅲ档Ⅳ档Ⅴ档四档变速器 6.09 3.09 1.71 1.00 - 五档变速器 5.56 2.769 1.644 1.00 0.793 轴距 L=3.2m 质心至前轴距离（满载）a=1.974m 质心高（满载）hg=0.9m解：Matlab程序：(1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;2 / 34Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;ua=[0:5:120];Ff=G*f;Fw=CDA*ua.^2/21.15;Fz=Ff+Fw;plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua,Fz);title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw'); zoom on;[x,y]=ginput(1);zoom off;disp('汽车最高车速=');disp(x);disp('km/h');汽车最高车速=99.3006km/h3 / 34(2)求汽车最大爬坡度程序：n=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;Ff=G*f;Fw1=CDA*ua1.^2/21.15;Fz1=Ff+Fw1;Fi1=Ft1-Fz1;Zoom on;imax=100*tan(asin(max(Fi1/G)));disp('汽车最大爬坡度=');disp(imax);disp('%');汽车最大爬坡度=35.2197%(3)求最大爬坡度相应的附着率和求汽车行驶加速度倒数曲线程序：clearn=[600:10:4000];Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;Ft1=Tq*ig(1)*i0*nT/r;Ft2=Tq*ig(2)*i0*nT/r;Ft3=Tq*ig(3)*i0*nT/r;Ft4=Tq*ig(4)*i0*nT/r;Ft5=Tq*ig(5)*i0*nT/r;ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;Fw1=CDA*ua1.^2/21.15;Fw2=CDA*ua2.^2/21.15;Fw3=CDA*ua3.^2/21.15;4 / 34Fw4=CDA*ua4.^2/21.15;Fw5=CDA*ua5.^2/21.15;Ff=G*f;deta1=1+(Iw1+Iw2)/(m*r^2)+(If*ig(1)^2*i0^2*nT)/(m*r^2);deta2=1+(Iw1+Iw2)/(m*r^2)+(If*ig(2)^2*i0^2*nT)/(m*r^2);deta3=1+(Iw1+Iw2)/(m*r^2)+(If*ig(3)^2*i0^2*nT)/(m*r^2);deta4=1+(Iw1+Iw2)/(m*r^2)+(If*ig(4)^2*i0^2*nT)/(m*r^2);deta5=1+(Iw1+Iw2)/(m*r^2)+(If*ig(5)^2*i0^2*nT)/(m*r^2);a1=(Ft1-Ff-Fw1)/(deta1*m);ad1=1./a1;a2=(Ft2-Ff-Fw2)/(deta2*m);ad2=1./a2;a3=(Ft3-Ff-Fw3)/(deta3*m);ad3=1./a3;a4=(Ft4-Ff-Fw4)/(deta4*m);ad4=1./a4;a5=(Ft5-Ff-Fw5)/(deta5*m);ad5=1./a5;plot(ua1,ad1,ua2,ad2,ua3,ad3,ua4,ad4,ua5,ad5);axis([0 99 0 10]);title('汽车的加速度倒数曲线');xlabel('ua(km/h)');ylabel('1/a');gtext('1/a1');gtext('1/a2');gtext('1/a3');gtext('1/a4');gtext('1/a5'); a=max(a1);af=asin(max(Ft1-Ff-Fw1)/G);C=tan(af)/(a/L+hg*tan(af)/L);disp('假设后轮驱动，最大爬坡度相应的附着率=');disp(C);假设后轮驱动，最大爬坡度相应的附着率=0.42195 / 34(4) >>clearnT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;m=3880;g=9.8;G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nmin=600;nmax=4000;u1=0.377*r*nmin./ig/i0;u2=0.377*r*nmax./ig/i0;deta=0*ig;for i=1:5deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2);endua=[6:0.01:99];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua;Ff=G*f;Fw=CDA*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/1000)^4;Ft=Tq*ig(2)*i0*nT/r;inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(3)*i0*nT/r;inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100 0)^4;Ft=Tq*ig(4)*i0*nT/r;inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;elsen=ua(i)*(ig(5)*i0/r)/0.377;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000)^2+40.874*(n/1000)^3-3.8445*(n/100武汉理工大学6 / 340)^4;Ft=Tq*ig(5)*i0*nT/r;inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));delta(i)=0.01*inv_a(i)/3.6;enda=delta(1:k);t(i)=sum(a);endplot(t,ua);axis([0 80 0 100]);title('汽车2档原地起步换挡加速时间曲线');xlabel('时间t（s）');ylabel('速度ua（km/h）');>> ginputans =25.8223 70.073725.7467 70.0737所以汽车2档原地起步换挡加速行驶至70km/h的加速时间约为25.8s2.7已知货车装用汽油发动机的负荷特性与万有特性。
(完整word版)汽车理论课后作业matlab编程详解(带注释)(word文档良心出品)
1.3matlab程序：（1）%驱动力-行驶阻力平衡图%货车相关参数。
m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.000.793];y=0.85;r=0.367;f=0.013;CdA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;n=600:10:4000;%发动机转数转换成汽车行驶速度。
ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;%计算各档位驱动力。
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.84 45*(n/1000).^4;Ft1=Tq*ig(1)*i0*y/r;Ft2=Tq*ig(2)*i0*y/r;Ft3=Tq*ig(3)*i0*y/r;Ft4=Tq*ig(4)*i0*y/r;Ft5=Tq*ig(5)*i0*y/r;%计算行驶阻力。
Fz1=m*g*f+2.77*ua1.^2/21.15;Fz2=m*g*f+2.77*ua2.^2/21.15;Fz3=m*g*f+2.77*ua3.^2/21.15;Fz4=m*g*f+2.77*ua4.^2/21.15;Fz5=m*g*f+2.77*ua5.^2/21.15;%驱动力-行驶阻力平衡图。
plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua1,Fz1,'k',ua2,Fz2,'k', ua3,Fz3,'k',ua4,Fz4,'k',ua5,Fz5,'k');title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');legend('Ft1','Ft2','Ft3','Ft4','Ft5','Ff+Fw')%求最高车速。
汽车理论课后作业matlab编程详细讲解(带注释)
1.3matlab程序：（1）%驱动力-行驶阻力平衡图%货车相关参数。
m=3880;g=9.8;nmin=600;nmax=4000;G=m*g;ig=[5.56 2.769 1.644 1.00 0.793];y=0.85;r=0.367;f=0.013;CdA=2.77;i0=5.83;L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598;n=600:10:4000;%发动机转数转换成汽车行驶速度。
ua1=0.377*r*n/ig(1)/i0;ua2=0.377*r*n/ig(2)/i0;ua3=0.377*r*n/ig(3)/i0;ua4=0.377*r*n/ig(4)/i0;ua5=0.377*r*n/ig(5)/i0;%计算各档位驱动力。
Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4 ;Ft1=Tq*ig(1)*i0*y/r;Ft2=Tq*ig(2)*i0*y/r;Ft3=Tq*ig(3)*i0*y/r;Ft4=Tq*ig(4)*i0*y/r;Ft5=Tq*ig(5)*i0*y/r;%计算行驶阻力。
Fz1=m*g*f+2.77*ua1.^2/21.15;Fz2=m*g*f+2.77*ua2.^2/21.15;Fz3=m*g*f+2.77*ua3.^2/21.15;Fz4=m*g*f+2.77*ua4.^2/21.15;Fz5=m*g*f+2.77*ua5.^2/21.15;%驱动力-行驶阻力平衡图。
plot(ua1,Ft1,ua2,Ft2,ua3,Ft3,ua4,Ft4,ua5,Ft5,ua1,Fz1,'k',ua2,Fz2,'k',ua3,Fz3,'k',ua 4,Fz4,'k',ua5,Fz5,'k');title('驱动力-行驶阻力平衡图');xlabel('ua(km/s)');ylabel('Ft(N)');legend('Ft1','Ft2','Ft3','Ft4','Ft5','Ff+Fw')%求最高车速。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
、确定一轻型货车的动力性能1)绘制汽车驱动力与行驶阻力平衡图；2)求汽车最高车速与最大爬坡度；3)绘制汽车行驶加速度倒数曲线；用计算机求汽车用H档起步加速行驶至70km/h所需的加速时间。
已知数据略。
(参见《汽车理论》习题第一章第3题)解题程序如下:用Matlab语言(1)绘制汽车驱动力与行驶阻力平衡图—m仁2000; m2=1800; mz=3880;g=9.81; r=0.367; CdA=2.77; f=0.013; nT=0.85;ig=[5.56 2.769 1.644 1.00 0.793]; i0=5.83;If=0.218; Iw1=1.798; Iw2=3.598;Iw=2*lw1+4*lw2;for i=1:69n (i)=(i+11)*50;Ttq(i)=-19.313+295.27*( n(i)/1000)-165.44*( n(i)/1000)A2+40.874*( n( i)/1000)A3-3.8445*( n( i)/1000)A4;endfor j=1:5for i=1:69Ft(i,j)=Ttq(i)*ig(j)*i0* nT/r;ua(i,j)=0.377*r* n( i)/(ig(j)*i0);Fz(i,j)=CdA*ua(i,j)A2/21.15+mz*g*f;endendplot(ua,Ft,ua,Ff,ua,Ff+Fw)title(' 汽车驱动力与行驶阻力平衡图');xlabel('ua(km/h)');ylabel('Ft(N)');gtext( 'Ft1')gtext( 'Ft2' )gtext( 'Ft3' )gtext( 'Ft4' )gtext( 'Ft5' )gtext( 'Ff+Fw')(2)求最大速度和最大爬坡度for k=1:175n1(k)=3300+k*0.1;Ttq(k)=-19.313+295.27*( n1(k)/1000)-165.44*( n1(k)/1000)A2+40.874* (n 1(k)/1000)A33.8445* (n 1(k)/1000)A4;Ft(k)=Ttq(k)*ig (5)*i0* nT/r;ua(k)=0.377*r* n1(k)/(ig(5)*i0);Fz(k)=CdA*ua(k)A2/21.15+mz*g*f;E(k)=abs((Ft(k)-Fz(k)));endfor k=1:175if(E(k)==mi n(E))disp(' 汽车最高车速=');disp(ua(k));disp('km/h');endendfor p=1:150n2(p)=2000+p*0.5;Ttq(p)=-19.313+295.27*( n2(p)/1000)-165.44*( n2(p)/1000)A2+40.874*( n2(p)/1000) 人3-3.8445*( n2(p)/1000)M;Ft(p)=Ttq(p)*ig(1)*iO* nT/r;ua(p)=0.377*r* n2(p)/(ig(1)*i0);Fz(p)=CdA*ua(p)A2/21.15+mz*g*f;af(p)=as in ((Ft(p)-Fz(p))/(mz*g));endfor p=1:150if(af(p)==max(af))i=ta n(af(p));disp(' 汽车最大爬坡度=');disp(i);endend汽车最高车速=99.0679km/h汽车最大爬坡度=0.3518(3)计算2档起步加速到70km/h所需时间for i=1:69n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*( n(i)/1000)-165.44*( n(i)/1000)A2+40.874*( n( i)/1000)A3-3.8 445*( n( i)/1000)A4;endfor j=1:5for i=1:69deta=1+lw/(mz*「A2)+lf*ig(j)A2*i0A2* nT/(mz*「A2);ua(i,j)=0.377*r* n( i)/(ig(j)*i0);a(i,j)=(Ttq(i)*ig(j)*i0* nT/r-CdA*ua(i,j)A2/21.15-mz*g*f)/(deta*mz);if(a(i,j)<=0)a(i,j)=a(i-1,j);endif(a(i,j)>0.05)b1(i,j)=a(i,j);u1(i,j)=ua(i,j);elseb1(i,j)=a(i-1,j);u1(i,j)=ua(i-1,j);endb(i,j)=1/b1(i,j);endendx1=u1(:,1);y 1= b(:,1);x2=u1(:,2);y2=b(:,2);x3=u1(:,3);y3=b(:,3);x4=u1(:,4);y4=b(:,4);x5=u1(:,5);y5=b(:,5);Plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5);title(' 加速度倒数时间曲线');axis([O 120 0 30]);xlabel('ua(km/h)');ylabel('1/aj');gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext('1/a5')加速度倒載时间曲线for i=1:69A=ua(i,3)-ua(69,2); if (A<1 &A>0)j=i;endB=ua(i,4)-ua(69,3);if(B<2&B>0)k=i;endif(ua(i,4)<=70)m=i;endend t=ua(1,2)*b(1,2); for p1=2:69t1(p1)=(ua(p1,2)-ua(p1-1,2))*(b(p1,2)+b(p1-1,2))*0.5; t=t+t1(p1);endfor p2=j:69t2(p2)=(ua(p2,3)-ua(p2-1,3))*(b(p2,3)+b(p2-1,3))*0.5; t=t+t2(p2);endfor p3=k:mt3(p3)=(ua(p3,4)-ua(p3-1,4))*(b(p3,4)+b(p3-1,4))*0.5; t=t+t3(p3);endt=t+(ua(j,3)-ua(69,2))*b(69,2)+(ua(k,4)-ua(69,3))*b(69,3) +(70-ua(m,4))*b(m,4);tz=t/3.6;disp(' 加速时间=');disp(tz);disp('s');加速时间=29.0585s二、计算与绘制题 1 中货车的1）汽车功率平衡图；2 ）最高档与次高档的等速百公里油耗曲线。
已知数据略。
（参见《汽车理论》习题第二章第7 题）解题程序如下: 用Matlab 语言m1=2000; m2=1800; mz=3880; g=9.81;r=0.367; CdA=2.77; f=0.013; nT=0.85;ig=[5.56 2.769 1.644 1.00 0.793];i0=5.83; If=0.218; Iw1=1.798; Iw2=3.598;n1=[815 1207 1614 2012 2603 3006 3403 3804];Iw=2*Iw1+4*Iw2;nd=400; Qid=0.299;for j=1:5for i=1:69 n(i)=(i+11)*50;Ttq(i)=-19.313+295.27*( n(i)/1000)-165.44*( n(i)/1000)人2+40.874*( n(。
/1000)人3-3.8 445*( n( 0/1000)人4;Pe(i)=n(i)*Ttq(i)/9549;ua(i,j)=0.377*r*n(i)/(ig(j)*i0);Pz(i,j)=(mz*g*f*ua(i,j)/3600.+CdA*ua(i,j)A3/76140.)/nT;endend plot(ua,Pe,ua,Pz);title(' 汽车功率平衡图)');xlabel('ua(km/h)');ylabel('Pe,Pz(kw)');gtext('I') gtext('II') gtext('III') gtext('IV') gtext('V') gtext('P 阻')for j=1:5for i=1:8Td(i)=-19.313+295.27*(n 1(i)/1000.0)-165.44*( n1(i)/1000.0)A 2+40.874*( n1(i)/10 00.0)人3-3.8445*( n1(i)/1000.0)A4;Pd(i)=n1(i)*Td(i)/9549;u(i,j)=0.377* n1(i)*r/(ig(j)*i0);end end b(1)=0.17768*Pd(1)A4-5.8629*Pd(1)A3+72.379*Pd(1)A2-416.46*Pd(1)+1326.8; b(2)=0.043072*Pd(2)A4-2.0553*Pd(2)A3+36.657*Pd (2)A2-303.98*Pd(2)+1354.7; b(3)=0.0068164*Pd(3)A4-0.51184*Pd (3)A3+14.524*Pd(3)A2-189.75*Pd (3)+1284.4; b(4)=0.0018555*Pd (4) A4-0.18517*Pd(4)A3+7.0035*Pd(4)A2-121.59*Pd(4)+1122.9; b(5)=0.00068906*Pd(5)A4-0.091077*Pd (5)^3+4.4763*Pd (5) A2-98.893*Pd(5)+1141.0;b(6)=0.00035032*Pd (6)A4-0.05138*Pd (6) A3+2.8593*Pd (6)A2-73.714*Pd(6)+1051.2; b(7)=0.00028230*Pd(7)A4-0.047449*Pd(7)A3+2.9788*Pd(7)A2-84.478*Pd(7)+1233.9; b(8)=-0.000038568*Pd(8)A40.00075215*Pd(8)A3+0.71113*Pd(8)A245.291*Pd(8) +1129.7; u1= u(:,1)'; u2=u(:,2)： u3=u(:,3)'; u4=u(:,4)'; u5=u(:,5)';B 仁 polyfit(u1,b,3);B2=polyfit(u2,b,3);B3=polyfit(u3,b,3);B4=polyfit(u4,b,3);B5=polyfit(u5,b,3);for q=1:69bh(q,1)=polyval(B1,ua(q,1));bh(q,2)=polyval(B2,ua(q,2));汽车功率平衡图 40 60ua/(krrVh) 80 m 120w o 90 70 Gn 5n 40 3n °0 20 MW Jbh(q,3)=polyval(B3,ua(q,3));bh(q,4)=polyval(B4,ua(q,4));bh(q,5)=polyval(B5,ua(q,5));endfor i=1:5for q=1:69Q(q,i)=Pz(q,i)*bh(q,i)/(1.02*ua(q,i)*7.05); end endplot(ua(:,4),Q(:,4),ua(:,5),Q(:,5));title(' 四档五档等速百公里油耗图');xlabel('ua(km/h)'); ylabel('Qs(L/100km)');三、改变1.3题中轻型货车的主减速器传动比，做出i0为5.17、5.43、 5.83、6.17、 6.33 时的燃油经济性—加速时间曲线，讨 论不同i0值对汽车性能的影响 Matlab 程序：m1=2000 ;m2=1800 ;m=3880 ;四档和五橙等速百公里油耗曲线40353025201510 020 40 60 B0 100 120ua/(km/h)r0=0.367 ;gt=0.85 ;f=0.013 ;CDA=2.77 ;i0=5.83 ;If=0.218 ;Iw1=1.798 ;Iw2=3.598 ;Ig5=[5.56 2.769 1.644 1.00 0.793];Ig0=[5.17 5.43 5.83 6.17 6.33];B=[1326.8 -416.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.00028230;1129.7 -45.291 0.71113 -0.00075215 -0.000038568;];n=[815 1207 1614 2012 2603 3006 3403 3804];for i=1:5for k=1:8ua(i,k)=0.377*0.367*n(k)/(Ig0(i)*Ig5(5));Ttq(i)=-19.313+295.27.*( n(i)/1000)-165.44.*( n( i)/1000)92+40.874.*( n( i)/1000)9 3-3.8445.* (n (。