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

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

>>A=size([5 5 6;6 6 6;7 7 8]); >>B=size([1 2;3 4;5 6]); >>C=A ’*B
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.用round 函数四舍五入对数组[2.48 6.39 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, 6]
B. array = [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]
C. array =6
D. array =10.
10.i=2; a=2i; b=2*i; 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
二、填空题（每空2分，共38分）
1. 已知A=[0 9 6;1 3 0];B=[1 4 3;1 5 0];写出下列各指令运行的结果。

A &
B 的运行结果ans= ；
A >
B 的运行结果ans= 。

2.产生4阶全0方阵的命令为 ________ ；产生3阶全1方阵的命令为
_____________ 。

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

4. 标点符号 _____ 可以使命令行不显示运算结果， _____ 用来表示该
行为注释行。

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

6.为了使两个plot的图形在同一个坐标显示，可以使用 ______ 命令
进行图形保持；可以使用 _________ 命令为图形添加网格。

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

8.在matlab中，执行命令dms2mat(rad2dms（pi）)，结果为_________
______________________________；函数dms2rad(mat2dms( )) 的功能是___________________________________________。

9.在语句末添加逗号或分号，区别是_______________________
_________________________________________________。

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

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

12.在当前路径下有一个数据文件my.txt,若要打开此文件并进行读
取数据，则执行命令为fid=_______________________________。

三、简短程序阅读题（每题6分，共24分）
1、请修改下面的程序，让程序中没有for循环语句也能完成任务。

（本题6分）
A=[1 2 3 4;5 6 7 8;9 10 11 12];
[r c]=size(A);
for i=1:r
for j=1:c
if (A(i,j)>9 | A(i,j)<4)
A(i,j)=0;
end
end
end
disp(A)
修改为：
２
2.写出执行以下代码后C，D，E的值 (6分)
A=[1,2,3;4:6;7:9];
C=[A;[10,11,12]],
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;
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:2π]均匀的取50个点，构成向量t
(2)在同一窗口绘制曲线y1=sin(2*t-0.3); y2=3cos(t+0.5)；要求
y1曲线为红色点划线，标记点为圆圈；y2为蓝色虚线，标记点为星号。

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

（本题8分）
function [Lv,V,d0]=szwtjpc(L,A,W,P) %______________________ nl=length(L); %_________________________________________ [ma, na]=size(A);
%___________________________________________________ nw=length(W); if nl~=na
disp(' 观测值个数与改正数个数不相等! ') return
end
if ma~=nw
warning('所列条件方程有误') return end
[mp,np]=size(P); if mp~=np
warning('权阵应为方阵') return end
Q=inv(P);
%_______________________________________________ W=(A*L+W);
%_______________________________________________ V=-Q*A'*inv(A*Q*A')*W
%_______________________________________________ Lv=L+V/1000
%_______________________________________________ d0=sqrt(V'*P*V/nw)
%_______________________________________________
A C
E。