空中三角电模拟平差器的构造及精度

一、空中三角平差公式及其改化在空中三角测量中,习惯采用下列平差公式：δh=C_0+C_1x+C_2y+C_3x~2+C_4xy式中：δh——高程改正数x、y——各点平面座标C_0…C_4——倾斜及扭曲各项系数     

等时角法同时测定经纬度、方位角

本文提出利用三星在十几分钟内连续通过同一时圈（等时角）瞬间,测定其间水平角,即可同时确定经纬度和方位角,而不需要知道测站经纬度、方位角的任何概值。如果只测方位角和纬度,则仅要秒表或普通手表,而不需天文表和收报机。从理论推导和实际观测结果,都证实本法用于方位角观测特别有利,其精度高、观测易、原理和计算简单。拱极星方位角θ_3公式有：ctgθ_3=K·ctgβ_1+（K-1）ctgβ_2式中：系数K=(tgδ_1+tgδ_3)/(tgδ_1+tgδ_2),而δ_i为三星赤纬已知量;β_1、β_2为三星间两个水平角观测量。天体方位角与水平角间误差关系式,有mθ_a=±1/(60)mβ。观测试验成果列于文末,算例见附表,从而清楚地见到本法精度及计算过程。     

补助点设在任意点的位置,从高斯投影带到相邻带及从兰孛氏割圆锥投影到高斯投影的坐标变换公式(续二,完)

（C）实用之部C_1.Cassini曲线Cassini曲线表示（ρ~4-2ρ~2cos2ψ+1）~(1/2)=常数的曲线,见图21,此处ρ=2r/s。图20内设有一点a,其坐标为（ξ）=12km,（η）=12km,s=24km,则r=2~(1/2)·12km,故ρ=2~(1/2),此处ψ=45°。由ρ=2~(1/2)=1.414（按图21内ρ的比例尺）及ψ=45°,得图21的A点,在A点读得：（ρ~4-2ρ~2cos2ψ+1）~(1/2)=5~(1/2)。C_2.Hunger公式及Schroeder-Kastner公式的应用（a）Hunger公式经改进后写成下列各式（参考（A）内式（18）,（19））：     

试评60分制圆分度光电检验仪的精度

本文试就研制单位多次对一块存档光栅度盘(64800线/φ300毫米)和一块存档光学度盘(720线/φ230毫米)检测的一万多数据中,整理分析出五种检测方法,进行了精度计算和比较。建议评价这一类型高精度检验仪,检测方法宜采用全组合比较法,从而计算仪器测量一次所取得检定值的测量中误差m仪为主精度指标,它的计算式为: m仪=±(m_γ~2+τ~2)~1/2 并以仪器测量一次观测值中误差m_γ为辅助精度指标,它的计算式为: m_γ=±([VV]/(s-1)(s-2))~1/2 以上两式中: τ为仪器标准器的直径全中误差: V为观测值的改正数: s为检测的直径数,本文例s为18。至于宜采用多少条直径为最有代表性,尚需进一步研究。     

大三角网的矩阵分组平差

本文论述了用矩阵法作大三角网分组平差的问题。作者的目的在于简化平差计算。首先,通过适空地排列法方程组,建立了一个分组平差的方案。此法把三角网按典型图形来进行分组,从而使法方程中主对角线上子矩阵的结构具有一定的规律性。作者用普通累代法同样得到了分组累进式K=b_0+b_1+b_2+…+b_s+…,并进一步得到了它的吉德尔累代法的算式。用公式1/P_F=[ff]+[qφ]来计算平差值函数之权倒数。其次,为了加速分组累进式的收敛性,作者找得了在累代矩阵B存在一小最大特征值λ_1,且|λ_1|＜1或二个最大特征植λ_1=-λ_2,且|λ_1|＜1时的收敛公式。第三,作者建议用汉森法或加边矩阵法以求逆矩阵A_(11)~(-1),同时还给出了一些用加边法求逆矩阵的计算用表。最后是一个21个法方程的算例。该例子由作者一人计算完成仅化去了16小小时。如果由三人同时进行计算,那末在8小时内是一定可以完成的。用此法进行计算速度是很快的,且计算精度还不受到限制。因而此法极适用于用手摇计算机或快速计算机作大三角网的平差计算。     

局部地形改正快速计算的GPU并行的棱柱法

在重力归算中,局部地形改正在重力勘探、地壳结构分析和大地水准面计算等领域有着重要意义,但严格棱柱体积分公式计算效率低,而快速计算公式则会降低计算精度。本文利用CUDA并行编程平台,提出一种地形格网重新编码和严格棱柱体积分八分量拆解方法,实现了基于CPU+GPU异构并行技术的严格棱柱体积分计算地形改正快速并行算法,克服了GPU各个线程计算任务分配和线程计算超载问题,解决了局部地形改正的高分辨率、高精度严密公式的快速计算难题。通过试验,在显卡型号为Tesla V100的计算机上进行4°×6°范围,积分半径40'和分辨率1'的局部地形改正计算仅需1.5 s;分辨率10″的局部地形改正计算仅需14.6 min;进行分辨率3″的地形改正计算耗时45.7 h,而传统串行算法则难以完成计算。在保证微伽级以上计算精度的条件下,计算加速比最高达到850倍以上,有效缩短了计算耗时,提高了计算效率。本文还依据上述并行算法对全国范围地形改正量进行计算。结果表明,我国地形改正量普遍低于80 mGal（1 Gal=10^-2 m/s²）,平均值1.83 mGal,最大值达到196 mGal。