p=[-3.8445 40.874 -165.44 295.27 -19.313];
Tq=polyval(p,n2);
ua=0.377.*r.*n1./(ig1.*i0);
Ft1=Tq(1,:).*i0.*ig1(1,:).*eta./r;
Ft2=Tq(2,:).*i0.*ig1(2,:).*eta./r;
Ft3=Tq(3,:).*i0.*ig1(3,:).*eta./r;
Ft4=Tq(4,:).*i0.*ig1(4,:).*eta./r;
Ft5=Tq(5,:).*i0.*ig1(5,:).*eta./r;
ua1=[0:max(ua(5,:))];
Fw=cda.*ua1.^2./21.15;
Ff=9.8*3880*f;
Fz=Fw+Ff;
grid on
box on
hold on
plot(ua(1,:),Ft1,'k','linewidth',2);
plot(ua(2,:),Ft2,'--b','linewidth',2);
plot(ua(3,:),Ft3,'-.g','linewidth',2);
plot(ua(4,:),Ft4,':m','linewidth',2);
plot(ua(5,:),Ft5,'linewidth',2);
plot(ua1,Fz,'r','linewidth',2);
xlabel('ua/(km/h)'),ylabel('F/N'),title(‘汽车驱动力行驶阻力平衡图’);
legend('Ft1','Ft2','Ft3','Ft4','Ft5','Fz');
%fff=@(n)0.13.*(0.377.*0.367.*n/(0.793*5.83)).^2-
(0.377.*0.367.*n/(0.793*5.83))+494.312;
% fff=@(n)n.^2-50000;
%fff=@(n)2.77./21.15.*(0.377.*0.367./(0.793.*5.83).*n).^2-
0.377.*0.367./(0.793*5.83).*n+3880*9.8*0.013;
fff=@(n)0.000117.*n.^2+494.312-0.793.*5.83.*0.85./0.367.*(-
19.313+295.27.*(n/1000)-165.44.*(n/1000).^2+40.874.*(n/1000).^3-
3.8445.*(n/1000).^4);
[x,favl]=fzero(fff,4000);
uamax=0.377.*r.*x/(0.793.*i0)
计算结果;
5 / 24
评论0
最新资源