【正文】 入A點(diǎn)的Y坐標(biāo)=)。 double AY = (())。 (請輸入A點(diǎn)的高程=)。 double HA = (())。 (請輸入B點(diǎn)的X坐標(biāo)=)。 double BX = (())。 (請輸入B點(diǎn)的Y坐標(biāo)=)。 double BY = (())。 (請輸入B點(diǎn)的高程=)。 double HB = (())。 (請輸入N點(diǎn)的X坐標(biāo)=)。 double NX = (())。 (請輸入N點(diǎn)的Y坐標(biāo)=)。 double NY = (())。 //計(jì)算M到A的方位角 double α0 = (MX, MY, AX, AY)。 //計(jì)算B到B的方位角 double α1 = (BX, BY, NX, NY)。 Listdouble 坐標(biāo)增量X集合 = new Listdouble()。 Listdouble 坐標(biāo)增量Y集合 = new Listdouble()。 Listdouble 高差集合 = new Listdouble()。 Listdouble 水平距離平方集合 = new Listdouble()。 Listdouble 方位角集合 = new Listdouble()。 int i = 1。 double 坐標(biāo)增量X和 = 0。 double 坐標(biāo)增量Y和 = 0。 double 坐標(biāo)增量X絕對值和 = 0。 double 坐標(biāo)增量Y絕對值和 = 0。 double 水平距離和 = 0。 double 高差和 = 0。 double 水平距離平方和 = 0。 double α = α0。 do { (請輸入{0}點(diǎn)到{1}點(diǎn)的斜距D直接回車結(jié)束輸入=, i, i + 1)。 string str = ()。 if (str == ) break。 double D = (str)。 (請輸入{0}點(diǎn)的水平角左角為正,右角為負(fù)=, i)。 double β = ((()))。 (請輸入{0}點(diǎn)到{1}點(diǎn)的垂直角=, i, i + 1)。 double 垂直角 = ((()))。//并化為弧度 (請輸入{0}的儀器高=, i)。 double 儀器高 = (())。 (請輸入{0}的反光鏡高=, i + 1)。 double 反光鏡高 = (())。 i++。 //計(jì)算高差 double 水平距離。 double h = (D, 垂直角, 儀器高, 反光鏡高, out 水平距離)。 (水平距離 * 水平距離)。 (h)。 水平距離平方和 += 水平距離 * 水平距離。 高差和 += h。 //計(jì)算坐標(biāo)增量 double dx, dy。 α = (0, 0, β, 水平距離, α, out dx, out dy)。 (α)。 (dx)。 (dy)。 坐標(biāo)增量X和 += dx。 。 坐標(biāo)增量Y和 += dy。 。 坐標(biāo)增量X絕對值和 += (dx)。 坐標(biāo)增量Y絕對值和 += (dy)。 水平距離和 += 水平距離。 } while (true)。 (請輸入最后一個連接角左角為正,右角為負(fù)=)。 double βn = ((()))。 //推算B到N的方位角 double αBN = (方位角集合[],βn)。 //計(jì)算方位角閉合差 double 方位角閉合差 = αBN α1。 double 方位角改正數(shù) = 方位角閉合差 / ( + 1)。 //計(jì)算高差閉合差 double 高差閉合差 = HA + 高差和 HB。 //計(jì)算單位長度改正數(shù),三角高程的權(quán)與距離的平方成反比 double 單位長度改正數(shù) = 高差閉合差 / 水平距離平方和。 //先對角度進(jìn)行平差計(jì)算,并重新計(jì)算坐標(biāo)增量 ()。 ()。 坐標(biāo)增量X和 = 0。 坐標(biāo)增量Y和 = 0。 坐標(biāo)增量X絕對值和 = 0。 坐標(biāo)增量Y絕對值和 = 0。 for (int j = 0。 j 。 j++) { 方位角集合[j] += 方位角改正數(shù) * (j + 1)。 //計(jì)算坐標(biāo)增量 double 平距 = (水平距離平方集合[j])。 double dx = 平距 * (方位角集合[j])。 double dy = 平距 * (方位角集合[j])。 (dx)。 (dy)。 坐標(biāo)增量X和 += dx。 。 坐標(biāo)增量Y和 += dy。 。 坐標(biāo)增量X絕對值和 += (dx)。 坐標(biāo)增量Y絕對值和 += (dy)。 } //計(jì)算閉合差 double 閉合差X = AX + 坐標(biāo)增量X和 BX。 double 閉合差Y = AY + 坐標(biāo)增量Y和 BY。 //計(jì)算單位長度改正數(shù) double X坐標(biāo)單位長度改正數(shù) = 閉合差X / 坐標(biāo)增量X絕對值和。 double Y坐標(biāo)單位長度改正數(shù) = 閉合差Y / 坐標(biāo)增量Y絕對值和。 //計(jì)算導(dǎo)線相對精度 int f = (int)(水平距離和 / (閉合差X * 閉合差X + 閉合差Y * 閉合差Y))。 (導(dǎo)線的相對精度為1/{0},f)。 //計(jì)算每點(diǎn)坐標(biāo) double X = AX。 double Y = AY。 double H = HA。 for (int j = 0。 j 。 j++) { H += 高差集合[j] + 單位長度改正數(shù) * 水平距離平方集合[j]。 X += 坐標(biāo)增量X集合[j] + X坐標(biāo)單位長度改正數(shù) * (坐標(biāo)增量X集合[j])。 Y += 坐標(biāo)增量Y集合[j] + Y坐標(biāo)單位長度改正數(shù) * (坐標(biāo)增量Y集合[j])。 (第{0}點(diǎn)的X坐標(biāo)={1},Y坐標(biāo)={2},高程H={3},j+2,X,Y,H)。 } (最后一個點(diǎn)的坐標(biāo)和高程應(yīng)該與B點(diǎn)的坐標(biāo)和高程相同!)。 } } class 導(dǎo)線 { // static public double DEG(double ang) { int fuhao = (int)(ang / (ang))。 ang = (ang)。 int d = (int)ang。 int m = ((int)(ang * 100)) d * 100。 double s = ang * 10000 m * 100 d * 10000。 return ((d + m / + s / ) * fuhao) / * 。 } //計(jì)算方位角,返回弧度值 public static double 方位角(double x1, double y1, double x2, double y2) { double deltaX = x2 x1。 double deltaY = y2 y1。 double angle = * 。 if ((deltaX) ) { angle = (deltaY, deltaX)。 } if (angle 0) { angle += 。 } if (deltaY ) { angle += 。 } return angle。 } //高差計(jì)算,返回高差 public static double 高差(double 斜距, double 垂直角, double 儀器高, double 反光鏡高, out double 水平距離) { double h = 斜距 * (垂直角) + 斜距 * (垂直角) / * 斜距 * (垂直角) * (1 ) + 儀器高 反光鏡高。 水平距離 = 斜距 * (垂直角)。 return h。 } //計(jì)算坐標(biāo),返回已知點(diǎn)到計(jì)算點(diǎn)的方位角 public static double 坐標(biāo)(double x0, double y0, double 左角, double 水平距離,double 已知方位角,out double x,out double y) { double 方位角 = 已知方位角 + 左角 + 。 //將方位角調(diào)整到0到2π之間 if (方位角 = * 2) 方位角 = * 2。 if (方位角 ) 方位角 += * 2。 x = x0 + 水平距離 * (方位角)。 y = y0 + 水平距離 * (方位角)。 return 方位角。 } //根據(jù)后視邊的方位角與左角,計(jì)算前進(jìn)邊的方位角 public static double 方位角(double 后視邊方位角, double 左角) { double 方位角 = 后視邊方位角 + 左角 + 。 //將方位角調(diào)整到0到2π之間 if (方位角 = * 2) 方位角 = * 2。 if (方位角 ) 方位角 += * 2。