求助，一个汽车行驶加速度倒数曲线拟合程序

大树YJS

求助，一个汽车行驶加速度倒数曲线拟合程序，请大侠们帮帮忙，急用啊，急用，谢谢了，谢谢了
下面是自己写得一个，但运行出事错的曲线，也请大侠们帮忙改改
n_min=600;n_max=4000;
m=2305;
eT=0.88;
r=0.322;
f=0.013;
i0=4.55;
ig=[3.5 2.16 1.47 1 0.79];
g=9.8;
CDA=2.4;
n=(n_min:10:n_max);
N=n./1000;
Tq=-19.313+295.27.*N-165.44.*N.^2+40.874.*N.^3-3.8445.*N.^4;
Ft1=Tq.*ig(1).*i0.*eT./r;
Ft2=Tq.*ig(2).*i0.*eT./r;
Ft3=Tq.*ig(3).*i0.*eT./r;
Ft4=Tq.*ig(4).*i0.*eT./r;
Ft5=Tq.*ig(5).*i0.*eT./r;
s=2*pi*r*0.001;
Ua1=n./ig(1)./i0.*60.*s;
Ua2=n./ig(2)./i0.*60.*s;
Ua3=n./ig(3)./i0.*60.*s;
Ua4=n./ig(4)./i0.*60.*s;
Ua5=n./ig(5)./i0.*60.*s;
G=m*g;
ua=(0:1:340);
Fw=CDA.*ua.^2./21.15;
D1=(Ft1-Fw)./G;
D2=(Ft2-Fw)./G;
D3=(Ft3-Fw)./G;
D4=(Ft4-Fw)./G;
D5=(Ft5-Fw)./G;
a11=(g.*(D1-f))./1.22;
a22=(g.*(D2-f))./1.09;
a33=(g.*(D3-f))./1.05;
a44=(g.*(D4-f))./1.03;
a55=(g.*(D5-f))./1.02;
AA1=(1./a11);
AA2=(1./a22);
AA3=(1./a33);
AA4=(1./a44);
AA5=(1./a55);
plot(Ua1,AA1,Ua2,AA2,'r',Ua3,AA3,'b',Ua4,AA4,'g',Ua5,AA5,'r')

yang10804050232

>> n_min=600;n_max=4000;
m=2305;
eT=0.88;
r=0.322;
f=0.013;
i0=4.55;
ig=[3.5,2.16,1.47,1,0.79];
g=9.8;
CDA=2.4;
n=(n_min:10:n_max);
Tq=(-19.313+295.27.*n-165.44.*n.^2+40.874.*n.^3-3.8445.*n.^4)./1000;
Ft1=Tq.*ig(1).*i0.*eT./r;
Ft2=Tq.*ig(2).*i0.*eT./r;
Ft3=Tq.*ig(3).*i0.*eT./r;
Ft4=Tq.*ig(4).*i0.*eT./r;
Ft5=Tq.*ig(5).*i0.*eT./r;
s=2*pi*r*0.001;
Ua1=n./ig(1)./i0.*60.*s;
Ua2=n./ig(2)./i0.*60.*s;
Ua3=n./ig(3)./i0.*60.*s;
Ua4=n./ig(4)./i0.*60.*s;
Ua5=n./ig(5)./i0.*60.*s;
G=m*g;
ua=(0:1:340);
Fw=CDA.*ua.^2./21.15;
D1=(Ft1-Fw)./G;
D2=(Ft2-Fw)./G;
D3=(Ft3-Fw)./G;
D4=(Ft4-Fw)./G;
D5=(Ft5-Fw)./G;
a11=(g.*(D1-f))./1.22;
a22=(g.*(D2-f))./1.09;
a33=(g.*(D3-f))./1.05;
a44=(g.*(D4-f))./1.03;
a55=(g.*(D5-f))./1.02;
AA1=(1./a11);
AA2=(1./a22);
AA3=(1./a33);
AA4=(1./a44); AA5=(1./a55);
plot(Ua1,AA1,Ua2,AA2,'r',Ua3,AA3,'b',Ua4,AA4,'g',Ua5,AA5,'r')

yang10804050232

不过画出来的图形好像还是不对