缓和曲线计算小程序(自编)
CASIO任意曲线坐标计算程序及说明
个人测量资料汇集韦习耀2005年8月起缓和曲线（ZH到HY）中桩坐标CASIO4800/4850计算程序A＂ZHX＂：B＂ZHY＂：C＂LS＂：Q：R：K：M＂DK ZH＂↙LB1 1：｛L｝：L＂DK＂：X=A＋（（L－M）－（L－M）5÷（40R2C2））cosQ-(（L－M）3÷(6RC))KsinQ◢Y=B＋(（L－M）-（L－M）5÷(40R2C2))sinQ＋(（L－M）3÷(6RC))KcosQ ◢Goto 1说明：（A，B）为ZH点坐标，LS为缓和曲线长，Q为ZH点放位角，R 为半径，左转K为-1，右转K为1，M为ZH点里程，L为所求点的里程。
缓和曲线（HY到HZ）中桩坐标CASIO4800/4850计算程序A＂HZX＂：B＂HZY＂：C＂LS＂：Q：R：K：M＂DK HZ＂↙LB1 1：｛L｝：L＂DK＂：X=A＋（（－（L－M））－（－（L－M））5÷（40R2C2））cosQ-(（－（L－M））3÷(6RC))KsinQ◢Y=B＋(（－（L－M））-（－（L－M））5÷(40R2C2))sinQ＋(（－（L－M））3÷(6RC))KcosQ◢Goto 1说明：（A，B）为HZ点坐标，LS为缓和曲线长，Q为HZ点放位角，R 为半径，左转K为-1，右转K为1，M为HZ点里程，L为所求点的里程。
另外注意，这个程序是把HZ点当作ZH点，倒过来算，方位角Q 和K均应注意其值。
圆曲线（HY到YH）中桩坐标CASIO4800/4850计算程序A＂HYX＂：B＂HYY＂：C＂LS＂：Q：R：K：M＂DK HY＂↙LB1 1：｛L｝：L＂DK＂：X=A＋2Rsin（90（L－M）÷（πR））cos（Q ＋90（（L－M）＋C）÷（πR）K）◢Y=B＋2Rsin（90（L－M）÷（πR））sin（Q＋90（（L－M）＋C）÷（πR）K）◢Goto 1说明：（A，B）为HY点坐标，LS为缓和曲线长，Q为ZH点放位角，R 为半径，左转K为-1，右转K为1，M为HY点里程，L为所求点的里程。
不等长缓和曲线、竖曲线、直线4850计算器程序
一、曲线坐标计算主程序：QXJSProg”1“：Q”K-ZH”:F”ZHX=”:G”ZHY=”:N”FANG,WEI,JIAO=”:I”(Y1,Z-1)”Lbl 1:Fixm:｛KDJ｝:K”LICHENG=”:D”WAIYI=(Z-,Y+)”:J”JIAJIAO=”:Goto 2Lbl 2:K<Q=>Goto 3:≠>Goto 4Lbl 3:”K-OVER”:Goto1Lbl 4:K>Q+L=>Goto 3: ≠>Goto 5Lbl 5:K≤Q+B=>Goto 6: ≠>Got o 7Lbl 6:S=K-Q:W=B:Prog”2”:O=N+I×90S2÷(ΠRB):O<0=>O=O+360⊿O>360=>O=O-360⊿O”FANGWEI=”X=Ucos N-Vsin N+F+D×cos(O+J)Y=Vcos N+Usin N+G+D×sin (O+J)Goto 1Lbl 7:K≤Q+L-C=>Goto 8:≠>Goto 9Lbl 8:S=K-Q:O=N+I×90(2S-B) ÷(ΠR):O<0=>O=O+360⊿O>360=>O=O-360⊿O”FANGWEI=”Prog”3“：X= Ucos N-Vsin N+F+D×cos(O+J)Y=Vcos N+Usin N+G+D×sin (O+J)Goto 1Lbl 9:S=Q+L-K:W=C: Prog”2“:E=N+IA+180:O=E-I×90S2÷(ΠRC)+180:O<0=>O=O+360⊿O>360=>O=O-360⊿O”FANGWEI=”V=-VX=(U-H) cos E-Vsin E+F+T×cosN+D×cos(O+J)Y=Vcos E+(U-H)sin E+G+T×sinN+D×sin(O+J)Goto 1子程序：（共3个）1、A:R:B”L01”:C”L02”M”M1”=B÷2-B3÷(240×R2)J”M2”=C÷2-C3÷(240×R2)P”P1”=B2÷(24×R)-B4÷(2688×R2)P”P2”=C2÷(24×R)-C4÷(2688×R2)T”T1”=M+(R+P)tan(A÷2)+(Z-P) ÷sinAH”T2”=J+(R+Z)tan(A÷2)+(P-Z) ÷sin AL=B÷2+ C÷2+ARΠ÷180E=(R+(Z+P) ÷2)÷cos(A÷2)-R2、 U=S-S5÷(40R2W2)+S9÷(3456R4W4)-S13÷(599040R6W6):V=I×(S3÷(6RW)-S7÷(336R3W3)+S11÷(42240R5W5))3、 U=Rsin(180(S-B÷2) ÷ΠR)+M:V=I×(R×(1-cos(180(S-B÷2) ÷ΠR))+P)备注：本程序可计算不等长缓和曲线及一般曲线中桩及边桩坐标。
4800P计算器圆曲线、缓和曲线通用程序清单
4800P通用程序清单圆曲线程序：—YQX（拼音）O“Xδ”：H“Yδ”：G“а0”：M“μв”：P“YP,ZP＝+1,-1”：R“R”： L“L C”： C“L h”: Q“ZYD.ZH”：LbI 1：{KD}：S＝K－Q：S＜0=＞S＝0:⊿ S＞L=＞S＝L:⊿V＝S÷R×180÷π：N＝(M－G)÷2：N＜0=＞N＝-N:⊿N＞90=＞N＝180-N:⊿U＝C÷2R×180÷π：C＝0=＞E＝R÷CosN:≠＞E＝(R＋C2÷24R)÷CosN:⊿ P＝1=＞F＝M＋(180－2N)÷2：T＝G＋270：W＝T＋U＋V：Z＝R－D:≠＞F＝M－(180－2N)÷2：T＝G＋90：W＝T－U－V:Z＝R＋D:⊿ Fix 3：X“X”＝O＋E×CosF＋Z×CosW◢Y“Y”＝H＋E×SinF＋Z×SinW◢Goto 1缓和曲线程序：—HHQX（拼音）O“Xδ”：H“Yδ”： G“а0”: Z“Bp g”：P“P=1,-1,2,-2”： R“R”： C“L h”： Q“QSD.ZH”：LbI 0：{KD}：S＝K－Q：J＝Abs P：J＝1=＞Goto 1:≠＞J＝2=＞Goto 2:≠＞P“P=1，-1，2，-2；Ran#”◢⊿ Goto 0：LbI 1：S＜0=＞S＝0:⊿ S＞C=＞S＝C:⊿ Goto 3：LbI 2：S≥0=＞S＝0:≠＞S＝Abs S:⊿ S＞C=＞S＝C:⊿ Goto 3：LbI 3：V＝S2÷(2×R×C)×180÷π： E＝S∧3÷(6×R×C)：F＝S－S∧5÷(40×R2×C2)：Goto 4：LbI 4：P＝-1=＞Goto 5:≠＞P＝1=＞Goto 6:≠＞P＝-2=＞Goto 7:≠＞P＝2=＞Goto 8:⊿LbI 5：U＝G＋270：W＝G－V＋90：Goto 9：LbI 6：U＝G＋90：W＝G＋V＋90：F＝-F：Goto 9：LbI 7：U＝G＋270：W＝G－V＋270：Goto 9：LbI 8：U＝G＋90：W＝G＋V＋270：F＝-F：Goto 9：LbI 9：L＝C÷2－C∧3÷(240×R2)：Z=AbsZ：T＝(R＋C2÷(24×R))×Tan(Z÷2)：M＝O＋(L＋T)×Cos(G＋180):N＝H＋(L＋T)×Sin(G＋180): Fix 3：X“X”＝M＋E×CosU－F×SinU＋D×CosW◢Y“Y”＝N＋E×SinU＋F×CosU＋D×SinW◢Goto 0：程序注译圆曲线程序---YQX（拼音）Xδ？--------------------输入交点纵坐标。
道路设计缓和曲线参数宏代码
道路设计缓和曲线参数宏代码1. 问题描述在道路设计中，为了提高行车的安全性和舒适性，常常需要设置缓和曲线。
缓和曲线是一种具有特定曲率变化规律的曲线，在转向时使车辆更平稳地通过道路的转弯部分。
为了快速且准确地确定缓和曲线的参数，我们可以使用宏代码来进行计算。
2. 缓和曲线的参数缓和曲线的参数包括缓和曲线的长度、曲率变化的速率以及曲线的半径。
通过调整这些参数，可以使得缓和曲线适应不同的道路设计需求。
2.1 缓和曲线的长度缓和曲线的长度是指曲线的起始点到终点的直线距离，通常以米为单位。
长度的选取应该考虑到车辆的速度和转向的舒适度，较高的速度和较大的转向角度需要更长的缓和曲线长度，以确保车辆安全平稳地通过转弯。
2.2 曲率变化的速率曲率变化的速率是指缓和曲线上曲率从起始点到终点逐渐变化的程度。
曲率变化的速率越大，缓和曲线上的曲率变化越快，车辆的转向更加迅速。
曲率变化的速率影响到车辆的进入和离开缓和曲线的过程，合理的速率可以提高车辆通过转弯时的舒适度。
2.3 缓和曲线的半径缓和曲线的半径是曲线在某一点上的曲率半径，通常以米为单位。
缓和曲线的半径决定了曲线的弯曲程度，较小的半径意味着更加急剧的转弯。
在设置缓和曲线的半径时，需要考虑到道路的设计需求和车辆的转向能力，以保证车辆能够安全通过道路。
3. 宏代码实现下面是一个示例的宏代码实现，用于计算缓和曲线的参数。
代码使用Python编写，并使用了适应车辆速度和转向角度的缓和曲线设计方法。
import mathdef calculate_curve_parameters(v, theta):"""计算缓和曲线的参数参数:v: 车辆行驶速度，单位 m/stheta: 车辆转向角度，单位弧度返回值:curve_length: 缓和曲线的长度，单位米curvature_rate: 曲率变化的速率curve_radius: 缓和曲线的半径，单位米"""# 根据车辆速度计算曲线长度curve_length = 0.5 * v ** 2 / 4# 根据转向角度计算曲率变化的速率curvature_rate = 2 * theta / curve_length# 根据曲率变化的速率计算曲线半径curve_radius = 1 / curvature_ratereturn curve_length, curvature_rate, curve_radius# 示例用法v = 10 # 车辆速度为10 m/stheta = math.pi / 6 # 车辆转向角度为30°curve_length, curvature_rate, curve_radius = calculate_curve_parameters(v, the ta)print("缓和曲线的长度:", curve_length, "米")print("曲率变化的速率:", curvature_rate)print("缓和曲线的半径:", curve_radius, "米")4. 结论本文介绍了道路设计中缓和曲线的参数以及如何使用宏代码计算这些参数。
缓、圆曲线程序
缓和曲线、圆曲线计算程序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：起点（终点）X坐标P“Y0”=B-SIN E×F：起点（终点）Y坐标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:△L1的前进方位角V=U+90：L2的前进方位角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输入02左偏为负，右偏为正。
适合隧道测量用的各种小程序
曲线坐标计算程序：Deg：Norm1↓（设置计算器角度单位为度及小数点显示为小数模式）主程序（程序名称：QXZBJS）R ：K：W：C↓（定义变量）K÷2- K3÷240R2→Z[1]：K2÷24R- K4÷2384R3→Z[2] ：(Z[2]+R)tan(Abs(W)÷2)+ Z[1]→Z[3] ：K+ Abs(W)×πR÷180→Z[4]↓Lbl 1↓Abs(D-C)→I↓If D＜C：Then D-C→X：0→Y ：0→Q ：Goto 4：IfEnd↓If I> Z[4]：Then Z[4]-I→Z[5]：0→Z[6] ：0→Q ：Goto 3：IfEnd↓If I≤Z[4]-K：ThenI→J：Else Z[4]-I→J：IfEnd↓If J≥K：Then Goto 2：IfEnd↓(J2×180÷(2πRK))°→Q↓（以度数为单位赋值给Q）J- J5÷(40R2K2)→X：J3÷(6RK)- J7÷(336R3K3)→Y↓If I＞Z[4]-K：Then X→Z[5]：Y→Z[6]：Goto 3：IfEnd↓Goto 4↓Lbl 2↓((I- 0.5K)×180RπR)°→Q↓（以度数为单位赋值给Q）Rsin(Q)+ Z[1]→X：R(1-cos(Q)) + Z[2]→Y↓Goto 4↓Lbl 3↓Z[3](1+cos(W))- Z[5]cos(W)- Z[6]sin(Abs(W))→X↓(Z[3]- Z[5]) sin(Abs(W)) + Z[6]cos(W)→Y↓Abs(W)-Q→Q↓Goto 4↓Lbl 4↓If W＜0：Then -Q→Q：-Y→Y：IfEnd↓X-Esin(Q)→X：Y+Ecos(Q)→Y↓Goto 5↓Lbl 5↓F：G：ZX- Z[3]→I：Y→J↓F+Icos(Z)-Jsin(Z)→X↓G+Isin(Z)+Jcos(Z)→Y↓Goto 6↓Lbl 6↓“X=” ：X◢“Y=” ：Y◢Pol(X-A, Y-A) ↓“I=” ：I◢If J＞0：Then“J=” ：J◢“J=” ：J+360◢IfEnd↓子程序：左线（QL）Lbl 0↓？A：？B↓Goto 1↓Lbl 1↓？D：？E↓If D>220249.506：Then Goto 2：IfEnd↓219215.613→C:1450→R:280→K:-29°47′23″→W↓54117.378→F:94561.972→G: 298°7′23″→Z↓Goto3↓Lbl 2↓If D>222040.807：Then “WU XIAO” ◢Goto 1：IfEnd↓221331.463→C:3000→R:0.00000000000001→K:-13°32′51″→W↓54060.248→F:92598.282→G: 268°20′1″→Z↓Goto3↓Lbl 3↓Prog“QXZBJS” ↓Goto1↓RETURN↓R：圆曲线半径K：缓和曲线长度W：转角（左偏取负值，右偏取正值）C：ZH点里程Z：第一切线方位角A：置镜点X坐标B：置镜点Y坐标F：交点X坐标G：交点点Y坐标D、E：待求点里程和偏距路线为分离式路线时，右线子程序与左线一样，根据右线的直曲转角表输入相应的参数即可导线测量:DXCLLbl 0↓Deg：Norm1↓（设置计算器角度单位为度及小数点显示为小数模式）?A：?B：“A1”→Z[8]:“B1”→Z[9] ↓Pol(A-Z[8], B-Z[9]) ↓If J>0：Then“ J”：J→Z[10]◢Else“J”: J+360→Z[10]◢IfEnd↓Goto 1↓Lbl 1↓?V:?L↓Z[10]+V-180→T↓“X”:A+Lcos(T) →X◢“Y”:B+Lsin(T) →Y◢Goto 1↓A、B：分别为置镜点X坐标和Y坐标“A1” ：“B1” 分别为后视点X坐标和Y坐标V：夹角L：边长缓和曲线的已知坐标反算里程和偏距程序主程序（程序名称：HHQXFS）Lbl 0↓Deg：Norm1↓（设置计算器角度单位为度及小数点显示为小数模式）Z[11] ：Z[12] ：Z[13] ：Z[14] ：M ：N ：R ：K：C ：Z↓（定义变量）Goto 1↓Lbl 1Pol(Z[13]-Z[11] ，Z1[4]-Z[12]) ↓I→Z[15] ↓If J>0：Then J→Z[16]：Else J+360→Z[16]：IfEnd↓Pol(X-Z[11] ，Y-Z[12]) ↓I→Z[17] ↓If J>0：Th en J→Z[18]：Else J+360→Z[18]：IfEnd↓Z[17]×cos(Abs(Z[16]-Z[18])) →Z[19]：Z[17]×sin(Abs(Z[16]-Z[18]))→Z[20]↓Z[19]×K÷Z[15]→SS- S5÷(40R2K2)→Z[21]：S3÷(6RK)- S7÷(336R3K3)→Z[22]↓Goto 2↓Lbl 2↓→Z[23]：tan-1 (Z[22] ÷Z[21])→Z[24]↓If N=1：Then Z+(-1)M×Z[24]→W：Else Z+180- (-1)M×Z[24]→W：IfEnd↓Z[11]+ Z[23] cosW→Z[25] ：Z[12]+ Z[23] sinW→Z[26]↓→Z[27]↓→Z[18]↓→Z[29]↓Goto 3↓Lbl 3↓If Z[18]> Z[16]：Then- Z[20]- Z[29]→E：Else Z[20]- Z[29]→E：IfEnd↓“D”：C-(-1) N×S+ →D◢（计算显示反算出的里程）其中为修正系数If M=1：Then“E”：E◢Else “E”：-E→E◢IfEnd↓（计算显示反算出的偏距）子程序：（如当我们想反算慕家梁隧道左线第二缓和曲线时）Lbl 0↓54102.075→Z[11] ↓（HZ点X坐标）94035.965→Z[12] ↓（HZ点Y坐标）54101.208→Z[13] ↓（HY点X坐标）94315.848→Z[14] ↓（HY点Y坐标）1→M：2→N：1450→R：280→K：220249.506→C：268°20′01″→Z↓?X：?Y （待反算点X坐标和Y坐标）Go to 1↓Lbl 1↓Prog“HHQXFS”（引用主程序“HHQXFS”）Goto 0↓Return↓（引用主程序后必须有的结束语句）变量含义：Z[11]、Z[12]、Z[13]、Z[14]：分别为缓和曲线的起终点X坐标和Y坐标；M ：缓和曲线左偏M=1，右偏M=2N ：计算第一缓和曲线时M=1，第二缓和曲线时M=2R ：圆曲线半径K：缓和曲线长C：缓和曲线的起点或终点里程（第一缓和曲线时为起点里程，反之为终点里程）；Z：第一切线方位角在实际工作中，当项目变换后，只需根据新项目中的曲线坐标表中给定参数改变子程序中的这些变量（Z[11] ：Z[12] ：Z[13] ：Z[14] ：M ：N ：R ：K：C ：Z）即可圆曲线反算程序：YQXFSLbl 0↓Deg：Norm1↓（设置计算器角度单位为度及小数点显示为小数模式）Z[30] ：Z[31] ：Z[32] ：Z[33] ：R ：L：K：Z：C↓（定义变量）Goto 1↓Lbl 1↓（180K÷(2 πR)）°→Z[34]↓Z+ ×Z[34]+ ×90→Z[35]↓Z[30]+R cos(Z[35])→Z[36]↓Z[31]+Rsin(Z[35])→Z[37]↓→Z[38]↓Pol(Z[30]- Z[36]，Z[31]- Z[37])↓J→Z[39]↓Pol(Z[32]- Z[36]，Z[33]- Z[37])↓J→Z[40]↓Pol(Z[X]- Z[36]，Y- Z[37])↓J→Z[41]↓Abs(Z[41]- Z[39])×(L-2K)÷Abs(Z[40]- Z[39]) →S↓“D”:C+S →D◢Z[38]-R→E↓If M=1：Then“E”：E◢：Else“E”：-E→E◢IfEnd↓子程序：（如当我们想反算慕家梁隧道左线圆曲线时）Lbl 0↓53993.248→Z[30]:94775.117→Z[31]↓54101.208→Z[32]↓94315.848→Z[33]↓1450→R:1033.893→L:298°7′23″→Z：219495.613→C↓?X：?Y↓Goto 1↓Lbl 1↓Prog“YQXFS”（引用主程序“YQXFS”）Goto 0↓Return↓（引用主程序后必须有的结束语句）Z[30]、Z[31]、Z[32]、Z[33]：分别为圆曲线的起终点X坐标和Y坐标；R：圆曲线半径L：曲线长K：缓和曲线长Z：第一切线方位角高程 GCLbl 0↓Deg：Norm1↓Z[41] ：Z[42] ：Z[43] ：C：R↓（定义变量）Z[41]-Z[42]→Z[44]↓R×Abs(Z[44])÷200→T↓Goto 1↓Lbl 1↓If D＜C-T：Then Goto2：IfEnd↓If D≤C：Then Goto3：IfEnd↓If D≤C+T：Then Goto4：IfEnd↓Lbl 2↓“H”：Z[43]-(C-T)×Z[41]÷100→H◢↓Goto 6↓Lbl 3↓(D+T-C) ÷2R→Y↓If Z[44]＜0：ThenY→Y：Else-Y→Y:IfEnd↓“H”：Z[43]-(C-T)×Z[41]÷100+Y→H◢↓Goto 6↓Lbl 4↓(C+T-D) ÷2R→Y↓If Z[44]＜0：ThenY→Y：Else-Y→Y:IfEnd↓“H”：Z[43]+(D-C)×Z[42]÷100+Y→H◢↓Got o 6↓Lbl 5↓“H”：Z[43]+(D-C)×Z[42]÷100→H◢↓Goto 6↓Lbl 6↓子程序：（左线）Lbl 0↓?D↓Goto 1↓Lbl 1↓If D＞221443.412：Then Goto2：IfEnd↓-2.2→Z[41]: -0.5→Z[42]:525.222→Z[43]:219640→C:45000→R↓Goto 3↓Lbl 2↓If D＞222000：Then “WU XIAO” ◢Goto 0：IfEnd↓-0.5→Z[41]: -2.473→Z[42]:514.972→Z[43]:221690→C:25000→R↓Goto 3↓Lbl 3↓Prog“GC”（引用主程序“HHQXFS”）Goto 0↓Return↓（引用主程序后必须有的结束语句）Z[41]、Z[42]：W 、WZ[43] ：交点高程C：交点里程桩号R：数曲线半径直线反算程序：（ZXFS）“A”?→Z[45]:“B”?→Z[46]↓:？C:？ZLbl 1↓？X：？Y↓Pol(X-Z[45], Y-Z[46])J-Z→Z[47]↓C+Icos(Z[47])→D↓“D”:D◢“E”:Isin(Z[47])→E◢Goto 1↓放样程序：（FY）Deg：Norm1↓（设置计算器角度单位为度及小数点显示为小数模式）?A：?B↓Goto 1↓Lbl 1↓?X：?Y↓Pol(X-A, Y-B)↓“I”:I◢If J>0：Then“ J”：J◢Else“J”: J+360◢IfEnd↓Goto 1↓添加额外变量：例如想添加100个额外变量100→FUNCTION→｛PROG｝→SHIFT→(Dimz) →EXE当显示器上显示“Done”时，表示添加了指定数量的额外变量。
非对称缓和曲线坐标计算程序
非对称缓和曲线坐标计算程序CASIO fx-4800P QXZB曲线坐标计算CASIO4800 QXZB可计算不等缓和曲线、圆曲线上的任意中、边桩坐标：该程序适用于计算器 CASIO fx-4800P，可计算与线路中心成任意夹角的缓和曲线、圆曲线中、边桩坐标及待测点方位角和距离。
1、DK(JD)？输入交点桩号2、X(JD)？输入交点坐标X3、Y(JD)？输入交点坐标Y4、T1？输入第一切线长(如果只有一条切线两者都输入一致)5、T2？输入第二切线长(如果只有一条切线两者都输入一致)6、FWJ？输入直线方位角(ZH→JD)7、A？输入转角：左转为负，右转为正8、R？输入圆曲线半径9、LS1？输入第一缓和曲线长(如果只有一条缓和曲线两者都输入一致)10、LY？输入圆曲线长(L-LS1-LS2)11、LS2？输入第二缓和曲线长(如果只有一条缓和曲线两者都输入一致)12、X(ZJD)？输入置镜点坐标X13、Y(ZJD)？输入置镜点坐标Y14、JSDK？输入前视点里程15、PL？输入偏距16、PA？输入偏角程序下载地址：/blog/post/QXZB-4800.html评价答案好:18不好:1原创:18非原创:0菲メ帆ぅ回答采纳率:52.8% 2010-06-02 17:37满意答案好评率：57%（for Casio-fx4850）扩展变量操作（15个）：Defm 15←┚ ( O为字母、0为数字)J-PQX (平面数据输入，自行切换到J-JSMS)Defm 15←┚A“JD” B“JDX” C“JDY” F“FWJ” O“A0：Z-，Y+” RE“LS1” K“LS2”：E<1=>E=1E-9⊿K<1=>K=1E-9⊿Z[1]=EE÷24R-E∧4÷2688RRR:Z[2]=E÷2-EEE÷240RR:X=(EE-KK)÷24R÷sin Abs O :“T1=”:Z[3]=(R+Z[1])tan(Abs O÷2)+Z[2]-X◢“T2=”:Z[4]=(R+KK÷24R-K∧4÷2688RRR)tan(Abs O ÷2)+K÷2-KKK÷240RR+X◢“L=”:L=Abs OπR÷180+(E+K)÷2◢J=tan-1((R+Z[1])÷(Z[3]-Z[2]):“E=”:X=(R+Z[1])÷sin J-R◢X=A-Z[3]:Y=X+E:E<1=>“ZY=”:X◢≠=> “ZH=”:X◢“HY=”:Y◢⊿ “QZ=”:Y =X+(L-K-E)÷2+E◢Y=X+L-K:X=X+L:K<1=> “YZ=”:X◢≠=> “YH=”:Y◢“HZ=”:X◢⊿Prog“J-JSMS” ←┚J-JSMS（放样模式主程序）“1-ZS, FY,3-FS”: Lb1 0:{Z}: Z“MS”≤2=>Goto 1: ≠=> Z≤4=> Goto2⊿⊿←┚Lb1 1 :Q“X-YIQI”S“Y-YIQI”T“X-HOUSI”U“Y-HOUSI”:{PDW}:P“F=”D“BZ=”W“BJ=”: Prog“JP”: Prog“FY”: Goto 1←┚Lb1 2 :{XYW}:XYW“BJ”: Prog“JF”: “P=”: P◢“BZ=”: D◢JP（平曲线正算子程序）Fixm←┚Lb1 1 :J=F:X= B-Z[3]cos F :Y= C-Z[3]sin F :G=E←┚P≤A-Z[3]=>I=A-Z[3]-P:M=-I:N=0:H=F+W:Goto 5←┚≠=>P≤A-Z[3]+E=>I=P-A+Z[3]:H=90II÷REπ:O<0=>H=-H⊿H=H+W+F:Goto 3←┚≠=>P≤A-Z[3]+L-K=>I=P-A+Z[3]-E:Goto 4: ≠=>Goto 2⊿⊿←┚Lb1 2 :X= B+Z[4]cos(F+O :Y= C+Z[4]sin(F+O :J=F+O+180:G=K←┚P≤A-Z[3]+L=>I=A-Z[3]+L-P:H=90II÷RKπ:O>0=>H=-H⊿H=H+J+W+180:Goto 3 ←┚≠=>I=P-A+Z[3]-L:M=-I:N=0:H=J+W+180:Goto 5 ←┚Lb1 3 :M=I-I∧5÷40RRGG:N=III÷6GR-I∧7÷336RRRGGG:Goto 5 ←┚Lb1 4 :H=(E+2I)×90÷πR:M=R sin H+Z[2]:N=R(1-cos H)+Z[1]:O<0=>H=-H⊿H=J+H+W: Goto 5←┚Lb1 5 :P≤A-Z[3]+L-K=> O<0=>N=-N⊿≠=>O>0=>N=-N⊿⊿Goto 6←┚Lb1 6 :X =X+M cos J-N sin J+D cos H:Y =Y+M sin J+N cos J+D sin H“X=”: X ◢“Y=”: Y ◢JF（平曲线反算子程序）Fixm←┚U=X:V=Y:D=0:J=F-W:P=A+(Y-C)cos J-(X-B)sin J←┚Lb1 1:Prog“JP”:J=H-180:I=(V-Y)cos J-(U-X)sin J:Abs I<1E-4=>Goto 2:≠=>P=P+I:Goto 1⊿←┚Lb2: D=(V-Y)÷sinHFY(放样计算子程序)Pol(T-Q,U-S):“HOUSI D=”:I◢J<0=>J=J+360⊿J>360=>J=J-360⊿“HOUSI J=”:J→DMS◢Pol(X-Q,Y-S):“I=”:I◢J<0=>J=J+360⊿J>360=>J=J-360⊿“J=”:J→DMS◢一、程序简介本套程序是共有2个主程序，3个子程序。
坐标计算实例(缓和曲线)
缓和曲线逐桩坐标计算（转载）摘要：利用一缓和曲线算例，通过数学分析，推导出缓和曲线逐桩坐标计算公式，此公式可作为道路测设中的范例来运用，有很强的指导意义。
关键词：缓和曲线、公式、逐桩坐标一、引言道路建设中，由于受地形或地质影响，经常需要改变线路方向，为满足行车要求，往往要用曲线把两条直线连接起来。
曲线的构成形式无外乎圆曲线和缓和曲线，本文以河北省沿海高速某曲线段为例推导出缓和曲线的逐桩坐标计算公式，以方便图纸的审核，满足施工放样的需求。
本公式具有良好的操作性，方便施工、提高精度，可作为道路测设中的范例运用。
二、公式推导1 、实例数据河北省沿海高速公路一缓和曲线（如图）：AB 段为缓和曲线段，A 为ZH 点，B 为HY 点，R B=800m ；A 点里程为NK0+080 ，切线方位角为θA=100 ° 00 ′ 24.1 ″，坐标为X A=4355189.493,Y A=476976.267 ；B 点里程为NK0+158.125 ，切线方位角为θB=102 ° 48 ′ 15.6 ″，坐标为X B=4355174.669 ，Y B=477052.964 ，推求此曲线段内任意点坐标。
2 、公式推导及实例计算方法一：弦线偏角法1 ）公式推导由坐标增量的计算方法我们不难理解，求一点坐标可以根据其所在直线的方位角以及直线上另一点的坐标和距待求点的距离。
所以我们可以利用ZH 点，只要知道待求点距ZH 点的距离（弦长S ）和此弦与ZH 点切线方位角的夹角（转角a ），即可求出该点坐标。
根据回旋线方程C=RL ，用B 点数据推导出回旋线参数：C=RL S=800*78.125=62500 （L S为B 点至ZH 点的距离）设待求点距ZH 点距离为L因回旋线上任意点的偏角β0=L2/2RL S, 且转角a=β0/3 ，可得该点转角a 。
（曲线左转时a 代负值）。
根据缓和曲线上的弧弦关系S=L-L5/90R2L S2，可以求出待求点至ZH 点的弦长。
casio5800缓和曲线、圆曲线、直线程序
LbI 0:4→Dimz:Cls"1.SZ→XY""2.XY→SZ"?→N:Cls"XO"?U: "YO"?V: "SO"?O:"FO"?G: "LS"?H: "RO"?P:"RN"?R:?Q1÷P→C：(P-R)÷(2HPR) →D：180÷π→EN=2=>Goto 2LbI 1:Cls:?SS=0=>Goto 0 (注:当S输入0时,程序会由LbI 0开始重新运行) ?Z:Abs(S-O)→WProg"SUB1":Cls"XS=":X◢"YS=":Y◢"FS=":F-90▼DMS◢Goto 1LbI 2:Cls:?XX=0=>Goto 0 (注:当X输入0时,程序会由LbI 0开始重新运行) ?Y:X→I:Y→JProg"SUB2":O+W→S:Cls"S=":S◢"Z=":Z◢Goto 2SUB1(子程序名)0.1739274226→A:0.3260725774→B0.0694318442→K:0.330009482→L1-L→F:1-K→MG+QEKW(C+KWD)→Z[1]G+QELW(C+LWD)→Z[2]G+QEFW(C+FWD)→Z[3]G+QEMW(C+MWD)→Z[4]Acos(Z[1])→XX+Bcos(Z[2])→XX+Bcos(Z[3])→XX+Acos(Z[4])→XU+WX→XAsin(Z[1])→YY+Bsin(Z[2])→YY+Bsin(Z[3])→YY+Asin(Z[4])→YV+WY→YG+QEW(C+WD)+90→FX+Zcos(F)→XY+Zsin(F)→YSUB2(子程序名)G-90→T(Y-V)cos(T)-(X-U)sin(T)→WAbs(W)→W:0→ZLbI 4:Prog"SUB1"T+QEW(C+WD)→L(J-Y)cos(L)-(I-X)sin(L)→ZIf Abs(Z)＜0.000001:Then 0→Z :Prog"SUB1"(J-Y)÷sin(F)→ZElse W+Z→W:Goto 4:IfEnd子程序（平面线形数据库）PM-SJ-1÷P→ CG-90→TLbl 1ifS ≥500（线元起点里程）AndS＜769.256（线元止点里程）：Then19942.837 →U（线元起点X坐标）：28343.561 →V（线元起点Y坐标）：500 →O（线元起点里程）：125°16°31°→G（线元起点方位角）：269.265 →H（线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥769.256（线元起点里程）AndS＜806.748（线元止点里程）：Then19787.34 →U（线元起点X坐标）：28563.378→V（线元起点Y坐标）：769.256→O（线元起点里程）：125°16°31°→G （线元起点方位角）：37.492→H（线元长度）：10^（45）→P（线元起点曲率半径）：221.75→R （线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥806.748（线元起点里程）AndS＜919.527（线元止点里程）：Then19766.566→U（线元起点X坐标）：28594.574→V（线元起点Y坐标）：806.748→O（线元起点里程）：120°25°54.07°→G（线元起点方位角）：112.779→H（线元长度）：221.75→P（线元起点曲率半径）：221.75→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥919.527（线元起点里程）AndS＜999.812（线元止点里程）：Then19736.072 →U（线元起点X坐标）：28701.893 →V（线元起点Y坐标）：919.527→O（线元起点里程）：91°17°30.63°→G（线元起点方位角）：80.285→H（线元长度）：221.75→P（线元起点曲率半径）：9579.228→R（线元终点曲率半径）：-1→Q（线元左右偏标志：左负右正）：IfEnd←┘ifS ≥999.812（线元起点里程）：Then19744.038 →U（线元起点X坐标）：28781.659→V（线元起点Y坐标）：999.812→O（线元起点里程）：80°40°50°→G（线元起点方位角）：100→H （线元长度）：10^（45）→P（线元起点曲率半径）：10^（45）→R（线元终点曲率半径）：0 →Q（线元左右偏标志：左负右正）：IfEnd←┘使用说明1、规定(1) 以道路中线的前进方向（即里程增大的方向）区分左右；当线元往左偏时，Q=-1；当线元往右偏时，Q=1；当线元为直线时，Q=0。
在excel vba 中计算缓和曲线的方法
在excel vba 中计算缓和曲线的方法在Excel VBA中计算缓和曲线的方法可以通过使用内置的函数和公式来实现。
缓和曲线是一种平滑的曲线，通常用于在两个点之间连接线条。
它可以用于绘制折线图、曲线图以及其他图表类型。
以下是在Excel VBA中计算缓和曲线的方法：1. 使用平滑函数：Excel VBA中有几种平滑函数可以用来计算缓和曲线，如LOESS函数和Spline函数。
这些函数会基于数据点生成平滑曲线，然后可以将其用于绘制图表。
2. 使用线性插值函数：线性插值可以使用Excel VBA中的LINEST函数来完成。
可以将已知点的坐标作为输入，并使用线性插值来计算缓和曲线的每个点的坐标。
3. 使用曲线拟合函数：曲线拟合是一种根据已知数据点来生成缓和曲线的方法。
Excel VBA中有多种曲线拟合工具，如多项式拟合、指数拟合和幂函数拟合等。
可以选择适合数据的拟合类型，并使用拟合函数来计算缓和曲线。
除了以上方法，还可以通过手动计算缓和曲线来实现。
该方法需要一定的数学知识和计算能力，但可以通过编写VBA代码来自动执行计算。
手动计算缓和曲线的方法包括：4. 使用Bezier曲线：Bezier曲线是一种通过控制点来生成平滑曲线的方法。
可以使用Excel VBA中的Bezier函数来计算缓和曲线。
5. 使用B样条曲线：B样条曲线是一种基于控制点和节点来生成平滑曲线的方法。
可以使用Excel VBA中的B样条曲线函数来计算缓和曲线。
以上是在Excel VBA中计算缓和曲线的方法。
可以根据需要选择适合的方法，然后编写代码来自动计算并绘制缓和曲线。
(完整word版)缓和曲线计算程序说明
缓和曲线测设数据计算程序说明（一）缓和曲线测设数据计算程序设计框图根据《测量学》(第四版）§10-4中所介绍的计算公式，设计“缓和曲线测设数据计算程序”（Ecurve。
lsp）。
其步骤如下：(1）设计“缓和曲线测设数据计算程序”的数据输入对话框（文件名为:Ecurve。
dcl）,用以输入道路中线的已知数据（交点桩号和坐标、前后转点坐标），设计数据（圆曲线半径、缓和曲线长度、曲线分段长度、道路半宽）等.（2) 设计“从编辑框取值”子程序，将字符型数据变换为实数型。
（3）应用“方位角函数"按道路中线点坐标计算中线方位角和交点处的道路偏角;根据交点桩号、偏角和曲线的设计数据计算曲线元素（切线长、曲线长、外矢距）,曲线主点桩号和坐标。
（4) 在独立坐标系中计算缓和曲线细部点及其边线点、圆曲线细部点及其边线点的独立坐标,然后变换为大地坐标；在CAD屏幕绘制曲线的主点和细部点.(5) 文件输出曲线测设所需要的主点和细部点的坐标数据。
图1为缓和曲线测设数据计算程序设计框图.图1 缓和曲线测设数据计算程序框图（二）缓和曲线坐标计算程序中变量的标识符规定1。
设计数据S2 —注记字高(单位：毫米）;PNjd —交点桩号（以后PN代表其余各桩桩号的开头）；JD —交点点位；Xj, Yj - 交点坐标；FD， BD —前视转点，后视转点的点位（F-fore,B—back）;Xf， Yf， Xb，Yb —前视转点,后视转点的坐标;R —圆曲线设计半径;Ls —缓和曲线设计长度；LN, N —缓和曲线分段间距，分段数；CN —圆曲线整桩间距;D —道路半宽；W —道路设计宽度。
2。
道路转向偏角和曲线计算元素Ab，bA,Af，fA —交点至前视转点、后视转点的正、反方位角（A-azimuth）;Am —内矢距方向的方位角;LR —判断路线左偏或右偏（1—左偏，2—右偏）；A —路线偏角 ；m， p —增设缓和曲线后的曲线增长，曲线内移距；B0 - 缓圆点切线方位角βo；ZH，HY，QZ,YH,HZ —曲线主点：直缓、缓圆、曲中、圆缓、缓直点的点位；ZH1，HY1,QZ1，YH1，HZ1 —直缓、缓圆、曲中、圆缓、缓直点的注字点位；Aq, Ah, Az - QZ、HY、YH点的注字方向的方位角；Xe，Ye — YH或HY点的独立坐标（在转换为大地坐标前）；L, Lc - 曲线总长度，圆曲线长度；TT， EE —切线长,外矢距;HY0，YH0 — HY点、YH点在切线上的投影点位。
应用Excel编制缓和曲线坐标计算程序
应用Excel编制缓和曲线坐标计算程序一、引言随着道路事业的飞速发展，各种先进的测量仪器和现代化的测设手段已被广泛采用。
曲线测设是道路建设放样的重要组成部分之一。
现在大多数道路工程曲线较多，加密桩放样随之增加，这就需要批量计算数据，但这往往给工程的数据处理带来许多麻烦。
兼于这些特点，曲线坐标计算非常适合在Excel 中进行。
应用Excel编制缓和曲线坐标计算程序和目前常用的程序相比，应用Excel进行缓和曲线坐标计算既不需要专门的程序和软件，也不需要用到计算机程序语言，直接应用曲线计算公式和Excel的功能就可以完成计算，速度较快，灵活性较强。
二、计算过程应用Excel编制缓和曲线坐标计算的过程实际上就是输入数据和输入公式的过程，关键是充分利用Excel所具有的函数功能，使计算变的更为简便而又直观。
1、已知数据输入首先，根据（表1）的提示将已知数据第一缓和曲线长、第二缓和曲线长、半径、偏角、初始切线方位角、直缓点里程、转向（左偏为负，右偏为正）、直缓点坐标、置镜点坐标、后视点坐标以及需要计算的里程一一输入相对应的单元格。
2、曲线元素计算如表1所示（1）、计算第一切线增长值m1在A7中输入“=A4/2-A4^3/240/C4^2”（2）、计算第二切线增长值m1在B7中输入“=B4/2-B4^3/240/C4^2”（3）、计算内移量p1、p2在C7中输入“=A4^2/24/C4-A4^4/2688/C4^3”在D7中输入“=B4^2/24/C4-B4^4/2688/C4^3”表1（4）计算切线长T1、T2在E7中输入“=A7+(C4+C7)*TAN(RADIANS((INT(D4)+INT((D4-INT(D4))*100)/60+(D4*100-INT(D4*100))*100/3600))/2)+(D7-C7)/SIN(RADIANS((INT(D4)+INT((D4-INT(D4))*100)/60+(D4 *100-INT(D4*100))*100/3600)))”在F7中输“=B7+(C4+D7)*TAN(RADIANS((INT(D4)+INT((D4-INT(D4))*100)/60+(D4*100-INT(D4*100))*100/3600))/2)+(C7-D7)/SIN(RADIANS((INT(D4)+INT((D4-INT(D4))*100)/60+(D4*100-IN T(D4*100))*100/3600)))”（5）计算曲线长L。