mathlab汽车动力匹配最小传动比的优化程序.docx
《mathlab汽车动力匹配最小传动比的优化程序.docx》由会员分享，可在线阅读，更多相关《mathlab汽车动力匹配最小传动比的优化程序.docx（29页珍藏版）》请在冰豆网上搜索。
mathlab汽车动力匹配最小传动比的优化程序
1、利用C曲线进行最小传动比的优化：
（1）加速时间：
；
而：
；；；；；
（1.1）加速度-速度曲线
1档：
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=8.4;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.22;
a=1/（s*m）*（Ft-Ff-Fw）;
p=polyfit（Ua,a,3）;
Ua1=3:
0.5:
11;
a1=polyval（p,Ua1）;
>>plot（Ua,a,'o',Ua1,a1,'-'）
2档：
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=5.516;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.09;
a=1/（s*m）*（Ft-Ff-Fw）;
p=polyfit（Ua,a,3）;
Ua1=5:
0.5:
17;
a1=polyval（p,Ua1）;
>>plot（Ua,a,'o',Ua1,a1,'-'）
3档：
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=3.616;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.04;
a=1/（s*m）*（Ft-Ff-Fw）;
p=polyfit（Ua,a,3）;
Ua1=8:
0.5:
26;
a1=polyval（p,Ua1）;
>>plot（Ua,a,'o',Ua1,a1,'-'）
4档：
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=2.357;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.02;
a=1/（s*m）*（Ft-Ff-Fw）;
p=polyfit（Ua,a,3）;
Ua1=12:
0.5:
40;
a1=polyval（p,Ua1）;
>>plot（Ua,a,'o',Ua1,a1,'-'）
5档：
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=1.539;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.01;
a=1/（s*m）*（Ft-Ff-Fw）;
p=polyfit（Ua,a,3）;
Ua1=18:
0.5:
60;
a1=polyval（p,Ua1）;
>>plot（Ua,a,'o',Ua1,a1,'-'）
6档：
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=1;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.005;
a=1/（s*m）*（Ft-Ff-Fw）;
p=polyfit（Ua,a,3）;
Ua1=25:
0.5:
92;
a1=polyval（p,Ua1）;
>>plot（Ua,a,'o',Ua1,a1,'-'）
（1.2）加速度倒数-速度曲线
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=8.4;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.22;
a=1/（s*m）*（Ft-Ff-Fw）;
b=a.^（-1）;
p=polyfit（Ua,b,3）;
Ua1=3:
0.5:
11;
b1=polyval（p,Ua1）;
>>plot（Ua1,b1）
holdon;
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=5.516;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.09;
a=1/（s*m）*（Ft-Ff-Fw）;
b=a.^（-1）;
p=polyfit（Ua,b,3）;
Ua1=5:
0.5:
17;
b1=polyval（p,Ua1）;
>>plot（Ua1,b1）
holdon;
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=3.616;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.04;
a=1/（s*m）*（Ft-Ff-Fw）;
b=a.^（-1）;
p=polyfit（Ua,b,3）;
Ua1=8:
0.5:
26;
b1=polyval（p,Ua1）;
>>plot（Ua1,b1）
holdon;
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=2.357;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.02;
a=1/（s*m）*（Ft-Ff-Fw）;
b=a.^（-1）;
p=polyfit（Ua,b,3）;
Ua1=12:
0.5:
40;
b1=polyval（p,Ua1）;
>>plot（Ua1,b1）
holdon;
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=1.539;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.01;
a=1/（s*m）*（Ft-Ff-Fw）;
b=a.^（-1）;
p=polyfit（Ua,b,3）;
Ua1=18:
0.5:
60;
b1=polyval（p,Ua1）;
>>plot（Ua1,b1）
holdon;
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=1;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*T;
n=[800,1000,1200,1400,1500,1600,1800,2000,2200,2400];
Ua=（0.377*r/（ig*i0））*n;
G=225400;
F=0.0076+0.000056*Ua;
C=0.9;
A=7;
Fw=C*A*Ua.^2/21.15;
m=23000;
s=1.005;
a=1/（s*m）*（Ft-Ff-Fw）;
b=a.^（-1）;
p=polyfit（Ua,b,3）;
Ua1=25:
0.5:
92;
b1=polyval（p,Ua1）;
>>plot（Ua1,b1）
（2）i0=5.476时的加速时间
closeall
T=[556,661,742,793,796,788,773,754,725,685];
r=0.545;
ig=8.4;
i0=5.476;
j=0.894;
Ft=（ig*i0*j/r）*
升级会员 