大地测量学编程实习报告--大地主题解算.docx

大地测量学编程实习报告 ——大地主题解算 学院： 测绘学院 专业： 测绘工程 班级： 2014级2班 学号： 2014301610339 姓名： 何杨 实习目的 1.?提高运用计算机语言编程开发的能力；? 2.?加深对大地主题解算计算公式及辅助参数的理解并掌握计算步骤；? 3.?通过编程语言实现大地主题解算。 二、注意事项 1.?计算所需变量多，容易混淆； ?2.?正反算函数的编写； ?3.?函数调用；? ?弧度与角度之间的转化。 程序框图 B01 = S * Math.Cos(A12) / M;L01 = S * Math.Sin(A12)/ (N * Math.Cos(B1));A01= S *Math.Sin(A12)*Math.Tan(B1)/N; B01 = S * Math.Cos(A12) / M; L01 = S * Math.Sin(A12)/ (N * Math.Cos(B1)); A01= S *Math.Sin(A12)*Math.Tan(B1)/N; 开始 输入e2,c,L1,B1,A12,S 正算 输入e2,c,L1,B1,A12,S Am = A12 + 0.5 * A01; Am = A12 + 0.5 * A01; Bm = B1 + 0.5 * B01; B02 = B1 + B02;L02 = L1 + L02; B02 = B1 + B02; L02 = L1 + L02; A02 = A12 + A02+PI; Yes 输出B2,L2,A21 输出B2,L2,A21 结束Am = A12 + 0.5 * A02;Bm = B1 + 0.5 * B02; 结束 Am = A12 + 0.5 * A02; Bm = B1 + 0.5 * B02; 反算 输出S,A12，A21S = U/ Math.Sin(Am);A12 = Am - 0.5 * dA;A21 = Am + 0.5 * dA + Math.PI;判断Am的象限 Am=arctan(U/V)U=SsinAm = r01 * dL + r21 * dB * dB * dL + r03 * dL * dL * dL;V=ScosAm = S10 * dB + S12 * dB * dL * dL + S30 * dB * dB * dB;开始输入e2,c,L1,B1,L2,B2Am = A12 + 0.5 * A01;Bm = B1 + 0.5 * B01; 输出S,A12，A21 S = U/ Math.Sin(Am); A12 = Am - 0.5 * dA; A21 = Am + 0.5 * dA + Math.PI; 判断Am的象限 Am=arctan(U/V) U=SsinAm = r01 * dL + r21 * dB * dB * dL + r03 * dL * dL * dL; V=ScosAm = S10 * dB + S12 * dB * dL * dL + S30 * dB * dB * dB; 开始 输入e2,c,L1,B1,L2,B2 Am = A12 + 0.5 * A01; Bm = B1 + 0.5 * B01; 源代码 调用的函数： public double Mjisuan(double B)//M计算 { double M; M = a * (1 - e2) *Math.Pow((1 - e2 * Math.Sin(B) * Math.Sin(B)),-1.5); return M; } public double Njisuan(double B)//N计算 { double N; N = a * Math.Pow((1 - e2 * Math.Sin(B) * Math.Sin(B)), -0.5);