MATLAB知识的简单运用.doc

实验04051001 MATLAB语言基础操作成绩报告成绩实验目的1) 熟悉MATLAB的运行环境2) 掌握MATLAB的矩阵和数组的运算3) 掌握MATLAB符号表达式的创建4) 熟悉符号方程的求解实验内容（任选6题）1. 利用rand等函数产生下列矩阵：产生一个均匀分布在（-5，5）之间的随机阵（502），要求显示精度为精确到小数点后一位（精度控制指令为format）。format bankM=10*rand(50,2)-52. 在一个已知的测量矩阵T（100100）中，删除整行数据全为0的行，删除整列数据全为0的列（判断某列元素是否为0方法：检查T(: , i) .* (T(: , i)是否为0）。T = rand(100,100);T(3,:) = 0;T(:,4) = 0;row,col = size(T);r = ; c = ;for i = 1:rowif(sum(T(i,:).*T(i,:)=0) r(size(r)+1) = i;endendif(size(r)T(r,:) = ;endfor j = 1:colif(sum(T(:,j).*T(:,j)=0) c(size(c)+1) = j;endendif(size(c)T(:,c) = end 3. 假设汽车系在下列各年度的人口统计如表所示年份类 别大一新生学士毕业生2002102822003120100200498942005105972006121110试用一个二维矩阵STU表示上述数据，并请计算出下列各数值：（1） 汽车系在20022006年之间的每年平均新生、毕业生人数。（2） 5年来汽车系共有多少毕业生？（3） 在哪几年，新生数目比毕业生多？（4） 5年来每年的毕业生对新生的比例平均值为何？STU = 2002 102 82 2003 120 100 2004 98 94 2005 105 97 2006 121 110;Ave_Stu = round(mean(STU(:,2,3) Sum_gra = sum(STU(:,3)Ave_Stu = 109 97Sum_gra = 483 STU_radio = mean(STU(:,3)./STU(:,2) STU_radio = 0.8859 4. 完成下列矩阵运算：（1） 使用randn产生一个（1010）的矩阵A（2） 计算B=（A+A）/2。请注意，B一定是一个对称矩阵（3） 计算矩阵B的特征向量e1 e2 e10（使用函数eig）（4） 验证在i不等于j的情况下，ei和ej的内积必定为0A=randn(10,10)B=(A+A)/2V,D=eig(B)5. 下列 MATLAB 语句用于画出函数 在 0, 10 的值。 x = 0:0.1:10; y = 2*exp(-0.2*x); plot(x, y); 利用 MATLAB的M-文件编辑器创建一个新的 M-文件，输入上述代码，保存为文件test1.m。然后在命令窗口中输入 test1 执行这个文件。看得到什么结果？ x = 0:0.1:10; y = 2*exp(-0.2*x); plot(x, y); 6. 熟悉MATALB帮助的使用：1) 通过以下两种方式得到关于 exp函数的帮助 （a） 在命令窗口中输入 help exp 命令；（b） 运用帮助空间窗口。2) 使用lookfor 命令查找出以10为底的对数函数，并计算1、10、1000的对数值。help exp EXP Exponential. EXP(X) is the exponential of the elements of X, e to the X. For complex Z=X+i*Y, EXP(Z) = EXP(X)*(COS(Y)+i*SIN(Y). See also EXPM1, LOG, LOG10, EXPM, EXPINT. Overloaded methods: lti/exp codistributed/exp Reference page in Help browser doc exp lookfor log10 log10 - Common (base 10) logarithm. LOG10 Common (base 10) logarithm. LOG10(X) is the base 10 logarithm of the elements of X. Complex results are produced if X is not positive. Class support for input X: float: double, single See also log, log2, exp, logm Overloaded methods: codistributed/log10 Reference page in Help browser doc log10log10(1)log10(10)log10(1000) ans = 0ans = 1ans = 3 7. 假设u=1和v=3，用 MATLAB 符号计算功能计算下列表达式的值：a. b.c.d.syms n v;a = 4*u/(3*v)b = 2*v(-2)/(u+v)c = v3/(v3-u3)d = 4*pi*v2/3 8. a = 3, A = 4, b = a2, B = b2-1, c = a+A-2B, C = a+B+2c，求C。（使用符号函数solve）S = solve(a-3,A-4,b-a2,B-b2+1,c-a-A+2*B,C-a-B-2*c);disp(A,B,C,c,b,a)disp(S.A,S.B,S.C,S.b,S.a)9. 求解常微分方程x=-x+x+1，x(0)=1，x(0)=0（使用符号函数dsolve）S = dsolve(D2x=Dx=Dx+x+1, Dx(0) = 1, x(0) = 0) 10. 应用 MATLAB 帮助窗口查找显示当前目录的命令。a) MATLAB 启动时的当前目录是什么？ D:My DocumentsMATLABb) 创建一个新的目录mynewdir，并将其改变为当前目录。然后打开M-文件编辑窗口，增加以下语句： % create an input array from -2*pi to 2*pit = -2*pi:pi/10:2*pi; % calculate|sin(t)| x=abs(sin(t); %plot result plot(t, x); 把此文件以 test2.m 为文件名保存，然后在命令窗口中输入 test2 执行此文件。 % create an input array from -2*pi to 2*pit = -2*pi:pi/10:2*pi; % calculate|sin(t)| x=abs(sin(t); %plot result plot(t, x); c) 关闭绘图窗口，返回到原始目录，然后在命令窗口中输入 test2。看有何情况发生，为什么。? Undefined function or variable test2.因为Matlab命令窗口要执行的 test2已不存在于当前目录下思考题1. MATLAB中，数组与矩阵在表示与应用上有哪些区别。二维数组相当于矩阵，所以矩阵是数组的子集。数组运算是指数组对应元素之间的运算,也称点运算.矩阵的乘法、乘方和除法有特殊的数学含义，并不是数组对应元素的运算，所以数组乘法、乘方和除法的运算符前特别加了一个点。矩阵是一个二维数组，所以矩阵的加、减、数乘等运算与数组运算是一致的。但有两点要注意：(1)对于乘法、乘方和除法等三种运算，矩阵运算与数组运算的运算符及含义都不同：矩阵运算按线性变换定义，使用通常符号；数组运算按对应元素运算定义，使用点运算符；(2)数与矩阵加减、矩阵除法在数学是没有意义的，在MATLAB中为简便起见，定义了这两类运算 实验04051002 MATLAB科学计算及绘图操作成绩报告成绩实验目的1) 熟悉MATLAB所提供的常用数值计算的函数(方程（组）的求解、插值、拟合)；2) 掌握MATLAB二维图形绘制命令及其图形控制（plot、loglog、contour、polar等）；3) 熟悉MATLAB三维图形绘制命令及其图形控制(mesh、surf等)。实验内容1. 求f(x) = 4x6-x+x3-95在（0，100）上最大值、最小值，根。（使用函数fminbnd、roots）yx1 = (x) 4*x6-x+x3-95;Xmin,ymin=fminbnd(yx1,-10,100)yx2 = (x) (4*x6-x+x3-95)*(-1);Xmax,ymax=fminbnd(yx2,-10,100)Ymax=-ymax;P= 4 0 0 1 0 -1 -95V=roots(p) 2. 请用梯形法、辛普森法分别计算积分值 （trapz、quad）采用trapzformat longd = 0.01; x = 0 : d : 10;y = sqrt (x.2+x+1);Itrapz = d * trapz(x,y) 采用quad fx = sqrt(x.2+x+1);Ic = quad(fx,0,1) %Ic = quad(sqrt(x.2+x+1),0,1)3. ，用两种方法求函数的根，并求其极值与零点。yx=(x)x3+(x-0.8)2/(x+1.25)3-5*(x+1/x); xi1=fzero(yx,100)S=solve(x3+(x-0.8)2/(x+1.25)3-5*(x+1/x),x); xi2=double(S)yx1=(x)(x3+(x-0.8)2/(x+1.25)3-5*(x+1/x)*(-1);x0,y0=fminbnd(yx,0,100);x1,y1=fminbnd(yx1,0,100);x0ymin=y0x1ymax=abs(y1) 4. 计算二重积分 （使用函数dblquad）format long S1=dblquad(inline(x.2 + y.2 + x*y + 2*x + y + 1),0,2,0,1)%采用匿名函数表示被积函数 S1 = 11.333333333333334 5. 求解常微分方程x=-x+x+1，x(0)=1，x(0)=0（使用函数ode45，需建立M函数）%function ydot = DyDt( t,y )ydot = y(2);-y(2)+y(1)+1;tspan=0,30;y0=1;0;tt,yy=ode45(DyDt,tspan,y0);plot(tt,yy(:,1)xlabel(t),title(x(t) 6. 矩阵M = 1, 2, 6; 4, 2, 7; 8, 9, 3，求M的LU分解，QR分解，对角阵，特征值分解。（使用函数lu、qr、svd、eig）M=1,2,6; 4,2,7; 8,9,3 L,U=lu(M) Q,R = qr(M) S,V = svd(M) V,D = eig(M) 7. a = 3, A = 4, b = a2, B = b2-1, c = a+A-2B, C = a+B+2c，求C。S=solve(a-3,A-4,b-a2,B-b2+1,c-a-A+2*B,C-a-B-2*c);C=S.C 8. 要求使用两种方法求解（建立矩阵A、x、b，然后求解方程Ax=b）。A = rand(4,4);b = rand(4,1);x1 = Abx2 = inv(A)*b 9. 建立一个55矩阵，分别用polyval函数和polyvalm函数将矩阵代入函数f(x) = 1.35+0.0668x+0.436x2+0.695x3 ，计算出结果并比较。M = rand (5,5);p = 0.695 0.436 0.0668 1.35; p1 = polyval (p,M) p2 = polyvalm (p,M) 10. 已知t = an2 + bn，测得对应数据如下：（多项式插值interp1，多项式拟合polyfit）t = 0, 20, 40, 60, 80, 100, 120, 140, 160, 183.5;n = 0, 1153, 2045, 2800, 3466, 4068, 4621, 5135, 5619, 6152;通过多项式拟合试求a和b的值。format bankt = 0,20,40,60,80,100,120,140,160,183.5;n = 0,1153,2045,2800,3466,4068,4621,5135,5619,6152;t1=1:0.5:183.5;n1=interp1(t,n,t1,spline);p = polyfit(t, n,2)a=p(1)b=p(2) 11. 北京内燃机总厂生产的492Q发动机，由试验测得的转矩特性为：发动机转速n（rmin-1）10001500200025003000350038004000发动机转矩T（N m）135.33147.10152.98156.91147.10138.27133.37125.53试求出转矩特性（T = f(n)）的多项式表示，并作出发动机的外特性图（要有适当的标注）。format bankn = 1000 1500 2000 2500 3000 3500 3800 4000;T = 135.33 147.10 152.98 156.91 147.10 138.27 133.37 125.53;p = polyfit (n0, te0, 5) ; n = 1000:50:4000te=poly2str (p,n) plot (n,te,Linewidch,3) ;title(北京492Q 发动机外特性);xlabel(发动机转速 (r/min);ylabel(发动机转矩 (N*m);grid onaxis onbox on 12. 假设汽车系在2004、2005、2006年的人员组成如表所示：年份类 别大一新生学士毕业生2002102822003120100200498942005105972006121110a) 请用bar3命令来画出上述数据的立体长条图b) 请用各种你可以想到的方法，在上述图形中加入各种解释和说明的文字。Year=2002:2006;grad1=102,10,98,105,121grad4=82,100,94,97,110;bar3(yea,grad1,grad4)title(学生人数年份)；思考题1. MATLAB求多项式的根是用什么方法，与传统方法相比有何优点MATLAB的多项式求根算法是采用roots(c)指令，行量c为多项式系数。它是借助于所谓的伴随矩阵特征值计算实现的，与传统的方法相比该算法最稳定、可靠。2. 画出横坐标在（-15，15）上的函数的曲线，应该使用什么命令。 ezplot(sin(x),-15,15) 3. 请思考网络线有什么作用，为什么要对图形进行标注。网格线可以增加图形的可读性。标注使图形表达的信息更加清楚，便于读图。实验04051003 MATLAB综合实例编程操作成绩报告成绩实验目的1) 了解Windows界面编程的基本概念和方法掌握MATLAB程序设计的方法；2) 熟悉MATLAB/GUI的基本特点；掌握MATLAB/GUI编制的基本步骤；3) 掌握MATLAB/Simulink的使用方法和基本步骤；4) 将MATLAB应用到所学专业。实验内容说明：本次实验为“MATLAB编程综合实例”，不同专业的同学应结合自己的专业进行实验。例如，车辆工程相关专业的同学可对汽车性能（汽车理论课程所研究的性能）进行仿真，材料工程专业同学可结合工程热力学、传热学等课程的内容参照进行。1. 建立二自由度汽车SIMULINK模型。更改参数，观察不同参数下前轮角阶跃输入下的汽车稳态响应。1） MATLAB指令窗口下输入指令simulink，打开simulink编辑器并新建一个mdl。2） 分别将两个拉入编辑窗口，并命名为k1、k2。3） 将拉入编辑窗口并双击打开属性窗口，将“List of signs”改为“+”，命名为k1+k2，再将k1、k2与之相连。4） 如步骤2）创建B，将拉入编辑窗口，再将k1+k2、B与之相连。5） 分别将、拉入编辑窗口，更改相应属性并连线6） 如步骤2）创建a、b，将拉入2个编辑窗口并双击打开属性窗口，分别将“Expression”改为“u*u”， “u*u”。7） 采用如上方法，得到最终框图。8） 点击运行仿真，结束后双击“Scope”与“Scope1”观察曲线9） 更改k1、k2等参数，运行仿真，不同参数下前轮角阶跃输入下的汽车稳态响应曲线。2. 结合汽车理论知识，针对汽车某项性能设计一仿真软件，要求可根据输入参数绘制相应曲线。（选做一题）选第一题1) 动力性计算，参数见汽车理论习题1.3设计完成以下工作：a） 可进行整车参数输入；b） 计算最高车速、加速时间（2档起步070km/h），最大爬坡度；c） 绘制驱动力阻力平衡图、汽车加速时间曲线、汽车爬坡度图。1）实现按钮功能的程序21function pushbutton1_Callback(hObject, eventdata, handles)global nmin nmax m1 m2 m3 r nt f CdA i0 If Iw1 Iw2 L a hg ig1 ig2 ig3 ig4 ig5;nmin=600;nmax=4000;m1=2000;m2=1800;m3=3880;r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;ig1=5.56;ig2=2.769;ig3=1.644;ig4=1.00;ig5=0.793;set(handles.edit1,string,num2str(nmin);set(handles.edit2,string,num2str(nmax);set(handles.edit3,string,num2str(m1);set(handles.edit4,string,num2str(m2);set(handles.edit5,string,num2str(m3);set(handles.edit6,string,num2str(r);set(handles.edit7,string,num2str(nt);set(handles.edit8,string,num2str(f);set(handles.edit9,string,num2str(CdA);set(handles.edit10,string,num2str(i0);set(handles.edit11,string,num2str(If);set(handles.edit12,string,num2str(Iw1);set(handles.edit13,string,num2str(Iw2);set(handles.edit14,string,num2str(L);set(handles.edit15,string,num2str(a);set(handles.edit16,string,num2str(hg);set(handles.edit17,string,num2str(ig1);set(handles.edit18,string,num2str(ig2);set(handles.edit19,string,num2str(ig3);set(handles.edit20,string,num2str(ig4);set(handles.edit21,string,num2str(ig5);2）动力性计算程序function pushbutton2_Callback(hObject, eventdata, handles)global nmin nmax m1 m2 m3 r nt f CdA i0 If Iw1 Iw2 ig1 ig2 ig3 ig4 ig5 L a hg;nmin=600;nmax=4000;m1=2000;m2=1800;m3=3880;r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;ig1=5.56;ig2=2.769;ig3=1.644;ig4=1.00;ig5=0.793;L=3.2;a=1.947;hg=0.9;n=linspace(600,4000,1000); Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;Ft1=(Tq*ig1*i0*nt)/r;Ft2=(Tq*ig2*i0*nt)/r;Ft3=(Tq*ig3*i0*nt)/r;Ft4=(Tq*ig4*i0*nt)/r;Ft5=(Tq*ig5*i0*nt)/r;Ff=m3*9.8*f;ua=linspace(0,200,1000);ua1=0.377*r*n/ig1/i0;ua2=0.377*r*n/ig2/i0;ua3=0.377*r*n/ig3/i0;ua4=0.377*r*n/ig4/i0;ua5=0.377*r*n/ig5/i0;F1=Ff+CdA*ua1.2/21.15;F2=Ff+CdA*ua2.2/21.15;F3=Ff+CdA*ua3.2/21.15;F4=Ff+CdA*ua4.2/21.15;F5=Ff+CdA*ua5.2/21.15;v=95:0.01:100; FD5=interp1(ua5,Ft5,v,spline);FF5=interp1(ua5,F5,v,spline);for i=1:501; if FD5(i)FF5(i); uamax=v(i); break endendset(handles.edit1,string,num2str(uamax); %输出最高车速的结果m1=2000; %开始计算最大加速时间m2=1800;m3=3880;r=0.367;nt=0.85;f=0.013;CdA=2.77;io=5.83;If=0.218;Iw1=1.798;Iw2=3.598;g=9.8;L=3.2;a=1.947;hg=0.9;ig=5.56 2.769 1.644 1.00 0.793;n=600:1:4000;Tq=-19.313+295.27*n/1000-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;Ft1=Tq*ig(1)*io*nt/r;Ft2=Tq*ig(2)*io*nt/r;Ft3=Tq*ig(3)*io*nt/r;Ft4=Tq*ig(4)*io*nt/r;Ft5=Tq*ig(5)*io*nt/r;ua1=0.377*r*n/(io*ig(1);ua2=0.377*r*n/(io*ig(2);ua3=0.377*r*n/(io*ig(3);ua4=0.377*r*n/(io*ig(4);ua5=0.377*r*n/(io*ig(5);ua=0:1:130;Fz=m3*9.8*f+CdA*ua.2/21.15;Pe1=Ft1.*ua1./(3600*nt);Pe2=Ft2.*ua2./(3600*nt);Pe3=Ft3.*ua3./(3600*nt);Pe4=Ft4.*ua4./(3600*nt);Pe5=Ft5.*ua5./(3600*nt);Pe=Tq.*n/9550;Pem=max(Pe);Pz=Fz.*ua./(3600*nt);d1=1+(Iw1+Iw2)/(m3*r2)+If*ig(1)*r2*io2*nt/(m3*r2);d2=1+(Iw1+Iw2)/(m3*r2)+If*ig(2)*r2*io2*nt/(m3*r2);d3=1+(Iw1+Iw2)/(m3*r2)+If*ig(3)*r2*io2*nt/(m3*r2);d4=1+(Iw1+Iw2)/(m3*r2)+If*ig(4)*r2*io2*nt/(m3*r2);d5=1+(Iw1+Iw2)/(m3*r2)+If*ig(5)*r2*io2*nt/(m3*r2);a1=(Ft1-m3*g*f-CdA*ua1.2/21.15)/(d1+m3);a2=(Ft2-m3*g*f-CdA*ua1.2/21.15)/(d2+m3);a3=(Ft3-m3*g*f-CdA*ua3.2/21.15)/(d3+m3);a4=(Ft4-m3*g*f-CdA*ua4.2/21.15)/(d4+m3);a5=(Ft5-m3*g*f-CdA*ua5.2/21.15)/(d5+m3);aj1=1./a1;aj2=1./a2;aj3=1./a3;aj4=1./a4;aj5=1./a5;ua1m=max(ua1);ua2m=max(ua2);ua3m=max(ua3);ua4m=max(ua4);ua5m=max(ua5);UA=ua1;ua2;ua3;ua4;ua5;AD=aj1;aj2;aj3;aj4;aj5;t(1)=0; ii=2; for uu=10:1:70 if uuua1m&uuua2m&uuua3m&uuua4m&uuua5m; aci=interp1(UA(5,:),AD(5,:),uu,spline); t(ii)=t(ii-1)+aci*1; end ii=ii+1; end tmax=t(end); %act(1) u=linspace(10,70,62); y=40+max(a5)*t; set(handles.edit3,string,num2str(tmax); %输出最大加速时间结果 nmin=600; %开始计算最大爬坡度nmax=4000;m1=2000;m2=1800;m3=3880;r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;ig1=5.56;ig2=2.769;ig3=1.644;ig4=1.00;ig5=0.793;L=3.2;a=1.947;hg=0.9;n=linspace(nmin,nmax,1000);Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;Ft1=(Tq*ig1*i0*nt)/r;Ft2=(Tq*ig2*i0*nt)/r;Ft3=(Tq*ig3*i0*nt)/r;Ft4=(Tq*ig4*i0*nt)/r;Ft5=(Tq*ig5*i0*nt)/r;ua=linspace(0,200,1000);ua1=0.377*r*n/ig1/i0;ua2=0.377*r*n/ig2/i0;ua3=0.377*r*n/ig3/i0;ua4=0.377*r*n/ig4/i0;ua5=0.377*r*n/ig5/i0;Ff=m3*9.8*f;Fw=CdA*ua.2/21.15;Fw1=CdA*ua1.2/21.15;Fw2=CdA*ua2.2/21.15;Fw3=CdA*ua3.2/21.15;Fw4=CdA*ua4.2/21.15;Fw5=CdA*ua5.2/21.15;Fi1=Ft1-Ff-Fw1;Fi2=Ft2-Ff-Fw2;Fi3=Ft3-Ff-Fw3;Fi4=Ft4-Ff-Fw4;Fi5=Ft5-Ff-Fw5;i1=tan(asin(Fi1./(m3*9.8)*100%;i2=tan(asin(Fi2./(m3*9.8)*100%;i3=tan(asin(Fi3./(m3*9.8)*100%;i4=tan(asin(Fi4./(m3*9.8)*100%;i5=tan(asin(Fi5./(m3*9.8)*100%;imax=max(i1);set(handles.edit2,string,num2str(imax); %输出最大爬坡度结果3）绘制汽车驱动力阻力平衡图程序 nmin=600;nmax=4000;m1=2000;m2=1800;m3=3880;r=0.367;nt=0.85;f=0.013;CdA=2.77;i0=5.83;If=0.218;Iw1=1.798;Iw2=3.598;ig1=5.56;ig2=2.769;ig3=1.644;ig4=1.00;ig5=0.793;L=3.2;a=1.947;hg=0.9;n=linspace(nmin,nmax,1000);Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;Ft1=(Tq*ig1*i0*nt)/r;Ft2=(Tq*ig2*i0*nt)/r;Ft3=(Tq*ig3*i0*nt)/r;Ft4=(Tq*ig4*i0*nt)/r;Ft5=(Tq*ig5*i0*nt)/r;ua=linspace(0,200,1000);ua1=0.377*r*n/ig1/i0;ua2=0.377*r*n/ig2/i0;ua3=0.377*r*n/ig3/i0;ua4=0.377*r*n/ig4/i0;ua5=0.377*r*n/ig5/i0;Ff=m3*9.8*f;Fw=CdA*ua.2/21.15;Fz=Fw+Ff;h=plot(ua1,Ft1,y,ua2,Ft2,b,ua3,Ft3,r,ua4,Ft4,c,ua5,Ft5,k,ua,Fz,g);title(汽车驱动力-阻力平衡图);xlabel(速度ua/(km/h);ylabel(力F/N );legend(h,Ft1,Ft2,Ft3,Ft4,Ft5,Ff+Fw);4）汽车加速时间曲线图程序m1=2000;m2=1800;m3=3880;r=0.367;nt=0.85;f=0.013;CdA=2.77;io=5.83;If=0.218;Iw1=1.798;Iw2=3.598;g=9.8;L=3.2;a=1.947;hg=0.9;ig=5.56 2.769 1.644 1.00 0.793;n=600:1:4000;Tq=-19.313+295.27*n/1000-165.44*(n/1000).2+40.874*(n/1000).3-3.8445*(n/1000).4;Ft1=Tq*ig(1)*io*nt/r;Ft2=Tq*ig(2)*io*nt/r;Ft3=Tq*ig(3)*io*nt/r;Ft4=Tq*ig(4)*io*nt/r;Ft5=Tq*ig(5)*io*nt/r;ua1=0.377*r*n/(io*ig(1);ua2=0.377*r*n/(io*ig(2);ua3=0.377*r*n/(io*ig(3);ua4=0.377*r*n/(io*ig(4);ua5=0.377*r*n/(io*ig(5);ua=0:1:130;Fz=m3*9.8*f+CdA*ua.2/21.15;Pe1=Ft1.*ua1./(3600*nt);Pe2=Ft2.*ua2./(3600*nt);Pe3=Ft3.*ua3./(3600*nt);Pe4=Ft4.*ua4./(3600*nt);Pe5=Ft5.*ua5./(3600*nt);Pe=Tq.*n/9550;Pem=max(Pe);Pz=Fz.*ua./(3600*nt);d1=1+(Iw1+Iw2)/(m3*r2)+If*ig(1)*r2*io2*nt/(m3*r2);d2=1+(Iw1+Iw2)/(m3*r2)+If*ig(2)*r2*io2*nt/(m3*r2);d3=1+(Iw1+Iw2)/(m3*r2)+If*ig(3)*r2*io2*nt/(m3*r2);d4=1+(Iw1+Iw2)/(m3*r2)+If*ig(4)*r2*io2*nt/(m3*r2);d5=1+(Iw1+Iw2)/(m3*r2)+If*ig(5)*r2*io2*nt/(m3*r2);a1=(Ft1-m3*g*f-CdA*ua1.2/21.15)/(d1+m3);a2=(Ft2-m3*g*f-CdA*ua1.2/21.15)/(d2+m3);a3=(Ft3-m3*g*f-CdA*ua3.2/21.15)/(d3+m3);a4=(Ft4-m3*g*f-CdA*ua4.2/21.15)/(d4+m3);a5=(Ft5-m3*g*f-CdA*ua5.2/21.15)/(d5+m3);aj1=1./a1;aj2=1./a2;aj3=1./a3;aj4=1./a4;aj5=1./a5;ua1m=max(ua1);ua2m=max(ua2);ua3m=max(ua3);ua4m=max(ua4);ua5m=max(ua5);UA=ua1;ua2;ua3;ua4;ua5;AD=aj1;aj2;aj3;aj4;aj5;t(1)=0; ii=2; for uu=10:1:70 if uuua1m&uuua2m&uuua3m&uuua4m&uuua5m; aci=interp1(UA(5,:),AD(5,:),uu,spline); t(ii)=t(ii-1)+aci*1; end ii=ii+1; end tmax=t(end); u=linspace(10,70,62); y=40+max(a5)*t; h=plot(t,u,g,t,y,m,linewidth,2.5); title(汽车加速时间曲线); xlabel(时间t/s); ylabel(速度ua/(km/h); legend(h,挡原地起步换挡加速,直接挡加速) grid on;5）汽车爬坡度图程序nmin=600;nmax=4000;m1=2000;m2=1800;m3=3880;r=0.367;nt=0.85;f=0