由本站论坛坛友　lsl47477948　提供，谢谢！！


以下是我自己编的测量程序,本人不是搞测量专业的,内有不足之处望大家多多指正,缓和曲线和曲线公式是我自己看着书上编的,不知对不对,希望大家支持! 注意π是圆周率
*********ZB-ZS**********                    [文件名:坐标-正算]
Lb1 1
"-----ZB--ZS-----"                                 ［标题，如文件名］ {A,B,Q,L}                                                      ［变量］
A"ZJD-X"B"ZJD-Y":X=A+L"JU-LI"cosQ"FWJ" ◢   [ZJD-X/Y置镜点X/Y轴/JU-LI距离(长度)/FWJ方位角]     Y=B+LsinQ◢ O=2                                        ［为下步选择时默认值］
{O}                                         [选择是否边线坐标计算O为字母]
"---BIAN-XIAN---"：O"1.YES   2.NO"：O≠1=>Goto 1:△  
P=Q:Prog "BX":Goto 1                                    [BX边线]
**********ZB-FS**********                   [文件名:坐标-反算]
LB1 1
"-----ZB--FS-----"                         ［标题，如文件名］
{A,B,C,D}                                   [ZJD-X/Y置镜点X/Y轴/FXD-X/Y方向点X/Y轴L(m)两点距离(长度)/∫开平方]A"ZJD-X"B"ZJD-Y"C"FXD-X"D"FXD-Y":X=C-A:Y=D-B:Prog "FWJ":L"L(m)"=∫(X^2+Y^2)◢
Q"FWJ="◢
goto 1
**********ZB-FY**********                    [文件名:坐标-放样]
LB1 1
"-----ZB--FY-----"
{A,B,C,D}                                [ZJD-X/Y置镜点X/Y轴/FXD-X/Y方向点X/Y轴]
A"ZJD-X"B"ZJD-Y"C"FXD-X"D"FXD-Y":X=C-A:Y=D-B  
prog "FWJ":P=Q {E,F} X=E"DCD-X"-A:Y=F"DCD-Y"-B                    [DCD-X/Y待测点X/Y轴数据]  prog "FWJ"
L"L(m)"=∫(X^2+Y^2)◢                        [待测点-置镜点距离/∫开平方]
N"PJ"=Q-P◢                                  [置镜点-后视方向点=>待测点偏角]
GOTO 1                                       [正数顺时针,负数逆时针]
**********QX-HQ**********                            [文件名:曲线-缓曲]
LB1 1
" HUAN-HE-QU-XIAN"                     ［标题：缓和曲线，前面加一空格］
{A,B,C,D,R,M}
A"ZH-X"B"ZH-Y"C"JD-X"D"JD-Y"R"QX-R"M"HQ-LS"  X=C-A:Y=D-B                                [QX-R曲线半径/HQ-LS缓和曲线长度] prog "FWJ" LB1 2{G}                                            [选择曲线方向Y向右/Z向左]
G"QFX 1.You 2.Zuo":G=1=>goto 4:△
G=2=>goto 3:△
goto 2
LB1 3
G=-1
LB1 4
{L,O}
J=L"HQ-L"^2÷(6RM)×(180/3.14159)G                [HQ-L直缓点到缓和曲线上任一点长度]
L>M=>goto 4:△                               [输入任一点长度大于缓和曲线长时返回］
H=L-(L^5÷(40R^2M^2))
I=H÷COS(Abs(J))
E"HQ-X"=A+Icos(Q+J)◢       [HQ-X/Y缓和曲线上任一点X/Y坐标]
F"HQ-Y"=B+Isin(Q+J)◢
"----XUAN-ZHE----":O"1.=>BX  2.=>YQX"     [BX选择是否进入缓和曲线上任一点边线计算]
O=1=>goto 5:△                               [YQX选择是否进入圆曲线计算]［O字母］
O=2=>goto 6:△ goto 4
LB1 5
P=Q+3J:X=E:Y=F Prog "BX"
goto 4
LB1 6                                    [进入圆曲线]
L<M=>goto 4:△                     [必须使缓和曲线上任一点L长度=缓和曲线长度M,否则返回]
LB1 7
"--YUAN-QU-XIAN--"
{L}
K=(360L"YQ-L")÷(2Rπ)G                      [YQ-L缓圆点-圆曲线上任一点长度]
L>(2R3.14159)=>goto 7:△  
I=∫(2R^2-2R^2cos(Abs(K)))                   [∫开方符号Abs取绝对值符号]
X"YQ-X"=E+Icos(Q+3J)+(K/2))◢                [YQ-X/Y圆曲线上任一点坐标值]
Y"YQ-Y"=F+Isin(Q+3J)+(K/2))◢
O=2                                          [O字母]
{O}                                          [BX选择是否进入圆曲线上任一点边线计算]
"---BIAN--XIAN---":O"1.YES  2.NO":O=1=>goto 8:△
goto 7
LB1 8
P=Q+3J+K
prog "BX"
goto7
**********QX-YQ**********                    [文件名:曲线-圆曲］
LB1 1
"--YUAN-QU-XIAN--"
{A,B,C,D,R}
A"HY-X"B"HY-Y"C"JD-X"D"JD-Y"R"YQ-R"
X=C-A:Y=D-B
prog "FWJ"
LB1 2
{G}
G"QXFX=> 1.Y  2.Z":G=1=>goto 4:△
G=2=>goto 3:△ goto 2
LB1 3
G=-1
LB1 4
{L}
K=(360L"YQ-L")÷(2Rπ)G                   ［YQ－L圆曲线上任一点到ZY点距离]
L>(2Rπ)=>goto 4:△
I=∫(2R^2-2R^2cos(Abs(K)))
X"YQ-X"=A+Icos(Q+K/2)◢
Y"YQ-Y"=B+Isin(Q+K/2)◢
O=2
{O}
"---BIAN--XIAN---":O"1.YES  2.NO":O=1=>goto 5:△
goto 4
LB1 5
P=Q+K
prog "BX"
goto 4
以下是上面文件运行时需要调用的文件[文件名要改时必须和上面所用到的一致]**********BX**********                   [文件名：边线]
"--BIAN-XIAN-XY--"
N=90
{N,Z}
S"ZUO-X"=X+Z"ZUO-L"cos(P-N"PJ")◢          [ZOU-X/Y左边线X/Y坐标值/ZUO-L左边线长度]
T"ZUO-Y"=Y+Zsin(P-N)◢                    [PJ计算边线时所用偏角,默认下90度]
U=Z
{U}
V"YOU-X"=X+U"YOU-L"cos(P+N"PJ")◢    [YOU-X/Y右边线X/Y坐标值/YOU-L右边线长] W"YOU-Y"=Y+Usin(P+N)◢
**********FWJ**********                  [文件名：方位角]
X≠0=>goto 2:△X=0=>goto 4:△
LB1 2
Q=tan-1(Y/X) X>0=>goto 3:△X<0=>Q=Q+180:goto 5:△
LB1 3
Y>0=>goto 5:△Y≤0=>Q=Q+360:goto5:△
LB1 4
Y=0=>Q=0:goto 5:△Y>0=>Q=90:goto 5:△Y<0=>Q=270:goto 5:△
LB1 5                                    [最好不要去掉这个空目录名]