1、%不同轉(zhuǎn)速下的燃油消耗率與扭矩的曲線(xiàn)擬合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; %轉(zhuǎn)換矩陣格式Be1=interp1(Ttq1,be1,T1,spline); %n=1400r/min時(shí)燃油消耗率與扭矩的曲線(xiàn)擬合be2=222.0,221.7,235.4,226.5,230.5,236.8,249.1,276.1,407.9,48

2、7.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,

3、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:38

4、4/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.

5、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

6、=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=GB; %A為6*1矩陣n,Ttq=meshgrid(1400:2800,100:600); %生成n-Ttq平面上的自變量“格點(diǎn)”矩陣be=A(1)+n.*A(2)+Ttq*A(3)+n.2*A(4)

7、+n.*Ttq*A(5)+Ttq.2*A(6); %5011401Pe=Ttq.*n/9550; %5011401%外特性實(shí)驗(yàn)數(shù)據(jù)擬合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); % 5011401矩陣ii=find(Ttqh); %確定超

8、出邊界的“格點(diǎn)”下標(biāo) %1be(ii)=NaN; %強(qiáng)制為非數(shù)Pe(ii)=NaN; %強(qiáng)制為非數(shù)%繪制等燃油消耗率曲線(xiàn)和等功率曲線(xiàn)三維擬合圖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(等燃油消耗曲線(xiàn)和等功率曲線(xiàn)的三維擬合圖)%繪制邊界線(xiàn)(外特性曲線(xiàn))subplot(1,2,2);plot(n0,Ttqw_N,LineWidth,2);a

9、xis(1400,2800,100,550);xlabel(n(r/min);ylabel(Ttq(N*m);title(萬(wàn)有特性曲線(xiàn));hold on;%繪制等油耗率曲線(xiàn)的二維圖B=contour(n,Ttq,be,11); %畫(huà)等位線(xiàn)，并給出標(biāo)識(shí)數(shù)據(jù)clabel(B); %把“等位值”沿等位線(xiàn)隨機(jī)標(biāo)識(shí)hold on;%繪制等功率曲線(xiàn)的二維圖P=contour(n,Ttq,Pe,11); %畫(huà)等位線(xiàn)，并給出標(biāo)識(shí)數(shù)據(jù)clabel(P); %把“等位值”沿等位線(xiàn)隨機(jī)標(biāo)識(shí)legend(等油耗曲線(xiàn),等功率曲線(xiàn),外特性曲線(xiàn))hold off%利用mesh作原始曲面N=3500:500:8000;T=1

10、.3:-1.2/6:0.1;B=xlsread(F:MatlabRanJia.xls,revise5);%mesh(N,T,B)%colormap;%colorbar;%xlabel(轉(zhuǎn)速n/ r*min-1), ylabel(p平均有效壓力/ Mpa), zlabel(燃油消耗率z/ kg/(kw*h)%hidden off 透明網(wǎng)孔%colormap; 表面小塊著色%colorbar;%surf(,FaceColor,red,EdgeColor,none);%camlight left; 增加光源%lighting phone 照明方式%view(-15,65) 改變視角（方位角，仰角）%

11、利用interrp2矩陣插值函數(shù)作優(yōu)化曲面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(轉(zhuǎn)速n/ r*min-1), ylabel(p平均有效壓力/ Mpa), zlabel(燃油消耗率z/ kg/(kw*h)%外特性實(shí)驗(yàn)數(shù)據(jù)擬合T3=0.90,0.95,0.99,1.06,1.12,1.18,1.24,1.30,1.22,1.11;T4=

12、interp1(N,T3,N1,spline);W=repmat(T4,241,1); %平鋪成38191矩陣jj=find(T2W); %確定超出邊界的“格點(diǎn)”下標(biāo)B1(jj)=NaN;%畫(huà)外特性曲線(xiàn)plot(N1,T4,LineWidth,2);axis(3500,8000,0.1,1.3);hold on%畫(huà)等油耗線(xiàn)v=450,460,470,480,490,500,510,525,540,560;A1,h=contour(N2,T2,B1,v,b:); %N2,T2限制在X,Y軸上的范圍，10為等高線(xiàn)條數(shù)clabel(A1,h,manual);hold on%畫(huà)等功率線(xiàn)A2,h=contour(N2,T2,P1,500:100:900,k:); %N2,T2限制在