FX-4850计算器曲线计算程序2011.721

合集下载

Casio 4850计算器数据库型全线高程计算程序

CASIO4850计算器公路测量计算程序程序设计:魏加训2009.3.10Casio 4850计算器数据库型全线任意位置高程计算程序一、主程序：HLbl 0: X “BM-H”：U “HS”：“SXG=”: U+X◢输入水准点高程及后视，显示视线高程Lbl 1：{ P K }: “QS”K “K××+×××”: Prog“SJK”↙（换路修改此处为对应数据库名即可）C= C÷100: D= D÷100:↙D>C=> W=1:≠>W = -1⊿K<A=> I= -C:≠> I=D⊿↙Abs (A-K):L> R Abs(D-C)÷2 =>H=Z+IL: "H(Z-SJ)=": H◢Goto2↙≠>H= Z+IL+W(R Abs(D-C)÷2- Abs (A-K))2÷2÷R: "H(Z-SJ)=": H◢Lbl 2: {ONJP}:O “JGC(h)” : N“Z”:J “HP”:N=Abs N : B=H-O+N J÷100: "H(C) =" : U+X-P :Pause 0 : “H(S)=”: B◢“⊿h=”: B-U-X+P◢输入前视后显示实测高、设计高及高差（+填，-挖）Goto 1二、竖曲线数据库格式：K<下一竖曲线起点里程=>C=本竖曲线变坡点之前坡度(%前分子有正负): D=本竖曲线变坡点之后坡度(%前分子有正负): R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙三、数据库样例: SJKK<下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负):D=本竖曲线变坡点之后坡度(%前分子有正负):R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙K<下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负):D=本竖曲线变坡点之后坡度(%前分子有正负):R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙…………………………………………………………………………………………………………………..(注：如一条路有多个竖曲线要素则按顺序排列依次添加入数据库S.Z中)K≤下一竖曲线起点里程=> C=本竖曲线变坡点之前坡度(%前分子有正负): D=本竖曲线变坡点之后坡度(%前分子有正负): R=本竖曲线半径: A=本竖曲线变坡点桩号: Z=本竖曲线变坡点高程: Goto 1⊿↙Lbl 1↙四、输入与显示说明BM-H ？输入水准点的高程HS？输入水准点上的后视读数SXG= ×××显示视线高程“QS”？输入所求点前视读数K××+×××? 输入所求的桩号H(Z-SJ)= ×××显示所求桩号的路面中桩设计高程JGC(h) ? 输入结构层厚度或从设计路面高程点处垂直向下的厚度Z ？输入所求点距中线的边距(注:不分正负号均可，有中央分隔带的输入时需自行扣除中央分隔带半宽)HP？输入所求桩号处的横坡(%前分子有正负).H(C) = 显示计算得出的所求点的实测高程H(S) = 显示计算得出的所求点的设计高程⊿h= 显示设所求点的计与实测的高差（+填，-挖）注：（①当Z和JGC(h)为0时，则H=×××为路面中桩设计高程等于前面显示的H(S) ；②当Z为0，JGC(h)不为0时，则H=×××为路面垂直向下至某层的中桩设计高程；③当Z、JGC(h)、HP三者同时不为0时，则H=×××为路面垂直向下至某层的边桩对应Z位置的设计高程；④当Z不为0、JGC(h)为0、HP不为0时，则H=×××为路面边桩对应Z位置的设计高程；）附：数据库中用到的竖曲线基本计算公式：竖曲线起点桩号=本竖曲线变坡点桩号+本竖曲线之切线长T竖曲线切线长T=RAbs（i2-i1）÷2。

4850全线连线计算程序

CASIO4850全线坐标正反算程序主程序：ZHANG LIANGLbI Ø ：Deg:T"WJM":C=2Ø+4Ø(T-1):T=Ø:N"1.SR,2.JS":N=2=> Goto A⊿M"JZ-XX":O"QD-ZH":U"X-Ø ":V"Y-Ø ":G"FWJ ":Lbl B:{HRQP}: P"QD-R"：H"XX-CD":R"ZD-R":Q").-1,(.+1,--.Ø":W=H:Z =Ø:Prog"S1":Prog"S4":U=X:V=Y:O=O+H:G=F-9Ø:P=R:"…SR…NEXT…!”:Isz T: Dsz M:Goto B⊿Lbl A：Fixm:{XYSZN}:T=Ø:N"1.JS-XY,2.XY-FS": N=2=>Goto 2⊿S"JS-ZH=":Z"D=":Lbl 1:S＞Z［8T+1+C］=>Isz T:Goto 1⊿Prog "S3":W=S-O: Prog "S1": "JS-X=":X：Pause Ø:"JS-Y=":Y◢Goto AΔLbl 2: X"X…":Y"Y…":I=X:J=Y:Lbl 3: Prog "S3":W=Abs((Y-V)cos(G-9Ø)-(X-U)sin(G-9Ø:W＞H=> Isz T:Goto 3⊿Prog "S2":"FS-ZH…":S=O+W:Pause Ø:"D…":Z ◢Goto A⊿子程序S1A=.1184634425:B=.2393143352:N=.2844444444:K=. Ø4691ØØ77:L=.23Ø7653449:D=Q r W2(P-R)÷2HPR:E=Q r WP-1:Z[1]=G+KE+K2D:Z[2]=G+LE+L2D:Z[3]=G+.5E+.25D:Z[4]=G +(1-L)2D+E(1-L:Z[5]=G +(1-K)2D+E(1-K:F=G+E+D+9Ø：X=U+ZcosF+W(AcosZ[1]+BcosZ[2]+NcosZ[3]+BcosZ[4]+AcosZ[5:Y=V+ZsinF+W(AsinZ[1]+BsinZ[2]+NsinZ[3]+Bsin Z[4]+Asin Z[5：子程序S2W=Abs((Y-V)cos(G-9Ø)-(X-U)sin(G-9Ø:Z=Ø:Lbl 4:Prog "S1":L=G-9Ø+E+D:Z=(J-Y)cosL-(I-X)sinL:AbsZ<.ØØ1=>Goto5:≠>W=W+Z:Goto 4⊿Lbl 5:Z=Ø:Prog "S1":Z=(J-Y)÷sinF:子程序S3H=Z［8T+C+6:O=Z［8T+C+1］-H:U=Z［8T+C+2:V=Z［8T+C+3:G=Z［8T+C+4:P=Z［8T+C+5:R=Z［8T+C+7: Q=Z［8T+C+8:子程序S4Z［8T+C+1］=O+H: Z［8T+C+2］=U: Z［8T+C+3］=V:Z［8T+C+4］=G: Z［8T+C+5］=P:Z［8T+C+6］=H:Z［8T+C+7］=R: Z［8T+C+8］=Q: 结束程序编制说明：1.本程序是为CASIO-fx4850p编制的。

CASIO fx—4850P计算器高程控制网按间接平差计算程序

则法方程式为ＣＸ＋Ｌ＝０Ｌ
解上式得ｘ＝一ｃＬＬ（）３
Ｍｅ：ｅ：Ｎ：ｉ：＝Ｔ＋Ｎ：＝０５Ｎ＋１）ｌＤｇＴＭＦｘＡｍＷ．Ｎ（：
式中：～观测值改正数向量，Ａ一误差方程
＝
单位权中误差：＝￣（ＰＶＵ／［Ｖ］÷（ｎ—ｎ）待ｒ）
Ｌｌ：［Ｍ＋］＝＝ｚＩ＝ｚＰ一［＋Ｍ＋ｂＦｚＡ＋ＪＰ［］［］ｚＡ２
ＪｌＬｌＧ： ≠ＭｊＧｔＪ＝０ＩｂＪｏｏＤ： ≠ＡｊＧｔＣｏｏ
将程序中划线部分内容删除，节省程序容量２５字可３
其法方程式为：
ＡＰＸ＋ＡＰＡＬ＝０（）２
令
Ｃ＝ＡＰＡ，ＬＬ＝ＡＰＬ
节，中问结果在计算结束后从内存中提出查看。３１主程序．
ＧＣＷＪＰＣＪ
Ｚ［＋Ｍ＋ＩＩＡ６］ ≠Ｍ：ｏｏＱＩ：ｔ＝０／￣Ｇ
ＬｌＩｚＩ “ ＬＮ）＝” Ｚ『＋Ｎ＋８＋ＩＩｂｓ：Ｌ（Ｘ：ＡＭ］ ≠Ｎ
ｊＧｏＩｔＸｏ＝０：ｂＩｚＩ “ Ｗ）＝” ＺＩ＋ＬｌＹ：ｓ：Ｃ（Ａ
求点高程中误差Ｍｕｉ（）＝±Ｕ￣ｃ／
２程序功能．程序可进行高程控制网按间接平差计算，计算数
据在扩大储存变量ｚＥ］中按一维存放，了节省内ｊ为存单元，解算中采用了内存单元叠加的办法，即由误

4850计算器计算程序

ZY
切线元素各弧长所对园心角(度） 1.137681652 4.957400286 8.77711892 12.59683755 16.41655619 20.23627482 24.05599346 25.8028 切线方位角施工坐标X 施工坐标Y施工转测量坐标X 施工转测量坐标Y (度） 99.17336673 7635.349 3210.668 100.3110484 104.130767 107.9504856 111.7702043 115.5899229 119.4096415 123.2293602 124.9761667 5.957 25.925 45.777 65.427 84.786 103.768 122.289 130.583 0.059 1.122 3.513 7.221 12.230 18.518 26.056 29.911 7634.341 7630.108 7624.583 7617.790 7609.758 7600.525 7590.131 7585.003 3216.539 3236.082 3255.300 3274.107 3292.419 3310.156 3327.239 3334.812
5.9569 25.957 45.957 65.957 85.957 105.957 125.957 135.103
坐标计算模板主编：黄绍敏曲线测设元素与主点桩号计算（范本）
要素计算值(m)
切线方位角 99.17336673 JD点X坐标 7624.394 JD点Y坐标 3278.506 偏向(左-右+)) 1
曲线长Li
方位角(度） 9.1734 10.3110 14.1308 17.9505 21.7702 25.5899 29.4096 33.2294 34.9762

CASIO-fx4850多交点道路平曲线程序

CASIO-fx485Ø道路综合放样程序6
№
PQX FANG YANG(平曲线放样主程序)（544步）
1
Deg:T“WENJIAN MING”: Q=1Ø5 (T-1) :LbI A:｛J｝: J“1SURU-6YIZB-7FY-8FS”
2
=7 =>Goto F⊿J=8=> Goto 7⊿J=6=> Goto 5⊿LbI 9:﹛ZL﹜: Z“JIGE JIAODIAN”:
8
“SQX QD=”:W =Z［7T+Q+6］-S◢“SQX ZD=”:D =Z［7T+Q+6］+S◢
9
Prog“A”:Isz T:Y=T:X=T:Dsz Z: Goto 8: O“LU KUANG”N“PO DU”: Fixm: Prog“AI”:
1Ø
LbI G:﹛VE﹜: V“YIQI H”E“MIAN XIAQU h”:
6
D= G:﹛YXUHGD﹜: Prog“AF”:“SHURU ZB”: Z［7T+Q+13］=Y“X3”:Z［7T+Q+14］
7
=X“Y3”:Prog“AL”:Prog“AG”:Prog“AH”:Prog“A”:Isz T:Y=T:X=T:Dsz Z: Goto 8:
8
LbI 5:fixm:﹛NOKL﹜: Prog“H”:Prog“AI”:
2
M=6 => Goto 5⊿M=8 => Goto 3⊿LbI 6:Fixm :｛ZL｝: Z“JIZU XIANXIN”: Z＞14 => Goto 6⊿
3
Z［12Ø+Q］=Z: T=3: L“1-JISUAN 2-SHURU”=2 => Goto 4⊿Fixm :｛WGHD｝: Prog“D”:

4850计算程序

ZDZBZBJS（匝道中桩、边桩坐标计算）LbI 0：{EG}：A″XA″：B″YA″：C″CA″：D″1/RA″：E″1/RB″：F″DKA″：G″DKB″LbI 1：{HLR}：H″DKI″：L″DL″：R″DR″：H>G＝＞Goto 3⊿P=（E-D）÷Abs（G-F）：Q= Abs（H-F）：I=P×Q：T=D+IJ=C+（I+2D）Q×90÷π◢M=C+（I÷4+2D）×Q×45÷（2π）：N=C+（3I÷4+2D）Q×135÷（2π）K=C+（I÷2+2D）Q×45÷πX=A+Q÷12×（Cos C+4（Cos M+ Cos N）+2 Cos K+ Cos J）◢Y=B+Q÷12×（Sin C+4（Sin M+ Sin N）+2 Sin K+ Sin J）◢U″XL″=X+L Cos（J-90）◢V″YL″=Y+L Sin（J-90）◢W″XR″=X+R Cos（J+90）◢Z″YR″=Y+R Sin（J+90）◢A=X：B=Y：C=J：D=T：F=H：Goto 1LbI 3：A=X：B=Y：D=E：F=G：C=J：Goto 0说明：XA？、YA？——ZH或HZ的X、Y坐标；CA？——切线方位角；1/RA？——起算点的曲率：直线段代0，曲线段代+1/R（右转）或-1/R（左转）；1/RB？——代终点曲率：即±1/R；DKA？——代起算点里程；DKB？——代待求段的终点里程；DKI？——代所求点里程；DL？——代左边桩距中桩的距离；DR？——代右边桩距中桩的距离；——是回车键符号。

（本程序适合于任何线形的计算）（2）ZXZBZJS（直线段中桩、边桩坐标计算）G″KA″A″XA″B″YA″D″XB″E″YB″：POl（D-A，E-B）：J＜0＝＞J=J+360：⊿F=J◢LbI 0：{KSP}：H=K-G：L=Abs H：Rec（L，F）：X=A+I◢Y=B+J◢Rec（S，F-90）：M″XL″=X+I◢N″YL″=Y+J◢Rec（P，F+90）：V″XR″=X+I◢W″YR″=Y+J◢Goto 0说明：KA——起算点里程 XA、YA——起算点坐标 XB、YB——后视点坐标K ——待求点里程 X、Y——待求点中桩坐标 S——左侧桩距中桩距离XL、YL——左侧桩坐标 P——右侧桩距中桩距离 XR、YR——右侧桩坐标（3）DBXMJJS（多边形面积计算）G=0：A=X：B=Y：C=A：D=B：LbI 0：{XY}：X＜0＝＞G=G+AD-BC：A= AbsG/2：A″A=″◢A=0.0015A：A″AA=″◢⊿G=G+XD-YC：C=X：D=Y：Goto 0说明：X、Y——待求点坐标 A——面积（m2） AA——面积（亩）测点坐标输入完毕后在X中输入-1。

Fx4850线元法测量程序

Fx4850最精简的万能坐标计算程序(线元法）1、JS-XY(计算主程序)Norm↙Z“KM”：D：N“PJ”:Prog“D”↙P=(S-R)÷T:Q= Z-C:I=PQ↙J=F+90Q(I+2R)÷π:J≥360=>J=J-360: ≠>J<0=>J=J+360⊿↙L=F+45Q(I÷4+2R)÷(2π):○=F+135Q(3I÷4+2R)÷(2π):K=F+45Q(I÷2+2R)÷π↙U=A+Q÷12×(cosF+4(cosL+cos○)+2 cos K+COS J):V =B+Q÷12×(sin F+4(sin L+ sin ○)+2 sin K+ sin J): W=J+N:X=U+D cosW:Y=V+D sin W:Fix 3: “XD=”:X▲“YD=”:Y↙本程序可以计算任何线型，包括直线、缓和曲线、圆曲线、不等半径曲线等，由其是匝道各曲线计算。

注：0为数字（零），○为字母（噢）KM-----计算点桩号（输入）D-------计算边桩距离左负，右正，中0（输入）PJ------路线右偏角，正交边桩为90（输入）XD------计算所得X坐标（计算得出）YD-----计算所得Y坐标（计算得出）FW---所求点桩号的切线方位角（计算得出）2、子程序：D（数据源）Z<803.961=>A=46670.1：B=13269.591:C=700：F=84°50′23″：R=0：S=0：T=103.961：Goto 0: ≠>Z<1085.902=> A=46679.45：B=13373.13:C=803.961：F=83°33′23″：R=0：S=0：T=281.941：Goto 0: ≠>Z<1363.366=> A=46711.091：B=13653.29:C=1085.902:F=83°33′23″：R=-1÷1000：S=-1÷1000：T=277.464：Goto 0:LbI 0:Fixm:A:B:C:F:R:S:T: ↙Z<本段线元的终点桩号=>A=起点X坐标：B=起点Y坐标:C=线元起点桩号：F=起点切线方位角：R=线元起点半径：S=线元终点半径：T=本段线元长度：Goto 0: ≠>Z<本段线元的终点桩号=>A=起点X坐标：B=起点Y坐标:C=线元起点桩号：F=起点切线方位角：R=线元起点半径：S=线元终点半径：T=本段线元长度：Goto 0: ≠>Z<？=>A=？：B=？:C=？：F=？°？′？″：R=？：S=？：T=？：Goto 0:Z<？=>A=？：B=？:C=？：F=？°？′？″：R=？：S=？：T=？：Goto 0:LbI 0:Fixm:A:B:C:F:R:S:T: ↙说明：字程序Z〈后面为每段线元的终点桩号，可以把你标段所有线元段都输入：其程序中的字母：A、B为起点坐标；C为起点桩号；F为起点方位角：R、S为起终点半径（以倒数输入，如R60，1÷60）分左右转如左转为-1÷60,右转为1÷60，直线为0；T为线元长度。

卡西欧4850计算器程序

卡西欧fx-4850P万能放线程序主程序：ZBJS↙◢“◢：H◢↙注：↙↙↙↙注：此程序可把全线路要素输入数据库，能执行直线、高程放线正反算。

卵形曲线。

卡西欧fx-4850P直线、各种曲线任意放中线(匝道)1：主程序：TYQXJS“1，KS=＞XY”：“2，XY=＞KS”：N：U“X0”：V“Y0”：o“K0”：G“FWJ”：H“LS”：P“R1”：R“R2”：Q“Y+1；0；Z-1”：C=P-1：D=（P-R）÷（2HPR）：E=180÷π：N=1=＞Goto 1:≠＞Goto 2↙◢“Goto 1↙2：正算子程序：S U B1A=0.1739274226:B=0.3260725774:Z=0.0694318442:L=0.33 00094782:F=1-L:M=1-Z:X=U+W(Acos(G+QEZW(C+ZWD))+Bcos (G+QELW(C+LWD))+Bcos(G+QEFW(C+FWD))+Acos(G+QEMW(C+M WD))):Y=V+W(Asin(G+QEZW(C+ZWD))+Bsin(G+QELW(C+LWD)) +Bsin(G+QEFW(C+FWD))+Asin(G+QEMW(C+MWD))):F=G+QEW(C +WD)+90:X=X+Scos F :Y=Y+SsinF↙3坐标反算子程序 SUB2T=G-90：W=Abs((Y-V)cosT-(X-U)sinT) :S=0:Lbl 0: Prog “SUB1”：L=T+QEW（C+WD）:S=（J-Y）cos L-(I-X)sin L : Abs S＜10-6=＞Goto 1: ≠＞W=W+S:Goto 0↙Lbl 1:S=0:Prog“SUB1”:S=(J-Y)÷sin F↙注：o→代表：欧 0→代表零使用说明：1规定：⑴本程序可执行直线、圆曲线、缓和曲线（完整或非完整型）、复曲线、卵形曲线的任意里程中边桩坐标的正反算。

4850曲线计算程序

ZDXLbI 0Prog “FWJ”I”S”▲LbI 1｛ALP｝︰A”J”︰L︰E=X+Lcos (A+J) ︰F=Y+Lsin (A+J) ︰”X”︰E▲“Y”︰F▲P=1=＞Goto 0︰≠=＞Goto 1︰⊿FWJFix 4｛XYMN｝︰X”X0”︰Y”Y 0”︰M”X1”︰N”Y1”︰P ol(M-X,N-Y)FixmJ＜0=>J=J+360︰⊿”A”︰J→DMS▲“S”︰IHYQXLbI 3Fixm﹛K﹜︰Prog “DAT”Prog “QJS”“X”:X▲“Y”:Y▲Goto 3SQXFixmLbl 6Prog “DAT1”Prog “SJS”“H”:X▲Goto 6SDDMDefm 1LbI 0Prog “DAT1”Prog “SJS”Z[1]=XProg “DAT”Prog “QJS”（此处计算坐标必须是隧道中心点坐标）A=X:B=Y:C=D:F=Z[1]｛Z｝︰Z”DM”FixmLbI 1｛VGEH｝︰G”X”︰H”Y”︰E”H”︰V”JING”M=(G-A)cos C+(H-B)sin C︰L=(H-B)cos C-(G-A)sin C︰L=Abs L︰N=E-F︰”K-GAI”︰M▲Prog “DAT2”{W}:W=1=＞Goto 0︰≠=＞Goto 1︰⊿JSKDefm 3LbI 0﹛K﹜︰Prog “DAT”G=Z[3]Prog “QJS”A=X:B=Y:C=D:FixmZ[3]=G:G=Z[1]:H=Z[2]:｛GH｝︰G”X”︰H”Y”Z[1]=G:Z[2]=H:M=(G-A)cos C+(H-B)sin C:”K-GAI”︰M▲K=K+M:”LICH”:K▲L=(H-B)cos C-(G-A)sin C:”BIAN”︰L▲Goto 0QJSS=B²÷24R︰M=C²÷24R︰N=B÷2-B^3÷240R²︰O=C÷2-C^3÷240R²︰T=(R+M-(R+S)cos J)÷sin J+N︰U=(R+S-(R+M)cos J)÷sin J+O︰L=JRπ÷180-B÷2-C÷2+B+C︰H=K-I︰V=L-C︰H≤V=＞W=T︰Z=A+180︰≠=＞W=U︰Z=A+PJ︰⊿E=X+UcosZ︰F=Y+Usin ZH>0=＞H<B=＞Prog”H”︰⊿⊿H≥B=＞H≤V=＞Prog “y”︰⊿⊿H>V=＞H<L=＞H= L-H︰Prog “H”︰H=K-I︰Y=-Y︰D=2A-D+PJ︰⊿⊿H≤0=＞X=H︰Y=0︰D=A︰⊿H≥L＝＞X=L-H︰Y=0︰D=A+PJ︰⊿H≤V=＞O=A︰≠＝＞O=A+PJ +180︰⊿M=E+Xcos O-YsinO︰N=F+Xsin O+Ycos O﹛GQ﹜︰G”B”︰G=0=＞Q= 0︰≠=＞Q”J”︰⊿X=M+Gcos (D+Q) ︰Y=N+Gsin (D+Q) YX=N+Rsin (((2H-B)÷2R)×180÷π) ︰Y=P (S+R(1-cos (90(2H-B)÷πR))) ︰D=A+P (90(2H-B)÷Rπ)HX=H-H^5÷40R²B²+H^9÷3456R^4B^4︰Y=(H^3÷6RB-H^7÷336R^3B^3+H^11÷42240R^5B^5)P︰D=A+P(90H²÷RBπ)SJSW=C-B︰L=WR︰W＜0=＞G=-1︰L= -L︰≠＝＞G=1︰⊿T=L÷2︰L=WR︰E=T²÷2÷R︰H=K-I︰H≤0=＞F= BH︰⊿H>0=＞F=CH︰⊿Abs H<T=＞F=F+(T-A bsH)²÷2÷RG︰⊿X=F+AJDFixmO=0︰V≠1=＞Goto 2︰⊿Pol(N-P,L-Q) ︰J<S=>Goto 3︰⊿J>T=>Goto 3︰⊿o=1︰U=I-R︰U<0=>”R-WAI”︰≠=>”R- NEI”︰⊿U=Abs U︰U▲Goto 3Lbl 2U=R²-(N-P)²︰U<0=> “WUJIE”Goto 3︰⊿U=√U+Q ︰Pol(N-P,U-Q) ︰J<S=>Goto 3 ︰⊿J>T=>Goto 3 ︰⊿o=1 ︰U=L-U ︰U<0=>”H-WAI “︰≠=>”H-NEI”︰⊿U=Abs U ︰U▲Lbl 3DATK<135200＝＞X=1408︰Y=0︰C=280：B=280︰R=6000︰J=23º51º50º︰A=0︰I=131624.36︰P=-1︰Goto 0︰⊿K<137000＝＞X=4475.753︰Y= -1357.128︰A=336º08º10º︰J=20º29º59º︰B=310︰C=310︰R=5000︰I=13 5290.64︰P=1︰Goto0︰⊿K<158950＝＞X=158082.391︰Y= 10000︰A=11º16º41º︰J=11º16º41º︰B=260︰C=260︰R=7000︰I=15726 5.73︰P=-1︰Goto0︰⊿K<166710＝＞X=16096.745︰Y= 10000︰A=0︰J=13º34º14º︰B=280︰C=280︰R=6000︰I=164 591.784︰P=1︰Goto0︰⊿K≥166710＝＞X=14549.512︰Y= 7357.2034︰A=176º15º7.8º︰J=13º42º56.5º︰B=260︰C=260︰R=7000︰I=1671 80.249︰P=1︰Goto0︰⊿Lbl 0DAT1{K}:K<159649=>A=16.692 ︰B=0.003 ︰C=-0.003 ︰R=20000 ︰I=157950 ︰Goto 4 ︰⊿K<161960=>A=8.772︰B=-0.006 ︰C=0︰R=20000 ︰I=161900 ︰G oto 4 ︰⊿K<164558=>A=8.772︰B=0 ︰C=0︰R=20000 ︰I=164500 ︰Goto 4 ︰⊿K<165700=>A=15.152 ︰R=20000︰B=0.0058 ︰C=-0.002︰I=16 5600 ︰Goto 4 ︰⊿K<166750=>A=12.939︰B=-0.002 ︰C=0.003︰R=20000 ︰I=166700 ︰Goto 4 ︰⊿K<168287=>A=17.439︰B=0.003 ︰C=-0.0057︰R=20000 ︰I=16 8200 ︰Goto 4 ︰⊿Lbl 4DAT2Z=330=> Goto A︰⊿Z=320=> Goto B︰⊿Lbl AP=2.27︰Q=0︰R=7.01︰S=0︰T=60︰Prog “JD”:o=1=>Goto Z︰⊿P=2.27︰Q=0︰R=7.01︰S=60︰T=107º53º55º︰Prog “JD” :o=1=>Goto Z︰⊿P=1.217︰Q=3.273︰R=3.57︰S=107º53º55º︰T=157º35º24º︰Prog “JD”:o=1=>Goto Z︰⊿Lbl BP=2.27︰Q=0︰R=7.01︰S=0︰T=60︰Prog“JD”:o=1=>Goto Z︰⊿P=2.27︰Q=0︰R=7.01︰S=60︰T=107º53º55º︰Prog “JD” :o1=>Goto Z︰⊿P=1.217︰Q=3.273︰R=3.57︰S=107º53º55º︰T=157º35º24º︰Prog “JD”:o=1=>Goto Z︰⊿Labl Z程序使用说明此程序适应于CASIO 4850计算器ZDX程序（支导线计算）X0，Y0测站点坐标X1，Y1后视点坐标A 反算方位角S 反算距离J 观测左角L 观测距离X，Y 计算坐标FWJ程序（方位角计算）X0，Y0测站点坐标X1，Y1后视点坐标A 反算方位角S 反算距离HYQX程序（基本型曲线计算）K 计算里程B 计算边距J 右偏角（90为右侧，270为左侧，支持任意角）X，Y 计算坐标SQX程序（竖曲线计算）K 计算里程H 中桩高程SDDM程序（隧道断面计算）K 计算断面所在里程，半径越小，里程越需要准确，7000半径里程偏差十米，误差约1cm，直线无误差。

casio4850全线坐标计算程序

该程序是本人自己编的适用于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 ZP≤YH点桩号＝＞Prog“HYQ.SUB”：Goto ZP≤HZ点桩号＝＞Prog“HQ2.SUB”：Goto Z（以上各字母含义详（HQ1.SUB、HYQ.SUB、HQ2.SUB）坐标计算子程序说明）该主程序是简单的示例说明，具体的应根据图纸的线型组合，随机应变的做顺序调整，并将其输入Casio4850计算器。

4850线路计算程序

4850线路计算程序1：已知线外任意点坐标，求对应线路里程在缓和曲线上，要计算任意里程的法线方向及任意宽度的边线坐标，非常简单。

但要计算任意一个已知坐标点，是对应哪一个里程法线方向上的点，就有一些困难。

很难推导一个这样的计算公式。

唯一的方法“渐进”，如果手工计算这可不是一个好方法。

但在有CASIO系列可编程计算器，如：FX-4500的情况下就变的非常简单了。

亦可用于直线和圆曲线的计算。

首先在缓和曲线上任选一点A为起始点，计算该点的坐标和切线方位角，通过坐标反算求起始点A与计算点B的方位角和距离，B点肯定对应A点切线方向上有一个垂足C点，把三点看成一个直角三角形，通过解直角三角形计算AC的距离，当该距离大于某一数值，如0。

001m，A点里程加AC的距离等于C点的里程，回到开始重新进入新一轮的计算，如果AC的距离小于某一规定值，则计算C点的里程与BC的距离即可。

求对应线路里程程序：主程序QLC (已知坐标求里程)Lb1 0：｛LＤE｝：Prog XH：Goto 0子程序：XH (循环)L1 Lb1 1L2 Norm：Prog LYYD：L3 PO1(D-X,E-Y)：W≤0=> W=W+360⊿L4 Z=W-I：A=V×cos Z：L=L+AL5 Abs A≥0.001=>Goto 1：≠=>B=V×sinZ：Fix 3：“FXJL=” ◢L6 L:Fix3:“DYLC=” ◢程序中字母代表D 任意点X坐标，E 任意点Y坐标，DYLC 对应里程，FXJL 中线法线距离。

程序中有坐标反算功能。

使用方法：只需输入计算点坐标、和较为接近的桩号。

桩号越接近计算速度越快2：逐桩坐标计算2.1编制方法：线路坐标程序是按照平曲线为单元，直线部分归属在曲线两端的方法，把整段路线分段装进数据库，根据桩号判断采用数据通过共用程序，进行任意点的坐标计算，在坐标转换示意土，第一直线段，是通过方位角和距离直接计算大地坐标，第一缓和曲线和圆曲线段，是先计算任意点切线支距和方位角然后转换大地坐标，第二缓和曲线段和直线段是先计算任意点切线支距和方位角。

4580圆曲线综合计算程序

fx--4850P圆曲线计算套件程序SJDYQX(实交单圆曲线主程序)Mcl: R:Z:H:V:W←Prog "ZBFS"←K=J:P=I←Prog "ZBFS"←R>800⇒S=20:≠>S=10△R<150⇒S=5△ G=J-K:G<0⇒G=G+180:≠>G=G-180△Q=Abs( G÷2):H=H÷2←Prog "YQXZD"←XJDYQX(虚交单圆曲线主程序)Mcl: R:Z:H:V:W←Prog "ZBFS"←K=J:P=I←Prog "ZBFS"←S=I：L=J：G=J-K:G<0⇒G=G+180:≠>G=G-180△E:F:Pol(E-C,F-D):J<0⇒J=J+360△E=J-L+180:E>360⇒E=E-360△E<0⇒E=E+180:≠>E=E-180△G=G+E:Q=Abs( G÷2):S=Ssin(Abs E)÷sin (180-2Q):"X.GJ":A=A-Scos K ▲"Y.GJ":B=B-Wsin K ▲ H=H÷2:P=P+S：R>800⇒S=20:≠>S=10 △ R<150⇒S=5　△ ←Prog "YQXZD"←YQXZD(圆曲线主点计算)"T":T=Rtan Q ▲ "L":L=2RQπ÷180▲"E":E=R(1÷cos Q-1)▲N=(Int(Z÷S))S:V≠0⇒U=(tan-1(W÷V))÷V △←Lbl 1： G<0⇒F=-90:"ZZJ":2Q▲≠>F=90:"YZJ":2Q ▲ △ "ZY.ZH":Z+P-T▲"QZ.ZH":Z+P-T+L÷2▲ "YZ.ZH":Z+P-T+L▲ "ZD.ZH":Z+P-2T+L+I ▲Z[1]=A+Tcos K:"ZY.X":Z[1]+(w÷2)cos(K-F)▲Z[2]=B+Tsin K:"ZY.Y":Z[2]+(w÷2)sin(K-F)▲ "QZ.X":A+(E+W÷2)cos(K-F+G÷2)▲ "QZ.Y":B+(E+W÷2)sin(K-F+G÷2)▲:Z[3]=A+Tcos J:"YZ.X":Z[3]+(w÷2)cos(J+F) ▲ Z[4]=B+Tsin J:"YZ.Y":Z[4]+(w÷2)sin(J+F) ▲Lbl 2:X=A+(R+E)cos(K-F+G÷2):Y=B+（R+E）sin (K-F+G÷2):M=K+F：H=0⇒Goto 5 △ G>0⇒Z[5]=0:Z[6]=W:≠>Z[5]=W:Z[6]=0 △"X.QZ":A+Pcos K+Hcos (K+90)▲ "Y.QZ":B+Psin K+Hsin (K+90)▲"X.QY":A+Pcos K+Hcos (K-90)▲ "Y.QY":B+Psin K+Hsin (K-90)▲"ZY.XZ":Z[1]+(H+Z[5])cos(K+90) ▲ "ZY.YZ":Z[2]+(H+Z[5])sin (K+90) ▲ "ZY.XY":Z[1]+(H+Z[6])cos(K-90)▲ "ZY.YY":Z[2]+(H+Z[6])sin (K-90) ▲Lbl 3:Z[7]=X+(R+H)cos(K+F+G÷2):Z[8]=Y+(R+H)sin(K+F+G÷2):Z[9]=X+(R-H-W)cos(K+F+G÷2):Z[10]=Y+(R-H-W)sin (K+F+G÷2):G>0⇒"QZ.XZ":Z[7] ▲ "QZ.YZ":Z[8] ▲ "QZ.XY":Z[9] ▲ "QZ.YY":Z[10] ▲:≠>"QZ.XZ":Z[9]▲"QZ.YZ":Z[10] ▲ "QZ.XY":Z[7] ▲ "QZ.YY":Z[8] ▲ △ ←Lbl 4:"YZ.XZ":Z[3]+(H+Z[5])cos(J-90)▲"YZ.YZ":Z[4]+(H+Z[5])sin(J-90)▲"YZ.XY":Z[3]+(H+Z[6])cos(J+90)▲ "YZ.YY":Z[4]+(H+Z[6])sin(J+90)▲"X.ZZ":A+Icos J+Hcos(J-90) ▲ "Y.ZZ":B+Isin J+Hsin (J-90)▲"X.ZY":A+Icos J+Hcos(J+90) ▲ "Y.ZY":B+Isin J+Hsin (J+90)▲Lbl 5：Prog "JKHHJS"←YQXZB(圆曲线坐标计算)Lbl 1: N=N+S:N>Z+L⇒Goto 4 △ ←Lbl 2: "ZH.P":N ▲ G>0⇒P=N-Z:≠>P=Z-N △ ○=180P÷R÷π:"X.P":X+(R-W÷2)cos (M+○)▲"Y.P":Y+(R-W÷2)sin (M+○)▲Lbl 3:H=0⇒Goto 1 △Z[7]=X+(R+H)cos(M+○):Z[8]=Y+(R+H)sin (M+○):Z[9]=X+(R-H-W)cos(M+○):Z[10]=Y+(R-H-W)sin(M+○):G>0⇒"X.ZP": Z[7]▲ "Y.ZP":Z[8] ▲ "X.YP":Z[9] ▲ "Y.YP": Z[10] ▲≠>"X.ZP":Z[9] ▲ "Y.ZP":Z[10] ▲ "X.YP": Z[7] ▲ "Y.YP":Z[8] ▲ △Goto 1 △ ←Lbl 4:N=N-S ←Lbl 5: Prog "JKHHJS" ←JKHHJS(加宽缓和计算含直线段)Z[12]=Z[12]+1:Z[12]=2⇒Goto 5△ C=A:D=B:A=A+Pcos K:B=B+Psin K:K=K+180:W=0⇒Goto 1:≠>"X.HQ":Z[1]=A+(P-T-V)cos K ▲"Y.HQ":Z[2]=B+(P-T-V)sin K ▲ "X.HQZ":Z[1]+Hcos (K-90) ▲"Y.HQZ":Z[2]+Hsin (K-90) ▲ "X.HQY":Z[1]+Hcos(K+90) ▲"Y.HQY":Z[2]+Hsin(K+90)▲ "X.HZ":Z[1]=C+(T+V)cos J ▲ "Y.HZ":Z[2]=D+(T+V)sin J ▲"X.HZZ":Z[1]+Hcos (J-90)▲ "Y.HZZ":Z[2]+Hsin (J-90)▲"X.HZY":Z[1]+Hcos (J+90)▲ "Y.HZY":Z[2]+Hsin (J+90)▲ △←Lbl 1: N=N+S:N≥Z+P-T⇒Goto 3 △ "ZH.P":N ▲ N>Z+P-T-V⇒○=(N-(Z+P-T-V))tan (U(N-(Z+P-T-V))):≠>○=0 △ H≠0⇒Goto 2:≠>"X.P":A+(N-Z)cos K ▲ "Y.P":B+(N-Z)sin K ▲ △ Goto 1 △←Lbl 2: G>0⇒Z[7]=A+(N-Z)cos K:"X.P":Z[7]+(○÷2)cos(K+90) ▲Z[8]=B+(N-Z)sin K:"Y.P":Z[8]+(○÷2)sin(K+90) ▲"X.ZP":Z[7]+Hcos(K-90) ▲ "Y.ZP":Z[8]+Hsin (K-90 )▲"X.YP":Z[7]+(H+○)cos (K+90)▲"Y.YP":Z[8]+(H+○)sin (K+90)▲≠>Z[7]=A+(N-Z)cos K:"X.P": Z[7]+(○÷2)cos(K-90)▲ Z[8]=B+(N-Z)sin K:"Y.P":Z[8]+(○÷2)sin(K-90)▲"X.ZP":Z[7]+(H+○)cos(K-90) ▲ "Y.ZP":Z[8]+(H+○)sin(K-90 )▲ "X.YP":Z[7]+Hcos (K+90)▲"Y.YP":Z[8]+Hsin (K+90)▲△ Goto 1 △← Lbl 3:N=N-S:Z=Z+P-T←Lbl 4: Prog "YQXZB"←Lbl 5: Z=Z+L:A=C:B=D：V≠0⇒U=(tan-1(W÷V))÷V：≠>U=0 △←Lbl 6:N=N+S:N≥Z+I-T⇒Goto 8 △"ZH.P":N ▲N<Z+V⇒○=(Z+V-N)tan (U(Z+V-N)):≠>○=0 △ H≠0⇒Goto 7:≠>"X.P":A+(N-Z+T)cos J ▲"Y.P":B+(N-Z+T)sin J ▲ △ Goto 6 △ ←Lbl 7: G>0⇒Z[7]=A+(N-Z+T)cos J:"X.P":Z[7]+(○÷2)cos(J+90) ▲Z[8]=B+(N-Z+T)sin J :"Y.P":Z[8]+(○÷2)sin (J+90) ▲"X.ZP":Z[7]+Hcos (J-90) ▲ "Y.ZP":Z[8]+Hsin(J-90) ▲ "X.YP":Z[7]+(H+○) cos (J+90) ▲ "Y.YP":Z[8]+(H+○)sin (J+90)▲ ≠>Z[7]=A+(N-Z+T)cos J:"X.P":Z[7]+(○÷2)cos(J-90) ▲ Z[8]=B+(N-Z+T)sin J:"Y.P":Z[8]+(○÷2)sin (J-90)▲ "X.ZP":Z[7]+(H+○)cos(J-90)▲ "Y.ZP":Z[8]+(H+○)sin(J-90)▲"X.YP":Z[7]+Hcos(J+90)▲"Y.YP":Z[8]+Hsin(J+90)▲ △ Goto 6 △←Lbl 8:"END"▲←X.ZZ、Y.ZZ=终点边左（x，y）X.ZY、Y.ZY=终点边右（x，y）ZY.XZ、ZY.YZ=直圆边左（ｘ，ｙ）ZY.XY、ZY.YY=直圆边右（ｘ，ｙ）QZ.XZ、QZ.YZ=曲中边左（x，y）QZ.XY、QZ.YY=曲中边右（x，y）YZ.XZ、YZ.YZ=圆直边左（ｘ，ｙ）YZ.XY、YZ.YY=圆直边右（ｘ，ｙ）X.P、Y.P =所求点中（x，y）XP.Z、YP.Z=所求边左桩点（x，y）XP.Y、YP.Y=所求边右桩点（x，y）单独显示内容：X.GJ、Y.GJ= 公交点（x，y）（在虚交单圆曲线中显示）。

CASIOfx4850万能计算程序.

CASIO fx4850万能计算程序[kaixin100sc正式版]卡西欧计算器程序2009-10-13 17:59:51 阅读236 评论0 字号：大中小经本人反复测试与验证,该程序可说是一个绝无仅有的万能程序,对于匝道当然也非常适用,匝道在立交工程中常常用到，由于其采用线型多样，计算复杂，因此一份适用的计算程序是不可缺少的，本文提供一份万能匝道计算程序，该程序简单，容易上手，输入参数少。

一、说明在高速公路和城市建设中常常运用到互通式立交工程，互通式立交由于曲线元的复杂性和多样性，给施工中的测量定位工作带来困难，手工计算已经不能满足要求，时下CASIO系列计算器在土木工程中已经得到普遍运用，为便于对立交工程匝道的准确测量定位，本人根据多年的测量工作经验，编写出一份适合CASIO fx-4850匝道计算程序，该程序界面好，能计算出匝道各种线元的中线坐标和边线坐标。

二、程序清单QXZDSJK (程序运行子程序,可计算n条线路)T=1=>Prog “SJK1”：Prog “QXZD”：Goto 0⊿(回车)T=n=>Prog “SJKn”：Prog “QXZD”⊿(回车)Lbi 0SJK (子程序--如果有多条线路，可改程序名为SJK1...到SJKn){H}：H“DKI”(回车)H≤曲线元终点里程=>A=曲线元起点坐标：B=曲线元终点坐标：C=曲线元起点正切线方位角：D=曲线元起点曲率：E=曲线元终点曲率：F=曲线元起点桩号：G=曲线元终点桩号：Goto 0⊿(回车)H≤曲线元终点里程=>A=曲线元起点坐标：B=曲线元终点坐标：C=曲线元起点正切线方位角：D=曲线元起点曲率：E=曲线元终点曲率：F=曲线元起点桩号：G=曲线元终点桩号：⊿(回车)Lbl 0(说明：一条线路中有多少个线元就要添加多少个进入数据库SJK中）QXZD (程序名)A"XA":B"YA": C"CA": D"RA":E"RB":F"DKA":G"DKB":Deg (回车)H>G=>{H}:Goto 0:≠>H<F=>{H}:Goto0:ΔL=(H-F)÷4:I=180L÷πD:U=90LL(D-E)÷πDE(G-F):T=C+I+U: W=C+2I+4U:K=C+3I+9U:O=C+4I+16U: ｛ZQ｝:Z"BJ(+0+)" :P=O+Q"JIAJIAO(-0+)":X=A+L(cosC+4(cosK+cosT)+2cosW+cosO)÷3+ZcosP:Y=B+L(sinC+4(sinK+sinT)+2sinW+sinO)÷3+Zsin P (回车)"X=":X▲"Y=":Y▲I=0:J=0:Pol（X-M"X0"，Y-N"Y0"）：J<0＝＞J=J+360：Δ“FWJ=”:J←DMS◢“I=”：I : Lbl 0三、变量及符号说明说明XA：YA：CA――――曲线元起点X、Y坐标及起点正切线方位角RA：RB――――――曲线元起点及终点曲率半径(左偏负，右偏正)DKA：DKB：―――曲线元起点和终点桩号DKI――――――――曲线元中待求点桩号Z―――――――――边桩距中线平距X：Y―――――――待求点中边桩桩号的X，Y坐标X0，Y0――――――置镜点的X，Y坐标FWJ――――――――待防样点的方位角L=―――――――――待防样点到置镜点的距离 Q"[JIA(-0+)]"--------------中点切线与两边桩连线的夹角(左负中0右正)四、示例龙门互通立交D匝道线型及相关数据如下(里面卖些关子,就是曲线元的半径没有给足,你会计算吗?)：点号桩号ＸY 方位角起点半径，止点半径QD, K0000 , 2777862.564 ,544381.888 ,295"56"01″1291,1291YH, K0081.721 ,2777900.604 ,544309.576 ,299"33"37″1291,250HY,K 122.038 ,2777921.779 ,544275.286 ,305"04"30″250,250YH, K 277.085,2778043.376 ,544183.130 ,340"36"33″250,10^45GQ,K 317.085 ,2778081.751 ,544171.883 ,345"11"34″10^45,-135HY,K 367.085 ,2778129.138 ,544156.172 ,334"34"57″-135,-135ZD, K 432.211 ,2778179.092 ,544115.375 ,306"56"31″五、计算步骤1、QD－YH曲线元，该段为圆曲线运行程序，按照提示输入曲线元QD纵横坐标及方位角，半径1291,QD和YH点桩号，再输入该曲线元内所求点桩号即可显示该中线坐标和切线方位角，可以输入YH点桩号K0 81.721验证计算出的结果，之后会提示输入“Z”值，即边线宽度(不分左右,均为正值)，"Q"值既中线与边桩连线的夹角(左负右正中0)。

(CASIO-4850计算器)曲线任意里程中边桩坐标正反算专用程序

曲线任意里程中边桩坐标正反算(CASIO-4850计算器)专用程序一、程序功能及原理1.功能说明：本程序由一个主程序(TYQXjs)和四个子程——正算子程序(SUB1)、反算子程序(SUB2)等构成，可以根据曲线段——直线、圆曲线、缓和曲线（完整或非完整型）的线元要素（起点坐标、起点里程、起点切线方位角、线元长度、起点曲率半径、止点曲率半径）及里程边距或坐标，对该曲线段范围内任意里程中边桩坐标进行正反算。

本修改版程序既可实现正算全线贯通，亦可实现反算全线贯通。

本程序也可以在CASIO fx-4800P计算器运行，但由于容量较小，不建议使用。

2．计算原理：利用Gauss-Legendre 5点通用(可查阅<测绘工程>2006年第二期，由武汉大学花向红所著---道路平面测设计算的通用公式)公式正算线路中边桩坐标、线外测点至曲线元起点和终点的垂距的符号是否相异（即Dca×Dcb<0=>该测点在其线元内）进行判断并利用该线元要素反算中桩里程、支距，最后计算出放样数据。

二、源程序1.主程序(TYQXjs)Lbl 0：Defm 50:"1.SZ => XY,2.XY => SZ"：N：N=1=>Goto 1：≠>Goto 2Δ←Lbl 1：{SZ}：SZS≤Z[18]=>J=1：Prog “DA T1”：Goto 3ΔS≤Z[26]=>J=2：Prog “DA T1”：Goto 3ΔS≤Z[34]=>J=3：Prog “DA T1”：Goto 3ΔS≤Z[42]=>J=4：Prog “DA T1”：Goto 3Δ……………S≤Z[8（N+1）+2]=>J=n：Prog “DAT1”：Goto 3ΔLbl 3：W=Abs(S-O)：Prog "SUB1"："XS=":X=X◢"YS=":Y=Y◢”FWI=”: F=F-M :F→DMS◢Goto 4←┘Lbl 2：{XY}：XY：Z[4]=X：Z[5] =Y：N=0Lbl A：Isz N：N=J: Prog" DA T1": Z[11]=U：Z[12]=V: Z[9]=G:A= Z[9]-M:J=J+1: Prog "DA T1": Z[13]=U: Z[14]=V: Z[10]=G:B= Z[10]-M:Prog "ZX1" ：Z[6]×Z[7]<0=>J=N：Prog DA T1" ：Goto BΔA=Z[9]+M：B=Z[10]+M：Prog "ZX1" ：Z[6]×Z[7]<0=>J=N：Prog "DA T1" ：Goto B：≠> Goto AΔLbl B：Prog "SUB2"："S=":S=O+W◢"Z=":Z=Z◢Goto 2Lbl 4:J=0:I=Pol(X-Z[1],Y-Z[2]):F=J:F<0=>F=F+360Δ”DIST=”: I◢”FW=”:F→DMS◢Goto 12. 正算子程序(SUB1)A=0.1184634425:B=0.2393143352:N=0.2844444444:K=0.0469100770:L=0.2307653449:Z[3]=0.5:X=U+W(Acos (G+QEKW(C+KWD))+Bcos(G+QELW(C+LWD))+Ncos(G+QEZ[3]W(C+Z[3]WD))+Bcos(G+QE(1-L)W(C+(1-L)WD))+Acos(G+QE(1-K)W(C+(1-K)WD))): Y=V+W(Asin(G+QEKW(C+KWD))+Bsin(G+QELW(C+LWD))+ Nsin(G+QEZ[3]W(C+Z[3]WD))+Bsin (G+QE(1-L)W(C+(1-L)WD))+Asin (G+QE(1-K)W(C+(1-K)WD)))：M”ANG=”: F=G+QEW(C+WD)+M：X=X+ZcosF：Y=Y+ZsinF3. 反算子程序(SUB2)M”ANG=”:T=G-M：W=Abs((Y-V)cosT-(X-U)sinT)：Z=0：Lbl 0：Prog "SUB1"：L=T+QEW(C+WD)：Z=(Z[5]-Y)cosL-(Z[4]-X)sinL：AbsZ<1E-6=>Goto1：≠>W=W+Z：Goto 0Δ←┘Lbl 1：Z=0：Prog "SUB1"：Z=(Z[5]-Y)÷sinF4. 垂距计算子程序(ZX1)Z[6]=(Z[5]-Z[12])COS A-(Z[4]-Z[11])SIN AZ[7]=(Z[5]-Z[14])COS B-(Z[4]-Z[13])SIN B5.曲线元要素数据库：DA T1Lbl1:J=1=>U=Z[8]:V=Z[9]:O=Z[10]:G=Z[11]:H=Z[12]:P=Z[13]:R=Z[14]:Q=Z[15]ΔJ=2=>U=Z[16]:V=Z[17]:O=Z[18]:G=Z[19]:H=Z[20]:P=Z[21]:R=Z[22]:Q=Z[23]ΔJ=3=> U=Z[24]:V=Z[25]:O=Z[26]:G=Z[27]:H=Z[28]:P=Z[29]:R=Z[30]:Q=Z[31]ΔJ=4=> U=Z[32]:V=Z[33]:O=Z[34]:G=Z[35]:H=Z[36]:P=Z[37]:R=Z[38]:Q=Z[39]Δ……………………………J=N=> U=Z[8N]:V=Z[8N+1]:O=Z[8N+2]:G=Z[8N+3]:H=Z[8N+4]:P=Z[8N+5]:R=Z[8N+6]:Q=Z[8N+7]Δ(注：如有多个曲线元要素继续添加入数据库DAT1中)Lbl 2 ：U"X0"：V"Y0"：O"S0"：G"F0"：H"LS"：P"R0"：R"RN"：Q：C=1÷P：D=(P-R)÷(2HPR)：E=1 80÷π三、使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线路往左偏时，Q=-1；当线路往右偏时，Q=1；当线路为直线时，Q=0。

CASIO4850公路线元法积分公式计算器程序V11版

请使用WORD的Times New Roman字体，若使用其它字体，字母O与数字0是区别不出的.切记切记!积分公式计算的线元法CASIO4850计算器程序第十一版1 (主程序)(使用时请运行这个)L1 LbI 0:{ QDV}L2 U<4=> Prog “4”△L3 U=5=> Prog “2”△L4 U=6=> Prog “2”△L5 Int U=7=> D=0:V=0△L6 U=8=>{M}: PSW:D=(B+M“01”P)÷Cos Abs(90-AbsV): Q=W+S√(1-M2)△L7 U>5=> Prog “4” △L8 LbI2:IntU=9=> PS: Prog“F”:A=R+90: Prog“9”:Goto 2△L9 LbI3:U=10=>PSWB: Prog“F”: Pol ( P-W,S-B):J<0=>J=J+360:△A=J: Prog“9”:Goto 3△L10 QDVHEFACZG: Frac Y=0=> Prog “3”: △O: Prog “10”L11 U=2=>T=T+ Int (-5+10 Ran# )÷1000: K=K+ Int (-5+10Ran# )÷1000: Prog “8”: Prog “5”△L12 U=3=> I= -S+2S Ran# :J=√(S2-I2): I= Int I÷1000:J= Int (-J+2J Ran# )÷1000:T+I: Pause 0:K+J: Pause 0:I: Pause 0:J◢√(I2+J2) ◢△L13LbI 1: IntU=7=> {DV}: Prog“10”: S=D:D=B÷2÷Sin Abs V:J=R:U=7.1=>D=B÷2:J=R+V-90: △ Prog “6”: Prog “8”: Prog “5”:D=2D:J=R-180: U=7.1=> J=R+V+90: △ Prog “6”: Prog “8”: Prog “5”: V≤0=>V=V+180: ≠>V=V-180: △D=S:Goto 1△L14 U<5=>Q=Q+B: △Goto 010（子程序）（为共用计算中桩/加转角/用正算计算边桩/反算/显示的过程）L1 Prog “F”:U=11=>R=0:△J= R+V: Prog “6” : Prog “8” : Prog “5”9 （子程序）（为点垂桩垂距/涵交桩交角的共同运行部分）L1 Prog“7”: Prog“8”: Abs(J-R)<1=>Q=Q-I:≠> Q=Q+I:△I<0.001=>T=P:K=S: Prog “8”: “OK”: Q: Pause 0:I: Pause 0: J=J-R:J>180=>J=-J+180:△J→DMS◢U=9.1=> Prog “11”: △△ Prog “4”7 （子程序）（两点斜直线交点部分）(也可作为两点斜直线交点主程序独立运行)L1E=Tan A:F= Tan R:Y=(EP-S-FT+K) ÷(E-F): Pause 0:O=E(Y-P)+S5 (子程序)( 显示部分)L1T: Pause 0:K: Pause 0:I: Pause 0:J→DMS◢F (子程序) (中桩计算部分)L1 RAD:I= Abs(Q-H):J= Abs(G-H):R=Aπ/180+CI+I2(Z-C)/2/J:T=E+∫(Cos (Aπ/180+CX+X2(Z-C)/2/J),0,I,4):K= F+∫(Sin (Aπ/180+CX+X2(Z-C)/2/J),0,I,4):Q>G =>R=R+π:△DEG:R=180R/πF2 (子程序) (中桩计算部分)(使用时把F2改为F;同时把F改为F1)(这个子程序计算速度快,但在缓曲长到2或3百米且半径小到6或7十米时,计算误差有几厘米)L1 Defm 7L1 I=Q-HL2 J= (Z-C)I÷(G-H)L3 Z[1]=A+90(J÷8+2C)(I÷8)÷πL4 Z[2]=A+90(J÷4+2C)(I÷4)÷πL5 Z[3]=A+90(3J÷8+2C)(3I÷8)÷πL6 Z[4]=A+90(J÷2+2C)(I÷2)÷πL7 Z[5]=A+90(5J÷8+2C)(5I÷8)÷πL8 Z[6]=A+90(6J÷8+2C)(6I÷8)÷πL9 Z[7]=A+90(7J÷8+2C)(7I÷8)÷πL10 R=A+90(J+2C) I÷πL11 T=E+ AbsI(Cos A+4(CosZ[1]+ CosZ[3]+ CosZ[5]+ CosZ[7])+2(CosZ[2]+ CosZ[4]+ CosZ[6])+ CosR)÷24L12 K=F+ AbsI(Sin A+4(SinZ[1]+ SinZ[3]+ SinZ[5]+ SinZ[7])+2(SinZ[2]+ SinZ[4]+ SinZ[6])+ SinR)÷248(子程序) (座标反算部分)(也可作为座标反算主程序独立运行, 方位角按RCL R SHIFT °′″ 显示)L1J=TK: Pol ( T-Y,K-O): Pause 0:J<0=>J=J+360△6(子程序) (边桩和座标正算计算部分)(也可作为座标正算主程序独立运行)L1TK:T=T+ Rec (D,J): Pause 0:K=K+J2(子程序) (推算开挖/填筑桩部分){平原微丘一般推算1-2次,山岭重丘一般4-5次,能达到5厘米以内的精度,看你的熟练程度} L1 {WSPB}：I= -W+ S◢L2 D= Abs IP+B◢L3 {D}:D3(导线点数据库)(数据库中的数字请按你所在公路的数据替换)L1 Y=21=>Y=36640.666:O=65111.666△L2 ……(省略号表示类同上行, 并非真在程序中输入省略号. 注意: 第一个Y为导线点编号, 第二个Y为导线点X座标,O为导线点Y座标, 数据库中的数字请按你所在公路的数据替换. 但注意: 导线点编号必须没字母也没小数;导线点X座标则必须有小数,若刚好没有小数,请加上0.0001并不影响精度)4(平曲线数据库) (数据库中的数字请按你所在公路的数据替换)L1 Q≤444.978=>H=374.978:E=36771.643:F=64929.942:A=151°29′38″:C=0: G=444.978:Z=-600-1: Goto 1△L2 Q≤811.819=>H=444.978: E=36710.800:F=64964.534:A=148°09′06″:C=-600-1: G= 811.819:Z=-600-1: Goto 1△L3 Q≤881.819=> H=811.819:E=36475.596:F=65238.599:A= 113°07′15″:C=-600-1: G=881.819:Z=0: Goto 1△L4 ……(省略号表示类同上行,并非真在程序中输入省略号)L5 LbI 1 (注意:这行始终是平曲线数据库的最后一行, 之前行都是相同的) (为更好的理解这个平曲线数据库的结构,请对照上面数据看下图)11 （子程序）（双圆垂轴隧道断面开挖超欠测量部分）A:实测高程 H:设计高程I:实测宽度 B:超欠值(设棱镜半径0.03)L1 {AHI}:A≤H=>B=110:≠>A≤H+2.75+1.056=>B=√((A-(H+1.905))2+(I+1.813)2) +0.03-4.5 : ≠> B =√ ((A-(H+2.75))2+I2 )+0.03-2.5:△△B◢平曲线数据库的数据图示:(图画得不好,现丑了.图示是对称缓曲,但线元法程序对任意公路线形都适用的)ZH HZ 0+374.978X=36771.643 0Y=64929.942 HY151°29′38″ 0+444.9780 X=36710.800Y=64964.534-600-1 X=36475.596Y=65238.599113°07′15″-600-1使用说明:{本程序是针对CASIO4850计算器的，其它的计算器请作相应修改}{有方框的字为一个字符，不是几个字母，为各种计算状态名或函数名，如RAD与Cos等等，∫为积分符号，L1\ L2\ L3\L4等为程序行号不要在程序中输入, ……省略号表示类同上行,并非真在程序中输入省略号}{如有其它疑问或进行技术交流，请打詹德军的电话015918402837欢迎交流!}运行主程序 1 (按1 FILE进入)U?(判断进入何种计算状态)1为数据库输入原始数据计算中边桩，2为数据库输入计算放线测量记录表，3为数据库输入计算平面位置检查记录表，4为人工输入原始数据计算中边桩，5为人工输入推算开挖/填筑桩状态，6为数据库输入推算开挖/填筑桩状态，7为数据库输入计算斜交斜做涵洞/单跨桥七点放样状态，7.1为数据库输入计算斜交正做涵洞七点放样状态，(斜交斜做涵洞的伸缩缝处的放样也可用)8为数据库输入计算锥坡放样，9为数据库输入反算线外点垂桩垂距，9.1为数据库输入双圆垂轴隧道断面开挖超欠测量，10为数据库输入求线外两点交路线的桩号与夹角，11为计算曲线桥平行布桩的桩/柱放样状态(因为D=0时就是中桩,D≠0时就是边桩,所以不再单独设立中桩或边桩的计算状态,方便程序进一步开发)（U=1）数据库输入原始数据计算中边桩1.U？1(判断进进入何种计算状态)2.Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3.D？中桩到边桩的距离（求中桩输入0）4.V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5.Y?置仪点X座标（也可以输入导线点编号，因为此时使用了导线点数据库）6.O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7.=所求桩号X座标（注:四个结果同时显示于屏幕）=所求桩号Y座标=所求桩号放样距离=所求桩号放样方位角8．B？步长（步长就是间隔长度，当要计算的是大量有相同距离间隔且相同中边距或都是中桩时，步长就是好东西了，计算完第一个桩后所有桩，都只需闭着眼睛连按四次EXE，就是下个桩要计算的所有结果了）9．以后循环第2/3/4/7步（U=2）计算放线测量记录表1．U？2(判断进进入何种计算状态)2．Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3．D？中桩到边桩的距离（求中桩输入0）4．V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5．Y?置仪点X座标（也可以输入导线点编号，因为此时使用了导线点数据库）6．O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7．=所求桩号计算X座标（注:四个结果同时显示于屏幕）=所求桩号计算Y座标=所求桩号计算放样距离=所求桩号计算放样方位角8．=所求桩号复测检查X座标（注:四个结果同时显示于屏幕）=所求桩号复测检查Y座标=所求桩号复测检查放样距离=所求桩号复测检查放样方位角9．B？步长（步长就是间隔长度）10．以后循环第2/3/4/7/8步（U=3）计算平面位置检查记录表1．U？3(判断进进入何种计算状态)2．Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）3．D？中桩到边桩的距离（求中桩输入0）4．V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)5．Y?置仪点X座标（也可以输入导线点编号，因为此时使用了导线点数据库）6．O?置仪点Y座标（当X座标输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）7．=所求桩号设计X座标（注:四个结果同时显示于屏幕，但后两个数据是用不上的） =所求桩号设计Y座标=所求桩号放样距离=所求桩号放样方位角8．？所允许的最大误差（以毫米为单位，比如高速公路的中线允许的最大误差为20MM，就输入20）9．=所求桩号检测X座标（注:四个结果同时显示于屏幕）=所求桩号检测Y座标=所求桩号X轴偏差=所求桩号Y轴偏差10．=所求桩号双轴偏差11．B？步长（步长就是间隔长度）12．以后循环第2/3/4/7/9/10步（U=4）人工输入原始数据计算中边桩1．Q？所求桩号（计算第二个桩开始Q已经自动加了步长，不用手工输入，但仍然会显示你看，并可修改成别的桩号）2．D？中桩到边桩的距离（求中桩输入0）3．V？中桩到边桩方位角与中桩前进切线方位角间的夹角(注意左转为负值,角度按需要以度分秒输入任意输入)4．H?起算点桩号5．E?起算点X座标6．F?起算点Y座标7．A?起算点方位角(以度分秒输入)8．C?起算点曲率(即半径分之一，直线为0，注意左转为负值)9．Z?终算点曲率(即半径分之一，直线为0，注意左转为负值)10.G?终算点桩号11.Y?置仪点X座标12.O?置仪点Y座标13.=所求桩号X座标（注:四个结果同时显示于屏幕）=所求桩号Y座标=所求桩号放样距离=所求桩号放样方位角14.？步长（步长就是间隔长度）15.以后循环第1/2/3/13步（U=5）人工输入原始数据推算开挖/填筑桩状态(平原微丘一般推1-2次，山岭重丘4-5次达到5厘米的精度，要看你的熟练程度)1．W?实测大地高程2．S?变坡点设计高程3．=变坡点设计高程与实测大地高程之高差(正为填方,负为挖方)4．P?填方或挖方的坡比(若1:1.5侧输1.5,无论填挖坡比全为正)5．B?中桩到变坡点宽度6．=推算的中桩到开挖/填筑桩距离7. D?采用的中桩到开挖/填筑桩距离7．Q?所求开挖/填筑桩桩号8．V?中桩到边桩的方位角与中桩的前进切线方位角间的夹角(注意左转为负值,只能输入右边90或左边-90，因为开挖/填筑桩总是垂直的)9．H?起算点桩号10.E?起算点X座标11.F?起算点Y座标12.A?起算点方位角(以度分秒输入)13.C?起算点曲率(即半径分之一，直线为0，注意左转为负值)14.Z?终算点曲率(即半径分之一，直线为0，注意左转为负值)15.G?终算点桩号16. Y?置仪点X座标17. O?置仪点Y座标18.=推算开挖/填筑桩桩号X座标（注:四个结果同时显示于屏幕）=推算开挖/填筑桩桩号Y座标=推算开挖/填筑桩桩号放样距离=推算开挖/填筑桩桩号放样方位角12．以后循环第1/2/3//4/5/6/7/8/18步（U=6）数据库输入原始数据推算开挖/填筑桩状态(平原微丘一般推1-2次，山岭重丘4-5次达到5厘米的精度，要看你的熟练程度)1．W?实测大地高程2．S?变坡点设计高程3．=变坡点设计高程与实测大地高程之高差(正为填方,负为挖方)4．P?填方或挖方的坡比(若1:1.5侧输1.5,无论填挖坡比全为正)5．B?中桩到变坡点宽度6．=推算的中桩到开挖/填筑桩距离7．D?采用的中桩到开挖/填筑桩距离8．Q?所求开挖/填筑桩桩号9．V?中桩到边桩的方位角与中桩的前进切线方位角间的夹角(注意左转为负值,只能输入右边90或左边-90，因为开挖/填筑桩总是垂直的)10．Y?置仪点X座标（也可以输入导线点编号，因为此时使用了导线点数据库）11．O?置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）12．=推算开挖/填筑桩桩号X座标（注:四个结果同时显示于屏幕）=推算开挖/填筑桩桩号Y座标=推算开挖/填筑桩桩号放样距离=推算开挖/填筑桩桩号放样方位角 12．以后循环第1/2/3//4/5/6/7/8/11步推算开挖/填筑桩的图示(以有二级填方边坡的填筑线为例,其它类推)中桩S变坡点高程中桩至该变坡点距离 B地面线（U=7/7.1）为计算斜交斜做涵洞单跨桥 / 斜交正做涵洞七点放样状态(U=7/7.1并能避免曲线内单跨桥涵进出水口长度不一至)(U=7同时适用于,曲线桥扇形布桩的盖梁放样; 曲线桥扇形布桩的支座放样)1.U? 7/7.1进入计算斜交斜做涵洞单跨桥/ 斜交正做涵洞七点放样状态2.Q? 涵洞桩号3.Y?置仪点X座标（也可以输入导线点编号，因为此时使用了导线点数据库）4. O?置仪点Y座标（当Y输入了导线点编号，则此行不显示，因为数据库已经装入了导线点的X与Y座标）5.=涵洞中桩X座标（注:四个结果同时显示于屏幕）=涵洞中桩Y座标=涵洞中桩放样距离=涵洞中桩放样方位角6.D?右侧涵斜长7.V? 右侧转角(进流水轴与中桩的前进切线方位角间的夹角,左负右正)8.= 右侧桩X座标（注:四个结果同时显示于屏幕）=右侧桩Y座标=右侧桩放样距离=右侧桩放样方位角9.B?正的跨径长度(程序自算斜跨径长度)10. =右侧前进侧台身脚X座标（注:四个结果同时显示于屏幕）=右侧前进侧台身脚Y座标=右侧前进侧台身脚放样距离=右侧前进侧台身脚放样方位角11. =右侧后退侧台身脚X座标（注:四个结果同时显示于屏幕）=右侧后退侧台身脚Y座标=右侧后退侧台身脚放样距离=右侧后退侧台身脚放样方位角12.D?左侧涵斜长13.V? 左侧转角(出流水轴与中桩的前进切线方位角间的夹角,左负右正,程序自动按180度转换了,其实不用输入,只是显示观看)14.= 左侧桩X座标（注:四个结果同时显示于屏幕）=左侧桩Y座标=左侧桩放样距离=左侧桩放样方位角15. =左侧前进侧台身脚X座标（注:U=7.1时为后退侧结果，四个结果同时显示于屏幕）=左侧前进侧台身脚Y座标=左侧前进侧台身脚放样距离=左侧前进侧台身脚放样方位角16. =左侧后退侧台身脚X座标（注: U=7.1时为前进侧结果，四个结果同时显示于屏幕）=左侧后退侧台身脚Y座标=左侧后退侧台身脚放样距离=左侧后退侧台身脚放样方位角17.以后循环第6-8/10-17步(当然这种循环只是重新检查而已,按两次AC/ON退出)计算涵洞单跨桥七点放样图示（注: U=7.1时左侧的15步与16步顺序互换）前进方向10步15步8步右侧左侧14步 5步11步16步切线路线（U=8）为数据库输入原始数据计算锥坡放样1．U?8为进入数据库输入原始数据计算锥坡放样2．P？椭圆单侧斜长轴3．S？椭圆单侧斜短轴(当从近桥0点往远桥1点推桩号增加时，S输入负值。