《数学实验》期末总复习题

电子科技大学成都学院二零壹零至二零壹壹学年第一学期

数樊实验课程考试题A(120分钟)闭卷考试时间：

成绩构成：平时(50%)+期末(50%)=总成绩(100%)

一二三合计

注意：请同学们将答案填写在答题纸上，否则无效。

一、单项选择题(共40分)

1、三阶幻方又称为九宫图，提取三阶幻方矩阵对角元并构造对角阵用()

(A)diag(magic(3));(B)diag(magic);

(C)diag(diag(magic(3)));(D)diag(diag(magic))<>

2、符号计算与一般数值计算有很大区别，它将得到准确的符号表达式。在MATLAB命令窗口中键入命令

symsxyl=x0；y2=sqrt(x)；int(yl-y2,x,0J)屏幕显示的结果是()

(A)yl=xA(l⑵(B)ans=1/3；(C)y2=xA2；(D)ans=-1/3

3、MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,0,pi)功能是()

(A)绘出函数f在[0,27]图形；(B)计算函数f在[0,2%]的积分；

(C)计算旋转曲面所围的体积；(D)计算旋转曲面的表面积。

4、用赋值语句给定x数据，计算47sin(3+2x)+e3In3对应的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)+eA3*log(3))(D)sqr(7sin(3+2x)+eA31og(3))

5、MATLAB语句strcat(int2str(2010),年是1s,年)的功能是()

(A)将数据2010转换为符号；(B)将数据2010与符号变量合并；

(C)将几个符号变量合并为一个；(D)将符号变量转换为数值变量；

6、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪"，命令k=mod(2008,12)+1的结果是()

(A)k指向第二动物牛；(B)k指向第三动物虎；

(C)k指向第四动物兔；(D)k指向第五动物龙。

7、下面有关MATLAB函数的说法，哪一个是错误的()

(A)函数文件的第一行必须由function开始，并有返回参数，函数名和输入参数；

(B)MATLAB的函数可以有多个返回参数和多个输入参数；

(C)在函数中可以用nargout检测用户调用函数时的输入参数个数；

(D)如果函数文件内有多个函数，则只有第一个函数可以供外部调用：

8、MATLAB命令x=rand(10,l)生成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(-2:2)返回数据中()

(A)x是行向量，y是列向量；(B)x是列向量，y是行向量；

(C)y是行元素相同的矩阵；(D)x是行向量相同的矩阵

10、设a,b,c表示三角形的三条边，表达式a+b>c|a+c>b|b+c>a,()

(A)是三条边构成三角形的条件；(B)是三条边不构成三角形的条件；

(C)构成三角形时逻辑值为真；(D)不构成三角形时逻辑值为假

二、程序阅读理解(24分)

1、用plot命令绘制平面曲线

v0=515;alpha=45*pi/l80;

T=2*v0*sin(alpha)/9.8;%第二行

t=T*(0:16)/16;

x=v0*t*cos(alpha);

y=vO*t*sin(alpha)-0.5*9.8*t.A2;

plot(x,y,x,y,，r*，)

Xmax=x(17)

8000

6000

4000

2000

0

00.511.522.53

(1)对下面有关程序的功能的说法确切的是()

(A)以515为初速度的抛射曲线的绘制；

(B)以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(l-x,A2);%第三行

meshz(x,y,z),axisoff

colormap([001])

view(-47,56),holdon

x1=cos(t);y1=sin(t);z1=abs(sin(t));

plot3(xl,yl,zl,'ro')；

(1)下面有关程序的功能的说法确切的是(

(A)绘圆柱面N+y2=i,"+z2=1的交线；

(B)绘圆柱面12+$=1,/+22=1所围区域的边界曲面；

(C)绘圆柱面N+y2=i,N+z2=l的交线及所围区域的边界曲面；

(D)绘圆柱面》2+$=i,#+z2=1的交线及所围区域的边界曲面的上半部分。

(2)关于第三行语句褚误的解释是()

(A)z是矩形域上曲顶柱面高度值；(B)z是与y同型的矩阵；

(C)z是圆域上曲顶柱面高度值；(D)z是与x同型的矩阵

3、非负函数y=f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴

的旋转曲面，其体积V=可："(x)F办.

symsabx

f=exp(a*x)*sin(b*x);

fl=subs(f,a,-0.2);

f2=subs(fl,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])

axisoff

view(-17,54)

(1)关于程序的功能确切的说法()

(A)计算曲线段f(x)=exp(ax)sin(bx),绕X轴旋转的旋转曲面体积

(B)计算曲线段f(x)=exp(-0.2x)sin(0.5x),0<x<2^绕X轴旋转的旋转曲面体积

(C)计算曲线段f(x)=exp(ax)sin(bx),()4％<2左绕X轴旋转的旋转曲面体积并绘图

(D)计算曲线段f(x)=exp(-0.2x)sin(0.5x),0<x<27绕X轴旋转的旋转曲面体积

并绘图

(2)由第9行至第十一行语句可得旋转曲面的方程.)

(A)f(x)=Jy+z+r7(B)f(x)=Jy2+z2

(C)=f(\lx2+z2)(D)y=/(Vx2+z2+/2)

三、程序填空(36分)

1、中国农历60年一大轮回，按天干“甲乙丙丁戊己庚辛壬癸”和地支“子丑寅卯辰巳午未申酉成亥”循环排

列而成。已知20XX年是农历己丑年，通过简单计算可以找出年份与天干/地支对应的规律。下面数学实验

程序对输入年份，计算并输出字符串农历纪年。填空完善程序。(8分)

functioncalendar=year(year)

ifnargin==0,year=2009;end

Sl=?甲乙丙丁戊已庚辛壬癸，；

S2d子丑寅卯辰巳午未申酉戍亥，；

k1=mod(year-4,l0)+1;①;%定位天干序数

sl=Sl(kl);

k2=mod(year-4,12)+1②；%定位地支序数

s2=S2(k2);

calendar=strcat(int2str(year),'年是',s1,s2,'年')

2.对于二重积分dxdy,。是直线y=x-2和抛物线y2=x所围区域。下面数学实验程序的

D

功能是利用MATLAB符号计算方法将二重积分处理为二次积分进行计算(选择先对X积分后对y积分的

秩序)，然后将计算结果的符号表达式转换为双精度实数。(8分)

symsxy;

f=x*y*y;

xl=y*y；x2=y+2①

Sl=int(f,x,xl,x2);

S2=int(Sl,y,-l,2);

S=②

3.某年A、B两城镇人各有1000人，以后每年A镇人口10%迁往B镇;B镇人口15%迁往A镇。则有迁

移矩阵LL有两个互异特征值;I],属于两个特征值的特征向量分别为：ai，。2,记初始人口分布

Xo=flOOO,lOOOf,则存在不全为零的数C1，62使得乂0=，1。1+。2。2(解此方程组求出口，C2)。下面

n

实验程序主要功能是利用LX0=C\Z"a\+。2丸2”&2计算人口变化。(8分)

n=input(*inputn:-);

A=[0.9,0.15;0.1,0.85];

X0=[1000;1000];

[P,D]=eig(A);

C=double(S2)①

cl=C(I);c2=C(2);

alfal=P(:,l);

alfa2=P(:,2);

lamdal=D(l,l);

lamda2=D(2,2);

Xn=P\X0;②

4、长征三号甲运载火箭提供给探月卫星的初始速度不足以将卫星送往月球轨道。为提高到奔月速度，中国

航天工程师使用了卫星变轨技术。数学实验程序根据变轨中轨道周期和近地点距离数据，利用开普列第二

定律模拟计算计算卫星飞行的最大速度。填空完善下面实验程序。(12分)

R=6378;Time=[l6,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*lamdalAn*alfal+c2*lamda2An*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%)

二三合计

注意：请同学们将答案填写在答题纸上，否则无效。

一、单项选择题(共40分)

1＞下面有关MATLAB变量名和函数名的说法中，错误的说法是()

(A)变量名的第一个字符必须是一个英文字母

(B)变量名可由英文字母、数字和下划线混合组成

(C)变量名不得包含空格和标点，但可以有下连字符

(D)变量名和函数名对于英文的大小使用没有区别

2、符号计算与一般数值计算有很大区别，它将得到准确的符号表达式。在MATLAB命令窗口中键入命令

symsxyl=xA2；y2=sqrt(x)；int(yl-y2,x,0,l)屏幕显示的结果是()

(A)yl=xA(l⑵(B)ans=l/3；(C)y2=xA2；(D)ans=-l/3

3、MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,O,pi)功能是()

(A)绘出函数f在［0,27］图形；(B)计算函数f在［0,27］的积分；

(C)计算旋转曲面所围的体积；(D)计算旋转曲面的表面积。

4、用赋值语句给定x数据，计算gsin(3+2x)+e31n3对应的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)+eA3*log(3))(D)sqr(7sin(3+2x)+eA310g(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指向第五动物龙。

7、MATLAB命令A=magic(3)创建3阶幻方矩阵，求A的特征值绝对值最小用()

(A)min(abs(eig(A)));(B)min(eig(abs(A)));(C)min(eig(A));(D)min(abs(A));

8、MATLAB命令x=rand(10,1)生成10个随机数，将它们从大到小排序，使用()

(A)y=sort(x);z=y(10:l);(B)[y,H]=sort(x);z=y(n);

(C)y=sort(x);z=y(10:-1;1);(D)[y,II]=sort(x);z=x(II);

9、MATLAB命令[x,y]=meshgrid(l:3);H=l./(x+y-l)产生的矩阵H是()

-11r-123-「123'■11/21/3

(A)222(B)123(C)233(D)1/21/31/4

3331233451/31/41/5

10、设a,b,c表示三角形的三条边，表达式a+b>c|a+c>b|b+c>a,()

(A)是三条边构成三角形的条件；(B)是三条边不构成二角形的条件;

(C)构成三角形时逻辑值为真；(D)不构成三角形时逻辑值为假

二、程序阅读理解(24分)

1、用plot命令绘制平面曲线

v0=515;alpha=45*pi/l80;

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.A2;

plot(x,y,x,y/r*，)

Xmax=x(17)

8000

6000

4000

2000

0

00.511.522.53

(1)对下面有关程序的功能的说法确切的是()

(A)以515为初速度的抛射曲线的绘制；

(B)以515为初速度的抛射曲线的绘制以及计算射程；

(C)以515为初速度以45度为发射角的抛射曲线的绘制以及计算射程;

(D)以515为初速度以45度为发射角的抛射曲线的绘制。

(2)第二行语句的功能是()

(A)计算射程；(B)计算飞行时间；

(C)计算航点的横坐标；(D)计算航点的纵坐标。

2、关于“牟合方盖”的实验程序如下

h=2*pi/100;t=0:h:2*pi;

1-0:0.05:1;x=r'*cos(t);y=r'*sin(t);

z=sqrt(l-x.A2);%第三行

meshz(x,y,z),axisoff

colormap([001])

view(-47,56),holdon

x1=cos(t);yl=sin(t);z1=abs(sin(t));

plot3(xl,yl,zl,,ro,)；

(1)下面有关程序的功能的说法确切的是()

(A)绘圆柱面N+y2=1,N+Z2=1的交线；

(B)绘圆柱面N+y2=I,x2+z2=1所围区域的边界曲面；

(C)绘圆柱面N+y2=I,x2+z2=1的交线及所围区域的边界曲面；

(D)绘圆柱面H+y2=i,H+z2=1的交线及所围区域的边界曲面的上半部分。

(2)关于第三行语句错误的解释是()

(A)z是矩形域上曲顶柱面高度值；(B)z是与y同型的矩阵；

(C)z是圆域上曲顶柱面高度值；(D)z是与x同型的矩阵

3、非负函数y=f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴

的旋转曲面，其体积.=»r"(切2小

symsabx”

f=exp(a*x)*sin(b*x);

fl=subs(f,a,-0.2);

f2=subs(fl,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=f*ones(size(t));%第九行

y=rf*cos(theta);%第十行

z=rr*sin(theta);%第十一行

mesh(z,y,x)

colormapdO00])

axisoff

view(-17,54)

(3)关于程序的功能确切的说法()

(A)计算曲线段f(x)=exp(ax)sin(bx),0＜兀42万绕X轴旋转的旋转曲面体积

(B)计算曲线段f(x)=exp(-0.2x)sin(0.5x),0MxM2万绕X轴旋转的旋转曲面体积

(C)计算曲线段f(x)=exp(ax)sin(bx),0«左427绕X轴旋转的旋转曲面体积并绘图

(D)计算曲线段f(x)=exp(-0.2x)sin(0.5x),04x42%绕X轴旋转的旋转曲面体积

并绘图

(4)由第到行至第十一行语句可得旋转曲面的方程.)

(A)/(x)=Jy2+z2+fr(B)f(x)=Jy2+z2

(C)y=/(V%2+z2)(D)y=/(A/X2+z2+t2)

三、程序填空(36分)

1、判断闰年条件有二

①能被4整除，但不能被100整除:

②能被4整除，又能被400整除。

yeai^inputC'inputyear:=,);

nl=year/4;

n2=year/100;

n3=year/400;

ifnl===fix(nl)&n2~=fix(n2),，不是闰年，)

_____________________①

dispC是闰年)

elseifn1==fix(nI)&n3==fix(n3)

dispC是闰年)

else

disp②

end

2.对于二重积分JJAy2dxdy,。是直线y=x-2和抛物线/=》所围区域。下面数学实验程序的

D

功能是利用MATLAB符号计算方法将二重积分处理为二次积分进行计算

(选择先对X积分后对y积分的秩序)，然后将计算结果的符号表达式转换

为双精度实数。(8分)

symsxy;

f=x*y*y;

x1=y*y;x2=y+2

Sl=int(f,x,xl,x2);

S2=int(Sl,y,-l,2);

S=double(S2)

__________________________②

3.某年A、B两城镇人各有1000人，以后每年A镇人口10%迁往B镇;B镇人口15%迁往A镇。则有迁

移矩阵L,L有两个互异特征值为，/U，属于两个特征值的特征向量分别为：ai，&2,记初始人口分布

r

X()=[1000,1000],则存在不全为零的数ci，C2使得X°=ciai+c2a2(解此方程组求出ci，c2)o下面

n

实验程序主要功能是利用LX0=ci九"a1+C2/l2〃a2计算人口变化。(10分)

n=input(finputn:-);

A=[0.9,0J5;0.1,0.85];

X0=[1000;I000];

[P,D]=eig(A);

C=P\X0;①

cl=C(l);c2=C(2);

alfal=P(:,l);

alfa2=P(:,2);

lamdal=D(l,l);

lamda2=D(2,2);

Xn=cl*lamdalAn*alfal+c2*lamda2An*alfa2

_________________________②

4、为了进入地月转移轨道，嫦娥一号卫星进行了四次变轨调速度。第一次变轨从16小时初始轨道进入16

小时轨道，第二次卫星进入24小时轨道，第三次卫星进入48小时轨道，第四次卫星进入116小时地月转

移轨道。上面小时数并不是准确轨道周期，变轨目的是将速度从10.3(km/s)逐渐提高到约10.9(km/s)。下面

数学实验程序是在区间[10.3,10.9]上插入线性等分点，即每个轨道的最大速度以等差数列出现，然后近似计

算出每个轨道的周期参数。(10分)

functionsatel1()

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=(H-h)/2;①；％计算短半轴数据

E2=(c./a).A2;

L=2*pi*a.*(1-E2/4-3*E2.A2/64)

formatbank

Vmax=linspace(10.3,10.9,5)

S=pi*a*b/Times②;%根据最大速度计算每秒钟扫过的面积

Times二a.*b.*pi./S;

myTimes=Times/3600

电子科技大学成都学院二零壹零至二零壹壹学年第二学期

数学实验课程考试题A(120分钟)闭卷考试时间：

成绩构成：平时(50%)+期末(50%)=总成绩(100%)

二三合计

注意：请同学们将答案填写在答题纸上，否则无效。

一、单项选择题(共40分)

1、下面有关MATLAB变量名和函数名的说法中，正确的说法是()

(A)变量名的第一个字符必须是一个英文字母

(B)变量名不可由英文字母、数字和下划线混合组成

(C)变量名可以包含空格、标点和下连字符

(D)变量名和函数名对于英文的大小使用没有区别

2、MATLAB命令x=5:2:100将创建等差数列，该数列是()

(A)以5为初值的97个数,；(B)以100为终值的97的个数；

(C)以99为终值的48个数；(D)以100为终值的48个数。

3、theta=linspace(0,2*pi,100);r=cos(5*theta);polar(theta,r；k，)功能是()

(A)绘四叶玫瑰线；(B)绘五叶玫瑰线；(C)绘心脏线；(D)绘十叶玫瑰线。

4、MATLAB绘三维曲面需要构建网格数据，语句［x,y］=meshgrid(-2:2)返回数据中()

(A)x是行向量，y是列向量；(B)x是列向量，y是行向量；

(C)y是行元素相同的矩阵；(D)x是行向量相同的矩阵

5、十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪"，命令k=mod(2011-4,12)+1的结果是()

(A)k指向第二动物牛；(B)k指向第三动物虎；

(C)k指向第四动物兔；(D)k指向第五动物龙。

6、MATLAB命令symsx;f=sin(x);int(f*f,x,0,pi)的功能是()

(A)计算函数f在［0,2*pi］图形；(B)计算函数协2在［0,2*pi］的积分；

(C)计算旋转曲面所围的体积；(D)计算旋转曲面的表面积。

7、己知A是一个nXn的可逆矩阵，下面正确的说法是()

(A)rank(A)求A的行列式；(B)inv(A)求A的秩；

(C)sum(A)求A的列和；(D)［P,D］=eig(A)求A的特征值。

8、MATLAB计算正态分布随机变量分布函数的方法是()

(A)binocdf(x,n,p);(B)normpdf(x,mu,p);(C)binopdf(x,n,p);(D)normcdf(xzn,p)o

9、v0=515m/s;515*515*sin(2*alpha)"8的功能是计算关于抛射体问题的()

(A)发射角；(B)最大高度；(C)最大飞行时间；(D)最大射程。

10、绘制绕Z轴旋转的曲面方程是()

(A)/(%)=\］y2+z2(B)f(y)=y/z2+x2

22

(C)f(z)=Jy2+x。(D)z=yly+x

二、程序阅读理解（24分）

1.数学实验程序如下：

h=439;H=2384;R=6400;

a=(h+H+2*R)/2;c=(H-h)/2;

el=c/a;b=sqrt(a*a-c*c);

symse2t

f=sqrt(1-e2*cos(t)A2);

ft=subs(f,e2,el*el);

S=int(ft,0,pi/2);

L=4*a*double(S)；

V=L/(114*60)

sl=pi*a*b/(l14*60);%第十行

Vmax=2*sl/(h+R)

Vmin=2*sl/(H+R)

（1）实验程序的运行后，将显示的数据是（)

（A）卫星运行的近地速度、远地速度；(B)卫星运行的近地速度、远地速度和平均速度;

（C）卫星运行时向径每秒扫过的面积：(D)卫星运行的平均速度数据

（2）第十行语句的功能是（）

（A）计算卫星运行的最小速度；(B)计算卫星运行时向径每秒扫过的面积；

（C）计算卫星运行的最大速度；(D)计算卫星运行轨道的平均速度

2.数学实验程序如下：

xy=[-4-4;4-4;44;-44;-4-4];

A=[cos(pi/24)-sin(pi/24);sin(pi/24)cos(pi/24)];

x=xy(:,l);y=xy(:,2);

line(x,y),pause(l)

fork=l:30

xy=89*xy*A,;%第六行

x=xy（:,l）;y=xy（:,2）;

line（x,y）,pause（l）

end

（1）第六行语句中的矩阵A'的功能是（）

（A）正交矩阵；

（B）作用于矩阵xy时，其效果是将xy逆时针旋转为pi?24；

（C）作用于矩阵xy时，其效果是将xy旋转为pi>24,并将xy压缩为0.89；

（D作用于矩阵xy时，其效果是将xy顺时针旋转为pi/24。

（2）对下面有关程序的功能的说法确切的是（）

（A）边长为4以原点为中心的正方形旋转pi/24重复30次并绘图；

（B）边长为8以原点为中心的正方形旋转pi/24并将其边长压缩0.89,重复30次并绘图；（

（C）边长为4以原点为中心的正方形顺时针旋转pi/24并将其边长压缩0.89,重复30次并绘图；

（D）边长为8以原点为中心的正方形逆时针针旋转pi/24并将其边长压缩0.89,重复30次并绘图。

3、Viviani体是圆柱体（x-R/2）2+y2<R2/4被球面x2+y2+z2=R2所割立体。下面的数学实验程序功自

是取R=2求体积上半部分，先利用符号处理重积分并转换为数值数据，再用蒙特卡罗方法计算体做对比。

数学实验程序如下:

symsxy;

f=sqrt(4-xA2-yA2);

y1=sqrt(2*x-xA2);y2=sqrt(2*x-xA2);

Sl=int(f,y,y1,y2);S2=int(S1,x,0,2)

V=double(S2)

P=rand(10000,3)；

X=2*P(:1);Y=2*P(:,2);Z=2*P(:,3);

II=find((X-l).A2+Y.A2<=1&Z<=sqrt(4-X.A2-Y.A2));

Vl=8*length(II)/10000

(1)符号计算所用的积分公式是()

⑻丫=同：然）4一/723

(C)V=\+^^-x2-y2dy

J-yj2x-x~

（2）蒙特卡罗方法选用的随机点变化范围的立方体区域是（）

(A)Q={(x,y,z)|xe(0,2),yG(0,2),zG(0,2)}；

(B)。={(x,y,z)|%£(0,2),ye(-1,1),zG(0,2)}

(C)Q={(x,y,z)10<x<2,0<y<1,0<z<2}

(D)Q={(x,y,z)|—2<x<2,0<y<1,0<z<2}

三、程序填空（24分）

1、判断闰年条件有二

①能被4整除，但不能被100整除;

②能被4整除，又能被400整除。

year=input（'inputyear:=*）;

nl=year/4;

n2=year/100;

n3=year/400;

if_______________________________

dispC是闰年）

elseif___________________________②

dispC是闰年》

else

dispC不是闰年）

end

2.对于二重积分“dxdy,。是由两条抛物线y=&和抛物线y=V所围区域。下面数学实验程

D

序的功能是利用MATLAB符号计算方法将二重积分处理为二次积分进行计算（选择先对x积分后对y积

分的秩序），然后将计算结果的符号表达式转换为双精度实数。

symsxy;

f=x*sqrt(y);

xl=yA2

x2=一①

Sl=int(f,x,xl,x2);

S2=____________②

S=double（S2）

3.汽车紧急刹车问题数据拟合实验

V:2025303540455055606570

T:202841537293118149182221266

V表示刹车时汽车行驶速度(英里/小时),T表示刹车后汽车滑行距离(英尺)

分别做二次多项式和三次多项式数据拟合实验，并绘出数据拟合曲线的图，计算出残差平方和，完成如下

实验程序填空

v=[2025303540455055606570]*1.609;

T=[202841537293118149182221266]*.3048;

P2==polyfit(v,T,2)

T2=polyval(P2,v);

R2=sum((T-T2).A2)

figure(2),plot(v,T；*\v,T2)

P3=®

T3=polyval(P3,v);

R3=②

figure⑶,plot(v,T,'*',v,T3)

四、数学建模(12分)

每年A镇的人口10%迁往B镇;B镇的人口15%迁往A镇.假设某年A、B两镇人口各有110人和90人.假

设两镇人口流动正常的话，问两年后两镇人口数量分布如何?寻找方案使每年人口正常流动而A镇和B镇

的人口数量不增不减。

电子科技大学成都学院二零壹零至二零壹壹学年第一学期

数学实验课程考试题B(120分钟)闭卷考试时间：

成绩构成：平时(50%)+期末(50%)=总成绩(100%)

一二三合计

注意：请同学们将答案填写在答题纸上，否则无效。

一、单项选择题(共40分)

1、下面有关MATLAB变量名和函数名的说法中，错误的说法是()

(A)变量名的第一个字符必须是一个英文字母

(B)变量名可由英文字母、数字和下划线混合组成

(C)变量名不得包含空格和标点，但可以有下连字符

(D)变量名和函数名对于英文的大小使用没有区别

2、符号计算与一般数值计算有很大区别，它将得到准确的符号表达式。在MATLAB命令窗口中键入命令

symsxyl=xA2；y2=sqrt(x)；int(yl-y2,x,0,l)屏幕显示的结果是()

(A)yl=xA(l/2)(B)ans=l/3；(C)y2=xA2；(D)ans=-l/3

3、MATLAB命令symsx;f=sin(x);V=pi*int(f*f,x,O,pi)功能是()

(A)绘出函数f在[0,2万]图形；(B)计算函数f在[0,2%]的积分；

(C)计算旋转曲面所围的体积；(D)计算旋转曲面的表面积。

4、用赋值语句给定x数据，计算J；sin(3+2x)+e31n3对应的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)+eA3*log(3))(D)sqr(7sin(3+2x)+eA31og(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指向第五动物龙。

7、MATLAB命令A=magic(3)创建3阶幻方矩阵，求A的特征值绝对值最小用()

(A)min(abs(eig(A)));(B)min(eig(abs(A)));(C)min(eig(A));(D)min(abs(A));

8、MATLAB命令x=rand(10,l)生成10个随机数，将它们从大到小排序，使用()

(A)y=sort(x);z=y(10:l);(B)[y,II]=sort(x);z=y(II);

(C)y=sort(x);z=y(10:-1;1);(D)[y,II]=sort(x);z=x(H);

9、MATLAB命令[x,y]=meshgrid(l:3);H=l./(x+y-l)产生的矩阵H是()

-iir-123--123-■11/21/3

(A)222(B)123(C)233(D)1/21/31/4

3331233451/31/41/5

10、设a,b,c表示三角形的三条边，表达式a+b>c|a+c>b|b+c>a,()

(A)是三条边构成三角形的条件；(B)是三条边不构成三角形的条件;

(C)构成三角形时逻辑值为真；(D)不构成三角形时逻辑值为假

二、程序阅读理解(24分)

1、用plot命令绘制平面曲线

v0=515;alpha=45*pi/l80;

T=2*v0*sin(alpha)/9.8;%第二行

t=T*(0:16)/16;

x=v0*t*cos(alpha);

y=v0*sin(alpha)・0.5*9.8*t.八2;

plot(x,y,x,y,,r*>)

Xmax=x(17)

8000

6000

4000

2000

00.511.522.53

(1)对下面有关程序的功能的说法确切的是()

(A)以515为初速度的抛射曲线的绘制；

(B)以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:l;x=r'*cos(t);y=r'*sin(t);

z=sqrt(l-x.A2);%第三行

meshz(x,y,z),axisoff

colortnap([001])

view(-47,56),holdon

x1=cos(t);y1=sin(t);z1=abs(sin(t));

plot3(xl,yl,zl,'ro')；

(1)下面有关程序的功能的说法确切的是()

(A)绘圆柱面x2+y2=l,x2+z2=1的交线；

(B)绘圆柱面x2+y=i,x2+z2=1所围区域的边界曲面；

(C)绘圆柱面N+y2=i,%2+z2=1的交线及所围区域的边界曲面；

(D)绘圆柱面x2+y2=i,%2+z2=1的交线及所围区域的边界曲面的上半部分。

(2)关于第三行语句错误的解释是()

(A)z是矩形域上曲顶柱面高度值；(B)z是与y同型的矩阵；

(C)z是圆域上曲顶柱面高度值；(D)z是与x同型的矩阵

3、非负函数y=f(x)在有限区间上的图形为上半平面的一条曲线，曲线绕x轴旋转时，产生以x为对称轴

的旋转曲面，其体积

symsabx

f=exp(a*x)*sin(b*x);

fl=subs(f,a,-0.2);

f2=subs(fl,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=tF*ones(size(t));%第九行

y=r**cos(theta);%第十行

z=r'*sin(lhela);%第^^一行

mesh(z,y,x)

colormap([000])

axisoff

view(-17,54)

(5)关于程序的功能确切的说法()

(A)计算曲线段f(x)=exp(ax)sin(bx),()《xw2乃绕X