第十二讲大地坐标系与大地极坐标系的关系(二)大地问题正反解的定义归化纬度与大地纬度的关系(会推)理解并记住贝塞耳大地问题解算公式的三个投影条件和解算步骤理解并记住贝塞耳微分方程(即dS与dσ、dl与dλ的关系)请思考野外测量可以获得哪些观测值,点的最终水平坐标和高程怎么得到?12)(11,LB)(22,LB1A2ASNL1,B1,S,A1L2,B2,A2大地问题正解大地问题反解L1,B1,L2,B2S,A1,A212)(11,LB)(22,LB1A2A2AS18022AAN21tantaneBu内容回顾ReviewMmusincossinMmu222sincos1cos4、贝塞尔微分方程的解:S与σ的关系cos1cos12222duedldueadsu2mM0902090A220QPP球面三角形mekdMkbdscos'sin1220220sin1dMkbdsS24cos2sin2cossinMCMBAbS)(Sdueads22cos1Mmu222sincos1cos4、贝塞尔微分方程的解:S与σ的关系442421281161464341kCkkBkkAmek222cosmekdMkbdscos'sin1220220sin1dMkbdsS24cos2sin2cossinMCMBAbS)(Sdueads22cos1Mmu222sincos1cos4、贝塞尔微分方程的解:S与σ的关系24cos2sin2cossin''MMS)(S442421281''''814''''64741''''kACkkABkkbbA1u1mM0901090A1111sincossincos/tantanAumAuM24cos2sin2cossinMCMBAbS110QPP球面三角形4、贝塞尔微分方程的解:S与σ的关系duAdcossinumM090A090umAcossinsindumd2cossin4、贝塞尔微分方程的解:l与λ的关系cos1cos12222duedldueads220QPP球面三角形duedl22cos1dumd2cossindumuedl222cossincos1dueemddl)cos82(sin242024200)cos82(sindueemdllMmu222sincos1cosMml2cossin''''sin''''mekkeeke22222222cos''''16''161821')(lMml2cossin''''sin'''')(l4、贝塞尔微分方程的解:l与λ的关系111,,,ASLB222,,ALBN)L,(BP111)L,(BP2221AS2A2A四、贝塞尔大地问题正解计算解算步骤:1)将椭球面元素投影到球面上121tan1tanBeu11costantanAuM11sincossinAum24cos2sin2cossinMMS四、贝塞尔大地问题正解计算111,,,ASLB222,,ALB2)解算球面三角形MmAAcostan'tantan22MAutancostan2212222111tansintansintantansintansintanAuMmAuMm解算步骤:1)将椭球面元素投影到球面上四、贝塞尔大地问题正解计算111,,,ASLB222,,ALB2A3)将球面元素换算到椭球面上222tan'1tanueB24cos2sin'2cossin''sinMMmllLL122)解算球面三角形解算步骤:1)将椭球面元素投影到球面上四、贝塞尔大地问题正解计算111,,,ASLB222,,ALBA1在第Ⅰ象限:m、M、λ1均在第I象限λ2与(M+σ)同象限A2在IV或III象限四、贝塞尔大地问题正解计算111,,,ASLB222,,ALB四、贝塞尔大地问题正解计算111,,,ASLB222,,ALBN)L,(BP111)L,(BP2221AS2A2211,,,LBLBSAA,,21五、贝塞耳大地问题反解计算解算步骤:1)将椭球面元素投影到球面上222121tan1tantan1tanBeuBeu?lMml2cossin''''sin''''l0022210coscoscossinsincosuuuu00210110sinsincoscossincossinuuAum0121001cossintancossintanuuuA01010tansinsintanAmuM五、贝塞耳大地问题反解计算2211,,,LBLBSAA,,212)解算球面三角形2122costancoscossintanuuuA3)将球面元素换算到椭球面上24cos2sin2cossin1MMS解算步骤:1)将椭球面元素投影到球面上五、贝塞耳大地问题反解计算2211,,,LBLBSAA,,211、大地问题正解已知试编程求L2、B2、A2?(以度分秒的形式输出到屏幕,保留小数点后3位;使用我们课堂上所给的公式,系数由克拉索夫斯基椭球参数a,α计算)注:这部分内容安排在课程结束到计算中心统一实习,不过,希望大家现在就要抽空开始编写程序,否则会来不及;写的程序要填入我们的程序框架中,即大地测量学基础程序框架.cpp和subroutine.h。算例mSABL2.15000000,33.000010022.000035,11.000090o1o1o1上机编程计算实习:第三部分1、大地问题正解要求1:从文件中读入已知数据,文件名及文件内容为:SUBL.txt:6378245.0298.3350000.22900000.111000000.3315000000.2要求2:每个功能模块单独编写一个子函数。文件名aα=1:298.3B1L1A1S上机编程计算实习:第三部分2、大地问题反解已知试编程求S、A1、A2?(以米、度分秒的形式输出到屏幕,保留小数点后3位;使用我们课堂上所给的公式,系数由克拉索夫斯基椭球参数a,α计算)注:这部分内容安排在课程结束到计算中心统一实习,不过,希望大家现在就要抽空开始编写程序,否则会来不及;写的程序要填入我们的程序框架中,即大地测量学基础程序框架.cpp和subroutine.h。算例69.209230,33.049521522.000035,11.000090o2o2o1o1BLBL上机编程计算实习:第三部分2、大地问题反解要求1:从文件中读入已知数据,文件名及文件内容为:SUSA.txt:6378245.0298.3350000.22900000.11-302920.9642155904.333要求2:每个功能模块单独编写一个子函数。文件名aα=1:298.3B1L1B2L2上机编程计算实习:第三部分贝塞尔微分方程的解•1、如何简化和正确判断贝塞尔大地问题解算的象限判断问题?(可收集网上资料)•2、大地问题正反解的实际编程计算