2014-2015测量程序设计基础B答案

１淮海工学院2014-2015学年第二学期测量程序设计基础模拟试卷B(闭卷）参考答案与评分标准题号一二三四总分核分人分值20381824100得分一、选择题（每小题2分，共20分）1.运行以下命令后，变量C的值是____B____。A=size([555;666;777]);B=size([12;34;56]);C=A’*BA.6699B.6969C.9669D.69962.如果想使matlab程序段暂行些时间，可以使用函数____D__。A.breakB.tic/tocC.timeD.pause3.下列哪个函数式M文件的定义是不合法的CA.abcde3B.xyzw_3C.abcdef-3D.x3yz4.下列哪条指令是求方阵行列式的命令CA.invB.diagC.detD.eig5.在循环结构中跳出循环，执行循环后的命令，则使用BA.returnB.breakC.continueD.keyboard6.清空Matlab工作空间内所有变量的命令是CA.clcB.clsC.clearD.clf7.用floor函数可以对数组[2.486.393.93-8.52]取整，结果为BA.[263-8]B.[263-9]C.[274-9]D.[374-8]8.已知a=2:2:11,b=2:6，下面的运算表达式中，出错的为CA.a'*bB.a.*bC.b*aD.a-b9.下面的程序执行后array的值为___A______fork=1:11ifk5return;elsearray(k)=k;endendA.array=[1,2,3,4,5]B.array=[1,2,3,4,5,6,7,8,9,10,11]C.array=5D.array=11.10.j=2;a=2*j;b=2j;c=2*sqrt(-1);程序执行后；a,b,c的值分别是________B__________A.a=4,b=4,c=2.0000iB.a=4,b=2.0000i,c=2.0000iC.a=2.0000i,b=4,c=2.0000iD.a=2.0000i,b=2.0000i,c=2.0000i２二、填空题（每空2分，共38分）1.已知A=[096;130];B=[143;150];写出下列各指令运行的结果。A|B的运行结果ans=[111;110]；A==B的运行结果ans=[000;101]。2.产生5阶单位阵的命令为_eye(5)；产生5阶全0方阵的命令为__zeros(5)。3.A=eye(2,6);b=size(A);c=length(A)；则b和c的值分别为[2,_6_]和6__。4.若命令后面带有逗号，则表示该命令行执行后__显示结果_，符号__%___用来表示该行为注释行。5.用if判断语句判断80＞x≥40,在MATLAB中if语言后的判断应写为_80x&x=40。6.为了使两个plot的图形不同时显示，可以使用holdoff命令清除第一个图形；可以使用_gridoff命令为图形取消网格线。7.MATLAB的工作空间中有三个变量a,b,c，写出把变量b保存到文件myfile.mat中的指令savemyfileb；写出把myfile.mat文件中的变量b读取到MATLAB工作空间内的指令loadmyfileb。8.在matlab中，执行命令dms2mat(rad2dms（2*pi）)，结果为[360____,0,_0_]。9.识别是函数式M文件还是命令式M文件，就看_有无函数定义行_。10.在数字测量技术中，通常要绘制等高线图，这时需要调用函数为__contour_。11.在matlab中，函数fid=fopen（filename,permission）,若fid大于零，则表明_文件已被打开_。12.执行命令st=fclose(fid)后，若st=-1,则表明___________________文件没有被关闭_。三、简短程序阅读题（每题6分，共18分）1、阅读下面的程序，写出程序运行的结果。（本题6分）A=[12345;678910;1112131415];[rc]=size(A);fori=1:rforj=1:cif(A(i,j)8|A(i,j)5)A(i,j)=0;endendenddisp(A)结果为：（要求写成矩阵形式）000056780000000完全写对得6分，否则酌情给分2.写出执行以下代码后C，D，E的值(写成矩阵形式)(6分)A=[1,2,34;5:8;9:2:16];C=[A;17:20],D=C(1:2,[23])E=D(2,[12])３C=123456789111315171819202分D=23674分E=676分3.写出执行以下代码后，MATLAB命令窗口上显示的x矩阵的值(6分)x=[0,1,0,2,0,3,0,4，0,5];fork=1:10ifx(k)==0x(k)=k;elsex(k)=k^2+1;endenddisp(x)运行此程序段，结果为（写成矩阵的形式，否则不得分）：153175377659101完全写对得6分，否则酌情给分四、程序设计题（本题满分24分）1.求解以下线性方程组，要求写出全部7行程序代码（7分）2x1-3x2+x3+2x4=8x1+3x2+x4=6x1-x2+x3+8x4=17x1+x2-2x3+2x4=5A=[2-312;1301;1-118;71-22];b=[8;6;1;5];x=A\bx1=x(1)x2=x(2)x3=x(3)x4=x(4)每条命令写对得1分2．绘制函数曲线，要求写出程序代码(8分)(1)在区间[0,3π]均匀的取200个点，构成向量t(2)在同一窗口绘制曲线y1=4cos(2*t+0.5);y2=3sin(t-0.3)；要求y1曲线为绿色点划线，标记点为圆圈；y2为红色虚线，标记点为星号。(1)t=linspace(0,3*pi,200);2分(2)y1=4*cos(2*t+0.5);4分y2=3*sin(t-0.3);6分plot(t,y1,’g-.o’,t,y2,’r--*’)8分４3.A和B是已知高程的水准点，C、D、E是待定点，A、B高程、观测高差和相应的水准路线长度见下表，为求各高差的平差值，利用间接平差进行编程，请根据题意和程序中上、下语句的关系，注释相应语句的意思。（本题9分）functionjiesuan()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]);%_从6行4列的矩阵中取第1和第2列_l=temp(:,3);S=temp(:,4);%_矩阵的第4列为6条水准路线的长度_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/1000disp('精度评定')n=6;%观测值的个数t=2;%必要观测数delta=sqrt(V'*P*V/(n-t))%_计算单位权中误差的估值_每答对一条注释得1分。线路编号观测高差（m）线路长度（km）已知高程（m）123456+1.360+2.010+0.364+1.013+0.658-0.3572.22.63.33.53.46.1HA=5.015HB=6.016ACEBDh1h6h7h3h5h2h4