2014-2015(2)测量程序设计基础B

淮 海 工 学 院

2014 - 2015 学年第二学期 测量程序设计基础模拟 试卷B(闭卷）

一 二 三 四 总分 核分人 题号 分值 20 38 18 24 100 得分 一、选择题（每小题2分，共20分） 运行以下命令后，变量C的值是________。

>>A=size([5 5 5;6 6 6;7 7 7]); >>B=size([1 2;3 4;5 6]); >>C=A’*B

A. ??99???96??9?66?? B. ????96?? C. ??6???69?? D. ???69???96?? ?如果想使matlab程序段暂行些时间，可以使用函数______。 A. break B. tic/toc C.time D. pause

3. 下列哪个函数式M文件的定义是不合法的

A. abcde3 B. xyzw_3 C. abcdef-3 D. x3yz 4.下列哪条指令是求方阵行列式的命令 A. inv B. diag C. det D. eig

5.在循环结构中跳出循环，执行循环后的命令，则使用 A. return B. break C. continue D. keyboard 6.清空Matlab工作空间内所有变量的命令是 A. clc B. cls C. clear D. clf

7.用floor函数可以对数组[2.48 6.39 3.93 -8.52]取整， 结果为

A. [2 6 3 -8] B. [2 6 3 -9] C. [2 7 4 -9] D. [3 7 4 -8]

8.已知a=2:2:11, b=2:6，下面的运算表达式中，出错的为 A. a'*b B. a .*b C. b*a D. a-b 9.下面的程序执行后array的值为_________ for k=1:11 if k>5 return; else

array(k) = k; end end

A. array = [1, 2, 3, 4, 5] B. array = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10,11] C. array =5 D. array =11.

10.j=2; a=2*j; b=2j; c=2*sqrt(-1); 程序执行后；a, b, c的 值分别是__________________ A.a=4, b=4, c=2.0000i B.a=4, b=2.0000i, c=2.0000i C.a=2.0000i, b=4, c=2.0000i D. a=2.0000i, b=2.0000i, c=2.0000i

１ 1. 2.

二、填空题（每空2分，共38分）

1. 已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。

A | B 的运行结果ans=[__________;___________]； A == B 的运行结果ans=[__________;___________]。

2. 产生5阶单位阵的命令为 _________ ；产生5阶全0方阵的命令为 _______________ 。

3. A=eye(2,6);b=size(A);c=length(A)；则b和c的值分别为[2 ,___]

和 __。

4. 若命令后面带有逗号，则表示该命令行执行后___________，

符号 _____ 用来表示该行为注释行。

5. 用if判断语句判断80＞ x ≥40,在MATLAB中if语言后的判断应写为 ___________ 。

6.为了使两个plot的图形不同时显示，可以使用 _______ 命令清除第一个图形；可以使用 ___________ 命令为图形取消网格线。 7.MATLAB的工作空间中有三个变量a, b, c，写出把变量b保存到文件myfile.mat中的指令 ______________ ；写出把myfile.mat文件中的变量b读取到MATLAB工作空间内的指令 _______________ 。 8.在matlab中，执行命令dms2mat(rad2dms（2*pi）)，结果为[________,0 ,______]。

9.识别是函数式M文件还是命令式M文件，就看_____________。 10. 在数字测量技术中，通常要绘制等高线图，这时需要调用函数为 ________________。

11.在matlab中，函数fid=fopen（filename,permission）,若fid大于零，则表明_______________。

12.执行命令st=fclose(fid)后，若st=-1,则表明_________________ ______________。

三、简短程序阅读题（每题6分，共18分） 1、阅读下面的程序，写出程序运行的结果。

（本题6分）

A=[1 2 3 4 5;6 7 8 9 10; 11 12 13 14 15]; [r c]=size(A); for i=1:r for j=1:c

if (A(i,j)>8 | A(i,j)<5) A(i,j)=0; end end end disp(A)

结果为：（要求写成矩阵形式）

2. 写出执行以下代码后C，D，E的值(写成矩阵形式) (6分)

A=[1,2,3 4;5:8;9:2:16]; C=[A;17:20], D=C(1:2,[2 3]) E=D(2,[1 2])

２C = D = E =

3.写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 (6分) x=[0,1,0,2,0,3,0,4，0,5]; for k=1:10 end disp(x)

运行此程序段，结果为（写成矩阵的形式）： 四、程序设计题（本题满分24分）

1.求解以下线性方程组，要求写出全部的7行程序代码（7分）

2x1- 3x2+ x3+2x4=8 x1+3x2+ x4=6 x1- x2+ x3+8x4=17x1+ x2-2x3+2x4=5

2．绘制函数曲线，要求写出程序代码(8分) (1)在区间[0,3π]均匀的取200个点，构成向量t

(2)在同一窗口绘制曲线y1=4cos(2*t+0.5); y2=3sin(t-0.3)；要求 y1曲线为绿色点划线，标记点为圆圈；y2为红色虚线，标记点为星号。(1)

(2)

if x(k)==0 x(k)=k; else

x(k)=k^2+1;

３ end

3.A和B是已知高程的水准点，C、D、E是待定点，A、B高程、观测高差和相应的水准路线长度见下表，为求各高差的平差值，利用间接平差进行编程，请根据题意和程序中上、下语句的关系，注释相应语句的意思。（本题9分） C h1h6 E h3 Ah5h7

h2B Dh4

线路编号 观测高差（m） 线路长度（km） 已知高程（m） 1 +1.360 2.2 HA=5.015 2 +2.010 2.6 HB=6.016 3 +0.364 3.3 4 +1.013 3.5 5 +0.658 3.4 6 -0.357 6.1

function jiesuan( )

disp('水准网间接平差示例') disp('已知高程')

Ha = 5.015 % 已知点高程 Hb = 6.016 % 已知点高程

disp('观测高差') %____________________________ L = [1.360;2.010; 0.364; 1.013; 0.658; -0.357] ;

fid=fopen(‘E:\\szjj.txt’,’r’) %______________________

temp=fscanf(fid,’%f%f%f%f’ ,[4,6])

%上面命令的意思是将文件中的数据读入一个4行6列的矩阵 fclose(fid); %________________________________ temp=temp’

B=temp(:,[1,2]); %______________________________ l=temp(:,3);

S=temp(:,4); %_________________________________ S=S’;

disp(‘C是单位线路公里数，S是线路长度') C = ones(1,6);

P = C./S ; %_______________________________ P = diag(P) % ______________________________ disp(‘参数的解’)

x=inv(B'*P*B)*B'*P*l;

disp('误差V(mm), 各待定点的高程平差值L1（m）') V=B*x-l; % ______________________________ L1 = L + V/1000 disp('精度评定')

n = 6; % 观测值的个数 t = 2; % 必要观测数

delta = sqrt(V'*P*V/(n-t)) %_________________________

４