k1=3*(k-1)+1;k2=3*k;
Sk=S(k1:k2)
Show=strcat('你选择了----',Sk)
N=2000;R=1+fix(6*rand(4,N));
x=sum(R);y=mod(x,12)+1;
II=find(y==k);
%第十行语句
n=length(II);Fn=n/N
(1) 当用户调用函数程序时，没有输入数据，则程序运行后将显示 2000 次随机实验( B )
有限位实数表示 3 2 ，就会犯下错误。数学实验程序验证了这个事实，程序运行后误差如右
文本框所示
a=2^(1/3); D=1; for k=1:8
D=D*10; b=fix(a*D)/D; V(k)=b^3; end error=V’-2
error=-2.7200e-001 -4.6875e-002 -4.3830e-003 -1.0024e-004 -4.9998e-006 -2.3761e-007 -2.3761e-007 -4.7121e-008
V=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 )
5
å 1．MATLAB 命令 A=rand(5,5)；创建 A
=
(aij)5´ 5 ，求 m
ax
j
|aij |用( A )
i= 1
(A) max(sum(abs(A))); (B) max(sum(abs(A’))); (C) max(sum(A))); (D) sum(max(A));
2．MATLAB 命令 x=[1,2,4,5,9];mean(x)，的计算结果是( B )
图 1 Vivinai 问题
2
+ 2x- x2
2
1+ 1- y2
蝌 蝌 (A) V = dx
4 - x 2 - y 2dy (B) V = dy
4 - x 2 - y 2dx
0
- 2x- x2
0
1- 1- y2
+ 2x- x2
ò (C) V =
4 - x 2 - y 2dy
- 2x- x2
+ 2x- x2
7．MATLAB 命令 syms x; f=sin(x); V=pi*int(f*f,x,0,pi)功能是( C )
(A) 绘出函数 f 在[0,2 p ]图形；
(B) 计算函数 f 在[0,2 p ]的积分；
(C) 计算旋转曲面所围的体积；
(D) 计算旋转曲面的表面积。
8．十二属相为“鼠牛虎兔龙蛇马羊猴鸡狗猪”，命令 k=mod(2008,12)+1 的结果是( D )
(A) k 指向第二动物牛；
(B) k 指向第三动物虎；
(C) k 指向第四动物兔；
(D) k 指向第五动物龙。
9．MATLAB 命令[x,y]=meshgrid(1:3);H=1./(x+y-1)产生的矩阵 H 是( D )
轾犏1 1 1 (A) 犏犏2 2 2
犏犏臌3 3 3
轾犏1 2 3 (B) 犏犏1 2 3
数学实验程序计算。
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 )
(A) 目标函数系数； (B) 等式约束系数； (C) 不等式约束系数； (D) 等式约束常向量
填空：
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= 0.9*xy*A' ①;
(A) S n =
2
k=
D
1
k
；(B)
Sn
=
D k ；(C) S n =
k= 1
2
k=
D
2
k
；(D)
Sn
=
Dk
k= 2
三、程序填空(10 分)
1．二阶正交矩阵作用于某一向量时，其效果是将该向量旋转，旋转解为 a (逆时针旋转为正)。
把一个以原点为中心的正三角形旋转 p / 50 ，并缩小 90%，迭代 33 次创建图 3。完成程序
下面的数学实验程序功能是取 R=2 求体积上半部分，先利用符号处理重积分并转换为数值
数据，再用蒙特卡罗方法计算体做对比。数学实验程序如下：
syms x y;
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)
ò (D) V = - 2x- x2
(2) 蒙特卡罗方法选用的随机点变化范围的立方体区域是( B )
4 - x 2 - y 2dx
(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)}
)
。则多边形面积可由二阶行列式求和计算，数
学实验程序如下：
data=[-1,-1;1,-1;1,1;0,0;-1,1]; n=size(data,1);Sk=0;pk1=data(1,:); for k=2:n
pk=data(k,:);Dk=det([pk1;pk]); Sk=Sk+Dk;pk1=pk; end pk=data(1,:);Dk=det([pk1;pk]); Sk=Sk+Dk; Sn=0.5*Sk (1) 程序中所用的二阶行列式是( C )
(2) 程序中变量 b 存放的数据是( D )
(A) 将 a 的小数点后第 k 位减 1 所得； (B) 将 a 的小数点 k 位后按四舍五入所得；
(C) 将 a 的小数点后第 k 位增 1 所得； (D) 将 a 的小数点 k 位后截断舍去所得。
2．Viviani 体是圆柱体(x - R / 2)2 + y 2 £ R 2 / 4 被球面 x2 + y2 + z2 = R 2 所割立体。
(B) 变量名可由英文字母、数字和下划线混合组成
(C) 变量名不得包含空格和标点，但可以有下连字符
(D) 变量名和函数名对于英文的大小使用没有区别
二、程序阅读题 (40 分)
1．传说古希腊曾流行瘟疫，人们为消除灾难求助于神。神说：把神庙中黄金祭台增容一倍，
可消除瘟疫。当立方体祭台尺寸放大一倍后，瘟疫仍然流行。人们才知道体积并不是扩大了 两倍。这个古希腊难题被称为倍立方体问题，在人类还没有认识到无理数时，企业界企图用
(A) 4
(B) 4.2 (B) 4.5 (D) 21
3．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);
x=xy(:,1);
y=
xy(:,2)
②;
line(x,y)
1 0.8 0.6 0.4 0.2
0 -0.2 -0.4 -0.6
-1 -0.8 -0.6 -0.4 -0.2 0 0.2 0.4 0.6 0.8 1
图 3 旋转三角形
end
2．长征三号甲运载火箭提供给探月卫星的初始速度不足以将卫星送往月球轨道。为提高到
(C) 统计 2000 次随机实验中，游戏者可能是第 k 个星座的次数；
(D) 统计 2000 次随机实验中，游戏者可能是第 k 个星座的频数。
5．一个平面多边形由它的
n
个顶点确定，将顶点按逆时针方向排列为：
பைடு நூலகம்
P
(x ,y
111
)
，...，
P
n
(xn
,yn
)。将第(n+1)个顶点设为
P
(x ,y
111
犏犏臌1 2 3
轾犏1 2 3 (C) 犏犏2 3 3
犏犏臌3 4 5
轾犏1 1 / 2 1 / 3 (D) 犏犏1 / 2 1 / 3 1 / 4
犏犏臌1 / 3 1 / 4 1 / 5
10．下面有关 MATLAB 变量名和函数名的说法中，错误的说法是( D )
(A) 变量名的第一个字符必须是一个英文字母
(1) 程序中循环控制变量 k 从 1 变量 8，而变量 D=10k 的作用是( C )
(A) 将 a 的小数点向右移 D 位取整；
(B) 将 a 的小数点向右移 D 位取整后再向左移 D 位；
(C) 将 a 的小数点向右移 k 位取整后再向左移 k 位；
(D) 将 a 的小数点向左移 k 位取整后再向右移 k 位；
(A) 游戏者可能是白羊座的频率； (B) 游戏者可能是金牛座的频率；
(C) 游戏者可能是双子座的频率； (D) 游戏者可能是其它星座的频率。