线元法路线计算程序(1602) - mr主线
线元法
线元法1：ZHY主程序“1.SZ→XZ”：“2.XY→SZ”：If N=1：Then goto 1 ：Else goto 2：lfEnd Lbl 1：‘S=’?S：“Z=”?Z：Prog〝QXSJK〞: 1÷P→C :（P－R）÷（2HPR）→D : 180÷π→E: Abs（S－O）→W:Prog〝SUB 1〞:〝XS〞: X→X ▲〝YS=〞: Y→Y▲〝FS=〞: F－90→F▲Prog〝GAO〞: goto 1LBI 2 :〝X=〞?X :〝Y=〞?Y :X→I:Y→J：Prog〝QXSJK〞: 1÷P→C :（P－R）÷（2HPR）→D : 180÷π→E: :Prog〝SUB 2〞: 〝S〞: O+W →S▲〝Z=〞: Z▲Prog〝GAO〞:〝SC〞?C〝R=〞:5.89-√((C-G)2+Z2) ▲goto 2“L=”：?L :199.965-(L-49040)*0.006+0.00 →H ◢设计高程R-√（P²+G²）2：SUB1（缓和曲线正算子程序，不能独立运行）0.1739274226→A : 0.3260725774→B : 0.0694318442→K : 0.3300094782→L : 1－L→F : 1－K→M : U+W（A Cos（G+QEKW （C+KWD））+B Cos（G+QELW（C+LWD））+ B Cos（G+QEFW （C+FWD））+A Cos（G+QEMW（C+MWD）））→X :V+W（A Sin （G+QEKW（C+KWD））+B Sin（G+QELW（C+LWD））+ B Sin （G+QEFW（C+FWD））+A Sin（G+QEMW（C+MWD）））→Y : G+QEW （C+WD）+90→F : X+Z cos（F）→X : Y+Z Sin（F）→Y3 : SUN 2（缓和曲线反算子程序，不能独立运行）G－90→T : Abs（（Y－V）Cos（T）－（X－U）Sin（T））→W : 0→Z : LBI 0 : Prog〝SUB 1〞: T+QEW（C+WD）→L :（J－Y）Cos（L）－（I－X）Sin （L）→Z : If Abs（Z）< 1×10－6 : Then Goto 1 : ELse W+Z→W : Goto 0 : IfendLBI 1 : 0→Z : Prog〝SUB 1〞:（J－Y）÷Sin（F）→Z4 : QX SJK （缓和曲线数据库子程序，不能独立运行）If S≥7000 And S<8552.052 : Then 39351.594→U : 72418.7097→V : 7000→O :257.8746719→G : 2000→H : 1×1045→P : 1×1045→R : 0→Q : ELse If S≥8552.052 And S<8752.052 : Then 39025.584→U : 70901.283→V : 8552.052→O: 257.8746719→G : 200→H : 1×1045→P : 1800→R : 1→Q : ELse If S≥8752.052 And S<9900.413 : Then 38987.2071→U : 70705.0275→V : 8752.052→O : 261.0577708→G : 1148.361→H : 1800→P : 1800→R : 1→Q : ELse If S≥9900.413 : Then 39170.3263→U : 69590.9925→V : 9900.413→O : 297.6112367→G : 200→H : 1800→P : 1×1045→R : 1→Q : Ifend : Ifend : Ifend : Ifend : Return注意：0 : 表示零。
线元法线路坐标正反算程序
经苦心钻研，奋战多日，终于编写出了代码短，速度快，精度高，功能全的线路坐标正反算程序，欢迎试用并提出宝贵意见。
功能简介及特点：1、选用高斯-勒让德公式作计算内核，保证精度，模块化设计，便于扩充功能。
2、线元数据可自动从数据库调用，也可手工输入。
3、可管理多条线路，如里程不在线路或线元范围，将警告里程偏大、偏小。
4、边桩计算设计为导线式递推方式，可用于由一个中桩推出结构物所有角点坐标。
5、反算实现了智能化操作，只需输入线路号（或手工输线元资料）、坐标，不需近似里程，即可自动从起点向后开始试算出里程、位置，如对算出里程、位置表示怀疑，还可以让计算器从终点起再向前试算下一个可能的位置（匝道、回头曲线同一坐标可能会有一个以上结果）。
第三次及以后试算才要求输入近似里程。
6、程序代码规范简洁，便于阅读、理解。
完整程序清单:ZFS %正反算主程序B=.1739274226:C=.5-B:Lbl 1:U"0 ZS 1 FS"=0=>Prog "ZS":≠>U=1=>Prog"FS":≠>Goto 1ZS %正算子程序{K}:Prog"ZZ":I=0:{I}:I"L"≠0=>"Prog"WY":≠>Prog"ZB"FS %反算子程序{KVW}:V"XC"W"YC":Lbl 2:Prog "ZZ":I=V-S:J=W-T:Pol(I,J: J=J-F:K=K+Rec(I,J:AbsI<1m=>Prog"WZ":≠>Goto 2ΔM=0:{M}:M"0 NEXT"=0=>U=U+1:Goto 2:≠>U=1ZZ %高斯法中桩子程序(4节点)Prog"XL":M=K-L:O=(P-R)÷2PQR:D=.0694318442:E=.3300094782:F=1:G=1-E:H=1-D:I=5:Lbl 1:C[I]=A+MrC[I](1÷P+OMC[I]:Dsz I:Goto 1:S=X+M(BcosD+CcosE+CcosG+BcosH:T=Y+M(BsinD+CsinE+CsinG+BsinHWY %外移点计算子程序Lbl 1:J=90:{J}:J=F+J"<":F=J:S=S+Rec(I,J:T=T+J: Prog"ZB":I=0:{I}:I"L"≠0=>Goto 1WZ %位置显示子程序"KJ":K:Pause 1:J◢ZB %坐标显示子程序"XY":S:Pause 1:T◢YC %异常处理子程序U=1=>K=L:U=2ΔU=3=>K=M:U=4ΔU=5=>{K}:U=4ΔK<L=>"<<!":Z=1ΔK>M=>">>!":Z=1DL %断链处理子程序"DL":K=L:I>0=>K=L+Q-------------------------以上为程序运算部分,以下为数据库部分-------------------------------XL %线路数据库选择子程序Lbl 1:Z=0:N"0 SD"=0=>Prog"0"△N=1=>Prog"1"△N=2=>Prog"2"△...有几条线路仿上行格式输几行Z=1=>{NLXYOPQRK}:Goto 10 %手工输入子程序L"K0"XYAQ"LS"P"R0"R"RN":M=L+Q:Prog"YC"1 %线路一数据库子程序①Lbl B:L=线路起点里程:M=线路终点里程:Prog"YC":Z=1=>Goto EΔ②Q=线元长:P=起点半径:R=终点半径:K≤L+Q=>X=起点X坐标:Y=起点Y坐标:A=起点方位角:Goto EΔL=L+Q:③......④Q=短链长:K<L+Q=>Prog "DL":Goto BΔL=L+Q:⑤Q=线元长:P=起点半径:R=终点半径:K≤L+Q=>X=起点X坐标:Y=起点Y坐标:A=起点方位角:Goto EΔL=L+Q:⑥......⑦Q=线元长:P=起点半径:R=终点半径:X=起点X坐标:Y=起点Y坐标:A=起点方位角:Lbl E2 %线路二数据库子程序输入要求和线路一相同。
线元法万能坐标计算程序
线元法万能坐标计算程序(适用于CASIO fx-9750GⅡ计算器)论文/：本论文仅供学习交流使用，本站仅作合理转载，原作者可来邮要求删除论文。
摘要：我国公路建设事业正处于一个高速发展的时期，在公路工程施工过程中，施工技术人员经常要使用全站仪、水准仪进行施工放样、高程测量，在测量过程中，手工计算速度慢，失误率高，工作效率极低。
利用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Ⅱ编程函数计算器即可。
5800公路测量程序22版(7条路线)
FX5800计算器公路测量常用程序集2.2 版一、程序功能主要功能：采用线元法与交点法相结合计算多条线路坐标正反算，可算任意复杂线型及立交匝道，包括C型，S型、卵型、回头曲线等；极坐标放样，全线路基边坡开挖口及坡脚放样计算、路基任意点设计高程、横坡、设计半幅宽度.隧道欠超挖放样计算。
新版本优化：1、调整一次显示结果；2、交点法中考虑了不对称缓和曲线；3、修改原版本部分地方笔误.(红色为修改处)二、源程序1.总主程序（1、坐标计算放样程序2、坐标反算程序；3、高程计算查阅程序；4、路基半幅标准宽度查阅程序；5、路基边坡及开挖口放样程序；6、路基标准距离放样；7、桥梁锥坡计算放样程序；8、极坐标计算程序；9、隧道超欠挖计算程序）运行后输入1~9数子则选择1至9的程序，返回时，在桩号输入-1,返回选择选择计算类型。
输入-2,返回选择线路。
坐标计算中输入-3,则显示本段曲线要素。
程序名：0ZCXLbl 0: 24→DimZ：Norm 2:1→A：＂A:XY=1，ZD=2 ，GC=3，GD=4， BP=5，FM=6，ZP=7，JS=8，SD=9＂?A:A=1=>Goto 1:A=2=>Goto 2:A=3=>Goto 3:A=4=>Goto 4:A=5=>Goto 5:A=6=>Goto 6:A=7=>Goto 7:A=8=>Goto 8:A=9=>Goto 9Lbl 1:Prog＂DX＂:LbI A:Prog＂QX＂:90→B:＂PJ1＂?B:B→C:＂PJ2＂?C:B→Z[1]:C→Z[8]:Lbl B:1→F:＂KM＂?Z:Z= -1=>Goto 0:Z=-2=>Goto A:Z=-3=>GotoX:Prog＂KM＂:?D:Porg＂THB＂:0→L:＂L0＂?L:Z[2]+Z[1]-Z[8] →E:X+L cos(E)→X:Y+Lsin(E) →Y:Prog＂XY＂:Prog＂JS＂:Goto BLbl 2:2→F:90→Z[1]:Prog＂QX＂: Lbl C:＂DKLC＂?Z:Z=-1=>Goto 0:Z=-2=>Goto 2:Z=-3=>Goto X:Prog＂KM＂: ＂X0＂?X: ＂Y0＂?Y:Prog＂THB＂:Porg＂ZD＂:Goto CLbI 3:Prog＂QX＂:→B:＂H-B＂?B:B→Z[9]:＂SXG＂?L:Lbl D:＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto 3:?D:＂ZS＂?P: Porg＂H＂:Fix 3:＂SJG=＂:Lcoate 6,4,H-Z[9]:＂SCG=＂:Lcoate 6,4,F:＂I=＂:Locate 6,4,I×100: If U<0:Then ＂W=＂: Locate 6,4, U: Else＂T=＂: Locate 6,4, U : IfEnd◢Goto DLbl 4:Prog＂QX＂:LbI E:＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto 4:?D:Prog＂GD＂:Fix 3:＂ SJGD=＂:Locate7,4,L: Goto ELbl 5:Prog＂QX＂:0.5→B:＂TH-GD＂?B:B→Z[23]:Lbl F:2→F:90→Z[1]:＂KM＂?Z:Z=-1 =>Goto 0:Z=-2 =>Goto 5:Prog＂KM＂:＂X0＂?X:＂Y0＂?Y:0→M:＂M0＂?M: M→Z[4]:Prog＂FBX＂:Fix 2: P→D ：Abs（D）-S→O: ＂ L0=＂:Locate 6,4,O:Prog＂ZD＂: Z[4]→T:＂ TW=＂:Lcoate 6,4,T◢ Goto FLbl 6:Prog＂DX＂:LbI G:Prog＂QX＂:LbI H:1→F:90→Z[1]:＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto G:Prog＂KM＂:?D:Prog＂THB＂:Prog＂XY＂:Prog＂JS＂:Prog＂H＂:0→M: ＂M0＂?M:Fix 2:H-M→T:＂ TW=＂:Locate 6,4,T◢ Goto HLbl 7:Prog＂DX＂:LbI I:Prog＂QX＂: ＂Z0＂?Z:0→D: ＂LD:Z-,Y+＂? D:Abs(D) →R:＂LR＂?R: D→Z[6]:R→Z[7]:Z→Z[8]:LbI J:0→Y: ＂L0:DZH+,XZH-＂?Y:Y=-1=>Goto 0:Y=-2=>Goto I:Prog＂ZP＂:Goto JLbl 8:Prog＂DX＂:Lbl K:?X:X=-1=>Goto 0:X=-2=>Goto 8:?Y:Prog＂JS＂:Goto KLbl 9:Porg＂DX＂:Lbl L:Prog”QX”:Lbl M:2→F:90→Z[1]: ＂KM＂?Z:Z=-1=>Goto 0:Z=-2=>Goto L:Z=-3=>Goto X: Prog＂SD1＂: 0.5→A:＂CQHD＂?A:“X0＂?X:＂Y0＂?Y:0→M:＂M0＂?M:M→Z[8]:A→Z[15]:Prog＂SD＂: Goto M: Lbl X:Mat F◢ Goto 02、次程序:路基开挖边线及填方坡脚线放样程序（输入填方放宽值、大概桩号及测量坐标、地面标高计算出偏移距离、桩号、距中距离、填挖高度）程序名：FBXProg“THB”:D→P:Prog“GD”:If D<0:Then 0.75-L→D:Goto H:Else L-0.75→D:Goto H:IfEnd:Lbl H:Prog“H”:H-0.03-Z[4] →Z[4]: Prog “W-1”:If Z[4]<0:Then –Z[4]→G:Goto W:Else Z[4] →G:Goto T:IfEnd:Lbl W: Z[10]+Z[11] +Z[12] +Z[13]→A: Z[10]+Z[11]+Z[12]→B:Z[10]+Z[11]→C: If G>A:Then Goto 1:Else If G>B:Then Goto 2:Else If G>C:Then Goto 3:Else If G>Z[10]:Then Goto 4:Else Goto 5:IfEnd:Lbl 1: L+Z[14]+Z[15]+Z[16]+Z[17]+Z[18]+(G-A)×Z[9]+Z[13]×Z[8]+Z[12]×Z[7]+Z[11]×Z[6]+Z[10]×Z[5]→S:Goto Z:Lbl 2: L+Z[14]+Z[15]+Z[16]+Z[17]+(G-B)×Z[8]+Z[12]×Z[7]+Z[11]×Z[6]+Z[10]×Z[5]→S:Goto Z:Lbl 3:L+Z[14]+Z[15]+Z[16]+(G-C)×Z[7]+Z[11]×Z[6]+Z[10]×Z[5]→S:Goto Z: Lbl 4:L+Z[14]+Z[15]+(G-Z[10])×Z[6]+Z[10]×Z[5]→S:Goto Z:Lbl 5:L+Z[14]+G×Z[5] →S:Goto Z:Lbl T:L+Z[23]→T:If G>Z[21]:Then T+Z[22] +Z[21]×Z[19]+(G-Z[21])×Z[20]→S:Goto Z:Else T+G×Z[19]→S:Goto Z: IfEnd:Lbl Z3.次程序3.隧道超欠挖值计算放样程序（输入隧道线路，大概桩号、输入衬砌厚度、测量三维坐标，计算准确桩号及位置、计算欠超挖值）程序名：SDProg “KM”:Prog“THB”: D→Z[13]:0→D:Prog“H”:Cos(Z[10])×Z[4] →E: H+Z[5] →Z[5]: H+Z[7] →Z[7]: E+Z[5] →Z[10]:If Z[8]>Z[10]:Then Goto R:Else Goto S:IfEnd:Lbl R:√((Z[13]-Z[9]) 2+(Z[8]-Z[5]) 2 )- Z[4]-Z[15] →L:Goto L:Lbl S:Z[4]-Z[6] →S:√(S 2-(Z[7]-Z[5]) 2 )×S÷Abs(S)→T:Abs(Z[13]-Z[9])-T→T:√(T 2+(Z[8]-Z[7]) 2 )- Z[6]-Z[15] →L:Goto L:Lbl L: Fix 2: “L0=”: Locate 6,4,L:Prog“ZD”:“H0=”: Locate 6,4,Z[8]-H◢4：极坐放样计算程序(计算放样点至置仪点方位角及距离)程序名：JS(一式)X：Y：Z[11]→K:Z[12]→L:Pol(X-K, Y-L):IF J<0:Then J+360→J:IFEnd:Int(J)+0.01Int(60Frac(J))+0.006Frac(60Frac(J)) →J: Fix 4: “FWJ=”: Locate 6,4,J:Fix 3:“JL=”:Locate 6,4,I◢JS(二式)X：Y：Z[11]→K:Z[12]→L:Pol(X-K, Y-L):IF J<0:Then J+360→J:IFEndFix 3:“JL=”:Locate 6,4,I: J▲DMS◢5．路基宽度子程序：路基标准半幅宽度计算程序(对于设计有加宽渐变的有用，如路基宽度无变化，则把此程序直接输入半幅宽度值至L)程序名： GD1→S: Prog “G-1”:Z-C→E: E ×(B-A)/S+A→L:6．坐标计算次程序程序名：THBIf Q>1:Then Goto J: Else 1÷P→C:(P-R)÷(2HPR)→S:180÷π→E:IfF=1:Then Abs(Z-O) →W:Prog "A":Ｇoto 2:Else X→I:Y→J:Prog ＂B＂:O＋Ｗ→Z:D→D:Ｇoto 2: IfEnd:LbI J: If F=1:Then Prog "Z":Ｇoto 1:Else Prog ＂ZX＂:Ｇoto 2: IfEnd: LbI 1:I+D×COS(Z[2]+Z[1]) →X: J+D×Sin(Z[2]+Z[1]) →Y: LbI 27. 线元法正算子程序(Ａ)程序名：A0.1184634425→A: 0.2393143352→B：0.2844444444→N：0.0→K：0.2307653449→L：0.5→M：U+W(Acos(G+QEKW(C+KWS))+Bcos(G+QELW(C+LWS))+Ncos(G+QEMW(C+MWS))+Bcos(G+Q E(1-L)W(C+(1-L)WS))+Acos(G+QE(1-K)W(C+(1-K)WS)))→X:V+W(Asin(G+QEKW(C+KWS))+Bsin(G+QELW(C+LWS))+Nsin(G+QEMW(C+MWS))+Bsin(G+Q E(1-L)W(C+(1-L)WS))+Asin(G+QE(1-K)W(C+(1-K)WS))) →Y:G+QEW(C+WS)→Z[2]:X+Ｄcos（Z[2]+Z[1]）→Ｘ: Y+Ｄsin（Z[2]+Z[1]）→Ｙ8. 线元法反算子程序(Ｂ)程序名： B G-90→T：Abs((Y-V)cos(T)-(X-U)sin(T))→W：0→D：Lbl 0：Prog "A"：T+QEW(C+WS)→L：(J-Y)cos(L)-(I-X)sin(L)→D：IF Abs(D)<0.01:Then Goto 1：Else W+D→W:Goto 0: IfEnd: Lbl 1:0→D：Prog "A"：(J-Y)÷sin(Z[2]+90) →D9. 交点法正算子程序(Z)程序名：ZH2÷R÷24-H∧(4)÷2688÷R∧(3)→A:H÷2-H∧(3)÷240÷R2→B:（（H2-N2）÷24÷R）÷Sin(Abs(P))-((H∧(4)-N∧(4))/2688/R∧(3))÷Sin(Abs(P))→E:(R+A)tan(Abs(P)÷2)+B-E→T:P÷Abs(P)→W:0→M:H→C:If Z≤O -T:Then Z-O→S:G→Z[2]:Goto 2: IfEnd:If Z≤O -T+H:Then Z-O +T→S:Prog “HX”:G+WK→Z[2]:Goto 4:IfEnd:If Z≤O -T+πR×Abs(P)÷180+H÷2-N÷2: Then 180(Z-O +T-0.5H)÷R÷π→S : A+R(1-Cos(S))→B : H÷2-H∧(3)÷240÷R2+Rsin(S)→A:R→M:G+WS→Z[2]:Goto 4: IfEnd:O -T+πR×Abs(P)÷180+H÷2+N÷2-Z→S:(R+N2÷R÷24-N∧(4)÷2688÷R∧(3))tan(Abs(P)÷2)+N÷2-N∧(3)÷240÷R2+E→T :N→H:Prog “HX”:G+P→S:S-WK→Z[2]:U+(T-A)Cos(S)-WBSin(S)→I:V+(T-A)Sin(S)+WBcos(S)→J:Goto 3:Lbl 4:U+(A-T)cos(G)-WBsin(G)→I:V+(A-T)Sin(G)+WBcos(G)→J: Goto 3:Lbl 2:U+Scos(Z[2])→I:V+Ssin(Z[2])→J: LbI 3:C→H10. 交点法缓和段转化子程序(HX)程序名：HXS-S∧(5)÷40÷R2÷H2+S∧(9)÷3456÷R∧(4)÷H∧(4)→A:S∧(3)÷6÷R÷H-S∧(7)÷336÷R∧(3)÷H∧(3)+S∧(11)÷42240÷R∧(5)÷H∧(5)→B:90S2÷π÷R÷H→K:RH÷S→M11. 交点法反算子程序(ZX)程序名：ZXZ:0→D:Lbl 0:Prog“Z”:Pol(X-I,Y-J):J-Z[2]→J:Isin(J)→S:Icos(J) →I:If Abs(I)<0.1:Then Z+I→Z:S→D:Goto 2:Else Goto 1: LbI 1:If M=0:Then Z+I→Z:Goto 0:Eles Pol(M-WS,I):(JMπ)/180→I:Z+I→Z:Goto 0:IfEnd:LbI 212．高程计算子程序（H）程序名：HProg “S-1”:N-I→A:A÷Abs(A→W:R×Abs(A÷2→T:If Z<C-T:Then G+I(Z-C→H:Goto 0:IfEnd:If Z<C+T:Then G-TI+(I+(Z-C+T)W÷2÷R)×(Z-C+T→H:Goto 0: IfEnd:N→I: G+I(Z-C→H:Lbl 0:If D=0:Then H→H: I→I:Goto I:Else Prog “I”:H+V→H:Abs(L-13.高程超高计算程序（I）程序名：IIf Z[3]=1:Then Prog“I-1”:Goto 1:Else If Z[3]=2:Then Prog“I-2” :Goto1: Else If Z[3]=3:Then Prog“I-3” :Goto 1: Else If Z[3]=4:Then Prog“I-4” :Goto 1: Else If Z[3]=5:Then Prog“I-5” :Goto 1: Else IfZ[3]=6:Then Prog“I-6” :Goto 1: Else Prog“I-7”:Goto 1: IfEnd:Lbl 1: If W=1:Then Goto Z:Else Goto X: IfEnd :Lbl Z:If S=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×((N-M)×(Z-C)÷S+M)→V:Goto 2:IfEnd:Lbl X:If S=0:Then Abs(D)×M→V:Goto 2:Else Abs(D)×(((3((Z-C)÷S)2-2((Z-C)÷S)∧(3))×(N-M))+M)→V:Goto 2:IfEnd:Lbl 2:Abs(D)→E:V÷E→I:I(E-K)→V14．导线点子程序（DX）程序名：DXZ[11]→K:Z[12]→L:“XZ”?K:“YZ”?L:K→Z[11]:L→Z[12]15．线路选择子程序(线路选择输0时,则输曲线参数)（QX）程序名：QX1→A:“LX:XX-1,XA-2,XB-3,XC-4,XD-5,XE-6,XF-7”?A:A→Z[3]:If A=0:Then Goto 1:Else Goto 2: IfEnd:Lbl 1:“JDO”?O:“XO”?U:“YO”?V:“FW”?G:“LS”?H:?P:?R:“J=2,3;X=(0,1,-1)”?Q:If Q=2:Then H→N:Goto 2:Else If Q=3:Then“LS”?N:Lbl 216．线路中线元段判别子程序（KM）程序名：KMIf Z[3]=0:Then Goto 0:Else If Z[3]=1:Then Prog “XX-1”:Goto 0: Else If Z[3]=2:Then Prog “XA-2”:Goto 0: Else If Z[3]=3:Then Prog “XB-3”:Goto 0: Else If Z[3]=4:Then Prog “XC-4”:Goto 0: Else IfZ[3]=5:Then Prog “XD-5”:Goto 0: Else If Z[3]=6:Then Prog “XE-6”:Goto 0: Else Prog “XF-7”:Goto 0: IfEnd: LbI 017．坐标计算参数转化子程序（F）程序名：FMat F[1,1]→O:Mat F[1,2]→U:Mat F[1,3]→V:Mat F[1,4]→G:Mat F[1,5]→H:Mat F[1,6]→P:Mat F[1,7]→R: Mat F[1,8]→Q:Q=2=>H→N:Q=3=>Mat F[1,9]→N18.桥梁锥坡放样子程序（ZP）程序名：ZP√(Z[6]2（1-Y2/Z[7]2）) →X:Z[8] →Z:Prog“GD”:If D>0:Then L+X→D:Else –L-X→D:IfEnd:Z[8]+Y→Z:D:1→F:90→Z[1]:Prog“KM”:Prog“THB”:Prog“XY”:Prog“JS”19．显示子程序（XY）和（ZD）①程序名：XYFix 3:＂XS=＂:Locate 6,4,Ｘ：＂YS=＂:Locate 6,4,Ｙ：②程序名：ZDFix 3:＂KM=＂:Locate 6,4,Z：＂ D=＂:Locate 6,4,Z：20．数据子程序(附后示例)①-1程序名：X X-1（线路1坐标计算要素程序）If Z≤26615.555:Then [[25900,11587.421,1847.983, 101。
线元法路线计算程序
线元法路线计算程序线元法（LE法）是一种用于计算电力系统潮流分布的方法，它将电力系统抽象成节点和支路的网络，通过对节点和支路进行编号，可以建立节点电压和潮流分布之间的方程，进而求解电力系统中各节点的电压和潮流分布。
下面是一个用于计算线元法路线的程序。
1.定义节点和支路：首先，我们需要对电力系统的节点和支路进行定义。
节点可以是发电站、变电站或负荷节点；支路可以是输电线路或变压器。
每个节点和支路都需要有一个唯一的编号，以便在后续的计算中进行引用。
2.建立节点电压方程：根据电力系统的KCL（电流平衡方程），我们可以得到节点电压方程。
每个节点的电压方程可以表示为：V(i) = Σ{(V(j) - V(i)) / Z(ij)}，其中V(i)表示第i个节点的电压，V(j)表示第j个节点的电压，Z(ij)表示第i个节点到第j个节点的支路阻抗。
3.建立支路潮流方程：根据每个支路的电流平衡方程，我们可以得到支路潮流方程。
每个支路的潮流方程可以表示为：I(ij) = (V(i) - V(j)) / Z(ij)，其中I(ij)表示从第i个节点到第j个节点的支路电流，V(i)表示第i个节点的电压，V(j)表示第j个节点的电压，Z(ij)表示第i个节点到第j个节点的支路阻抗。
4.解线性方程组:将节点电压方程和支路潮流方程组合成一个线性方程组，我们可以通过求解线性方程组，得到电力系统中各节点的电压和潮流分布。
5.输出结果:根据求解的节点电压和支路潮流，我们可以将结果输出，以便进行分析和评估。
下面是一个基本的线元法路线计算程序的伪代码示例：```Input: 节点和支路的定义，节点电压和支路阻抗的初始值Output: 节点电压和支路潮流的计算结果1.建立节点电压方程和支路潮流方程-初始化节点电压和支路潮流的初始值-根据节点电压方程和支路潮流方程2.解线性方程组-使用数值计算方法求解线性方程组，得到节点电压和支路潮流的计算结果3.输出结果-将节点电压和支路潮流的计算结果输出，以便进行分析和评估```这是一个简化的线元法路线计算程序的框架，具体实现时需要根据具体的电力系统结构和算法细节进行调整和优化。
4800、4850线元法单心圆隧道三维(坐标正反计算、高程)程序
2010-4-2修改版1. XH（此为主程序）Lbl 0：”1.LC=>XY”：“2.XY=>LC”：“3.ZHZL=>GC”：{V}：V=1=>GOTO 1△V=2=> GOTO 2 △V=3=>GOTO 3△:≠> GOTO 0△Lbl 3：{HZ}：H”ZH”：Z”XLZJ,-Z+Y”：Z”SDZF，-Z+Y”=Z-0.125▲Prog”SQXBG”：GOTO 0：Lbl 1 ：{L}：L”ZH”：L>173000=>L<174661.96=>GOTO 4△△GOTO 0△Lbl 4：{QJ}：L”ZH”=L：Q”XLZJ,-Z+Y”：Prog “QXYS”：Q=0=>J=0：≠>J=90△Prog “ZSZB”：Q “SDZJ,-Z+Y”=Q-0.125▲ X”X=“▲ Y”Y=“▲O”FWJ=“▲H=L：Z=Q：Prog “SQXBG”：GOTO 0：Lbl 2：{MR}：M”XO”：R”YO”：M>3845505.273 =>M<3846506.099=>R>499371.832 =>R<500352.224 => GOTO 5△△△△GOTO 0△Lbl 5：L：Q=0：J=0：M”XO”=M：R”YO”=R：L”ZH” =173300：Prog “FSZB”：L”LC=“▲Q”JL=“▲Q “SDZJ，-Z+Y”=Q-0.125▲H=L：Z=Q：Prog “SQXBG”：GOTO 02.正算坐标ZSZBH=( L-S ) / 4：F=90/π：U=HHF（1/T-1/I）/（K-S）：D=2HF/ I：O=C+4D+16U：P=O+J ： E=C+ D+ U： W=C+2D+4U： G=C+3D+9U：X=A+AbsH/3*(cosC+4(cosG+cosE)+2cosW+cosO)+Qcos P ：Y=B+AbsH/3*(sinC+4(sinG+sinE)+2sinW+sinO)+Qsin P3.反算坐标：FSZBLbl 0：Prog “QXYS”：Prog “ZSZB”：Z=O-90：P= (R-Y)cosZ-(M-X) sinZ ：L=L+P：AbsP≥0.001=> GOTO 0 ：≠> GOTO 1 △Lbl 1：Q= (R-Y)cosO-(M-X) sinO4. 曲线元要素数据库：QXYSL≥S=>L＜K=> S=**：A=**：B=**：C=**：I=**：K=**T=**⊿⊿←┘L≥S=>L＜K=> S=**：A=**：B=**：C=**：I=**：K=**T=**⊿⊿←┘L≥S=>L＜K=> S=**：A=**：B=**：C=**：I=**：K=**T=**⊿⊿←┘L≥S=>L＜K=> S=**：A=**：B=**：C=**：I=**：K=**T=**⊿⊿←┘L≥S=>L＜K=> S=**：A=**：B=**：C=**：I=**：K=**T=**⊿⊿←┘……………………………L≥S=>L＜K=> S=**：A=**：B=**：C=**：I=**：K=**T=**⊿⊿←┘(注：如有多个曲线元要素继续添加入数据库QXYS中)5高程计算主程序：SQXBG(后有修改说明){X}:X”SCGC”：H”ZH”：Z”ZJ,-Z+Y”：M”R”：N”CS”：Prog＂SJK＂：F=C-D：T=Abs(RF÷2)：R=R AbsF÷F：H≤B-T=＞K=0：≠=＞H≥B+T=＞K=0：C=D：≠=＞K=H-B+T⊿⊿G”XLZG”=A-（B-H）C-K^2÷2R▲Z≥-0.125=＞Z=Z-0.125：X>100=＞J”YO1XGC”=G+N▲J”Y GCFSKD”= Abs√(M^2-(X-（G+N）)^2) ▲J”Y KDCQ,+C,-Q”=J- Abs (Z+0.125) ▲J”Y SBSJGC”=G+N+√(M^2-(Z+0.125)^2) ▲J”Y GCCQ,C+,Q-”=X-J ▲Goto1：≠=＞J”Y SBSJGC”=G+N+√(M^2-(Z+0.125)^2) ▲Goto1：≠=＞Z=Z-0.125：X>100=＞J”Z O1XGC”=G+N▲J”ZGCFSKD”= Abs√(M^2-(X-（G+N）)^2) ▲J”ZKDCQ,+C,-Q”=J-Abs(Z+0.125) ▲J”ZSBSJGC”=G+N+√(M^2-(Z+0.125)^2) ▲J”ZGCCQ,C+,Q-”=X-J ▲Goto1：≠=＞J”ZSBSJGC”=G+N+√(M^2-(Z+0.125)^2) ▲Goto1：Lbi16高程计算主程序子程序：SJKH≤第二竖曲线起点桩号=>A=第一竖曲线交点高程：B=第一竖曲线交点桩号：C=第一竖曲线前坡度：D=第一竖曲线后坡度：R=第一竖曲线半径：Goto1⊿H≤第三竖曲线起点桩号=>A=第二竖曲线交点高程：B=第二竖曲线交点桩号：C=第二竖曲线前坡度：D=第二竖曲线后坡度：R=第二竖曲线半径：Goto1⊿H≤第四竖曲线起点桩号=>A=第三竖曲线交点高程：B=第三竖曲线交点桩号：C=第三竖曲线前坡度：D=第三竖曲线后坡度：R=第三竖曲线半径：Goto1⊿………………………继续添加要素：Goto1：Lbi1说明：第一部分坐标部分（1、2、3、4）V=1进入坐标正算V=2进入坐标反算V=3进入单独的高程计算当V不等于1、2、3时，则返回程序，要求再次输入V值。
公路、铁路曲线坐标计算程序
公路、铁路曲线坐标计算程序一、主程序：M（M是文件名，下同）L1 A”CX”B”CY”C”HX”D”HY”E”ZHX”F”ZHY”Q”HZX”L”HZY”G”ZHDK”H”ZHJD°”V”A°”W”L0”K”L1”RL2 Pol(C-A,D-B):M=JL3 Lbl 2:J=90:{NOJZ}:Z”DK”N”Z/!/Y(1/2/3)”:N=2=>Goto 1⊿O”M”J”A°”L4 Lbl 1:P=Z-G:S=E:T=F:U=H:I=H:FixmL5 P≤0=>Prog ”3”:≠=>P≤W=>Prog ”1”:≠=>P≤W+K=>Prog ”2”:I=90W/π/R+180π-1R-1(P-W)⊿⊿P≤W+K=>Prog ”4”:I=H+I⊿⊿L6 P>W+K=>P=2W+K-P:S=Q:T=L:P>0=>Prog ”1”:X=-X:U=H+V:Prog ”4”: I=U-I:≠=>U=H+V:I=U:P=-P:Prog ”3”⊿⊿L7 N≠2=>P=O:U=I-J:N=3=>P=-P⊿S=X:T=Y:Prog ”3”⊿L8 Pol(X-A,Y-B):J=J-M:J<0=>J=J+360⊿L9 J:”°°°=”◢I:”S=”◢X:”X=”◢Y:”Y=”◢Goto 2二、子程序：1L1 I=90P2(πRW)-1L2 X=P-Px y5(40R2W2)-1L3 Y=PIπ/540三、子程序：2L1 Y=180π-1R-1(P-0.5W)L2 X=RsinY+W/2- Wx y3/240/R2L3 Y=W/24/R+R-RcosY四、子程序：3L1 X=S+PcosUL2 Y=T+PsinU五、子程序：4L1 V<0=>Y=-Y:I=-I⊿L2 S=S+XcosU-YsinUL3 Y=T+XsinU+YcosUL4 X=S说明：1、该程序适用于计算器 CASIO fx-4800P。
fx9750卡西欧9750公路线元法计算程序
新版程序把线元法和交点法已经集成在一个模块中了，用户只需修改JD程序和ZA程序中的数据部分即可,其余不需作任何的改动。
2.因为每条路高程计算不尽相同，且比较复杂，现在可利用PC机EXCEL计算好打印成表格带到工地上使用，所以本版程序未对线路高程序进行专门的编程计算，而是利用统计计算模式中来输入桩号（第一列X）及左、右高程（第二、三列Y，Freq），这种输入数据的方式最为直观，易发现错误，也易修改，输入完毕后运行S 程序对数据按桩号进行排序，在程序中通过调用GG程序来进行内插计算，SG=-1得左标高，SG=1得右标高（若SG输入0,则可进行一般的线性内插计算）。
3.在JD程序和XY程序中，先将一个计算单元的数据置入矩阵F中（1行8列或1行9列），这样程序可读性极好。
4．相比原CASIO4850程序操作习惯，作了一点小小的改动，测站坐标存在Z[10]，N中，X坐标原存在M中容易被误操作修改，而设计标高存在M中，这样易于修改，因为CASIO5800没有IN，OUT功能，很不方便。
4.程序利用Z[2]变量值来判断是采用交点法还是线元法模型计算，Z[2]＝0为线元法，否则为交点法。
一、PQX程序：计算中边桩坐标及近似的桩号反算，在运行模式直接调用。
① Z*10+→S:”XO”?S:S→Z*10+:”YO”?N:Prog “AU”② Lbl 2:?L:Prog “Z”:Prog “E”:1n→O:90→S③ Lbl 4:”JJ”?S:”YC”?O:SO=0 =>Goto 2‘原来lbl 后没有标号４的。
④ O=-1 =>Goto 6⑤ “X,Y”:R+OCos(Z+S)→X▲U+OSin(Z+S)→Y▲Prog“D”:Goto 4⑥ Lbl 6:Z*7+→X:Z*8+→Y:”XF”?X:”YF”?Y:XY=0 =>Goto 4⑦ X→Z*7+:Y→Z*8+:Pol(X-R,Y-U+1p):Z+S-J→J:”YC,DL,L”:ICos(J)→O▲ISin(J)→I▲L+I▲Goto 6二、P程序：在程序中提供一个自由运算的模式。