汽车理论MATLAB程序

武汉理工大学

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选 其中的一种进行整车性能计算）：

1）绘制汽车驱动力与行驶阻力平衡图。

2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。

3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h的加速时间。

轻型货车的有关数据：

汽油发动机使用外特性的Tq-n曲线的拟合公式为

nn2n3n4Tq??19.313?295.27()?165.44()?40.874()?3.8445()

1000100010001000式中，Tq为发动机转矩（N?m）;n为发动机转速（r/min）。 发动机的最低转速nmin=600r/min,最高转速nmax=4000r/min。 装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径

%% 汽车动力性计算（自己编的动力性计算程序，供大家计算动力性时参考，具体参数大家根据所给程序对应输入，并对坐标轴数值按需要进行修改） clc; clear; close all;

%%根据所给发动机数据拟合外特性曲线（发动机数据按照你所得到的数据进行输入）

n_test=[500 600 700 800 900 1000 1100 1200 1300 1400 1500 1600 1700 1800 1900 2000 2100 2200];

T_test=[975 1108 1298 1496 1546 1620 1670 1785 1974 1974 1970 1889 1829 1748 1669 1700 1524 1105]; figure(1)

plot(n_test,T_test,'g'); hold on grid on

%p=polyfit(n_test,T_test,7); p=polyfit(n_test,T_test,2); n=[450:1:2200]; Ttq=polyval(p,n); plot(n,Ttq,'k');

xlabel('发动机转速n(r/min)'); ylabel('

设计目的：车牌定位系统的目的在于正确获取整个图像中车标的区域，并识别出车标。 程序效果：

程序实现：

STEP1：输入待处理的原始图像: 程序：

1 2 3 4

clear;clc;close all;

%Step1 获取图像 装入待处理彩色图像并显示原始图像 Scolor = imread('1.jpg');%imread函数读取图像文件 subplot(3,4,1);imshow(Scolor),title('原始图像')

输出：

原始图像

SETP2：图像的灰度化:

彩色图像包含着大量的颜色信息，不但在存储上开销很大，而且在处理上也会降低系统的执行速度，因此在对图像进行识别等处理中经常将彩色图像转变为灰度图像，以加快处理速度。由彩色转换为灰度的过程叫做灰度化处理。选择的标准是经过灰度变换后，像素的动态范围增加，图像的对比度扩展，使图像变得更加清晰、细腻、容易识别。 程序：

1 %将彩色图像转换为黑白并显示

2 Sgray = rgb2gray(Scolor);%rgb2gray转换成灰度图

3 subplot(3,4,2);imshow(Sgray),title('原始黑白图像');

输出：

原始黑白图像

STEP3：对原始图像进行开操

汽车传动系统，多目标优化（粒子群法）matlab原程序代码，供参考学习！可以根据实际情况进行更改进行运算。原码：

function apso

% 参数设置定义全局变量

global lamda1 lamda2 m ua_max eta_T r G f alpha Cd A rou K Ttq_max Fz fai ge_ne_pe du

lamda1 = 0.2; % 动力性发挥程度加权因子； lamda2 = 0.8; % 经济性加权因子； m = 1092; % 整车质量（kg）； ua_max = 50; % 最大车速(km/h)； eta_T = 0.9; % 传动系的传动效率； r = 0.3; % 车轮半径(m)；

g = 9.8; % 重力加速度（g*m/s^2） G = m*g; % 汽车重力G=mg，(N)； f = 0.015; % 汽车的滚动阻力系数； alpha = 25*pi/180; % 道路坡度角-->弧度； Cd = 0.32; % 空气阻力系数；

A = 1.5; % 迎风面积，即汽车行驶方向的投影

汽车理论中期作业

孙野 20081268

1.3（1）绘制汽车驱动力与行驶阻力平衡图 选用5挡变速器进行整车性能计算

发动机转速与汽车行驶速度之间的关系：ua?0.377发动机使用外特性的Tq-n曲线拟合公式：

Tq??19.313?295.27(n1000)?165.44(n1000)?40.874(2rnigi0

n1000)?3.8445(3n1000)

4汽车驱动力：Ft?Tqigi0?Tr

CDA21.152行驶阻力：Ff?Fi?Fw?Fj?Gf?Gi?

用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; for i=1:69

ua??mdudt?mgf?CDA21.15ua

2n(i)=(i+11)*50;

Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)^2+

2009 届 成 绩

海南大学机电工程学院

汽车工程系

汽 车 理 论 课 程 设 计

题 目： 汽车动力性的仿真 学 院： 机 电 工 程 学 院 专 业： 09级交通运输 姓 名： 黄 生 锐 学 号： 20090504 指导教师：

编号 名 称 课程设计论文件数 页数 编号 名 称 Matlab编程源程序 件数 1 页数 1 2

设计任务书 1 1 3 2012年6月20日

汽车理论课程设计任务书

姓名 黄生锐 学号 20090504 汽车动力性的仿真 专业 09交通运输 课程设计题目 内容摘要： 本设计的任务是对一台Passat 1.8T手动标准型汽车的动力性能进行仿真。采用MATLAB编程仿真其性能，其优点是：一是能过降低实际成本，提高效率；二是获得较好的参数模拟，对汽车动力性能提供理论依据。 主要任务： 根据该车的外形、轮距、轴距、最小离地

1. 用MATLAB命令产生如下信号，并绘出波形图

(1) 2e?1.5tu(t)

t=-2*pi:pi/50:2*pi; a=t>=0;

f(a)=2*exp(-1.5*t(a)); f(~a)=0; plot(t,f)

（2）2sin(2?t??4)

t=-2*pi:pi/50:2*pi; y=2*sin(2*pi*t+pi/4); plot(t,y)

(3) f(t)???2?00?t?1

t?0,t?1t=-5:0.01:5;

y=2*rectpuls(t-0.5,1); plot(t,y) axis equal

2. 已知信号f(t)的波形如图所示，试用MATLAB命令画出f(t?2)、

f(3t)、f(?t)、f(?3t?2)的波形图。

f(t)1?201t

function y=p(t) y1=(t<0&t>=-2);

y2=(-t+1).*(t<=1&t>=0); y=y1+y2

t=-5:0.01:5; y=p(t-2); plot(t,y) axis equal

t=-5:0.01:5; y=p(3*t); plot(t,y) axis equal

t=-5:0.01:5; y=p(-t); plot(t,y) axis equa

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

汽车理论大作业

同济大学汽车学院

xxx xxxxxx 2013.12

%%0. 作业要求

%%1. 发动机外特性曲线

nmin=600; nmax=4000;

n=nmin:0.01:nmax; m=n/1000;

Tq=-19.313+295.27*m-165.44*m.^2+40.874*m.^3-3.8445*m.^4; figure; plot(n,Tq);

title('发动机外特性曲线 Tq-n'); xlabel('n/(r/min)'); ylabel('Tq/Nm');

Pe=(Tq.*n)/9550; figure; plot(n,Pe);

title('发动机外特性曲线 Pe-n'); xlabel('n/(r/min)'); ylabel('Pe/kw');

%%2. 驱动力曲线

nmin=600; nmax=4000;

n=nmin:0.01:nmax; m=n/1000;

Tq=-19.313+295.27*m-165.44*m.^2+40.874*m.^3-3.8445*m.^4; r=0.367; elta=0.85; i0=5.83;

ig_4=[6.09 3.09 1.71 1.00];

ig_

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