电子科技大学《数学实验》2008-2009学年期末试题(含答案)

电子科技大学二零零八到二零零九学年第二学期期末考试《数学实验》课程考试题A卷(120分钟) 考试形式：闭卷考试日期：2009年7月8日

一、单项选择题(20分)

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

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

(C) diag(diag(magic(3))); (D) diag(diag(magic))。

2、MATLAB命令P=pascal(3)将创建三阶帕斯卡矩阵，max(P)的计算结果是( )

(A) 1 2 3 (B) 1 2 1 (C) 3 6 10 (D) 1 3 6

3、命令J=*1;1;1+**1,2,3+;A=j+j’-1将创建矩阵( )

(A)

123

234

345

⎡⎤

⎢⎥

⎢⎥

⎢⎥

⎣⎦

; (B)

234

345

456

⎡⎤

⎢⎥

⎢⎥

⎢⎥

⎣⎦

(C)

123

123

123

⎡⎤

⎢⎥

⎢⎥

⎢⎥

⎣⎦

(D)

111

222

333

⎡⎤

⎢⎥

⎢⎥

⎢⎥

⎣⎦

4、data=rand(1000,2);x=data(:,1);y=data(:,2);II=find(yx.^2);的功能是( )

(A) 统计2000个随机点中落入特殊区域的点的索引值；

(B) 统计1000个随机点落入特殊区域的点的索引值；

(C) 模拟2000个随机点落入特殊区域的过程；

(D) 模拟1000个随机点落入特殊区域的过程。

5、MATLAB计算二项分布随机变量分布律的方法是( )

(A) binocdf(x,n,p); (B) normpdf(x,mu,s); (C)binopdf(x,n,p); (D) binornd(x,n,p)。

6、MATLAB命令syms e2;f=sqrt(1-e2*cos(t)^2);S=int(f,t,0,pi/2)功能是()

(A) 计算f(x)在[0,pi/2]上的积分；(B) 计算f(t)不定积分符号结果；

(C) 计算f(x)积分的数值结果；(D) 计算f(t)定积分的符号结果。

7、y=dsolve(‘Dy=1/(1+x^2)-2*y^2’,’y(0)=0’,’x’);ezplot(y)的功能是( )

(A) 求微分方程特解并绘图；(B) 解代数方程(C) 求定积分；(D)求微分方程通解。

8、X=10000 ;0.5*asin(9.8*X/(515^2))的功能是计算关于抛射体问题的()

(A) 十公里发射角；(B) 十公里飞行时间；(C)最大飞行时间；(D)最大射程。

9、theta=linspace(0,2*pi,100) ;r=cos(4*theta) ;polar(theta,r,’k’)功能是()

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

10、北京和纽约的经度分别是：东经118和西经76，根据经度差计算时差用()

(A) fai1=118;fai2=-76;Dfai=(fai1+fai2)/24; (B) fai1=118;fai2=-76;Dfai=(fai1+fai2)/15;

(C) fai1=118;fai2=-76;Dfai=(fai1-fai2)/24; (D) fai1=118;fai2=-76;Dfai=(fai1-fai2)/15。

二、程序阅读题(40分)

1、直方图功能是将数据分为n个类，统计各个类的数据量并绘图。借用现有的直方图命令hist，编写新直方图程序如下。

function m=myhist(data,n)

if nargin==1,n=7;end

Xmin=min(data);Xmax=max(data);h=(Xmax-Xmin)/n;

m=hist(data,n)/length(data)/h;

t=linspace(Xmin,Xmax,n+1);

II=1:4:4*n-3;JJ=1:n;

x(II)=t(JJ);y(II)=zeros(1,n);

x(II+1)=t(JJ);y(II+1)=m;

x(II+2)=t(JJ+1);y(II+2)=m;

x(II+3)=t(JJ+1);y(II+3)=zeros(1,n);

plot(x,y,'k')

(1) 变量data存放了1000个数据，在命令窗口调用myhist(data)的结果是( )

(A) 只绘数据的直方图而不显示被分类后各类的数据量；

(B) 只显示被分类后各类的数据量而不绘数据的直方图；

(C) 既绘数据直方图也显示被分类后各类的数据量；

(D) 根据默认值在数据范围内插入七等分点绘直方图。

(2) 关于新直方图绘图程序下面说法不正确的是( )

(A)h是n等分直方图中小区间长度；(B) 修改程序最后一行可绘红色直方图；

(C) 直方图中所有小矩形面积之和为1；(D) 直方图中所有小矩形的高度和为1。

2、3n+1问题反映一个数学猜想：对任一自然数n，按如下法则进行运算：若n为偶数，则将n除2，若n为奇数，则将n乘3加1。重复这种操作，结果终会为1。实验程序如下。function [k,N]=threeN(n)

if nargin==0,n=5;end

k=1;N=n;

while n~=1

r=rem(n,2);

if r==0

n=n/2;

else

n=3*n+1;

end

N=[N,n];k=k+1;

end

(1)在MATLAB命令窗口中直接调用threeN运行结果为( ) %5 16 8 4 2 1

(A)只显示k的最后数值为6；(B) 只显示k的最后数值5；

(C) 同时显示k和N的数据；(D) 仅显示N的所有数据。

(2)实验程序运行过程中( )

(A) 输入变量n不发生改变；(B)N是记录数据变化的一维数组；

(C) N记录每次数据变化的单个数据；(D)n是记录数据变化的一维数组。

3、将半径为r的球体(密度1

ρ)置入水中，球体将浮出水面一定高度。程序如下：

<

function [h,Rou]=highNu(r)

if nargin==0,r=10;end

Rou=0.3:0.1:1;

N=length(Rou);

for k=1:N

rouk=Rou(k);

P=[1,-3*r,0,4*r^3*rouk];