Æû³µÀíÂÛmatlab×÷Òµ ÏÂÔر¾ÎÄ

a2=(Ft2-F2)./(Q(2)*m); t2=trapz(ua2,1./a2)/3.6;

ua3=uamax(2):uamax(3);

n3=ua3*Ig0(k)*Ig5(3)/(0.377*0.367);

Ttq3=-19.313+295.27*(n3/1000)-165.44*(n3/1000).^2+40.874*(n3/1000).^3-3.8445*(n3/1000).^4;

Ft3=Ttq3*Ig0(k)*Ig5(3)*0.85/0.367;

F3=0.013*3880*9.8+2.77*ua3.^2/21.15; a3=(Ft3-F3)./(Q(3)*m); t3=trapz(ua3,1./a3)/3.6;

ua4=uamax(3):70;

n4=ua4*Ig0(k)*Ig5(4)/(0.377*0.367);

Ttq4=-19.313+295.27*(n4/1000)-165.44*(n4/1000).^2+40.874*(n4/1000).^3-3.8445*(n4/1000).^4;

Ft4=Ttq4*Ig0(k)*Ig5(4)*0.85/0.367;

F4=0.013*3880*9.8+2.77*ua4.^2/21.15; a4=(Ft4-F4)./(Q(4)*m); t4=trapz(ua4,1./a4)/3.6; t(k)=t2+t3+t4; end

f=polyfit(Qz,t,2);

QQ=Qz(1):0.001:Qz(5); b=polyval(f,QQ) plot(QQ,b,Qz,t,'*')

gtext('5.17') gtext('5.43') gtext('5.83') gtext('6.17') gtext('6.33')

title('ȼÓ;­¼ÃÐÔ-¶¯Á¦ÐÔÇúÏß') xlabel('ȼÓ;­¼ÃÐÔ/[L(100km)-1]') ylabel('t/s')

ËÄ¡¢Ò»ÖÐÐÍ»õ³µ×°ÓÐÇ°ºóÖƶ¯Æ÷·Ö¿ªµÄË«¹Ü·Öƶ¯Ïµ£¬ÆäÓйزÎÊýÈçÏ£º

ÖÊÐÄÖÁÖÊÁ¿ÔØºÉ £¨kg£© ¿ÕÔØ ÂúÔØ 1£© 2£©

4080 9290 hg/m Àëa/m 0.845 1.170 3.950 3.950 2.100 2.950 Êý¦Â 0.38 0.38 ÖÊÐĸßÖá¾àL/m Ç°Öá¾à·ÖÅäϵÖƶ¯Á¦¼ÆËã²¢»æÖÆÀûÓø½×ÅϵÊýÇúÏߺÍÖƶ¯Ð§ÂÊÇúÏß

ÇóÐÐÊ»³µËÙUa£½30km/h£¬ÔÚ?£½0.80·ÃæÉϳµÂÖ²»±§ËÀµÄ

Öƶ¯¾àÀë¡£¼ÆËãʱȡÖƶ¯Ïµ·´Ó¦Ê±¼ä?2'£½0.02s£¬Öƶ¯¼õËÙ¶ÈÉÏÉýʱ¼ä?2''£½0.02s¡£ 3£©

ÇóÖƶ¯ÏµÇ°²¿¹Ü·Ëð»µÊ±Æû³µµÄÖƶ¯¾àÀës£¬Öƶ¯Ïµºó²¿

¹Ü·Ëð»µÊ±Æû³µµÄÖƶ¯¾àÀës'¡£ Matlab³ÌÐò£º

m1=4080;hg1=0.845;a1=2.100; m2=9290;hg2=1.17;a2=2.95; beta=0.38;L=3.95; z=0:0.05:1

gf1=beta.*z*L./(L-a1+z*hg1); gf2=beta.*z*L./(L-a2+z*hg2); gr1=(1-beta).*z*L./(a1-z*hg1); gr2=(1-beta).*z*L./(a2-z*hg2); g=z;

for i=1:21

if (z(i)<0.3&z(i)>0.15); g3(i)=z(i)+0.08;

end

if(z(i)>=0.3);

g3(i)=0.38+(z(i)-0.3)/0.74; end end

z1=0.15:0.01:0.3; g4=z1-0.08;

plot(z,gf1,'-.',z,gf2,z,gr1,'-.',z,gr2,z,g,z,g3,'xk',z1,g4,'x') axis([0 1 0 1.2])

title('ÀûÓø½×ÅϵÊýÓëÖƶ¯Ç¿¶ÈµÄ¹ØϵÇúÏß') xlabel('Öƶ¯Ç¿¶Èz/g') ylabel('ÀûÓø½×ÅϵÊýg') gtext('¿Õ³µÇ°Öá') gtext('¿Õ³µºóÖá') gtext('ÂúÔØÇ°Öá') gtext('ÂúÔغóÖá') gtext('ECE·¨¹æ')

C=0:0.05:1;

Er1=(a1/L)./((1-beta)+C*hg1/L)*100; Ef=(L-a2)/L./(beta-C*hg2/L)*100; Er=(a2/L)./((1-beta)+C*hg2/L)*100; plot(C,Er,C,Ef,C,Er1) axis([0 1 0 100]) title('Ç°ºó¸½×ÅЧÂÊÇúÏß') xlabel('¸½×ÅϵÊýC') ylabel('Öƶ¯Ð§ÂÊ£¨%£©') gtext('ÂúÔØ') gtext('Ef') gtext('Er') gtext('¿ÕÔØ') gtext('Er')