计算程序.docx
- 文档编号:5950329
- 上传时间:2023-01-02
- 格式:DOCX
- 页数:13
- 大小:23.31KB
计算程序.docx
《计算程序.docx》由会员分享,可在线阅读,更多相关《计算程序.docx(13页珍藏版)》请在冰豆网上搜索。
计算程序
工程测量(闽浙)浦南线2006~2008计算器程序归档
水准测量平差程序(SZJS)
D=A”HO”:
B”HN”:
C”N”:
T=0:
I=Z
LBI1:
G=Z[I]:
S=Z[I+C]:
{GS}:
I:
Pause1:
“H”:
D+D+G”DH”:
T=T+S:
Z[I]=G:
Z[I+C]=S:
I
Goto1⊿
“FH”:
F=D-B◢
I=1:
D=A
LBI2:
I:
Pause1
“H”:
D=D+Z[I]-FZ[I+C]/T◢
I
Goto2:
≠=>”END”⊿
备注:
HO为起点高程,HN为终点高程,S测站距离或测站数N高差个数DH测站高差FH高差闭合差H待求高差
水准测量平差程序(SZPC)
P:
D:
A:
B:
Defm8
D=D+1:
N=0:
F=0:
M=0
LBI0
N=N+1
{CKL}
Z[2N-1]=C:
F=F+C
Z[2N]=K:
M=M+K
N<0=>Goto0⊿
P>16=>“W=”:
W=0.012SQRTP≠=>“W=”:
w=0.02SQRTL>⊿
“F=”:
F=F+A-B◢
ABSF
≠=>GotoE⊿
N=0:
G=A
LBI1
“N=”:
N=N+1◢
“H=”:
G=G+Z[2N-1]+FZ[2N]◢
N
G-B
LBIE
备注:
A起点高程,B终点高程,N测段计数,C、Z[2N]观测高差,K、Z[2N-1]测段路线长度(Km)或测站数,F路线闭合差,G待求点高程,D未知水准点的数量,W在山地每公里测站数P〉=16,或高速公路和国家一级公路线长L的闭合差,Defm8是按3个未知点设置,其数值为2(n+1).
附和导线平差程序(DXPC)
P:
A:
B:
I:
J:
K:
L:
Defm12
P=P+2:
N=0:
M=A
W=60SQRTP/3600
LBI0
N=N+1
{C}
Z[2N-1]=C:
M+M+C
M>180=>GotoM=M-180”≠=>M+M+180⊿
N
Goto0⊿
“F=”:
F=M-B→DMS◢
ABSF
≠=>GotoE⊿
N=0:
M=A
LBI1
N=N+1”M=M+Z[2N-1]+F
M>180=>M=M-180:
≠=>M=M+180◢⊿
“M=”:
M→DMS◢
Z[2N-1]=M
N
Goto1⊿
N=0:
M=0”G=0:
H=0:
P=P-1
LBI2
N=N+1
(D):
Z[2N]=D:
M=M+D
X=DcosZ[2N-1]:
Y=DsinZ[2N-1]
G=G+X:
H=H+Y
Z[2N-1]=X:
Z[2N]=Y
N “fx=”: G=G+I-K◢ “fy=”: H=H+J-L◢ “f=”: T=SQRTG^2+H^2◢ “1/K”: Q=M/T◢ Q>2000=>G=-G/M: H=-H/M≠=>GotoE⊿ N=0: X=I: Y=J LBI3 “N”: N=N+1◢ D=SQRTZ[2N-1]^2+Z[2N]^2 “OX”: V=Z[2N-1]+DG◢ “OY”: U=Z[2N]+DH◢ “NX”: X=X+V◢ “NY”: Y+Y+U◢ Z[2N-1]=X: Z[2N]=Y N Goto3⊿ LBIE 备注: P未知点点数,A起点方位角,B终点方位角,I、J起点坐标K、L中点坐标C导线点观测左角,D导线边边长,X、Y未知点平差值,Defm12是按照4个未知点设定的,其数值为2(n+2),SQRT为数学根号. 竖曲线高程计算(SQXJS) Defm5: “1.ZTS”: “2.FLS”: N: A“BP”: D“BPH”: I: J: : Q“SR”: F: C: T“BT”=QAbs(J-I)/2◢ Z[1]“BE”=T^2/2/Q◢ LBI0: {S}: S<=A=>P=I: ≠=>P=J⊿Abs(S-A)>T=>Z[2]=0: ≠=>Z[2]=F⊿ K=D-P(A-S)-Z[2](T-Abs(A-S))^2/2/Q◢ C=0=>Z[5]=0: Goto1: ≠=>Z[5]=1: Goto4⊿ LBI1: H: R“PR”: O“LS”: U: Y“HY”: R<550=>W=0.06: M=170: ≠=>R<710=>W=0.05: M=150: ≠=>R<960=>W=0.04: M=130: ≠=>R<1410=>W=0.03: M=105: ≠=>R<2500=>W=0.02: M=85: ≠=>R>=2500=>W=0.02: M=0: Goto2⊿⊿⊿⊿⊿ LBI2: M>=O=>L=O: Goto3: ≠=>L=M: Goto3⊿ LBI3: S<=(Y-L)=>C=0.02H: G=-C: Goto4⊿ S<=Y=>C=H(0.02-(S-Y+L)(0.02+W)/L): G=C: AbsC<=0.02=>G=-0.02H: ⊿Goto4⊿ S<=(Y+U)=>C=-HW: G=C: Goto4⊿ S<=(Y+U+L)=>C=H(0.02-(Y+U+L-S)(0.02+W)/L): G=C: AbsC<=0.02=>G=-0.02H: ⊿Goto4⊿S>(Y+U+L)=>C=0.02H: G=-C: Goto4: LBI4: N=1=>X“PL”=(1+H)C/2-(H-1)G/2◢ V“PR”=(1+H)G/2-(H-1)C/2◢ ≠=>N=2=>X“PL”=C◢ V“PR”=C◢Goto5⊿⊿ LBI5: {Z}: Z[3]“SH”=K-(AbsE-E)X/2+(AbsE+E)V/2◢ Z≠0=>Z[4]“GH”=Z[3]-(AbsZ-Z)X/2+(AbsZ+Z)V/2◢Goto5: ≠=>C=Z[5]C: ”END”⊿ LBI0⊿ 备注: N=1为整体式,N=2分离式,A竖曲线交点里程,D竖曲线交点高程,I前坡度,J后坡度,Q竖曲线半径,F凹(-1)凸 (1),C整体式超高侧或分离式左侧横坡度,G对应C另一侧横坡度,T竖曲线切线长,K设计高程,P竖曲线计算坡度,H平曲线转向值,R平曲线半径,O对称型缓和曲线长,U圆曲线长,Y缓圆点里程,W最大超高横坡度,M超高渐变段曲线长,X左侧计算横坡值,V右侧横坡计算值,E设计高程距中桩距离(左负右正),Z距中桩偏距(左负右正),Z[1]竖曲线外失距,Z[2]判断待求里程是否在切线长绝对值以内,Z[3]路基中桩高程,Z[4]距中桩偏距Z的高程,Z[5]判断横坡是否在超高渐变段内。 平曲线坐标程序(TYQXJS) 二、源程序 1.主程序(TYQXJS) "1.SZ=>XY": "2.XY=>SZ": N: U"X0": V"Y0": O"S0": G"F0": H"LS": P"R0": R" RN": Q: C=1÷P: D=(P-R)÷(2HPR): E=180÷π: N=1=>Goto1: ≠>Goto2Δ←┘ Lbl1: {SZ}: SZ: W=Abs(S-O): Prog"SUB1": X"XS"=X◢ Y"YS"=Y◢ F"FS"=F-90◢ Goto1←┘ Lbl2: {XY}: XY: I=X: J=Y: Prog"SUB2": S"S"=O+W◢ Z"Z"=Z◢ Goto2 2.正算子程序(SUB1) A=0.1739274226: B=0.3260725774: K=0.0694318442: L=0.3300094782: F=1-L: M=1-K: X=U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW (C+FWD))+Acos(G+QEMW(C+MWD))): Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+ QELW(C+LWD))+Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))): F=G+QEW(C+ WD)+90: X=X+ZcosF: Y=Y+ZsinF 3.反算子程序(SUB2) T=G-90: W=Abs((Y-V)cosT-(X-U)sinT): Z=0: Lbl0: Prog"SUB1": L=T+QEW(C+ WD): Z=(J-Y)cosL-(I-X)sinL: AbsZ<1E-6(0.000001)=>Goto1: ≠>W=W+Z: Goto0Δ←┘ Lbl1: Z=0: Prog"SUB1": Z=(J-Y)÷sinF 三、使用说明 1、规定 (1)以道路中线的前进方向(即里程增大的方向)区分左右;当线元往左偏时,Q=-1;当线元往右偏时,Q=1;当线元为直线时,Q=0。 (2)当所求点位于中线时,Z=0;当位于中线左铡时,Z取负值;当位于中线中线右侧时,Z取正值。 (3)当线元为直线时,其起点、止点的曲率半径为无穷大,以10的45次代替。 (4)当线元为圆曲线时,无论其起点、止点与什么线元相接,其曲率半径均等于圆弧的半径。 (5)当线元为完整缓和曲线时,起点与直线相接时,曲率半径为无穷大,以10的45次代替;与圆曲线相接时,曲率半径等于圆曲线的半径。 止点与直线相接时,曲率半径为无穷大,以10的45次代替;与圆曲线相接时,曲率半径等于圆曲线的半径。 (6)当线元为非完整缓和曲线时,起点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径。 止点与直线相接时,曲率半径等于设计规定的值;与圆曲线相接时,曲率半径等于圆曲线的半径。 2、输入与显示说明 输入部分: 1.SZ=>XY 2.XY=>SZ N? 选择计算方式,输入1表示进行由里程、边距计算坐标;输入2表示由坐标反算里程和边距。 X0? 线元起点的X坐标 Y0? 线元起点的Y坐标 S0? 线元起点里程 F0? 线元起点切线方位角 LS? 线元长度 R0? 线元起点曲率半径 RN? 线元止点曲率半径 Q? 线元左右偏标志(左偏Q=-1,右偏Q=1,直线段Q=0) S? 正算时所求点的里程 Z? 正算时所求点距中线的边距(左侧取负,值右侧取正值,在中线上取零) X? 反算时所求点的X坐标 Y? 反算时所求点的Y坐标 显示部分: XS=××× 正算时,计算得出的所求点的X坐标 YS=××× 正算时,计算得出的所求点的Y坐标 FS=××× 正算时,所求点对应的中线点的切线方位角 S=××× 反算时,计算得出的所求点的里程 Z=××× 反算时,计算得出的所求点的边距 1E-6=0.000006 平曲线坐标正、反算计算程序(PQX) prog"PQX" {mn}: M"X0": N"Y0": Defm1: Z[1]=0 Lbl1: {ABGHVWK}: K"JL": G"JX": H"JY": V"I0": W"J0": A"R": B"LS": Fixm P=B^2/24/A: Q=B/2-B^3/240/A^2: T=(A+P)tan(AbsW/2)+Q: C=K-T: D=C+B: F=D+π*A*AbsW/180: E=F-B w<0=>s=-1: ≠=>s=1△ Lbl2: {L}: L: L=0=>GOTO1△ Lbl6: L≤C=>O=K-L: R=G-OcosV: U=H-OsinV: Z=V: GOTO3△ L≥F=>O=L-F+T: Z=V+W: R=G+OcosZ: U=H+OsinZ: GOTO3△ L prog"XY": Z=V+SR: R=G+(I-T)cosV-SJsinV: U=H+(I-T)sinV+JScosV: GOTO3△ L>E=>O=F-L: prog"XY": O=V+W: Z=O-SR: R=G+(T-I)cosO-SJsinO: U=H+(T-I)sinO+SJcosO: GOTO3△ O=180(L-D+B/2)/π/A: I=AsinO+Q: J=P+A(1-cosO): Z=V+SO: R=G-(T-I)cosV-SJsinV: U=H-(T-I)sinV+SJcosV Lbl3: Z[1]=0.002=>GOTO5: ≠=>O=0△R"XL"◢U"YL"◢Z"ZL"◢ Lbl4: {O}: O"YC": O=0=>GOTO2△ O=-1=>GOTO5△ X=R-OsinZ: Y=U+OcosZ: pol(X-M,Y-N: J<0=>J=J+360△ X"XZ"◢Y"YZ"◢J"A0"◢I"DD"◢GOTO4 Lbl5: {XY}: X"XF": X=-1=>Z[1]=0: GOTO4△Y"YF": pol(X-R,Y-U: O=Icos(J-Z): AbsO≤0.001=>O"DL"◢L+O◢O=Isin(J-Z): O"YC"◢GOTO5: ≠=>Z[1]=0.002: L=L+O: GOTO6△ prog"XY" U=AB: I=O-O^5/40/U^2: J=O^3/6/U-O^7/336/U^3: R=90O^2/π/U 变量名称说明: M"X0": N"Y0"为测站点坐标 K"JL": G"JX": H"JY": V"I0": W"J0": A"R": B"LS"分别为交点桩号、坐标、起始方位角、转角(左-右+)、半径、缓长 L、R"XL"、U"YL"、Z"ZL"分别为输入桩号、求出中桩坐标、即时方位角 O"YC": X"XZ"◢Y"YZ"◢J"A0"◢I"DD"分别为输入边桩距离(左-右+)、求出边桩坐标、边桩点到测站点的方位角和平距 O"YC"=0返回计算中桩坐标 O"YC"=-1进入计算桩号过程(输入坐标X"XF"、Y"YF"分别为求出桩号L+O,及至中桩边距O"YC",O"DL"为桩号计算的误差距) O"YC"=约0.00001可求中桩点到测站点方位角和平距 X"XF"=-1返回计算坐标过程 计算桩号过程时,需多次确认输入的坐标X"XF"、Y"YF"直到出现O"DL"桩号误差距较小时即可得出准确的桩号 =>前不用冒号 ≠=>前有冒号 △前后不用冒号 测量计算程序(CLJS) 1.直线段任意点坐标的计算(ZXZB) ABSTE X=A+(T-S)*cosE▲ Y=B+(T-S)*sinE▲ X+Z*cos(E+J)▲ Y+Z*sin(E+J) 注: A、B为直线上任意点坐标,S为已知点里程,T为待求点里程,E为前进方位角,J为偏角(左为-,右为+),Z为偏距 2.坐标反算(ZBFS) ABCD {AB}: A"X1": B"Y1": {CD}: C"X2": D"Y2": Pol(C-A,D-B): D-B<0=〉E=J+360: Goto3▲≠〉E=J: Goto3▲LbI3: I: E▲ 注: A、B测站点坐标,C、D为后视点坐标,E为测站至后视点方位角,I为两点间距离 3.夹角计算(JJ) ABCDEF Pol(C-A,D-B): K=J▲ Pol(E-A,F-B): L=J▲ N=L-K▲ N<0=>N=N+360▲≠>N▲ 注: A、B测站点坐标,C、D为后视点坐标,E、F为待求点坐标,N为夹角 4.竖曲线计算(SQX) BDFIJRL T=0.5*R*Abs(J-I): A=B-T: C=B+T: L≤A=〉H=D-I(B-L): Goto1△⊿ L≤B=〉H=D-I(B-L)-F(L-A)^2/2/R: Goto1△⊿ L≤C=〉H=D+J(L-B)-F(C-L)^2/2/R: Goto1△⊿ L≥C=〉H=D+J(L-B): ⊿Lbl1: H◢◢(H+Z)◢ 注: B为坡脚交点里程,D为坡脚交点高程,I为前一曲线坡度,J为后一曲线坡度,R为竖曲线半径,F为凹曲线为-1,凸曲线为+1,L为待求点里程,H为待求里程高程,Z为外矢距 曲线要素 α1=aretanI α2=aretanJ W=α1-α2 T=R*Tan(W/2) Z=R*(sec(W/2)-1) 5.圆曲线计算(YQX) ABFRHST O=(T-S)/R*90/π: D=2R*SIN*(O): E=F+HO: X=A+D*COS(E)◢ Y=B+D*SIN(E)◢ V=F+2HO: ◢ X+Z*COS(V+J)◢ Y+Z*SIN(V+J)◢ 注: AB为ZY(HY)点坐标,F为切线方位角,O为偏角,D为弦长,R为半径,E为计算方位角,H为曲线转向(左为-1,右为1),S为ZY(HY)点的里程,T为待求点里程,V为待求点前进方位角,J为偏角(左为-,右为+),Z为偏距,SORT是数学根号. 曲线要素: 切线长T=R*Tan(α/2) 曲线长L=R*α*π/180º(α以度为单位) 外矢距E=R*(sec(α/2)-1) 切曲差D=2T-L R--圆曲线半径,α--转向角 桩号计算 ZY里程=JD里程-T QZ里程=ZY里程+L/2 YZ里程=QZ里程+L/2 YZ里程=JD里程+T-D(检核公式) 6.缓和曲线计算(HQX) ABRLFHST K=Abs(S-T): N=K-K^5/40/R^2/L^2: O=K^3/6/R/L: D=SQRT(N^2+O^2): I=K^2/6/R/L*180/PI(): E=F+HI: X=A+D*cosE▲ Y=B+D*sinE▲◢ V=F+3HI: ◢ X+Z*cos(V+J)▲◢ Y+Z*sin(V+J)◢ 注: AB为ZH(HZ)点坐标,F为切线方位角,K为待求点至已知点弧长,N、O为坐标增量,L为缓和曲线长,D为待求点至已知点弦线长,E为计算方位角,I为偏角,R为半径,H为曲线转向(左为-1,右为1),S为起ZH(HZ)点的里程,T为待求点里程,V为待求点前进方位角,J为偏角(左为-,右为+),Z为偏距,SORT是数学根号,PL()为圆周率(派),Abs为绝对值,当运算YH点至HZ点缓和曲线段时S为HZ点里程,AB为HZ点坐标,F等于线路前进方位角+180,所求得的V为线路前进方位角+180。 曲线要素: 切线长T=(R+P)*Tan(α/2)+m 曲线长L=R*(α-2β0)*π/180º+2l0=R*α*π/180º+l0 外矢距E=(R+P)*sec(α/2)-R 切曲差D=2T-L α=3β0 曲线参数β0=l0/2/R(弧度)=l0*180º/2/R/π(度) 半径变更率C=R*l0 内移量P=l02/24/R-l04/2688/R3≈l02/24/R 切垂距m=l0/2-l03/240/R2≈l0/2 R--圆曲线半径,α--转向角,l0--缓和曲线的弧长 桩号计算 ZH里程=JD里程-T HY里程=ZH里程+l0 QZ里程=HY里程+(L/2-l0) YH里程=QZ里程+(L/2-l0) HZ里程=YH里程+l0 HZ里程=JD里程+T-D(检核公式) 7.卵形曲线计算(LQX) ABFLRQHST C=RQL/Abs(R-Q): K=ABS(S-T): M=C/R U=(MK+K^2/2)/C*180/PI() N=K-K^3(M^2/3+MK/4+K^2/20)/2/C^2+K^3(M^4/5+M^3K/3+3M^2K^2/14+MK^3/16+K^4/144)/24/C^4: O=K^2(M+K/3)/2/C-K^4(M^3/2+3M^2K/5+MK^2/4+K^3/24)/12/C^3+K^6(M^5/3+5M^4K/7+5M^3K^2/8+5M^2K^3/18+MK^4/16+K^5/176)/240/C^5: D=SQRT(N^2+O^2): E=F+H*TAN-1(O/N): X=A+D*COS(E)◢ Y=B+D*SIN(E)◢ V=F+HU: X+Z*cos(E+J)▲ Y+Z*sin(E+J) 注: AB为曲线起点坐标,大半径半径方向点(HY或YH)的X、Y、切线方位角,C、M为曲线参数,N、O坐标增量,D为曲线起点只待求点弦线长,K为曲线起点只待求点弧长,L为曲线长,R为大半径,Q为小半径,H为曲线转向(左为-1,右为1),S为起点的里程,T为待求点里程,E为计算方位角,F曲线起点方位角,V为待求点前进方位角,Z为偏距,J为偏角。 8.坐标计算及后方交会 "1.ZXZB": "2.HFJH": V: A"X1"B"Y1"C"X2"D"Y2": Pol((C-A),(D-B)): D-B<0=>F=J+360: Goto0: ≠=>F=J: Goto0: Lbl0: I: F◢ V=1=>Goto1: ≠=>Goto2 Lbl1: {STZE}: L=Abs(S-T): X=ALcosF◢ Y=B+LsinF◢ X+Zcos(F+E)◢ Y+Zsin(F+E)◢ Goto1 Lbl2: {HOMN}: E=cos-1((I2+M2-N2)/2/I/M): Q=cos-1((I2+N2-M2)/2/I/N): K=(180-O-Q-E)/3◢ P=F+H(E+K): X=A+McosP◢ Y=B+MsinP◢ Goto2 注: AB为第一坐标,CD为第二坐标,I为两点间距,F为方位角,S为第一座标点里程,T为待求点里程,Z为偏距,H为曲线转向(左为负,右为正)O为待求点到两已知坐标交角,M为第一边距,N为第二边距,E为第一偏角,Q为第二偏角,K为内角补偿数,P为计算方位角。 任意多边形面积计算程序 Lbl0: M: N: E=M: F=N: S=0↙ Lbl1: {XY}: S=S-(XN-YM)/2: X=E=>Y=F=>"S=": S◢Goto0⊿⊿↙ M=X: N=Y: Goto1 动行时,M、N输入第一点的坐标,以后提示X、Y则依次按顺时针或逆时针方向输入各点的坐标,当所有点输入完成后,还需要再次输入第一点的坐标,这样就会自动显示计算结果。 计算结果。 不论顺时针或逆时针方向输入,都有可能导致计算的结果为负,取其绝对值即可。 三维断面面积体积测量CASIO4850 F1A3 L1T=0: H=0: WQ L2Lbl0: {NEFG}: NEFG: A=0: P=A: D=G: K=0: M=1 L3Lbl1: {XYZ}: XYZ: I=POI(X-E,Y-F): J≤0=>J=J+360⊿J≤180=>U=1: ≠>U=-1⊿C=UI: K=K+(A+C)(Z-G)/2: A=C: G=Z: M=M+1 L4M L5T≠1=>GOTO2: ≠>L=G-H: V=(R+K+√(RK))(Q-H)/3◢⊿W=W+V◢ L6Lbl2: R=K: H=G: T=1: {Q}: Q: GOTO0 2、说明 (1)、程序可自动计算每一断面面积,当进行到第二个断面时就会显示出1~2
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算 程序