卡西欧万能路线坐标计算程序(值得收藏)
算土方?算坐标?简单容易：https://www.360docs.net/doc/9117927166.html,/[/B]（道路填挖、河塘清淤回填算量、出图等）。
CASIO万能坐标程序
作者（马超）
程序特点：
真正的全线贯通坐标计算，在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，扩充变量数据库，无需修改程序内容；傻瓜操作，适用初级用户。
一、程序：MC
W“1.JS 2.SZ”:W=1=>Z[2]=0:Goto 1 ΔW=2=> O“KOU LING”:O≠123456=>O=0: “OUT” ◢Goto 4ΔO=0: Z[1]=0:Goto 0←┘
Lbi 0←┘
”N0.”：Z[1]+1 ◢←┘
Z[1]=0=>{ABCREFGU}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G” :
Z[Z[1]×8+3]=A:Z[Z[1]×8+4]=B:Z[Z[1]×8+5]=C:Z[Z[1]×8+6]= R-1:Z[Z[1]×8+7]=
E-1:Z[Z[1]×8+8]=F: Z[Z[1]×8+9]=F+G: Z[Z[1]×8+10]=U: “NEXT”◢Isz Z[1]: Goto
0ΔZ[1]=1=>D=Z[9]:Z=0:Z[2]=0:GOTO 2ΔD=Z[(Z[1]-1)×8+9]:Z=0:Z[2]=Z[1]-1:GOTO 2←┘
Lbi A: Z[Z[1]×8+3]=X:Z[Z[1]×8+4]=Y:Z[Z[1]×8+5]=J: Z[Z[1]×8+8]=D: {REGU}:R“Ro”：E“RN”: G“LS”:U“G”: Z[Z[1]×8+6]=R-1 :Z[Z[1]×8+7]=E-1: Z[Z[1]×8+9]=D+G:
Z[Z[1]×8+10]=U: “NEXT”◢Isz Z[1]: Goto 0←┘
Lbi 1←┘
{DZ}:D：Z：Z[2]=0:Goto 2←┘
Lbi 2←┘
Z[2]>Z[1]=>GoTo 4ΔD≤Z[Z[2]×8+9]=>A=Z[Z[2]×8+3]:B=Z[Z[2]×8+4]: C
=Z[Z[2]×8+5]:R=Z[Z[2]×8+6]: E=Z[Z[2]×8+7]: F=Z[Z[2]×8+8]: G=Z[Z[2]×8+9]:
U=Z[Z[2]×8+10]: Goto3ΔIsz Z[2]:Goto 2
Lbi 3←┘
W=2 =>N=8：≠>N=5ΔP＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）÷N：S=90Q÷π：J＝
C+(NPQ+2UR)NS：L=1←┘
X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Cos (C+(LPQ+2UR)LS）,L,1,(N-1))）+ZCos（J+90)←┘
Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Sin (C+(LPQ+2UR)LS）,L,1,(N-1))）+ZSin（J+90):W=2=>GOTO AΔ
Z=0=>“X”:X:Pause 0: “Y” :Y◢Goto 1Δ
Z<0=>“XL”:X:Pause 0: “YL”:Y◢Goto 1Δfx4850 ①
Z>0=>“XR”:X:Pause 0: “YR”:Y ◢Goto 1 ←┘
Z=0=> X “X” ◢Y “Y”◢Goto 1Δ
Z<0=> X “XL” ◢Y “YL”◢Goto 1Δfx4800 ②
Z>0=> X “XR” ◢Y “YR”◢Goto 1 ←┘
Lbi 4←┘
二、说明
a、编制说明
本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。通过对N=？进行修改，可对辛普森公式进行任意等分进行运算。（注：N为不小于2的整数，N越大精度越高，计算速度越慢；N越小精度越低，运算速度越快，一般曲线取N=4就能满足精度要求，在能满足精度的情况下尽量N取小值，已获得最佳运算速度，不要盲目的追求精度）。
b、程序操作说明
程序分为两部分：1.JS为计算，2.SZ为设置。
1、首先用Shift+Defm键对计算器内存变量进行扩充，扩充变量数为8×X+2（X为曲线元段数，变量数视内存情况尽量大些），运行程序，选2进行曲线要素设置，为防止误操作在正确输入口令123456时方可进行设置，否则显示OUT跳到程序尾（LBI4）。在执行程序中将第一段曲线元要素按规律输入到扩充内存变量中（A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G”），以后曲线仅需输入R“R0”、E“RN” :G“LS”:U“G”，其他参数自动计算得出。在显示NEXT时为提示是否输入下一曲线，按EXE继续输入，如不需再输入则退出程序。
2、选1.JS进行坐标运算，D输入桩号，Z输入左右距离（负为左，正为右，0为中）程序首先自动判断其在哪一线元内，并把其线元要素调出进行计算，即可对全线进行坐标计算。
c、变量说明
X0：Y0：F0――――曲线元起点X、Y坐标及起点正切线方位角
R0：RN――――――曲线元起点及终点半径
D0：LS：Q ――――曲线元起点桩号、路线长度及线路左右偏标志（左=-1,直线=0,右=1）
D ――――――――曲线元中待求点桩号
Z ――――――――计算边桩距中线平距，左边输入负值，右边输入正值，中桩输入0
X：Y （中）
XL:YL （左）
XR:YR （右）―――待求点的X，Y坐标
扩充变量：Z[Z[1]×8+3]: Z[Z[1]×8+4]: Z[Z[1]×8+4]: Z[Z[1]×8+6]: Z[Z[1]×8+7]:
Z[Z[1]×8+8]: Z[Z[1]×8+9]: Z[Z[1]×8+10]: 分别为各线元X0:Y0:F0;R0:RN:D0:LS:G,
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0，坐标显示X Y；当位于中线左侧时，Z取负值,坐标显示XL YL,；当位于中线右侧时，Z取正值,坐标显示XR YR。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次方代替。
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径等
于设计规定的值；与圆曲线相接时，曲率半径等于圆曲线的半径。
(7)当两段缓和曲线相连时，相连处要素均输为10的45次方（可以看作缓直+直缓）。
本程序可在4800 4850上运行，只需将输出部分进行修改。（4850① 4800②）[/B]
郑石高速路面NO.9标余官营互通E匝道
序号Xo Yo Fo Ro RN D0 LS G
1 744383.11
2 510406.456 43。34"00.0" 7000 7000 0 170 -1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
85 744445.055 510464.662 170 744507.700 510522.112
2 744507.700 510522.112 42。10"30.8" 7000 500 170 116.071 1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
228.036 744550.193 510561.637 286.071 744590.136 510603.713
3 744590.136 510603.713 49。18"02.3" 500 500 286.071 299.282 1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
435.712 744669.415 510729.969 585.353 744707.945 510873.987
4 744707.94
5 510873.987 83。35"44.8" 500 无穷大585.353 60 1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
615.353 744710.545 510903.872 645.353 744712.247 51093
3.823
CASIO4800、4850坐标正反计算程序
2008年1月27日更新，加入了坐标反算功能，坐标正算加入斜角计算功能。
程序特点：
真正的全线贯通坐标正反计算！在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出，解决了主线坐标计算无法获得第二段及其以后曲线元起点参数的问题；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，参数存储采用扩充变量数据库，无需修改程序内容；多功能采用单程序编程，苊馄捣钡饔米映绦颍 岣咴怂闶涠龋簧倒喜僮鳎 乇鹗视贸跹д摺?
一、程序：MC
W“1.ZS 2.FS 3.SZ”:W=1=>Z[2]=0:V=0：Goto 1 ΔW=2=> Goto 4ΔW=3=>
O“KOU LING”:O≠123456=>O=0: “OUT” ◢Goto DΔO=0: V=0:Z[1]=0:Goto 0←┘
Lbi 0←┘
”N0.”：Z[1]+1 ◢←┘
Z[1]=0=>{ABCREFGU}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G” :
Z[Z[1]×8+3]=A:Z[Z[1]×8+4]=B:Z[Z[1]×8+5]=C:Z[Z[1]×8+6]= R-1:Z[Z[1]×8+7]= E-1:Z[Z[1]×8+8]=F: Z[Z[1]×8+9]=F+G: Z[Z[1]×8+10]=U: “NEXT”◢Isz Z[1]: Goto
0ΔZ[1]=1=>D=Z[9]:Z=0:Z[2]=0:GOTO 2ΔD=Z[(Z[1]-1)×8+9]:Z=0:Z[2]=Z[1]-1:GOTO
2←┘
Lbi A←┘
Z[Z[1]×8+3]=X:Z[Z[1]×8+4]=Y:Z[Z[1]×8+5]=J: Z[Z[1]×8+8]=D: {REGU}:R“Ro”：E“RN”: G“LS”:U“G”: Z[Z[1]×8+6]=R-1 :Z[Z[1]×8+7]=E-1: Z[Z[1]×8+9]=D+G: Z[Z[1]×8+10]=U:
“NEXT”◢Isz Z[1]: Goto 0←┘
Lbi 1←┘
{DZO}:D：Z：O“RJ”：Z[2]=0:Goto 2←┘
Lbi 2←┘
V≠1=>Z[2]>Z[1] =>GoToDΔΔD≤Z[Z[2]×8+9]=> A=Z[Z[2]×8+3]:B=Z[Z[2]×8+4]: C =Z[Z[2]×8+5]:R=Z[Z[2]×8+6]: E=Z[Z[2]×8+7]: F=Z[Z[2]×8+8]: G=Z[Z[2]×8+9]:
U=Z[Z[2]×8+10]: Goto3ΔIsz Z[2]:Goto 2
Lbi 3←┘
W=3 =>N=8：≠>N=4ΔP＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）÷N：S=90Q÷π：J＝
C+(NPQ+2UR)NS：L=1←┘
X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Cos (C+(LPQ+2UR)LS),L,1,(N-1)))+ZCos（J+ O）←┘
Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Sin (C+(LPQ+2UR)LS),L,1,(N-1)))+ZSin（J+ O）:V=1=>Goto6ΔV=2=>Goto9ΔV=3=> GOTO
CΔW=3=>GOTO AΔ
Z=0=>“X”:X:Pause 0: “Y” :Y◢Goto 1Δ
Z<0=>“XL”:X:Pause 0: “YL”:Y◢Goto 1Δfx4850 ①
Z>0=>“XR”:X:Pause 0: “YR”:Y ◢Goto 1 ←┘
Z=0=> X “X” ◢Y “Y”◢Goto 1Δ
Z<0=> X “XL” ◢Y “YL”◢Goto 1Δfx4800 ②
Z>0=> X “XR” ◢Y “YR”◢Goto 1 ←┘
Lbi 4←┘
{MH} ：M“X” ：H“Y” ：Z[2]=0：GOTO 5←┘
Lbi 5←┘
V=1：D= Z[Z[2]×8+9]：Z=0 ：O=90：GOTO 2←┘
Lbi 6←┘
K=（（H -B）Cos（C-90）-（M-A）Sin（C-90））×（（H -Y）Cos（J-90）-（M-X）Sin（J-90））：
K≤0=> Goto 7ΔIsz Z[2]：Goto5←┘
Lbi 7←┘
D=F+Abs((H -B)Cos(C-90)-(M-A)Sin（C-90）)：D>G=> Isz Z[2]: Goto5ΔGoto 8←┘
Lbi 8←┘
V=2 ：GOTO 3←┘
Lbi 9 ←┘
K=(H -Y)Cos(J-90)-(M-X)Sin(J-90)：AbsK
Lbi B←┘
V=3 ：Z=0 ：Goto 3←┘
Lbi C←┘
Z=(H-Y) ÷Sin(J+90)：“D”:D:Pause 0: “Z” Z◢
(Z=(H-Y) ÷Sin(J+90)：D“D” ◢“Z” Z◢)4800输出
GOTO 4←┘
Lbi D←┘
二、说明
a、编制说明
本程序是运用复化辛普生公式根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距，对该
曲线段范围内任意里程中边桩坐标进行计算，以及对卡西欧扩充变量的灵活应用，实现了真正意义上的的全线贯通及曲线要素输入程序化（在不修改程序内容的情况下可通过运行程序输入任意多段曲线元要素）。通过对N=？进行修改，可对辛普森公式进行任意等分进行运算。（注：N为不小于2的整数，N越大精度越高，计算速度越慢；N越小精度越低，运算速度越快，一般曲线取N=4就能满足精度要求，在能满足精度的情况下尽量N取小值，已获得最佳运算速度，不要盲目的追求精度）
b、程序操作说明
程序分为两部分：1.JS为计算，2.SZ为设置。
1、首先用Shift+Defm键对计算器内存变量进行扩充，扩充变量数为8×X+2（X为曲线元段数，变量数视内存情况尽量大些），运行程序，选2进行曲线要素设置，为防止误操作在正确输入口令123456时方可进行设置，否则显示OUT跳到程序尾（LBI4）。在执行程序中将第一段曲线元要素按规律输入到扩充内存变量中（A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G”），以后曲线仅需输入R“R0”、E“RN” :G“LS”:U“G”，其他参数自动计算得出。在显示NEXT时为提示是否输入下一曲线，按EXE继续
输入，如不需再输入则退出程序。
2、选1.JS进行坐标运算，D输入桩号，Z输入左右距离（负为左，正为右，0为中）程序首先自动判断
其在哪一线元内，并把其线元要素调出进行计算，即可对全线进行坐标计算。
c、变量说明
X0：Y0：F0――――曲线元起点X、Y坐标及起点正切线方位角
R0：RN――――――曲线元起点及终点半径
D0：LS：Q―――曲线元起点桩号、路线长度及线路左右偏标志（左=-1,直线=0,右=1）
D――――――――曲线元中待求点桩号
Z―――――――――计算边桩距中线平距，左边输入负值，右边输入正值，中桩输入0
RJ ：坐标正算右夹角
X：Y （中）
XL:YL （左）
XR:YR （右）―――――――待求点的X，Y坐标
扩充变量：Z[Z[1]×8+3]: Z[Z[1]×8+4]: Z[Z[1]×8+4]: Z[Z[1]×8+6]: Z[Z[1]×8+7]: Z[Z[1]×8+8]: Z[Z[1]×8+9]: Z[Z[1]×8+10]: 分别为各线元X0:Y0:F0;R0:RN:D0:LS:G,
(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，
Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
(2) 当所求点位于中线时，Z=0，坐标显示X Y；当位于中线左侧时，Z取负值,坐标显示XL YL,；
当位于中线右侧时，Z取正值,坐标显示XR YR。
(3) 当线元为直线时，其起点、止点的曲率半径为无穷大，以10的45次方代替。（对于特别长的
直线应适当提高10的次方数，以防止以大半径曲线误算，建议以10的90次方计算）
(4) 当线元为圆曲线时，无论其起点、止点与什么线元相接，其曲率半径均等于圆弧的半径。
(5) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。止点与直线相接时，曲率半径为无穷大，以10的45次方代替；与圆曲线相接时，曲率半径等于圆曲线的半径。
(6) 当线元为非完整缓和曲线时，起点与直线相接时，曲率半径等于设计规定的值；与圆曲线相接
时，
本程序可在4800 4850上运行，只需将输出部分进行修改。（4850① 4800②）
郑石高速路面NO.9标余官营互通E匝道
序号Xo Yo Fo Ro RN D0 LS G
1 744383.11
2 510406.456 43。34"00.0" 7000 7000 0 170 -1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
85 744445.055 510464.662 170 744507.700 510522.112
2 744507.700 510522.112 42。10"30.8" 7000 500 170 116.071 1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
228.036 744550.193 510561.637 286.071 744590.136 510603.713
3 744590.136 510603.713 49。18"02.3" 500 500 286.071 299.282 1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
435.712 744669.415 510729.969 585.353 744707.945 510873.987
4 744707.94
5 510873.987 83。35"44.8" 500 无穷大585.353 60 1
中点桩号中点X 中点Y 终点桩号终点X 终点Y
615.353 744710.545 510903.872 645.353 744712.247 510933.823
计算时间测试（秒），计算精度满足1mm时测试每段曲线元的终点
4等分
正算：2.3 2.4 2.5 2.6
反算：9.6 12.4 16.5 17.0
5等分
正算：2.5 2.7 2.8 2.9
反算：11.8 13.8 18.7 19.1
CASIO4800、4850标高程序
程序特点：
全线标高数据一次性程序化输入，扩充变量数据库，无需修改程序内容；全线贯通计算，标高计算时将全线数据分为4等分，首先判断待求点在哪个区间，再循环选择，提高运算时间；傻瓜操作，适用初级用户。
一、程序：MC
W“1.BG 2.SZ ”：W=1=> Goto 1ΔO“KOU LING”：O≠123456=>O=0：Goto ５ΔO=0：Z[1]=0:
Goto 0←┘
Lbi 0←┘
”No.”：Z[1]+1◢
{ZHPR}：Z“D0” ：H“BG” ：P“PD” ：R“R0” ：Z[Z[1]×4+3]=Z：Z[Z[1]×4+4]=H：Z[Z[1]×4+5]=P：Z[Z[1]×4+6]=R：“NEXT”◢Isz Z[1]：Goto 0←┘
Lbi 1 ←┘
Z[2]= Z[1]←┘
{D}：D≤Z÷4+1)×4+3]=>Z[2]= Int(Z [2]÷4) ：Goto 2ΔD≤Z
÷2+1)×4+3]=>Z[2]=Int(Z[2]÷2)：Goto 2ΔD≤Z×.75+1)×4+3]=>Z[2]=Int (Z[2]×.75)：
Goto 2ΔGoto 2←┘
Lbi 2 ←┘
D
Lbi 3←┘
Z[Z[2]×4+6]≠0=>G=Z[Z[2]×4+4]+(D-Z[Z[2]×4+3])×Z[Z[2]×4+5]÷100+(D-Z[Z[2]×4
+3]) ^2÷2÷Z[Z[2]×4+6]：Goto
4ΔG=Z[Z[2]×4+4]+(D-Z[Z[2]×4+3])×Z[Z[2]×4+5]÷100：Goto 4←┘
Lbi 4 ←┘
{Z}：Z：G=G-0.02Abs(Z):“G”: G◢（4850输出）
{Z}：Z：G=G-0.02Abs(Z): G“G”◢（4800输出）
Goto 1 ←┘
Lbi 5 ←┘
“OUT” ◢
程序操作说明
程序分为两部分：1.BG为标高计算，2.SZ为参数设置。
1、首先用Shift+Defm键对计算器内存变量进行扩充，扩充变量数为4×X+2（X为标高参数段数，变量数视内存情况尽量大些），运行程序，选2进行标高参数设置，为防止误操作在正确输入口令123456时方可进行设置，一组参数输入完毕后显示“NEXT”，如输入完毕请退出，按“ＥＸＥ”继续。
2、选1.ＢＧ进行标高运算，D输入桩号，Z输入距离，程序首先自动判断其待求点在哪一段参数内，并
将参数调出进行计算，即可对全线进行标高计算。
c、变量说明
Z“D0”―――――――起点桩号（非变坡点）
H“BG”―――――――起点标高（非变坡点）
P“PD”―――――――前坡度（％）
R“R0”―――――――竖曲线半径（凸曲线为负值；凹曲线为正值，直坡为０）
公路测量卡西欧5800万能程序
一、前言本程序是《CASIO fx-5800P计算与道路坐标放样计算》中道路坐标放样计算程序的升级改进版本。原道路坐标放样计算程序只基于道路的单个基本型曲线，有效计算范围仅包括平曲线部分和前后的两条直线段，使用时需要输入平曲线设计参数，无坐标反算桩号功能。改进后的程序名称为：道路中边桩坐标放样正反算程序（全线贯通），增加了可实现全线贯通的数据库功能和坐标反算桩号功能，主要是： 1．使用道路平面数据库子程序，可将一段或若干段道路的交点法格式平面参数（可容易从直线、曲线及转角表中获得）以数据库子程序形式输入计算器，程序在计算时省却了输入原始数据的麻烦； 2．坐标正算方面，输入桩号即可进行道路的中、边桩坐标计算，若输入了测站坐标，还可同时计算全站仪极坐标放样数据（拨角和平距）； 3．坐标反算方面，输入平面坐标，即可计算对应的桩号和距中距离（含左右信息）； 4．对于存在断链的道路，可分段分别编写数据库子程序，然后在主程序中添加一个路段选择的功能即可实现（可参照立交匝道程序中匝道的选择）。程序的特点： 1．可进行中桩坐标的正、反算，程序代码简洁，便于阅读和改写； 2．主程序通过调用数据库子程序，省却了使用时输入平面参数的繁琐； 3．使用数据库子程序，换项目只需改写数据库子程序，程序通用性强。二、道路示例项目基本资料基本资料同《CASIO fx-5800P计算与道路坐标放样计算》第6章HY高速公路第2合同段（合同段起止桩号： K4+800~K9+600）。这里摘取直线、曲线及转角表资料如下.
.
. 三、程序代码 .
. .
. .
5800计算器程序下载
阅览室馆友我的图书馆 kaixin100 | | 分享
U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos( G+QEFW (C+FWD))+Acos(G+QEMW(C+MWD)))→X： V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+ QEFW(C+FWD))+Asin(G+QEMW(C+MWD))) →Y： G+QEW(C+WD)+90→F：X+Zcos（F）→X：Y+Zsin（F）→Y 反算子程序SUB2 G-90→T (Y-V)cosT-(X-U)sin(T) →W Abs(W)→W:0→Z Lbl6：Prog "SUB1" T+QEW(C+WD) →L：(J-Y)cos（L）-(I-X)sin（L）→Z IF Abs（Z）<1E-6：Then0→Z:Prog "SUB1":(J-Y)÷sin（F）→Z:Else W+Z→W：Goto6:IfEnd 数据库子程序SUB0 Goto 1(线元可输入多条，分离式可在前多加一位，匝道一样。例：左幅为K129+500,右幅输线元参数里程为1129+500，其他不变，前面 1为任意数字，计算机便于区分) Lbl 1：IF S<线元终点里程：Then@@@→O(线元起点里程) ：@@@ →U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G(线元起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEnd IF S<线元终点里程：Then@@@→O(线元起点里程) ：@@@ →U(线元起点X坐标)：@@@→V(线元起点Y坐标)：@@@→G线元(起点计算方位角)：@@@→P(线元起点半径)：@@@→R(线元止点半径)：@@@→H(线元长度)：@@@→Q(线元左、右偏标志，左偏-1，右偏1，直线为0)：Return：IfEnd 。。。。。。。。。。。。。。。。。。。。。。。。。。。 一程序功能 本程序由一个主程序(ZBJS)和3个子程——正算子程序(SUB1)、反 算子程序( SUB2) 、数据库子程序（SUB0）构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线 元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲 率半径、止点曲 率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标 进行正反算。另 外也可以将本程序中核心算法部分的两个子程序移植到其它相关的 程序中，用于对曲 线任意里程中边桩坐标进行正反算。本程序也可以在CASIO fx-4500P计算器及CASIO fx-4850P计算器上运行。 二、使用说明
5800计算器全线坐标计算放样程序(修改第三版)
5800计算器全线坐标计算放样程序（修改第三版） 5800计算器全线坐标计算放样程序（修改版） “XLZBJSCX” ◢ LB1 0 ↙ CLS : FIX 4 : 30→DIM Z ↙ “XHS="?G(后视点X):"YHS="?L(后视点Y):"XZJ="?M(置镜点X):"YZJ="?N(置镜点Y)0l(G-M,L-N):"DH=":I(后视距)◢J<0=>J+360→J:"FH=":J→DMS◢(后视 方位角) LB1 1 ↙ “K=”?K ◢（计算里程） IF K<本段曲线终点里程 AND K≥上段曲线终点里程：THEN 本段终点里程→Z[1] : 上段曲线终点里程→Z[2] ：1→0 （注：左偏曲线输入-1→0,右偏曲线输入1→0）: 偏角→A：半径→R : 第一缓和曲线→Z[6] : 第二缓和曲线 →Z[7] : 交点X→B :交点Y→C : 小里程向交点方位角→E : 交点向大里程方位角→F : G0T0 2 : IFEND↙ …………（曲线段分段输入） 补充直线段输入如下 IF K<本段直线终点里程 AND K≥本段直线起点里程:THEN 1→0:本段直线终点里程→Z[2]:终点X→Z[16]:终点Y→Z[11]:方位角→E:G0T0 4:IFEND LB1 2 ↙（曲线要素计算） Z[6]/2- Z[6]^3/(240*R^2)+ Z[6]^5/(34560*R^4) →Z[8] ↙（M1） Z[7]/2- Z[6]^3/(240*R^2)+ Z[7]^5/(34560*R^4) →Z[9] ↙（M2） Z[6]^2/(24*R)- Z[6]^4/(2688*R^3) →Z[10] ↙（P1） Z[7]^2/(24*R)- Z[7]^4/(2688*R^3) →Z[11] ↙（P2） π*A*R/180+0.5*( Z[3]+ Z[2])→W ↙（曲线总长） 90* Z[6]/(R*π) →Z[14] ↙（第一缓和曲线总偏角） 90* Z[7]/(R*π) →Z[15] ↙（第二缓和曲线总偏角,可以省略） Z[8]＋（R+Z[10])TAN(A/2)-(Z[11]-Z[11] )/SIN A→Z[11]↙ (切线T1) Z[9]＋（R+Z[12])TAN(A/2)+(Z[10]-Z[12] )/SIN A→Z[12]↙ (切线T2) B+ Z[12]*C0S (E+180)→ Z[13] ↙（ZH点X） C+ Z[12]*SIN(E+180)→ Z[15] ↙（ZH点Y） Z[1]-S→Z[3] ↙ (ZH点里程) Z[3]+ Z[6]→Z[4] ↙ (HY点里程) Z[1]- Z[7]→Z[5] ↙ (YH点里程) G0T0 3 ↙ LB1 3 ↙(判断里程点与曲线关系) IF K≤Z[3] AND K> Z[2] : THEN G0T0 4 : IFEND ↙ IF K≤Z[4] AND K> Z[3] : THEN G0T0 5 : IFEND ↙ IF K≤Z[5] AND K> Z[4] : THEN G0T0 6 : IFEND ↙
CASIO fx-5800P实用工程测量程序
一、QXFY 辛甫森公式放样程序 1. “X0”? U：“Y0”?V 2. “XA”? A：“Y A”? B：“CA”? C：“1÷RA”?D：“1÷RB”?E： “KA”?F：“KB”? G 3. Lb1 1：“KI”?H：“JJ”?L：“Y+Z－”?R 4. If H>G Or H 5800全线任意坐标计算程序 1. 正算主程序(ZHCX) （不运行） 8→DimZ 1÷P→Z[4 ]：(P-R)÷(2HPR)→D: 180÷π→E “Z=”?Z:”YJJ=”?A:Abs(S-O)→W 0.26→Z[1 ]: 0.74→B: 0.02→K: 0.82→Z[3 ]: 1-Z[3 ]→F:1-K→Z[2 ] U+W(Z[1 ]cos(G+QEKW(Z[4 ]+KWD))+Bcos(G+Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+Bcos(G+QEFW (Z[4 ]+FWD))+ Z[1 ]cos(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→X: V+W(Z[1 ] sin (G+QEKW(Z[4 ]+KWD))+B sin(G+ Z[3 ]QEW(Z[4 ]+ Z[3 ]WD))+B sin(G+QEFW (Z[4 ]+FWD))+ Z[1 ] sin(G+ Z[2 ]QEW(Z[4 ]+ Z[2 ]WD)))→Y： G+QEW(Z[4 ]+WD)→F：X+Zcos(F+A)→X：Y+Zsin(F+A)→Y:If F≧360:Then F-360→F:IfEnd ”X=”:X→X◢ ”Y=”:Y→Y◢ If F﹤0:Then F+360→F:IfEnd ”QX FWJ=”:F▼DMS◢ “C=1=>XX: C=2=>XZ”: ”C=”?C: ”QHJU=”?L: If C=1：Then Goto 1:Else Goto 2: IfEnd 可以计算斜交斜做或斜交正做的桥涵坐标 Lbi 1 X+L cos(F)→X：Y+Lsin(F)→Y: Goto 3 Lbi 2 X+L cos(F+A-90)→X：Y+Lsin(F+A-90)→Y: Goto 3 Lbi 3 “QH-X=”: X →X◢ “QH-Y=”: Y →Y◢ Prog “FY” 2 . 参数子程序（直接运行） M(主线) 一条线路一个名称 “S=”?S If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd … … If S≦线元终点：Then 线元起点X值→U: 线元起点Y值→V:线元起点切线方位角→G:线元起点桩号→O:线元长度→H：线元起点半径→P：线元终点半径→R：（左偏-1，或右偏 1）→Q:Goto 1:IfEnd Lbi 1 Prog “ZBJS” 3. 放样程序（FY）（不运行） “X0=”?M:“Y0=”？N Pol((X-M, Y-N) 线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器) 论文https://www.360docs.net/doc/9117927166.html,/：本论文仅供学习交流使用，本站仅作合理转载，原作者可来邮要求删除论 文。 摘要：我国公路建设事业正处于一个高速发展的时期，在公路工程施工过程中，施工技术人员经常要使用全站仪、水准仪进行施工放样、高程测量，在测量过程中，手工计算速度慢，失误率高，工作效率极低。利用CASIO fx-9750GⅡ编程函数计算器强大的内存（可诸存63000个字符）和编程功能，编写各种计算程序，能够在2秒钟内计算出施工放样、桩点坐标等施工过程中的各项数据资料，同时也使我们有更多的时间去挑战更富有创造性的工作。 关键词：坐标放线线元测量程序 1、前言 本程序采用Gauss-Legendre（高斯-勒让德）五节点公式作内核,计算速度（太约2秒）适中，计算精度很高。在此之前，本人曾用过以下公式作内核：①积分公式simpson法②双重循环复化高斯2节点③高斯-勒让德3节点④求和公式复化simpson法⑤双重循环复化simpson法⑥高斯-勒让德4节点，⑦高斯-勒让德5节点，经过测试③计算最快，⑦代码稍长但计算速度只比③⑥稍慢，精度最高，可满足线元长小于1/2πD 的所有线形的精度要求。⑦作内核分别计算圆曲线长1/4πD、1/2πD、3/4πD、πD处的精度,1/4πD时偏差为0.001mm,1/2πD时偏差为0.55m m,3/4πD时偏差为31.63mm，πD时偏差为968mm，偏差按半径倍数增大，如线元长大于1/2πD（1/2圆周长）时，可将其拆分二个或多个线元单位，以确计算保精度。 2、程序特点 事先将所有的平曲线交点的线元要素诸存到计算器内，测量时只输桩号、边距等程序会自动寻找各类要素，一气呵成地完成施工测量任务，中途不需人工转换各类要素数据，本程序可诸存几百条线路的要素数据，计算时可按需选择线路编号进行测量。测量时不需查阅及携带图纸，仅一台CASIO fx-9750GⅡ编程函数计算器即可。 本程序含一个主程序:3XYF，五个子程序:GL（公式内核）、QD（线路选择）、XL（线路要素判断）、GF（坐标反算）、File 1 （要素存放的串列工作簿）。可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、终点里程、起点曲率半径、止点曲率半径）及里程边距或坐标，对该线元段范围内任意里程中边桩坐标进行正反算。 3、计算公式及原理 如图：BC 间为一曲线元，曲线元上任一点的曲率随至B 点的弧长作线性变化。设起点B 的曲率为KA ，终点C 的曲率为KB ，R 为曲线半径。±表示曲线元的偏向，当曲线元左偏时取负号，当曲线元右偏时取正号,直线段以1的45次方代替（即半径无穷大）。 式中：αΑ＝起始方位角l ＝p 点到B的距离lS=曲线总长αp＝p 点切线方位角 R1=R5=0.118463442528095 ,R2 = R4 = 0.239314335249683 , R3 = 0.28444444444444 V1=1-V5= 0.046910070 ,V 2= 1-V4 = 1 0.2307653449 V3= 0.5 利用上面公式及CASIO fx-9750GⅡ编程函数计算器可编写下列计算程序。 4、程序清单 （1）、3XYF（主程序） "1→XY2→FS"?→V:V=1=>Goto 1:V=2=>Goto 2↙（选择计算功能） Lbl 1:File 1:”XLn”?→S:Prog “QD”↙（选择线路） 5800计算器公路全线坐标正、反算计算程序FX5800全线贯通万能正、反算程序（一体化、超好用、短小、易懂） FX5800计算器的积分程序（正反算、全线贯通、新线路）终极版 ZHUCHENGXU主程序 "1.ZS,2.FS" ?→Q←┘输入1正算，输入2反算 “NEW=0,OLD≠0”?Z←┘ IfZ=0:Then“X0=”?A:“Y0=”?B:“C0=”?C:“1/R0=”?D:“1/RI=”?E:“SP=”?F:“EP=”?G:Ifend:Q=2=>Goto 2←┘ Lbl1 :“KM=,<0 Stop”?H:H<0=>Stop:“PJ=”?O:“PY=”?L←┘ LblZ:Z=1=> Prog“01”：Z=2=> Prog“02”←┘选择数据库文件，可增加 H- F→X:0.5(E-D)÷(G-F)→N←┘ C+(XD+NX2)*180÷π→P：P<0=>P+360→P：P>360=>P-360→P←┘- A+∫(cos(C+(XD+NX2)*180÷π),0,X)+Lcos(P+O)→U←┘ B+∫(sin(C+(XD+NX2)*180÷π),0,X)+Lsin(P+O)→V←┘ Q=2=>Goto4：Cls：Fix 3←┘ "Xn=":Locate4，1，U："Yn=": Locate5，2，V：“FWJ=”：PDMS◢ Norm 2：Cls：Goto1←┘ Lbl2:“XD=,<0,STOP”？R：R<0=>Stop:“YD=”？S←┘ “KMDG=”？H ：90→O:0→L：GotoZ←┘（H线路范围内的任意桩号） Lbl4:Pol（R-U，S-V）：J<0 => J+360→J←┘ Whileabs（Icos（J-P））≤0.001:P-J>180=> J+360→J： P-J<-180=> P+360→P:If P-J>0:then -I→L:else I→L :Ifend:Goto3: Whileend：H+Icos（J-P）→H:GotoZ←┘ Lbl3:Cls：Fix 3←┘ 程序：MC W“1.JS 2.SZ”:W=1=>Z[2]=0:Goto 1 ΔW=2=> O“KOU LING”:O≠123456=>O= 0: “OUT” ◢Goto 5ΔO=0:Defm 42:Z[1]=0:Goto 0←┘ Lbi 0←┘ {ABCREFGU}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G” ←┘ Z[Z[1]*8+3]=A:Z[Z[1]*8+4]=B:Z[Z[1]*8+5]=C:Z[Z[1]*8+6]=1÷R:Z[Z[1]*8+7]=1÷E:Z [Z[1]*8+8]=F: Z[Z[1]*8+9]=F+G: Z[Z[1]*8+10]=U←┘ A=0=> Z[2]=0:Goto 1ΔIsz Z[1]: Goto 0←┘ Lbi 1←┘ {DZ}:D：Z：Z[2]=0:Goto 2←┘ Lbi 2←┘ D≤Z[Z[2]*8+9]=>A=Z[Z[2]*8+3]:B=Z[Z[2]*8+4]: C =Z[Z[2]*8+5]:R=Z[2]*8+6]: E= Z[Z[2]*8+7]: F=Z[Z[2]*8+8]: G=Z[2]*8+9]: U=Z[Z[2]*8+10]: Goto3:ΔIsz Z[2]:Got o 2 Lbi 3←┘ P＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）：I＝PQ：J＝C+90 Q（I+2UR）/π：J<0=>J=J+360ΔM＝C+45 Q（I÷4+2UR）÷2π：N＝C+135Q（3I÷4+2UR）÷2π：K＝C+45Q（I÷2+ 2UR）÷π:L=C+45Q(I÷8+2UR) ÷4π: S=C+135Q(3I÷8+2 UR) ÷4π: T=C+225Q(5I÷8+2UR) ÷4π: H=C+315Q(7I÷8+2UR) ÷4π←┘ Lbi 4 X＝A+Q(Cos C+4（Cos L+ Cos S +Cos T+Cos H）+2(Cos M + Cos N+Cos K)+Cos J) ÷24+ZCos（J+90)←┘ Y＝B+Q(Sin C+4（Sin L + Sin S + Sin T+ Sin H）+2(Sin M + Sin N+ Sin K) + Sin J) ÷24 +ZSin（J+90)←┘ Z=0=>“X”:X:Pause 0: “Y” :Y◢Goto 1Δ Z<0=>“XL”:X:Pause 0: “YL”:Y◢Goto 1Δfx4850 Z>0=>“XR”:X:Pause 0: “YR”:Y ◢Goto 1 ←┘ 卡西欧FX5800---辛普森公式(万能公式) 复化辛普森公式 1.Lbl 0：“XA=”?A：“YA=”?B: “CA=”?C：“1/RA=”?D：“1/RB=”?E：“DKA=”?F:“DKB=”?G 2.Lbl 1:“DKI=”?H：“DL=”?O：“DR=”?R：IF H＞G:THEN Goto0 IFEND 3.(E-D)/Abs(G-F)→P:Abs(H-F) →Q: P×Q→I:D+I→T 4.C+(I+2D)×Q×90/π→J 5. C+(I/4+2D)Q×45/(2π) →M: C+(3I/4+2D)Q×135/(2π) →N 6. C+(I/2+2D)Q×45/π→K 7. A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12 →X 8.B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12 →Y 9.”FW=”:J▲DMS ▲ 10. “X=”:X▲ 11. “Y=”:Y▲ 12.“XL=”: X+Ocos(J-90) →U ▲ 13.“YL=”: Y+Osin(J-90) →V▲ 14.“XR=”:X+Rcos(J+90)→ W▲ 15.“YR=”:Y+Rsin(J+90)→ Z▲ 16.Goto 1 程序结束 程序说明： A- 曲线元起点A的坐标； B- 曲线元起点B的坐标； C- 曲线元起点A的切线坐标方位角； F- 曲线元起点A的里程； G- 曲线元起点B的里程； H- 曲线上待求点i的里程； D- 曲线元起点A的曲率； E- 曲线元终点B的曲率； XL-左边线点位X坐标； YL-左边线点位Y坐标； XR-右边线点位X坐标； YR-右边线点位Y坐标； X- 中线点位纵坐标； Y- 中线点位横坐标； DL-左边线距中线平距； DR-右边线距中线平距； 该程序需要输入的数据为： （1）.曲线元起点A的坐标及切线坐标方位角，计算器上用“XA”，“YA”，“CA”显示； （2）.曲线元起点A和B的曲率，计算器上用I÷RA，I÷RB显示（曲线左偏时取“-”）； （3）.曲线元起点A和终点B的里程，计算器上用“DKA”，“DKB”显示；（4）.输入待求点里程和该点距左右的水平距离，计算器上用“DKI”，“DL”，“DR”显示； 每算完一个待求点的中线及边线坐标，程序又让输入下一点的“DKI”，“DL”，“DR” 当输入的“DKI”大于“DKB”时，此时输入下一个曲线元起点的曲率和里程，即可计算下一个曲线中线及边线点位坐标。 使用该程序应注意事项； 该程序以前进方向为有意识，不可倒退计算；缓和段和圆曲线段应分开计算在计算圆曲线时应记下缓和 曲线尾的坐标方位角即“J”的角度；在计算第一段缓和曲线时曲率“1÷RA”输入0；在计算第二段缓和曲 线时“1÷RB” 输入0 。 （5）第一个0为零。 1线路中线和边线点位坐标的万能通用程序 Lb1?: { E G }:A“X?”:B“Y?”:C“FWJ”:D“1/R-QD”: E“1/R-ZD”:F“QD-ZH”:G“ZD-ZH”： Lb1 1:{ H O }:H“JS-ZH”:O“JZ-JL”:H > G Goto 2 Δ Prog 1: X“X-JS”=X+Ocos(J+9?)▲Y“Y-JS”=Y+Osin(J+9?)▲Goto1Δ Lb1 2:L=H:H=G:“WARING>JS-FW”▲“SR…SUJU”▲ Prog 1: A=X:B=Y:D=E:F=G:C=J:H=L: Got o ?Δ 子程序1: P= (E-D)/ABS(G-F):Q=ABS(H-F):I=PQ： J=C+9?Q(I+2D)/π： M=C+45Q(I/8+D)/π: N=C+135Q(3I/8+D)/π K=C+45Q(I/2+2D)/π X=A+Q(cosC+4(cosM+cosN)+2cosK+cosJ)/12: Y=B+Q(sinC+4(sinM+sinN)+2sinK+sinJ)/12: 注释：A“X?”:B“Y?”:-起点的xy坐标 C“FWJ”:-起点的方位角 D“1/R-QD”:起点的 1/半径，如直线为1/e50 E“1/R-ZD”:终点的 1/半径，如半径100为1/100（如左转为负，右转为正）“QD-ZH”:G“ZD-ZH”：起点、终点的桩号 H“JS-ZH”:O“JZ-JL”:输入计算的桩号，距中距离、如大于计算范围则显示 “WARING>JS-FW”▲“SR…SUJU”▲只需继续输入下一点的终点的 1/半径和终点桩号即可连续计算下一线型。 X“X-JS”、y“y-JS”显示计算点的xy坐标 卡西欧5800道路坐标测量程序（通俗易懂版） 主程序（自己给取个名吧） 2→Dimz “1.SZ＝﹥XY 2.XY＝﹥SZ”◢ “N=”?N If N=1:Then Goto 1:Else Goto 2:IfEnd Lbl 1 “S=”?S “Z=”?Z “J=”?T If S﹤0:Then Prog “SYL”:Else Prog “SYR”:IfEnd(可拿掉) 1÷P→C:(P-R) ÷(2HPR)→D:180÷∏→E:Abs(Abs(S)-O)→W:Prog “SUB1” “QJ=”:Z[1]◣DMS◢ “X=”:X◢ “Y=”:Y◢ Goto 1 Lbl 2 “S=”?S If S<0:Then Prog “SYL”:Else Prog “SYR”:IfEnd(可拿掉) 90→T 1÷P→C:(P-R) ÷(2HPR)→D:180÷∏→E “X=”?X “Y=”?Y X→I:Y→J Prog “SUB2” O+W→S（红色加粗加下画线为字母，常规为零） “S=”:S◢ Z→Z “Z=”:Z◢ Goto 2 子程序（SUB1） 0.1739274226→A：0.3260725774→B:0.0694318442→K:0.3300094782→L:1-L→F:1-K→M:U+W(Acos(G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Bcos(G+QEFW(C +FWD))+Acos(G+QEMW(C+MWD)))→X V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+Bsin(G+QEFW(C+F WD))+Asin(G+QEMW(C+MWD))) →Y G+QEW(C+WD)→Z[1] 万能坐标计算公式 X=起点x+(待求点桩号-起点桩号)*【0.1184634425*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.046910077+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.046910077^2)+0.2393143352*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.2307653449+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.2307653449^2)+0.2844444444*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.5+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.5^2)+0.2393143352*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.7692346551+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.7692346551^2)+0.1184634425*COS(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.953089923+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.953089923^2))+边距*COS(偏角+起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)】 Y=起点Y+(待求点桩号-起点桩号)*【0.1184634425*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.046910077+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.046910077^2)+0.2393143352*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.2307653449+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.2307653449^2)+0.2844444444*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.5+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.5^2)+0.2393143352*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.7692346551+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.7692346551^2)+0.1184634425*SIN(起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)*0.953089923+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)*0.953089923^2))+边距*SIN(偏角+起始方位角/180*3.14159265+起点曲率*(待求点桩号-起点桩号)+0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号)*(待求点桩号-起点桩号)】 切线方位角A=起始方位角+(终点曲率*(待求点桩号-起点桩号) +0.5*(终点曲率-起点曲率)/(终点桩号-起点桩号)*(待求点桩号-起点桩号) ^2)*180/3.14159265 如何设置手持GPS相关参数及全国各地坐标转换参数一、如何设置手持GPS相关参数 （一）手持GPS的主要功能 手持GPS，指全球移动定位系统，是以移动互联网为支撑、以GPS智能手机为终端的GIS系统，是继桌面Gis、WebGis之后又一新的技术热点。目前功能最强的手持GPS，其集成GPRS通讯、蓝牙技术、数码相机、麦克风、海量数据存储、USB/RS232端口于一身，能全面满足您的使用需求。 主要功能：移动GIS数据采集、野外制图、航点存储坐标、计算长度、面积角度（测量经纬度，海拔高度）等各种野外数据测量；有些具有双坐标系一键转换功能；有些内置全国交通详图，配各地区地理详图，详细至乡镇村落，可升级细化。 （二）手持GPS的技术参数 因为GPS卫星星历是以WGS84大地坐标系为根据建立的，手持GPS单点定位的坐标属于WGS84大地坐标系。WGS84坐标系所采用的椭球基本常数为：地球长半轴a=6378137m；扁率F=1／298.257223563。 常用的北京54、西安80及国家2000公里网坐标系，属于平面高斯投影坐标系统。北京54坐标系，采用的参考椭球是克拉索夫斯基椭球，该椭球的参数为：地球长半轴a=6378245m；扁率F=1／298.2。西安80坐标系，其椭球的参数为：地球长半轴a=6378140m；扁率F=1／298.257。国家2000坐标系，其椭球的参数为：地球长半轴a=6378137m；扁率F=1／298.298.257222101。 （三）手持GPS的参数设置 要想测量点位的北京54、西安80及国家2000公里网高精度坐标数据，必须学习坐标转换的基础知识，并分别科学设置手持GPS的各项参数。 首先，在手持式GPS接收机应用的区域内(该区域不宜过大)，从当地测绘部门收集1至两个已知点的北京54、西安80或国家2000坐标系统的坐标值；然后在对应的点位上读取WGS84坐标系的坐标值；之后采用《万能坐标转换》软件，可计算出DX、DY、DZ的值。 将计算出的DX、DY、DZ三个参数与DA、DF、中央经线、投影比例、东西偏差、南北偏差等六个常数值输入GPS接收机。将GPS接收机的网格转换为“UserGrid”格式，实际测量已知点的公里网纵、横坐标值，并与对应的公里网纵、横坐标已知值进行比较，二者相差较大时要重新计算或查找出现问题的原因。详细过程可查看《万能坐标转换》软件的【手持GPS参数设置】界面。 （四）自定义坐标系统（User）投影参数的确定 1、自己观测计算 新机拿到手之后，供应商都给提供一个投影参数，这对于要求不高的一般用户来说基本可以满足工作需要，而对于一些专业用户来说，就要自己来测算参数。一般型号的导航型手持GPS自定义坐标系统（User）投影参数设置界面都提供了五个变量（△X、△Y、△Z、△A、△F）需要设置，而实际工作中，后两个参数（△A、△F）针对某一坐标系统来说为固定参数（北京54坐标系△A=-108、△ F=0.0000005），无需改动，需要自己测算的参数主要为前三个（△X、△Y、△Z）， 线路任意点坐标计算、及任意点对应桩号，左右偏距计算程序卡西欧4800、4850系列计算器测量计算程序 一、字母含义； K: 表示拟计算位置的线路桩号; H: 表示计算位置距路线中心线的偏距，左偏为正，右偏为负。 T: 各段线路上作为起算点处的切线方位角。 S: 拟计算点到起算点的曲线长。 L: 在圆曲线上表示曲线长，在缓和曲线上表示缓和段长度，在直线上为零。 R: 表示曲线半径，左偏为正，右偏为负 E、F: 起算点的坐标值。 M”X1”N”Y1”: 已知点坐标，求其对应位置桩号及左右偏距。 二、程序 1、坐标计算(COORD) {K,H}:KH“L+，R-”:Prog “DATA”:”(X,Y)=”: X=X+HSinW :Pause 1:Y=Y-HcosW: 2、坐标反算线路桩号(FS ZH) Fix 4:M”X1”N”Y1”:LbI 1:Prog “DATA”:PoI (M-X,N-Y):Fixm:J<0=>J=J+360:≠>J=J⊿ Abs(Sin(W-J ))=1=>”K=”:K: Pause 1 : “L+,R- =”: H=ISin(W-J):≠>K=K+Icos(W-J): Goto 1 3、子程序Prog “DATA” K<（第１段与第２段线路分界处的路线桩号）＝> T=(第１段起算点处的切线方位角值):Ｓ＝Ｋ－（第１段起算点处的路线桩号）:Ｌ＝（在圆曲线上等于Ｓ；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第1段起算点的X坐标值）:F=（第1段起算点的Y 坐标值） ≠ >K<(第2段与第3段线路分界处的路线桩号):=>T=(第2段起算点处的切线方向角值):S=K-（第2段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第2段起算点的X坐标值）F=（第2段起算点的Y坐标值）≠ >K<(第3段与第4段线路分界处的路线桩号):=>T=(第3段起算点处的切线方向角值):S=K-（第3段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第3段起算点的X坐标值）F=（第3段起算点的Y坐标值）…… ≠ >K<(第n-1段与第n段线路分界处的路线桩号):=>T=(第n-1段起算点处的切线方向角值):S=K-（第n-1段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第n-1段起算点的X坐标值）F=（第n-1段起算点的Y坐标值） ≠>T=(第n段起算点处的切线方向角值):S=K-（第n段起算点处的路线桩号）:L=（在圆曲线上等于S；在缓和曲线上等于缓和段长度；在直线上为零。）:R=（曲线半径）:E=（第n段起算点的X坐标值）F=（第n段起算点的Y坐标值）⊿⊿⊿⊿⊿⊿⊿⊿(⊿共n-1个) L ≠ 0 = > C =90S^2/ЛRL :≠> C=0:⊿ L=S => V=2RsinC : Q=T-C : W=T-2C ≠> L=0 =>V=S : Q=T : W=T :≠> Pol(S-S^5/40R^2L^2+S^9/3456R^4L^4,S^3/6RL-S^7/336R^3L^3+S^11/42240R^5L^5):Fixm: V=I:Q=T-J:W=T-CΔΔ X=E+VcosQ: 手持GPS参数设置及全国各地坐标转换参数 ／298.2。西安80坐标系，其椭球的参数为：地球长半轴 a=6378140m；扁率F=1／298.257。国家2000坐标系，其椭球的参数为：地球长半轴a=6378137m；扁率F=1／298.298.257222101。 （三）手持GPS的参数设置 要想测量点位的北京54、西安80及国家2000公里网高精度坐标数据，必须学习坐标转换的基础知识，并分别科学设置手持GPS的各项参数。 首先，在手持式GPS接收机应用的区域内(该区域不宜过大)，从当地测绘部门收集1至两个已知点的北京54、西安80或国家2000坐标系统的坐标值；然后在对应的点位上读取WGS84坐标系的坐标值；之后采用《万能坐标转换》软件，可计算出DX、DY、DZ的值。 将计算出的DX、DY、DZ三个参数与DA、DF、中央经线、投影比例、东西偏差、南北偏差等六个常数值输入GPS接收机。将GPS接收机的网格转换为“UserGrid”格式，实际测量已知点的公里网纵、横坐标值，并与对应的公里网纵、横坐标已知值进行比较，二者相差较大时要重新计算或查找出现问题的原因。详细过程可查看《万能坐标转换》软件的【手持GPS参数设置】界面。 （四）自定义坐标系统（User）投影参数的确定 1、自己观测计算 新机拿到手之后，供应商都给提供一个投影参数，这对于要求不高的一般用户来说基本可以满足工作需要，而对于一些专业用户来说，就要自己来测算参数。一般型号的导航型手持GPS 自定义坐标系统（User）投影参数设置界面都提供了五个变量（△X、△Y、△Z、△A、△F）需要设置，而实际工作中，后两个参数（△A、△F）针对某一坐标系统来说为固定参数（北京54坐标系△A=-108、△F=0.0000005），无需改动，需要自己测算的参数主要为前三个（△X、△Y、△Z），一般称为三参数。 2、经验坐标 三参数对于非专业人员大多采用经验坐标，可别人的成果。 已知坐标点校正GPS的误差 （1）用GPS去测量已知坐标点得到坐标XGPS和YGPS； （2）计算两者的差值：△X＝XGPS－X已知△Y＝YGPS －Y已知 （3）计算FALSE′E′（东西偏差）和FALSE′N′（南北偏差）东西偏差＝500000－△X南北偏差＝0－△Y （4）更改GPS参数中的FALSE′E′（东西偏差）和FAL 坐标表示方法： 教你如何通过Excel VBA编写测量坐标计算程序 发布日期： 摘要：认识VBA、理解VBA，并利用Office Excel VBA编写测量坐标计算程序。 关键词：Excel VBA程序坐标编写 了解：VBA是什么？简单的说就是一种自动化语言，它可以使常用的程序自动化,可以创建自定义的解决方案。可以用Excel的宏语言来使Excel自动化运行等……Microsoft让它开发出来的应用程序共享一种通用的自动化语言——Visual Basic For Application(V BA)，可以认为VBA是非常流行的应用程序开发语言Visual Basic的子集，事实上VBA 是VB应用程序的版本，尽管存在有些不同VBA和VB在结构上仍然十分相似。如果你已经了解VB会发现学习VBA非常快。相应的学完VBA会给学习VB打下坚实的基础。 理由：选择Excel VBA编程的理由是因为它的计算功能非常强大，是现今任何编程计算器无法逾越的。它运用范围广，计算速度快，计算精度高，合理化显示等。或许很多测量人员对Excel VBA还有些陌生，主要是大家寄托于计算器、电脑、手机PDA等系列软件使用。Excel VBA对于大多数测量人员而没有系统学过计算机语言程序设计的人群来讲有一定含糊，不过只要有基本数学知识、测量常识和逻辑理解的人，都能通过Excel VBA编写设计出称心如意的测量程序。 目标：基于Excel VBA的测量坐标计算程序的设计目标是将繁琐计算过程转入到计算机中，利用程序语言的重复性原理，在计算机中可将坐标计算得出更精确的结果，使坐标计算更加可靠。最终目标是让用户可以通过Excel VBA自行完成坐标计算程序设计。 认识：学习VBA到底需要什么基础和了解些什么？ 学习VBA需要认识英文字母、一般的单词（如：函数所用的过程）、数学基础知识、测量常识、逻辑性思维即可。 在VBA中需要了解VBA的过程、变量、属性、方法、事件、语句等。 Excel VBA程序可以分为“录制宏、自定义函数”，由于录制宏编写计算类程序它限制了计算涵式过程，而无法达到自定义数据直接运算的目的，所以大家可以通过按钮式点击进行自定义函数过程（还可以通过窗体定义过程）。 基本常识： 算土方?算坐标?简单容易：https://www.360docs.net/doc/9117927166.html,/[/B]（道路填挖、河塘清淤回填算量、出图等）。 CASIO万能坐标程序 作者（马超） 程序特点： 真正的全线贯通坐标计算，在曲线元要素输入时仅需要输入第一段全部曲线元要素，后面曲线元要素除起点半径、终点半径、曲线长、转向需输入外其他要素均从前一曲线按辛普森8等分计算得出；辛普森公式任意等分，满足所有精度要求；全线曲线元数据一次性程序化输入，扩充变量数据库，无需修改程序内容；傻瓜操作，适用初级用户。 一、程序：MC W“1.JS 2.SZ”:W=1=>Z[2]=0:Goto 1 ΔW=2=> O“KOU LING”:O≠123456=>O=0: “OUT” ◢Goto 4ΔO=0: Z[1]=0:Goto 0←┘ Lbi 0←┘ ”N0.”：Z[1]+1 ◢←┘ Z[1]=0=>{ABCREFGU}:A“X0”:B“Y0”:C“F0”:R“R0”:E“RN”:F“D0”:G“LS”:U“G” : Z[Z[1]×8+3]=A:Z[Z[1]×8+4]=B:Z[Z[1]×8+5]=C:Z[Z[1]×8+6]= R-1:Z[Z[1]×8+7]= E-1:Z[Z[1]×8+8]=F: Z[Z[1]×8+9]=F+G: Z[Z[1]×8+10]=U: “NEXT”◢Isz Z[1]: Goto 0ΔZ[1]=1=>D=Z[9]:Z=0:Z[2]=0:GOTO 2ΔD=Z[(Z[1]-1)×8+9]:Z=0:Z[2]=Z[1]-1:GOTO 2←┘ Lbi A: Z[Z[1]×8+3]=X:Z[Z[1]×8+4]=Y:Z[Z[1]×8+5]=J: Z[Z[1]×8+8]=D: {REGU}:R“Ro”：E“RN”: G“LS”:U“G”: Z[Z[1]×8+6]=R-1 :Z[Z[1]×8+7]=E-1: Z[Z[1]×8+9]=D+G: Z[Z[1]×8+10]=U: “NEXT”◢Isz Z[1]: Goto 0←┘ Lbi 1←┘ {DZ}:D：Z：Z[2]=0:Goto 2←┘ Lbi 2←┘ Z[2]>Z[1]=>GoTo 4ΔD≤Z[Z[2]×8+9]=>A=Z[Z[2]×8+3]:B=Z[Z[2]×8+4]: C =Z[Z[2]×8+5]:R=Z[Z[2]×8+6]: E=Z[Z[2]×8+7]: F=Z[Z[2]×8+8]: G=Z[Z[2]×8+9]: U=Z[Z[2]×8+10]: Goto3ΔIsz Z[2]:Goto 2 Lbi 3←┘ W=2 =>N=8：≠>N=5ΔP＝U（E-R）÷Abs（G-F）：Q＝Abs（D-F）÷N：S=90Q÷π：J＝ C+(NPQ+2UR)NS：L=1←┘ X＝A+Q÷6×(Cos C+Cos J +4∑（Cos (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Cos (C+(LPQ+2UR)LS）,L,1,(N-1))）+ZCos（J+90)←┘ Y＝B+Q÷6×(Sin C+Sin J +4∑（Sin (C+((L+0.5)PQ+2UR)×(L+.5)S),L,0,(N-1)）+2∑（Sin (C+(LPQ+2UR)LS）,L,1,(N-1))）+ZSin（J+90):W=2=>GOTO AΔ Z=0=>“X”:X:Pause 0: “Y” :Y◢Goto 1Δ Z<0=>“XL”:X:Pause 0: “YL”:Y◢Goto 1Δfx4850 ① Z>0=>“XR”:X:Pause 0: “YR”:Y ◢Goto 1 ←┘ Z=0=> X “X” ◢Y “Y”◢Goto 1Δ Z<0=> X “XL” ◢Y “YL”◢Goto 1Δfx4800 ② Z>0=> X “XR” ◢Y “YR”◢Goto 1 ←┘ Lbi 4←┘5800简单全线坐标计算程序
线元法万能坐标计算程序
FX5800计算器公路全线坐标正、反算计算程序
CASIO万能坐标计算程序
卡西欧FX5800线路曲线坐标计算程序
线路中线和边线点位坐标的万能通用程序
卡西欧5800坐标计算程序(通俗易懂版)
坐标万能计算公式
手持GPS三参数计算及各地坐标转换经验参数
线路任意点坐标计算、及任意点对应桩号,左右偏距计算程序.tmp
手持GPS参数设置及全国各地坐标转换参数
教你如何通过EXCEL VBA编写测量坐标计算程序
卡西欧万能路线坐标计算程序(值得收藏)