（2）、创建word对象Set wordapp = CreateObject("Word.Application")（2分）
（3）、设置对象可见wordapp.Visible = True（1分）
（4）、释放对象Set worddoc = Nothing（1分）
得分
阅卷人
三、编程题（每题7分，共14分）
Private Sub F1_Click()
Dim dM As String, i%
Dim n%, s%, x1#, y1#, h1, h2, gy
Dim A() As Double, L() As Double, cs(1 To 6) As Double
Dim At#(), Naa#(), W#()
山东交通学院期末考试测量程序设计课程试卷 （A）卷2010——2011 学年第 1 学期 第 1 页 共 3 页
题号
一
二
三
四
五
六
七
总分
审核
得分
得分
阅卷人
一、填空题（每空2分，共40分）
1、VB采用面向对象的程序设计，对象是系统基本运行实体,对象具有属性、事件和方法。
2、VB的代码存储在模块中，在VB中提供了三种类型的模块：窗体模块、标准模块和类模块。
4、简述单像空间后方交会的程序设计步骤？
答：（1）、读入原始数据内方位元素，像点坐标和地面点坐标；（1分）
（2）、确定外方位元素初值，线元素用像控点平均坐标，角元素初值为零；（2分）
（3）、利用已知值和近似值组误差方程式；（1分）
(4)、解算外方位元素改正数；（1分）
(5)、改正数是否小于给定限差，结果输出；（1分）
（2）、假设区域内高程异常与大地坐标存在曲面数学关系，
利用控制点上的高程异常序列和相应的大地坐标组成多项式拟合方程组。（2分）
（3）、根据最小二乘原理求解拟合方程系数，得到校正多项式函数参数，利用所求参数计算待定高程点的高程异常，高程异常与GPS观测得到的大地高计算GPS观测点的正常高。（2分）
（4）、程序：
2、简述水准网间接平差程序设计的计算流程。
答：(1) 按文件方式进行输入, 约定数据格式；（1分）
(2) 高程近似值计算；（1分）
(3) 列立观测值的误差方程；（1分）
(4) 解误差方程并求高程平差值并进行精度估计；（2分）
(5) 计算结果输出。（1分）
山东交通学院期末考试测量程序设计课程试卷 （A）卷2010——2011 学年第 1 学期 第 2 页 共 3 页
For i = 1 To n
Input #1, dh, x1#, y1#, h1, h2
A(i, 1) = 1: A(i, 2) = x1#: A(i, 3) = y1#
A(i, 4) = x1# * x1#: A(i, 5) = y1# * x1#: A(i, 6) = y1# * y1#
L(i) = h2 - h1（2分）
Dim c(1 To m, 1 Top) As Double，a(1 To m, 1 To n) As Double, b(1 To n, 1 To p) As Double（1分）
For i = 1 To m
For j = 1 To p
For k = 1 To n
c(i, j) = c(i, j) + a(i, k) * b(k, j)（3分）
8、高斯投影坐标反算过程中，一般采用迭 代方法计算底点纬度 。
9、在圆曲线和直线之间增设缓和曲线后，整个曲线发生了变化，为了保证缓和曲线和直线相切，圆曲线应均匀地向圆心方向内移一段距离 ，计算公式为 ，切线也应相应地增长 ，计算公式为 。
10、set语句的作用是将对象引用赋给变量，createobject语句的作用是创建并返回一个对 ActiveX 对象的引用 。
1、已知 和 两点的坐标，试编制计算AB边坐标方位角的子程序（以度表示即可）。
答：Public sub fwjjs(x1,y1,x2,y2,a)
Dim x1#,y1#,x2#,y2#,dx#,dy#（2分）
dx=x(2)-x(1)
dy= y(2)-y(1)+0.0000001
s=sqr(dx*dx+dy*dy)（2分）
得分
阅卷人
二、简答题（每题6分，共36分）
1、简述程序设计的基本步骤。
（1) 问题描述：按照一定的逻辑顺序将所要解决的实际问题描述清楚；（1分）
（2) 算法设计：首先对算法的分析，设计解决问题的算法；（2分）
（3) 程序编码：按照算法，编写出满足要求的程序代码；（1分）
（4) 程序测试：程序测试是为了发现错误而执行程序的过程。（2分）
6、在默认情况下，声明的静态数组下标下界从0开始，设定数组的默认下界为1的语句为Option Base 1。
7、VB 使用工程来管理构成应用程序的所有不同的文件，一个工程主要包括工程文件（.vbp ）、窗体文件（.frm ）、标准模块文件（.bas ）、类模块文件（.cls ）、资源文件（.res ）。
3、对于给定字符串A$=“abcdefghijk”，测试字符串A$长度的语句是len(A$)，截取中间第5-8个字符“efgh”的语句是Mid$(A$,5,4)。
4、编程序进行间接平差的计算过程中，高程控制网一般选取待定点的高程为参数，平面控制网一般选取待定点的坐标为未知参数。
5、在定义数组时已确定了数组的大小，称为静态数组。在定义数组时未给出数组的大小，省略了括号中的下标，称动态数组。
'读取控制点数据
Open "d:\sj.txt" For Input As #1
Open "d:\jg.txt" For Output As #2
Input #1, n, s
ReDim A(1 To n, 1 To 6) As Double, L(1 To n) As Double
ReDim At(1 To 6, 1 To n), Naa(1 To 6, 1 To 6), W(1 To 6)
3、怎么利用高斯投影正、反算公式间接进行换带计算？
答：首先将某投影带内已知点的平面坐标（x1, y1），按高斯投影坐标反算公式求得其在椭球面上的大地坐标（B, L）；（2分）然后根据纬度和所需换算的投影带的中央子午线经度，计算该点在新投影带内的经差，（2分）再按高斯投影坐标正算公式计算该点在新投影带内的高斯平面坐标（x2, y2）。（2分）至此，就完成了高斯投影坐标的换带计算问题。
a=180-Sgn(dy)*90-Atn(dx/dy)* 57.29577951308°（3分）
E程试卷 （A）卷2010——2011学年第1学期 第3页 共3页
2、设A是m×n矩阵，B是n×p矩阵，编制矩阵相乘的子程序过程。
答：Public Sub MatrixMulti(a, b, c)
Next k（1分）
Next j（1分）
NextI（1分）
End Sub
得分
阅卷人
四、综合题（10分）
详细叙述二次曲面拟合法GPS高程计算的原理和参数计算程序（可直接调用矩阵计算子过程）并进行注释。
答：（1）、在拟合区域内选取部分均匀分布的GPS/水准点作为控制点，根据GPS观测的的大地高和水准测量得到的正常高计算控制点的高程异常.（1分）
End Sub
Next i
'计算转换参数
MatrixTrans A, At '求系数阵的转置矩阵
Matrix_Multy Naa, At, A '求AtA
kk = MRinv(Naa)
Matrix_Multy W, At, L '法方程常数向量
Matrix_Multy cs, Naa, W（2分）
Print #2, cs(1), cs(2), cs(3), cs(4), cs(5), cs(6)'六参数求解完成（1分）
5、写出用四参数进行平面直角坐标转换的公式，并解释四参数含义。
答：四参数进行平面直角坐标转换的公式包括两个平移参数，一个尺度比参数，一个旋转角参数，（2分）公式为：
（4分）
6、简述VB调用WORD的方法和主要语句（定义对象变量、创建word对象、设置对象可见、释放对象）?
答：（1）、定义对象变量Dim worddoc As Object（2分）