测量程序设计的基本函数.pptVIP

测量程序设计 第1章 常用计算及辅助函数 §1-1 概述 一、课程主要内容 第一章 常用辅助函数 第二章 法方程式的解算 第三章 常用大地测量计算 第四章 高程网平差（水准网平差） 第五章 平面控制网平差（导线网平差） 第六章 GPS基线网平差（二维约束平差） 第一章 常用计算及辅助函数 二、软件开发的一般步骤 （1）需求分析 （2）系统设计 （3）编写代码 （4）系统测试 （5）系统维护 三、作业问题 第一章 常用计算及辅助函数 §1-2角度格式的相互化算 一、度、分、秒转换弧度 为方便和规则起见，通常在输入角度时采用度分秒连写的形式，如“53°3′12.78″”输入时为“53.031278”。分、秒的整数部分都要补齐两位。 函数名称： Double dms_rad (Double a) 第一章 常用计算及辅助函数 private double dms_radc(double a) { //提取符号 double sign = (a 0) ? -1.0 : 1.0; a = System.Math.Abs(a); //提取整度 int dd = (int)(System.Math.Floor(a)); a = a - dd; //if (a 0.0) { dd = dd - 1; a = a + 10000; } //提取角度值整分、秒 int mm = (int)System.Math.Floor(a * 100.0); a = (a * 100.0 - mm) * 100.0; //if (a0.0){dd=dd-1;a=a+10000;} a = sign * (dd * 3600.0 + mm * 60.0 + a) / const1.p2; return a; } 第一章 常用计算及辅助函数 二、弧度转换为度分秒 public double rad_dmsc(double a) { //提取符号 double sign = (a 0) ? -1.0 : 1.0; a = System.Math.Abs(a); a = a * 57.2957795130823; //换成度的形式 int d = (int)(System.Math.Floor(a)); a = a - d; int m = (int)System.Math.Floor(a * 60.0); a = (a - m / 60.0) * 3600.0; if (a + 0.0001 60) { a = 0; m = m + 1; } if (m = 60) { m = 0; d = d + 1; } a = d + m / 100.0 + a / 10000.0; return a; } 第一章 常用计算及辅助函数 §1.3 坐标方位角的计算 一、C#中函数 ATAN 和 ATAN2的区别 二．坐标计算方位角 注意问题： 与数学上坐标系的关系 考虑特殊情况，例如：北坐标差等于0的情况 单位问题 不能大于360°或2π 角度是弧度表示的，我们通常用度分秒表示。 第一章 常用计算及辅助函数 public double xyazimuth(double x1, double y1, double x2, double y2) { double cc = 1, cx, cy; cx = x2 - x1; cy = y2 - y1; if (cx == 0) if (cy 0) cc = const1.P / 2.0; else cc = const1.P * 1.5; else cc