计算程序文档格式.docx
- 文档编号:18941803
- 上传时间:2023-01-02
- 格式:DOCX
- 页数:13
- 大小:23.31KB
计算程序文档格式.docx
《计算程序文档格式.docx》由会员分享,可在线阅读,更多相关《计算程序文档格式.docx(13页珍藏版)》请在冰豆网上搜索。
G=G+Z[2N-1]+FZ[2N]◢
D=>
Goto1⊿
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:
J:
K:
L:
Defm12
P=P+2:
M=A
W=60SQRTP/3600
{C}
Z[2N-1]=C:
M+M+C
M>
180=>
GotoM=M-180”≠=>
M+M+180⊿
P=>
Goto0⊿
F=M-B→DMS◢
F=-F/P:
GotoE⊿
N=N+1”M=M+Z[2N-1]+F
M=M-180:
M=M+180◢⊿
“M=”:
M→DMS◢
Z[2N-1]=M
Goto1⊿
M=0”G=0:
H=0:
P=P-1
LBI2
(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
PGoto2⊿
“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≠=>
X=I:
Y=J
LBI3
“N”:
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◢
Goto3⊿
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”:
:
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⊿
H:
R“PR”:
O“LS”:
U:
Y“HY”:
R<
550=>
W=0.06:
M=170:
710=>
W=0.05:
M=150:
960=>
W=0.04:
M=130:
1410=>
W=0.03:
M=105:
2500=>
W=0.02:
M=85:
R>
=2500=>
M=0:
Goto2⊿⊿⊿⊿⊿
=O=>
L=O:
Goto3:
L=M:
Goto3⊿
LBI3:
=(Y-L)=>
C=0.02H:
G=-C:
=Y=>
C=H(0.02-(S-Y+L)(0.02+W)/L):
G=C:
AbsC<
=0.02=>
G=-0.02H:
⊿Goto4⊿
=(Y+U)=>
C=-HW:
=(Y+U+L)=>
C=H(0.02-(Y+U+L-S)(0.02+W)/L):
⊿Goto4⊿S>
(Y+U+L)=>
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:
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÷
D=(P-R)÷
(2HPR):
E=180÷
π:
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:
SUB2"
S"
=O+W◢
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:
L=T+QEW(C+
WD):
Z=(J-Y)cosL-(I-X)sinL:
AbsZ<
1E-6(0.000001)=>
Goto1:
W=W+Z:
Goto0Δ←┘
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"
N"
Defm1:
Z[1]=0
Lbl1:
{ABGHVWK}:
K"
JL"
JX"
JY"
I0"
W"
J0"
A"
B"
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<
s=-1:
s=1△
Lbl2:
{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:
L<
O=L-C:
XY"
Z=V+SR:
R=G+(I-T)cosV-SJsinV:
U=H+(I-T)sinV+JScosV:
L>
E=>
O=F-L:
O=V+W:
Z=O-SR:
R=G+(T-I)cosO-SJsinO:
U=H+(T-I)sinO+SJcosO:
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}:
YC"
O=0=>
GOTO2△
O=-1=>
GOTO5△
X=R-OsinZ:
Y=U+OcosZ:
pol(X-M,Y-N:
J<
J=J+360△
XZ"
◢Y"
YZ"
◢J"
A0"
◢I"
DD"
◢GOTO4
Lbl5:
{XY}:
XF"
X=-1=>
Z[1]=0:
GOTO4△Y"
YF"
pol(X-R,Y-U:
O=Icos(J-Z):
AbsO≤0.001=>
DL"
◢L+O◢O=Isin(J-Z):
◢GOTO5:
Z[1]=0.002:
L=L+O:
GOTO6△
U=AB:
I=O-O^5/40/U^2:
J=O^3/6/U-O^7/336/U^3:
R=90O^2/π/U
变量名称说明:
为测站点坐标
分别为交点桩号、坐标、起始方位角、转角(左-右+)、半径、缓长
L、R"
、U"
、Z"
分别为输入桩号、求出中桩坐标、即时方位角
分别为输入边桩距离(左-右+)、求出边桩坐标、边桩点到测站点的方位角和平距
=0返回计算中桩坐标
=-1进入计算桩号过程(输入坐标X"
、Y"
分别为求出桩号L+O,及至中桩边距O"
,O"
为桩号计算的误差距)
=约0.00001可求中桩点到测站点方位角和平距
=-1返回计算坐标过程
计算桩号过程时,需多次确认输入的坐标X"
、Y"
直到出现O"
桩号误差距较小时即可得出准确的桩号
=>
前不用冒号
前有冒号
△前后不用冒号
测量计算程序(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:
E▲
A、B测站点坐标,C、D为后视点坐标,E为测站至后视点方位角,I为两点间距离
3.夹角计算(JJ)
ABCDEF
K=J▲
Pol(E-A,F-B):
L=J▲
N=L-K▲
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:
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
α=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:
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)):
F=J+360:
Goto0:
F=J:
Lbl0:
F◢
V=1=>
Goto1:
Goto2
{STZE}:
L=Abs(S-T):
X=ALcosF◢
Y=B+LsinF◢
X+Zcos(F+E)◢
Y+Zsin(F+E)◢
Goto1
{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◢
AB为第一坐标,CD为第二坐标,I为两点间距,F为方位角,S为第一座标点里程,T为待求点里程,Z为偏距,H为曲线转向(左为负,右为正)O为待求点到两已知坐标交角,M为第一边距,N为第二边距,E为第一偏角,Q为第二偏角,K为内角补偿数,P为计算方位角。
任意多边形面积计算程序
M:
E=M:
F=N:
S=0↙
S=S-(XN-YM)/2:
X=E=>
Y=F=>
S="
S◢Goto0⊿⊿↙
M=X:
N=Y:
动行时,M、N输入第一点的坐标,以后提示X、Y则依次按顺时针或逆时针方向输入各点的坐标,当所有点输入完成后,还需要再次输入第一点的坐标,这样就会自动显示计算结果。
计算结果。
不论顺时针或逆时针方向输入,都有可能导致计算的结果为负,取其绝对值即可。
三维断面面积体积测量CASIO4850
F1A3
L1T=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<
N=>
GOTO1⊿K=K+(C+P)(D-Z)/2◢
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文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 计算 程序