您的位置：360文档中心› 汽车理论图形MATLAB程序

汽车理论图形MATLAB程序

汽车理论图形MATLAB程序
汽车理论图形MATLAB程序

功率平衡图

m=1230;g=9.8;

ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75;

r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;

np=6000;Pemax=83;

%绘制汽车驱动力与行驶阻力平衡图

for i=1:56;

n=500:100:6000;

Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3);

Tq(i)=9549*Pe(i)/n(i);

end

for j=1:5

for i=1:56

Ft(i,j)=Tq(i)*ig(j)*i0*yt/r;

ua(i,j)=0.377*r*n(i)/(ig(j)*i0);

Fz(i,j)=m*g*f+CD*A*(ua(i,j)^2)/21.15;

end

end

figure

plot(ua,Ft,ua,Fz);

title('汽车驱动力与行驶阻力平衡图');

xlabel('ua(km/h)');

ylabel('Ft(N)');

text(20,6700,'Ft1');

text(40,4000,'Ft2');

text(50,2800,'Ft3');

text(80,2000,'Ft4');

text(100,1600, 'Ft5');

text(100,800,'Ff+Fw');

for k=1:56;

n=500:100:6000;

Pe(k)=Pemax*(n(k)/np+(n(k)/np)^2-(n(k)/np)^3);

Tq(k)=9549*Pe(k)/n(k);

Ft(k)=Tq(k)*ig(4)*i0*yt/r;

ua(k)=0.377*r*n(k)/(ig(4)*i0);

Fz(k)=m*g*f+CD*A*(ua(k)^2)/21.15;

E(k)=abs((Ft(k)-Fz(k)));

end

[Emin,kmin]=min(E);

Umax=ua(kmin)

加速度曲线图

m=1230;g=9.8;

ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75;

r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;

np=6000;Pemax=83;

deta=[1.33 1.12 1.07 1.05 1.04];

for i=1:56;

n=500:100:6000;

Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3);

Tq(i)=9549*Pe(i)/n(i);

end

for j=1:5

for i=1:56

Ft(i,j)=Tq(i)*ig(j)*i0*yt/r;

ua(i,j)=0.377*r*n(i)/(ig(j)*i0);

Fz(i,j)=m*g*f+CD*A*(ua(i,j)^2)/21.15;

a(i,j)=(Ft(i,j)-Fz(i,j))/(deta(j)*m);

b(i,j)= 1./a(i,j);

end

end

x1=ua(:,1);y1=b(:,1);

x2=ua(:,2);y2=b(:,2);

x3=ua(:,3);y3=b(:,3);

x4=ua(:,4);y4=b(:,4);

x5=ua(:,5);y5=b(:,5);

figure

plot(x1,y1,x2,y2,x3,y3,x4,y4,x5,y5);

axis([0 185 0 15]);

title('汽车行驶加速度倒数曲线');

xlabel('ua(km/h)');

ylabel('1/a');

汽车二挡原地起步加速时间

nT=0.9;r=0.31;f=0.0177;CD=0.31;A=2.2;;i0=3.75;

m=3880;g=9.8;Pemax=83;np=6000;

G=m*g; ig=[3.615 2.053 1.393 1.031 0.837];

nmin=600;nmax=6000;

u1=0.377*r*nmin./ig/i0;

u2=0.377*r*nmax./ig/i0;

deta=[1.33 1.12 1.07 1.05 1.04];

ua=[1:0.01:180];N=length(ua);n=0;Tq=0;Ft=0;inv_a=0*ua;delta=0*ua; Ff=G*f;

Fw=CD*A*ua.^2/21.15;

for i=1:N

k=i;

if ua(i)<=u2(2)

n=ua(i)*(ig(2)*i0/r)/0.377;

Tq=9549*(Pemax*(n/np+(n/np)^2-(n/np)^3))/n;

Ft=Tq*ig(2)*i0*nT/r;

inv_a(i)=(deta(2)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

elseif ua(i)<=u2(3)

n=ua(i)*(ig(3)*i0/r)/0.377;

Tq=9549*(Pemax*(n/np+(n/np)^2-(n/np)^3))/n;

Ft=Tq*ig(3)*i0*nT/r;

inv_a(i)=(deta(3)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

elseif ua(i)<=u2(4)

n=ua(i)*(ig(4)*i0/r)/0.377;

Tq=9549*(Pemax*(n/np+(n/np)^2-(n/np)^3))/n;

Ft=Tq*ig(4)*i0*nT/r;

inv_a(i)=(deta(4)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

else

n=ua(i)*(ig(5)*i0/r)/0.377;

Tq=9549*(Pemax*(n/np+(n/np)^2-(n/np)^3))/n;

Ft=Tq*ig(5)*i0*nT/r;

inv_a(i)=(deta(5)*m)/(Ft-Ff-Fw(i));

delta(i)=0.01*inv_a(i)/3.6;

end

a=delta(1:k);

t(i)=sum(a);

end

plot(t,ua);

axis([0 80 0 100]);

title('汽车2档原地起步换挡加速时间曲线'); xlabel('时间t（s）');

ylabel('速度ua（km/h）');

功率平衡图

m=1230;g=9.8;

ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75;

r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;

np=6000;Pemax=83;

deta=[1.33 1.12 1.07 1.05 1.04];

for j=1:5

for i=1:56;

n=500:100:6000;

Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3);

Tq(i)=9549*Pe(i)/n(i);

ua(i,j)=0.377*r*n(i)/(ig(j)*i0);

Pz(i,j)=(m*g*f*ua(i,j)/3600+CD*A*ua(i,j)^3/76140)/yt;

end

end

plot(ua,Pe,ua,Pz);

title('功率平衡图');

xlabel('ua(km/h)');

ylabel('P(kW)');

text(40,90,'一档Pe');

text(80,90,'二档Pe');

text(110,90,'三档Pe');

text(150,90,'四档Pe');

text(200,90,'五档Pe');

text(100,25,'Pz');

%绘制最高档和次高档的等速百公里油耗曲线

n0=[815,1207,1614,2012,2603,3006,3403,3804];

for j=1:5

for i=1:8

Pd(i)=Pemax*(n0(i)/np+(n0(i)/np)^2-(n0(i)/np)^3);

Td(i)=9549*Pd(i)/n0(i);

u(i,j)=0.377*r*n0(i)/(ig(j)*i0);

end

end

b(1)=1326.8-416.46*Pd(1)+72.379*Pd(1)^2-5.8629*Pd(1)^3+0.17768*Pd(1)^4;

b(2)=1354.7-303.98*Pd(2)+36.657*Pd(2)^2-2.0553*Pd(2)^3+0.043072*Pd(2)^4;

b(3)=1284.4-189.75*Pd(3)+14.524*Pd(3)^2-0.51184*Pd(3)^3+0.0068164*Pd(3)^4;

b(4)=1122.9-121.59*Pd(4)+7.0035*Pd(4)^2-0.18517*Pd(4)^3+0.0018555*Pd(4)^4;

b(5)=1141.0-98.893*Pd(5)+4.4763*Pd(5)^2-0.091077*Pd(5)^3+0.00068906*Pd(5)^4;

b(6)=1051.2-73.714*Pd(6)+2.8593*Pd(6)^2-0.05138*Pd(6)^3+0.00035032*Pd(6)^4;

b(7)=1233.9-84.478*Pd(7)+2.9788*Pd(7)^2-0.047449*Pd(7)^3+0.00028230*Pd(7)^4;

b(8)=1129.7-45.291*Pd(8)+0.71113*Pd(8)^2+0.00075215*Pd(8)^3-0.000038568*Pd(8)^4; u1=u(:,1)';

u2=u(:,2)';

u3=u(:,3)';

u4=u(:,4)';

u5=u(:,5)';

B1=polyfit(u1,b,3);

B2=polyfit(u2,b,3);

B3=polyfit(u3,b,3);

B4=polyfit(u4,b,3);

B5=polyfit(u5,b,3);

for k=1:56

bh(k,1)=polyval(B1,ua(k,1));

bh(k,2)=polyval(B2,ua(k,2));

bh(k,3)=polyval(B3,ua(k,3));

bh(k,4)=polyval(B4,ua(k,4));

bh(k,5)=polyval(B5,ua(k,5));

end

for i=1:5

for k=1:56

Q(k,i)=Pz(k,i)*bh(k,i)/(1.02*ua(k,i)*7.10);

end

end

plot(ua(:,4),Q(:,4),ua(:,5),Q(:,5));

axis([0 185 0 20]);

title('最高档（五档）和次高档（四档）等速百公里油耗图'); xlabel('ua(km/h)');

ylabel('Qs(L/100km)');

text(90,10,'四档');

text(125,10,'五档');

爬坡度图

m=1230;g=9.8;

G=m*g;

ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75;

r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;

np=6000;Pemax=83;

for i=1:56;

n=500:100:6000;

Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3);

Tq(i)=9549*Pe(i)/n(i);

end

for j=1:5

for i=1:56

Ft(i,j)=Tq(i)*ig(j)*i0*yt/r;

ua(i,j)=0.377*r*n(i)/(ig(j)*i0);

Fz(i,j)=m*g*f+CD*A*(ua(i,j)^2)/21.15;

F=Ft-Fz;

I=tan(asin(Ft/G))*100;

end

end

plot(ua,I);

hold on;

title('汽车的爬坡度图');

xlabel('ua(km/h)');

ylabel('i(%)');

text(50,55,'一档');

text(50,35,'二档');

text(100,22,'三档');

text(150,15,'四档');

text(200,13,'五档');

动力因数图

m=1230;g=9.8;

G=m*g;

ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75;

r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2;

np=6000;Pemax=83;

%绘制汽车驱动力与行驶阻力平衡图

for i=1:56;

n=500:100:6000;

Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3);

Tq(i)=9549*Pe(i)/n(i);

end

for j=1:5

for i=1:56

Ft(i,j)=Tq(i)*ig(j)*i0*yt/r;

ua(i,j)=0.377*r*n(i)/(ig(j)*i0);

Fw(i,j)=CD*A*(ua(i,j)^2)/21.15;

F=Ft-Fw;

D=F/G;

end

end

plot(ua,D);

hold on;

axis([0 190 0 0.6]);

title('汽车的动力特性图');

xlabel('ua(km/h)');

ylabel('动力因数)');

text(40,0.5,'一档');

text(60,0.3,'二档');

text(100,0.2,'三档');

text(140,0.1,'四档');

text(140,0.04,'五档');

Matlab图形绘制经典案例

Matlab图形绘制经典案例 1、 三维曲线 >> t=0:pi/50:10*pi; >> plot3(sin(2*t),cos(2*t),t) >> axis square >> grid on

2、一窗口多图形>> t=-2*pi:0.01:2*pi; >> subplot(3,2,1)

>> plot(t,sin(t)) >> subplot(3,2,2) >> plot(t,cos(t)) >> subplot(3,2,3) >> plot(t,tan(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,4) >> plot(t,cot(t)) >> axis([-pi pi -100 100]) >> subplot(3,2,5) >> plot(t,atan(t)) >> subplot(3,2,6) >> plot(t,acot(t))

3、图形样式、标注、题字(也可以利用菜单直接Insert) >> x=0:pi/20:2*pi;

>> plot(x,sin(x),'b-.') >> hold on >> plot(x,cos(x),'r--') >> hold on >> plot(x,sin(x)-1,'g:') >> hold on >> plot(x,cos(x)-1) >> xlabel('x'); >> xlabel('x轴'); >> ylabel('y轴'); >> title('图形样式、标注等'); >> text(pi,sin(pi),'x=\pi'); >> legend('sin(x)','cos(x)','sin(x)-1','cos(x)-1'); >> [x1,y1]=ginput(1) %利用鼠标定位查找线上某点的值x1 = 2.0893 y1 = -0.5000 >> gtext('x=2.5') %鼠标定位放置所需的值在线上

汽车理论课后习题Matlab程序

1.3确定一轻型货车的动力性能（货车可装用 4挡或5挡变速器，任选 其中的一种进行整车性能计算）： 1） 绘制汽车驱动力与行驶阻力平衡图。 2） 求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3） 绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用 2档起步加速行 驶至70km/h 的车速一时间曲线，或者用计算机求汽车用 2档起步加速行驶至 70km/h 的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n 曲线的拟合公式为 19.313 295.27(孟o )165,44(^)2 40-874(^)3 环45為4 式中，Tq 为发动机转矩（N?m ） ;n 为发动机转速（r/min ）。 发动 机的最低转速n min =600r/min,最高转速n max =4000r/min 。 装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m 传动系机械效率 n =0.85 滚动阻力系数 f=0.013 空气阻力系数泌风面积 C D A=2.77m 2 解：Matlab 程序： (1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=[600:10:4000]; Tq=-19.313+295.27*( n/1000)-165.44*(门/1000)八2+40.874*(门/1000)八3-3.8445*( n/10 00).A 4; m=3880;g=9.8; nmi n=600; nm ax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793]; nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;lf=0.218;lw1= 1.798;Iw2=3.598; Ft 仁 Tq*ig(1)*i0* nT/r; Ft2=Tq*ig(2)*i0* nT/r; 主减速器传动比 飞轮转动惯量 二前轮转动惯量 四后轮转动惯量 i 0=5.83 l f =0.218kg?m I w1=1.798kg?nf I w2=3.598kg?m ig(数据如下表) 轴距 质心至前轴距离（满载） 质心高（满载） a=1.974m hg=0.9m

MATLAB绘图功能大全

Matlab绘图 强大的绘图功能是Matlab的特点之一，Matlab提供了一系列的绘图函数，用户不需要过多的考虑绘图的细节，只需要给出一些基本参数就能得到所需图形，这类函数称为高层绘图函数。此外，Matlab 还提供了直接对图形句柄进行操作的低层绘图操作。这类操作将图形的每个图形元素（如坐标轴、曲线、文字等）看做一个独立的对象，系统给每个对象分配一个句柄，可以通过句柄对该图形元素进行操作，而不影响其他部分。 本章介绍绘制二维和三维图形的高层绘图函数以及其他图形控制函数的使用方法，在此基础上，再介绍可以操作和控制各种图形对象的低层绘图操作。 一、二维绘图 二维图形是将平面坐标上的数据点连接起来的平面图形。可以采用不同的坐标系，如直角坐标、对数坐标、极坐标等。二维图形的绘制是其他绘图操作的基础。 （一）绘制二维曲线的基本函数 在Matlab中，最基本而且应用最为广泛的绘图函数为plot，利用它可以在二维平面上绘制出不同的曲线。 1．plot函数的基本用法

plot函数用于绘制二维平面上的线性坐标曲线图，要提供一组x 坐标和对应的y坐标，可以绘制分别以x和y为横、纵坐标的二维曲线。plot函数的应用格式 plot(x,y) 其中x,y为长度相同的向量，存储x坐标和y坐标。 例51 在[0 , 2pi]区间，绘制曲线 程序如下：在命令窗口中输入以下命令 >> x=0:pi/100:2*pi; >> y=2*exp(-0.5*x).*sin(2*pi*x); >> plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 注意：指数函数和正弦函数之间要用点乘运算，因为二者是向量。 例52 绘制曲线 这是以参数形式给出的曲线方程，只要给定参数向量，再分别求出x,y向量即可输出曲线： >> t=-pi:pi/100:pi; >> x=t.*cos(3*t); >> y=t.*sin(t).*sin(t); >> plot(x,y) 程序执行后，打开一个图形窗口，在其中绘制出如下曲线 以上提到plot函数的自变量x,y为长度相同的向量，这是最常见、最基本的用法。实际应用中还有一些变化。

汽车理论课后习题Matlab程序

1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选 其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。 2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n 曲线的拟合公式为 234 19.313295.27()165.44()40.874() 3.8445()1000100010001000 q n n n n T =-+-+- 式中，Tq 为发动机转矩（N?m）;n 为发动机转速（r/min ）。 发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。 装载质量 2000kg 整车整备质量 1800kg 总质量 3880kg 车轮半径 0.367m 传动系机械效率 ηt =0.85 滚动阻力系数 f =0.013 空气阻力系数×迎风面积 C D A =2.77m 2 主减速器传动比 i 0=5.83 飞轮转动惯量 I f =0.218kg?m 2 二前轮转动惯量 I w1=1.798kg ?m 2 四后轮转动惯量 I w2=3.598kg?m 2

变速器传动比ig(数据如下表) 轴距L=3.2m 质心至前轴距离（满载）a=1.974m 质心高（满载）hg=0.9m 解：Matlab程序： (1) 求汽车驱动力与行驶阻力平衡图和汽车最高车速程序： n=[600:10:4000]; Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3 .8445*(n/1000).^4; m=3880;g=9.8;nmin=600;nmax=4000; G=m*g; ig=[5.56 2.769 1.644 1.00 0.793];nT=0.85;r=0.367;f=0.013;CDA=2.77;i0=5.83; L=3.2;a=1.947;hg=0.9;If=0.218;Iw1=1.798;Iw2=3.598; Ft1=Tq*ig(1)*i0*nT/r; Ft2=Tq*ig(2)*i0*nT/r; Ft3=Tq*ig(3)*i0*nT/r; Ft4=Tq*ig(4)*i0*nT/r;

第二讲 MatLab图形绘制功能Final

第二讲 MatLab图形绘制功能 教学目的：熟练掌握二维曲线作图；了解隐函数绘图；掌握三维曲线和三维曲面作图的方法和技巧；了解MATLAB图形修饰处理方法；学会一些简单动画制作手段． 教学重点：二维曲线的作图, 极坐标和参数作图; 三维曲线和曲面作图. 教学难点：点运算在作图中的应用; 曲线和曲面作图的差别; 隐函数作图.通过解说，演示以及举例子克服难点难关． 一、二维平面图形 表1基本绘图函数

hold on 命令用于在已画好的图形上添加新的图形. plot是绘制一维曲线的基本函数，但在使用此函数之前，我们需先定义曲线上每一点的x及y 坐标.下例可画出一条正弦曲线： x=0:0.001:10; % 0到10的1000个点的x坐标 y=sin(x); % 对应的y坐标 plot(x,y); % 绘图 Y=sin(10*x);plot(x,y,'r:',x,Y,'b') % 同时画两个函数

若要改变颜色，在坐标对后面加上相关字串即可： x=0:0.01:10; plot(x,sin(x),'r') 若要同时改变颜色及图线型态（Line style），也是在坐标对后面加上相关字串即可：plot(x,sin(x),'r*') 用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围：

axis([0,6,-1.5,1]) MATLAB 也可对图形加上各种注解与处理： xlabel('x 轴'); % x 轴注解 ylabel('y 轴'); % y 轴注解 title('余弦函数'); % 图形标题 legend('y = cos(x)'); % 图形注解 gtext('y = cos(x)'); % 图形注解 ,用鼠标定位注解位置 grid on; % 显示格线 fplot 的指令可以用来自动的画一个已定义的函数分布图，而无须产生绘图所需要的一组数据作为变数.其语法为fplot('fun',[xmin xmax ymin ymax])，其中 fun 为一已定义的函数名称，例如 sin , cos 等等；而 xmin , xmax , ymin , ymax 则是设定绘图横轴及纵轴的下限及上限. 以下的例子是将一函数 f(x)=sin(x)/x 在-20> fplot('sin(x)./x',[-20 20 -0.4 1.2]) 例 画椭圆12 322 22=+y x . a = [0:pi/50:2*pi]'; %角度 π20-

汽车理论图形MATLAB程序

功率平衡图 m=1230;g=9.8; ig=[3.615 2.053 1.393 1.031 0.837]; i0=3.75; r=0.31;yt=0.9;f=0.017;CD=0.31;A=2.2; np=6000;Pemax=83; %绘制汽车驱动力与行驶阻力平衡图 for i=1:56; n=500:100:6000; Pe(i)=Pemax*(n(i)/np+(n(i)/np)^2-(n(i)/np)^3); Tq(i)=9549*Pe(i)/n(i); end for j=1:5 for i=1:56 Ft(i,j)=Tq(i)*ig(j)*i0*yt/r; ua(i,j)=0.377*r*n(i)/(ig(j)*i0); Fz(i,j)=m*g*f+CD*A*(ua(i,j)^2)/21.15; end end figure plot(ua,Ft,ua,Fz); title('汽车驱动力与行驶阻力平衡图'); xlabel('ua(km/h)'); ylabel('Ft(N)'); text(20,6700,'Ft1'); text(40,4000,'Ft2'); text(50,2800,'Ft3'); text(80,2000,'Ft4'); text(100,1600, 'Ft5'); text(100,800,'Ff+Fw'); for k=1:56; n=500:100:6000; Pe(k)=Pemax*(n(k)/np+(n(k)/np)^2-(n(k)/np)^3); Tq(k)=9549*Pe(k)/n(k); Ft(k)=Tq(k)*ig(4)*i0*yt/r; ua(k)=0.377*r*n(k)/(ig(4)*i0); Fz(k)=m*g*f+CD*A*(ua(k)^2)/21.15; E(k)=abs((Ft(k)-Fz(k))); end [Emin,kmin]=min(E); Umax=ua(kmin)

MatLab图形绘制功能

MatLab & 数学建模 第二讲 MatLab图形绘制功能 一、二维平面图形 基本绘图函数 hold on 命令用于在已画好的图形上添加新的图形 plot是绘制一维曲线的基本函数，但在使用此函数之前，我们需先定义曲线上每一点的x及y座标。下例可画出一条正弦曲线： x=0:0.001:10; % 0到10的1000个点的x座标

y=sin(x); % 对应的y座标 plot(x,y); % 绘图 Y=sin(10*x); plot(x,y,'r:',x,Y,'b') % 同时画两个函数 ?若要改变颜色，在座标对後面加上相关字串即可：x=0:0.01:10;

plot(x,sin(x),'r') 若要同时改变颜色及图线型态（Line style），也是在坐标对後面加上相关字串即可： plot(x,sin(x),'r*') 用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围 axis([0,6,-1.5,1])

MATLAB也可对图形加上各种注解与处理： xlabel('x轴'); % x轴注解 ylabel('y轴'); % y轴注解 title('余弦函数'); % 图形标题 legend('y = cos(x)'); % 图形注解 gtext('y = cos(x)'); % 图形注解 ,用鼠标定位注解位置 grid on; % 显示格线 fplot的指令可以用来自动的画一个已定义的函数分布图，而无须产生绘图所须 要的一组数据做为变数。其语法为fplot('fun',[xmin xmax ymin ymax])，其中fun 为一已定义的函数名称，例如sin, cos等等；而xmin, xmax, ymin, ymax则是设定绘图横轴及纵轴的下限及上限。 以下的例子是将一函数 f(x)=sin(x)/x 在-20> fplot('sin(x)./x',[-20 20 -0.4 1.2])

汽车理论1.3和2.7matlab编程答案

孙野 20081268 1.3（1）绘制汽车驱动力与行驶阻力平衡图 选用5挡变速器进行整车性能计算 发动机转速与汽车行驶速度之间的关系：0 377.0i i rn u g a = 发动机使用外特性的T q -n 曲线拟合公式： 4 32)1000 (8445.3)1000(874.40)1000(44.165)1000(27.295313.19n n n n T q -+-+-= 汽车驱动力：r i i T F T g q t η0= 行驶阻力：2 215 .2115.21a D a D j w i f u A C mgf dt du m u A C Gi Gf F F F F +=++ +=+++δ 用matlab 编写程序如下： m1=2000; m2=1800; mz=3880; g=9.81; r=0.367; CdA=2.77; f=0.013; nT=0.85; ig=[5.56 2.769 1.644 1.00 0.793]; i0=5.83; If=0.218; Iw1=1.798; Iw2=3.598; %输入已知参数 Iw=2*Iw1+4*Iw2; for i=1:69 n(i)=(i+11)*50; Ttq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000)^2+40.874*(n(i)/1000)^3-3.8445*(n(i)/1000)^4; end %从600~4000rpm 油拟合公式计算发动机转 距 for j=1:5 for i=1:69 Ft(i,j)=Ttq(i)*ig(j)*i0*nT/r; ua(i,j)=0.377*r*n(i)/(ig(j)*i0); Fz(i,j)=CdA*ua(i,j)^2/21.15+mz*g*f;

汽车理论课后题matlab程序

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*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000) .^4; %驱动力行驶阻力平衡图 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; plot(Ua,Ft); hold on; end Ff=G*f; ua=0:0.1:max(Ua); Fw=CdA*ua.^2/21.15; plot(ua,(Ff+Fw)); title('驱动力-行驶阻力平衡图'); xlabel('Ua/(km/h)');ylabel('Ft/N'); gtext('Ft1'),gtext('Ft2'),gtext('Ft3'),gtext('Ft4'),gtext('Ft5'),gtext('Ff+Fw') [x,y]=ginput(1); disp('汽车的最高车速');disp(x);disp('km/h'); %最大爬坡度及最大爬坡度时的附着率 Ua=0.377*r*n/5.56/i0; Ft=Ttq*5.56*i0*eff/r; Fw=CdA*Ua.^2/21.15; i=tan(asin((Ft-(Ff+Fw))/G)); disp('汽车的最大爬坡度');disp(max(i)); C=max(i)/(a/L+hg/L*max(i)); disp('克服最大爬坡度时的附着率');disp(C); %加速度倒数曲线 figure; for ig=[5.56,2.769,1.644,1.00,0.793] Ua=0.377*r*n/ig/i0; q=1+Iw/(m*r^2)+If*ig^2*i0^2*eff/(m*r^2); Ft=Ttq*ig*i0*eff/r; Fw=CdA*Ua.^2/21.15; as=(Ft-(Ff+Fw))/q/m; plot(Ua,1./as); hold on; end axis([0 98 0 10]); title('行驶加速度倒数曲线');xlabel('Ua/(km/h)');ylabel('1/a'); gtext('1/a1'),gtext('1/a2'),gtext('1/a3'),gtext('1/a4'),gtext('1/a5'); %加速时间曲线

MATLAB在“汽车理论”中的简单应用

MATLAB在 汽车理论 中的简单应用 汽车理论 可以使用MATLAB作为科学计算工具的课程三本文档的目的在于向之前没有接触过MATLAB的同学介绍其在 汽车理论 中的简单应用,更多关于MATLAB的知识需要大家结合其他课程,通过进一步的学习和研究来掌握三 需要说明的是: 1)本文档默认读者已掌握C语言和线性代数基本知识三 2)红色字体的语句可以直接复制到MATLAB中进行运算三语句中用到的符号(引号内部分除外)均为英文状态下输入三本文档所有语句均已在MATLAB7.8.0(R2009a)版本下调试过三 3)例程编写样式纯属个人习惯,大家不必拘泥三 4)如对本文档及‘汽车理论“第6版中的附录C一一个学生的 汽车理论 课程MATLAB习题编程思路有任何问题,欢迎发送邮件一起讨论三邮件地址:wang-d07@ https://www.360docs.net/doc/1b382357.html,三 1.MATLAB使用指南 打开MATLAB后出现的是软件的主界面三几个窗口中最重要的是Command Window,用户可以在里面输入程序,编写的函数的结果也在这里显示三 首先对MATLAB进行配置三如果是在自己的计算机上编程,建议先建一个属于自己的workspace,用来保留用户的工作环境三选择 file Save Workspace As ,在合适的目录下保存用户自己的workspace(.mat文件)三这样,以后用户编程的环境就默认为这个workspace了三 虽然可以直接在Command Window里编写程序,但这样的程序无法被保存和发布三选择 file New Blank M File ,可以新建空白的MATLAB函数文件(.m文件),在这里编写的程序可以被保存和发布三编写好的程序要经过编译之后才能执行三选择 Tools Save File and Run (也可以直接按键),进行程序的编译和执行,结果显示在Command Window中三如果程序有错误,MATLAB也会提示程序运行到哪一行命令终止了,并指出可能的错误类型,然后用户就开始了漫长的debug三 如果编写了子程序和主程序,那么需要先对子程序进行编译(此时不会产生任何结果),只有编译通过后再对主程序进行编译,才会显示结果三 当然,也可以选择新建 Function M File ,它直接提供一个函数模板三(个人感觉没有太大意义) 明白该如何操作MATLAB后,接下来介绍MATLAB中与C语言不同的变量单位 矩阵三 2.矩阵及其运算 MATLAB全称是Matrix Laboratory(矩阵实验室),因此矩阵是MATLAB中最基本的运算单位,熟练掌握矩阵的知识对于灵活运用MATLAB有很大帮助三 1

汽车理论1-4MATLAB编程

汽车理论1-4MATLAB 编程 0102030405060708090100 5000 10000 15000 汽车驱动力与阻力平衡图 u a /km.h -1 F /N 1020 3040 5060708090100 024 68 10 12 14 加速度倒数-速度曲线图 u 1/a

010203040 5060708090100 10 20 30 40 50 60 u/(km/h) P /k W 汽车功率平衡图 10 203040 5060708090100 1214 16 18 20 22 24 最高档等速百公里油耗曲线 Ua/(km/h) Q s /L

23 2425 2627 2829 1213 14 15 16 17 18 燃油积极性-加速时间曲线 燃油经济性(qs/L) 动力性--原地起步加速时间 (s t /s ) 源程序： 《第一章》 m=3880; g=; r=; x=; f=; io=; CdA=; If=; Iw1=; Iw2=; Iw=Iw1+Iw2; ig=[ ]; %变速器传动比 L=; a=; hg=;

n=600:1:4000; T=+*n/*(n/1000).^2+*(n/1000).^*(n/1000).^4; Ft1=T*ig(1)*io*x/r;%计算各档对应转速下的驱动力 Ft2=T*ig(2)*io*x/r; Ft3=T*ig(3)*io*x/r; Ft4=T*ig(4)*io*x/r; u1=*r*n/(io*ig(1)); u2=*r*n/(io*ig(2)); u3=*r*n/(io*ig(3)); u4=*r*n/(io*ig(4)); u=0:130/3400:130; F1=m*g*f+CdA*u1.^2/;%计算各档对应转速下的驱动阻力 F2=m*g*f+CdA*u2.^2/; F3=m*g*f+CdA*u3.^2/; F4=m*g*f+CdA*u4.^2/; figure(1); plot(u1,Ft1,'-r',u2,Ft2,'-m',u3,Ft3,'-k',u4,Ft4,'-b',u1,F1,'-r',u2,F2,'-m',u3,F3,'-k',u4,F4,'-b','LineWidth',2) title('汽车驱动力与阻力平衡图'); xlabel('u_{a}/^{-1}') ylabel('F/N') gtext('F_{t1}') gtext('F_{t2}') gtext('F_{t3}') gtext('F_{t4}') gtext('F_{f}+F_{w}') %由汽车驱动力与阻力平衡图知，他们无交点，u4在最大转速时达到最大 umax=u4(3401) Ft1max=max(Ft1); imax=(Ft1max-m*g*f)/(m*g) disp('假设是后轮驱动'); C=imax/(a/L+hg*imax/L) % 附着率 delta1=1+(Iw1+Iw2)/(m*r^2)+If*ig(1)*r^2*io^2*x/(m*r^2); delta2=1+(Iw1+Iw2)/(m*r^2)+If*ig(2)*r^2*io^2*x/(m*r^2); delta3=1+(Iw1+Iw2)/(m*r^2)+If*ig(3)*r^2*io^2*x/(m*r^2); delta4=1+(Iw1+Iw2)/(m*r^2)+If*ig(4)*r^2*io^2*x/(m*r^2); a1=(Ft1-F1)/(delta1*m); %加速度 a2=(Ft2-F2)/(delta2*m); a3=(Ft3-F3)/(delta3*m); a4=(Ft4-F4)/(delta4*m); h1=1./a1; %加速度倒数 h2=1./a2; h3=1./a3; h4=1./a4;

汽车理论课后习题MATLAB编程

汽车理论MATLAB编程 （第1-4章） 第一章 m=3880; g=9.8; r=0.367; x=0.85; f=0.013; io=5.83; CdA=2.77; If=0.218; Iw1=1.798; Iw2=3.598; Iw=Iw1+Iw2; ig=[6.09 3.09 1.71 1.00]; %变速器传动比 L=3.2; a=1.947; hg=0.9; n=600:1:4000; T=-19.313+295.27*n/1000-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1000).^4; Ft1=T*ig(1)*io*x/r;%计算各档对应转速下的驱动力 Ft2=T*ig(2)*io*x/r; Ft3=T*ig(3)*io*x/r; Ft4=T*ig(4)*io*x/r; u1=0.377*r*n/(io*ig(1)); u2=0.377*r*n/(io*ig(2)); u3=0.377*r*n/(io*ig(3)); u4=0.377*r*n/(io*ig(4)); u=0:130/3400:130; F1=m*g*f+CdA*u1.^2/21.15;%计算各档对应转速下的驱动阻力 F2=m*g*f+CdA*u2.^2/21.15; F3=m*g*f+CdA*u3.^2/21.15; F4=m*g*f+CdA*u4.^2/21.15; figure(1); plot(u1,Ft1,'-r',u2,Ft2,'-m',u3,Ft3,'-k',u4,Ft4,'-b',u1,F1,'-r',u2,F2,'-m',u3,F3,'-k',u4,F4,'-b','LineWidth',2) title('汽车驱动力与阻力平衡图'); xlabel('u_{a}/km.h^{-1}') ylabel('F/N') gtext('F_{t1}') gtext('F_{t2}')

MatLab图形绘制功能

MatLab图形绘制功能MatLab图形绘制功能例子解读 一、二维平面图形 基本绘图函数 命令含义 plot 建立向量或矩阵各队队向量的图形 loglog x、y轴都取对数标度建立图形 semilogx x轴用于对数标度，y轴线性标度绘制图形semilogy y轴用于对数标度，x轴线性标度绘制图形title 给图形加标题 xlabel 给x轴加标记 ylabel 给y轴加标记 text 在图形指定的位置上加文本字符串 gtext 在鼠标的位置上加文本字符串 grid 打开网格线 plot绘图函数的叁数 字元颜色字元图线型态 y 黄色 . 点 k 黑色 o 圆 w 白色 x x b 蓝色 + + g 绿色 * * r 红色 - 实线

c 亮青色 : 点线 m 锰紫色 -. 点虚线 -- 虚线 hold on 命令用于在已画好的图形上添加新的图形 plot是绘制一维曲线的基本函数，但在使用此函数之前，我们需先定义曲线上每一点的x及y座标。下例可画出一条正弦曲线: x=0:0.001:10; % 0到10的1000个点的x座标 y=sin(x); % 对应的y座标 plot(x,y); % 绘图 Y=sin(10*x); plot(x,y,'r:',x,Y,'b') % 同时画两个函数 , 若要改变颜色，在座标对後面加上相关字串即可: x=0:0.01:10; plot(x,sin(x),'r') 若要同时改变颜色及图线型态(Line style)，也是在坐标对後面加上相关字串 即可: plot(x,sin(x),'r*') 用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围 axis([0,6,-1.5,1]) MATLAB也可对图形加上各种注解与处理: xlabel('x轴'); % x轴注解 ylabel('y轴'); % y轴注解 title('余弦函数'); % 图形标题 legend('y = cos(x)'); % 图形注解 gtext('y = cos(x)'); % 图形注解 ,用鼠标定位注解位置 grid on; % 显示 格线

汽车理论课后题matlab程序

汽车理论课后题matlab程序

————————————————————————————————作者：————————————————————————————————日期: ?

1.3 n＝600：1：４000; r=0．367; i０=５.83; eff＝0.85; f＝0.０13； ｍ＝3880；g=9.8; G=m*g; CdA=2.77；ａ=１.94７; hg=０.9；L＝3.2； Iw1=1．79８; Iw２＝3.598; Ｉｗ=Iw１+Iw２; Iｆ=0.218； Tｔｑ＝-19.313+295.２7*n/1000－1６5.44*（ｎ／10０0).^2+40.874*(n/1０00).^3－３.８445*（ｎ/１00０)．^4； %驱动力行驶阻力平衡图 for ｉｇ=[5.56,2.７69,１.644,１.０0,0.793] Uａ=0.377*r＊n/ig/ｉ0； Ft=Tｔq*ｉg＊i０*eff/r； ｐlｏt(Ua,Ft)； ｈold on; ｅｎｄ Ff=G*f; ua=0：０.1:max（Uａ); Ｆｗ＝CdＡ*ｕa．^２/２1.15; plot(ua，(Ff+Ｆw)); title('驱动力-行驶阻力平衡图'）; xlａbel(＇Ua／(ｋm/ｈ）'）;ylab ｅl(＇Fｔ/Ｎ'）; gtext('Ｆt1')，gｔext('Fｔ２'),ｇｔext（'Ft3')，gｔexｔ('Ft4＇），ｇteｘt('Fｔ5')，gtext(＇Fｆ+Fw') [x，y]=ginput(1); disp('汽车的最高车速')；disp(x);dｉsp(＇km/h'); %最大爬坡度及最大爬坡度时的附着率 Ua=0.37７*r*ｎ／5.56/ｉ０； Ft=Ｔtq*５.56*i0*eff/r; Ｆw=CｄA*Uａ.^２/21.15； i=tａｎ（asiｎ((Ft-(Ｆf+Fｗ))/Ｇ)); disp（'汽车的最大爬坡度');diｓp（max(ｉ）); Ｃ=maｘ(i)/(a/L＋hg/L＊maｘ(i))； diｓp（＇克服最大爬坡度时的附着率'）;disｐ(C); %加速度倒数曲线 fｉguｒe; ｆoｒ ig=[5.56,２.７６９,1.644，1．00，0．７93] Ua=0．３７７*ｒ*ｎ/ig／i0; q=1+Iｗ/(m*ｒ＾2)+If＊ig^2*i０^２＊eff/(m＊r^2); Ft=Tｔq*ig*i0*ｅfｆ／r; Ｆw=CdA*Ua.^2／２1.15; ａs=(Ｆt-（Fｆ+Ｆw）)/q/ｍ; ｐlｏt(Ua,１．/as)； hold on; ｅnd axiｓ([0 98 0 １0]); titlｅ(＇行驶加速度倒数曲线')；ｘlabel(＇Ua/(kｍ/h）');ylａｂel（'１/a')； gteｘｔ('１/a1'),gｔｅxt('１/a2'),ｇtext('1／ａ3＇),ｇtext('１/a4')，ｇtext('1/

MatLab图形绘制功能口令及代码

第二讲MatLab图形绘制功能 一、二维平面图形 基本绘图函数 hold on 命令用于在已画好的图形上添加新的图形 plot是绘制一维曲线的基本函数，但在使用此函数之前，我们需先定义曲线上每一点的x及y座标。下例可画出一条正弦曲线： x=0:0.001:10; % 0到10的1000个点的x座标

y=sin(x); % 对应的y座标 plot(x,y); % 绘图 Y=sin(10*x); plot(x,y,'r:',x,Y,'b') % 同时画两个函数 若要改变颜色，在座标对後面加上相关字串即可：x=0:0.01:10;

plot(x,sin(x),'r') 若要同时改变颜色及图线型态（Line style），也是在坐标对後面加上相关字串即可： plot(x,sin(x),'r*') 用axis([xmin,xmax,ymin,ymax])函数来调整图轴的范围 axis([0,6,-1.5,1])

MATLAB也可对图形加上各种注解与处理： xlabel('x轴'); % x轴注解 ylabel('y轴'); % y轴注解 title('余弦函数'); % 图形标题 legend('y = cos(x)'); % 图形注解 gtext('y = cos(x)'); % 图形注解 ,用鼠标定位注解位置 grid on; % 显示格线 fplot的指令可以用来自动的画一个已定义的函数分布图，而无须产生绘图所须 要的一组数据做为变数。其语法为fplot('fun',[xmin xmax ymin ymax])，其中fun 为一已定义的函数名称，例如sin, cos等等；而xmin, xmax, ymin, ymax则是设定绘图横轴及纵轴的下限及上限。 以下的例子是将一函数 f(x)=sin(x)/x 在-20> fplot('sin(x)./x',[-20 20 -0.4 1.2])

汽车理论matlab作业

确定一轻型货车的动力性能。 绘制汽车行驶加速度倒数曲线；用计算机求汽车用□档起 的加速时间。 (1)绘制汽车驱动力与行驶阻力平衡图 m1=2000; m2=1800; mz=3880; g=9.81; r=0.367; CdA=2.77; f=0.013; nT=0.85; ig=[5.56 2.769 1.644 1.00 0.793]; i0=5.83; If=0.218; Iw1=1.798; Iw2=3.598; Iw=2*lw1+4*lw2; for i=1:69 n (i)=(i+11)*50; Ttq(i)=-19.313+295.27*( n(i)/1000)-165.44*( n(i)/1000)人2+40.874*( n(i)/1000)人3-3.8445* (n (i)/1000 )人4; end for j=1:5 for i=1:69 Ft(i,j)=Ttq(i)*ig(j)*i0* nT/r; ua(i,j)=0.377*r* n(i)/(ig(j)*i0); Fz(i,j)=CdA*ua(i,j)A2/21.15+mz*g*f; end end pl ot(ua,Ft,ua,Ff,ua,Ff+Fw) title('汽车驱动力与行驶阻力平衡图 '); xlabel('ua(km/h)'); ylabel('Ft(N)'); gtext( 'Ft1') 1) 绘制汽车驱动力与行驶阻力平衡图; 2) 求汽车最高车速与最大爬坡度; 3) 步加速行驶至 70km/h 所需 已知数据略。 (参见《汽车理论》习题第一章第 3题) 解题程序如下 :用Matlab 语言

Matlab绘图功能

第五章Matlab绘图功能 5.1 二维图形的绘制 5.1.1 常用的二维图形绘图函数 基本的二维绘图函数有 plot ——绘制2维曲线； title ——给图形加标题； grid ——显示网格线； xlabel ——给x轴加标记； ylabel ——给y轴加标记； text ——在坐标图中加入文字注释。 π的曲线图。 例：画出函数x =，其中x从0到π2步进100 yπ2 sin / X=0:pi/100:2*pi; Y=sin(X); plot(X,Y); % 作图 grid on; % 网格线显示，若该为grid off则不显示网格 ylabel('y=sin 2\pi x'); % Y轴标注，可以有汉字 xlabel('x'); % X轴标注，可以有汉字 title('function plot y=sin 2\pi x'); % 图标题 text(0.5,sin(0.5),'\leftarrow sin 2 \pi 0.5'); % text()可以在指定坐标处写文字标注 text(2.3,sin(2.3),'\leftarrow sin 2 \pi 2.3'); % 所有标注中均可使用汉字 % 对于特殊符号，如希腊字母，箭头等需要采用LaTeX格式 结果如图5.1 所示。

图5.1 基本的二维绘图函数用法 5.1.2 图形的线型和颜色控制 在命令plot的高级用法中，可以设置作图的线型，标记类型，线和标记的颜色，粗细等特征。用命令doc LineSpec和doc plot可以查询详细的帮助文档。 线型的定义如下： - solid line (default) 实线 -- dashed line 虚线 : dotted line 虚点连线 -. dash-dot line 点划线 常用标记的定义为： + plus sign 十字标记 o circle 小圈标记 * asterisk 星号标记 . point 黑点标记

汽车理论matlab

汽车理论 1.3 确定一轻型货车的动力性能（货车可装用4挡或5挡变速器，任选其中的一种进行整车性能计算）： 1）绘制汽车驱动力与行驶阻力平衡图。 2）求汽车最高车速，最大爬坡度及克服该坡度时相应的附着率。 3）绘制汽车行驶加速度倒数曲线，用图解积分法求汽车用2档起步加速行驶至70km/h 的车速－时间曲线，或者用计算机求汽车用2档起步加速行驶至70km/h 的加速时间。 轻型货车的有关数据： 汽油发动机使用外特性的Tq-n 曲线的拟合公式为 234 19.313295.27()165.44()40.874() 3.8445()1000100010001000 q n n n n T =-+-+- 式中，Tq 为发动机转矩（N ?m ）;n 为发动机转速（r/min ）。 发动机的最低转速n min =600r/min,最高转速n max =4000r/min 。

装载质量2000kg 整车整备质量1800kg 总质量3880kg 车轮半径0.367m 传动系机械效率ηt=0.85 滚动阻力系数f=0.013 空气阻力系数×迎风面积C D A=2.77m2 主减速器传动比i0=5.83 飞轮转动惯量I f=0.218kg?m2 二前轮转动惯量I w1=1.798kg?m2 四后轮转动惯量I w2=3.598kg?m2 变速器传动比ig(数据如下表) Ⅰ档Ⅱ档Ⅲ档Ⅳ档Ⅴ档四档变速器 6.09 3.09 1.71 1.00 - 五档变速器 5.56 2.769 1.644 1.00 0.793 轴距L=3.2m 质心至前轴距离（满载）a=1.974m 质心高（满载）hg=0.9m MATlab程序如下： n=600:4000; r=0.367; nt=0.85; f=0.013; CA=2.77; io=5.83; m=3880; g=9.8; Tq=-19.313+295.27.*n./1000-165.44.*(n./1000).^2+40.874.*(n./1000).^3-3.8445.*(n./1000).^4; ig1=5.56; n= linspace(600,4000,40); Ua1=0.377.*r.*n./(ig1.*io); Ft1=(-19.313+295.27.*n./1000-165.44.*(n./1000).^2+40.874.*(n./1000).^3-3.8445.*(n./1000).^4). *ig1.*io.*nt./r; plot(Ua1,Ft1); ig2=2.769; Ua2=0.377.*r.*n./(ig2.*io); Ft2=(-19.313+295.27.*n./1000-165.44.*(n./1000).^2+40.874.*(n./1000).^3-3.8445.*(n./1000).^4). *ig2.*io.*nt./r; ig3=1.644; Ua3=0.377.*r.*n./(ig3.*io); Ft3=(-19.313+295.27.*n./1000-165.44.*(n./1000).^2+40.874.*(n./1000).^3-3.8445.*(n./1000).^4). *ig3.*io.*nt./r; ig4=1.00; Ua4=0.377.*r.*n./(ig4.*io); Ft4=(-19.313+295.27.*n./1000-165.44.*(n./1000).^2+40.874.*(n./1000).^3-3.8445.*(n./1000).^4). *ig4.*io.*nt./r; ig5=0.793; Ua5=0.377.*r.*n./(ig5.*io); Ft5=(-19.313+295.27.*n./1000-165.44.*(n./1000).^2+40.874.*(n./1000).^3-3.8445.*(n./1000).^4). *ig5.*io.*nt./r; Fr=m.*g.*f+(CA/21.15).*Ua5.^2;

相关主题
相关文档
最新文档