电子科大成都学院数学实验复习题

2011年数学实验复习题一、选择题1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用(C)(A)diag(magic(3));(B)diag(magic);(C)diag(diag(magic(3)));(D)diag(diag(magic))。2、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(y<sqrt(x)&y>x.^2);的功能是(B)(A)统计2000个随机点中落入特殊区域的点的索引值；(B)统计1000个随机点落入特殊区域的点的索引值；(C)模拟2000个随机点落入特殊区域的过程；(D)模拟1000个随机点落入特殊区域的过程。3、MATLAB计算二项分布随机变量分布律的方法是(C)(A)binocdf(x,n,p);(B)normpdf(x,mu,s);(C)binopdf(x,n,p);(D)binornd(x,n,p)。4、MATLAB命令symse2;f=sqrt(1-e2*cos(t)^2);S=int(f,t,0,pi/2)功能是(D)(A)计算f(x)在[0,pi/2]上的积分；(B)计算f(t)不定积分符号结果；(C)计算f(x)积分的数值结果；5、y=dsolve(‘Dy=1/(1+x^2)2*y^2’,’y(0)=0’,’x’);ezplot(y)的功能是(A)(A)求微分方程特解并绘图；(B)解代数方程(C)求定积分；(D)求微分方程通解。(D)计算f(t)定积分的符号结果。6、X=10000?;0.5*asin(9.8*X/(515^2))的功能是计算关于抛射体问题的(A)(A)十公里发射角；(B)十公里飞行时间；(C)最大飞行时间；(D)最大射程。7、theta=linspace(0,2*pi,100)?;r=cos(4*theta)?;polar(theta,r,’k’)功能是(D)(A)绘四叶玫瑰线；(B)绘三叶玫瑰线；(C)绘心脏线；(D)绘八叶玫瑰线。A=a)5max|a|å´ij558、MATLAB命令A=rand(5,5)；创建，求用(A)ji=1(A)max(sum(abs(A)));(B)max(sum(abs(A’)));(C)max(sum(A)));(D)sum(max(A));9、MATLAB命令x=[1,2,4,5,9];mean(x)，的计算结果是(B)(A)4(B)4.2(B)4.5(D)2110、MATLAB命令x=rand(10,1)生成10个随机数，将它们从大到小排序，使用(C)(A)y=sort(x);z=y(10:1);(B)[y,II]=sort(x);z=y(II);(D)[y,II]=sort(x);z=x(II);(C)y=sort(x);z=y(10:-1;1);11、MATLAB命令roots([1,0,0,-1])的功能是(D)x+1=03(A)产生向量[1,0,0,1]；(B)求方程的根；x-13x-1=03(C)求多项式的值(D)求方程的根。12、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));13、命令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)14、MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,0,2*pi)功能是(C)p(A)绘出函数f在[0,2]图形；p(B)计算函数f在[0,2]的积分；(C)计算旋转曲面所围的体积；(D)计算旋转曲面的表面积。15、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令k=mod(2008,12)+1的结果是(D)(A)k指向第二动物牛；(C)k指向第四动物兔；(B)k指向第三动物虎；(D)k指向第五动物龙。16．下面有关MATLAB变量名和函数名的说法中，错误的说法是(D(A)变量名的第一个字符必须是一个英文字母)(B)变量名可由英文字母、数字和下划线混合组成(C)变量名不得包含空格和标点，但可以有下连字符(D)变量名和函数名对于英文的大小使用没有区别17、在MATLAB命令窗口中，键入命令symsx；int(x*sin(x))。结果是（A）（A）ans=sin(x)-x*cos(x)；（B）ans=cos(x)+x*sin(x)；（C）ans=sin(x)-cos(x)；（D）ans=-1/2*cos(x)*sin(x)+1/2*x18、在MATLAB命令窗口中键入命令A=[142;312;615];det(A(1:2,2:3).*A(2:3,2:3))。结果是（C）（A）ans=-143（B）ans=60（C）ans=36（D）ans=-1919、MATLAB命令x=3:2:100将创建等差数列，该数列是（C）（A）以3为初值的98个数,；（B）以100为终值的98的个数；（C）以99为终值的49个数；（D）以3为初值的97个数。20、在MATLAB命令窗口中输入命令data=[4123131424];y=hist(data,4)，结果是（C）（A）y=4123；（B）y=3232；（C）y=3223；（D）y=421121、MATLAB语句[x,y]=meshgrid(-2:2)的数据结果中（D）（A）x是行向量，y是列向量；（B）x是五行五列的矩阵；（C）x是行元素相同的矩阵；（D）x是列向量相同的矩阵22、MATLAB的语句colormap(001)（D）（A）将三维网面图确定为红色；（B）将三维网面图确定为绿色；（C）将三维网面图确定为蓝色；（D）语句使用格式错误23、火炮发射炮弹的初始速度和发射角为已知，由此可估算出炮弹在空中的飞行时间Tfly，使用语句Tspan=Tfly*(0:20)/20，将获得一些数据，下面不正确的说法是（D）A）Tspan为包括发射时刻在内的炮弹在空间飞行的21个不同的飞行时刻；B）Tspan中任意两个相邻数据之差的绝对值相等；C）Tspan包含了21个数据，第一个数据为0，最后一个数据为Tfly；D）Tspan是一个等差数列，公差为Tfly/2124、北京和纽约的经度分别是：东经118和西经76，根据经度差计算时差用(D)(A)fai1=118;fai2=-76;Dfai=(fai1+fai2)/24;(C)fai1=118;fai2=-76;Dfai=(fai1-fai2)/24;(B)fai1=118;fai2=-76;Dfai=(fai1+fai2)/15;(D)fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。25、用MATLAB随机产生60个1到365之间的正整数，应该使用下面的哪一条命令(B)A）fix(365*rand(1,60))；C）1+fix(364*rand(1,60))；B）1+fix(366*rand(1,60))；D）1+fix(365*rand(1,60))二、程序阅读1、关于“牟合方盖”的实验程序如下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）下面有关程序的功能的说法确切的是（）2222（A）绘圆柱面x+y=1,x+z=1的交线；2222（B）绘圆柱面x+y=1,x+z=1所围区域的边界曲面；2222（C）绘圆柱面x+y=1,x+z=1的交线及所围区域的边界曲面；2222（D）绘圆柱面x+y=1,x+z=1的交线及所围区域的边界曲面的上半部分。（2）关于第三行语句错误的解释是（（A）z是矩形域上曲顶柱面高度值；（B）z是与y同型的矩阵；（C）z是圆域上曲顶柱面高度值；（D）z是与x同型的矩阵）2、非负函数y=f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴的旋转曲面，ò其体积V=p2p[f(x)]dx.20symsabxf=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)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）关于程序的功能确切的说法（(A)计算曲线段f(x)=exp(ax)sin(bx),）绕X轴旋转的旋转曲面体积0xp(B)计算曲线段f(x)=exp(-0.2x)sin(0.5x),(C)计算曲线段f(x)=exp(ax)sin(bx),（D）计算曲线段f(x)=exp(-0.2x)sin(0.5x),绕X轴旋转的旋转曲面体积0xp绕X轴旋转的旋转曲面体积并绘图0xp绕X轴旋转的旋转曲面体积并绘图0xp（2）由第九行至第十一行语句可得旋转曲面的方程（）(A)（C）（B）f(x)=y+z+t222(D)y=f(x+z)2y=f(x+z+t)22223．数学实验程序如下h=439;H=2384;R=6400;a=(h+H+2*R)/2;c=(H-h)/2;e1=c/a;b=sqrt(a*a-c*c);symse2tf=sqrt(1-e2*cos(t)^2);ft=subs(f,e2,e1*e1);S=int(ft,0,pi/2);L=4*a*double(S)；V=L/(114*3600)；%第九行s1=pi*a*b/(114*60);Vmax=2*s1/(h+R)Vmin=2*s1/(H+R)（1）实验程序的运行后，将显示的数据是（）（A）卫星轨道的周长数据；（B）卫星运行的近地速度和远地速度；（C）卫星运行时向径每秒扫过的面积；（D）卫星运行的平均速度数据（2）第九行语句的功能是（）（A）计算卫星运行的最小速度；（B）计算卫星运行时向径每秒扫过的面积;（C）计算卫星运行的最大速度；（D）计算卫星运行轨道的平均速度x+y+z=R(x-R/+y?R/422222224、Viviani体是圆柱体被球面所割立体。下面的数学实验程序功能是取R=2求体积上半部分，先利用符号处理重积分并转换为数值数据，再用蒙特卡罗方法计算体做对比。数学实验程序如下：symsxy;f=sqrt(4-x^2-y^2);y1=-sqrt(2*x-x^2);y2=sqrt(2*x-x^2);S1=int(f,y,y1,y2);S2=int(S1,x,0,2)V=2*double(S2)P=rand(10000,3);X=2*P(:,1);Y=2*P(:,2);Z=2*P(:,3);II=find((X-1).^2+Y.^2<=1&Z<=sqrt(4-X.^2-Y.^2));图1V1=8*length(II)/10000(1)符号计算所用的积分公式是(A)V=dx4-x-ydy2+2x-x2蝌22(A)0xx-2-2(B)=4-x-y蝌2+2x-x2-2x-x2V220V=4-x-ydyV=xydx4--ò+2-2ò+2-2xx2222(C)(D)-2x-x2xx-2-2(2)蒙特卡罗方法选用的随机点变化范围的立方体区域是(B)x,y,z)|xyz?(A)(B)(C)(D)；x,y,z)|xy(1,1),z?x,y,z)|0<x<<y<1,0<z<x,y,z)|0<x<<y<<z<5、某厂生产两种产品，产一吨甲产品用A资源3B资源4m；产一吨乙产品用A资源2吨，B资源6mC资源7个单33位。一吨甲产品和乙产品分别价值7万元和5万元，三种资源限制分别为90200m和210个单位。生产两种产品使总3价值最高的生产方案可用数学实验程序计算。C=[-7,-5];A=[32;46;07];b=[90;200;210];Aeq=[];Beq=[];e0=[0,0];e1=[inf,inf];[x,fval]=linprog(C,A,b,Aeq,beq,e0,e1);(1)程序中变量C表示()(A)目标函数系数；(B)等式约束系数；(C)不等式约束系数；(D)等式约束常向量(2)程序中变量A表示()(A)等式约束矩阵；(B)不等式约束矩阵；(C)决策变量的值；(D)目标函数的最大值三、程序填空1、中国农历60年一大轮回，按天干“甲乙丙丁戊已庚辛壬癸”和地支“子丑寅卯辰巳午未申酉戍亥”循环排列而成。已知2009年是农历已丑年，通过简单计算可以找出年份与天干/地支对应的规律。下面数学实验程序对输入年份，计算并输出字符串农历纪年。填空完善程序。functioncalendar=year(year)ifnargin==0,year=2009;endS1=’甲乙丙丁戊已庚辛壬癸’;S2=’子丑寅卯辰巳午未申酉戍亥’;k1=mod(S1-4,10)+1s1=S1(k1);①;②;%定位天干序数%2010庚%定位地支序数k2=mod(S2-4,12)+1s2=S2(k2);calendar=strcat(int2str(year),’年是’,s1,s2,’年’)2．对于任意正整数n，如果n只能被1和它自身整除，则称这个数为素数(或质数)。判素数程序的算法思想是试商法，即用2，3,……，（n-1）去除n，如果能被这些数中一个整除，则n不是素数，否则是素数。完成下面填空。n=input('inputn:=');fork=2:nifmod(n,k)==0,break,end①endifk<ndisp('不是素数')elsedisp('是素数')end3．汽车紧急刹车问题数据拟合实验V2025303540455055606570T202841537293118149182221266V表示刹车时汽车行驶速度(英里/小时),T表示刹车后汽车滑行距离(英尺)分别做二次多项式和三次多项式数据拟合实验，并绘出数据拟合曲线的图，计算出残差平方和，完成如下实验程序填空v=[2025303540455055606570]*1.609;T=[202841537293118149182221266]*.3048;P2=①T2=polyval(P2,v);R2=sum((T-T2).^2)figure(2),plot(v,T,'*',v,T2)P3=polyfit(v,T,3);T3=polyval(P3,v);R3=②figure(3),plot(v,T,'*',v,T3)4、红、绿两队从相距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;while(B-A)>0.2iff==1tk=(B-A)/(vb+vc);elsetk=(B-A)/(vc+va);①;%计算A与C相遇时间endA=A+va*tkB=B-vb*tkf=-f;②;③;%计算A点位置%计算B点位置k=k+1;end5、二阶正交矩阵作用于某一向量时，其效果是将该向量旋转，旋转解为/50(逆时针旋转为正)。把一个以原点为中心的正三角形旋转，并缩小90%，迭代33次创建图3。完成程序填空：bata=[1/2;7/6;11/6;15/6]*pi;x=cos(bata);y=sin(bata);1line(x,y)xy=[x,y];alfa=pi/50;A=[cos(alfa),-sin(alfa);sin(alfa),cos(alfa)];fork=1:33xy=00.9*xy*A'xy(:,2);x=xy(:,1);y=line(x,y)②;图301end6、长征三号甲运载火箭提供给探月卫星的初始速度不足以将卫星送往月球轨道。为提高到奔月速度，中国航天工程师使最大速度。填空完善下面实验程序。R=6378;Time=[16,15.63,23.3,50.5,225]*3600;h=[200,600,600,600,600];H=[51000,51000,71000,128000,370000];a=(h+H+2*R)/2;c=(H-h)/2;①b=sqrt(a.*a-c.*c);S=pi*a.*b./Time;②Vmax=2*S./(R+h)8、为了进入地月转移轨道，嫦娥一号卫星进行了四次变轨调速度。第一次变轨从16小时初始轨道进入16小时轨道，第二次卫星进入24小时轨道，第三次卫星进入48小时轨