【doc】三维导线平差计算程序
三维导线平差计算程序 铡绘技术l98.年第2期总第3日期
三维导线平差计算程序
龙一鸣
(长江艟逋局一处删量队)
近年来由于光电测距在弼量工作中迅速广泛的应用,使传统的测量方法和对成果
的
处理也发生了相应的改变.在控制网测量中,以往是将平面控制和高程控制分开来
做,而今
用光电测距仪来测量控翮网时,则可一次完成对控制网点位的三维坐标的测量亡
作,它大大
地提高工效和减轻劳动强度,深受欢迎.本文根据目前有关三维导线平差公式.编
制了
BASIC语言(PB-700)程序.
一
,
计算公式
六西格玛计算公式下载结构力学静力计算公式下载重复性计算公式下载六西格玛计算公式下载年假计算公式
设有一条三维附合导线,如图】,A,B
为已知三维坐标点,其中共有r/--1个待定
点,共观测了条边的斜距(经改化过
的);个垂直角7'和+1个水平角(左
角)卢.
图
?+w;0
(一)坐标增量及离羞由下式计算
?-Sc.s-c..
?=s-COSY-?sinaj}(1)
AH'=Stsinri
式中:口t——坐标方位角. 在实际怍【业时『义高I与镜葛,基本相 等,其J与,的差值在O.1m以内.有的使 ,等于,,这样r=T,如果J不等于,可 按下式计算出r值.
sjnr—sjn+,--J(2) 式中:r一计算得的垂直角
1"--一垂直角的观澳f值.
(二)条件方程戒如下
妻訾st一妾?H—cost一砉cr一, +(一.)dm+t~30
害一妾?圩n矾+砉c一B +(y一一',1)drn+w.0 妻1蠹叼州H~)一.
(3)
.
,
/?I_J_____??-__?l_l?lllf_J
础
含
点
日
1989年第2期总第j5期测绘技术 一,'—一
:f:——瓤边的改正数}
一——垂直角?改正敛;
——
水平角的改正数J
dm——尺度改正数1
w=一l
=:
Y—r?1
=
,彳自一日l
p=206265.
(三)观测值权的确定
各观泌值的精度可由方位角的闭合差,垂直角及斜距的对向观测互差来确定.它们
的权为
PI.
P一讲;
P一m;I—m;"一
m;'Pm.
式中:——水平角的测角中误差} r一
垂直角的测角中误差J
——
测距的中误差.
(四)法方程式为
[]+【]+['K.+[a,dK+一o
[++]Ka4-[一.
[]十[].+[].+[ePd]+[?]d+. []+[]+[]+[dPd]+[AHJd+=o
[Ax]K2+[Ay~Ka+[AH'IK4=0
(五)改正数F々计算
=
(訾一+As/-/K)
r一(一?日lcos'2,AH'sin矾K3+AHetgrlK)
f,I=(一(一一)2+(n—I)3) .
1
S=(+I)?(1+din) 式中.——经过尺度改正后的斜距平差值. (六)最弱点K的坐标平差值及K点高程平差值的权
函数
excel方差函数excelsd函数已知函数 2 f x m x mx m 2 1 4 2拉格朗日函数pdf函数公式下载
式
(4)
(5)
一,.Jl;{r
JJIl
rPP
,???l
矗
一p
12测绘技术1989年筻2期总赘35期 =
耋警一耋^H—s.一妻cy一+c一m" ,=
耋'.一As+害c一+c一m
=
委.+妾咆+
(七)最弱点h的点位中误差及h点的 高程中误差按下式计算
m.?me'l?立
mr皇?m.?
点的点位中误差:
??~,mj+mr:
h点的高程中误差.
mn=?mP1
H
式中:m.一单位权中误差(水平角中误 差).
m.一??P]
=,符号及使用
说明
关于失联党员情况说明岗位说明总经理岗位说明书会计岗位说明书行政主管岗位说明书
N——导线的边
MS一测边中误差,阻毫米为单位| MT一垂直角的测角中误差,以秒为 单位;
MB一一水平角的测角中误差,以秒为 单位|
s(I)——斜边距离,以米为单位I W(I)——闭合差,角度闭台差以秒为 单位,坐标闭合差以毫米为单 位I
B(I)一_水平角}
T(I)——按(2)计算的垂直角| R(I)一一坐标方位角|
H(I)一高程,以米为单位}
X(D一张坐标,以米为单位} (6)
Y(I)一,横坐标,以米为单位, o206265
置入数据的顺序:
700DATA导线边数N,铡边中误差 MS,垂直角的测角中误差MT,水平角测 角中误差MB.
710DATA番点号的字符,由A点至
B点.
720DATA置入巳知数据:A点的
纵,横坐标,起端的坐标方位角R(0),A 点的高程H(0),B点的纵,横坐标.末端 的坐标方位角RB,B点的高程HB. 730DATA置入观测数据:水平角
B(1),垂直角T(I),斜距S(I),………… B(N+1),T(N),S(N). 数据置入后,必须校对后方能启动运 算,运算中发出鸣声,
表
关于同志近三年现实表现材料材料类招标技术评分表图表与交易pdf视力表打印pdf用图表说话 pdf
示计算到某一步 序可打巳知数据,观测数据,平黄后
的坐标,边长,坐标方位角,待定点的高 程,单位权中误差f水平角中误差J,最弱 点的点位中误差及高程中误差.如果需再打 印平差的成果,再启动205行,即可又打印 成果表.如另外计算其它导线时,必须清除 自700行以后的数据,再置入新数据后启动 运算.
OO4"
三,源程序(磁带SONYB面159,185) 5CLEAR:READN.MS,MT,MB: K=INT(N/2):P一206.265
10DIMA(8,3?N+1),B(N+1).
C(9,9),G(N),H(N),K(5),L(N),M(N),
_?1I??__l_{J1{?}?J l__【
i989年第2期总第35期测绘技术ii P(3?N+1)
15DIMRCN+1),S(N),T(N), W(4),X(N),Y(N),S(N)
20FORI0TON:READF$(I): NEXTI:READX(0),Y(O).R(O),H(O),
XB,YB,RB,HB
24LPRINT"N="jN}'1"AB(7); "Ms一JMS;i11i/1",TAB(J8);ML一":
MTj""JCHR$(217);
27LPRINTTAB(29)j"Mb=;MB; ,CHR$(247)
30LPRINT…,F$(O);TAB(7): X2;"X(0)TAB(25);V2="JY(O) 32LPRINTTAB(7);R2:;R(O): TAB(25)J"H2=;H(0)
35LPRINTF$(N);TAB(7):"X6 一
"
;XB}TAB(25)YB:=",YB 37LPRINTTAB(7);"Rb=:RB; TAB(25);H6;:HB
40A=R(O):GOSUBd00:R(0)一B:
A-=RB:GOSUB400:RB—B
45i-"ORl=1TON:READB(I) .
T(1),S(I);L=上+S(I):EXTl:READB (N+1)
48LPR【NT…
50FORI一1TOL\:A=B(I):
GOSUB400:B(I)=B
55LPRINT;…B}I;=":A:
TAB(2o)S"J一"S(f):NEXTI
60A=B(N+]):GOSUB400:BrN
+1)=B;LPRINT;BN_1;"::A 62LPRINT
65FORI=1TON;A=ABST") GOSUB400:LPRINT…;TjI;:":
T(I)
67T(I):=SGNT(I)?B;NEXTI: GOSUB600
70W(1)=(R(N十I)一RB)?3600
W(2)=(X(N)一XS)?lea;W(3):(V (N)一YB)?lEa
75S=SQR(W(3)一2+W(2)2):
W(4)=(E{(N)一HB)?lE3
80LPRINTUSING"#{#{. }";";
Wx="Jw(2);m/21";TAB(20)}"Wy =
"
jW(3);"m?】
85LPRINTUSING"#}?.jW2
一;w(1);"";CHR$(247);TAB(2O);
Wh=;W(4);illm
9OLPRINT"1/Fs=1/tIKT(L? 1E3/S)
95FORI=1T0N:A=l/S(I):A (2,I)一L(I)?A:A(3,I)一M(I)?A:A(4, I)=G(I)?A
100IFI=KTHENFORJ=2TO 4:A(J+4,I)=A(J,I):NEXTJ 105P(I)MS^2/MB^2:NEXTr
1lOFORl=N+1TO2?NrBI—
N:A=G(B)/P:A(2,r)篇一ACOSR(B):
A(3,I)=一A?SINR(B)
112A(4,I)=A/TANT(B)
115lFB<=KTHENFORJ=2TO 4:A(J?4,I)=A.(J,I):NEXTJ 120P(I)=MT2/MB^2:NEXI l25A=0:FORI=2*N+lTO3*N:
(VB—Y(A)l/P:A A(1,r)=1:A(2.I)一
(3,I)=(XB—X(A))/P
l3OrF(r一2*N)<一KTHENA(6.
I)=一(V(K)一Y(A))/P:A(7.I)=(X(K) 一
X(A))/P
135P(I)=1:A=A+l;BEEP:NEXT I:A(1,3*N+1)=1;P(3*N+1)一1
140FORI=1TO8}F广)RJ=1TO
8:FORE=1TO3*N1
145C(I,J)=C(I,J)+A(I,E)*A(J, E)*P(E);NEXTE}NEXTJ;NEXTI 150FORr=1TO4;C(I.9)=W(I): NEXTJ
152A—N;B=0FORI一1TO2C
(5.2+B)=(X(A)一X(O))/1E3:C(5,3+ B)=(Y(A)一Y(0))/1E3
14测鐾技术1980年第2期总第35媳
155C(5,4+B)(H(A)一H(0))/
1E3:B=4:A=-K:NEXTI
157FORI=2TO8:C(I,5)一:C(5,I)
:BEEP:NEXTI
160FOR1=:1TO5:D—C(I.
I):FOR
j—I+1TO9:IFJ9rHENU=U,C
(I,)一2/D
IG5C(I,J)一,C(I,J)/DNEX],
】:FORE=1TO9:IFE=ITHEN175 170C=C(E.I):FORJI+1TO9: C(E,J)=C(E,J)+C*C(I,J):NEXTJ 175NEXTE:NEXTI:FORI=1TO 5:K(I)=C(r,9):BEEP:NEXTI 180FORI=1TON:FORJ=1TO
4:S(I)=S(I)+A(J,I)*K(J)*P(I)/1E3 185T(I)=T(I)+A(J,N+I)*K(J)* P(N+I)/3600
190B(I)=B(I)+A(J,2*N+I)*K(J) /3600:NEXTJ:NEXTI:B(N+1)B(N +1)+K(1)
200IORI一1TON:S(I)=S(I)*(1 +K(5)*IE一6):NEXTI:GOSUB600 205FORr=0TON:LPRINT….
FS(I)
210LPRINTUSING#}###}}. ##}:
TAB(2);X=";X(I);TAB(20);Y一jY(r)
2121FI—NTH匡N225
2l5DS(I+1)*C0S1,(1):LIRI NTUSING#####.}##;TAB(2);D=; D;TAB(20):
220A=R(I+1):(OSUB500:LjRj NTUSlNG"}##.
{###";R一;B,TAB
(2){"H一:H(!)
225NEXTI:LPRlNTUSING##{ #}##";TAB(2);H一:H(N)
230O—SQR(U/3):X=O*SQRC(6 6):Y=O*SQRC(7,7):H—O~+SQRC(8.8) 235M—SQR(X~2?Y^2):LPRINT
USING"#{.
#:,"m.一+一;0;
CHR$(247)
240IPRINT…j"M";K:=;USI
NG"###.
#'';h{;Ifl-'lI,mh==";EI;"r? 245END
400Ij=:fNrA1/GOINT(100~FRA CA)+1/36-~RAC(J09*~FRACA):RET URN
500B—INTA+1E一2.1NT(60{'】A
CA){GE-3.FRAC(60aFRACA):RE'l" URN
600FORI=1TON+1:R(I)一R(I
—
I)+B(I)+180
605IFR(r)>360TliENR(I)一R(I)
一
360:GOTO6O5
607IFi-二N'1Tt-]ENRETURN
610A=S(I)*COST(I):L(I)=A*CO
SR(I):M(r)一A*SINR(1):G(I)一S(I)*
SINT(I)
615I{I—HI一1)+Ij:XfI;X(I一 1)+LrI]:Y(I):==Y(I—1),M(I):NEXTI
四,算倒,
图z
(一)起算薮据厦州台差 N一6kfs一5mmM
A4lXa一6096.882 Ra一2d2.20'222 242Xb一6312.187 Rb=:24220222 B一54112
B2=138.2917 B3=1670}|
B4—178561
B5=!13484l B6=88.
1336
B7:=159.17 2.5秒川)一5秒
0一4396863
L}a一227618
Y1]=:蚰6488,I lIb一237.
934
Sl255.
548
S2=]38.
336
S3166995
S4=【78813
S5一】82.
767
S6=182635
..
;
1989年第2期总第35期测绘技术 双边平角后方交会的平差及精度分析 严伯铎
(山东省城乡建设勘察院)
随着电子速测仪的广泛应用,各种边一 角交会定点法,已成为
工程
路基工程安全技术交底工程项目施工成本控制工程量增项单年度零星工程技术标正投影法基本原理
测量中经常被采 用前一种方法.笔者巩?文献[1]中介绍过单 边单角后方交会定点的方法,谈法虽具有布 点灵活方便的优点,但在应H{中尚存在一定 限制,且缺少计算校核.本文再介绍一种双 边单角后方交会定点的方沽,文中首先讨论 几种平差计算韵方法,然后推寻出点位中误 差的计算式,同时又计算并绘示r待定点点 T1=一5.27055
T2=一1.14007
T3=一0.10275
T4:0.03276
T5—7.43422
T5=4.
11035
Wx=0.7ramWy=57.0ram
W.=".0秒Wh=4.0ram 1/Fs=1/19395
(=)平差成果硬精度评定
A4l
X=6096.882Y4396.863 D一254.371R:116. 31372
H一22761800
F5l
X=5983.
275Y--4524.
455
D138.
293R=75.
0~554
tt=203.
34093
F52
X=6019.
035Y=4758.
044
位误差椭圆的分布状况图及等误差曲线图, 并加以精度的讨论与分析,最后给出应用中 的几点说明供参考.
一
,平差计算的基本方法
?图l,A,B为已知点,P为待定点, 在P点上'侧量了两条边s=P和s:一PB, 并观测了-尸点上的水平角r,欲确定P点的坐
标.这样由于存在多余观测,便产生平差问 F53
F54
}_55
B42
D=166.981R=B2.04576 H=20036284
X=6097.215
D一178.799
H=19985440
X一5183.848
D=181.094
}{=20003.385
Y=4905.593
R=6101061
Y一5O62.002
R=354.49469
X=6364.205Y=5045.682 D=I82.134R263.03232 H=224.60950
X=6342187Y=4864.884 H237.9340
IYIO~+一4.6秒
M3=9.8ram
rah=2.
6ram