fx4800P-2章
CASIO fx-4800p 用户说明书
目錄
使用圓括號的計算 .................................................................... 百分比計算 ............................................................................... 2-2 角度測量單位 ............................................................................... 2-3 三角函數和反三角函數 ................................................................. 2-4 對數函數和指數函數 .................................................................... 2-5 雙曲函數和反雙曲函數 ................................................................. 2-6 其他函數 ....................................................................................... 2-7 坐標變換 ....................................................................................... 2-8 排列和組合 ................................................................................... 2-9 分數 .............................................................................................. 2-10 工學記法計算 ............................................................................... 2-11 小數位數﹑有效位數﹑及指數記法 .............................................. 2-12 使用存儲器的計算 ........................................................................ 獨立存儲器 ............................................................................... 變量存儲器 ............................................................................... 49 50 50 51 51 52 53 54 55 56 57 58 59 59 59
卡西欧4800P计算器程序
环运算语句J 子程序名称H=X+DcosG◣ 公式运算、数据显示语I=Y+DsinG◣ 公式运算、数据显示语T=X—EcosG◣ 公式运算、数据显示语U=Y—EsinG◣ 公式运算、数据显示语最后计算器状态设定语句是大家最容易忽视的。
如果将单位进行预设那么计算器就会默认其使用单位在进行下一单位换算时要一定要进行单位转换，否者会使计算结果错误。
在显示屏幕的左下角可以清楚地发现小提示符号：如D代表度为现在的缺省单位、R代表弧度为现在的缺省单位、G 代表梯度为现在的缺省单位。
首页┊ GPS测量┊工程测量┊遥感测量┊地籍测量┊控制测量┊数字测图┊地图制图┊数据处理┊测绘仪器┊软件应用┊房产测量┊地理信息┊测量管理┊程序开发┊理论文章┊其他┊ 文章更新列表┊游客发布论文所属分类：程序开发来源：作者：更新日期：2006-8-25 阅读次数：156 文章更新列表查看公告列表00-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＞501▲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）Z、待求点的桩号。
K、待求点在左幅输0右幅输1。
B、待求点距中距离。
Goto 5前的桩号为起止点或断链点桩号。
P T÷100T-Int T）÷.6X=Q+Rec（L,T）:Y=U+J:“X=”:X:Pause 0:“Y=”:Y▲Lbl 1:Prog“55”:Prog“YS”:Z＞48980 oto 4△Prog“FQ”: Goto ZLbl 4:Prog“57-1”:Prog“YS”:Z＞50354.96=＞Goto 5△Prog“FQ”Z＞50167.16 =＞Z og“FQ”: Goto ZLbl 7:Prog“60”:Prog“YS”:Z＞52546.8=＞Goto 8△Prog“FQ”: Goto ZLbl 8:Prog“61”:Pro B△Prog“FQ”: Goto ZLbl B:“CX ERROR”▲ Goto 0Lbl Z:Z=Intg（1000Z+.5）÷1000:B=Intg（1000B+.5）÷100号为第2个曲线组合的ZH点桩号，以后类推。
“第1章 绪论”CASIOfx-4800P计算程序
“第1章 绪论”CASIOfx-4800P 计算程序§1.1 参考椭球元素计算测量中，高程的基准面是大地水准面，地表某点距离大地水准面的垂直距离称为该点的高程H ；大地坐标系的基准是参考椭球，用大地经度L 和大地纬度B 表示点在参考椭球表面的位置，用大地方位角A 表示球面上边长的方向。
定义参考椭球大小的元素值为：a ——长半轴b ——短半轴f ——扁率 aba f −=e ——第一偏心率 2222a b a e −=e ′ ——第二偏心率 2222b b a e −=c ——极曲率半径 ba c 2=我国采用过的两个参考椭球元素值以及GPS 测量使用的参考椭球元素值列于表1-1。
表1-1 参考椭球元素值坐标系名称 1954北京坐标系 1980西安坐标系 WGS-84坐标系(GPS 用)参考椭球名称 克拉索夫斯基椭球 IUGG —1975椭球 IUGG —1979椭球a 6378245 6378140 6378137b6356863.01877 6356755.28815 6356752.31424066 f 1:298.3 1:298.257 1:298.2572235632e 0.006693421622966 0.006694384999588 0.006694379990132e ′ 0.006738525414684 0.006739501819473 0.00673949674227 c6399698.90178 6399596.65199 6399593.62601(1) 数学模型通过参考椭球面上任一点的法线，可以作无数个法截面，法截面与参考椭球面的交线称为法截线，球面上不同纬度及方向法截弧线的曲率半径是不相同的。
1) 子午圈曲率半径M通过参考椭球旋转轴南北极的法截面(也称子午面)与参考椭球面相截形成的闭合圈称为子午圈，其曲率半径M 的计算公式为：3Vc M =(1-1) 式中21η+=V(1-2) B e cos ′=η(1-3)2) 卯酉圈曲率半径N与该点子午面相垂直的法截面与参考椭球面相截形成的闭合圈称为卯酉圈，其曲率半径N 的计算公式为：Vc N =(1-4)3) 平均曲率半径R2V c MN R == (1-5)由于V >1，所以有M R N >>。
卡西欧计算器程序fx-4800p
Casio Fx-4800P计算器测量计算程序抚顺市市政建设(集团)有限公司第一工程处：张和光Casio Fx-4800P一、整机概况Fx-4800P计算器是日本Casio公司的产品，它除了具备函数型计算器的全部功能之外，还具备如下特殊功能：1、二进制和八进制字组显示；2、二进制、八进制、十进制和十六进制的转换和加、减、乘、除运算。
3、进行f（x）函数的输入与积分的计算；4、进行标准差和回归计算。
Fx-4800P计算器具有同时显示运算表达式和计算结果的四行式显示的特点。
该计算器具有记存4500个程序步的记忆储存装置，计算器储存了程序之后，只要您输入数据，再按一下EXE键，计算器就会用数据进行程序执行,对于随数据而变化的反复计算的运用十分方便。
Fx-4800P计算器具备条件指令转移和非条件指令转移的能力，还具备逻辑判断能力，因此它可以进行测量工作中比较复杂的计算，当然其它专业的计算工作也同样可以解决。
（1）Fx-4800P计算器操作方法比较简单、易学，程序编排灵活、易懂，可称之谓各专业技术人员的物美、价廉、得心应手的计算工具。
每一台Fx-4800P计算器都具备一本说明书，该计算器各部键的功能以及操作方法，只要认真按照说明书各章节学习，很快就能掌握和运用，现将我在生产实践中应用的部分计算程序介绍给测绘界的同行以便共同探讨。
（2）二、程序库1．坐标反算：程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C ,F-D:I”S1-2=”◢J≤O=>J=J+360△J”A1-2=”操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→Y2→EXE→EXE→EXES1-2：计算得出的距离；A1-2：计算得出的角度。
（按shift°′″转换为60进制的角度）注：此程序可循环计算。
（3）2．坐标正算：程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X 2”=C+Rec(S,R◢Y” Y2”=D+J操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离→EXE→输入A1-2角度（例125°31′23.25″）→EXE→EXE注：此程序可循环计算。
CASIO_fx4800P与fx4850P说明书
次并显示计算结果，按 键删除显示后，计算结果不被保存；而 fx-4800P 可以完整地显示用户输
入的表达式，按 键显示的计算结果被保存在 Ans 存储器中，按 键删除显示后，可以按 键调出最近一次表达式的计算结果。还可以在完成表达式的计算后，立即按重演键 或 调出
并显示最近执行的表达式，再应用 或 键移动光标、输入新字符覆盖当前光标处的字符、应
第 1 章 CASIO fx-4800P 与 fx-4850P 的操作方法与编程技巧
§1.1 fx-4800P 的特点
众所周知，CASIOfx-4500PA 可编程函 数计算器早在 80 年代末就已经在我国工程 界 得 到 了 广 泛 的 应 用 ， fx-4800P 是 比 fx-4500PA 更高一个级别的产品(见图 1-1)， 它采用硬翻盖设计，除了具有 fx-4500PA 的 全部功能以外，还主要改进了下列特性：
30，但π 2 不可以输入为 2，4 3 显然也不可以输入为 43。
③ 计算器遵循括号内的表达式、A、B 型函数，×、÷、＋、－的运算规则，而在 A、B 型函
数内的表达式，省略了×号的乘法运算级别最高。对于有同级别的运算，遵循由左至右的运算规则。
如输入表达式 2π sin 30 最简单的按键方式为 2
。
③ 2500 加上其 15%(答案为 2875)。按键：2500 15
。
④ 3500 减去其 25%(答案为 2625)。按键：3500 25
。
⑤ 168、98 及 734 的和减去其 20%(答案为 800)。
按键操作：168 98 734
20
。
⑥ 原重量为 543 克，将其重量加上 345 克后，总重量为原重量的百分之几(答案为 163.54%)?
4800P放样程序(2)
1.)TYQXJS(程序名)Defm 2:“1.SZ=>XY”: “1.XY=>SZ”:N:N=1=>Goto 1≠>Goto 2⊿←Lb1 1:｛SZ｝:SZ←S≤每段线终点桩号=>J=1：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=2：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=3：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=4：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=5：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=6：Prog“DAT1”: Goto 3⊿←S≤每段线终点桩号=>J=7：Prog“DAT1”: Goto 3⊿←Lb1 3:W=Abs(S-O): Prog“SUB1”:X“XS”=X▲Y“YS”=Y▲Z[5]=X：Z[6]=Y：F“FS”=F-90▲Lb12:Prog“DAT1”: ｛XY｝:XY:I=X:J=Y:Prog“SUB2”:S“S”=O+W▲Z“Z”=Z▲Goto 22.)SUB1(程序名)A=0.1739274226:B=0.3260725774:K=0.0694318442:L=0.3300094782:F=1-L:Z[1]=1-K:X=U+W(A COS(G+QEKW(C+KWD))+B COS(G+QELW(C+LWD))+B COS(G+QEFW(C+FWD))+A COS(G+QE Z[1]W(C+Z[1]WD))):Y=V+W(A SIN(G+QEKW(C+KWD))+B SIN(G+QELW(C+L WD))+B SIN(G+QEFW(C+FWD))+A SIN(G+QEZ[1]W(C+Z[1]WD))):M=90:F=G+QEW(C+WD)+M:X=X+Z COS F:Y=Y+Z SIN F←3.)SUB2(程序名)M=90:T=G-M:W=Abs((Y-V)COS T-(X-U)SIN T):Z=0:Lb10:Prog“SUB1”:L=T+QEW(C+WD):Z= (J-Y)COS L-(I-X) SIN L: Abs Z＜1E-6=> Goto 1:≠>W=W+Z: Goto 0⊿←Lb1 1:Z=0: Prog“SUB1”:Z=(J-Y)÷SIN F4.)DAT1(程序名)J=1=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=2=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=3=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=4=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=5=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=6=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿J=7=>U=起点X坐标：V=起点Y坐标：O=起点桩号：G=起点方位角：H=线形长度：P=XXXXXX ：R=XXXXXX：Q=XXXX⊿……………………………………………………U“X0”：V“Y0”：O“S0”：G“F0”：H“LS”：P“R0”：R“RN”：Q：C=1÷P：D=（P-R）÷（2HPR）：E=180÷∏←说明：1．Q-------以道路中线前进方向，直线输0，左拐曲线输-1，右拐曲线输12．P------ 起点曲率R----- 终点曲S-----桩号Z----边桩距离左为负，右为正，0为中线坐标，N输入1时为正算，输入2时为反算；M为正交或斜交，正交输90，斜交输设计数M在程序内容里改3．当线段为直线时，起点、止点的曲率半径输入1E45；当线段为圆曲线时，起点、止点的曲率半径均输入圆的半径；4．当线段为完整缓和曲线时，起点与直线相接时，起点曲率半径输1E45，与圆曲线相接时，起点曲线半径等于圆曲线半径；止点与直线相接时，止点曲率半径输1E45，与圆曲线相接时，止点曲率半径等于圆曲线半径5．当线段为非完整缓和曲线时，起点与直线相接时，起点曲率半径等于设计规定值，与圆曲线相接时，曲线半径等于圆曲线半径；止点与直线相接时，终点曲率半径等于设计规定值，与圆曲线相接时，曲线半径等于圆曲线半径6．本程序可以连续计算整条路的坐标，但要把整条路的要素输入计算机，也可以计算一部分。
FX4800P测量程序集
CASIO fx 系列编程计算器测量程序集编写: 王劲松第一部分：CASIO fx-4800P 编程计算器测量程序集一．SHI JU CE TU （视距测图）1． 用途：该程序用于“经纬仪视距法碎部测量”时，计算测站点至碎部点间的平距及碎部点的高程。
2． 程序清单：C “H0”：I ：A “UP ”：B “DOWN ”：V “MIDDLE ”：L ：E=90-L ↙D=100Abs （A-B ）（cosE ）2◢ H=（C+I ）+DtanE-V3． 程序说明：H0——测站点高程；I ——测站仪器高；UP ——塔尺上丝读数；DOWN ——塔尺下丝读数；MIDDLE ——塔尺中丝读数；L ——经纬仪竖盘读数盘左读数；D ——测站点至碎部点间平距；H ——碎部点的高程。
5．该程序若在fx-3950中的程序清单为：？→M ：？→A ：？→B ：？→C ：？→D ：100（A-B ）（cos(90-D)）2→X ◢ M+Xtan(90-D)-C →Y说明：M ——将测站高程H0+仪器高I ；A ——上丝读数；B ——下丝读数；C ——中丝读数；D ——竖盘盘左读数L 。
X ——测站点至碎部点的平距；Y ——碎部点的高程。
二．FAN SUAN （反算）1．用途：该程序是“坐标反算公式”程序，用于根据两点1、2的平面坐标（X1，Y1）、（X2，Y2），计算两点间的距离D12及坐标方位角A12。
2． 程序清单：G “X1”：B “Y1”：C “X2”：H “Y2”：E=C-G ：F=H-B ↙Pol （E ，F ）：I “D12”◢ ∆+=⇒<3600J J J J “A12”3． 程序说明：X1，Y1——点1的平面坐标；X2，Y2——点2的平面坐标； D12——点1至点2间平距；A12——点1至点2坐标方位角。
由于编程时用到了极坐标函数POL （），故在程序运行中，会出现I=？XXX （数字）及J=？XXX （数字）这是函数POL计算出的两个量，不能改变，直接回车默认即可。
卡西欧 fx-4800 P计算器面积计算程序及其应用
Ａｒ ｅ ａ Ｃａ ｌ ｃ ｕｌ ａ ｔ ｉ ｏ ｎ Ｐｒ ｏ ｃ ｅ ｄｕｒ ｅ ａ ｎ ｄ Ａｐｐ ｌ ｉ ｃ ａ ｔ ｉ ｏ ｎ ｏ ｆ Ｃａ ｓ ｉ ｏ ｆ ｘ － ４８ ０ ０ Ｐ Ｃａ ｌ ｃ ｕｌ ａ ｔ ｏ ｒ
卡 西欧 ｆ ｘ一 ４ ８ ０ ０ Ｐ计 算 器 面 积 计 算 程 序 及 其 应 用
何 如 伟
（ 石屏 县 采 伐 林 场 ， 云 南 石屏 ６ ６ ２ ２ ０ １ ）
摘要： 在 Ｒ Ｔ Ｋ—Ｇ Ｐ Ｓ系统 、 全站 仪 未普 及 ， Ｇ Ｐ Ｓ又达 不 到 小 面积 测 量精 度 的 当下 ， 计 算 器 配合 罗盘 仪、 测 绳进 行 小面积 测量 还具 有一 定意 义 。文章根 据 闭合 导线 、 极 坐标 、 三 角形 等 多边形测 量原 理 ，
第３ ９卷 第 １期 ２ ０ １ ４年 ２月
ｄ ｏ ｉ ： １ ０ ． ３ ９ ６ ９ ／ ｊ ． ｉ ｓ ｓ ｎ ． １ ６ ７ １ ・ ３ １ ６ ８ ． ２ ０ １ ４ ． ０ ２ ． ０ ０ ３
林
业
调
查
规
划
Ｖ０ １ ． ３ ９ ＮＯ ． １ Ｆｅ ｂ． ２ ０ １ ４
Ｆ ｏ ｒ ｅ ｓ ｔ Ｉ ｎｖ ｅ ｎ ｔ ｏ ｒ ｙ ａ ｎ ｄ Ｐｌ ａ ｎ ｎ ｉ ｎ ｇ
编 制 了卡 西欧 ｆ ｘ一 ４ ８ ０ ０ Ｐ与 ４ ８ ５ ０ Ｐ计 算 器面 积计 算 程 序 ， 并 对其 应 用 方 法进 行 阐述 。 应 用 结果 表 明， 程序 简单 实用 ， 不懂程 序 的人 按 输入代 号 和操作 步骤 就会 计 算 。 关键 词 ： 卡 西欧 ｆ ｘ 一 ４ ８ ０ ０ Ｐ计 算 器 ； 面积计 算程序 ； 闭合导 线测 量 ； 极 坐标 测量
卡西欧计算器程序fx-4800p
Casio Fx-4800P计算器测量计算程序抚顺市市政建设(集团)有限公司第一工程处：张和光Casio Fx-4800P一、整机概况Fx-4800P计算器是日本Casio公司的产品，它除了具备函数型计算器的全部功能之外，还具备如下特殊功能：1、二进制和八进制字组显示；2、二进制、八进制、十进制和十六进制的转换和加、减、乘、除运算。
3、进行f（x）函数的输入与积分的计算；4、进行标准差和回归计算。
Fx-4800P计算器具有同时显示运算表达式和计算结果的四行式显示的特点。
该计算器具有记存4500个程序步的记忆储存装置，计算器储存了程序之后，只要您输入数据，再按一下EXE键，计算器就会用数据进行程序执行,对于随数据而变化的反复计算的运用十分方便。
Fx-4800P计算器具备条件指令转移和非条件指令转移的能力，还具备逻辑判断能力，因此它可以进行测量工作中比较复杂的计算，当然其它专业的计算工作也同样可以解决。
（1）Fx-4800P计算器操作方法比较简单、易学，程序编排灵活、易懂，可称之谓各专业技术人员的物美、价廉、得心应手的计算工具。
每一台Fx-4800P计算器都具备一本说明书，该计算器各部键的功能以及操作方法，只要认真按照说明书各章节学习，很快就能掌握和运用，现将我在生产实践中应用的部分计算程序介绍给测绘界的同行以便共同探讨。
（2）二、程序库1．坐标反算：程序步骤：C”X1”:D”Y1”:E”X2”:F”Y2”:Fixm:Pol(E-C ,F-D:I”S1-2=”◢J≤O=>J=J+360△J”A1-2=”操作过程：ZBFS→EXE→输入X1值→EXE→输入Y1值→EXE→输入X2值→EXE→Y2→EXE→EXE→EXES1-2：计算得出的距离；A1-2：计算得出的角度。
（按shift°′″转换为60进制的角度）注：此程序可循环计算。
（3）2．坐标正算：程序步骤：C”X1”:D”Y1”:S”S1-2”:R”A1-2”:Fixm:X”X 2”=C+Rec(S,R◢Y” Y2”=D+J操作过程：ZBZS→EXE→输入X1值→EXE→输入Y1值→EXE→输入S1-2距离→EXE→输入A1-2角度（例125°31′23.25″）→EXE→EXE注：此程序可循环计算。
fx4800编程计算器程序二
6. 圆曲线段 文件名：YQX (COMP) 程式： L”K0”:O”X0”:P”Y0”:W”ALF”: Lbl 1:{K}: J=(-1)^N (K-L)÷R×180÷π:D=2Rsin((-1)^N J÷2): X”XZ”=O+Dcos(W+J÷2) ◢ Y”YZ”=P+Dsin(W+J÷2) ◢ {B}:S”XL”=X-B”LL”cos(W+J+Q) ◢ T”YL”=Y-Bsin(W+J+Q) ◢ {C}:U”XR”=X+C”LR”cos(W+J+Q) ◢ V”YR”=Y+Csin(W+J+Q) ◢ Goto 1 注：若没有直接 HY 点方位角，则 ALFHY=ALFZH±Ls/2/R×180/π,(左转-，右转 +)。
3.缓和曲线曲线要素： 程序步骤： B=0°1718.87′﹡L/R◢ X=L-L＾3/40/ R2◢ Y=L2/6/ R◢ P=Y-R(1-cos B◢ Q= X-Rsin B◢ T=(R+P)tan(A/2)+Q◢ E=(R+P)(cos(A/2))-1-R◢ Z= R(A-2B)л/180+2L◢ J=2T-Z◢ D=X-Y/tan B 操作过程： HQQXYS→EXE→输入 L 值（即缓和曲线总长）→EXE→输入 R 值（即圆曲线半径） →EXE→得β角度→EXE→得 Xh 值→EXE→得 Yh 值→EXE→得 P 值→EXE→得 Q 值 →EXE→输入 A 角（例 125°31′23.25″）→EXE→得 T 值→EXE→得 E 值→EXE →得 Z 值→EXE→得 J 值→EXE→得 D 值 注：此程序可循环计算。
CASIO-fx4800P型计算器连续计算线路高程、坐标
CASIO-fx4800P型计算器连续计算线路高程、坐标及放样程序【简述】为了适应现场快速、准确、灵活的放样要求，作者根据线路施工测量特点，运用计算器的编程功能，使整条线路的施工测量计算数据能连续快速的计算，只要把整条线路各个交点处的计算要素输入子程序中，你便拥有整条线路的“数据库”。
在忙碌的施工现场，你不必携带大量的施工图和资料来查找平面、纵断面的设计要素，更不必担心设计要素的输入错误，从而使得计算和放样轻松、方便、快捷。
【计算范围】坐标部分包括：各等级公路和高速路的直线段、圆曲线段、加对称缓和曲线的圆曲线段中、边桩坐标及放样计算。
高程部分包括：各等级公路和高速路的直线段、竖曲线段、缓和超高段中、边桩高程及放样计算。
一、坐标及放样程序【起算数据】：交点里程桩号、交点坐标、前直线方位角、交点转角、圆曲线半径、缓和曲线长（一）主程序：XYFY2Fix4M"Xc=":P"Yc=" ＇如须放样输入测站坐标值,否则输入0Lbl0{HS}:H"K0+0":S"B=(-B,0,B)":Prog"A" ＇变量输入和声明S=0=>E=0:≠>E=90 ＇求边桩坐标时的偏角值判定L=πRN÷180+V ＇含有缓曲线的曲线总长T=(V÷2-V^3÷(240R^2))+(R+(V^2÷(24R)-V^4÷(2688R^3)))tan(N÷2) ＇切线长A=Q-T:B=A+V:D=A+L:C=D-V ＇A-ZH(ZY)点桩号,B-HY点桩号,C-YH点桩号,D-HZ(YZ)点桩号Rec(T,F+180) ＇求ZH(ZY)点的坐标增量Z[1]=W+I:Z[2]=K+J ＇ZH(ZY)点坐标值Rec(T,F+GN) ＇求HZ(YZ)点的坐标增量Z[3]=W+I:Z[4]=K+J ＇HZ(YZ)点的坐标值Lbl1U=S÷(Abs(S+10E-19)):S=AbsS ＇左右偏角每件判定H<A=>Goto2: ≠> ＇第一直线上坐标计算判定H<B=>Goto3: ≠> ＇第一缓和曲线上坐标计算判定H<C=>Goto4: ≠> ＇净圆曲线上坐标计算判定H<D=>Goto5: ≠> Goto6⊿⊿⊿⊿ ＇第二缓和曲线和第二直线上坐标计算判定Lbl2Rec(Q-H,F+180) ＇第一直线上中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标值计算Rec(S,F+180-(180-E)U) ＇第一直线上与中桩同一断面的边桩增量计算X=X+I:Y=Y+J ＇边桩坐标值计算S=SU ＇边长条件判定Goto8Lbl3 ＇进入第一缓曲线段计算Z=H-A ＇待求点至ZH点的距离O=90Z^2÷(πRV) ＇所求点缓曲线对应的圆心角X=Z-Z^5÷(40R^2V^2)+Z^9÷(3456R^4V^4) ＇缓曲线上ZH点至待求的支距X坐标Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇缓曲线上ZH点至待求点的支距Y坐标Lbl7Rec(X,F)X=Z[1]+I:Y=Z[2]+JRec(Z,F+90G) ＇把中桩支距坐标转换成测量统一坐标X=X+I:Y=Y+JRec(S,F+O G+EU) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SU ＇边长条件Goto8Lbl4 ＇进入主圆计算部分Z=H-A-V÷2 ＇曲线上待求点到HY点的里程O=180Z÷(Rπ) ＇待求点方位角所对应的圆心角X=Rsin O+(V÷2-V^3÷(240R^2))Z=R(1-cos O)+(V^2÷(24R)-V^4÷(2688R^3)) ＇主圆部分待求点支距坐标XY Goto7Lbl5 ＇进入第二缓和曲线计算Z=D-H ＇曲线上待求点到HZ点的里程O=90Z^2÷(πRV) ＇待求点方位角所对应的圆心角X=Z-Z^5÷(40Z R^2V^2)+^9÷(3456R^4V^4)Z=Z^3÷(6RV)-Z^7÷(336R^3V^3)+Z^11÷(42240R^5V^5) ＇待求点支距坐标XY Rec(X,F+GN+180)X=Z[3]+I:Y=Z[4]+J ＇把中桩支距坐标转换成测量统一坐标Rec(Z,F+GN+180-90G)X=X+I:Y=Y+JRec(S,F+GN+180-O G-(180-E)U) ＇曲线边桩坐标增量计算X=X+I:Y=Y+J ＇坐标计算结果S=SUGoto8Lbl6 ＇进入第二直线段计算部分Rec(H-D+T,F+GN) ＇中桩坐标增量计算X=W+I:Y=K+J ＇中桩坐标计算结果Rec(S,F+GN+EU) ＇边桩坐标增量计算X=X+I:Y=Y+J ＇边桩坐标计算结果S=SUGoto8Lbl8X=X◢ ＇显示坐标计算结果Y=Y◢M≠0=>Goto9: ≠>Goto0 ⊿Lbl9 ＇进入放样部分计算Pol(X-M,Y-P) ＇增量计算J<0=>J=J+360: ≠>J=J ⊿J"A"=Int J+0.01Int(60Frac J)+0.006Frac(60Frac J) ◢ ＇方位角计算结果（此显示值为度分秒格式）I"L"=I◢ ＇极距计算结果Goto0∙(二)子程序:AH≤ZY2(ZH2)=>Q=JD1:W=X1:K=Y1:R=R1:F=F0:N=N1:V=Ls1:G=1(-1):=>H≤ZY3(ZH3)=>Q=JD2:W=X2:K=Y2:R=R2:F=F1:N=N2:V=Ls2:G=1(-1):=>H≤ZY4(ZH4)=>Q=JD3:W=X3:K=Y3:R=R3:F=F2:N=N3:V=Ls3:G=1(-1):=>H≤ZY5(ZH5)=>Q=JD4:W=X4:K=Y4:R=R4:F=F3:N=N4:V=Ls4:G=1(-1):=>……………………………H≤ZY n+1(ZH n+1)=>Q=JD n:W=X n:K=Y n:R=R n:F=F n-1:N=N n:V=Ls n:G=1(-1)注:H—待求点桩号ZY n+1(ZH n+1)—后曲线起点桩号 Q—本曲线交点桩号 W—本曲线交点X坐标K—本曲线交点Y坐标 R—本曲线内圆曲线半径 F—前直线方位角 N—本交点转角V—本曲线内缓和曲线长 G—线路转向符,左转角-1,右转角为+1二、高程及放样程序【起算数据】：变坡点桩号A、变坡点高程B、前后纵坡C D、竖曲半径R、横坡度F、最大超高横坡度G、ZH点桩号、YH点桩号、HZ点桩号、施工层厚度N、视线高W （一）主程序：GCFYW"SXG="：N ＇常量输入（视线高程、施工层厚度）Lbl 1｛LО｝:L"K0+0"：О"B=" :Prog"B" ＇变量输入及声明K=C-D ＇坡度差及竖曲线偏角K>0=>J=-1:≠>J=1⊿ ＇凹凸竖曲线判断条件T=Abs(RK)÷2 ＇竖曲线切线长,如需显示则在后面加上◢H=A-T ＇竖曲线起点ZY桩号,如需显示则在后面加上◢X=A+T ＇竖曲线终点YZ桩号,如需显示则在后面加上◢M=L-A ＇待求点至变坡点的距离L≤H=>P=B+MC-N-W: ≠> ＇前纵坡直线段中桩高程计算L≤A=>P=B+MC+J(T-AbsM)2÷(2R)-N-W: ≠> ＇前纵坡竖曲线中桩高程计算L≤X=>P=B+MD+J(T-AbsM)2÷(2R)-N-W: ≠> ＇后纵坡竖曲线中桩高程计算L>X=>P=B+MD-N-W⊿⊿⊿⊿⊿ ＇后纵坡直线段中桩高程计算E=0=>Goto 2 ＇直线部分中、边桩计算判断条件≠>L≤E=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>V=> Goto 2 ＇直线部分中、边桩计算判断条件≠>L>E=> Goto 3 ＇前缓曲超高部分中、边桩计算判断条件≠> L>U=> Goto 4⊿⊿⊿⊿ ＇后缓曲超高部分中、边桩计算判断条件Lbl 2O=0=>Z“H”=AbsP◢＇显示直线部分中、边桩高程计算结果≠>Y“HB”=Abs(P-OF)◢不需放样显示高程,否则显示放样水准尺计数⊿Goto 1Lbl 3S=（L-E）（F+G）÷（V-U）-F ＇前缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至零界面、零界面至全超高及全超高≠> Goto 7 ⊿⊿ 部分的计算判断条件Lbl 4S=（V-L）（F+G）÷（V-U）-F ＇后缓曲线超高横坡度计算S≤F=> Goto 5：≠>S≤G=> Goto 6 ＇超高起点(终点)至超高变坡零界面、零界面至全≠> Goto 7 ⊿⊿ 超高、全超高部分计算判断条件Lbl 5O=0=>Z“H1”=AbsP◢＇超高起点(终点)至超高变坡零界面中、边桩高程≠>Y“H-D1”=Abs(P-OF)◢计算结果显示Q“H-G1”=Abs(P+OS)◢⊿Goto 1Lbl 6O=0=>Z“H2”=AbsP◢＇超高变坡零界面至全超高中、边桩高程≠>Y“H-D2”=Abs(P-OS)◢计算结果显示Q“H-G2”=Abs(P+OS)◢⊿Goto 1Lbl 7O=0=>Z“H3”=AbsP◢＇全超高中、边桩高程计算结果显示≠>Y“H-D3”=Abs(P-OG)◢Q“H-G3”=Abs(P+OG)◢⊿Goto 1∙∙注：W“SXG”-视线高N-施工层厚度L"K0+0"-待求点桩号О"B="-中桩至边桩距离（二）子程序：BL≤ZY2=>A=JD1:B=H1:C=I0:D=I1: R=R1:F=I h:E=ZH1:U=YH1:V=HZ1:G=I max1 :≠>L≤ZY3=>A=JD2:B=H2:C=I1:D=I2: R=R2:F=I h:E=ZH2:U=YH2:V=HZ2:G=I max2 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠>L≤ZY4=>A=JD3:B=H3:C=I2:D=I3: R=R3:F=I h:E=ZH3:U=YH3:V=HZ3:G=I max3 :≠> ……………………L≤ZY n+1=>A=JD n:B=H n:C=I n-1:D=I n: R=R n:F=I h:E=ZH n:U=YH n:V=HZ n:G=I maxn注:ZY N+1:-下一个竖曲线起点桩号（如：变坡点2的计算范围为前竖曲线终点YZ1至后竖曲线起点ZY3，则在程序中“L≤ZY3”就应输入ZY3的里程桩。
卡西欧4800计算器测量放样程序汇编
卡西欧4800计算器测量放样程序汇编4800程序使用说明书程序代码术语：连续曲线数据块、独立曲线数据块、长短链功能：1、放样对应里程中心点2、放样对应里程旋转任意角度外移点（图A点）3、放样对应里程旋转任意角度外移点的外移点（图B点）4、获取放样点坐标及对应里程中心点沿线路前进方向切线方位角程序输入注意事项：1、常用程序文件放在前面，便于调用，一般先建立“CIRCLE”、“XY-AL”、“AL-XY”三文件，文件名千万不能输错。
2、加粗带下划线词组为计算器中的函数或符号，只能从计算器中调出此函数或符号，不能从键盘输入。
符号←表示回车键“EXE”坐标输入通用格式：屏幕显示OPP NAME?（置镜点名称）或CPP NAME?（后视点名称）：1：输入-1：手动输入坐标X,Y2：输入大于1数字：输入用数字所代表的控制点名。
如数据库中没有本数字所代表的点名，则系统显示SynERROR in…（文件名）。
曲线放样程序（CIRCLE）运行说明：步骤1：屏幕提示“LOAD?”（输入曲线要素）1.1：输入1（默认）：输入里程后将自动判断所在曲线并调用（此功能必须保证曲线数据库中有连续曲线数据块，具体详见数据库建立）。
1.2：输入-1：屏幕提示手动输入各曲线要素，ZJ（转角）、R（半径）、L0(缓和曲线长)、ZH KM(直缓里程)、JDX(曲线交点X坐标)、JDY(曲线交点Y坐标)、A0(起始直线边前进方向方位角)、1 R,-1 L(右偏输入1，左偏输入-1)。
1.3：输入大于1的数字：输入用数字代表的曲线名称，比如13，表示放样点位于曲线13 上（曲线名称在建立数据库是自己命名，但一定要大于1）。
如库中没有此名称，系统显示“Syn ERROR in CIRCLE”。
（此功能必须保证曲线数据库中有独立曲线数据块）步骤2：屏幕提示“OPP NAME?”（输入置镜点名称）详见坐标输入格式。
步骤3：屏幕提示“DK+M?”（输入放样点对应中心里程）如K15+002.35应输入15002.35。
卡西欧4800P计算器测量程序集-10页文档资料
福建林业职业技术学院工程系土木工程室高国兴CASIO fx-4500PA计算中桩、边桩坐标程序直线、圆曲线中桩及左、右边桩计算程序fx-4800pC“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：L“Lμo=”：K“D=”：M“Jμo=”：Q=9θL÷（πR）：S=2RsinQ ：O =C＋Scos（N-Q）◢P=D+Ssin（N-Q）◢W=E-2Q◢X=O+Kcos（N-2Q-M）◢Y=P+Ksin（N-2Q-M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－直线R为无穷大Lμo=置镜点到待求点距离D左、右边距离Jμo=前进方向左（+90），右（－90）O=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标缓和曲线中桩及左、右边桩计算程序Fx-4800pC“Xθ=”：D“Yθ=”：E“Wθ=”：N=Ans：R“R=”：F“Lθ=”：L“Lμo=”：K“D=”：M“Jμo=”←U=L－L^5÷（40R2F2）：V=L^3÷（6RF）－L^7÷（336R^3F^3）：Q=tan－1（V÷U）：S=√（U2+V2）←B=9θL2÷（πRF）：O=C+Scos（N－Q）◢P=D+Ssin（N－Q）◢W=E－B◢X=O+Kcos（N－B－M）◢Y=P+Ksin（N－B－M）◢注：Xθ=起算点坐标Yθ=起算点坐标Wθ=起算点方位角R=曲线半径，线路前进方向，左+，右－Lθ=缓和曲线长Lμo=置镜点到待求点距离D左、右待求距离Jμo=前进方向左（+90），右（－90）O=中桩坐标P=中桩坐标W=待求点方位角X=左、右边桩坐标Y=左、右边桩坐标放样程序fx-4800pRLI“XZ”：W：“YZ”：G“XH”：H“YH”←B=H－W ←C=G－I ←Prog“P”：U=A←K=√（B2＋C2）←C=X“X－”－I←B=Y“Y－”－W←Prog“P”←L“L”=√（B2＋C2）◢V=A－U：V<θ==>V=V＋36θ⊿V “V=”◢PA=tan－1（B÷C）←B≥θ==>C≥θ==>A=A：≠=>C<θ==>A=A＋18θ⊿⊿≠=>B<θ==> C≥θ==>A=A＋36θ: ≠=>C<θ==>A= A＋18θ⊿⊿⊿←XZ、YZ—置镜点X、Y坐标；XH、YH—后视点坐标；X-、Y-—前视点坐标；输出：R=水平角、L=置镜点至前视点水平距离。
正式-CASIO-fx-4800P计算器的使用-基础知识上课讲义
表1-3 功能菜单的意义
§1.2 CASIO fx-4800P的基本操作
(4) fx-4800P与普通科学计算器常规操作的主要区别
1) 键与 键 在普通科学计算器中，
键的功能是计算前
已输入的表达式；而对fx-4800P，在COMP模式
下，
键的功能是计算前已输入的表达式，在
PROG的编程模式下，
用户输入的表达式，按 键显示的计算结果被保存在Ans存储
器中，按 键删除显示后，可以按
键调出最近一次
表达式的计算结果。还可以在完成表达式的计算后，立即按重演
键 或 调出并显示最近执行的表达式，再应用 或
键移动光标、直接输入新的字符覆盖当前光标处的字符、应用
键删除字符、应用
键插入输入的字符可以编辑当前显示
模式菜单
§1.2 CASIO fx-4800P的基本操作
(1) 键盘区简介
3) 第三键盘区有4行5列共20个键，其键面功能主要是 数字和＋、－、×、÷四则运算。
每个按键一般有键面字符、键上部1~2个字符三种功能。
直接按一个键则输入该键面白色字符的功能；先按
键，再按一个键，则输入该键左上方橘色字符的功
§1.2 CASIO fx-4800P的基本操作
(3) 键与 键的使用
按 键，屏幕显示图1-2所示的模式菜单。fx-4800P 共有8种模式，键入模式前的数字就可以选中该模式，各 模式的意义见表1-2。
表1-2 模式菜单的意义
§1.2 CASIO fx-4800P的基本操作
(3) 键与 键的使用
(5) 基本操作
2) A型函数和B型函数
提 示
§1.2 CASIO fx-4800P的基本操作
卡西欧4800程序
支导线程序(fx4800P)ZDXProg”006”:A=J:Z”1—2”:X=U:Y=W:LbI0:{Q}:Q”Q=”:{L}:L”L=”:X=X+LCos(A+Q):Y=Y+LSin(A+ Q): X”X=”▲ Y”Y=”▲A=A+Q+180:Z=2=>X=U:Y=W:A=J△Goto0放样程序(fx4800P)FYProg”006”: LbI0:{X}:X”X=”:{Y}:Y”Y=”: Prog”007” :Goto0竖曲线程序(fx4800P)SHU--QXI“I=”:J”J=”:R”R=”:W=I-J:C=Abs(RW)/1000:T=C/2:E=T^2/(2R): C”L=”▲ T”T=”▲ E”E=”▲W<0=>K=1△W>0=>K=-1△D”BIAN—LI=”:Z”BIAN—H=”:A=Z+(-IT)/1000:LbI0:{L}:L”DK=”:L≤D-T=>H=Z+(D-L)(-I/1000)△L>D-T=>H=A+I(L-(D-T))/1000+K(L-(D-T))^2/(2R) △L≥D+T=>H=Z+J(L-D)/1000△H”H=”▲Goto0已知点坐标求对应里程和距离(fx4800P)1、主程序(DYLC)Defm42: AB”N1=”C”N2=”:R:Pro g”000” :Pro g”009”:E”PYL=”:K”Z—Y”: Z”ZH=”:M”ZH—X”:P”ZH—Y”:F:Lbl0:Prog”XH”:Goto02、子程序(XH)Lbl1:Pro g”QX”:{U}:U”YZ—X=”:{W}:W”YZ—Y=”:I=U-X:J=W-Y:Pol(I,J):J<0=>J=J+360△Z[37]=J-Z[34]:Z[38]=cosZ[37]×I:D=D+Z[38]:AbsZ[38]≥0.001=>Goto1:≠> D”DYLC=”△▲Z[35]=sinZ[37] ×I: Z[35] “FXJL=”▲已知3点坐标求圆心和半径(fx4800P)主程序(3DZB--QYX)A”AX=”B”AY=”C”BX=”D”BY=”E”CX=”F”CY=”: X=((F-B)(D^2-B^2+C^2-A^2)+(D-B)(B^2-F^2+A^2-E^2))/(2(C-A)(F-B)-2(E-A)(D-B)):X”X=”▲Y=((E-A)(C^2-A^2+D^2-B^2)+(C-A)(A^2-E^2+B^2-F^2))/(2(D-B)(E-A)-2(F-B)(C-A)):Y”Y=”▲R=√((A-X)^2+(B-Y)^2):R“R=”▲自由设站(fx4800P)(边角后方交会)HFJHLbI0:{A}:A”H1—X=”:{B}:B”H1—Y=”:{C}:C”H2—X=”:{D}:D”H2—Y=”:I=C-A:J=D-B:POL(I,J):J≤0=>J=J+360⊿ I:”S=”▲Z[42]=J: Prog”006”:Z[42]:”FW—Q=”▲{E}:E”L1=”:{F}:F”L2=”:{G}:G”Q=”:H=Cos-1((I2+ E2-F2)/(2IE)):K=Cos-1((I2+ F2-E2)/(2IF)):L=A+E×Cos(J+H) :M=B+E×Sin(J+H):N=C+F×Cos(J+180-K):O=D+F×Sin(J+180-K):X=(L+N)/2:Y=(M+O)/2:X:”X=”▲Y:”Y=”▲Goto0不对称缓和曲线放样数据计算(fx4800P)(带偏移量+中边桩)1 主程序(--QX—Z-B)Defm42:AB”N1=”C”N2=”:R:Pro g”000”: Pro g”009”:E”PYL=”:K”Z—Y”: Z”ZH=”:M”ZH—X”:P”ZH—Y” :F Lbl0:Pro g”QX”: Prog”006”:Prog”007”:Goto02 子程序(QX){D}:Z[1]=D”DK=”-Z:O=Z+Z[2] :N=(M+Z[11]CosF)+Z[12]Cos(F+KA) :S=(P+Z[11]SinF)+Z[12]Sin (F+KA):{L}:L”BIN—L”:{Q}:Q”BIN—Q”:Z[1]≤0=>Prog”001”:≠>Z[1]≤B=>Prog”002”:≠>Z[1]≤(Z[2]-C)=>Prog”003”:≠>Z[1]≤Z[2]=>Prog”004”:≠>Z[1]>Z[2]=>Porg”005”△△△△△ X”X=”▲ Y”Y=”▲Z[34]<0=>Z[34]=Z[34]+360:≠>Z[34]>360=> Z[34]=Z[34]-360△△Z[42]=Z[34]:Pro g”010”:Z[42]”QX—Q”▲(该点切线方位角)3 子程序（000）曲线要素计算Z[4]=B/2-B^3/(240R^2):Z[6]=B^2/(24R)-B^4/(2688R^3):Z[8]=90B/(πR）:Z[5]=C/2-C^3/(240R^2):Z[7]=C^2/(24R)-C^4/(2688R^3):Z[9]=90C/(πR) :Z[11]=(R+Z[6])tng(A/2)+Z[4]-（Z[6]-Z[7]）/sinA:Z[12]=(R+Z[7])tng(A/2)+Z[5]+(Z[6]-Z[7])/sinA:Z[10]=(R+(Z[6]+Z[7])/2)(1/cos(A/2)) -R:Z[2]=(A-Z[8]-Z[9])Rπ/180+B+C: Z[3]=Z[2]-B-C4 子程序（001）第一切线段计算X=M+Z[1]cosF:Y=P+Z[1]sinF:X=X+Lcos(F+Q):Y=Y+Lsin(F+Q): Z[34]=F:5 子程序（002）第一缓和曲线段坐标计算Z[28]=F:Z[27]=E:Z[13]=90B/（πR）：Z[14]=Z[13]/3：Z[15]=90Z[1]^2/(πRB):Z[16]=Z[14](Z[1]/B)^2 :Z[17]=RB:Z[18]=Z[1]-Z[1]^5/(40Z[17]^2)+Z[1]^9/(3456Z[17]^4):Z[19]=Z[1]^3/(6Z[17])-Z[1 ]^7/(336Z[17]^3)+Z[1]^11/(42240Z[17]^5):Z[20]=√(Z[18]^2+Z[19]^2):X=M+Z[20]cos(Z[28]+KZ[16]) :Y=P+Z[20]sin(Z[28]+KZ[16]):K<0=> Z[27]= Z[27]:≠>K≥0=> Z[27]=-1 Z[27]△△Z[23]= Z[27] (Z[1]/B):Z[21]= Z[28]+KZ[15]+90: X=X+Z[23]cosZ[21]:Y=Y+Z[23]sinZ[21]:Z[22]=Z[28]+KZ[15]+Q : X=X+LcosZ[22]:Y=Y+LsinZ[22]: Z[34]=Z[28]+KZ[15]6 子程序（003）圆曲线段坐标计算Z[27]=E:Z[28]=F:Z[13]=90B/(πR):Z[22]=Z[28]+K(Z[13]/3):Z[17]=RB:Z[18]=B-B^5/(40Z[17]^2)+B^9/ (3456Z[17]^4):Z[19]=B^3/(6Z[17])-B^7(336Z[17]^3+B^11/(42240Z[17]^5):Z[20]=√(Z[18]^2+Z[19]^2) :X=M+Z[20]cosZ[22]:Y=P+Z[20]sinZ[22]:Z[14]=Z[1]-B:Z[15]=180Z[14]/( πR):Z[23]= Z[28]+KZ[13] :Z[21]=2Rsin(Z[15]/2):Z[16]=Z[15]/2:X=X+Z[21]cos(Z[23]+KZ[16]):Y=Y+Z[21]sin(Z[23]+KZ[16]):Z[24]= Z[23]+KZ[15]:K<0=> Z[27]= Z[27]:≠>K≥0=> Z[27]=-1 Z[27]△△X=X+ Z[27]cos(Z[24]+90):Y=Y+ Z[27]sin(Z[24]+90):X=X+Lcos(Z[24]+Q):Y=Y+Lsin(Z[24]+Q): Z[34]＝Z[28]+KZ[13]+KZ[15]7 子程序（004）第二缓和曲线段坐标计算Z[15]=O-D:Z[22]=F+KA+180:Z[27]=E:K>0=>Z[27]=Z[27]:≠>K<0=> Z[27]=-1Z[27]△△Z[25]=Z[27]( Z[15]/C):Z[23]=-K:Z[13]=90C/(πR):Z[14]=Z[13]/3:Z[16]=90Z[15]^2/(πRC):Z[17]=Z[14](Z[15]/C)^2 :Z[18]=RC:Z[19]=Z[15]-Z[15]^5/(40Z[18]^2)+Z[15]^9/(3456Z[18]^4):Z[20]=Z[15]^3/(6Z[18])-Z[15]^7 /(336Z[18]^3+Z[15]^11/(42240Z[18]^5):Z[21]=√(Z[19]^2+Z[20]^2):Z[30]=N+Z[21]cos(Z[22]+Z[23]Z[17]): Z[31]=S+Z[21]sin(Z[22]+Z[23]Z[17]):Z[24]=Z[22]+Z[23]Z[16]+90: Z[32]= Z[30]+Z[25]cosZ[24]: Z[33]= Z[31]+Z[25]sinZ[24]:Z[26]=F+KA+Z[23]Z[16]+ Q : X= Z[32]+LcosZ[26]:Y= Z[33]+LsinZ[26]: Z[34]＝Z[22]+（-1K）Z[16]-1808 子程序(005)计算第二切线段坐标Z[13]=F+KA:Z[14]=D-O:X=N+Z[14]cosZ[13]:Y=S+Z[14]sinZ[13]:X=X+Lcos(Z[13]+Q):Y=Y+Lsin(Z[ 13]+Q): Z[34]=Z[13]9 程序(006)计算后视数据{U}:U”ZX=”:{W}:W”ZY=”:{G}:G”HX=”:{H}:H”HY=”:I=G-U:J=H-W:Pol(I,J): I”S=”▲Prog”008”: Z[42]=J: Prog”010”: Z[42]”F0=”▲10 子程序(007)计算放样数据T=J:Pol(X-U, Y-W): Prog”008”: I”FY--L=”▲V=J-T:V<0=>V=V+360:≠>V>360=>V=V-360△△Z[42]=V: Prog”010”: Z[42]”FY--Q=”▲11 子程序(008)角度判断J<0=>J=J+360:≠>J>360=>J=J-360△△12 子程序（009）曲线要素结果显示Z[11]”T1=”▲ Z[12]”T2=”▲ Z[10]”E=”▲ Z[2]”QX--L=”▲ Z[3]”YQX—L=”▲13 子程序(010)角度的度分秒输出Z[40]=60FracZ[42]:Z[41]=60FracZ[40]:Z[42]=IntZ[42]+IntZ[40]÷100+Z[41]÷10000。
卡西欧fx-4800P编程
一、本人编写整段曲线坐标程序一则。
可作为参考，从ZH-HY-YH-HZ可一次性计算出中心或切点方向或法线方法等坐标，当只有圆曲线时，L0输入0，只计算单圆周曲线。
示意图无法粘贴。
若有问题请回复。
本程序经调试计算准确。
缓和曲线、圆曲线计算程序A“JDX”：B“JDY”：C“R”：D“L0”：J“HR”：F“T”：G“QD”：H“ZD”：I“Z，Y，-，+”：T“E”：I>0=>Z=1：≠=>Z= -1:△LB1 1：{K，L，M}：K“ZHUANG”：L“L1”：M“L2”：K<G+D=>N“Li”=K-G：Y=1：E=J：GOTO 2：≠=>K≤H-D=> GOTO 3：≠=>N=H-K：E=J+I-180：Y=-1: GOTO 2：△△Lb1 2：O“X0”=A-COS E×F：P“Y0”=B-SIN E×F：Q=N-N^5÷（40×C2×D2）：R=N^3÷（6×C×D）：S=TAN-1（R÷Q）×YZ：U=Y Z N2×90÷（ДCD）+E：K ≤G+D=>U=U-90：≠=>U=U+90:△V=U+90：W“X1”=O+COS（E+S）×√（Q2+R2）+COS U×L+COS V×M▲X“Y1”=P+SIN（E+S）×√（Q2+R2）+SIN U×L+SIN V×M▲GOTO 1：Lb1 3：O=J+Z（0.5×Abs I+90）：P“X0”=Cos O×（C+T）+A：Q“Y0”=Sin O×（C+T）+B：R=O-180-Z（（H-G）-2D）×90÷Д÷C：S=（K-（G+D））×180÷Д÷C：U=R+ZS：W“X1”=P+COS U×（C+ZL）+COS（U+Z×90）×M▲X“Y1”=Q+SIN U×（C+ZL）+ SIN（U+Z×90）×M▲GOTO 1说明：JDX：交点X坐标 JDY：交点Y坐标R：半径 L0：缓和曲线长 HR：前进方位角T：切线长 QD：起点桩号 ZD：终点桩号Z，Y，-，+：偏角，左为负，右为正E：外矢距ZHUANG：待求点桩号L1：待求点法线方向距离，左侧为正，右侧为负，求中桩坐标时为0L2：待求点切线方向距离，左侧为正，右侧为负，求中桩坐标时为0本程序用于计算缓和曲线圆曲线缓和曲线缓和曲线为对称型缓和曲线，若无缓和曲线，计算单圆曲线时L0输入0卡西欧4800计算全站仪坐标程序！这个是我月个朋友编写的“卡西欧4800计算全站仪坐标程序”一共是5个程序。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第2章 在土木工程测量中的应用§2.1 参考椭球与高斯投影计算§2.1.1参考椭球元素计算测量中，高程的基准面是大地水准面，地表某点距离大地水准面的垂直距离称为该点的高程H ；大地坐标系的基准是参考椭球，用大地经度L 和大地纬度B 表示点在参考椭球表面的位置，用大地方位角A 表示球面上边长的方向。
定义参考椭球大小的元素值为：a ——长半轴b ——短半轴f ——扁率 aba f -=e ——第一偏心率 2222a b a e -= e ' ——第二偏心率 2222b b a e -=c ——极曲率半径 ba c 2= 我国采用过的两个参考椭球元素值以及GPS 测量使用的参考椭球元素值列于表2-1。
通过参考椭球面上任一点的法线，可以作无数个法截面，法截面与参考椭球面的交线称为法截线，球面上不同纬度及方向法截弧线的曲率半径是不相同的。
1) 子午圈曲率半径M通过参考椭球旋转轴南北极的法截面(也称子午面)与参考椭球面相截形成的闭合圈称为子午圈，其曲率半径M 的计算公式为：3V c M =(2-1) 式中21η+=V(2-2) B e cos '=η(2-3)2) 卯酉圈曲率半径N与该点子午面相垂直的法截面与参考椭球面相截形成的闭合圈称为卯酉圈，其曲率半径N 的计算公式为：Vc N =(2-4)3) 平均曲率半径R2Vc MN R == (2-5)由于V >1，所以有M R N >>。
4) 任意法截弧的曲率半径A R方位角为A 方向法截弧的曲率半径A R 的计算公式为ANR A 22cos 1η+=(2-6)当A =0时，M R A =；当A =90时，N R A =。
(2) 程序与案例分别在1954北京坐标系、1980西安坐标系和WGS-84坐标系中计算参考椭球表面上任意点的M 、N 、R 及A R 。
1)2) 程序程序名：MNR Lab 0 {Z}Z=54C=6399698.90178: E=0.006738525414684 : Z=80C=6399596.65199:E=0.006739501819473 : Z=84C=6399593.62601: E=0.00673949674227 : Goto 0:{BA} U=EcosB: V=(1+U 2)M=C ÷V^3 N=C ÷VR=C ÷V 2 P=N ÷(1+(UcosA)2)3) 操作步骤设=B 22°32′50″，=A 153°19′18″，分别计算其在1954北京坐标系、1980西安坐标系及WGS-84坐标系中的M 、N 、R 及A R 值。
将上述程序以MNR 的文件名输入计算器后，按键 5 2及选择程序MNR ，按键，屏幕提示542232 15319188084§2.1.2边长投影计算使用电磁波测距仪或全站仪观测得到地面某直线的水平距离D 需要经过两次投影才能在高斯平面坐标系中进行计算。
第一次投影是从地面投影到参考椭球面上，第二次投影是从参考椭球面投影到高斯平面上。
(1) 数学模型设投影到参考椭球面上的长度为0S '，则计算公式为 D R h H R h H D S n g m n g m ⎪⎭⎪⎬⎫⎪⎩⎪⎨⎧+-+-='220)( (2-7) 式中，m H 为边长D 高程，g h 为边长所在地区大地水准面对于参考椭球面的高度，n R 为边长中点方向上参考椭球面法截弧的曲率半径。
设参考椭球面上的长度0S '投影到高斯平面上的长度为0S ，则计算公式为 0222200242S R y R y S S m m m '⎭⎬⎫⎩⎨⎧∆++'= (2-8) 式中，m y 为直线两端点近似y 坐标的平均值，y ∆为直线两端点近似y 坐标的增量，m R 为参考椭球在直线中点的平均曲率半径。
当边长不长时，可以取边长中点方向法截弧曲率半径A R 代替n R 和 m R 。
测区内，可以认为g h 为一常数，如果未知，工程计算中可以取g h =0。
(2)、程序与案例分别在1954北京坐标系和1980西安坐标系中计算边长投影到高斯平面上的距离0S 。
1)2) 程序名：DS Lab 0{Z}Z=54C=6399698.90178: E=0.006738525414684:Z=80C=6399596.65199:E=0.006739501819473:Goto 0{BA}U=EcosB: V=(1+U 2):R= C ÷V(1+(UcosA)2) {DHGYW}J=D-((H+G)÷R-(H+G)2÷R 2)D K=(Y+W)÷2: L=Y-WS=J+(K 2÷2÷R 2+L 2÷24÷R 2)J 3) 操作步骤设=B 22°33′04″，=A 107°41′02″，1y =210154m ，2y =213577m ，D =3591.293m ，m H =35.4m ，g h =4.3m ，分别计算其在1954北京坐标系和1980西安坐标系中的高斯平面距离值。
将上述程序以DS 的文件名输入计算器后，按键 5 2及选择程序DS ，按键，屏幕提示及操2233041074102 3591.293§2.1.3高斯投影正算与子午线收敛角的计算将大地经纬度为L 、B 的点投影到高斯平面上并求出其坐标x 、y 称高斯投影正算，令γ为子午线收敛角。
(1) 数学模型⎥⎦⎤⎢⎣⎡+-+++-++=64244222)5861(7201)495(24121m t t m t m Nt X x ηη(2-9) ⎥⎦⎤⎢⎣⎡-++-++-+=522242322)5814185(1201)1(61m t t t m t m N y ηηη(2-10) []523420)2(12)231(60180m t m m t-++++=ηηπγ(2-11)0γ表示以十进制度为单位的子午线收敛角。
式中B t tan =(2-12) 180cos πBlm =(2-13)N 为卯酉圈曲率半径，计算公式见(2-4)，η的计算公式见(2-3)，X 为子午线弧长。
对于克拉索夫斯基椭球：BB B B B B X cos )sin 0039.0 sin 6973.0sin 9238.133sin 7799.32005(8611.1111347530+++-= (2-14)对于IUGG1975椭球：BB B B B B X cos )sin 0039.0 sin 6976.0sin 9602.133sin 8575.32009(0047.1111337530+++-= (2-15)(2)、程序与案例分别在1954北京坐标系和1980西安坐标系中进行高斯投影正算。
1) 变量对照表2) 程序程序名：GSJSZ Lab 0 {ZLB} S=sinBZ=54C=6399698.90178: E=0.006738525414684 :F=111134.8611B-(32005.7799S+133.9238S^3+0.6973S^5 +0.0039S^7)cosB: Goto 1:Z=80C=6399596.65199: E=0.006739501819473 :F=111133.0047B-(32009.8575S+133.9602S^3+0.6976S^5 +0.0039S^7)cosB:Goto 0 :Lab 1 U=EcosB: V=(1+U 2):N= C ÷V T=tanB: M=(cosB)L π÷180 X=F+NT(0.5M 2+24-1(5-T 2+9U 2+4U^4)M^4+720-1(61-58T 2+T^4)M^6) Y=N(M+6-1(1-T 2+U 2)M^3+120-1(5-18T 2+T^4+14U 2-58U 2T 2)M^5) R=T ÷π×(180M+60(1+3U 2+2U^4)M^3+12(2-T 2)M^5)提示常数π后的×不能省略。
3) 操作步骤设L =113°25′31.4880″=B 21°58′47.0845″，分别计算其在1954北京坐标系和1980西安坐标系中的高斯平面坐标。
将上述程序以GSJSZ 的文件名输入计算器后，按键 5 2及选择程序GSJSZ ，按键，屏幕提215847.0845542531.48880§2.1.4高斯投影反算与子午线收敛角的计算计算高斯平面坐标为x 、y 的点在参考椭球面上的大地经纬度L 、B 称高斯投影反算，令γ为子午线收敛角。
(1) 数学模型[]642422222200)459061(25.0)935(5.7901nt t n t t n t B B f f f f f f f ff +++-++-+-=ηηπη(2-16) []5423220)24285(5.1)21(30180cos 1n t t n t n B l f f f f f +++++-=ηπ(2-17) []5423220)352(12)1(60180nt t n t n t f f f f f +++-+-=ηπγ (2-18) 0f B 、0B 、0l 、0γ分别是以十进制度为单位的底点纬度、纬度、经度差和子午线收敛角，其余符号的计算公式如下：cyV N yn f f == (2-19)21f f V η+=(2-20) f f B e cos '=η(2-21)f f B t tan =(2-22)底点纬度0f B 的计算公式，对于克拉索夫斯基椭球为：654320)3(90000001939.0)3(70000021572.0 )3(50000485728.0)3(20004353257.0 )3(20057974044.0)3(30246825708.951111537259.27---+-+-----+=X X X X X X B f(2-23)对于IUGG1975椭球为：654320)3(40000001940.0)3(90000021576.0 )3(70000485835.0)3(90004354002.0 )3(60057985065.0)3(90248365772.951116228946.27---+-+-----+=X X X X X X B f (2-24)上式中的610/x X =，也即以千公里为单位的点的x 坐标。
(2)、程序与案例分别在1954北京坐标系和1980西安坐标系中进行高斯投影反算。