汽车理论编程题.docVIP

1.3 n=600:1:4000; r=0.367; i0=5.83; eff=0.85; f=0.013; m=3880; g=9.8; G=m*g; CdA=2.77; a=1.947; hg=0.9; L=3.2; Iw1=1.798; Iw2=3.598; Iw=Iw1+Iw2; If=0.218; Ttq=-19.313+295.27*n/1000-165.44*power(n,2)/1000000+40.874*power(n,3)/1000000000-3.8445*power(n,4)/1000000000000; for ig=[5.56 2.769 1.644 1.00 0.793]; Ua=0.377*r*n/ig/i0; Ft=Ttq*ig*i0*eff/r; Ff=m*g*f; Fw=CdA*power(Ua,2)/21.12; Fs=Ff+Fw; plot(Ua,Ft,Ua,Fs); hold on; end text(15,13000,F_t_1); text(20,7000,F_t_2); text(40,4000,F_t_3); text(55,2500,F_t_4); text(70,1900,F_t_5); text(105,2100,F_f+F_w); xlabel(u_a km/h