《电子科大成都学院211年数学实验复习题.docx》由会员分享，可在线阅读，更多相关《电子科大成都学院211年数学实验复习题.docx（9页珍藏版）》请在三一文库上搜索。
1、201年数学实验复习题一、选择题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.A2);的功能是(B )(A) 统计2000个随机点中落入特殊区域的点的索引值；(B) 统计1000个随机点落入特殊区域的点的索引值；(C) 模拟2000个随机点落入特殊区域的过程；(D) 模拟1000个随机点落入特殊区域的过
2、程。3、MATLA时算二项分布随机变量分布律的方法是(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)A2);S=int(f,t,0,pi/2)功能是(D)(A) 计算 f(x) 在 0,pi/2 上的积分；(B) 计算 f(t) 不定积分符号结果；(C) 计算 f(x) 积分的数值结果；(D) 计算 f(t) 定积分的符号结果。5、 y=dsolve( Dy=1/(1+xA2) - 2*yA2 , y(0)=
3、0 , x );ezplot(y) 的功能是( A )(A) 求微分方程特解并绘图；(B) 解代数方程(C) 求定积分；(D) 求微分方程通解。6、 X=10000 ;0.5*asin(9.8*X/(515A2) 的功能是计算关于抛射体问题的(A)(A) 十公里发射角；(B) 十公里飞行时间；(C) 最大飞行时间；(D) 最大射程。7、 theta=linspace(0,2*pi,100);r=cos(4*theta) ;polar(theta,r, k )功能是(D)8、MATLAB令 A=rand(5,5)A=(A) 绘四叶玫瑰线；(B) 绘三叶玫瑰线；(C) 绘心脏线；(D) 绘八叶玫瑰
4、线。a )5 5 max? |aj |用(a ) i=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);(C) y=sort(x);z=y(10:-1;1);(D) y,II=so
5、rt(x);z=x(II);11、MATLAB令 roots(1,0,0,-1)的功能是(D )(A) 产生向量1,0,0,1；(B)求方程 x +1=0 的根；(C) 求多项式 x -1 的值(D)3求方程 x - 1 = 0 的根。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);的结果是( B )13、命令 factor() 用于分解因式，syms x; f=4*xA3+9*xA2-30*x; factor
6、(diff(f) (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)功能是(C )14、MATLAB令 syms x; f=sin(x); V=pi*int(f*f,x,0,2*pi)(A)绘出函数f在0,2 p图形； (B)(C)计算旋转曲面所围的体积；(D)15、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”(A) k指向第二动物牛；(B) k(C) k指向第四动物兔；(D) k计算函数f在0,2 p 的积分；计算旋转曲面的表面积。,命令 k=mod(2008,12)+1 的结果是( D )指向第三
7、动物虎；指向第五动物龙。16.下面有关MATLA酸量名和函数名的说法中，错误的说法是(D )(A)变量名的第一个字符必须是一个英文字母(B)变量名可由英文字母、数字和下划线混合组成(C)变量名不得包含空格和标点，但可以有下连字符(D)变量名和函数名对于英文的大小使用没有区别17、在MATLA命令窗口中，键入命令 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)*sin(x)+1/2*x18、在 M
8、ATLA命令窗口中键入命令 A=1 4 2;3 1 2;6 1 5;det(A(1:2,2:3).*A(2:3,2:3)。结果是(C )(D) ans= -143(B) ans= 60(C) ans= 36( D) ans= -1919、MATLA命令x = 3: 2:100将创建等差数列，该数列是( C)(A)以3为初彳！的98个数，；(B)以100为终值的98的个数；(。以99为终值的49个数；(D)以3为初值的97个数。20、在 MATLA命令窗口中输入命令 data=4 1 2 3 1 3 1 4 2 4;y=hist(data,4),结果是(C)(A) y= 4 1 2 3 ; (B
9、) y=3 2 3 2 ; (C) y= 3 2 2 3; ( D) y= 4 2 1 121、MATLAB1句x,y=meshgrid(-2:2)的数据结果中(D )(A) x是行向量，y是列向量；(B) x是五行五列的矩阵；(0 x是行元素相同的矩阵；(D) x是列向量相同的矩阵22、MATLA的语句 colormap(0 0 1)(D )(A)将三维网面图确定为红色；(B)将三维网面图确定为绿色；(。将三维网面图确定为蓝色；(D)语句使用格式错误23、火炮发射炮弹的初始速度和发射角为已知，由此可估算出炮弹在空中的飞行时间Tfly ,使用语句Tspan=Tfly*(0:20)/20,将获得
10、一些数据，下面不正确的说法是( D )A Tspan为包括发射时刻在内的炮弹在空间飞行的21个不同的飞行时刻;B) Tspan中任意两个相邻数据之差的绝对值相等；0 Tspan包含了 21个数据，第一个数据为0,最后一个数据为Tfly ;D) Tspan是一个等差数列，公差为 Tfly/21 24、北京和纽约的经度分别是：东经 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=
11、(fai1-fai2)/24; (D) fai1=118;fai2=-76;Dfai=(fai1-fai2)/1525、用MATLAB机产生60个1到365之间的正整数，应该使用下面的哪一条命令(B)A fix(365*rand(1,60); B ) 1+fix(366*rand(1,60)0 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:0.05:1;x=r*cos(t);y=r*sin(t);z=sqrt(1-x.A2);%第三行meshz(x,y
12、,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+(B)绘圆柱面x2+y2=1,x2+(C)绘圆柱面x2+y2=1,x2+(D)绘圆柱面x2+y2=1,x2+(2)关于第三行语句错误的解释是(z2 = 1z2 = 1z2 = 1z2 = 1的交线；所围区域的边界曲面；的交线及所围区域的边界曲面；的交线及所围区域的边界曲面的上半部分。)(A) z是矩形域上曲顶柱面高度值;
13、(C) z是圆域上曲顶柱面高度值；(B) z是与y同型的矩阵;(D) z是与x同型的矩阵2、非负函数y =f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕 x轴旋转时，产生以x为对称轴的旋转2P2曲面，其体积 V = p Q f (x) dx.syms a b x f=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);
14、 %第十行z=r*sin(theta); %第 H一行mesh(z,y,x) colormap(0 0 0) axis off view(-17,54)(1) 关于程序的功能确切的说法()(A)计算曲线段 f(x)=exp(a x)sin(b x), (B)计算曲线段 f(x)=exp(-0.2x)sin(0.5x), (C) 计算曲线段 f(x)=exp(a x)sin(b x),(D)计算曲线段 f(x)=exp(-0.2x)sin(0.5x),(2) 由第九行至第十一行语句可得旋转曲面的方程(A)(B)(C)f ( x) (D)J y2 + z2 +t2223.数学实验程序如丫= f (
15、 x +z )h=439;H=2384;R=6400;绕X轴旋转的旋转曲面体积绕X轴旋思勺看x曲裾积0,绕x轴旋转的旋晶稀积2P图绕X轴旋转的旋转曲步)2py = f (. x2 + z2+t2)a=(h+H+2*R)/2;c=(H-h)/2;e1=c/a; b=sqrt(a*a-c*c);syms e2 tf=sqrt(1-e2*cos(t)A2);ft=subs(f,e2,e1*e1);S=int(ft,0,pi/2);V1=8*length(II)/10000(1)符号计算所用的积分公式是(A )y2dyL=4*a*double(S);V=L/(114*3600) ; % 第九行 s1=
16、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)计算卫星运行轨道的平均速度222 . .2.2.224、Viviani 体是圆柱体 (x- R/2) +y ? R /4被球面 x + y +z = R所割立体。下面的数学实验程序功能是取 R=2求体积上半部
17、分，先利用符号处理重积分并转换为数值数据，再用蒙 特卡罗方法计算体做对比。数学实验程序如下：syms x y;f=sqrt(4-xA2-yA2);y1=-sqrt(2*x-xA2);y2=sqrt(2*x-xA2);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).A2+Y.A2=1&Z=sqrt(4-X.A2-Y.A2);(B )+ -j2x- x，22(C) V = O 2x- x2 4 - x - y dy(2)蒙特卡罗
18、方法选用的随机点变化范围的立方体区域是(A) W=( x,y,z)|x挝(0,2), y (0,2), z? (0,2)；(B) W=( x, y, z)| x挝(0,2), y (-1,1),z? (0,2)(C) W=( x, y, z)|0x2),0 y 1,0z2(D) W=( x, y, z) | 0 x 2),0 y 2,0 z25、某厂生产两种产品，产一吨甲产品用A资源3吨、B资源4谓；产一吨乙产品用 A资源2吨，B资源6+，C资源7个单位。一吨甲产品和乙产品分别价值7万元和5万元，三种资源限制分别为90吨、200m3和210个单位。生产两种产品使总价值最高的生产方案可用数学实验
19、程序计算。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）;程序中变量C表示（）（A）目标函数系数；（B）等式约束系数；（C）不等式约束系数；（D）等式约束常向量程序中变量A表示（）（A）等式约束矩阵；（B）不等式约束矩阵；（C）决策变量的值；（D）目标函数的最大值三、程序填空1、中国农历60年一大轮回，按天干“甲乙丙丁戊已庚辛壬癸”和地支“子丑寅卯辰巳午未申酉戍亥”循环排列而成。已知2009年是农历已丑年，通过简单计算可以找出年份与天干/地支对应
20、的规律。下面数学实验程序对输入年份，计算并输出字符串农历纪年。填空完善程序。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（year）, 年是,s1,s2,年）2 .对于任意正整数n,如果n只能被1和它自身整除，则称这个数为素数（或质数）。判素数程序的算法思想是试商法，即
21、用2, 3,，（n-1 ）去除n,如果能被这些数中一个整除，则 n不是素数，否则是素数。完成下面填空。n=input（input n:=）;for k=2:nif mod（n,k）= 0,break,end endif k 0.2 if f=1tk=(B-A)/(vb+vc);elsetk= (B-A)/(vc+va);% 计算A与C相遇时间end A= A+va*tk；% 计算 A点位置B= B-vb*tk；% 计算B点位置 f=-f； k=k+1;endCL5、二阶正交矩阵作用于某一向量时，其效果是将该向量旋转，旋转解为(逆时针旋转为正)。把一个以原点 为中心的正三角形旋转 二 /50,并
22、缩小90%迭代33次创建图3。完成程序填空:bata=1/2;7/6;11/6;15/6*pi;-0.6-1-0.8-0.6-0.4-0.200.20.40.60.81图3旋转三角形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:33xy= 0.9*xy*A；x=xy(:,1);y= xy(:,2);line(x,y)end 6、长征三号甲运载火箭提供给探月卫星的初始速度不足以将卫星送往月球轨道。为提高到奔月速度，中国航天工程师使用了卫星
23、变轨技术。数学实验程序根据变轨中轨道周期和近地点距离数据，利用开普列第二定律模拟计算计算卫星飞行的最大速度。填空完善下面实验程序。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
24、小时轨道，第三次卫星进入48小时轨道，第四次卫星进入116小时地月转移轨道。上面小时数并不是准确轨道周期，变轨目的是将速度从10.3(km/s)逐渐提高到约10.9(km/s)。下面数学实验程序是在区间10.3,10.9上插入线性等分点，即每个轨道的最大速度以等差数列出现，然后近似计算出每个轨道的周期参数。填空完善程序。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.A2-c.A2)；%计算短半轴数据E2=(c./
25、a)A2;L=2*pi*a.*(1-E2/4-3*E2 A2/64) format bankVmax=linspace(10.3,10.9,5)S= 0.5*Vmax.*(R+h)；根据最大速度计算每秒钟扫过的面积Times=a.*b.*pi./S;myTimes=Times/36007、五月十二日以来汶川地区发生五级以上地震已经越过了25次，将地震数据整理，表示成n行三列的矩阵data,第二列为经度 X k ,第一列为为纬度y k ,第三列为震级d k。不同震发点和震级的地震都对周边地区产生影响，距离近则影响强烈，距离远则影响减弱。用地震影响曲面描述，其数学原理如下：以每次震发点的经纬度为中
26、 心构造函数5。数学实验程将这一函数离散化为矩阵并逐次累加，最后除以累加后的矩阵的最大值，可绘出地震影响曲面如图 序如下，请填空完善。load ( data.txt )110d=data(:,3);n=length(d); x=data(:,2); y=data(:,1);X,Y=meshgrid(100:02110,30:0235); Z=zeros(size(X); for k=1:n xk=x(k);yk=y(k); dk=d(k);Z=Z+dk*exp(-(X-xk).A2+(Y-yk)A2)/dk);end Maxz= max(Z);图5地震影响曲面Z=Z./Maxz;mesh(X,Y,Z) colormap(0 0 0)