1、系别 班次 学号 姓名 . 密封线以内答题无效电子科技大学成都学院二零壹零至二零壹壹学年第一学期 数学实验课程考试题A（120分钟） 闭卷 考试时间： 成绩构成：平时（50%）+期末（50%）=总成绩（100%） 一二三合计注意：请同学们将答案填写在答题纸上，否则无效。一、单项选择题（共40分）1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用( ) (A) diag(magic(3); (B) diag(magic); (C) diag(diag(magic(3); (D) diag(diag(magic)。2、符号计算与一般数值计算有很大区别，它将得到准确的符号表达式。在MAT
2、LAB命令窗口中键入命令syms x y1=x2；y2=sqrt(x)；int(y1-y2,x,0,1) 屏幕显示的结果是（ ）（A）y1 =x(1/2) （B）ans= 1/3； （C）y2 =x2； （D）ans= -1/33、MATLAB命令syms x; f=sin(x); V=pi*int(f*f,x,0,pi)功能是( )(A) 绘出函数f在0,2图形； (B) 计算函数f在0,2的积分；(C) 计算旋转曲面所围的体积； (D) 计算旋转曲面的表面积。4、用赋值语句给定x数据，计算对应的MATLAB表达式是（ ）（A）sqrt(7sin(3+2x)+exp(3)log(3) （B）
3、sqrt(7*sin(3+2*x)+exp(3)*log(3)（C）sqr(7*sin(3+2*x)+e3*log(3) （D）sqr(7sin(3+2x)+ e3log(3)5、MATLAB语句strcat(int2str(2010),年是, s,年)的功能是（ ）（A）将数据2010转换为符号； （B）将数据2010与符号变量合并；（C）将几个符号变量合并为一个； （D）将符号变量转换为数值变量；6、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令k=mod(2008,12)+1的结果是( )(A) k指向第二动物牛； (B) k指向第三动物虎；(C) k指向第四动物兔； (D) k指向第五动
4、物龙。7、下面有关MATLAB函数的说法，哪一个是错误的（ ）（A）函数文件的第一行必须由function开始，并有返回参数，函数名和输入参数；（B）MATLAB的函数可以有多个返回参数和多个输入参数；（C）在函数中可以用nargout检测用户调用函数时的输入参数个数；（D）如果函数文件内有多个函数，则只有第一个函数可以供外部调用；8、MATLAB命令x=rand(10,1)生成10个随机数，将它们从大到小排序，使用( )(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=
5、sort(x);z=x(II);9、MATLAB绘三维曲面需要构建网格数据，语句x,y=meshgrid(-2:2)返回数据中（ ）（A）x是行向量，y是列向量； （B）x是列向量，y是行向量；（C）y是行元素相同的矩阵； （D）x是行向量相同的矩阵10、设a,b,c表示三角形的三条边，表达式a+bc|a+cb|b+ca，（ ）（A）是三条边构成三角形的条件； （B）是三条边不构成三角形的条件；（C）构成三角形时逻辑值为真； （D）不构成三角形时逻辑值为假二、程序阅读理解（24分）1、用plot命令绘制平面曲线v0=515;alpha=45*pi/180;T=2*v0*sin(alpha)/9
6、.8; %第二行t= T*(0:16) /16;x=v0 *t *cos(alpha);y=v0 *t *sin(alpha) -0.5*9.8*t.2;plot(x,y,x,y,r*)Xmax=x(17)00.511.522.53x 10402000400060008000（1）对下面有关程序的功能的说法确切的是（ ）（A）以515为初速度的抛射曲线的绘制； （B）以515为初速度的抛射曲线的绘制以及计算射程；（C）以515为初速度以45度为发射角的抛射曲线的绘制以及计算射程； （D）以515为初速度以45度为发射角的抛射曲线的绘制。（2）第二行语句的功能是（ ）（A）计算射程；（B）计算飞
7、行时间；（C）计算航点的横坐标； （D）计算航点的纵坐标。2、关于“牟合方盖”的实验程序如下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),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 + y
8、2 = 1, x2 + z2 = 1所围区域的边界曲面；（C）绘圆柱面x2 + y2 = 1, x2 + z2 = 1的交线及所围区域的边界曲面；（D）绘圆柱面x2 + y2 = 1, x2 + z2 = 1的交线及所围区域的边界曲面的上半部分。（2）关于第三行语句错误的解释是（ ）（A）z是矩形域上曲顶柱面高度值； （B）z是与y同型的矩阵； （C）z是圆域上曲顶柱面高度值； （D）z是与x同型的矩阵 3、 非负函数 y =f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴的旋转曲面，其体积.syms a b xf=exp(a*x)*sin(b*x);f1=s
9、ubs(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(0 0 0)axis offview(-17,54)（1） 关于程序的功能确切的说法（ ）(A)计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积(B)计算曲线段 f(x)=exp(-0.2x)sin(0.5x
10、), 绕X轴旋转的旋转曲面体积 (C) 计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积并绘图（D）计算曲线段 f(x)=exp(-0.2x)sin(0.5x), 绕X轴旋转的旋转曲面体积并绘图（2） 由第九行至第十一行语句可得旋转曲面的方程（ ）(A) （B） （C） (D)三、程序填空（36分）1、中国农历60年一大轮回，按天干“甲乙丙丁戊已庚辛壬癸”和地支“子丑寅卯辰巳午未申酉戍亥”循环排列而成。已知2009年是农历已丑年，通过简单计算可以找出年份与天干/地支对应的规律。下面数学实验程序对输入年份，计算并输出字符串农历纪年。填空完善程序。（8分）func
11、tion calendar=year(year)if nargin=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对于二重积分，D是直线y= x 2和抛物线 y2 = x 所围区域。下面数学实验程序的功能是利用MATLAB符号计算方法将二重积分处理为二次积分进行计算（选择先对x积分后对y积分的秩序），然后
12、将计算结果的符号表达式转换为双精度实数。（8分）syms x y;f=x*y*y;x1=y*y;x2= y+2 S1=int(f,x,x1,x2);S2=int(S1,y,-1,2);S= 3某年A、B两城镇人各有1000人，以后每年A镇人口10%迁往B镇; B镇人口15%迁往A镇。则有迁移矩阵L，L有两个互异特征值l1，l2，属于两个特征值的特征向量分别为：a 1，a 2，记初始人口分布X0=1000，1000T，则存在不全为零的数c1，c2使得 X0 = c1a 1 + c2a 2（解此方程组求出c1，c2）。下面实验程序主要功能是利用LnX0 = c1l1na 1+ c2l2na2计算人
13、口变化。（8分）n=input(input n:=);A=0.9,0.15;0.1,0.85;X0=1000;1000;P,D=eig(A);C= double(S2) c1=C(1);c2=C(2);alfa1=P(:,1);alfa2=P(:,2);lamda1=D(1,1);lamda2=D(2,2);Xn= PX0; 4、长征三号甲运载火箭提供给探月卫星的初始速度不足以将卫星送往月球轨道。为提高到奔月速度，中国航天工程师使用了卫星变轨技术。数学实验程序根据变轨中轨道周期和近地点距离数据，利用开普列第二定律模拟计算计算卫星飞行的最大速度。填空完善下面实验程序。（12分）R=6378;Ti
14、me=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= c1*lamdaln*alfa1+ c2*lamda2n*alfa2 b=sqrt(a.*a-c.*c);S= sqrt(a.*a-c.*c); Vmax=2*S./(R+h)Vmin= S= 0.5*Vmax.*(R+h) 电子科技大学成都学院二零壹零至二零壹壹学年第一学期 数学实验课程考试题B（120分钟） 闭卷 考试时间： 成绩构成：平时（50%）+期末（50%）=总成绩（100%）
15、 一二三合计注意：请同学们将答案填写在答题纸上，否则无效。一、单项选择题（共40分）1、下面有关MATLAB变量名和函数名的说法中，错误的说法是( )(A) 变量名的第一个字符必须是一个英文字母(B) 变量名可由英文字母、数字和下划线混合组成(C) 变量名不得包含空格和标点，但可以有下连字符(D) 变量名和函数名对于英文的大小使用没有区别2、符号计算与一般数值计算有很大区别，它将得到准确的符号表达式。在MATLAB命令窗口中键入命令syms x y1=x2；y2=sqrt(x)；int(y1-y2,x,0,1) 屏幕显示的结果是（ ）（A）y1 =x(1/2) （B）ans= 1/3； （C）
16、y2 =x2； （D）ans= -1/33、MATLAB命令syms x; f=sin(x); V=pi*int(f*f,x,0,pi)功能是( )(A) 绘出函数f在0,2图形； (B) 计算函数f在0,2的积分；(C) 计算旋转曲面所围的体积； (D) 计算旋转曲面的表面积。4、用赋值语句给定x数据，计算对应的MATLAB表达式是（ ）（A）sqrt(7sin(3+2x)+exp(3)log(3) （B）sqrt(7*sin(3+2*x)+exp(3)*log(3)（C）sqr(7*sin(3+2*x)+e3*log(3) （D）sqr(7sin(3+2x)+ e3log(3)5、MATL
17、AB语句strcat(int2str(2010),年是, s,年)的功能是（ ）（A）将数据2010转换为符号； （B）将数据2010与符号变量合并；（C）将几个符号变量合并为一个； （D）将符号变量转换为数值变量；6、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令k=mod(2008,12)+1的结果是( )(A) k指向第二动物牛； (B) k指向第三动物虎；(C) k指向第四动物兔； (D) k指向第五动物龙。7、MATLAB命令A=magic(3)创建3阶幻方矩阵，求A的特征值绝对值最小用( )(A) min(abs(eig(A); (B) min(eig(abs(A); (C)min(
18、eig(A); (D) min(abs(A);8、MATLAB命令x=rand(10,1)生成10个随机数，将它们从大到小排序，使用( )(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);9、MATLAB命令x,y=meshgrid(1:3);H=1./(x+y-1)产生的矩阵H是( )(A) (B) (C) (D) 10、设a,b,c表示三角形的三条边，表达式a+bc|a+cb|b+ca， （ ）（A）是三条边构成三角形的条件； （B）是
19、三条边不构成三角形的条件；（C）构成三角形时逻辑值为真； （D）不构成三角形时逻辑值为假二、程序阅读理解（24分）1、用plot命令绘制平面曲线v0=515;alpha=45*pi/180;T=2*v0*sin(alpha)/9.8; %第二行t= T*(0:16) /16;x=v0 *t *cos(alpha);y=v0 *t *sin(alpha) -0.5*9.8*t.2;plot(x,y,x,y,r*)Xmax=x(17)00.511.522.53x 10402000400060008000（1）对下面有关程序的功能的说法确切的是（ ）（A）以515为初速度的抛射曲线的绘制； （B）以
20、515为初速度的抛射曲线的绘制以及计算射程；（C）以515为初速度以45度为发射角的抛射曲线的绘制以及计算射程； （D）以515为初速度以45度为发射角的抛射曲线的绘制。（2）第二行语句的功能是（ ）（A）计算射程；（B）计算飞行时间；（C）计算航点的横坐标； （D）计算航点的纵坐标。2、关于“牟合方盖”的实验程序如下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),axis offcolormap(0 0 1)view(-47,56),hold onx1=cos(t);
21、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的交线及所围区域的边界曲面的上半部分。（2）关于第三行语句错误的解释是（ ）（A）z是矩形域上曲顶柱面高度值； （B）z是与y同型的矩阵； （C）z是圆域上曲顶柱面
22、高度值； （D）z是与x同型的矩阵 3、 非负函数 y =f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴的旋转曲面，其体积syms a b xf=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(0 0 0)ax
23、is offview(-17,54)（3） 关于程序的功能确切的说法（ ）(A)计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积(B)计算曲线段 f(x)=exp(-0.2x)sin(0.5x), 绕X轴旋转的旋转曲面体积 (C) 计算曲线段 f(x)=exp(a x)sin(b x), 绕X轴旋转的旋转曲面体积并绘图（D）计算曲线段 f(x)=exp(-0.2x)sin(0.5x), 绕X轴旋转的旋转曲面体积并绘图（4） 由第九行至第十一行语句可得旋转曲面的方程（ ）(A) （B） （C） (D)三、程序填空（36分）1、判断闰年条件有二能被4整除，但不能被
24、100整除；能被4整除，又能被400整除。year=input(input year:=);n1=year/4;n2=year/100;n3=year/400;if n1=fix(n1) &n2=fix(n2), (不是闰年) disp(是闰年)elseif n1=fix(n1)&n3=fix(n3) disp(是闰年)else disp end2对于二重积分，D是直线y= x 2和抛物线 y2 = x 所围区域。下面数学实验程序的功能是利用MATLAB符号计算方法将二重积分处理为二次积分进行计算（选择先对x积分后对y积分的秩序），然后将计算结果的符号表达式转换为双精度实数。（8分）syms
25、x y;f=x*y*y;x1=y*y;x2= y+2 S1=int(f,x,x1,x2);S2=int(S1,y,-1,2);S= double(S2) 3某年A、B两城镇人各有1000人，以后每年A镇人口10%迁往B镇; B镇人口15%迁往A镇。则有迁移矩阵L，L有两个互异特征值l1，l2，属于两个特征值的特征向量分别为：a 1，a 2，记初始人口分布X0=1000，1000T，则存在不全为零的数c1，c2使得 X0 = c1a 1 + c2a 2（解此方程组求出c1，c2）。下面实验程序主要功能是利用LnX0 = c1l1na 1+ c2l2na 2计算人口变化。（10分）n=input(input n:=);A=0.9,0.15;0.1,0.85;X0=1000;1000;P,D=eig(A);C= PX0; c1=C(1);c2=C(2);alfa1=P(:,1);alfa2=P(:,2);lamda1=D(1,1);lamda2=D(2,2);Xn= c1*lamdaln*alfa1+ c2*lamda2n*alfa2 4、为了进入地月转移轨道，嫦娥一号卫星进行了四次变轨调速度。第
