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

１
淮 海 大 学
2014 - 2015 学年第一学期 测量程序设计基础 试卷（B 闭卷）
1.运行以下命令后，变量C 的值是________。

>>A=size([1 2 3;4 5 6;7 8 9]); >>B=size([1 2;3 4;5 6]); >>C=B ’*A
A. ⎪⎪⎭⎫ ⎝⎛6699
B. ⎪⎪⎭⎫ ⎝⎛6969
C. ⎪⎪⎭⎫ ⎝⎛9669
D. ⎪⎪⎭⎫
⎝⎛6996 2.matlab 的循环结构中，如果想让程序运行调出循环，
则可以使用函数______。

A. break
B. tic/toc
C.time
D. pause
3. 下列哪个变量的定义是合法的
A. abcd-3
B. Xyz?3
C. _abcdef
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.用fix 函数对数组[2.58 6.37 3.93 8.52]取整， 结果为
A. [2 6 3 8]
B. [2 6 4 8]
C. [2 6 4 9]
D. [3 7 4 9]
8.已知a=2:2:8, b=2:5，下面的运算表达式中，错误的为 (A) a/b (B) a .*b (C) a\b (D) a*b 9.下面的程序执行后array 的值为_________ for k=1:10 if k>=6 break; 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]
C. array =5
D. array =10.
10.j=2; a=2j; b=2*j; c=2*sqrt(-1); 程序执行后；a, b, c 的 值分别是__________________ A.a=4, b=4, c=2.0000i B.a=2.0000i, b=4, c=2.0000i C.a=4, b=2.0000i, c=2.0000i D. a=2.0000i, b=2.0000i, c=2.0000i
二、填空题（每空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阶全0方阵的命令为 ________ ；产生4阶全1方阵的命令为
_____________ 。

3.A=randn(2,5);b=size(A);c=length(A)；则b和c的值分别为
_____________和 ________ 。

4. 标点符号 _____ 可以使命令行显示运算结果， ________ 用来表示
该程序段为函数式M文件。

5. 用if判断语句判断90> x >50,在MATLAB中if语言后的判断应写为
________________________ 。

6.为了使后面用plot绘制的图形不和前面的图形在同一个坐标系内，可
以使用 ______ 命令；可以使用 _________ 命令取消图形的网格。

7.MATLAB的工作空间中有三个变量v1, v2, v3，写出把v1, v2保存到
文件my_data.mat中的指令 _______ ；写出把my_data.mat文件中的变量v1读取到MATLAB工作空间内的指令 _________ 。

8.在matlab中，函数dms2rad(mat2dms( ))的功能为_________
______________________________；函数dms2mat(rad2dms（）) 的功能是___________________________________________。

9.在语句末添加分号的作用是___________________________________。

10. 在数字测量技术中，通常要绘制等高线图，这时需要调用函数为
______________________。

11.在matlab中，函数fopen（）的功能是_________________
_________________________________________。

12.在当前路径下有一个文件Ex0805.txt,若在MATLAB环境中显示
的文件内容 ,则需要执行命令为:_______________________________。

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

（本题6分）
A=[1 2 3 4;5 6 7 8;9 10 11 12;13 14 15 16];
[r c]=size(A);
for i=1:r
for j=1:c
if (A(i,j)>13 | A(i,j)<9)
A(i,j)=0;
end
end
end
disp(A)
结果为(必须用矩阵表示):
２
2.写出执行以下代码后C，D，E的值 (6分)
A=[1,2,3;4:6;7:9;[10 11 12]];
C=[A;[13,14,15]],
D=C(1:3,[2 3])
E=C(2,[1 2])
3.写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值 (6分) x=[0,1,0,2,0,3,0,4];
for k=1:8
if x(k)==0
x(k)=k^2;
else
x(k)=2*k-1;
end
end
disp(x);
4.求解以下线性方程组，只要求写出程序代码（6分）
2x1- 3x2+ x3+2x4=8
x1+3x2+ x4=6
x1- x2+ x3+8x4=1
7x1+ x2-2x3+2x4=5
四、程序设计题（本题满分18分）
1、绘制函数曲线，只要求写出程序代码(10分)
(1)在区间[0:4π]均匀的取100个点，构成向量t
(2)在同一窗口绘制曲线y1=cos(2*t-0.5); y2=3sin(t+0.3)；要求
y1曲线为蓝色点划线，标记点为圆圈；y2为绿色虚线，标记点为星号。

３
2、A和B是已知高程的水准点，C、D、E是待定点，A、B高程、观测高差和相应的水准路线长度见下表，为求各待定点的高程平差值，利用条件平差进行编程，请根据题意和程序中上、下语句的关系，注释相应语句的意思。

（本题8分）
function jiesuan2( )
disp('水准网间接平差示例')
disp('已知高程')
Ha = 5.015
Hb = 6.016
disp('观测高差')
L = [1.359;2.009; 0.363; 1.012; 0.657; -0.357]
f id=fopen(‘E:\szjj.txt’,‘r’) %____________________________________ temp=fscanf(fid,‘%f%f%f%f’,*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))%___________________________
A
E
４。