第二章测量程序基本算法.pptVIP

* * * * * * * * * * * * * * * * 选择法排序（升序）的VB程序： for i1=1 to n-1 ‘n-1趟遍历 P=i1 for i2=i1+1 to n ‘1次遍历 if a(P)a(i2) then P=i2 end if next i2 temp=a(i1):a(i1)=a(P):a(P)=temp next i1 思考：选择排序降序算法怎样？ 2） 冒泡排序（升序） 基本思想：(将相邻两个数比较，小的调到前头) （1）有n个数（存放在数组a(n)中），第一趟将每相邻两个数比较，小的放到前头，经n-1次两两相邻比较后，最大的数已“沉底”，放在最后一个位置，小数上升“浮起”； （2）第二趟对余下的n-1个数（最大的数已“沉底”）按上法比较，经n-2次两两相邻比较后得次大的数； （3）依次类推，n个数共进行n-1趟比较，在第j趟中要进行n-j次两两比较。 冒泡排序（升序)VB程序： for i=1 to n-1 ‘n-1趟遍历 for j=1 to n-i ‘1次遍历 if a(j)a(j+1) then temp=a(j+1) a(j+1)=a(j) a(j)=temp end if next j next i 3） 合并法排序（将两个有序数组A、B合并成另一个有序的数组C，升序） 基本思想： （1）先在A、B数组中各取第一个元素进行比较，将小的元素放入C数组； （2）取小的元素所在数组的下一个元素与另一数组中上次比较后较大的元素比较，重复上述比较过程，直到某个数组被先排完； （3）将另一个数组剩余元素抄入C数组，合并排序完成。 两个有序数组A、B合并成一个有序数组C，VB实现的程序： Do While ia = UBound(A) And ib = UBound(B) 当A和B数组均未比较完 If A(ia) B(ib) Then C(ic) = A(ia)： ia = ia + 1 Else C(ic) = B(ib)： ib = ib + 1 End If ic = ic + 1 Loop Do While ia = UBound(A) A数组中的剩余元素抄入C数组 C(ic) = A(ia) ia = ia + 1 ic = ic + 1 Loop Do While ib = UBound(B) B数组中的剩余元素抄入C数组 C(ic) = B(ib) ib = ib + 1： ic = ic + 1 Loop 2.2.5 坐标方位角计算 坐标方位角 MN：是指坐标北到某一向量的角度，范围为0～360 o。 假若存在两点M(xm,ym)与N(xn,yn)，求取方位角 ……(1) ……(2) ……..(3) ………(4) 其中， 从(2)、(3)、(4)式计算求得的方位角都有一定的局限性，就是计算的结果不一定是真正的坐标方位角，都要进行象限判断后加上一定的数值后，才能变成真正的坐标方位角。 改进后的计算公式: 式中,sgn()为取符号函数,即 ， 通过以上分析，应用式(5)计算时,除第1象限结果 大于360度外,其它象限,计算出的方位角都是真方位角，即使在 大于360时, -360也为第一象限真方位角 。 ……..(5) a=Atan((YN-YM)/(XN-XM))+180*2(1+sgn(xn-xm))/2。注意：当xn，xm不相等时，可采此算法，当xn与xm相等时，在编程时需单独考虑。 若方位角采用改进后的式（5）来计算， VB程序为： Function Qfangweijiao(xn As Double, yn As Double, xm As Double, ym As Double) As Double Dim c As Double c = gn(xn, xm)