最新对卡西欧4800计算器编程

对卡西欧4800计算
器编程
对卡西欧4800计算器在编程中的浅见
作者: luzhmu 发布日期: 2005-10-08 查看数: 793 出自:
在计算能力飞速发展的今天，各种计算工具应运而生代替了以往我们大量、繁琐的分析计算工作，简化了计算步骤节省了计算时间。

在种类繁多的计算工具中适合测量工作的计算器也大有所在。

比如有一款SHAPPC-E500计算器，它相当于是基于BASIC语言的微型电脑，功能也相当齐全能编写较复杂的程序、内存也相当大有32KB、64 KB、128 KB、256 KB、几种型号，同时具有输入输出通讯设备。

但价格不菲，现在市场价格大约为2000元到3000元，个人经济承受能力有限。

通过几年的实际测量工作我个人为认为真正适合测量人员较为专业的计算器应首推基于科学计算器具有简单程序语言的casio４x００计算器。

它有很多优点：一、使用简单，最大的一个优点就是只要具有初中水平会运用数学公式对函数有一般的了解的人，能基本理解casio４x００的内装函数即可进行简单的编程。

非常适合测量初学者和数学程序爱好者学习，也适合专业的测量人员的使用和能力提高。

不像其它的专业一定要具有较高的专业技术水平才能进行电脑编程。

二、成本低廉，价格在400元左右。

三、携带方便，体积很小可随时放在口袋里随拿随用。

casio系列较好的编程型号有casio４５００（以下简称４５００）、casio ４８００（以下简称４８００），（好象近来还推出了一款casio4850）前者较内存小，只有１１０３个字节，能应付一些较为简单的公式计算和科学计算，但由于内存有限，对一些较复杂或子程序过多的程序就力不从心了，不能出色的完成测量任务。

４８００就比４５００有较大的改进，４８００内存达到的４５００个字节，而且显示屏是４５００的几倍大，能更准确的显示数据，内装函数字符一目了然。

且具有简单的人机对话功能，出现了菜单子菜单。

４８００还在４５００增加了啊佛加德罗常数、万有引力长常数、详见《操作说明书》。

４８００有很多很实用的功能如他有一个公式解答功能，其原理是用牛顿法解方程。

举个例子：有一方程式：a=2b-c求当a=2、c=5时的b值。

将该方程式存入公式存贮器中，：先按按解答键“SOLVE” a输入2 、c输入5，再按解答键“SOLVE”计算器就会显示：
b=3.5。

此项功能被称为自动解答功能。

同时它也是非常实用的，在实际工作中通常要有一个经常使用的小公式，可借助它来完成。

４８００使用的程序语言可以算做简单的BASIC语言，有的命令如GOTO（转移到）、PAUSE（暂停）就与BASIC语言的一模一样。

现在４８００的程序语言来说说。

其主要命令有：1、=>…… 条件转移成立符号，其用法相当于BASIC 中的IF……THEN（假设语句相当于假如……然后，IF相当于条件……THEN相当于结果）语句2、≠>……条件转移不成立符号，其用法相当于BASIC中的IF……ELSE语句通常二者连用，相当于BASIC中的IF……THEN……ELSE语句（它的英语形式一般为if a＞b then c＞d else if b＞a the……）3、⊿…… 条件转移结束符号，与=>和≠>配合使用，放在条件语句最后面。

4、LbI……标记命令。

用于将一段语句作转换标记。

后
可接字母、数字、符号，但不能超过两个字节，如不能用≥10的数字作行标，否则会出现出错信息。

5、Goto…… （条件）转移命令。

前面可加条件语句，与BASIC中的GOTO作用相同。

通常与LbI一起用，如果所转移的行号无效，则会显示：GO ERROR（详见说明书）出错信息6、Dsz……减量循环命令。

可减少未知数的数量。

7、Isz……增量循环命令。

8、Pause……暂停命令。

后可接0~9之间的整数n，可使某一数据显示n/2秒钟，然后继续运行下面的程序。

它被认为是一个语句。

9、Fixm……变量锁定命令。

该命令能使其所有变量值（A~Z）均当成定数处理。

当程序运行时，将不需要输入变数（“{}”内的变数除外），而是将存贮器中原有的数值来完成计算。

10、{}…… 变量输入命令。

只程序在循环使用时经常发生改变的数字，如里程、和宽度。

它的输入方式可以使很多种如{AB}{A,B}{A B}都可以。

注意“{”和“}”必需成对输入。

否则会出现Syn ERROR（详见说明书）出错信息。

11、=、≠、>、<、≤、≥……条件运算关系算子，常与Goto命令构成条件转换语句。

12、Prog……在正常情况是下打开程序的快捷键。

在编程过程中是运行子程序命令，后接子程序名（一定要加引号，且要注意空格，否则会出现Syn ERROR（详见说明书）的出错信息。

13、↓……换行，只保留计算过程不显示计算结果。

当不想对其换行时也可用：代替。

14、◣……数据显示命令。

该命令输入后会自动换行。

保留计算过程并显示计算结果。

有一条总原则即：①学会运用程序的语言，尽可能使程序变得简明扼要；我们编写程序应该尽可能地使程序变得简明扼要，能省略的要一定省略。

烦琐的语句过多的字节只能使计算器的运算速度变慢没有任何好处，而且相当站用内存。

学会节省字节和使用符号是相当重要的。

尤其要灵活运用计算器语句因为它会使你更多的节省字节达到预期效果。

比如下程序就灵活运用了 Dsz （减量循环命令）。

比如使其能输入10个数值，并计算10个数值的平均值。

一般程序求10个数字的平均值需要有11个数字的提示符号。

但学会灵活运用了 Dsz （减量循环命令）那么只要有三个就可以了，这样就大大节省了字节的占用。

例程序如下：A=10
C=0 Lbi 1 {B} C=B+C Dsz A Goto 1
C÷10 但要注意的是：如果你是初学者或你对程序的编程不熟练，首先一定要先按照你的思路把程序步骤一步一步的列好在确定它能正确的计算后在想办法对其进行精简修改，否则只会使程序出现过多的错误；②尽可能使程序所包含的子程序减少；子程序过多就会造成程序结构松散，有的计算器主程序需要三个或四子程序，过多的子程序只会引起程序之间紊乱、混淆。

子程序过多对在使用时查找也比较麻烦。

而且子程序过多如果其中某个环节出现错误很难发现其错误所在，在编写程序时要尽量的少编写子程序，即使要编写子程序时也要注意尤其在容易出错的地方要多加注意。

有弊就有利如果你对子程序了解得多了那么可以几个主程序合用一个之程序也到到了要求的减少程序的字节使程序更简化。

例程序如下：CX CD Lbi 1 Lbi
1Prog”V” Prog”V”B=L-(K-S) ◣ B=L+(K-S) ◣Goto 1 Goto
1 V Y=√A2B2+B2X2÷A③尽量少用或不使用扩充变数存贮器，如
A[1]、A[2]等：使用扩充存贮器是一个利少弊多的做法。

每扩充一个存贮器就要减少10个字节的容量，而每个扩充存贮器至少要占四个字节，比一个A~Z变量净增三个。

有时你会觉得变量存贮器不够用。

其实不尽然，一般程序变量数很少会超过26个，只是你不懂得去使用。

一般来说，两个相对独立的程序步骤之间根本不需要考虑变量重复问题。

针对某一个程序，只要不是固定变量（{}内的变量），也就是那些通过计算出来的用于下一步计算的数值。

我们就可以通过重复赋值来得到某些计算量。

反正在下一轮循环中该量是变化的。

明白了各种命令的含义和注意事项就可以编程了。

举例有公式如下：CX 程序名称Lbl 0↓ 起始标记命令语句QMNFJ↓ 数据输入语句（指公式循环运算时的不变量）{KDE}↓ 数据输入语句（指公式循环运算时变量）S=K—Q:G=F+J↓ 公式运算命令X=M+ScosF↓ 公式运算命令Y=M+SsinF↓ 公式运算命令Prog ＂j＂↓ 运行子程序命令Goto 0↓ 循环运算语句J 子程序名称H=X+DcosG◣公式运算、数据显示语
I=Y+DsinG◣公式运算、数据显示语T=X—EcosG◣公式运算、数据显示语U=Y—EsinG◣公式运算、数据显示语最后计算器状态设定语句是大家最容易忽视的。

如果将单位进行预设那么计算器就会默认其使用单位在进行下一单位换算时要一定要进行单位转换，否者会使计算结果错误。

在显示屏幕的左下角可以清楚地发现小提示符号：如D代表度为现在的缺省单位、R代表弧度为现在的缺省单位、G 代表梯度为现在的缺省单位。

首页┊GPS测量┊工程测量┊遥感测量┊地籍测量┊控制测量┊数字测图┊地图制图┊数据处理┊测绘仪器┊软件应
用┊房产测量┊地理信息┊测量管理┊程序开发┊理论文章┊其他┊文章更新列表┊游客发布论文
所属分类：程序开发来源：作者：更新日期：2006-8-25 阅读次数：156 文章更新列表查看公告列表
公路编程实例
ZBZS★坐标正算主程序Deg:Fixm Lbl 0:｛O｝:O“1,47264-50167.162,50150-50673.503,50700-54700”Lbl 1:｛ZKB｝:Z:K:B:O=1=＞Goto 2△O=2=＞Goto 3△O=3=＞Goto 4△Lbl 2:Z＜47264=＞Goto 5△Z＞50167.16=＞Goto 5△Z＜48708.52=＞Prog“55”: Goto 6△Z＜49310.50=＞Prog“56”: Go to 6△Z ＜49874.21=＞Prog“57”: Goto 6△Prog“57-1”: Goto 6△Lbl 3:Z＜50150=＞Goto 5△Z＞50673.5=＞Goto 5△Prog“58”: Goto 6△Lbl 4:Z＜50700=＞Goto 5△Z＞54700=＞Goto 5△Z＜51265.12=＞Prog“59”: Goto 6△Z＜52266.78=＞Prog“60”: Goto 6△Z＜53341.82=＞Prog“61”: Goto 6△Z＜53868.03=＞
Prog“62”: Goto 6△Prog“63”: Goto 6△Lbl 5:“ZH ERROR”▲Goto 0Lbl 6:Prog“DQXY”:“X=”:X:Pause 0:“Y=”:Y▲Pol（X-Q,Y-U）:J＜0=＞J=360+J △:I=Intg（1000I+.5）÷1000“L=”:I:Pause 0:“T=”:J→DMS▲Goto 1 说明：本程序适用于casio4850计算器。

把线路从断链处分为三段，O的取值，待求点在第一段输1，第二段输2，第三段输3。

Z、待求点的桩号。

K、待求点在左幅输0右幅输1。

B、待求点距中距离。

Goto 5前的桩号为起止点或断链点桩号。

Prog“xx”:Goto 6前的桩号为曲线组合的HZ点桩号。

Q、U为测站坐标，运行程序前先存入。

ZHFS★桩号反算主程序Deg:Fixm Lbl 0:｛TL｝:T:L:T=Int T+Frac 100T÷36+（Int 100T÷100T-Int T）÷.6 X=Q+Rec（L,T）:Y=U+J:“X=”:X:Pause 0:“Y=”:Y▲ Lbl 1:Prog“55”:Prog“YS”:Z＞48980=＞Goto 2△Prog“FQ” Z＜47264=＞Goto B△Goto Z Lbl 2:Prog“56”:Prog“YS”:Z＞49437.85=＞Goto 3△Prog“FQ”: Goto Z Lbl 3:P rog“57”:Prog“YS”:Z＞49874.21=＞Goto 4△Prog“FQ”: Goto Z Lbl 4:Prog“57-1”:Prog“YS”:Z＞50354.96=＞Goto 5△Prog“FQ”Z＞50167.16 =＞Z=Z-17.16△Goto Z Lbl
5:Pro g“58”:Prog“YS”:Z＞50931.36=＞Goto 6△Prog“FQ”Z＞50673.5 =＞Z=Z+26.5△Goto Z Lbl 6:Prog“59”:Prog“YS”:Z＞51953.51=＞Goto 7△Prog“FQ”: Goto Z Lbl 7:Prog“60”:Prog“YS”:Z＞52546.8=＞Goto 8△Prog“FQ”: Goto Z Lbl 8:Prog“61”:Prog“YS”:Z＞53477.04=＞Goto 9△Prog“FQ”: Goto Z Lbl
9:Prog“62”:Prog“YS”:Z＞54118.18=＞Goto A△Prog“FQ”: Goto Z Lbl A:Prog“63”:Prog“YS”:Z＞54700=＞Goto B△Prog“FQ”: Goto Z Lbl B:“CX ERROR”▲Goto 0 Lbl Z:Z=Intg（1000Z+.5）÷1000:B=Intg（1000B+.5）÷1000:“LC=”:Z:Pause 0:“BJ=”:B▲Goto 0 说明：T、L测点相对测站的方位角和距离，T的取值例：12°34′56″输为12.3456。

Q、U为测站坐标，程序运行前先存入。

Goto 2前的桩号为第2个曲线组合的ZH点桩号，以后类推。

“BJ＝”显示正
值则测点在线路右幅，显示负值则测点在线路左幅。

DQXY★单曲线坐标正算子程序S=.5Abs（C-A）:G=（C-A）÷2S C-A＜-180=＞S=180-S:G=-G △C-A＞180=＞S=180-S:G=-G△T=（R+N2÷24R）tan S+N÷2-N^3÷240R2:L=πRS÷90+N:K=0=＞K=-1△Z≤D-T=＞Goto 1△Z≥D-T+L=＞Goto 2△Z＜D-T+N=＞Goto 3△Z＞D-T+L-N=＞Goto 4△Lbl 0:M=A+G（90+90N÷πR+180（Z-D+T-N）÷πR）-180 X=E+Rec（（R+ N2÷24R）÷cos S,C+G（90-S））:Y=F+J X=X+Rec（R-GKB,M）:Y=Y+J: Goto 5 Lbl 1:X=E-Rec（D-Z,A）:Y=F-J:X=X+Rec（B,A+90K）:Y=Y+J: Goto 5 Lbl 2:X=E+Rec（Z-D+2T-L,C）:Y=F+J:X=X+Rec（B,C+90K）:Y=Y+J: Goto 5 Lbl 3:X=E-Rec（T,A）:Y=F-J:H=Z-D+T:Prog“XY”X=X+Rec（P,A+Gtan-1
（W÷V））:Y=Y+J X=X+Rec（B,A+90GH2÷πRN+90K）:Y=Y+J: Goto 5 Lbl 4:X=E+Rec（T,C）:Y=F+J:H=D-T+L-Z:Prog“XY”X=X+Rec（P,C+180-Gtan-1（W÷V））:Y=Y+J X=X+Rec（B,C-90GH2÷πRN+90K）:Y=Y+J: Goto 5 Lbl 5 XY★缓和曲线参数计算子程序V=H-
H^5÷40R2N2+H^9÷3456R^4N^4 W=H^3÷6RN-H^7÷336R^3N^3+H^11÷42240R^5N^5:P=√（V2+W2）YS★曲线要素计算子程序S=.5Abs（C-A）:G=（C-A）÷2S C-A＜-180=＞S=180-S:G=-G△C-A＞180=＞S=180-S:G=-G△T=（R+N2÷24R）tan S+N÷2-N^3÷240R2:L=πRS÷90+N Pol（X-E,Y-F）:J=J-C:Z=Icos J-T:Z=D-T+L+Z FQ★单曲线反求桩号子程序Z≥D-T+L=＞Z:B=Isin J:Goto A△D=D-T+L÷2:V=E+Rec（（R+N2÷24R）÷cos S,C+G （90-S））:W=F+J O=C+G（90-S）+180:O＞360=＞O=O-360△Pol（X-V,Y-W）:J＜0=＞J=360+J△Z=D+G（J-O）πR÷180:B=G（R-I）Z＜D-（L-
2N）÷2=＞Goto 1△Z＞D+（L-2N）÷2=＞Goto 2△Goto A Lbl 1:D=D-L÷2:E=E-Rec（T,A）:F=F-J:P=A-90:G=-G:Prog“HHXY”:Z=D+T: Goto A Lbl
2:D=D+L÷2:E=E+Rec（T,C）:F=F+J:P=C+90:Prog“HHXY”:Z=D-T: B=-B:Goto A Lbl A HHXY★反求缓和曲线长度子程序Pol（X-E,Y-F）:J=J-
P:X=Icos J:Y=Isin
J B=G÷46080R^6N^6:H=X÷3840R^5N^5:I=13G÷5760R^4N^4:J=GY÷384R^4N^4 K=X÷48R^3N^3:M=G÷15R2N2:O=GY÷8R2N2:P=X÷2RN:S=GY:Z=Y Lb 2:V=BZ^13+HZ^10-IZ^9+JZ^8-KZ^6+MZ^5-OZ^4+PZ2-GZ+S W=13BZ^12+10HZ^9-9IZ^8+8JZ^7-6KZ^5+5MZ^4-4OZ^3+2PZ-G T=Z-V÷W:Abs（T-Z）＜1E-5=＞Goto 3△Z=T:Goto 2 Lbl 3:V=G（T^3÷6RN-T^7÷336R^3N^3+T^11÷42240R^5N^5）W=T-T^5÷40 R2N2+T^9÷3456R^4N^4 V=X=＞W=Y=＞
B=0:Goto 5△△Pol（X-V,Y-W）:B=sin（J+90GT2÷πRN-90）:S=-1:B＞0=＞S=1△B=SI Lbl 5 55★ 55～63为平曲线参数置数子程序
D=48212.03:E=42650.6693:F=16021.6936:A=247°56′58″:C=151°56′02″:R=826.192:N=60 56★
D=49148.91:E=41425.5219:F=16674.9305:A=151°56′02″:C=185°56′02″:R=354.724:N=120 57★
D=49661.59:E=40908.3058:F=16621.1706:A=185°56′02″:C=151°54′14″:R=600:N=80 57-1★
D=50019.67:E=40582.613:F=16795.0457:A=151°54′14″:C=172°51′03″:R=570.059:N=8058★
D=50498.68:E=40087.8038:F=16857.11:A=172°51′03″:C=184°48′17″:R=1250:N=6059★
D=51112.95:E=39501.143:F=16807.7977:A=184°48′17″:C=206°30′23″:R=600:N=8060★
D=52113.01:E=38603.5949:F=16360.1705:A=206°30′23″:C=237°03′19″:R=400:N=10061★
D=52962.43:E=38138.5471:F=15642.5421:A=237°03′19″:C=193°17′44″:R=910:N=10062★
D=53677.91:E=37406.98:F=15469.6681:A=193°17′44″:C=228°55′55″:R=500:N=8063★
D=54344.49:E=36962.0015:F=14959.0062:A=228°55′55″:C=181°27′44″:R=400:N=80说明：线路起点K47＋264，止点K54＋700。

共十个交点，两处断链，第一处K50+ 167.16＝K50+150，长17.16M。

第二处K50＋673.5＝K50＋700，短26.5M。

置数子程序中，D交点桩号。

E、F交点坐标。

A、C交点前后边正方位角。

R曲线半径。

N缓和曲线长。

程序“HHXY”引用借鉴了徐光辉2000年7月5日编写的“《CASIO fx-4800P线路计算 6.0版》程序设计”中“反求缓曲长度子程序HQ”的计算式。

GCjs★路面高程计算主程序Lbl 0:｛O｝:O“1,47264-50167.162,50150-50673.503,50700-54700”Lbl 1:
｛Z:K:B｝:Z:K:B:Fixm 意义同坐标正算程序V=2:K=0=＞K=-1△B＜.75=＞B=.75△B＞11.25 =＞ Goto 5△ .75为1/2中央分隔带宽O=1=＞Goto 2△O=2=＞Goto 3△O=3=＞Goto 4△ 11.25为路面半幅宽度Lbl 2:Z＜47264=＞Goto 5△Z＞50167.16=＞Goto 5△Prog“A”: Prog“SQ”:
Prog“O”: Goto 6Lbl 3:Z＜50150=＞Goto 5△Z＞50673.5=＞Goto 5△Prog“B”: Prog“SQ”: Prog“P”: Goto 6Lbl 4:Z＜50700=＞Goto 5△Z＞54700=＞Goto 5△Prog“C”: Prog“SQ”:Prog“Q”: Goto 6 Lbl 5:“ZH ERROR”▲Goto 0Lbl 6:Prog“CG”:“H=”:H▲Goto 1 SQ★竖曲线计算子程序Z≤D-T=＞Goto 1△Z≥D+T=＞I=J: Goto 1△Lbl 0:H=G-（D-Z）I÷100+（Z-D+T）2÷2RN: Goto 2 Lbl 1:H=G-（D-Z）I÷100: Goto2 Lbl 2 CG★超高计算子程序（弯道未设加宽）P=2VS÷（V+W）Z＜M=＞Goto 1△Z＞M+L=＞Goto 1△Z≤M+P=＞F=Z-M:Goto 2△Z≥M+L-P=＞F=M+L-Z:Goto 2△Z＜M+S=＞F=Z-M:Goto 3△Z＞M+L-S=＞F=M+L-Z:Goto 3△Lbl 0:X=H+（B-.75）WK÷100:Y=H-（B-.75）WK÷100 E=-1=＞H=X:≠＞H=Y△Goto 5 Lbl 1:H=H-（B-75）V÷100: Goto 5 Lbl 2:X=H-（B-.75）V÷100:Y=H+（（V+W）F÷S-V）（B-.75）÷100 E=-1 =＞K=-1 =＞H=X: ≠＞H=Y△≠＞K=1=＞H=X: ≠＞H=Y△△Goto 5 Lbl 3:X=H+（（V+W）F÷S-V）（B-.75）K÷100 Y=H-（（V+W）F÷S-V）（B-.75）K÷100 E=-1=＞H=X: ≠＞H=Y△Goto 5 Lbl 5 说明：竖曲线参数置数子程序中，D竖曲线顶点桩号。

G顶点高程。

I竖曲线前段纵坡。

J竖曲线后段纵坡。

R竖曲线半径。

T切线长。

N凹曲线取1凸曲线取-1。

平曲线参数置数子程序中，M起超点桩号。

S超高缓和段长。

L超高曲线总长。

W最大单向横坡坡比。

E曲线左偏取-1右偏取1。

V直线段双向横坡坡比。

纵坡按设计给出的正负取值，横坡均按正数取值。

坡比取值示例：2.5%取2.5。

A★Z＜47635=＞Goto 1△Z＜48680.5=＞Goto 2△Z＜49200=＞Goto 3△Z＜49369.25=＞Goto 4△Lbl 0:D=50040:G=1025.21:I=-6:J=-3:R=4000:T=60.05:N=1: Goto 5 Lbl 1:D=47560:G=1135.61:I=-6:J=-
3:R=5000:T=75:N=1: Goto 5 Lbl 2:D=48640:G=1103.21:I=-3:J=-6:R=2700:T=40.5:N=-1: Goto 5 Lbl 3:D=49140:G=1073.21:I=-6:J=-3:R=4000:T=60:N=1: Goto 5 Lbl 4:D=49340:G=1067.21:I=-3:J=-6:R=1950:T=29.25:N=-1: Goto 5 Lbl 5 B★D=50220:G=1019.3:I=-3:J=-6:R=3000:T=45.04:N=-1
C★Z＜50957.5 =＞Goto 1△Z＜51125 =＞Goto 2△Z＜51805 =＞Goto 3△Z＜51953 =＞Goto 4△Z＜52546 =＞Goto 5△Z＜52727.44 =＞Goto 6△Z＜52950.32 =＞Goto 7△Z＜53800 =＞Goto 8△Lbl 0:D=54400:G=835.84:I=-2.5:J=-1.2:R=22000:T=143:N=1: Goto 9 Lbl 1:D=50920:G=978.89:I=-6:J=-
3:R=2500:T=37.5:N=1: Goto 9 Lbl 2:D=51080:G=974.09:I=-3:J=-6:R=3000:T=45:N=-1: Goto 9 Lbl 3:D=51760:G=933.29:I=-6:J=-3:R=3000:T=45:N=1: Goto 9 Lbl 4:D=51920:G=928.49:I=-3:J=-6:R=2200:T=33:N=-1: Goto 9 Lbl 5:D=52500:G=893.69:I=-6:J=-2:R=2300:T=46:N=1: Goto 9 Lbl
6:D=52680:G=890.09:I=-2:J=-5:R=2500:T=47.44:N=-1: Goto 9 Lbl 7:D=52900:G=877.34:I=-5:J=-3:R=3600:T=50.32:N=1: Goto 9 Lbl 8:D=53700:G=853.34:I=-3:J=-2.5:R=40000:T=100:N=1: Goto 9 Lbl 9 A、B、C 为竖曲线参数置数子程序O★Z＜48708.52 =＞Goto 1△Z＜49310.5 =＞Goto 2△Z＜49874.21 =＞Goto 3△Lbl 0:M=49874.21:S=80:L=288.41:W=3:E=1: Goto 4 Lbl 1:M=47264:S=60:L=1444.52:W=2:E=-1: Goto 4 Lbl 2:M=48980:S=120:L=330.5:W=5:E=1: Goto 4 Lbl 3:M=49437.85:S=80:L=436.36:W=3:E=-1: Goto 4 Lbl 4 P★M=50337.8:S=60:L=320.8:W=2:E=1 Q★Z＜51265.12 =＞Goto 1△Z＜52266.78 =＞Goto 2△Z＜53341.82 =＞Goto 3△Z＜53868.03 =＞Goto 4△Lbl 0:M=54118.18:S=100:L=431.41:W=5:E=-1: Goto 5 Lbl
1:M=50957.86:S=80:L=307.26:W=3:E=1: Goto 5 Lbl 2:M=51953.51:S=100:L=313.27:W=5:E=1: Goto 5 Lbl 3:M=52546.8:S=100:L=795.02:W=2:E=-1: Goto
5 Lbl 4:M=53477.04:S=80:L=390.99:W=4:E=1: Goto 5 Lbl 5 O、P、Q 为平曲线参数置数子程序--------------------------------------------------------------------------------
【目前共有0 条对该新闻的评论】【查看参与评论】
收集来自网上，版权归原来作者所有 ,如果侵犯您的版权，请来信告诉我们，我们会马上做出整理，谢谢Copyright © 2004 测量论
文中心 All Rights Reserved
(适用于无缓和曲线的平曲线)
本程序系本人精心编制而成,公路全线只需要一个数据库,欢迎各位测友下载试用！
知识共享是本人的宗旨。

联系方式： QQ 303223247漂人
更多程序请参/u/1221782002
附件(点击下载)：
主程序名:ZBJS(坐标计算)
J=E-D:J=Arg(cosJ+ⅰsinJ):H=AbsJ÷J:T=Rtan （JH÷2）:L= JHRπ÷180:
K＜C-T⇒I=C-K:F=0:Goto 1⊿
K＜C-T+L⇒W=K-C+T:F=180W÷Rπ:I=T-RsinF+(R-RcosF)ⅰ⊿
Lb1 1:S“PJU”:Z“PJIAO”:
“X=”:A-cos(D-HArgI)AbsI+Scos(D+HF+Z)◢
“Y=”:B-sin(D-HArgI)AbsI+Ssin(D+HF+Z)
数据库:SJK (全线交点要素一个数据库)
K＜YZ点桩号⇒A=本交点X坐标:B=本交点Y坐标:C=本交点桩号:
D=前直线方位角:E=后直线方位角:R=曲线半径:Goto 0⊿
……………
Lbl 0:P rog“ZBJS”
程序说明:
K? 输入桩号
PJU? 输入边桩至中桩的斜距
PJIAO? 输入边桩方向与路线前进方向的右夹角
(n是圆周率,如程序用于4800计算器，显示部分需要更改才能正常：“X=”：更改为
X“X=”=
“Y=”：更改为Y“Y=”=
谢谢发贴！平均用户评分：+5【测量空间】
大中小主程序名:ZBJS(坐标计算)
J=E-D:J=Arg(cosJ+ⅰsinJ):H=AbsJ÷J:
P=R+G2÷24R:Q=G÷2-G3÷240R2:M=JH÷2:
N=90÷Rπ:V=RG:T=PtanM+Q:L=M÷N+G:
U=C-T:S“PJU”:Z“PJIAO”:
K＜U=>I=C-K:F=0:Goto 1⊿
K＜U+G=>W=K-U:Lb1 2:F=NW2÷G:I=T-W+W5÷40V2+ⅰW3÷6V:Goto 1⊿
K＜U+L-G=>W=K-U-G:F=N(2W+G):I=T-RsinF-Q+(P-RcosF)ⅰ:Goto 1⊿
K＜U+L=>W=U+L-K:D=E+180:H=-H:S=-S:Goto 2⊿
Lb1 1:“X=”:A-cos(D-HArgI)AbsI+Scos(D+HF+Z)◢
“Y=”:B-sin(D-HArgI)AbsI+Ssin(D+HF+Z)
数据库:SJK(全线交点要素)
K＜HZ点桩号=>A=本交点X坐标:B=本交点Y坐标:C=本交点桩号:
D=前直线方位角:E=后直线方位角:R=曲线半径:G=缓和曲线长:Goto 0⊿……………
Lbl 0:Prog“ZBJS”
程序说明:
K? 输入桩号
PJU? 输入边桩至中桩的斜距
PJIAO? 输入边桩方向与路线前进方向的右夹角
(如程序用于4800，显示部分需要更改才能正常：
“X=”：更改为X“X”=；
“Y=”：更改为Y“Y”=
程序名：ZDZH(主点桩号计算)
Q=C“LS”/2-C3/240R2:P=C2/24R-C4/2384R3:T=(R+P)tan
(AbsJ“PJ”/2)+Q:L=RπAbsJ/180+C:B=Z“JDZH”-T:
“ZH=”:B◢
“HY=”:B+C◢
“QZ=”:B+L/2◢
“YH=”:B+L-C◢
“HZ=”:B+L
程序说明:
LS? 输入缓和曲线长
R? 输入半径
PJ? 偏角
JDZH? 交点桩号
程序名：WJPQX
A “PJ”:E:L“LS”:
N=1÷cos(0.5A):M=E÷(N-1):W=NL2÷(24N-24):Q=M2-4W:
Q≥0⇒“R=”:0.5(M+√Q)⊿
程序说明:本程序适用于外距控制单交点平曲线PJ? 输入路线转角
E? 输入外距
LS? 输入缓和曲线长
2006-03-2817:52:15
大中小程序名:LZXJD(两直线交点)
A"X1":B"Y1":C"X2":D"Y2":E"X3":F"Y3":G"X4":H"Y4":
N=C-A:O=D-B:P=G-E:Q=H-F:K=QN-OP:U=NB-OA:V=PF-QE:
"X=":(PU-NV)/K◢"Y=":(QU-OV)/K
程序说明:
X1? 输入第一点X坐标
Y1? 输入第一点Y坐标
X2? 输入第二点X坐标
Y2? 输入第二点Y坐标
X3? 输入第三点X坐标
Y3? 输入第三点Y坐标
X4? 输入第四点X坐标
Y4? 输入第四点Y坐
casio4850全线坐标计算程序
一、LXZB（坐标计算主程序）
LbI Z：｛P｝
A=×××：V=××°××′××″：C=×××：D=×××：
P≤直线的终点桩号＝＞Prog“ZX.SUB”：Goto Z （以上各字母含义详直线坐标计算子程序说明）
A=×××：V=××°××′××″：C=×××：D=×××：
R=×××：I=（1或-1）
P≤圆曲线终点桩号＝＞Prog“YQX.SUB”：Goto Z
（以上各字母含义详圆曲线坐标计算子程序说明）
A=×××：V=××°××′××″：C=×××：D=×××：
B=×××：W=××°××′××″：E=×××：F=×××：
R=×××：L=×××：I=（1或-1）
P≤HY点桩号＝＞Prog“HQ1.SUB”：Goto Z
P≤YH点桩号＝＞Prog“HYQ.SUB”：Goto Z
P≤HZ点桩号＝＞Prog“HQ2.SUB”：Goto Z
（以上各字母含义详（HQ1.SUB、HYQ.SUB、HQ2.SUB）坐标计算子程序说明）该主程序是简单的示例说明，具体的应根据图纸的线型组合，随机应变的做顺序调整，并将其输入Casio4850计算器。

二、ZX.SUB（直线坐标计算子程序）
G=P-A：T=V
“X=”：X=C+GcosT◢
“Y=”：Y=D+GsinT◢
｛Q｝：Q=1 =>Prog“BZ.SUB”Δ
说明：
G：所求桩号至起算点的距离
P：所求点的桩号
A：起算点的桩号
T：所求桩号的方位角
V：起始方位角
C、D：起算点的X、Y坐标
X、Y：显示所求桩号的坐标
Q：参数，为1时计算边标坐标，为0时不计算三、YQX.SUB（圆曲线坐标计算子程序）
G=P-A
H=Rsin((G÷R)×(180÷π))
J=(R(1-cos((G÷R)×(180÷π))))I
T=V+((G÷R)×(180÷π))I
“X=”:X=C+HcosV-JsinV ◢
“Y=”:Y=D+HsinV+JcosV◢
｛Q｝:Q=1 =>Prog“BZ.SUB”Δ
说明：
G：所求桩号至起算点的距离
P：所求点的桩号
A：起算点的桩号(指ZY点)
R：圆曲线半径
T：所求桩号的方位角
V：起始方位角(指ZY点的方位角)
C、D：起算点（ZY点）的X、Y坐标
X、Y：显示所求桩号的坐标
I：参数，（左偏为－1，右偏为＋1）
Q：参数，为1时计算边标坐标，为0时不计算
四、HQ1.SUB（第一缓和曲线坐标计算子程序）
G=P-A
H=G-G^5÷(40R2L2)+G^9÷（3456R^4L^4）-G^13÷（599040 R^6L^6）J=(G^3÷(6RL)-G^7÷(336R^3L^3)+G^11÷(42240R^5L^5))I
T=V+90G2÷(πRL)I
“X=”:X=C+HcosV-JsinV ◢
“Y=”:Y=D+HsinV+JcosV◢
｛Q｝:Q=1 =>Prog“BZ.SUB”Δ
说明：
G：所求桩号至起算点的距离
P：所求点的桩号
A：起算点的桩号(指ZH点)
R：曲线半径
L：缓和曲线长Ls
T：所求桩号的方位角
V：起始方位角(指ZH点的方位角)
C、D：起算点（ZH点）的X、Y坐标
X、Y：显示所求桩号的坐标
I：参数，（左偏为－1，右偏为＋1）
Q：参数，为1时计算边标坐标，为0时不计算
五、HYQ.SUB（带缓和曲线的圆曲线坐标计算子程序）
G=P-A
H=Rsin((L÷2÷R)×180÷π+(G-L)÷R×180÷π)+L÷2-L^3÷(240R2)
J=(R(1-cos((L÷2÷R)×180÷π+(G-L)÷R×180÷π))+L2÷(24R)-L^4÷(2688R3))I T=V+90L÷(πR)I+(180(G-L)÷(πR))I
“X=”:X=C+HcosV-JsinV ◢
“Y=”:Y=D+HsinV+JcosV◢
｛Q｝:Q=1 =>Prog“BZ.SUB”Δ
说明：
G：所求桩号至起算点的距离
P：所求点的桩号
A：起算点的桩号(指ZH点)
R：曲线半径
L：缓和曲线长Ls
T：所求桩号的方位角
V：起始方位角(指路线前进方向ZH点的方位角)
C、D：起算点（ZH点）的X、Y坐标
X、Y：显示所求桩号的坐标
I：参数，（左偏为－1，右偏为＋1）
Q：参数，为1时计算边标坐标，为0时不计算
六、HQ2.SUB（第二缓和曲线坐标计算子程序）
G=B-P
H=G-G^5÷(40R2L2)+G^9÷（3456R^4L^4）-G^13÷（599040 R^6L^6）
J=-(G^3÷(6RL)-G^7÷(336R^3L^3)+G^11÷(42240R^5L^5))I T=W-90G2÷(πRL)I
“X=”:X=E-HcosW+JsinW ◢
“Y=”:Y=F-HsinW-JcosW ◢
｛Q｝:Q=1 =>Prog“BZ.SUB”Δ
说明：
G：所求桩号至起算点的距离
P：所求点的桩号
B：起算点的桩号(指HZ点桩号)
R：曲线半径
L：缓和曲线长Ls
T：所求桩号的方位角
W：起始方位角(指路线前进方向HZ点的方位角)
E、F：起算点（HZ点）的X、Y坐标
X、Y：显示所求桩号的坐标
I：参数，（左偏为－1，右偏为＋1）
Q：参数，为1时计算边标坐标，为0时不计算
七、BZ.SUB（边桩坐标计算子程序）
｛O，S｝：
O“ZBK”：S“YBK”：“XZ=”：M=X-O cos(T+90) ◢“YZ=”：N=Y-O sin(T+90)◢
“XY=”：K=X+S cos(T+90) ◢
“YY=”：U=Y+S sin(T+90) ◢
说明：ZBK：中桩至左边桩的距离YBK：中桩至右边桩的距离XZ、YZ：显示所计算左边桩的X、Y坐标XY、YY：显示所计算右边桩的X、Y坐标
线路中边桩任意里程坐标正反算改编程序CASIO4850
Posted by 忘不了の糊涂 2006-10-5 23:45:00
含数据库、测站放样计算（第五次修改-正式版）
1.加编数据库及计算总调度程序,计算中不必逐项输入"线元要素",提高运算速度，避免现场忙中出错
2.可加入多条线路的数据库,内业输入数据,外业一目了然
3.计算直观,人性化
4.正算直接输入里程和边距,反算输入近似里程便可
5.增加了“计算点与测站点”的距离和方位角计算语句，方便直接放样
6.愿收获与大家共享,同时也希望大家提出心得和宝贵建议
7.核心计算程序摘自“yshf”
一.改动后的程序清单，增设数据库程序,可输入n条线路的数据库
1.(QXZDJS 计算总调度程序)
M=1=>Prog "SJK1":Prog "SUBSJK"△←┘
...........
M=n=>Prog "SJKn":Prog "SUBSJK"△←┘
2.(SJK1 数据库程序)
"1.SZ => XY"："2.XY => SZ"：｛N,S｝：N：S"DKI"←┘
S≤本线元终点里程＝＞U=本线元起点X坐标：V=本线元起点Y坐标：G=本线元起点正切线方位角：P=本线元起点曲率半径：R=本线元终点曲率半径：O=本线元起点里程：H=本线元终点里程：Q=本线元左直右偏向(0或1、-1):Goto 0Δ←┘
S≤下一线元终点里程=＞O=。

：Goto0Δ←┘LB1 0
3.(SUBSJK 运算主程序)
{UVGPROHQ}：U"XA"：V"YA"：G"FA"：P"RA"：R"RB"：O"DKA"：H"DKB "：Q←┘
Deg:S>H=>Goto AΔD=(P-R)÷2PR(H-O)：N=1=>Goto 1：≠>Goto 2Δ←┘
Lbl 1：{Z}：Z：W=Abs(S-O)：Prog "SUB1"："XS="：X◢
"YS="：Y◢
"QX-FWJ="：F=F-90←DMS◢
I=0：J=0：Pol（X-C"XC"，Y-E"YC"）：J<0＝＞J=J+360：Δ“FWJ=”:J←DMS ◢
“I=”：I◢
Goto A←┘
Lbl 2：{XY}：I=X：J=Y：Prog "SUB2"："S=":S=O+W◢
"Z=":Z◢
Lbl A
4. (SUB1 正算子程序)
A=0.1739274226：B=0.3260725774：K=0.0694318442：L=0.3300094782：X=U+W(Acos(G+180QKW(1÷P+KWD)÷π)+Bcos(G+180QLW(1÷P+LWD)÷π)+ Bcos(G+180(1-L)QW(1÷P+(1-L)WD)÷π)+A cos(G+180(1-K)QW(1÷P+(1-K)W D)÷π))：Y=V+W(Asin(G+180QKW(1÷P+KWD)÷π)+Bsin(G+180QLW(1÷P+L WD)÷π)+Bsin(G+180(1-L)QW(1÷P+(1-L)WD)÷π)+Asin(G+180(1-K)QW(1÷P+ (1-K)WD)÷π))：F=G+180QW(1÷P+WD)÷π+90：X=X+ZcosF：Y=Y+ZsinF
5.(SUB2 反算子程序)
T=G-90：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+18 0QW(1÷P+WD)÷π：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1÷10^6=>Goto1：≠>W=W +Z：Goto 0Δ←┘
Lbl 1：Z=0：Prog "SUB1"：Z=(J-Y)÷sinF
二. 使用说明：
1.把所有相关的”线元要素“依次输入”SJK“，如果有多条线路，可把“SJK”分成SJK1...SJKn等n条线路数据库，然后对应“计算总调度程序”中的M1...Mn备用
2.运算时直接调用“SJK”运行，程序提示输入里程“S”？时，正算直接输入待求点里程，反算输入所求点“近似“里程
3.如果在正算中途想进行反算，重新给”N“？赋值"2"便可
4.程序中“XC、YC”为测站坐标，“I＝”、“FWJ＝”为放样距离和方位角
5.程序中“QX-FWJ=”为计算点在中线上的正切线方位角，此方位角在新线路输入数据库时非常重要，我们在每输入完一段曲线元要素后，立即退出返回到“计算总调度程序(SUBSJK)”中，然后输入该段曲线元的终点里程，其计算出来的X，Y中桩坐标应该与设计图纸(当然是经过复查无误的图纸)上的该点中桩坐标一致，否则一定是哪个地方数据输入或者程序有误，这时候我们一定要把该计算点在中线上的正切线方位角“QX-FWJ=”记录下来为下一段曲线元要素数据输入服务，因为它就是下一段曲线元的起点正切线方位角，嘿嘿！
6.其它输入与"yshf"原程序一致
-----------------------------------------------------------------------
附件:(原yshf程序)
曲线任意里程中边桩坐标正反算(CASIO fx-4800P计算器)程序
一、程序功能
本程序由一个主程序(TYQXJS)和两个子程——正算子程序(SUB1)、反算子程序(
SUB2)序构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线
元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲
率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

另
外也可以将本程序中核心算法部分的两个子程序移植到其它相关的程序中，用于对曲
线任意里程中边桩坐标进行正反算。

本程序也可以在CASIO fx-4500P计算器及 CASIO fx-4850P计算器上运行。

特别申明：(1). 适用于弧长小于2倍半径的各种线元坐标正反算，精度优于1mm；
(2). 在引用该核心计算部分时，请注明来源。

二、源程序
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=>Goto 1：≠>Goto 2Δ←┘
Lbl 1：{SZ}：SZ：W=Abs(S-O)：Prog "SUB1"：X"XS"=X◢
Y"YS"=Y◢
F"FS"=F-90◢
Goto 1←┘
Lbl 2：{XY}：XY：I=X：J=Y：Prog "SUB2"：S"S"=O+W◢
Z"Z"=Z◢
Goto 2
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+LW D))+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：Lbl 0：Prog "SUB1"：L=T+Q EW(C+
WD)：Z=(J-Y)cosL-(I-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘
Lbl 1：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) 当线元为完整缓和曲线时，起点与直线相接时，曲率半径为无穷大，以1 0的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=×××反算时，计算得出的所求点的边距
四、算例
某匝道的由五段线元（直线+完整缓和曲线+圆曲线+非完整缓和曲线+直线）组
成，各段线元的要素（起点里程S0、起点坐标X0 Y0、起点切线方位角F0、线元长度
LS、起点曲率半径R0、止点曲率半径RN、线元左右偏标志Q）如下：
S0 X0 Y0 F0 LS R0 RN Q
500.000 19942.837 28343.561 125 16 31.00 269.256 1E45 1E4
5 0
769.256 19787.340 28563.378 125 16 31.00 37.492 1E45 221.7 5 -1
806.748 19766.566 28594.574 120 25 54.07 112.779 221.75 221.7 5 -1
919.527 19736.072 28701.893 91 17 30.63 80.285 221.75 9579.22 8 -1
999.812 19744.038 28781.659 80 40 50.00 100.000 1E45 1E4
5 0
1、正算
（注意：略去计算方式及线元要素输入，请自行根据所求点所在的线元输入线元要素）
S=700 Z=-5 计算得 XS=19831.41785 YS=28509.72590
S=700 Z=0 计算得 XS=19827.33592 YS=28506.83837
S=700 Z= 5 计算得 XS=19823.25398 YS=28503.95084
S=780 Z=-5 计算得 XS=19785.25749 YS=28575.02270
S=780 Z=0 计算得 XS=19781.15561 YS=28572.16358
S=780 Z= 5 计算得 XS=19777.05373 YS=28569.30446
S=870 Z=-5 计算得 XS=19747.53609 YS=28654.13091
S=870 Z=0 计算得 XS=19742.68648 YS=28652.91379
S=870 Z= 5 计算得 XS=19737.83688 YS=28651.69668
S=940 Z=-5.123 计算得 XS=19741. 59118 YS=28722.05802
S=940 Z=0 计算得 XS=19736.47687 YS=28722.35642
S=940 Z= 3.009 计算得 XS=19733.47298 YS=28722.53168
2、反算
X=19831.418 Y=28509.726 计算得 S=699.9999974 Z= -5 .00018 164
X=19827.336 Y=28506.838 计算得 S=699.9996493 Z= 0.000145 136
X=19823.25398 Y=28503.95084 计算得 S=699.9999985 Z= 5.0000 03137
X=19785.25749 Y=28575.02270 计算得 S=780.0000035 Z= -5 .000 001663
X=19781.15561 Y=28572.16358 计算得 S=780.0000025 Z=- 0.000 002979
X=19777.05373 Y=28569.30446 计算得 S=780.0000016 Z= 4.9999 9578
X=19747.536 Y=28654.131 计算得 S=870.0001137 Z= -4.99941 049
X=19742.686 Y=28652.914 计算得 S=870.0003175 Z=- 0.00041 814
X=19737.837 Y=28651.697 计算得 S=870.0002748 Z= 4.999808 656
X=19741.5912 Y=28722.0580 计算得 S=939.9999786 Z= -5.1230 24937 X=19736.4769 Y=28722.3564 计算得 S=939.9999862 Z=- 0.000027710 X=19733.4730 Y=28722.5317 计算得 S=940.00002 38 Z= 3.00898694
FX4800P常用测量程序集
1、由一个已知点坐标和一个已知方位角计算支导线坐标（程序名：ZBJS）
程序清单：
X：Y：T：Lb1 0：Fixm：{AS}：T=T+A-180：T〈0 =〉T=T+360⊿ T〉360 =〉T= T-360⊿T◢
X=X+Scos T◢
Y=Y+Ssin T◢
Goto 0
程序运行说明：
运行“ZBJS”，每运行一步按EXE键一次，显示X？时输入已知点X坐标，显示Y？时输入已知点Y坐标，显示T？时输入已知点由后视到前视的已知坐标方位角，运行一次本程序只需输入一次已知数据。

显示A？时输入第一个左转角值，显示S？时输入测站到前视的水平距离。

接下来计算器显示的T是测站到前视的坐标方位角，按 + 可将显示的十进制角度换算为度分秒。

如此继续可计算出所有支导线未知坐标。

2、由两个已知点坐标计算支导线坐标（程序名：ZBJS2）
程序清单：
M“X0”：N“Y0”：D“XH”：E“YH” ↲
Pol(M-D，N-E)：Fixm：T=J：X=M：Y=N：T〈0=〉T=T+360⊿T“TH-J=” ◢
I“S=”◢
Lb1 0:{AS}: T=T+A-180：T〈0 =〉T=T+360⊿ T〉360 =〉T =T-360⊿T◢
X=X+Scos T◢
Y=Y+Ssin T◢
Goto 0
程序运行说明：
运行后，显示X0，Y0时输入已知置镜点坐标，显示XH，YH时输入已知后视点坐标，其它同ZBJS的输入与显示。

3、由一个已知点坐标和一个已知方位角进行放射型导线坐标计算（程序名FSXZ BJS）
程序清单：
C“X0”：D“Y0”：E“T0”：Fixm：Lb1 1:{SA}↲
T=E+A-180：T〈0=〉T=T+360⊿ T〉360 =〉T=T-360⊿T◢
X=C+Scos T◢
Y=D+Ssin T◢
Goto 1
程序运行说明：。