MATLAB和Origin绘制发动机万有特性曲线程序及方法

合集下载

基于MATLAB语言的发动机特性研究

$AB ) C %,-
) ?++ $"96; -D-")+)". ))?"? )E>") .+?H. .)+"F !"9@ ﹒
$AB ) C %,-
- )-? - .F- /+> - ?)+ ) ++D ) -?/ ) .F> ) D?? ) E?>
)>)H> ).+HD ))/H. )).HF ))+H/ ))+HE ))+HF ))DHF ).+H-
* 本文利用()*+),强大的数据处理和三
维曲线绘图功能 % 提出了一种绘制发动机特性的新
"% #
!
发动机特性数学模型
对发动机特性及其与车辆动力传动系匹配进行
式中 %) 和 , 分别表示任意一条特性曲线的纵坐标与横坐标 /* 为多项式的阶次 /’+为多项式各项拟合系数 % 一般采用最小二乘法确定 *
根据以上原理 ! 用 !"#$%& 语言编制相应的程序解出回归方程的系数向量 ! 即可得到试验数据的回归方程 "
&+ & # $ &( ’ !"’ ! !$ ’ ( !#"!$"!#")!#"$"!$" *!#!$" % " * ’ $ &’,- ( %
- >++ $"96; !"9@ ﹒
$AB ) C%
()*+), 中的曲线拟合等相应的工具箱或 ./01234%./01560 %./014//0%./017/82 等函数可以非常好

基于MATLAB语言的发动机万有特性曲线的绘制

ＭＡＴＬＡＢ是世界上最广泛的工程计算语言，利用其强
Plotting of Engine Univer sal Char acter istics Cur ve Based on MATLAB
HUANG Meimei, ZHAO Zhiwei，JIN Hualei, JIA Yantao,SUN Haipeng
低。其实质是以二维的图形方式表达三维的信息，不直观
并且难以保证对数据进行深刻分析。
随着ＭＡＴＬＡＢ语言的广泛应用，因为其强大的数据处
理和三维曲线绘图功能，可进行工程计算、建模仿真和数
据分析处理等。本文则利用ＭＡＴＬＡＢ强大的功能，提出了
[3] 易雪梅，吴伶 . 用 MATLAB 语言绘制发动机万有特性的两种方法 [J]. 北京汽车, 2005, 5: 33- 35
[4] 李金辉，徐立友 . 基于 MATLAB 语言的发动机特性研究 [J]. 汽车科技, 2005, 3: 40- 42
参考文献
[1] 杨丽娟，赵丹平 . 基于 MATLAB 基础上的发动机万有特性曲线的建立 [J]. 汽车节能, 2010, 1: 32- 33
[2] 薛定宇，陈阳泉 . 基于 MATLAB/Simulink 的系统仿真与应用 [M]. 北京: 清华大学出版社 .2002
Key words: straight welded pipe, burr, broach, hydraulic system
!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

MATLAB和Origin绘制发动机万有特性曲线程序及方法

MATLAB及Origin绘制发动机万有特性曲线实例1.MATLAB绘制发动机万有特性实例（含程序代码）2.Origin绘制发动机万有特性曲线方法（操作步骤）一、MATLAB方法MTALB绘制发动机万有特性曲线需要建立M文件，并在M文件中将程序写入，运行即可，当然也可以通过调用excel数据来绘制万有特性曲线，下面是MATLAB绘制实例，可借鉴修改然后运行。

MATLAB程序如下：clcclear all%不同转速下的燃油消耗率与扭矩的曲线拟合be1=[222.8,220.4,232.4,228.5,227.8,232.6,248.5,245.9,272.4,329.7];Ttq1=[399.8,354.1,318.5,278.1,236.2,203.6,185.3,157.2,117.2,80.8];T1=80:320/9:400;%转换矩阵格式Be1=interp1(Ttq1,be1,T1,'spline');%n=1400r/min时燃油消耗率与扭矩的曲线拟合be2=[222.0,221.7,235.4,226.5,230.5,236.8,249.1,276.1,407.9,487.0];Ttq2=[409.1,365.7,328.3,284.1,243.7,203.2,164.3,123.9,83.5,39.7];T2=39:371/9:410;Be2=interp1(Ttq2,be2,T2,'spline');be3=[226.0,225.3,226.4,233.9,242.1,283.3,253.9,271.4,323.5,468.6];Ttq3=[408.3,368.3,328.3,289.0,244.4,208.8,167.7,132.1,89.5,46.1];T3=46:363/9:409;Be3=interp1(Ttq3,be3,T3,'spline');be4=[206.5,231.1,231.1,233.0,242.0,244.9,265.0,299.8,398.0,596.8];Ttq4=[425.6,380.3,332.7,290.9,244.4,205.1,160.2,114.5,68.8,30.7];T4=30:396/9:426;Be4=interp1(Ttq4,be4,T4,'spline');be5=[234.7,259.8,235.5,237.6,242.8,292.3,277.9,308.7,396.2,605.9];Ttq5=[420.7,379.6,334.6,291.6,244.4,202.8,157.5,116.0,74.1,37.8];T5=37:384/9:421;Be5=interp1(Ttq5,be5,T5,'spline');be6=[174.2,242.2,252.1,287.4,253.6,263.6,290.6,316.8,378.0,518.8];Ttq6=[404.6,360.5,322.7,283.0,243.3,205.5,162.1,124.7,86.8,52.4];T6=52:353/9:405;Be6=interp1(Ttq6,be6,T6,'spline');be7=[256.9,253.7,253.5,260.0,303.8,280.7,300.6,346.6,435.6,812.9];Ttq7=[378.0,344.7,310.3,264.3,226.1,186.8,154.2,115.3,76.3,34.1];T7=34:344/9:378;Be7=interp1(Ttq7,be7,T7,'spline');be8=[257.9,295.3,282.4,288.7,301.9,329.7,357.0,475.4,580.3,1080.1];Ttq8=[315.6,275.5,242.5,210.3,178.5,145.6,118.6,72.6,52.8,22.4];T8=22:294/9:316;Be8=interp1(Ttq8,be8,T8,'spline');B=[Be1';Be2';Be3';Be4';Be5';Be6';Be7';Be8'];N=[1400*ones(10,1);1600*ones(10,1);1800*ones(10,1);2000*ones(10,1);2200*ones(10,1);2400*ones (10,1);2600*ones(10,1);2800*ones(10,1)];Ttqn=[T1';T2';T3';T4';T5';T6';T7';T8'];G=[ones(80,1),N,Ttqn,N.^2,N.*Ttqn,Ttqn.^2];A=G\B;%A为6*1矩阵[n,Ttq]=meshgrid(1400:2800,100:600);%生成n-Ttq平面上的自变量“格点”矩阵be=A(1)+n.*A(2)++Ttq*A(3)+n.^2*A(4)+n.*Ttq*A(5)+Ttq.^2*A(6);Pe=Ttq.*n/9550;%外特性实验数据拟合Nw=[1403,1597,1797,1986,2102,2199,2303,2400,2507,2598,2700,2802]; Ttqw=[474,497,515,526,528.8,522.8,509.5,492.2,471.2,448.4,408.3,357.4]; n0=1400:2800;Ttqw_N=interp1(Nw,Ttqw,n0,'spline');h=repmat(Ttqw_N,501,1);ii=find(Ttq>h);%确定超出边界的“格点”下标be(ii)=NaN;%强制为非数Pe(ii)=NaN;%强制为非数%绘制等燃油消耗率曲线和等功率曲线三维拟合图subplot(1,2,1);mesh(n,Ttq,be);hold on;mesh(n,Ttq,Pe);axis([1000,3000,100,600,0,500]);hold on;xlabel('n(r/min)')ylabel('Ttq(N*m)')zlabel('Pe(KW) be(g/(KW*h))')title('等燃油消耗曲线和等功率曲线的三维拟合图')%绘制边界线(外特性曲线)subplot(1,2,2);plot(n0,Ttqw_N,'LineWidth',2);axis([1400,2800,100,550]);xlabel('n(r/min)');ylabel('Ttq(N*m)');title('万有特性曲线');hold on;%绘制等油耗率曲线的二维图B=contour(n,Ttq,be,11);%画等位线，并给出标识数据clabel(B);%把“等位值”沿等位线随机标识hold on;%绘制等功率曲线的二维图P=contour(n,Ttq,Pe,11);%画等位线，并给出标识数据clabel(P);%把“等位值”沿等位线随机标识legend('等油耗曲线','等功率曲线','外特性曲线')hold off二、Origin方法用origin软件绘制发动机万有特性曲线方法一、万有特性数据输入在excel中整理好发动机万有特性数据，主要包括发动机转速、扭矩、燃油消耗率及功率数据。

Matlab程序模拟汽车理论中的数据特性图1

车辆理论作业第三组第一大题确定一辆轻型货车的动力性能 1、绘制发动机的使用外特性曲线 n=600:1:4000 N=n/1000T=-19.313+295.27*n/1000-165.44*N.*N+40.874*N.*N.*N-3.8445*N.*N.*N .*NPe=T.*n/9550 %发动机的净功率 figure(1) plot(n,T) xlabel('转速n') ylabel('转矩T') figure(2) plot(n,Pe) xlabel('转速n') ylabel('功率Pe')5001000150020002500300035004000100110120130140150160170180转速n转矩T5001000150020002500300035004000010203040506070转速n功率P e2、绘制驱动力---行驶阻力平衡图 i0=6.17; %减速器的传动比 ig(1)=6.09 % 一档的传动比 ig(2)=3.09 ig(3)=1.71 ig(4)=1.00 nT=0.85 r=0.367 f=0.013 m1=3880 m2=9590 CdA=2.77 Ff1=m1*9.8*f Ff2=m2*9.8*fft1=T*ig(1)*i0*nT/r %一档所产生的驱动力 ft2=T*ig(2)*i0*nT/r ft3=T*ig(3)*i0*nT/r ft4=T*ig(4)*i0*nT/r Ft=[ft1;ft2;ft3;ft4]';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) ua=[ua1;ua2;ua3;ua4]';fw1=CdA*ua1.^2/21.15+Ff1 %一档的行驶阻力 fw2=CdA*ua2.^2/21.15+Ff1 fw3=CdA*ua3.^2/21.15+Ff1 fw4=CdA*ua4.^2/21.15+Ff1Fz1=[fw1;fw2;fw3;fw4]'; figure(3)plot(ua,Ft) %绘制各档的驱动力曲线hold onplot(ua,Fz1) %绘制行驶阻力曲线 hold onplot(ua4,Ff2) %绘制滚动阻力曲线 xlabel('ua/(km.h)')ylabel('F/N')title('汽车驱动力-行驶阻力平衡图')text(20,12000,'Ft1'); text(30,6000,'Ft2'); text(40,4000,'Ft3'); text(80,3000,'Ft4');text(90,2200,'（空载）Ff+fw'); text(105,1450,'满载Ff')102030405060708090200040006000800010000120001400016000ua/(km.h)F /N汽车驱动力-行驶阻力平衡图Ft1Ft2Ft3Ft4（空载）Ff+Fw满载Ff一档的驱动力很大，可用于爬坡且坡度较大；二档用于直接原地起步加速 3、绘制动力特性图D=(Ft-Fw)/(m1*9.8) %求动力因数f1=0.0076+0.000056*ua %一定车速条件下的滚动阻力 figure(4) plot(ua,D) hold onplot(ua,f1) text(20,0.3,'I') text(20,0.154,'II') text(40,0.1,'III') text(80,0.05,'IV') text(110,0.025,'f') title('汽车动力特性图') xlabel('ua/(km.h)') ylabel('D')1020304050607080900.050.10.150.20.250.30.350.40.45IIIIIIIVf汽车动力特性图ua/(km.h)DIf=0.218Iw1=1.798 Iw2=3.598die1=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(1).^2)/(m1*r^2);%一档的旋转质量换算系数die2=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(2).^2)/(m1*r^2); die3=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(3).^2)/(m1*r^2); die4=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(4).^2)/(m1*r^2); die=[die1;die2;die3;die4]';a1=9.8*D(:,1)/die1-9.8*f/die1; %求一档的加速度 a2=9.8*D(:,2)/die2-9.8*f/die2; a3=9.8*D(:,3)/die3-9.8*f/die3; a4=9.8*D(:,4)/die4-9.8*f/die4; aI=1./a1 %求一档的加速度倒数 aII=1./a2 aIII=1./a3 aIV=1./a4aa=[aI,aII,aIII,aIV]';figure(5) %绘制各档的加速度倒数曲线 plot(ua1,aI) hold on plot(ua2,aII) hold onplot(ua3,aIII) hold onplot(ua4,aIV)axis([5 99 0.01 10])10203040506070809012345678910aI aII aIIIaIV加速度倒数曲线ua/(km.h)1/a由图可知，加速度倒数没有交点，则可以在每档达到最大车速时换挡 nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=6.17;If=0.218; Iw1=1.798;Iw2=3.598;L=3.950;a=2.000;hg=0.9;m=3880;g=9.8; G=m*g; ig=[6.09 3.09 1.71 1.00 ]; nmin=600;nmax=4000; u1=0.377*r*nmin./ig/i0; u2=0.377*r*nmax./ig/i0;deta=0*ig; for i=1:4deta(i)=1+(Iw1+Iw2)/(m*r^2)+(If*(ig(i))^2*i0^2*nT)/(m*r^2); endua=[0: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:N k=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;else 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; enda=delta(1:k); t(i)=sum(a); endplot(t,ua,'b','LineWidth',2);axis([0 100 0 100]);title('四档汽车二档原地起步换挡加速时间曲线'); xlabel('时间t （s ）');ylabel('速度ua （km/h ）');1020304050607080901000102030405060708090100四档汽车二档原地起步换挡加速时间曲线时间t （s ）速度u a （k m /h ）第二大题n=600:1:4000 N=n/1000T=-19.313+295.27*n/1000-165.44*N.*N+40.874*N.*N.*N-3.8445*N.*N.*N .*NPe=T.*n/9550 figure(1) plot(n,T)text(1000,120,'转矩T') hold on plot(n,Pe)text(2000,20,'发动机净功率Pe')5001000150020002500300035004000020406080100120140160180转矩T发动机净功率Pei0=6.17; ig(1)=6.09 ig(2)=3.09 ig(3)=1.71 ig(4)=1.00 nT=0.85 r=0.367 f=0.013 m1=3880 m2=9590 CdA=2.77ua1=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); ua=10:1:100Pf=m1*9.8*f.*ua/(3600) Pw=CdA*ua.^3/(76140) figure(2) plot(ua1,Pe) hold on plot(ua2,Pe) hold on plot(ua3,Pe) hold onplot(ua4,Pe) hold onplot(ua,(Pf+Pw)/nT)010203040506070809010010203040506070评价：一二三档的后备功率都比较大，可在一定坡度的道路上直接启动加速；在最高档达到最高车速时，n0=[815 1207 1614 2012 2603 3006 3403 3884]B00=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7] B10=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291] B20=[72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113]B30=[-5.8629 -2.0533 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215]B40=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568] B0=spline(n0,B00,n) B1=spline(n0,B10,n) B2=spline(n0,B20,n) B3=spline(n0,B30,n) B4=spline(n0,B40,n)Pf3=m1*9.8*f.*ua3/(3600) Pw3=CdA*ua3.^3/(76140) Pf4=m1*9.8*f.*ua4/(3600) Pw4=CdA*ua4.^3/(76140) Ff1=m1*9.8*ffz3=CdA*ua3.^2/21.15+Ff1 fz4=CdA*ua4.^2/21.15+Ff1 %Pe3=fz3.*ua3./(3600*nT); %Pe4=fz4.*ua4./(3600*nT); Pe3=(Pf3+Pw3)/nTPe4=(Pf4+Pw4)/nTb3=B0+B1.*Pe3+B2.*Pe3.^2+B3.*Pe3.^3+B4.*Pe3.^4 b4=B0+B1.*Pe4+B2.*Pe4.^2+B3.*Pe4.^3+B4.*Pe4.^4Qs4=Pe4.*b4./(1.02*ua4*7.10*9.8) Qs3=Pe3.*b3./(1.02*ua3*7.10*9.8) figure(3)plot(ua4,Qs4) hold on plot(ua3,Qs3)01020304050607080901.41.61.822.22.42.62.8n=600:1:4000 N=n/1000T=-19.313+295.27*n/1000-165.44*N.*N+40.874*N.*N.*N-3.8445*N.*N.*N .*NPe=T.*n/9550 i0=6.17; ig(1)=6.09 ig(2)=3.09 ig(3)=1.71 ig(4)=1.00 nT=0.85 r=0.367 f=0.013 m1=3880 m2=9590 CdA=2.77ua1=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);n0=[815 1207 1614 2012 2603 3006 3403 3884]B00=[1326.8 1354.7 1284.4 1122.9 1141.0 1051.2 1233.9 1129.7] B10=[-416.46 -303.98 -189.75 -121.59 -98.893 -73.714 -84.478 -45.291]B20=[72.379 36.657 14.524 7.0035 4.4763 2.8593 2.9788 0.71113]B30=[-5.8629 -2.0533 -0.51184 -0.18517 -0.091077 -0.05138 -0.047449 -0.00075215]B40=[0.17768 0.043072 0.0068164 0.0018555 0.00068906 0.00035032 0.00028230 -0.000038568]B0=spline(n0,B00,n)B1=spline(n0,B10,n)B2=spline(n0,B20,n)B3=spline(n0,B30,n)B4=spline(n0,B40,n)Pf3=m1*9.8*f.*ua3/(3600) %三档的滚动阻力功率Pw3=CdA*ua3.^3/(76140) %三档的的空气阻力功率Pf4=m1*9.8*f.*ua4/(3600)Pw4=CdA*ua4.^3/(76140)Ff1=m1*9.8*fPe3=(Pf3+Pw3)/nT %三档的发动机净功率Pe4=(Pf4+Pw4)/nT %四档的发动机净功率b3=B0+B1.*Pe3+B2.*Pe3.^2+B3.*Pe3.^3+B4.*Pe3.^4b4=B0+B1.*Pe4+B2.*Pe4.^2+B3.*Pe4.^3+B4.*Pe4.^4Qs4=Pe4.*b4./(1.02*ua4*7.10) %四档的百公里油耗Qs3=Pe3.*b3./(1.02*ua3*7.10) %三档的百公里油耗ua=[25 40 50]ss=[50 250 250]fz=CdA*ua.^2/21.15+Ff1 %匀速条件下的阻力Pee=fz.*ua/(3600*nT);b=spline(b4,Pe4,Pee) %匀速条件下的燃油消耗率Qsy=Pee.*b./(1.02*ua*7.10)a=[0.2 0.2 -0.36]dt1=1/(3.6*0.2)dt2=1/(3.6*0.2)dt3=1/(3.6*(0.36))ua11=25:1:40ua22=40:1:50ua33=50:-1:25b11=spline(ua4,b4,ua11) %在加速条件下的各个速度节点的燃油消耗率b22=spline(ua4,b4,ua22)b33=spline(ua4,b4,ua33)If=0.218Iw1=1.798Iw2=3.598die4=1+(Iw1+Iw2)/(m1*r^2)+(If*i0^2*nT*ig(4).^2)/(m1*r^2);Pej1=(m1*9.8*f.*ua11/(3600)+CdA.*ua11.^3/(76140)+(die4*m1.*ua11/3600)*0.2)/nT %在三个加速条件下的各个速度节点Pej2=(m1*9.8*f.*ua22/(3600)+CdA.*ua22.^3/(76140)+(die4*m1.*ua22/3 600)*0.2)/nT %的发动机净功率Pej3=(m1*9.8*f.*ua33/(3600)+CdA.*ua33.^3/(76140)+(die4*m1.*ua33/3 600)*0.2)/nT %Qstj1=Pej1.*b11/(367.1*7.10) %在三中加速条件下的，各自的单位时间的燃油消耗量Qstj2=Pej2.*b22/(367.1*7.10)Qstj3=Pej3.*b33/(367.1*7.10)i11=size(Qstj1)i22=size(Qstj2)i33=size(Qstj3)i1=i11(2)i2=i22(2)i3=i33(2)Qt1=Qstj1(2:i1-1)Qt2=Qstj2(2:i2-1)Qt3=Qstj3(2:i3-1)q1=((Qstj1(1)+Qstj1(i1))*dt1./2+sum(Qt1)*dt1)/10q2=((Qstj2(1)+Qstj2(i2))*dt2./2+sum(Qt2)*dt2)/10q3=((Qstj3(1)+Qstj3(i3))*dt3./2+sum(Qt3)*dt3)/10Qall=(q1+q2+q3+Qsy(:,1)+Qsy(:,2)+Qsy(:,3))*100/1075q1 =5.8526q2 =4.9810 q3 =6.0186 Qall =2.7916第三大题(1)m0=4880; %空载质量m/kghg0=0.845; %空载质心高L0=3.95; %空载轴距a0=2.100; %空载质心至前轴距离m=9290; %满载质量m/kghg=1.170; %满载质心高L=3.95; %满载轴距a=2.950; %满载质心至前轴距离BB=0.38; %B为制动力分配系数g=9.8;G=m*g;b=L-a;G0=m0*g;b0=L0-a0;Fu1 =0:1000:100000;Fu2=1/2*(G/hg*(b^2+4*hg*L/G*Fu1).^0.5-(G*b/hg+2*Fu1));Fu02=1/2*(G0/hg0*(b0^2+4*hg0*L0/G0*Fu1).^0.5-(G0*b0/hg0+2*Fu1));F2=(1-B)*Fu1/B;plot(Fu1,Fu2) %绘制满载时的I曲线hold onplot(Fu1,Fu02) %绘制空载时的I曲线hold onplot(Fu1,F2); %%绘制beta曲线title('beta线和I曲线')text(30000,50000,'beta曲线')text(25000,10000,'空载时I曲线')text(25000,30000,'满载时I曲线')xlabel('Fu1')ylabel('Fu2')figure(2)fe=[0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9];for i=1:10Fxb1fk=0:100:120000;Fxb1rk=0:100:120000;Fxb2fk=(L0-fe(i)*hg0)/(fe(i)*hg0).*Fxb1fk-G0*b0/hg0;Fxb2rk=(-1*fe(i)*hg0)/(L0+fe(i)*hg0).*Fxb1rk+fe(i)*G0*a0/(L0+fe(i )*hg0);plot(Fxb1fk,Fxb2fk)hold onplot(Fxb1rk,Fxb2rk,'r')hold onendxlabel('Fxb1')ylabel('Fxb2')axis([0 120000 0 40000])title('空载时的f、r线组图')figure(3)for i=1:10Fxb1fk=0:100:220000;Fxb1rk=0:100:220000;Fxb2fk=(L-fe(i)*hg)/(fe(i)*hg).*Fxb1fk-G*b/hg;Fxb2rk=(-1*fe(i)*hg)/(L+fe(i)*hg).*Fxb1rk+fe(i)*G*a/(L+fe(i)*hg); plot(Fxb1fk,Fxb2fk)hold onplot(Fxb1rk,Fxb2rk,'r')hold onendxlabel('Fxb1')ylabel('Fxb2')axis([0 220000 0 80000])title('满载时的f、r线组图')012345678910x 104-2024681012141618x 104beta 线和I 曲线beta 曲线空载时I 曲线满载时I 曲线Fu1F u 2024681012x 10400.511.522.533.54x 104Fxb1F x b 2空载时的f 、r 线组图00.20.40.60.81 1.21.4 1.6 1.822.2x 105012345678x 104Fxb1F x b 2满载时的f 、r 线组图(2)-(3)m0=4880; %空载质量m/kghg0=0.845; %空载质心高L0=3.95; %空载轴距a0=2.100; %空载质心至前轴距离m=9290; %满载质量m/kghg=1.170; %满载质心高L=3.95; %满载轴距a=2.950; %满载质心至前轴距离BB=0.38; %B 为制动力分配系数g=9.8;G=m*g;b=L-a;G0=m0*g;b0=L0-a0;z=0:0.1:1; %制动强度fef=B*z*L./(b+z*hg);fef0=B*z*L0./(b0+z*hg0);fer=(1-B)*z*L./(a-z*hg);fer0=(1-B)*z*L0./(a0-z*hg0);figure(1)plot(z,fef,z,fef0,z,fer,z,fer0,z,z)text(0.8,1.5,'空车后轮')text(0.8,1.2,'满载后轮')text(0.7,0.6,'满载前轮')text(0.8,0.4,'空车前轮')figure(2)Ef=z./fef*100Ef0=z./fef0*100Er=z./fer*100Er0=z./fer0*100plot(fef,Ef,fer,Er)hold onplot(fef0,Ef0,fer0,Er0)axis([0 1 0 100])text(0.3,90,'Ef')text(0.8,90,'Er 满载')text(0.3,70,'Er 空载')00.10.20.30.40.50.60.70.80.9100.20.40.60.811.21.41.61.82空车后轮满载后轮满载前轮空车前轮00.10.20.30.40.50.60.70.80.910102030405060708090100Ef Er 满载Er 空载m0=4880;hg=0.845;l=3.950;a=2.100;b=1.850;beita=0.38;g=9.8;G0=m0*g ;m=9290;Hg=1.170;L=3.950;A=2.950;B=1.050;beita=0.38;g=9.8;G=m*g;t1=0.02;t2=0.02;ua0=30;fai=0.80;z=0:0.01:1.0;faifk=beita*z*l./(b+z*hg);faifm=beita*z*L./(B+z*Hg);fairk=(1-beita)*z*l./(a-z*hg);fairm=(1-beita)*z*L./(A-z*Hg);Efk=z./faifk*100;Efm=z./faifm*100;Erk=z./fairk*100;Erm=z./fairm*100;ak1=Erk(81)*g*fai/100;am1=Erm(81)*g*fai/100;Sk1=(t1+t2/2)*ua0/3.6+ua0^2/(25.92*ak1);Sm1=(t1+t2/2)*ua0/3.6+ua0^2/(25.92*am1);disp('空载时，汽车制动距离Sk1=');disp(Sk1);disp('满载时，汽车制动距离Sm1=');disp(Sm1);ak2=fai*g*a/(l+fai*hg);am2=fai*g*A/(L+fai*Hg);ak3=fai*g*b/(l-fai*hg);am3=fai*g*B/(L-fai*Hg);Sk2=(t1+t2/2)*ua0/3.6+ua0^2/(25.92*ak2);Sm2=(t1+t2/2)*ua0/3.6+ua0^2/(25.92*am2);Sk3=(t1+t2/2)*ua0/3.6+ua0^2/(25.92*ak3);Sm3=(t1+t2/2)*ua0/3.6+ua0^2/(25.92*am3);disp('空载时，前制动器损坏，汽车制动距离Sk2=');disp(Sk2);disp('满载时，前制动器损坏，汽车制动距离Sm2=');disp(Sm2);disp('空载时，后制动器损坏，汽车制动距离Sk3=');disp(Sk3);disp('满载时，后制动器损坏，汽车制动距离Sm3=');disp(Sm3);空载时，汽车制动距离Sk1=7.8668满载时，汽车制动距离Sm1=5.6354空载时，前制动器损坏，汽车制动距离Sk2=10.0061满载时，前制动器损坏，汽车制动距离Sm2=7.5854空载时，后制动器损坏，汽车制动距离Sk3=8.0879满载时，后制动器损坏，汽车制动距离Sm3=12.9629。

如何用matlab绘制电机效率map图或发动机万有特性曲线

如何用matlab绘制电机效率map图或发动机万有特性曲线前段时间写论文，需要绘制电机效率map图，其实和发动机万有特性曲线一样。

看了好多资料都不会，问问师兄也没具体画过。

困惑中查到貌似有几个软件可以画map图，由于我比较熟悉matlab，就选用它了，可是matlab也不知道咋画呀，我查看了matlab图形处理这一块，突然发现等高线图绘制，咦？？？这不就是高中地理学的吗？？？和map图万有特性图本质一样吗？？？就是contour函数啦，惊喜万分5.2.13 等值线图等值线图可用于绘制地理数据中的等高图、气象数据中的等势图等。

等值线图在二维图形中把第三维中相同大小的数据连接为等值线，一定程度上可以表示第三维的信息，同时等值线图相比三维图更容易观察数据之间的关系，被广泛的应用于各个领域。

MATLAB中提供了一系列的函数用于绘制不同形式的等高线图，其中包括：1．contour()函数contour()函数可用于绘制二维等值线图，函数的调用格式为：❑contour(z)：输入数据z为二维矩阵，绘制数据z的等值线，绘图时等值线的数量和数值根据矩阵z的数据范围自动确定。

❑contour(z,n)：绘制等值线图，设置等值线数目为n。

❑contour(z,v)：绘制等值线图，向量v设置等值线的数值。

❑contour(x,y,z)：绘制矩阵z的等值线图，输入参数x、y用于指定绘制的等值线图的坐标轴数据，同时输入数据x、y、z必须为大小相等的矩阵。

❑contour(x,y,z,n)：为指定坐标轴的等值线图设置等值线的数目n。

❑contour(x,y,z,v)：为指定坐标轴的等值线图设置等值线的数值v。

❑contour(...,LineSpec)：输入参数LineSpec用于设置等值线的线型。

❑[c,h] = contour(...)：返回contour()函数绘制的等高线图中的等值线的数值标签c和包含所有图形对象的句柄h；2．contourf()函数contourf()函数用于绘制带填充的二维等值线图。

基于MATLAB基础上的发动机万有特性曲线的建立

基于MATLAB基础上的发动机万有特性曲线的建立杨丽娟;赵丹平【摘要】发动机性能的好坏直接影响着整车运行的平顺性、安全性、稳定性等,要全面评价发动机性能,万有特性曲线则是一个很好的工具.万有特性曲线是以发动机转速为横坐标,以扭矩或平均有效压力为纵坐标,在坐标系内画出等燃油消耗率曲线和等功率曲线[1].绘制万有特性曲线的方法有很多种,MATLAB语言是其中之一.本文利用强大的MATLAB绘图工具,绘制了发动机的万有特性曲线,经分析,该方法是一个有效的精确度较高的方法.【期刊名称】《交通节能与环保》【年(卷),期】2010(000)001【总页数】3页(P32-33,48)【关键词】MATLAB;发动机;万有特性【作者】杨丽娟;赵丹平【作者单位】内蒙古工业大学,内蒙古,赤峰,010051;内蒙古工业大学,内蒙古,赤峰,010051【正文语种】中文0 引言发动机是汽车的动力源，发动机性能的好坏直接影响着整车的动力性与经济性。

汽车的运行工况是个随机的过程，受到很多因素的影响，如道路条件、交通流量、气候条件以及汽车自身技术性能的变化等等。

在所有的运行工况下，发动机都应能够与传动系实现最佳匹配，以使整车动力性、经济性、排放性和噪声污染等方面均处于最佳状态。

然而，对发动机性能的分析与研究是保证整车性能达到最佳的重要前提。

MATLAB语言是一个功能强大的仿真软件，可以完成复杂的数学运算，实现对动态系统的建模仿真等，在工程计算中应用非常广泛。

本文利用MATLAB的强大的绘图功能，建立了发动机的输出转矩模型、油耗模型及万有特性曲线。

通过输出转矩模型，由当前发动机节气门开度及转速既可得到与之相对应的发动机的转出转矩，从而为分析发动机性能奠定基础；发动机油耗模型反应了其有效燃油消耗率与转速和转矩之间的关系，发动机在不同工况下运行时，由此模型既可得到不同工况下发动机的比油耗，从而为分析整车燃油经济性提供数据支持；万有特性曲线是在由发动机转速和转矩构成的坐标系内，绘制出等油耗曲线、等功率曲线、外特性曲线等，通过万有特性曲线既可较全面地了解发动机在不同工况下的性能指标。

MATLAB的发动机万有特性曲线绘制方法程序

%不同转速下的燃油消耗率与扭矩的曲线拟合clear allbe仁[222.8,220.4,232.4,228.5,227.8,232.6,248.5,245.9,272.4,329.7];Ttq仁[399.8,354.1,318.5,278.1,236.2,203.6,185.3,157.2,117.2,80.8];T1=80:320/9:400;%转换矩阵格式Be仁in terp1(Ttq1,be1,T1,'spli ne');% n=1400r/mi n 时燃油消耗率与扭矩的曲线拟合be2=[222.0,221.7,235.4,226.5,230.5,236.8,249.1,276.1,407.9,487.0];Ttq2=[409.1,365.7,328.3,284.1,243.7,203.2,164.3,123.9,83.5,39.7];T2=39:371/9:410;Be2=i nterp1(Ttq2,be2,T2,'spl in e');be3=[226.0,225.3,226.4,233.9,242.1,283.3,253.9,271.4,323.5,468.6];Ttq3=[408.3,368.3,328.3,289.0,244.4,208.8,167.7,132.1,89.5,46.1];T3=46:363/9:409;Be3=i nterp1(Ttq3,be3,T3,'spl in e');be4=[206.5,231.1,231.1,233.0,242.0,244.9,265.0,299.8,398.0,596.8];Ttq4=[425.6,380.3,332.7,290.9,244.4,205.1,160.2,114.5,68.8,30.7];T4=30:396/9:426;Be4=i nterp1(Ttq4,be4,T4,'spl in e');be5=[234.7,259.8,235.5,237.6,242.8,292.3,277.9,308.7,396.2,605.9];Ttq5=[420.7,379.6,334.6,291.6,244.4,202.8,157.5,116.0,74.1,37.8];T5=37:384/9:421;Be5=i nterp1(Ttq5,be5,T5,'spl in e');be6=[174.2,242.2,252.1,287.4,253.6,263.6,290.6,316.8,378.0,518.8];Ttq6=[404.6,360.5,322.7,283.0,243.3,205.5,162.1,124.7,86.8,52.4];T6=52:353/9:405;Be6=i nterp1(Ttq6,be6,T6,'spl in e');be7=[256.9,253.7,253.5,260.0,303.8,280.7,300.6,346.6,435.6,812.9];Ttq7=[378.0,344.7,310.3,264.3,226.1,186.8,154.2,115.3,76334.1];T7=34:344/9:378;Be7=i nterp1(Ttq7,be7,T7,'spl in e');be8=[257.9,295.3,282.4,288.7,301.9,329.7,357.0,475.4,580.3,1080.1];Ttq8=[315.6,275.5,242.5,210.3,178.5,145.6,118.6,72.6,52.8,22.4];T8=22:294/9:316;Be8=i nterp1(Ttq8,be8,T8,'spl in e');B=[Be1';Be2';Be3';Be4';Be5';Be6';Be7';Be8'];N=[1400*o nes(10,1);1600*o nes(10,1);1800*o nes(10,1);2000*o nes(10,1);2200* on es(10,1);2400*o nes(10,1);2600*o nes(10,1);2800*o nes(10,1)];Ttq n=[T1';T2';T3';T4';T5';T6';T7';T8'];G=[o nes(80,1),N,Ttq n,N.A2,N.*Tt qn ,Ttq n. A2];A=G\B;%A 为6*1 矩阵[n,Ttq]=meshgrid(1400:2800,100:600);% 生成n-Ttq 平面上的自变量“格点”矩阵be=A(1)+n. *A (2)++Ttq*A(3)+n.A2*A(4)+n.*Ttq*A(5)+Ttq.A2*A(6);Pe=Ttq.* n/9550;%外特性实验数据拟合Nw=[1403,1597,1797,1986,2102,2199,2303,2400,2507,2598,2700,2802];Ttqw=[474,497,515,526,528.8,522.8,509.5,492.2,471.2,448.4,408.3,357.4]; n0=1400:2800; Ttqw_N=i nterp1(Nw,Ttqw ,n 0,'spli ne');h=repmat(Ttqw_N,501,1);ii=find(Ttq>h);%确定超出边界的“格点”下标be(ii)=NaN;%强制为非数Pe(ii)=NaN;%强制为非数%绘制等燃油消耗率曲线和等功率曲线三维拟合图subplot(1,2,1);mesh( n, Ttq,be);hold on;mesh( n,Ttq,Pe);axis([1000,3000,100,600,0,500]);hold on;xlabel(' n(r/mi n)')ylabel('Ttq(N*m)')zlabel('Pe(KW) be(g/(KW*h))')title('等燃油消耗曲线和等功率曲线的三维拟合图')%绘制边界线(外特性曲线)subplot(1,2,2);plot( nO,Ttqw_N,'Li neWidth',2);axis([1400,2800,100,550]);xlabel(' n(r/mi n)');ylabel('Ttq(N*m)');title('万有特性曲线');hold on;%绘制等油耗率曲线的二维图B=contour(n,Ttq,be,11);%画等位线，并给出标识数据clabel(B);%把“等位值”沿等位线随机标识hold on;%绘制等功率曲线的二维图P=contour(n,Ttq,Pe,11);%画等位线，并给出标识数据clabel(P);%把“等位值”沿等位线随机标识lege nd('等油耗曲线','等功率曲线','外特性曲线')hold off。

汽车发动机万有特性曲线簇计算机绘制算法

因而，三次Ｂ样条曲线矢量方程可表示
为Ｉ
ＣＩ．．（ｕ）＝Ｎ１．．（ｕ）Ｐｉｌ＋Ｎ２． ‘ （ｕ）Ｐｉ＋Ｎ３，４（ｕ）Ｐｉ十１＋Ｎ‘ ， ‘ ）Ｐ｛＋２
（４）
至此，所有的Ｍｊ均已求得。
２．２．２最优逼近曲线
对象，采用本文后面的传统Ｂ样条曲线拟和
方式．拟和的曲线在节点容易出现奇异性，从图形上看就是曲线两次经过同一节点处，即
其中，Ｘ一ｉ￣Ｘ＜Ｘｊ，ｈｊ＝】【ｊ－Ｘｊ一１，ｊ＝１，
２， … … ，Ｎ。
“ 打折，这样的情况必须避免。本文处理方法是以节点作为特征，多边形控翩点构造二阶导数连续的Ｂ样条曲线。若从空间ｎ＋１个顶点ＰｉＯ；０，１， …ｎ）中选取相邻的四个硬点，可构造出一段三次Ｂ样条曲线，其相应的基函数是：
油耗）。
ｘ轴代表转速，Ｙ轴代表扭矩，ｚ代表油
移动式试车线在我国单缸柴油机试车中的真正使用，仅有一年多的时问，其突出的优越性已赢得生产厂家的信赖，并在行业中引起很大反响。随着它在生产实践和使用过程中的不断改进和提高，必将显示出强大的生命力。在多缸汽油机和柴油机的出厂试验中，这种移动式试车线也必将会占有一席之地
・２１・
若ｉ，ｊ的高度值不等，则当ｆ＜Ｏ时，等值线必与棱边有唯一交点；当ｆ＞ｏ，两者无交点Ｉ而当ｆ一０时，等值线过节点ｉ或ｊ，见图２。值得注意的是，ｆ＝Ｏ的情形极易引起后面

MTLAB绘制发动机万有特性MAP详细程序

MTLAB绘制发动机万有特性MAP详细程序采用MATLAB来绘制发动机MAP，通过附表1和2，将这两个表放在同一个根目录下，并编写M文件（见后面附录程序），将M文件跟附表1和2同时放在同一个不带中文的根目录下，点击运行即可。

详细绘图程序见附录。

附表1：转速扭矩燃油消耗率2500 613 221.662501 549.4 222.472500 484.3 221.462500 401.8 227.392499 300.3 242.942500 205.5 272.412498 99.1 375.582300 649.9 212.042300 578.8 212.76附表2：1168 408.821200 477.51337 5001485 514.481500 5731543 618.921401 749.821600 716.251656 749.71760 759.661885 759.952000 5001．将附表1粘贴到Excel中，并将Excel名称更改为12.xlsx文件，并将该12.xlsx 文件存放在E盘根目录下，即E:\12.xlsx2．将附表2粘贴到另外一个Excel中，并将Excel名称更改为hua.xlsx文件，并将hua.xlsx文件存放在E盘根目录下，即E:\hua.xlsx3．在MTLAB中点击File，New，Scrip新建一个M文件4. 将附录后的程序粘贴进去后点击运行，即可生成发动机MAP图。

附录程序：（注意：程序和两个Excel存放在同一个跟目录，因为程序会调用者两张表格）clearclcA=xlsread('E:\12.xlsx');%%%表示12.xlsx文件存放在E盘根目录下x=A(:,1);y=A(:,2);z=A(:,3);xi=(linspace(min(x),max(x),100));size(xi)%%%% linspace是Matlab中的一个指令，用于产生x1,x2之间的N点行矢量。

基于MATLAB的发动机万有特性曲线绘制方法_2_图文(精)

【设计研究】基于 M AT LAB 的发动机万有特性曲线绘制方法周广猛 1, 郝志刚 2, 刘瑞林 1, 陈东 3, 管金发 1, 张春海4(1. 军事交通学院汽车工程系 , 天津 300161;2. 军事交通学院训练部 , 天津300161; 3. 军事交通学院基础部 , 天津 300161;4. 兰州军区军械汽车技工训练大队 , 陕西 710111摘要 :利用 MAT LAB 数学运算能力 , , , 有曲线直观明了 , 把等燃油消耗率曲线、 , 拟合程度较高。

关键词 ; :A文章编号 :1673-6397(2009 02-0034-03U niversal Characteristics Curve Plotting Method based on MAT LABZ H O U G uang -m eng 1,H A O Z hi -gang 2, L I U Rui -lin 1,CHE N D ong 3,G U A N Jin -fa 1,Z H A NG Chun -hai 4(1. Autom obile Engineering Department , Academy of Military T ransportation , T ianjin 300161,China ;2. T raining Department ,Academy of Military T ransportation , T ianjin 300161,China ;3. G eneral C ourse Department , Academy of Military T ransportation , T ianjin 300161,China ;4. Ordnance Mechanic T raining Brigade , Lan Zhou Theater , X i ’ an 710111,China Abstract :Taking advantage of MAT LAB mathematic operation , data from engine characteristic test was processed , the method is sim ple and credible , The universal characteristics curve plotted is intuitionistic and perspicuous ,and was in g ood fit with data g ot in test.K ey Words :MATLAB ;Universal Characteristics Curve ;Plot作者简介 :周广猛 (1984- , 男 , 山东邹城人 , 在读硕士研究生 , 主要研究方向为动力机械特殊环境适应性。

MATLAB曲线拟合在发动机特性拟合中的应用

MATLAB曲线拟合工具箱在发动机特性拟合中的应用摘要：详细介绍了MA TLAB的曲线拟合工具箱（Curve Fitting Toolbox）的GUI交互式图形界面和函数调用两种方式的使用方法，并以某柴油发动机稳态工况下的试验数据绘制万有特性图为例，介绍了利用曲线拟合工具箱进行曲面拟合的步骤和方法，并给出了关键代码。

对拟合结果的数理统计校验表明，拟合结果完全符合要求。

结果表明，MA TLAB的曲线拟合工具箱为发动机特性的曲面拟合提供了极大的便利。

主题词：发动机特性，曲面拟合，多元回归，曲线拟合工具箱中图分类号: TK421 文献标识码：AApplication of MATLAB Curve Fitting Toolbox on Engine Characteristic Maps FittingAbstract：The usage of MATLAB Curve Fitting Toolbox GUI mode and function mode is detailed introduced. The procedures and method of surface fitting employ Curve Fitting Toolbox are demonstrated with the universal characteristic map plotting of a diesel engine steady state test data. Key codes are provided also. The statistics results show that the fitting accuracy is very high. The results prove that MATLAB Curve Fitting Toolbox can provide much convenience to the surface fitting of engine characteristics.Key Words：Engine Characteristic Map, Surface Fitting, Multiple Regression, Curve Fitting Toolbox引言根据发动机稳态工况下的性能试验数据绘制转速特性、负荷特性和万有特性曲线的过程较为繁琐，工作量较大，目前一般应用计算机编制程序对大量的试验数据进行处理，从而得到万有特性曲线[1] ;尤其是MA TLAB 软件，因其强大的数值计算和绘图功能，越来越广泛的被应用于发动机万有特性的研究[2,3,4]。

基于MATLAB的发动机万有特性曲面拟合

２０１４年第３期第２０卷（总第１４８期）
基于ＭＡＴＬＡＢ的发动机万有特性曲面拟合
黄风清
（同济大学汽车学院，上海２０００９２）摘要基于软件建模的理论研究方法在整车开发中占据了越来越重要的地位。发动机作为整车中最为重要的一个零部件，发动机的万有特性如何通过数学模型更准确地表达出来，对整车模
Ａｂｓｔｒａｃｔ：Ａｒｅｓｅａｒｃｈｍｅｔｈｏｄｏｆｂｕｉｌｄｉｎｇｍａｔｈｅｍａｔｉｃａｌｍｏｄｅｌｓｂａｓｅｄｏｎｐｒｏｆｅｓｓｉｏｎｌａｓｏｆｔｗａｒｅｉｓ
柴油机设计与制造
Ｄｅｓｉｇｎ＆ＭａｎｕｆａｃｔｕｒｅｏｆＤｉｅｓｅｌＥｎｇｉｎｅ
ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．１６７１－０６１４．２０１４．０３．００３
ｏｆｇｒｅａｔｓｉｇｎｉｉｆｃａｎｃｅｔｏｔｈｅｐｒｅｃｉｓｉｏｎｏｆｖｅｈｉｃｌｅｃｌａｃｕｌａｔｉｏｎ．Ａｍｏｒｅａｃｃｕｒａｔｅｉｆｔｔｉｎｇｅｑｕａｔｉｏｎｆｏｒｅｎｉｇｎｅ
ｐｅｆｏｒｒｍａｎｃｅｍａｐｉｓｄｅｖｅｌｏｐｅｄｂｙｂｕｉｌｄｉｎｇｓｅｖｅｒｌａｄｉｆｆｅｒｅｎｔｍｏｄｅｌｓａｎｄａｄｏｐｔｉｎｇｔｈｅＭＡＴＬＡＢｉｆｔｔｉｎｇｃｌｃａｕｌａｔｉｏｎｆｕｎｃｔｉｏｎ．

MATLAB万有特性曲线

%不同转速下的燃油消耗率与扭矩的曲线拟合clear allbe1=[222.8,220.4,232.4,228.5,227.8,232.6,248.5,245.9,272.4,329.7];Ttq1=[399.8,354.1,318.5,278.1,236.2,203.6,185.3,157.2,117.2,80.8];T1=80:320/9:400; %转换矩阵格式Be1=interp1(Ttq1,be1,T1,'spline'); %n=1400r/min时燃油消耗率与扭矩的曲线拟合be2=[222.0,221.7,235.4,226.5,230.5,236.8,249.1,276.1,407.9,487.0];Ttq2=[409.1,365.7,328.3,284.1,243.7,203.2,164.3,123.9,83.5,39.7];T2=39:371/9:410;Be2=interp1(Ttq2,be2,T2,'spline');be3=[226.0,225.3,226.4,233.9,242.1,283.3,253.9,271.4,323.5,468.6];Ttq3=[408.3,368.3,328.3,289.0,244.4,208.8,167.7,132.1,89.5,46.1];T3=46:363/9:409;Be3=interp1(Ttq3,be3,T3,'spline');be4=[206.5,231.1,231.1,233.0,242.0,244.9,265.0,299.8,398.0,596.8];Ttq4=[425.6,380.3,332.7,290.9,244.4,205.1,160.2,114.5,68.8,30.7];T4=30:396/9:426;Be4=interp1(Ttq4,be4,T4,'spline');be5=[234.7,259.8,235.5,237.6,242.8,292.3,277.9,308.7,396.2,605.9];Ttq5=[420.7,379.6,334.6,291.6,244.4,202.8,157.5,116.0,74.1,37.8];T5=37:384/9:421;Be5=interp1(Ttq5,be5,T5,'spline');be6=[174.2,242.2,252.1,287.4,253.6,263.6,290.6,316.8,378.0,518.8];Ttq6=[404.6,360.5,322.7,283.0,243.3,205.5,162.1,124.7,86.8,52.4];T6=52:353/9:405;Be6=interp1(Ttq6,be6,T6,'spline');be7=[256.9,253.7,253.5,260.0,303.8,280.7,300.6,346.6,435.6,812.9];Ttq7=[378.0,344.7,310.3,264.3,226.1,186.8,154.2,115.3,76.3,34.1];T7=34:344/9:378;Be7=interp1(Ttq7,be7,T7,'spline');be8=[257.9,295.3,282.4,288.7,301.9,329.7,357.0,475.4,580.3,1080.1];Ttq8=[315.6,275.5,242.5,210.3,178.5,145.6,118.6,72.6,52.8,22.4];T8=22:294/9:316;Be8=interp1(Ttq8,be8,T8,'spline');B=[Be1';Be2';Be3';Be4';Be5';Be6';Be7';Be8'];N=[1400*ones(10,1);1600*ones(10,1);1800*ones(10,1);2000*ones(10,1);2200*ones(10,1);2400*ones (10,1);2600*ones(10,1);2800*ones(10,1)];Ttqn=[T1';T2';T3';T4';T5';T6';T7';T8'];G=[ones(80,1),N,Ttqn,N.^2,N.*Ttqn,Ttqn.^2];A=G\B; %A为6*1矩阵[n,Ttq]=meshgrid(1400:2800,100:600); %生成n-Ttq平面上的自变量“格点”矩阵be=A(1)+n.*A(2)++Ttq*A(3)+n.^2*A(4)+n.*Ttq*A(5)+Ttq.^2*A(6); %501×1401Pe=Ttq.*n/9550; %501×1401%外特性实验数据拟合Nw=[1403,1597,1797,1986,2102,2199,2303,2400,2507,2598,2700,2802]; Ttqw=[474,497,515,526,528.8,522.8,509.5,492.2,471.2,448.4,408.3,357.4]; n0=1400:2800;Ttqw_N=interp1(Nw,Ttqw,n0,'spline');h=repmat(Ttqw_N,501,1); % 501×1401矩阵ii=find(Ttq>h); %确定超出边界的“格点”下标 %155109×1be(ii)=NaN; %强制为非数Pe(ii)=NaN; %强制为非数%绘制等燃油消耗率曲线和等功率曲线三维拟合图subplot(1,2,1);mesh(n,Ttq,be);hold on;mesh(n,Ttq,Pe);axis([1000,3000,100,600,0,500]);hold on;xlabel('n(r/min)')ylabel('Ttq(N*m)')zlabel('Pe(KW) be(g/(KW*h))')title('等燃油消耗曲线和等功率曲线的三维拟合图')%绘制边界线(外特性曲线)subplot(1,2,2);plot(n0,Ttqw_N,'LineWidth',2);axis([1400,2800,100,550]);xlabel('n(r/min)');ylabel('Ttq(N*m)');title('万有特性曲线');hold on;%绘制等油耗率曲线的二维图B=contour(n,Ttq,be,11); %画等位线，并给出标识数据clabel(B); %把“等位值”沿等位线随机标识hold on;%绘制等功率曲线的二维图P=contour(n,Ttq,Pe,11); %画等位线，并给出标识数据clabel(P); %把“等位值”沿等位线随机标识legend('等油耗曲线','等功率曲线','外特性曲线')hold off%利用mesh作原始曲面N=3500:500:8000;T=1.3:-1.2/6:0.1;B=xlsread('F:\Matlab\RanJia.xls','revise5');%mesh(N,T,B)%colormap;%colorbar;%xlabel('转速n/ r*min^-^1'), ylabel('p平均有效压力/ Mpa'), zlabel('燃油消耗率z/ kg/(kw*h)')%hidden off 透明网孔%colormap; 表面小块着色%colorbar;%surf(,,,'FaceColor','red','EdgeColor','none');%camlight left; 增加光源%lighting phone 照明方式%view(-15,65) 改变视角（方位角，仰角）%利用interrp2矩阵插值函数作优化曲面N1=3500:50:8000;T1=1.3:-0.005:0.1;[N2,T2]=meshgrid(N1,T1);B1=interp2(N,T,B,N2,T2,'cubic');P1=T2.*N2/9.55;%figure;%surf(N1,T1,B1)%colormap;%colorbar;%xlabel('转速n/ r*min^-^1'), ylabel('p平均有效压力/ Mpa'), zlabel('燃油消耗率z/ kg/(kw*h)')%外特性实验数据拟合T3=[0.90,0.95,0.99,1.06,1.12,1.18,1.24,1.30,1.22,1.11];T4=interp1(N,T3,N1,'spline');W=repmat(T4,241,1); %平铺成381×91矩阵jj=find(T2>W); %确定超出边界的“格点”下标B1(jj)=NaN;%画外特性曲线plot(N1,T4,'LineWidth',2);axis([3500,8000,0.1,1.3]);hold on%画等油耗线v=[450,460,470,480,490,500,510,525,540,560];[A1,h]=contour(N2,T2,B1,v,'b:'); %N2,T2限制在X,Y轴上的范围，10为等高线条数clabel(A1,h,'manual');hold on%画等功率线[A2,h]=contour(N2,T2,P1,500:100:900,'k:'); %N2,T2限制在X,Y轴上的范围，10为等高线条数clabel(A2,h,'manual');plot(6000,1.114,'.','color','r');xlabel('转速—r/min');ylabel('扭矩—N ·m');legend('外特性曲线','等油耗线-g/(kW ·h)','等功率线-W') hold off% title('等油耗线');%xlabel('转速n/ r*min^-^1'), ylabel('p 平均有效压力/ Mpa')转速（r/min ）扭矩（N ·m ）。

MATLAB曲线拟合工具箱在发动机特性拟合中的应用

ａ０，口２， …，Ｏｋ一１）为模型中的待定系数；ｂｅ＝［１，ｎｅ，Ｍｅ，２，，， … ，１－１，］ × ● ：十 ● ：式中，（
ａｋ
一
１Ｉｌｅｋ一１
（，ｅ２， …，一。）为随机误差； Ⅳ 为试验数据点数；，为多项式的最高幂次；七为多项式的项数，（ｆ＋１）（１＋２）／２．
黄兵锋，等：ＭＡＴＬＡＢ曲线拟合工具箱在发动机特性拟合中的应用
采用拟合算法时应选择合适的最高幂次．虽然采用较高的多项式最高幂次将使得拟合曲线与试验点数值更逼近，但次数过高，容易出现龙格（Ｒｕｎｇｅ）现象，在试验点以外的数据上偏差可能很大；同时考虑到试验点数据受偶然误差等随机因素的影响，曲线拟合最高幂次一般为２～５，最大不超过７［８】．曲线拟合时，ＭＡＴＬＡＢ的曲线拟合工具箱允许最高幂次为９，曲面拟合时最高幂次为双５次（即ｎ，＾的最高幂次均为５）．
２．吉林大学汽车工程学院汽车仿真与控制国家重点实验室，吉林长春１３００２２；３．武汉科技大学信息科学与工程学院，湖北武汉４３００８１）摘要：文章阐述ＭＡＴＬＡＢ的曲线拟合工具箱（ＣｕｒｖｅＦｉｔｔｉｎｇＴｏｏｌｂｏｘ）的交互式图形界面和程序调用两种方式的使用方法，并以某柴油发动机稳态工况下的试验数据绘制万有特性图为例，介绍利用曲线拟合工具箱进行曲面拟合的步骤和方法，并给出关键代码，探讨利用该工具箱改进拟合效果和绘图质量．对拟合结果的数理统计校验表明拟合结果完全符合要求，ＭＡＴＬＡＢ的曲线拟合

用origin绘制万有特性的过程1

利用Origin绘制发动机万有特性曲线以绘制CYQD32Ti型柴油机万有特性曲线为例，介绍在Origin软件环境下绘制发动机万有特性曲线的过程。

首先，打开Origin软件，选择 File-Save Project As，命名为CYQD32Ti万。

一、输入原始数据。

先从EXCEL中处理好负荷特性试验数据、外特性数据和等功率数据，把外特性数据的转速、平均有效压力2列数据粘贴到表1中，把负荷特性的转速、平均有效压力、燃油消耗率3列数据粘贴到Origin软件的表2中，把等功率数据中的转速和各功率下的平均有效压力共11列数据粘贴到表3中。

结果见图1。

二、将负荷特性数据表的燃油消耗率Y列属性转化为Z，为生成为矩阵做准备。

点选Data2负荷特性数据表的燃油消耗率列，在右键菜单中选择Set as Z。

见图2。

三、将负荷特性数据转化为矩阵。

从菜单中选择Edit-Convert to Matrix-Random XYZ，在出现的Random XYZ Gridding对话框中选择OK按钮，生成Matrix1。

过程见图3。

四、利用矩阵生成等燃油消耗率曲线。

选中矩阵Matrix1，从菜单中选择Plot-Contour Plot-Coutour-B/W Lines+Labels，过程见图4。

等油耗曲线图Graph1见图5。

五、在等油耗线图中加入等功率线数据和外特性数据。

在Graph1窗口左上角的1处单击右键，在快捷菜单中选择Layer Contents。

在出现的Layer1对话框中将左框中与Data3 有关的数据添加到右侧，再把Data1的数据添加到右侧，和单击OK按钮。

过程见图6。

Graph1中已经加入等功率线和外特性曲线，生成万有特性曲线雏形，见图7。

六、设置图形外观，使其满足万有特性曲线报告要求。

1.在Graph1图形区域双击打开Plot Detail 对话框，设置等油耗线的间距。

见图8。

设置等油耗线标签的有效位数，见图9。

基于MATLAB的发动机万有特性曲面拟合

基于MATLAB的发动机万有特性曲面拟合黄风清【摘要】基于软件建模的理论研究方法在整车开发中占据了越来越重要的地位.发动机作为整车中最为重要的一个零部件,发动机的万有特性如何通过数学模型更准确地表达出来,对整车模型的计算精度有非常重要的意义.通过建立不同的数学模型,依靠MATLAB的曲面拟合计算功能,找到更为精确的发动机万有特性曲面拟合方程.【期刊名称】《柴油机设计与制造》【年(卷),期】2014(020)003【总页数】4页(P14-16,56)【关键词】发动机;数学模型;MATLAB【作者】黄风清【作者单位】同济大学汽车学院,上海200092【正文语种】中文基于软件建模的理论研究方法在整车设计开发中发挥了越来越重要的作用。

通过建立整车的模型，形成虚拟样机，可以更好地辅助设计师选择更优的设计参数、较准确地评估整车性能以及预测整车的各项参数。

整车油耗的模拟计算往往是整个模型中的重点与难点，因为影响整车油耗因素有：发动机本体的性能、传动效率、外部阻力、使用工况、人为操作习惯等。

针对如何提升模型对整车油耗计算的准确性，本文主要研究如何提高对发动机万有特性拟合的精度，以达到提高模拟计算精度的目的。

本文的研究对象为上柴动力SC9DK220G3型发动机。

该款发动机的供油系统为高压共轨，排放限值达到非道路用第Ⅲ阶段的要求，是上柴公司为非道路用第Ⅲ阶段排放打造的升级动力，具有油耗低、加速性好等优点。

其主要参数如表1所示。

发动机的万有特性曲线可以较全面地表示发动机的性能。

通常以发动机的转速为横坐标，发动机的扭矩为纵坐标，在图上可以表示出等油耗率曲线、等小时耗油量曲线、等功率曲线、等过量空气系数曲线、冒烟极限等。

3.1 基于多项式的数学模型根据发动机负荷特性定义，固定发动机转速，对表达式进行降维处理[1]，对其进行一元曲线拟合，其燃油消耗率be的表达式为：式（1）表示为发动机在转速ni下的负荷特性，若固定j，则式（1）中的Aij均为转速n的函数，其表达式可写为：将式（2）代入式（1），可以得到通用的计算公式其中i=(1,2,3…K)，j=(1,2,3…S)根据式（3）及图1燃油消耗率be与转速n和功率P的曲线关系，结合计算复杂程度，分别选取K=S=3进行拟合计算。

汽车万有特性曲线绘制

Ttq=xlsread('d:\Program Files\MATLAB\R2010b\bin\wanyoutexing1\wanyoutexing.xls','sheet1','b4:b83');
be=xlsread('d:\Program Files\MATLAB\R2010b\bin\wanyoutexing1\wanyoutexing.xls','sheet1','c4:c83');
ylabel('Ttq/N.m')
%%%%%%%%%%%%%%%%拟合外特性曲线% 1800 2000 2200 2400 2600 2800];
Ttq_max0=[399.8 409.1 408.3 425.6 420.7 404.6 378 315.6];
be1=A(1)+A(2)*n1+A(3)*Ttq1+A(4)*n1.^2+A(5)*n1.*Ttq1+A(6)*Ttq1.^2;
%%%%%%%%%%%%%绘制be-pme-n%%%%%%%%%%%%%
mesh(n1,Ttq1,be1)
title('be-Ttq-n三维曲面')
xlabel('n / (r/min)')
Z=[ones(size(n)),n,Ttq,n.^2,n.*Ttq,Ttq.^2];
%%%%%%%%%%%%%%%%求出回归参数%%%%%%%%%%%%%%%%%
A=Z\be;
n1=1400:2:2800;
Ttq1=20:0.1:426;
[Ttq1,n1]=meshgrid(Ttq1,n1);

发动机及驱动力曲线绘制Matlab绘图详解

资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
国内最大的维修资料库精通维修下载站
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/

资料来源：国内最大的维修资料库精通维修下载站：/
资料来源：国内最大的维修资料库精通维修下载站：/

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

打开origin，将excel中整理好的数据直接复制粘贴到Book1中即可，可以在左下方的信息栏对Book1进行重命名。

origin中表格操纵与excel中类似。

可以编辑数据的名称、单位、备注等信息，也可空着以后再绘制好的图表上修改。

同时选中燃油消耗率和功率数据两列，点右键选择Set as Z，也可以在Column菜单下点选Set as Z，如下图所示。

二、绘制万有特性曲线选中表格中所有数据列表，在绘图命令菜单Plot下绘制等高线命令Contour的颜色填充Color Fill选项，将出现图表窗口，如下图所示。

三、万有特性曲线图调整上一步完成的万有特性曲线只是一个雏形，与常见的还很不一样，需要进行调整细化。

在已绘制好的万有特性曲线图中，由于点选的是颜色填充的绘制方法，数据源有两组Z轴分量，相当于在一张图上绘制了两层，而等燃油消耗率曲线在等功率曲线的下方，被覆盖住无法看到。

因此，需要取消等功率曲线图层的颜色填充效果。

如图所示，在图片窗口左上角的1上右击，选择图层属性Layer Properties命令，在弹出的对话框中打开图层Layer1的下一级，选中转速、扭矩、功率曲线，并去掉Color Fill Control下Enabled之前的勾选，即可取消等功率曲线图层的颜色填充效果。

等燃油消耗率曲线还需要进行进一步的调整，才能变成最常见的样子。

需要调整曲线的层次间隔，让等值曲线分布的疏密合理，还要加上等值线的数据标注，以便于观看查阅。

同样是在图层属性对话框里，选中转速、扭矩、燃油消耗率曲线，在列标题或单元格编辑表格中，单击列标题或单元格可以完成曲线的层次、填充颜色、层次线型及添加数据标注的修改。

如图所示，单击Level栏表头，弹出Set Level设置层次对话框，先点击find min/max命令，找到燃油消耗率的最小值与最大值，并自动设置为等值曲线变化范围，选择线性变化，并选中增量increment选项，设置增量值为1，既让等燃油消耗率曲线从184.5g/kw.h开始，每隔1g/kw.h就绘制一条，直至367g/kw.h结束，一共183条。

单击Fill栏表头，在弹出的填充编辑对话框中，可以修改等值线之间的颜色变化，可以由一种颜色直接过渡到另一种，也可以中间插入其他颜色过渡，或者选择一种确定的颜色变化趋势，还可以为曲线之间添加纹理过渡趋势等。

单击标注Labels栏表头，弹出标注编辑对话框，选择show all命令显示所有等燃油消耗率曲线的数值标注。

点击apply应用后，效果如图所示。

调整等值线的疏密程度及平滑度，调整标注尺寸大小以使特性图变得更美观。

还是在曲线细节与图层属性对话框中，曲线密的地方选中曲线，单击Level下方的Delete命令即可删除，同样道理，在曲线疏的地方点击Insert即可插入新曲线。