汽车理论课程设计以及MATLAB培训

大学汽车理论课程设计以及MATLAB培训PPT课件( 90页)

MATLAB培训
出现，要求入PLP
MATLAB培训
这时在matlab7.0安装文件夹中找到 keygen文件，运行
MATLAB培训
出现，复制生产的密码，粘贴到MATLAB的PLP 中
MATLAB培训
然后Next，出现选Yes，然后Next
MATLAB培训
出现选Custom，然后Next
将setup为安装文件，序列号在sn文件里面
MATLAB6.5的安装
双击setup安装，选next
MATLAB6.5的安装
出现
打开sn 文件，复制序列号填入，点 next
MATLAB6.5的安装
出现
点yes
MATLAB6.5的安装
出现，在里面随便填几个字母
点next，出现，选no
出现
MATLAB培训
把matlab7.0安装文件夹都复制到硬盘，大约占用1.04G。建议不要复制到安装操作系统的分区，记住位置
MATLAB培训
点击右下角图标运行daemon，找到第一虚拟光驱文件，选中文件EFAML71，打开，MATLAB安装程序自动运行
Win7操作系统，直接安装可能会不兼容。
注：MATLAB界面可以更改，如果想回到默认界面，依次选中菜单：Desktop/ Desktop
Layout/Default
MATLAB基本操作
命令窗口
各个版本的界面可能有差别，但是一定有 Command命令窗口
MATLAB基本操作
用MATLAB进行计算有两种方法
可以在命令窗口里面直接输入命令，可以立即执行。用于对程序的调试，课程设计不建议用这种方法
安装过程中，可能提示安装第二张光盘，这时点右下角Deamon图标

汽车理论课程设计：基于Matlab的汽车动力性的仿真

2009 届海南大学机电工程学院汽车工程系汽车理论课程设计题目：汽车动力性的仿真学院：机电工程学院专业：09级交通运输姓名：黄生锐学号：20090504指导教师：编号名称件数页数编号名称件数页数1 课程设计论文 1 3Matlab编程源程序 12 设计任务书 12012年6月20日成绩汽车理论课程设计任务书姓名黄生锐学号20090504 专业09交通运输课程设计题目汽车动力性的仿真内容摘要：本设计的任务是对一台Passat 1.8T手动标准型汽车的动力性能进行仿真。

采用MATLAB编程仿真其性能，其优点是：一是能过降低实际成本，提高效率；二是获得较好的参数模拟，对汽车动力性能提供理论依据。

主要任务：根据该车的外形、轮距、轴距、最小离地间隙、最小转弯半径、车辆重量、满载重量以及最高车速等参数，结合自己选择的适合于该车的发动机型号求出发动机的最大功率、最大扭矩、排量等重要的参数。

并结合整车的基本参数，选择适当的主减速比。

依据GB、所求参数，结合汽车设计、汽车理论、机械设计等相关知识，计算出变速器参数，进行设计。

论证设计的合理性。

设计要求：1、动力性分析：1)绘制汽车驱动力与行驶阻力平衡图；2)求汽车的最高车速、最大爬坡度；3)用图解法或编程绘制汽车动力特性曲线4)汽车加速时间曲线。

2、燃油经济性分析：1) 汽车功率平衡图；完成内容：1．Matlab编程汽车驱动力与行驶阻力平衡图2．编程绘制汽车动力特性曲线图3．编程汽车加速时间曲线图4．课程设计论文1份汽车动力性仿真摘要本文是对Passat 1.8T 手动标准型汽车的动力性能采用matlab 编制程序，对汽车动力性进行计算。

从而对汽车各个参数做出准确的仿真研究，为研究汽车动力性提供理论依据，本文主要进行的汽车动力性仿真有：最高车速、加速时间和最大爬坡度。

及相关汽车燃油性经济。

关键词：汽车；动力性；试验仿真；matlab1. Passat 1.8T 手动标准型汽车参数功率Pe （kw ）转速n （r/min ）15 1000 36 1750 50 2200 66 2850 80 3300 90 4000 110 5100 1055500各档传动比主减速器传动比第1档 3.665 4.778第2档 1.999 第3档 1.407 第4档 1 第5档 0.472 车轮半径0.316（m ）传动机械效率0.91 假设在良好沥青或水泥路面上行驶，滚动阻力系数 0.014 整车质量1522kgC D A2.4m22. 最高车速汽车的最高车速是指汽车标准满载状态，在水平良好的路面(混凝土或沥青路面)上所能达到的最高行驶速度。

汽车理论习题答案及MATLAB编程

汽车理论课后习题答案（部分题目）教材：汽车理论（第5版）清华大学余志生主编机械工业出版社出版本文档包含习题答案：第一章：1.1 1.2 1.7 1.8第二章：2.2 2.3 2.4第三章第四章：4.1 4.2 4.3 4.5第五章：5.1 5.8 5.10 5.11 5.12 5.14 5.17第六章：6.1~6.5其中第三章 4.3 5.11 6.2 6.4 6.5所要求的MATLAB编程见本作者另一百度文库文档。

获取方法：点击左上角本作者的昵称，进入本作者的文库，《汽车理论课后习题答案—附MATLAB编程》第三章：解：（完整解题含MATLAB编程见另一文件）根据“货车最高车速为100km/h左右”，为满足五种主减速器与两种变速器匹配后车速均能达到100km/h，参照表中发动机转速，可选3006r/min、3403r/min、3804r/min三种。

结合万有特性图，考虑到动力性与燃油经济性的统一，选择转速3403r/min进行计算。

由于本人能力有限，暂难以完成EPA循环工况下的燃油消耗量计算，故以“货车采用最高档、90％负荷率工况下”的等速百公里油耗代替表示燃油经济性。

取汽油密度为0.7kg/L。

燃油经济性以等速百公里油耗Qs(L/100km)表示；动力性以2挡起步加速行驶至96.6km/h车速所用时间t(s)表示。

下面运用MATLAB编程实现计算和绘图，编程（m文件）列于最后，运行结果（4挡）如下：t=48.4235t=46.3511t=41.5557t=33.3796t=30.556490％负荷Pe=55.5647Qs4=25.103126.365628.307829.958730.7356图1加速度倒数-速度曲线图2燃油经济性-加速时间曲线由此可得出：主减速器传动比越大，动力性越好，燃油经济性越差。

与4挡变速器相比，采用5挡变速器时动力性和燃油经济性有所提高和改善。

（完整解题含MATLAB 编程见另一文件）图1图25.8某种汽车的质心位置、轴距和前后轮胎的型号已定。

汽车理论课程设计

汽车理论课程设计姓名学号专业班级指导教师年月日汽车理论课程设计任务书本次计算设计以小组为单位进行，每组计算两种车型(大型车、小型车)。

先进行手工计算并绘图（每个档位选取5-8个特征点），然后计算机编程（MATLAB）实现并绘图，最终打印计算说明书和程序。

具体设计要求如下：1.汽车动力性、经济性分析计算通过查阅收集有关资料，计算分析给定型号汽车的动力性能及燃油经济性，并绘出该车型的发动机外特性曲线，驱动力——行驶阻力平衡图，动力特性图，百公里油耗曲线。

根据计算结果，结合实际情况，分析该车型发动机参数和底盘性能参数匹配是否合理，并提出修改意见。

2.参数部分车型参数不完整，请查阅相关资料或用经验公式计算选取，并经手动计算分析后修正获得。

（说明书中注明出处）3.说明书说明书内容包括：任务书、目录、各车型参数、参数分析、计算、图表、结论、设计总结及心得体会。

车型一：解放CA1091载货汽车一、发动机CA6102（附表一）Nmax=99kw(相应转速3000r/min)Mmax=373N.m(相应转速1300r/min)二、整车参数：1.尺寸参数：全长L=7205mm，全宽B=2476mm，全高H=2436mm，轴距L1=4050mm,前轮距B1=1850mm，后轮距B2=1740mm.2.重量参数（附表二）3.性能参数：变速箱传动比i1=7.64，i2=4.835，i3=2.856，i4=1.895，i5=1.377，i6=1，i 倒=7.66。

主减速器比i0=6.33。

车轮：9.00-20。

三、使用数据：滚动阻力系数f=0.03；道路阻力系数：强度计算用Φ=1性能计算用Φ=0.8空气阻力系数：Cd=0.8；迎风面积：A=0.78×宽×高；最大速度：Vmax=90km/h；最大爬坡度：28%；传动系效率：η=0.9表一：发动机参数车型四：SATANA2000轿车一、发动机2VQS（附表一）Nmax=74kw(相应转速5200r/min)Mmax=155N.m(相应转速3800r/min)二、整车参数：1.尺寸参数：全长L=4680mm，全宽B=1700mm，全高H=1423mm，轴距L1=2656mm,前轮距B1=1414mm，后轮距B2=1422mm.2.重量参数（附表二）3.性能参数：变速箱传动比i1=3.45，i2=1.94，i3=1.29，i4=0.985，i5=0.8,主减速器比i0=4.1。

汽车理论matlab编程,汽车理论课后作业matlab编程详解带注释[10页]

汽车理论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编程

汽车理论课后习题答案—附MATLAB编程教材：汽车理论（第5版）清华大学余志生主编机械工业出版社出版本文档包含第三章 4.3 5.11 6.2 6.4 6.5所要求的MATLAB编程其它习题答案：第一章：1.1 1.2 1.7 1.8第二章：2.2 2.3 2.4第三章第四章：4.1 4.2 4.3 4.5第五章：5.1 5.8 5.10 5.11 5.12 5.14 5.17第六章：6.1~6.5见本作者另一百度文库文档。

获取方法：点击左上角本作者的昵称，进入本作者的文库，《汽车理论课后习题答案（含MATLAB编程）》第三章：m=3880;%总质量g=9.8;%重力加速度r=0.367;%车轮半径eta_t=0.85;%传动系机械效率f=0.013;%滚动阻力系数CDA=2.77;%空气阻力系数*迎风面积i0=[5.175.435.836.176.33];%主减速器传动比If=0.218;%飞轮转动惯量Iw1=1.798;Iw2=3.598;%两前轮/四后轮转动惯量Iw=Iw1+Iw2;ig4=[6.09 3.091.711];%4挡变速器传动比ig5=[5.56 2.7691.64410.793];%5挡变速器传动比n=600:1:4000;%发动机转速Tq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;%4挡发动机for i=1:5Ft1=Tq*ig4(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig4(2)*i0(i)/r;Ft3=Tq*ig4(3)*i0(i)/r;Ft4=Tq*ig4(4)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig4(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig4(2));u3=0.377*r*n/(i0(i)*ig4(3));u4=0.377*r*n/(i0(i)*ig4(4));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;delta1=1+(Iw+If.*ig4(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig4(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig4(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig4(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;figure(1)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=96.6;x1=[];x2=[];x3=[];x4=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];elseif u4(k)<=u4maxx3=[k];endendq2=max(x2);q3=max(x3);ua4=u4(q2:q3);a4=b4(q2:q3);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s=[s2s3s4];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs4(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs4=');disp(Qs4);tt=[48.423546.351141.555733.379630.5564];figure(2)plot(Qs4,tt,'*')hold onxi=25:0.001:31;pp=interp1(Qs4,tt,xi,'cubic');plot(xi,pp)title('燃油经济性-加速时间曲线');xlabel('燃油经济性--等速百公里油耗(L/100km)');ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')%5挡发动机for i=1:5Ft1=Tq*ig5(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig5(2)*i0(i)/r;Ft3=Tq*ig5(3)*i0(i)/r;Ft4=Tq*ig5(4)*i0(i)/r;Ft5=Tq*ig5(5)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig5(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig5(2));u3=0.377*r*n/(i0(i)*ig5(3));u4=0.377*r*n/(i0(i)*ig5(4));u5=0.377*r*n/(i0(i)*ig5(5));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;F5=m*g*f+CDA*u5.^2/21.15;delta1=1+(Iw+If.*ig5(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig5(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig5(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig5(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta5=1+(Iw+If.*ig5(5).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);a5=(Ft5-F5)/(delta5*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;b5=1./a5;figure(3)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4,u5,b5);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext('1/a5')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=max(u4);u5min=u4max;u5max=96.6;x1=[];x2=[];x3=[];x4=[];x5=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];endendq2=max(x2);ua4=u4(q2:y);a4=b4(q2:y);for l=1:3401if u5(l)<=u5minx3=[l];elseif u5(l)<=u5maxx4=[l];endendq2=max(x2);q3=max(x3);q4=max(x4);ua5=u5(q3:q4);a5=b5(q3:q4);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s5=trapz(ua5,a5);%五挡运行时间s=[s2s3s4s5];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;Pe5=Ft5.*u5./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs5(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs5=');disp(Qs5);tt=[93.973058.283444.678445.366744.9793];figure(2)plot(Qs5,tt,'*')hold onxi=30:0.001:100;pp=interp1(Qs5,tt,xi,'cubic');plot(xi,pp)%title('燃油经济性-加速时间曲线');%xlabel('燃油经济性--等速百公里油耗(L/100km)');%ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')运行结果：t=48.4235t=46.3511t=41.5557t=33.3796t=30.556490％负荷Pe=55.5647Qs4=25.103126.365628.307829.958730.7356图1加速度倒数-速度曲线图2燃油经济性-加速时间曲线4.3%空载(no load)—1；满载(full load)—2m1=4080;m2=9290;hg1=0.845;hg2=1.170;L=3.950;a1=2.100;a2=2.950;b1=L-a1;b2=L-a2;beta=0.38;%利用附着系数与制动强度的关系曲线z=0:0.01:1;phi_f1=L*beta.*z./(b1+z*hg1);%前轮利用附着系数(空载)phi_r1=L*(1-beta).*z./(a1-z*hg1);%空载时后轮利用附着系数(空载)phi_f2=L*beta.*z./(b2+z*hg2);%前轮利用附着系数(满载)phi_r2=L*(1-beta).*z./(a2-z*hg2);%后轮利用附着系数(满载)phi=z;%最理想情况figure(1);plot(z,phi_f1,'g--',z,phi_f2,'g-',z,phi_r1,'b--',z,phi_r2,'b-',z,phi,'k-'); axis([0,1,0,1]);grid on;box off;axis square;title('利用附着系数与制动强度的关系曲线');xlabel('制动强度z/g');ylabel('利用附着系数{\phi}');%ECE法规要求界限hold on;z1=0.2:0.01:0.8;z2=0.15:0.01:0.3;z3=0.3:0.01:1;phi1=(z1+0.07)/0.85;phi2=z2-0.08;phi3=z2+0.08;phi4=(z3-0.02)/0.74;plot(z1,phi1,'r-.',z2,phi2,'r-.',z2,phi3,'r-.',z3,phi4,'r-.'); legend('{\phi}_{f}(空载)','{\phi}_{f}(满载)','{\phi}_{r}(空载)',...'{\phi}_{r}(满载)','{\phi}=z','ECE法规要求界限');legend('Location','northwest');%制动效率曲线Ef1=z./phi_f1*100;Er1=z./phi_r1*100;Ef2=z./phi_f2*100;Er2=z./phi_r2*100;figure(2);plot(phi_f1,Ef1,'k',phi_r1,Er1,'k',phi_f2,Ef2,'b',phi_r2,Er2,'b'); axis([0,1,0,100]);grid on;box off;axis square;title('制动效率曲线');xlabel('附着系数{\phi}');ylabel('制动效率(%)');text(0.25,90,'E_f');text(0.6,88,'满载');text(0.8,90,'E_r');text(0.7,65,'空载');text(0.72,72,'E_r');程序运行结果如下：图1图2 5.11L=3.048;K=0.0024;i=20;ua=0:0.01:120;f=1/L.*(1./(3.6./ua+K.*ua./3.6));syms uF(u)=(1/L.*(1./(1./u+K.*u)))./i;plot(ua,f,'b-');xlabel('u_a(km/h)');ylabel('\omega_r/\delta)_s');F(22.35)运行结果：ans=0.1667386094155036.2f=0.1:0.1:100;Gq1=5.12*10^(-5)./f.^2;Gq2=2.02*10^(-3).*f./f;Gq3=7.98*10^(-2).*f.^2;figure(1);loglog(f,Gq1,'b-');box off;axis square;title('位移功率谱密度'); xlabel('f/Hz');ylabel('G_q(f)/(m^2s)');figure(2);loglog(f,Gq2,'b-');box off;axis square;title('速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''(f)/(m^2/s)');figure(3);loglog(f,Gq3,'b-');box off;axis square;title('加速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''_''(f)/(m^2/s^3)');运行结果如下：6.4f1=0.1:0.01:1;f2=1:0.01:10;G1=0.010106.*f1.^2;G2=0.010106.*f2./f2;loglog(f1,G1,'b-');hold onloglog(f2,G2,'b-');title('车身加速度的功率谱密度'); xlabel('f/Hz');ylabel('G_z_''_''(f)/m^-^1'); axis([0.1,10,0.0001,0.1]);运行结果如下：6.51.计算说明1）①幅频特性()122222114z q λζλγ⎡⎤-+⎢=⎢⎥∆⎣⎦，其中0ωλω=，()2222222111141λλγζλγλμμ⎡⎤⎡⎤⎛⎫⎛⎫∆=-+--+-⎢⎥⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎣⎦⎣⎦；()12222222211414z z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦；()1222222221414s s s s s p z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦，其中s s ωλω=；②均方根谱()()221~j 2z q H f f ωππ=；()()222~j 2z qH f f ωππ== ；()()2221~1j 2p qz z H f f z qωππ== ；③其他值()0.5360q q G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536110z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536220z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.5360a a G f df σ⎡⎤=⎢⎥⎣⎦⎰()()()()()()13622w 012222412.5362202412.5d 12.5 0.5d d 1d d 4a a a a a a W f G f f f G f f G f f G f f G f f f ⎡⎤=⎢⎥⎣⎦⎡⎤⎛⎫⎛⎫=+++⎢⎥ ⎪ ⎪⎝⎭⎢⎥⎝⎭⎣⎦⎰⎰⎰⎰⎰其中，()()()()()0.5 0.52 2441 412.512.512.536f f f W f f f f <<⎧⎪⎪<<⎪=⎨<<⎪⎪<<⎪⎩()aw w 020lg L a =2）见1）中计算式3）()()()22004q q q G f G G n n uωπ== ()()()1122363622d ~00d d dd f q q f q f H f G f f G f f q σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122d 12f q f γλπ⎡⎤=⎢⎥∆⎣⎦ ()()()1122363622d //~00d d dd F Gq q F G q F H f G f f G f f Gq σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122222d 1412F f Gq g λζλμπγ⎡⎤⎛⎫-+⎢⎥⎪+⎢⎥⎝⎭=⎢⎥∆⎢⎥⎢⎥⎣⎦其余见1）中计算式2.程序清单1）f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi*fs;lambda=omega./omega0;lambda_s=omega./omega_s;%计算并绘制幅频特性delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;figure(1)loglog(f,z1_q,'b-',0:0.1:100,(0:0.1:100)./(0:0.1:100),'g-');axis([0.1,100,0.1,10]);grid on ;title('z_1~q 的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_1}{q}|$$','Interpreter','latex');gtext('0:1');figure(2)loglog(f,z2_z1,'b-',0:0.1:1,(0:0.1:1)./(0:0.1:1),'g-',1:0.1:100,1./...(1:0.1:100),'g-');axis([0.1,100,0.01,10]);grid on;title('z_2~z_1的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_2}{z_1}|$$','Interpreter','latex'); gtext('0:1');gtext('-1:1');figure(3)loglog(f,p_z2,'b-');axis([0.1,100,0.01,10]);grid on;title('p~z_2的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{p}{z_2}|$$','Interpreter','latex');%计算并绘制均方根谱sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q; figure(4)loglog(f,sqrt_Gz1,'b-');grid on;title('车轮位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''1}(f)}$','Interpreter','latex'); figure(5)loglog(f,sqrt_Gz2,'b-');grid on;title('车身位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''2}(f)}$','Interpreter','latex'); figure(6)loglog(f,sqrt_Ga,'b-');grid on;title('传至人体的位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_a(f)}$','Interpreter','latex');%计算其它值sigma_q=sqrt(trapz(f,sqrt_Gq.^2));%路面不平度加速度均方根值sigma_z1=sqrt(trapz(f,sqrt_Gz1.^2));%车轮加速度均方根值sigma_z2=sqrt(trapz(f,sqrt_Gz2.^2));%车身加速度均方根值sigma_a=sqrt(trapz(f,sqrt_Ga.^2));%传至人体的加速度均方根值for i=1:Nif f(i)<=2W(i)=0.5;elseif f(i)<=4W(i)=f(i)/4;elseif f(i)<=12.5W(i)=1;elseW(i)=12.5/f(i);endendaw=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));%加权加速度均方根值a0=10^(-6);Law=20*log10(aw/a0);%加权振级format shortdisp('路面不平度加速度均方根值=');disp(sigma_q);disp('车轮加速度均方根值=');disp(sigma_z1);disp('车身加速度均方根值=');disp(sigma_z2);disp('传至人体的加速度均方根值=');disp(sigma_a);disp('加权加速度均方根值=');disp(aw);disp('加权振级=');disp(Law);2）%随fs变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=1.5:0.01:6;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;lambda=omega./omega0;for i=1:length(fs)omega_s=2*pi.*fs(i);lambda_s=omega./omega_s;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(fs,aw);grid ontitle('a_w随f_s的变化曲线');xlabel('f_s/Hz');ylabel('a_w/m·s^-^2');figure(2)plot(fs,Law);grid ontitle('L_aw随f_s的变化曲线');xlabel('f_s/Hz');ylabel('L_a_w/dB');%随zeta_s变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_ss=0.125:0.001:0.5;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi.*fs;lambda=omega./omega0;lambda_s=omega./omega_s;for i=1:length(zeta_ss)zeta_s=zeta_ss(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(zeta_ss,aw);grid ontitle('a_w随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('a_w/m·s^-^2');figure(2)plot(zeta_ss,Law);grid ontitle('L_aw随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('L_a_w/dB');3）%随f0变化f0=0.25:0.01:3;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;lambda_s=omega./omega_s;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(f0)omega0=2*pi*f0(i);lambda=omega./omega0;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(f0==1.5);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(f0,sm_z2,f0,sm_fd,f0,sm_FdG);axis([0.25,3,-25,10]);grid on;title('各响应量均方根值随f_0变化的曲线');xlabel('f_0/Hz');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随zeta变化f0=1.5;zeta0=0.125:0.001:0.5;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(zeta0)zeta=zeta0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(zeta0==0.25);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(zeta0,sm_z2,zeta0,sm_fd,zeta0,sm_FdG);axis([0.125,0.5,-4,4]);grid on;title('各响应量均方根值随\zeta变化的曲线');xlabel('\zeta');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');%随gamma变化f0=1.5;zeta=0.25;gamma0=4.5:0.1:18;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(gamma0)gamma=gamma0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(gamma0==9);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(gamma0,sm_z2,gamma0,sm_fd,gamma0,sm_FdG);axis([4.5,18,-6,6]);grid on;title('各响应量均方根值随\gamma变化的曲线');xlabel('\gamma');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随mu变化f0=1.5;zeta=0.25;gamma=9;mu0=5:0.1:20;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(mu0)mu=mu0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(mu0==10);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(mu0,sm_z2,mu0,sm_fd,mu0,sm_FdG);axis([5,20,-2,2]);grid on;title('各响应量均方根值随\mu变化的曲线');xlabel('\mu');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB');legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');3.结果分析1）路面不平度加速度均方根值=0.3523车轮加速度均方根值=0.2391车身加速度均方根值=0.0168传至人体的加速度均方根值=0.0161加权加速度均方根值=0.0100加权振级=80.0287分析：根据课本中表6-2，a w=0.0100<0.315，L aw=80.0287<110，故乘客没有不舒适。

matlab汽车经济性课程设计

matlab汽车经济性课程设计一、课程目标知识目标：1. 学生能理解并掌握MATLAB软件的基本操作和应用；2. 学生能运用MATLAB进行汽车经济性的相关数据分析；3. 学生了解并掌握汽车经济性评价的主要指标及其计算方法；4. 学生了解我国汽车经济性标准及政策。

技能目标：1. 学生能够运用MATLAB软件进行数据处理和图像绘制；2. 学生能够运用所学知识对汽车经济性进行简单评估和分析；3. 学生能够通过小组合作，完成汽车经济性课程设计任务。

情感态度价值观目标：1. 学生培养对汽车工程学科的兴趣和热情；2. 学生树立节能环保、绿色出行的意识，关注汽车产业发展；3. 学生养成合作、探究的学习习惯，培养创新思维和解决问题的能力。

课程性质：本课程为实践性较强的课程，结合MATLAB软件，让学生在实际操作中学习汽车经济性相关知识。

学生特点：高中年级学生对汽车有一定了解，具备基本的数据分析能力，对实际应用有较高的兴趣。

教学要求：教师需注重理论与实践相结合，引导学生主动参与，关注学生个体差异，提高学生的动手能力和团队协作能力。

通过课程学习，使学生达到上述设定的知识、技能和情感态度价值观目标，为后续相关课程打下坚实基础。

二、教学内容1. MATLAB软件基本操作与应用：包括数据类型、矩阵运算、函数编写和绘图功能等，涉及课本第二章内容。

2. 汽车经济性评价指标：介绍油耗、排放、动力性等指标，结合课本第三章相关内容。

3. 汽车经济性计算方法：讲解等速油耗、综合油耗的计算方法，参考课本第四章内容。

4. 汽车经济性数据分析：利用MATLAB对汽车经济性数据进行处理、分析和图像绘制，结合课本第五章实例。

5. 我国汽车经济性标准及政策：解读相关标准及政策，涉及课本第六章内容。

教学大纲安排：第一周：MATLAB软件基本操作与应用学习；第二周：汽车经济性评价指标介绍；第三周：汽车经济性计算方法学习；第四周：汽车经济性数据分析实践；第五周：我国汽车经济性标准及政策解读。

汽车理论第二版课程设计

汽车理论第二版课程设计1. 课程简介本课程是汽车理论的第二版课程设计，旨在通过理论知识和实践操作，深入了解汽车的构造、运行原理、维修和保养等方面的知识。

本课程将以课堂讲解和实操为主，为学员提供全面而深入的汽车知识。

2. 课程安排2.1 第一周课程目的了解汽车的基本概念和结构，初步认识发动机的工作原理。

课程内容•汽车的基本结构和分类•发动机的工作原理和分类•发动机的零部件和功能学习任务•理解汽车基本结构和分类•掌握发动机的工作原理和分类•熟悉发动机的零部件和功能2.2 第二周课程目的学习汽车的传动系统和底盘系统，了解其工作原理。

课程内容•传动系统的组成和工作原理•底盘系统的组成和工作原理•刹车系统的组成和工作原理学习任务•熟悉传动系统和底盘系统的组成和工作原理•理解刹车系统的组成和工作原理2.3 第三周课程目的学习汽车的电路系统和车身系统，了解其基本结构和工作原理。

课程内容•电路系统的组成和工作原理•车身系统的组成和工作原理•灯光系统的组成和工作原理学习任务•理解汽车电路系统的组成和工作原理•掌握车身系统和灯光系统的组成和工作原理2.4 第四周课程目的学习汽车的维修和保养，掌握相关知识和技能。

课程内容•汽车的日常维护和保养•常见故障的排除方法•汽车维修的基本流程和技能学习任务•掌握汽车日常维护和保养的方法•掌握常见故障的排除方法•学习汽车维修的基本流程和技能3. 实践操作本课程将安排实践操作环节，学员必须参加实践操作并通过相应的考核。

实践操作将包括汽车的基础维修和保养、故障排除和部分组装。

4. 课程评估本课程的评估方式包括课堂出勤率、实践操作、作业和期末考试，总成绩将按一定比例综合评定。

其中，实践操作和期末考试为必考核项，作业和课堂出勤率将作为参考。

5. 教学资料•《汽车理论第二版》教材•相关课件和PPT•汽车维修手册和资料•实践操作相关工具和设备6. 教学方法本课程将采用讲授、案例分析、实践操作等多种教学方法，注重理论与实践相结合，注重启发式教学和问题导向式教学。

汽车理论(余志生)matlab编程实验报告

汽车理论1.汽车驱动力-行驶阻力平衡图的绘制代码：>> n=600:1:4000;>>Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;>> i0=5.83;>> ig1=5.56;>> ig2=2.769;>> ig3=1.644;>> ig4=1;>> ig5=0.793;>> Ft1=Tq*ig1*i0*0.85/0.367;>> Ft2=Tq*ig2*i0*0.85/0.367;>> Ft3=Tq*ig3*i0*0.85/0.367;>> Ft4=Tq*ig4*i0*0.85/0.367;>> Ft5=Tq*ig5*i0*0.85/0.367;>> u1=0.377*n*0.367/(5.83*ig1);>> u2=0.377*n*0.367/(5.83*ig2);>> u3=0.377*n*0.367/(5.83*ig3);>> u4=0.377*n*0.367/(5.83*ig4);>> u5=0.377*n*0.367/(5.83*ig5);>> u=0:1:120;>> Ff=3880*9.8*0.013+2.77*u.^2/21.15;>> plot(u1,Ft1,u2,Ft2,u3,Ft3,u4,Ft4,u5,Ft5,u,Ff);2.最高车速、最大爬坡度及附着率：最高车速时Ft=Ff+Fw，由曲线图可知在5档时有最高车速：U max=(Ft5-3880*9.8*0.013*21.15/2.77)2; n max=4000r/min Umax=99.08km/h最大爬坡度：由曲线图可知最大Ft在一档，则计算一档时最大爬坡度为：i max=tan(arcsin F t−F f−F wG )=tan arcsin(F t1−F fG)由matlab计算得最大爬坡度i max=0.3522.在之前代码基础上输入：>> Ff1=3880*9.8*0.013+2.77*u1.^2/21.15; >> i=tan(asin((Ft1-Ff1)/(3880*9.8)));>> plot(u1,i);>> [m,index]=max(i);最大附着率：Φ=q a L +ℎg L ∗q =(0.35221.9473.2+0.93.2∗0.3522)≈0.498→0.53.加速度倒数曲线与加速时间：在以上代码继续输入matlab 代码：>>Ff2=3880*9.8*0.013+2.77*u2.^2/21.15;>>Ff3=3880*9.8*0.013+2.77*u3.^2/21.15;>>Ff4=3880*9.8*0.013+2.77*u4.^2/21.15;>> Z1=1+((1.798+3.598)+(0.218*5.83.^2*0.85*ig1.^2))/(3880*0.367.^2); %设置旋转质量换算系数>> Z2=1+((1.798+3.598)+(0.218*5.83.^2*0.85*ig2.^2))/(3880*0.367.^2); >> Z3=1+((1.798+3.598)+(0.218*5.83.^2*0.85*ig3.^2))/(3880*0.367.^2); >> Z4=1+((1.798+3.598)+(0.218*5.83.^2*0.85*ig4.^2))/(3880*0.367.^2); >> Z5=1+((1.798+3.598)+(0.218*5.83.^2*0.85*ig5.^2))/(3880*0.367.^2); >> a1=Z1*3880./(Ft1-Ff1); % 设置加速度倒数>> a2=Z2*3880./(Ft2-Ff2);>> a3=Z3*3880./(Ft3-Ff3);>> a4=Z4*3880./(Ft4-Ff4); % 由上题知5档最高车速99km/h ，因此重设转速、扭矩、5档驱动力和阻力>> n=600:1:3308;>> Tqa=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;>> Ft55=Tqa*ig5*i0*0.85/0.367;u55=0.377*n*0.367/(5.83*ig5);>> Ff55=3880*9.8*0.013+2.77*u55.^2/21.15;>> a5=Z5*3880./(Ft55-Ff55);>> plot(u1,a1,u2,a2,u3,a3,u4,a4,u55,a5);axis([0,100,0,14]); % 最高车速时Ft≈Ff，1a<14无限大，因此设置画图范围0<u a<100,0<1a积分法计算最大二档加速到70km/h的时间：t=26.2s脚本代码：Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;ig=[5.56,2.769,1.644,1,0.793];Z=[1.3829,1.1027,1.0429,1.0224,1.0179];nmin=600;nmax=4000;u1=0.377*0.367*nmin./(ig*5.83);u2=0.377*0.367*nmax./(ig*5.83);ua=[6:0.01:99];N=length(ua);Fw=2.77*ua.^2/21.15;for i=1:Nk=i;if ua(i)<=u2(2)n=ua(i)*(ig(2)*5.83/0.367)/0.377;Ft=Tq*ig(2)*5.83*0.85/0.367;IA(i)=(Z(2)*3880)/(Ft-3880*9.8*0.013-Fw(i));da(i)=0.01*IA(i)/3.6;elseif ua(i)<=u2(3)n=ua(i)*(ig(3)*5.83/0.367)/0.377;Ft=Tq*ig(3)*5.83*0.85/0.367;IA(i)=(Z(3)*3880)/(Ft-3880*9.8*0.013-Fw(i));da(i)=0.01*IA(i)/3.6;elseif ua(i)<=u2(4)n=ua(i)*(ig(4)*5.83/0.367)/0.377;Ft=Tq*ig(4)*5.83*0.85/0.367;IA(i)=(Z(4)*3880)/(Ft-3880*9.8*0.013-Fw(i));da(i)=0.01*IA(i)/3.6;elsen=ua(i)*(ig(5)*5.83/0.367)/0.377;Ft=Tq*ig(5)*5.83*0.85/0.367;IA(i)=(Z(5)*3880)/(Ft-3880*9.8*0.013-Fw(i));da(i)=0.01*IA(i)/3.6;enda=da(1:k);t(i)=sum(a);endplot(t,ua);速度-时间图像如下：取平均值得t=26.2s。

汽车理论课后习题答案及MATLAB编程

汽车理论课后习题答案—附MATLAB编程教材：汽车理论（第5版）清华大学余志生主编机械工业出版社出版本文档包含第三章 4.3 5.11 6.2 6.4 6.5所要求的MATLAB编程其它习题答案：第一章：1.1 1.2 1.7 1.8第二章：2.2 2.3 2.4第三章第四章：4.1 4.2 4.3 4.5第五章：5.1 5.8 5.10 5.11 5.12 5.14 5.17第六章：6.1~6.5见本作者另一百度文库文档。

获取方法：点击左上角本作者的昵称，进入本作者的文库，《汽车理论课后习题答案（含MATLAB编程）》第三章：m=3880;%总质量g=9.8;%重力加速度r=0.367;%车轮半径eta_t=0.85;%传动系机械效率f=0.013;%滚动阻力系数CDA=2.77;%空气阻力系数*迎风面积i0=[5.175.435.836.176.33];%主减速器传动比If=0.218;%飞轮转动惯量Iw1=1.798;Iw2=3.598;%两前轮/四后轮转动惯量Iw=Iw1+Iw2;ig4=[6.09 3.091.711];%4挡变速器传动比ig5=[5.56 2.7691.64410.793];%5挡变速器传动比n=600:1:4000;%发动机转速Tq=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4;%4挡发动机for i=1:5Ft1=Tq*ig4(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig4(2)*i0(i)/r;Ft3=Tq*ig4(3)*i0(i)/r;Ft4=Tq*ig4(4)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig4(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig4(2));u3=0.377*r*n/(i0(i)*ig4(3));u4=0.377*r*n/(i0(i)*ig4(4));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;delta1=1+(Iw+If.*ig4(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig4(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig4(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig4(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;figure(1)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=96.6;x1=[];x2=[];x3=[];x4=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];elseif u4(k)<=u4maxx3=[k];endendq2=max(x2);q3=max(x3);ua4=u4(q2:q3);a4=b4(q2:q3);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s=[s2s3s4];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs4(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs4=');disp(Qs4);tt=[48.423546.351141.555733.379630.5564];figure(2)plot(Qs4,tt,'*')hold onxi=25:0.001:31;pp=interp1(Qs4,tt,xi,'cubic');plot(xi,pp)title('燃油经济性-加速时间曲线');xlabel('燃油经济性--等速百公里油耗(L/100km)');ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')%5挡发动机for i=1:5Ft1=Tq*ig5(1)*i0(i)/r;%各档对应转速下的驱动力Ft2=Tq*ig5(2)*i0(i)/r;Ft3=Tq*ig5(3)*i0(i)/r;Ft4=Tq*ig5(4)*i0(i)/r;Ft5=Tq*ig5(5)*i0(i)/r;u1=0.377*r*n/(i0(i)*ig5(1));%各档对应转速下的车速u2=0.377*r*n/(i0(i)*ig5(2));u3=0.377*r*n/(i0(i)*ig5(3));u4=0.377*r*n/(i0(i)*ig5(4));u5=0.377*r*n/(i0(i)*ig5(5));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;F5=m*g*f+CDA*u5.^2/21.15;delta1=1+(Iw+If.*ig5(1).^2.*i0(i).^2.*eta_t)./(m.*r.^2);%汽车旋转质量换算系数delta2=1+(Iw+If.*ig5(2).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta3=1+(Iw+If.*ig5(3).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta4=1+(Iw+If.*ig5(4).^2.*i0(i).^2.*eta_t)./(m.*r.^2);delta5=1+(Iw+If.*ig5(5).^2.*i0(i).^2.*eta_t)./(m.*r.^2);a1=(Ft1-F1)/(delta1*m);%加速度a2=(Ft2-F2)/(delta2*m);a3=(Ft3-F3)/(delta3*m);a4=(Ft4-F4)/(delta4*m);a5=(Ft5-F5)/(delta5*m);b1=1./a1;%加速度倒数b2=1./a2;b3=1./a3;b4=1./a4;b5=1./a5;figure(3)subplot(2,3,i)plot(u1,b1,u2,b2,u3,b3,u4,b4,u5,b5);title('加速度倒数-速度曲线');xlabel('u')ylabel('1/a')gtext('1/a1')gtext('1/a2')gtext('1/a3')gtext('1/a4')gtext('1/a5')gtext(['i0='num2str(i0(i))]);%计算采用2挡起步加速至96.6km/h的原地起步加速时间u1min=min(u1);u1max=max(u1);u2min=u1max;u2min=min(u2);u2max=max(u2);u3min=u2max;u3max=max(u3);u4min=u3max;u4max=max(u4);u5min=u4max;u5max=96.6;x1=[];x2=[];x3=[];x4=[];x5=[];y=3401;for j=1:3401if u3(j)<=u3minx1=[j];endendq1=max(x1);ua3=u3(q1:y);a3=b3(q1:y);for k=1:3401if u4(k)<=u4minx2=[k];endendq2=max(x2);ua4=u4(q2:y);a4=b4(q2:y);for l=1:3401if u5(l)<=u5minx3=[l];elseif u5(l)<=u5maxx4=[l];endendq2=max(x2);q3=max(x3);q4=max(x4);ua5=u5(q3:q4);a5=b5(q3:q4);s2=trapz(b2,u2);%二挡运行时间s3=trapz(ua3,a3);%三挡运行时间s4=trapz(ua4,a4);%四挡运行时间s5=trapz(ua5,a5);%五挡运行时间s=[s2s3s4s5];t=sum(s)*1000/3600;%总时间disp('t=');disp(t);gtext(['t='num2str(t)'s']);Pe1=Ft1.*u1./3600;%各档对应转速下的功率Pe2=Ft2.*u2./3600;Pe3=Ft3.*u3./3600;Pe4=Ft4.*u4./3600;Pe5=Ft5.*u5./3600;endPemax=max([max(Pe1)max(Pe2)max(Pe3)max(Pe4)max(Pe4)]);%发动机最大功率disp('90％负荷');Pe=0.9*Pemax%90％负荷rou=0.7;Ua=0.377*r*3401./(i0.*ig4(4));B=[1233.9-84.4782.9788-0.0474490.00028230];%负荷特性曲线拟合公式系数，对应n=3403r/min for i=1:5b=B*[1Pe Pe^2Pe^3Pe^4]';%根据拟合公式计算比油耗Qs5(i)=Pe.*b./(1.02.*Ua(i).*rou.*g);%百公里油耗量enddisp('Qs5=');disp(Qs5);tt=[93.973058.283444.678445.366744.9793];figure(2)plot(Qs5,tt,'*')hold onxi=30:0.001:100;pp=interp1(Qs5,tt,xi,'cubic');plot(xi,pp)%title('燃油经济性-加速时间曲线');%xlabel('燃油经济性--等速百公里油耗(L/100km)');%ylabel('动力性--原地起步加速时间(s)');gtext('5.17')gtext('5.43')gtext('5.83')gtext('6.17')gtext('6.33')运行结果：t=48.4235t=46.3511t=41.5557t=33.3796t=30.556490％负荷Pe=55.5647Qs4=25.103126.365628.307829.958730.7356图1加速度倒数-速度曲线图2燃油经济性-加速时间曲线4.3%空载(no load)—1；满载(full load)—2m1=4080;m2=9290;hg1=0.845;hg2=1.170;L=3.950;a1=2.100;a2=2.950;b1=L-a1;b2=L-a2;beta=0.38;%利用附着系数与制动强度的关系曲线z=0:0.01:1;phi_f1=L*beta.*z./(b1+z*hg1);%前轮利用附着系数(空载)phi_r1=L*(1-beta).*z./(a1-z*hg1);%空载时后轮利用附着系数(空载)phi_f2=L*beta.*z./(b2+z*hg2);%前轮利用附着系数(满载)phi_r2=L*(1-beta).*z./(a2-z*hg2);%后轮利用附着系数(满载)phi=z;%最理想情况figure(1);plot(z,phi_f1,'g--',z,phi_f2,'g-',z,phi_r1,'b--',z,phi_r2,'b-',z,phi,'k-'); axis([0,1,0,1]);grid on;box off;axis square;title('利用附着系数与制动强度的关系曲线');xlabel('制动强度z/g');ylabel('利用附着系数{\phi}');%ECE法规要求界限hold on;z1=0.2:0.01:0.8;z2=0.15:0.01:0.3;z3=0.3:0.01:1;phi1=(z1+0.07)/0.85;phi2=z2-0.08;phi3=z2+0.08;phi4=(z3-0.02)/0.74;plot(z1,phi1,'r-.',z2,phi2,'r-.',z2,phi3,'r-.',z3,phi4,'r-.'); legend('{\phi}_{f}(空载)','{\phi}_{f}(满载)','{\phi}_{r}(空载)',...'{\phi}_{r}(满载)','{\phi}=z','ECE法规要求界限');legend('Location','northwest');%制动效率曲线Ef1=z./phi_f1*100;Er1=z./phi_r1*100;Ef2=z./phi_f2*100;Er2=z./phi_r2*100;figure(2);plot(phi_f1,Ef1,'k',phi_r1,Er1,'k',phi_f2,Ef2,'b',phi_r2,Er2,'b'); axis([0,1,0,100]);grid on;box off;axis square;title('制动效率曲线');xlabel('附着系数{\phi}');ylabel('制动效率(%)');text(0.25,90,'E_f');text(0.6,88,'满载');text(0.8,90,'E_r');text(0.7,65,'空载');text(0.72,72,'E_r');程序运行结果如下：图1图2 5.11L=3.048;K=0.0024;i=20;ua=0:0.01:120;f=1/L.*(1./(3.6./ua+K.*ua./3.6));syms uF(u)=(1/L.*(1./(1./u+K.*u)))./i;plot(ua,f,'b-');xlabel('u_a(km/h)');ylabel('\omega_r/\delta)_s');F(22.35)运行结果：ans=0.1667386094155036.2f=0.1:0.1:100;Gq1=5.12*10^(-5)./f.^2;Gq2=2.02*10^(-3).*f./f;Gq3=7.98*10^(-2).*f.^2;figure(1);loglog(f,Gq1,'b-');box off;axis square;title('位移功率谱密度'); xlabel('f/Hz');ylabel('G_q(f)/(m^2s)');figure(2);loglog(f,Gq2,'b-');box off;axis square;title('速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''(f)/(m^2/s)');figure(3);loglog(f,Gq3,'b-');box off;axis square;title('加速度功率谱密度'); xlabel('f/Hz');ylabel('G_q_''_''(f)/(m^2/s^3)');运行结果如下：6.4f1=0.1:0.01:1;f2=1:0.01:10;G1=0.010106.*f1.^2;G2=0.010106.*f2./f2;loglog(f1,G1,'b-');hold onloglog(f2,G2,'b-');title('车身加速度的功率谱密度'); xlabel('f/Hz');ylabel('G_z_''_''(f)/m^-^1'); axis([0.1,10,0.0001,0.1]);运行结果如下：6.51.计算说明1）①幅频特性()122222114z q λζλγ⎡⎤-+⎢=⎢⎥∆⎣⎦，其中0ωλω=，()2222222111141λλγζλγλμμ⎡⎤⎡⎤⎛⎫⎛⎫∆=-+--+-⎢⎥⎢⎥ ⎪ ⎪⎝⎭⎝⎭⎣⎦⎣⎦；()12222222211414z z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦；()1222222221414s s s s s p z ζλλζλ⎡⎤+⎢=⎢⎥-+⎣⎦，其中s s ωλω=；②均方根谱()()221~j 2z q H f f ωππ=；()()222~j 2z qH f f ωππ== ；()()2221~1j 2p qz z H f f z qωππ== ；③其他值()0.5360q q G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536110z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.536220z z G f df σ⎡⎤=⎢⎥⎣⎦⎰ ，()0.5360a a G f df σ⎡⎤=⎢⎥⎣⎦⎰()()()()()()13622w 012222412.5362202412.5d 12.5 0.5d d 1d d 4a a a a a a W f G f f f G f f G f f G f f G f f f ⎡⎤=⎢⎥⎣⎦⎡⎤⎛⎫⎛⎫=+++⎢⎥ ⎪ ⎪⎝⎭⎢⎥⎝⎭⎣⎦⎰⎰⎰⎰⎰其中，()()()()()0.5 0.52 2441 412.512.512.536f f f W f f f f <<⎧⎪⎪<<⎪=⎨<<⎪⎪<<⎪⎩()aw w 020lg L a =2）见1）中计算式3）()()()22004q q q G f G G n n uωπ== ()()()1122363622d ~00d d dd f q q f q f H f G f f G f f q σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122d 12f q f γλπ⎡⎤=⎢⎥∆⎣⎦ ()()()1122363622d //~00d d dd F Gq q F G q F H f G f f G f f Gq σ⎡⎤⎡⎤==⎢⎥⎢⎥⎣⎦⎢⎥⎣⎦⎰⎰ ，其中122222d 1412F f Gq g λζλμπγ⎡⎤⎛⎫-+⎢⎥⎪+⎢⎥⎝⎭=⎢⎥∆⎢⎥⎢⎥⎣⎦其余见1）中计算式2.程序清单1）f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi*fs;lambda=omega./omega0;lambda_s=omega./omega_s;%计算并绘制幅频特性delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;figure(1)loglog(f,z1_q,'b-',0:0.1:100,(0:0.1:100)./(0:0.1:100),'g-');axis([0.1,100,0.1,10]);grid on ;title('z_1~q 的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_1}{q}|$$','Interpreter','latex');gtext('0:1');figure(2)loglog(f,z2_z1,'b-',0:0.1:1,(0:0.1:1)./(0:0.1:1),'g-',1:0.1:100,1./...(1:0.1:100),'g-');axis([0.1,100,0.01,10]);grid on;title('z_2~z_1的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{z_2}{z_1}|$$','Interpreter','latex'); gtext('0:1');gtext('-1:1');figure(3)loglog(f,p_z2,'b-');axis([0.1,100,0.01,10]);grid on;title('p~z_2的幅频特性');xlabel('激振频率f/Hz');ylabel('$$|\frac{p}{z_2}|$$','Interpreter','latex');%计算并绘制均方根谱sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q; figure(4)loglog(f,sqrt_Gz1,'b-');grid on;title('车轮位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''1}(f)}$','Interpreter','latex'); figure(5)loglog(f,sqrt_Gz2,'b-');grid on;title('车身位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_{z''''2}(f)}$','Interpreter','latex'); figure(6)loglog(f,sqrt_Ga,'b-');grid on;title('传至人体的位移均方根谱');xlabel('激振频率f/Hz');ylabel('$\sqrt{G_a(f)}$','Interpreter','latex');%计算其它值sigma_q=sqrt(trapz(f,sqrt_Gq.^2));%路面不平度加速度均方根值sigma_z1=sqrt(trapz(f,sqrt_Gz1.^2));%车轮加速度均方根值sigma_z2=sqrt(trapz(f,sqrt_Gz2.^2));%车身加速度均方根值sigma_a=sqrt(trapz(f,sqrt_Ga.^2));%传至人体的加速度均方根值for i=1:Nif f(i)<=2W(i)=0.5;elseif f(i)<=4W(i)=f(i)/4;elseif f(i)<=12.5W(i)=1;elseW(i)=12.5/f(i);endendaw=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));%加权加速度均方根值a0=10^(-6);Law=20*log10(aw/a0);%加权振级format shortdisp('路面不平度加速度均方根值=');disp(sigma_q);disp('车轮加速度均方根值=');disp(sigma_z1);disp('车身加速度均方根值=');disp(sigma_z2);disp('传至人体的加速度均方根值=');disp(sigma_a);disp('加权加速度均方根值=');disp(aw);disp('加权振级=');disp(Law);2）%随fs变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=1.5:0.01:6;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;lambda=omega./omega0;for i=1:length(fs)omega_s=2*pi.*fs(i);lambda_s=omega./omega_s;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(fs,aw);grid ontitle('a_w随f_s的变化曲线');xlabel('f_s/Hz');ylabel('a_w/m·s^-^2');figure(2)plot(fs,Law);grid ontitle('L_aw随f_s的变化曲线');xlabel('f_s/Hz');ylabel('L_a_w/dB');%随zeta_s变化f0=1.5;zeta=0.25;gamma=9;mu=10;fs=3;zeta_ss=0.125:0.001:0.5;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega0=2*pi*f0;omega_s=2*pi.*fs;lambda=omega./omega0;lambda_s=omega./omega_s;for i=1:length(zeta_ss)zeta_s=zeta_ss(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;p_z2=((1+(2*zeta_s.*lambda_s).^2)./((1-lambda_s.^2).^2+(2*zeta_s.*...lambda_s).^2)).^0.5;sqrt_Gq=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u);sqrt_Gz1=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z1_q;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sqrt_Ga=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*p_z2.*z2_z1.*z1_q;for j=1:Nif f(j)<=2W(j)=0.5;elseif f(j)<=4W(j)=f(j)/4;elseif f(j)<=12.5W(j)=1;elseW(j)=12.5/f(j);endendaw(i)=sqrt(trapz(f,W.^2.*sqrt_Ga.^2));enda0=10^(-6);Law=20*log10(aw/a0);figure(1)plot(zeta_ss,aw);grid ontitle('a_w随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('a_w/m·s^-^2');figure(2)plot(zeta_ss,Law);grid ontitle('L_aw随\zeta_s的变化曲线');xlabel('\zeta_s');ylabel('L_a_w/dB');3）%随f0变化f0=0.25:0.01:3;zeta=0.25;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;lambda_s=omega./omega_s;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(f0)omega0=2*pi*f0(i);lambda=omega./omega0;delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(f0==1.5);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(f0,sm_z2,f0,sm_fd,f0,sm_FdG);axis([0.25,3,-25,10]);grid on;title('各响应量均方根值随f_0变化的曲线');xlabel('f_0/Hz');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随zeta变化f0=1.5;zeta0=0.125:0.001:0.5;gamma=9;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(zeta0)zeta=zeta0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(zeta0==0.25);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(zeta0,sm_z2,zeta0,sm_fd,zeta0,sm_FdG);axis([0.125,0.5,-4,4]);grid on;title('各响应量均方根值随\zeta变化的曲线');xlabel('\zeta');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');%随gamma变化f0=1.5;zeta=0.25;gamma0=4.5:0.1:18;mu=10;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(gamma0)gamma=gamma0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(gamma0==9);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(gamma0,sm_z2,gamma0,sm_fd,gamma0,sm_FdG);axis([4.5,18,-6,6]);grid on;title('各响应量均方根值随\gamma变化的曲线');xlabel('\gamma');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB'); legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');legend('Location','southeast');%随mu变化f0=1.5;zeta=0.25;gamma=9;mu0=5:0.1:20;fs=3;zeta_s=0.25;u=20;Gq_n0=2.56*10^(-8);n0=0.1;delta_f=0.2;N=180;f=delta_f*[1:N];omega=2*pi*f;omega_s=2*pi*fs;omega0=2*pi*f0;lambda_s=omega./omega_s;lambda=omega./omega0;Gq_f=4*pi^2*Gq_n0*n0^2*u;g=9.8;for i=1:length(mu0)mu=mu0(i);delta=((1-lambda.^2).*(1+gamma-1./mu.*lambda.^2)-1).^2+4*zeta^2.*lambda.^...2.*(gamma-(1./mu+1).*lambda.^2).^2;z1_q=gamma.*(((1-lambda.^2).^2+4*zeta.^2.*lambda.^2)./delta).^0.5;z2_z1=((1+4*zeta.^2.*lambda.^2)./((1-lambda.^2).^2+4*zeta.^2.*lambda.^2))....^0.5;sqrt_Gz2=4*pi^2.*f.*sqrt(Gq_n0*n0^2*u).*z2_z1.*z1_q;sigma_z2(i)=sqrt(trapz(f,sqrt_Gz2.^2));fd_q=gamma*lambda.^2./(2*pi.*f).*(1./delta).^0.5;sigma_fd(i)=sqrt(trapz(f,Gq_f.*fd_q.^2));Fd_Gq=2*pi.*f.*gamma./g.*(((lambda.^2./(1+mu)-1).^2+4*zeta.^2*lambda.^2)..../delta).^0.5;sigma_FdG(i)=sqrt(trapz(f,Gq_f.*Fd_Gq.^2));endm=find(mu0==10);sgm_z2=sigma_z2(m);sgm_fd=sigma_fd(m);sgm_FdG=sigma_FdG(m);sm_z2=20*log10(sigma_z2/sgm_z2);sm_fd=20*log10(sigma_fd/sgm_fd);sm_FdG=20*log10(sigma_FdG/sgm_FdG);plot(mu0,sm_z2,mu0,sm_fd,mu0,sm_FdG);axis([5,20,-2,2]);grid on;title('各响应量均方根值随\mu变化的曲线');xlabel('\mu');ylabel('\sigma_z_''_''_2,\sigma_f_d,\sigma_F_d_/_G/dB');legend('\sigma_z_''_''_2','\sigma_f_d','\sigma_F_d_/_G');3.结果分析1）路面不平度加速度均方根值=0.3523车轮加速度均方根值=0.2391车身加速度均方根值=0.0168传至人体的加速度均方根值=0.0161加权加速度均方根值=0.0100加权振级=80.0287分析：根据课本中表6-2，a w=0.0100<0.315，L aw=80.0287<110，故乘客没有不舒适。

《汽车理论》课程设计说明书-第1章

驶
F =F +F +F
t w i
j
知当加速阻力 F =0 时，可求出最大爬坡度。 i =arcsin[ F -（ F + F ）]/( m g)
i
max
t
j
w
0
F
= t
T i i ηt
tq g 5 0
r
CDA *U a ^ 2 21.15
F
F
j
= w
=
mgf
0
由上式
联合求出 i
max
3.4 3.4 汽车外特性图
F =F
i
j
=0,判断汽车在最高档
F
t
与
F +F
f a max
w
两
u
。
u
a max
。
⑵超车加速时间（50～70km/h）
它是指汽车从静止状态下，由第一挡起步，并以最大的加速强度（包括节气门全开和选择最恰当的换挡时机）逐步换至高挡后，达到某一预定的车速或距离所需要的时间。这里需要注意换挡点的确定：在加速度倒数曲线中，根据动力性换挡的原则：若无交点，尽可能的用低档行驶；若有交点，则在交点处换挡。
湖北汽车工业学院
Hubei Automotive Industries Institute
课程设计说明书
课程名称设计题目班号专业学生姓名指导教师( 指导教师(签字) 签字) 成绩学号
起止日期
年
月
日 --
年_
月
日
1 / 19
目录
一、设计任务及要求 …………………………………………………………3 1.1 课程设计的目的 ………………………………………………………3 1.2 课程设计的任务及思路 ………………………………………………3 二、界面及功能实现三、理论模型的建立 3.1 最高车速 u a max …………………………………………………………5 …………………………………………………………7 ………………………………………………………7

汽车理论课程设计

汽车理论课程设计姓名学号专业班级指导教师年月日汽车理论课程设计任务书本次计算设计以小组为单位进行，每组计算两种车型(大型车、小型车)。

先进行手工计算并绘图（每个档位选取5-8个特征点），然后计算机编程（MATLAB）实现并绘图，最终打印计算说明书和程序。

具体设计要求如下：1.汽车动力性、经济性分析计算通过查阅收集有关资料，计算分析给定型号汽车的动力性能及燃油经济性，并绘出该车型的发动机外特性曲线，驱动力——行驶阻力平衡图，动力特性图，百公里油耗曲线。

根据计算结果，结合实际情况，分析该车型发动机参数和底盘性能参数匹配是否合理，并提出修改意见。

2.参数部分车型参数不完整，请查阅相关资料或用经验公式计算选取，并经手动计算分析后修正获得。

（说明书中注明出处）3.说明书说明书内容包括：任务书、目录、各车型参数、参数分析、计算、图表、结论、设计总结及心得体会。

车型一：解放CA1091载货汽车一、发动机CA6102（附表一）Nmax=99kw(相应转速3000r/min)Mmax=373N.m(相应转速1300r/min)二、整车参数：1.尺寸参数：全长L=7205mm，全宽B=2476mm，全高H=2436mm，轴距L1=4050mm,前轮距B1=1850mm，后轮距B2=1740mm.2.重量参数（附表二）3.性能参数：变速箱传动比i1=7.64，i2=4.835，i3=2.856，i4=1.895，i5=1.377，i6=1，i 倒=7.66。

主减速器比i0=6.33。

车轮：9.00-20。

三、使用数据：滚动阻力系数f=0.03；道路阻力系数：强度计算用Φ=1性能计算用Φ=0.8空气阻力系数：Cd=0.8；迎风面积：A=0.78×宽×高；最大速度：Vmax=90km/h；最大爬坡度：28%；传动系效率：η=0.9表二：重量参数车型四：SATANA2000轿车一、发动机2VQS（附表一）Nmax=74kw(相应转速5200r/min)Mmax=155N.m(相应转速3800r/min)二、整车参数：1.尺寸参数：全长L=4680mm，全宽B=1700mm，全高H=1423mm，轴距L1=2656mm,前轮距B1=1414mm，后轮距B2=1422mm.2.重量参数（附表二）3.性能参数：变速箱传动比i1=3.45，i2=1.94，i3=1.29，i4=0.985，i5=0.8,主减速器比i0=4.1。

MATLAB软件在汽车理论教学中的应用

1.前言《汽车理论》是汽车类专业的一门必修专业课，其特点是理论性强，计算和分析工作多，教和学两方面都感觉难度较大。

现有的实验和实习效果有限，为学生提供的主要是感性的认识和对书本知识的验证，对抽象概念和复杂公式无能为力。

MATLAB是一个可视化的计算软件，功能强大、使用简单，被广泛地应用于科学和工程计算领域，包括汽车的设计和分析。

各大汽车公司、零部件企业的研究院或技术中心，在研究和分析工作中大量应用MATLAB实现。

在《汽车理论》教学中使用MATLAB为工具，利用其强大的计算功能和图形功能，可以方便地完成各种性能的计算；同时，利用MATLAB 的数值计算函数和Simulink软件，可以对《汽车理论》中复杂的过程进行仿真分析和求解。

这些计算和分析的结果都可以通过MATLAB提供的可视化手段呈现给学生，有助于清晰地阐释抽象的概念，并系统地仿真复杂的分析过程。

2.易学易用的计算工具MATLAB被称为“演草纸”式的计算机语言，在一般的科学计算中完全可以替代传统的计算机语言如C、Fortran、Basic等，矩阵运算更是MATLAB的拿手好戏。

MATLAB语法简单、界面友好、使用方便，既不用苦于指针等晦涩概念的掌握，也不必与难以操作的DOS环境打交道。

《汽车理论》前四章介绍了汽车动力性、经济性以及制动性的计算和校核，工作量较大，有时候还要不断的调整参数并反复计算。

这些性能的计算对于任何一种车型的设计和开发都是必不可少且极为重要的。

利用MATLAB编制程序，其简单的语法，能够使设计人员将精力集中于核心工作。

3.方便实用的图形功能《汽车理论》课程中需要绘制的曲线图形非常多，“驱动力-行驶阻力图”、“功率平衡图”、“动力因素图”、“百公里油耗曲线”、“C曲线”、“制动效率曲线”等等，数不胜数。

传统的程序设计语言如C、Fortran、Basic等，绘图功能的使用极为不便，程序员需要了解计算机的图形模式等复杂的软件知识，在绘图前还要做大量的设置和转换工作。

融入Matlab和Carsim联合仿真的汽车理论教学方法探讨

AUTOMOBILE EDUCATION | 汽车教育融入Matlab和Carsim联合仿真的汽车理论教学方法探讨程准南京林业大学江苏省南京市 210037摘要：传统的《汽车理论》课程教学方法以讲授法为主，而课程具有一定的理论深度且重难点多，学生在缺乏对汽车实体、系统动力学状态、测量信号等直接感知的情况下相关知识和能力的增进效果有限。

本文围绕融入Matlab和Carsim联合仿真的《汽车理论》课程新型教学方法进行了分析和探讨。

以汽车的制动性为例进行了教学设计，通过Matlab调节汽车参数及观测信号，通过Carsim观察汽车运行状态，以直接感知、实际训练和自主探究等多方面出发培养学生工程综合能力。

关键词：汽车理论　联合仿真　教学方法　直接感知1　引言《汽车理论》主要涉及与汽车动力学有关的汽车使用性能（主要包括动力性、燃油经济性、制动性、操纵稳定性、平顺性和通过性）[1]，是车辆相关专业的核心课程。

其各章节的知识内容包含了一定数量的动力学方程，知识点表征时多借助简图进行受力分析并且文字性描述较多。

传统的教学方法以课堂讲授为主并辅之以课后练习（主要采用习题的方式）。

学生学习时面对纯粹的理论公式、数据图表、受力简图和文字分析难以提高对知识的理解和应用程度。

在工程教育专业认证的背景下《汽车理论》是车辆相关专业“毕业要求”的关键课程支撑[2]。

这使得《汽车理论》课程有了更高的目标[3]，例如学生通过对这门课程需掌握车辆工程专业知识，并将相关知识和数学模型方法用于车辆工程专业问题解决方案的比较与综合；学生能基于相关科学原理和数学建模方法正确表达车辆工程领域的复杂工程问题；学生能够对实验结果进行分析和解释，通过信息综合解决复杂工程问题，并得到合理有效的结论等。

借助于仿真技术模拟再现汽车运动过程能够显著提升教学效果。

文献[4]阐述了ADAMS软件在《汽车理论》（主要介绍了汽车的制动性、操纵稳定性和平顺性）教学过程中的应用。