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

1、电子科大成都学院数学实验复习题 Company number【1089WT-1898YT-1W8CB-9UUT-92108】电子科大成都学院数学实验复习题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(yx.2);的功能是( B ) (A) 统计2000个随机点中落入特殊区

2、域的点的索引值； (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命令syms e2;f=sqrt(1-e2*cos(t)2);S=int(f,t,0,pi/2)功能是(D) (A) 计算f(x)在0,pi/2上的积分； (B) 计算f(t)不定积分符号结果

3、； (C) 计算f(x)积分的数值结果； (D) 计算f(t)定积分的符号结果。5、y=dsolve(Dy=1/(1+x2)-2*y2,y(0)=0,x);ezplot(y)的功能是( A ) (A) 求微分方程特解并绘图； (B) 解代数方程 (C) 求定积分； (D)求微分方程通解。6、X=10000;*asin*X/(5152)的功能是计算关于抛射体问题的(A) (A) 十公里发射角； (B) 十公里飞行时间； (C)最大飞行时间； (D)最大射程。7、theta=linspace(0,2*pi,100);r=cos(4*theta);polar(theta,r,k)功能是(D) (A)

4、 绘四叶玫瑰线； (B)绘三叶玫瑰线； (C)绘心脏线； (D) 绘八叶玫瑰线。8、MATLAB命令A=rand(5,5)；创建，求用( A )(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) (B) (D) 2110、MATLAB命令x=rand(10,1)生成10个随机数，将它们从大到小排序，使用( C )(A) y=sort(x);z=y(10:1); (B) y,II=sort(x);z=y(I

5、I);(C) y=sort(x);z=y(10:-1;1); (D) y,II=sort(x);z=x(II);11、MATLAB命令roots(1,0,0,-1)的功能是( D )(A) 产生向量1,0,0,1； (B) 求方程的根；(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()用于分解因式，syms x; f=4*x3+9*x2-30*x

6、; 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命令syms x; f=sin(x); V=pi*int(f*f,x,0,2*pi)功能是( C )(A) 绘出函数f在0,2图形； (B) 计算函数f在0,2的积分；(C) 计算旋转曲面所围的体积； (D) 计算旋转曲面的表面积。15、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令k=mod(2008,12)+1的结果是( D )(A) k指向第二动物牛； (B) k指向第三动物

7、虎；(C) k指向第四动物兔； (D) k指向第五动物龙。16下面有关MATLAB变量名和函数名的说法中，错误的说法是( D )(A) 变量名的第一个字符必须是一个英文字母(B) 变量名可由英文字母、数字和下划线混合组成(C) 变量名不得包含空格和标点，但可以有下连字符(D) 变量名和函数名对于英文的大小使用没有区别17、在MATLAB命令窗口中，键入命令syms x； 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)

8、*sin(x)+1/2*x 18、在MATLAB命令窗口中键入命令A=1 4 2;3 1 2;6 1 5;det(A(1:2,2:3).*A(2:3,2:3)。结果是（C ） （A）ans= -143 （B）ans= 60 （C）ans= 36 （D）ans= -19 19、MATLAB命令x = 3: 2: 100 将创建等差数列，该数列是（ C） （A）以3为初值的98个数,； （B）以100为终值的98的个数； （C）以99为终值的49个数； （D）以3为初值的97个数。20、在MATLAB命令窗口中输入命令data=4 1 2 3 1 3 1 4 2 4;y=hist(data,4)，

9、结果是（C） （A）y= 4 1 2 3；（B）y=3 2 3 2；（C）y= 3 2 2 3 ；（D）y= 4 2 1 121、MATLAB语句x,y=meshgrid(-2:2) 的数据结果中（D ） （A）x是行向量，y是列向量； （B）x是五行五列的矩阵； （C）x是行元素相同的矩阵； （D）x是列向量相同的矩阵22、MATLAB的语句colormap(0 0 1) （D ） （A）将三维网面图确定为红色； （B）将三维网面图确定为绿色； （C）将三维网面图确定为蓝色； （D）语句使用格式错误23、火炮发射炮弹的初始速度和发射角为已知，由此可估算出炮弹在空中的飞行时间Tfly，使用语句

10、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; (B) fai1=118;fai2=-76;Dfai=(fai1+fai2)/15; (C)

11、fai1=118;fai2=-76;Dfai=(fai1-fai2)/24; (D) fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。25、用MATLAB随机产生60个1到365之间的正整数，应该使用下面的哪一条命令 (B)A）fix(365*rand(1,60)； B）1+fix(366*rand(1,60)；C）1+fix(364*rand(1,60)； D）1+fix(365*rand(1,60)二、程序阅读1、关于“牟合方盖”的实验程序如下h=2*pi/100;t=0:h:2*pi;r=0:1;x=r*cos(t);y=r*sin(t);z=sqrt(1-x

12、.2); %第三行meshz(x,y,z),axis offcolormap(0 0 1)view(-47,56),hold onx1=cos(t);y1=sin(t);z1=abs(sin(t);plot3(x1,y1,z1,ro)； （1）下面有关程序的功能的说法确切的是（ ）（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的交线及所围区

13、域的边界曲面的上半部分。（2）关于第三行语句错误的解释是（ ）（A）z是矩形域上曲顶柱面高度值； （B）z是与y同型的矩阵； （C）z是圆域上曲顶柱面高度值； （D）z是与x同型的矩阵 2、 非负函数 y =f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴的旋转曲面，其体积.syms a b xf=exp(a*x)*sin(b*x);f1=subs(f,a,;f2=subs(f1,b,;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=

14、r*cos(theta); %第十行z=r*sin(theta); %第十一行mesh(z,y,x)colormap(0 0 0)axis offview(-17,54)（1）关于程序的功能确切的说法（ ）(A)计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积(B)计算曲线段 f(x)=expsin, 绕X轴旋转的旋转曲面体积 (C) 计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积并绘图（D）计算曲线段 f(x)=expsin, 绕X轴旋转的旋转曲面体积并绘图（2）由第九行至第十一行语句可得旋转曲面的方程（ ） (A) （B

15、） （C） (D)3数学实验程序如下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); syms e2 t f=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）卫星运行的近地速度和

16、远地速度； （C）卫星运行时向径每秒扫过的面积；（D）卫星运行的平均速度数据（2）第九行语句的功能是（ ） （A）计算卫星运行的最小速度； （B）计算卫星运行时向径每秒扫过的面积; （C）计算卫星运行的最大速度； （D）计算卫星运行轨道的平均速度4、Viviani体是圆柱体被球面所割立体。下面的数学实验程序功能是取R=2求体积上半部分，先利用符号处理重积分并转换为数值数据，再用蒙特卡罗方法计算体做对比。数学实验程序如下：syms x y;f=sqrt(4-x2-y2);y1=-sqrt(2*x-x2);y2=sqrt(2*x-x2);S1=int(f,y,y1,y2);S2=int(S1,x,

17、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);V1=8*length(II)/10000(1) 符号计算所用的积分公式是( A )(A) (B) (C) (D) (2) 蒙特卡罗方法选用的随机点变化范围的立方体区域是( B )(A) ；(B) (C) (D) 5、某厂生产两种产品，产一吨甲产品用A资源3吨、B资源4m3；产一吨乙产品用A资源2吨，B资源6m3，C资源7个单位。一吨甲产品和乙产品分别价值7万元和5万元，三种资源限制

18、分别为90吨、200m3和210个单位。生产两种产品使总价值最高的生产方案可用数学实验程序计算。C=-7,-5;A=3 2;4 6;0 7;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年一大轮回，按天干“甲乙

19、丙丁戊已庚辛壬癸”和地支“子丑寅卯辰巳午未申酉戍亥”循环排列而成。已知2009年是农历已丑年，通过简单计算可以找出年份与天干/地支对应的规律。下面数学实验程序对输入年份，计算并输出字符串农历纪年。填空完善程序。function calendar=year(year)if nargin=0, year=2009;endS1= 甲乙丙丁戊已庚辛壬癸;S2=子丑寅卯辰巳午未申酉戍亥;k1= mod(S1-4,10)+1 ; %定位天干序数%2010庚s1=S1(k1);k2= mod(S2-4,12) +1 ; %定位地支序数s2=S2(k2);calendar=strcat(int2str(yea

20、r),年是,s1,s2,年)2对于任意正整数n，如果n 只能被1和它自身整除，则称这个数为素数(或质数)。判素数程序的算法思想是试商法，即用2，3,，（n-1）去除n，如果能被这些数中一个整除，则n不是素数，否则是素数。完成下面填空。n=input(input n:=); for k=2:n if mod(n,k)= 0,break,end end if k if f=1 tk=(B-A)/(vb+vc); else tk= (B-A)/(vc+va); ; %计算A与C相遇时间 end A= A+va*tk ; %计算A点位置 B= B-vb*tk ; %计算B点位置 f=-f; k=k+1

21、;end5、二阶正交矩阵作用于某一向量时，其效果是将该向量旋转，旋转解为(逆时针旋转为正)。把一个以原点为中心的正三角形旋转，并缩小90%，迭代33次创建图3。完成程序填空：bata=1/2;7/6;11/6;15/6*pi;x=cos(bata);y=sin(bata);line(x,y)xy=x,y;alfa=pi/50; A=cos(alfa),-sin(alfa);sin(alfa),cos(alfa);for k=1:33 xy= *xy*A ; x=xy(:,1); y= xy(:,2) ; line(x,y)end6、长征三号甲运载火箭提供给探月卫星的初始速度不足以将卫星送往月球

22、轨道。为提高到奔月速度，中国航天工程师使用了卫星变轨技术。数学实验程序根据变轨中轨道周期和近地点距离数据，利用开普列第二定律模拟计算计算卫星飞行的最大速度。填空完善下面实验程序。R=6378;Time=16,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小时轨道，第

23、二次卫星进入24小时轨道，第三次卫星进入48小时轨道，第四次卫星进入116小时地月转移轨道。上面小时数并不是准确轨道周期，变轨目的是将速度从(km/s)逐渐提高到约(km/s)。下面数学实验程序是在区间,上插入线性等分点，即每个轨道的最大速度以等差数列出现，然后近似计算出每个轨道的周期参数。填空完善程序。function satel1()R=6378;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.2-c.2) ; %计算短半轴数据E2=(c./a).2;L=2*

24、pi*a.*(1-E2/4-3*E2.2/64)format bankVmax=linspace,5)S= *Vmax.*(R+h) ; %根据最大速度计算每秒钟扫过的面积Times=a.*b.*pi./S;myTimes=Times/36007、五月十二日以来汶川地区发生五级以上地震已经越过了25次，将地震数据整理，表示成n行三列的矩阵data，第二列为经度，第一列为为纬度，第三列为震级。不同震发点和震级的地震都对周边地区产生影响，距离近则影响强烈，距离远则影响减弱。用地震影响曲面描述，其数学原理如下：以每次震发点的经纬度为中心构造函数将这一函数离散化为矩阵并逐次累加，最后除以累加后的矩阵的最大值，可绘出地震影响曲面如图5。数学实验程序如下，请填空完善。load ()d=data(:,3);n=length(d); x=data(:,2);y=data(:,1);X,Y=meshgrid(100:110,30:35);Z=zeros(size(X);for k=1:n xk=x(k);yk=y(k); dk=d(k); Z=Z+dk*exp(-(X-xk).2+(Y-yk).2)/dk); endMaxz= max(Z); Z=Z./Maxz;mesh(X,Y,Z)colormap(0 0 0)