模拟试题-MATLAB

模拟试题-MATLAB一、单项选择题1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用(c)(A)diag(magic(3));(B)diag(magic);(C)diag(diag(magic(3)));(D)diag(diag(magic))。2、MATLAB命令P=pascal(3)将创建三阶帕斯卡矩阵，max(P)的计算结果是(d)(A)123(B)121(C)3610(D)1363、命令J=[1;1;1]*[1,2,3];A=j+j’-I将创建矩阵(a)(A)123234345;(B)234345456(C)123123123(D)1112223334、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(ysqrt(x)&yx.^2);的功能是(b)(A)统计2000个随机点中落入特殊区域的点的索引值；(B)统计1000个随机点落入特殊区域的点的索引值；(C)模拟2000个随机点落入特殊区域的过程；(D)模拟1000个随机点落入特殊区域的过程。5、MATLAB计算二项分布随机变量分布律的方法是(c)(A)binocdf(x,n,p);(B)normpdf(x,mu,s);(C)binopdf(x,n,p);(D)binornd(x,n,p)。6、MATLAB命令symse2;f=sqrt(1-e2*cos(t)^2);S=int(f,t,0,pi/2)功能是(d)(A)计算f(x)d[0,2/]上的积分；(B)计算f(t)的不定积分符号结果；(C)计算f(x)的积分的数年结果(double(S))；(D)计算f(t)定积分的符号结果。7、y=dsolve(‘Dy=1/(1+x^2)-2*y^2’,’y(0)=0’,’x’);ezplot(y)的功能是(a)(A)求微分言和特解并绘图；(B)解代数方程(C)求定积分；(D)求微分方程通解。8、X=10000;0.5*asin(9.8*X/(515^2))的功能是计算关于抛射体问题的(a)(实验一)(A)十公里发射角；(B)十公里飞行时间；(C)最大飞行时间；(D)最大射程。9．MATLAB命令A=rand(5,5)；创建55)(ijaA，求51||maxiijja用(a)(A)max(sum(abs(A)))（列和）;(B)max(sum(abs(A’)))(行和);(C)max(sum(A)));(D)sum(max(A));10．MATLAB命令x=[1,2,4,5,9];mean(x)(求平均值)，的计算结果是(b)(A)4(B)4.2(B)4.5(D)2111．MATLAB命令x=rand(10,1)生成10个随机数，将它们从大到小排序，使用(c)(A)y=sort(x);z=y(10:1);(B)[y,II]=sort(x);z=y(II);(C)y=sort(x);z=y(10:-1:1);(D)[y,II]=sort(x);z=x(II);12．MATLAB命令roots([1,0,0,-1])的功能是(d)(A)产生向量[1,0,0,1]；(B)求方程013x的根；(C)求多项式13x的值(D)求方程013x的根。13．MATLAB命令A=magic(3)创建3阶幻方矩阵，求A的特征值绝对值最小用(a)(A)min(abs(eig(A)));(B)min(eig(abs(A)));(C)min(eig(A));(D)min(abs(A));14．命令factor()用于分解因式，symsx;f=4*x^3+9*x^2-30*x;factor(diff(f))的结果是(b)(A)(x-1)*(2*x-5)(B)6*(x-1)*(2*x+5)(C)6*(x+1)*(2*x+5)(D)(x+1)*(2*x-5)15．MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,0,pi)功能是(c)(A)绘出函数f在[0,2]图形；(B)计算函数f在[0,2]的积分；(C)计算旋转曲面所围的体积；(D)计算旋转曲面的表面积。16．十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令k=mod(2008,12)+1的结果是(d)(A)k指向第二动物牛；(B)k指向第三动物虎；(C)k指向第四动物兔；(D)k指向第五动物龙。17．MATLAB命令[x,y]=meshgrid(1:3);H=1./(x+y-I)(逆矩阵)产生的矩阵H是(D)(A)333222111(B)321321321(C)543332321(D)5/14/13/14/13/12/13/12/1118．下面有关MATLAB变量名和函数名的说法中，错误的说法是(d)(A)变量名的第一个字符必须是一个英文字母(B)变量名可由英文字母、数字和下划线混合组成(C)变量名不得包含空格和标点，但可以有下连字符(D)变量名和函数名对于英文的大小使用没有区别19、北京和纽约的经度分别是：东经118和西经76，根据经度差计算时差用(d)(A)fai1=118;fai2=-76;Dfai=(fai1+fai2)/24;(B)fai1=118;fai2=-76;Dfai=(fai1+fai2)/15;(C)fai1=118;fai2=-76;Dfai=(fai1-fai2)/24;(D)fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。二、程序阅读题(40分)1．传说古希腊曾流行瘟疫，人们为消除灾难求助于神。神说：把神庙中黄金祭台增容一倍，可消除瘟疫。当立方体祭台尺寸放大一倍后，瘟疫仍然流行。人们才知道体积并不是扩大了两倍。这个古希腊难题被称为倍立方体问题，在人类还没有认识到无理数时，企业界企图用有限位实数表示32，就会犯下错误。数学实验程序验证了这个事实，程序运行后误差如右文本框所示a=2^(1/3);D=1;fork=1:8D=D*10;b=fix(a*D)/D;V(k)=b^3;enderror=V’-2(1)程序中循环控制变量k从1变量8，而变量D=10k的作用是(c)(A)将a的小数点向右移D位取整；(B)将a的小数点向右移D位取整后再向左移D位；(C)将a的小数点向右移k位取整后再向左移k位；(D)将a的小数点向左移k位取整后再向右移k位；(2)程序中变量b存放的数据是(d)(A)将a的小数点后第k位减1所得；(B)将a的小数点k位后按四舍五入所得；(C)将a的小数点后第k位增1所得；(D)将a的小数点k位后截断舍去所得。error=-2.7200e-001-4.6875e-002-4.3830e-003-1.0024e-004-4.9998e-006-2.3761e-007-2.3761e-007-4.7121e-0082、3n+1问题反映一个数学猜想：对任一自然数n，按如下法则进行运算：若n为偶数，则将n除2，若n为奇数，则将n乘3加1。重复这种操作，结果终会为1。实验程序如下。function[k,N]=threeN(n)ifnargin==0,n=5;endk=1;N=n;whilen~=1r=rem(n,2);ifr==0n=n/2;elsen=3*n+1;endN=[N,n];k=k+1;end(1)在MATLAB命令窗口中直接调用threeN运行结果为(A)(A)只显示k的最后数值为6；(B)只显示k的最后数值5；(C)同时显示k和N的数据；(D)仅显示N的所有数据。(2)实验程序运行过程中(B)(A)输入变量n不发生改变；(B)N是记录数据变化的一维数组；(C)N记录每次数据变化的单个数据；(D)n是记录数据变化的一维数组。3、关于“牟合方盖”的实验程序如下h=2*pi/100;t=0:h:2*pi;r=0:0.05:1;x=r'*cos(t);y=r'*sin(t);z=sqrt(1-x.^2);%第三行meshz(x,y,z),axisoffcolormap([001])view(-47,56),holdonx1=cos(t);y1=sin(t);z1=abs(sin(t));plot3(x1,y1,z1,'ro')；（1）下面有关程序的功能的说法确切的是（d）（A）绘圆柱面x2+y2=1,x2+z2=1的交线；（B）绘圆柱面x2+y2=1,x2+z2=1所围区域的边界曲面；（C）绘圆柱面x2+y2=1,x2+z2=1的交线及所围区域的边界曲面；（D）绘圆柱面x2+y2=1,x2+z2=1的交线及所围区域的边界曲面的上半部分。（2）关于第三行语句错误的解释是（a）（A）z是矩形域上曲顶柱面高度值；（B）z是与y同型的矩阵；（C）z是圆域上曲顶柱面高度值；（D）z是与x同型的矩阵4、非负函数y=f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴的旋转曲面，其体积symsabxf=exp(a*x)*sin(b*x);f1=subs(f,a,-0.2);f2=subs(f1,b,0.5);V=pi*int(f2*f2,x,0,2*pi)double(V)220[()]Vfxdxt=(0:20)*pi/10;theta=t;r=f2(t);x=t'*ones(size(t));%第九行y=r'*cos(theta);%第十行z=r'*sin(theta);%第十一行mesh(z,y,x)colormap([000])axisoffview(-17,54)（1）关于程序的功能确切的说法（d）(A)计算曲线段f(x)=exp(ax)sin(bx),绕X轴旋转的旋转曲面体积(B)计算曲线段f(x)=exp(-0.2x)sin(0.5x),绕X轴旋转的旋转曲面体积(C)计算曲线段f(x)=exp(ax)sin(bx),绕X轴旋转的旋转曲面体积并绘图（D）计算曲线段f(x)=exp(-0.2x)sin(0.5x),绕X轴旋转的旋转曲面体积并绘图（2）由第九行至第十一行语句可得旋转曲面的方程（b）(A)（B）（C）(D)三、程序填空(40分)1、中国农历60年一大轮回，按天干“甲乙丙丁戊已庚辛壬癸”和地支“子丑寅卯辰巳午未申酉戍亥”循环排列而成。已知2009年是农历已丑年，通过简单计算可以找出年份与天干/地支对应的规律。下面数学实验程序对输入年份，计算并输出字符串农历纪年。填空完善程序。functioncalendar=year(year)ifnargin==0,year=2009;endS1=’甲乙丙丁戊已庚辛壬癸’;S2=’子丑寅卯辰巳午未申酉戍亥’;k1=mod(year-4,10)+1①;%定位天干序数s1=S1(k1);k2=mod(year-4,12)+1②;%定位地支序数s2=S2(k2);calendar=strcat(int2str(year),’年是’,s1,s2,’年’)2、红、绿两队从相距100公里的地点同时出发相向行军。红队速度为10(公里/小时)，绿队速度为8(公里/小时)。开始时，通讯员骑摩托从红队出发为行进中的两队传递消息。摩托车的速度为60(公里/小时)往返于两队之间。每遇一队，立即回驶向另一队。当两队距离小于0.2公里时，摩托车停止，下面数学实验程序模拟计算摩托车跑了多少趟。请填空完善程序。functionk=moto(A,B)ifnargin==0,A=0;B=100;endva=10;vb=8;vc=60;f=1;k=0;02x02x02x02x222()fxyzt22()fxyz22()yfxz222()yfxztwhile(B-A)0.2iff==1tk=(B-A)/(vb+vc);elsetk=(B-A)/(va+vc)①;%计算A与C相遇时间endA=A+va*tk②;%计算A点位置