公路通用復化辛普森公式匝道點位坐標計算4800源程序
------------------杭浦高速臨平互通---------------------
本文利用的是計算公路匝道點位坐標的復化辛普森通用公式數學模型,集直線、圓曲線、回旋線通用，占字符內存較小，計算精度不限的程序
一、運行變量名稱說明：
V=1、2分別進入坐標計算、樁號反算
K1、K2-------曲線起點、終點里程
F0-----------曲線起點方位角
R1、 R2------曲線起點、終點半徑 (ρ左-右+，0為直線)
X0、Y0-------曲線起點、終點坐標
M------------求和累積次數n的2倍(偶數)，精度迭代次數
K------------曲線待求點里程
BP-----------求點左右偏距(左-右+)
ANG-------- -求點的右斜交角
X、Y---------曲線求得坐標
FW-----------待求點的即時切線方位角
XF、YF-------為需求樁號的點坐標
DL、K+O、LP分別為樁號誤差、求得樁號、左右偏距(左-右+)
當曲線的設計半徑較小時，為保證點位計算精度，M(即程序中n的2倍)的取值可適當的大些。M為偶數，直線時M=2即可，經計算M=16即可滿足半徑為60的小半徑曲線精度。
二、曲線計算程序名: Prog "CURVE"
Defm 4
V"V=1 2"
Lbl 0:{KLW}
Lbl 4
Q"OPT:M0AB1C2D3E4FH5G6I7J8CR9"=0=>Prog "M"△
Q=1=>Prog "AB"△
Q=2=>Prog "C"△
Q=3=>Prog "D"△
Q=4=>Prog "E"△
Q=5=>Prog "FH"△
Q=6=>Prog "G"△
Q=7=>Prog "I"△
Q=8=>Prog "J"△
Q=9=>Prog "CR"△
A"K1"B"K2"C"F0"D"R1"E"R2"F"X0"G"Y0"
D≠0=>I=1/D:≠=>I=D△E≠0=>J=1/E:≠=>J=E△
AbsD+AbsE=0=>M=2:≠=>M=16△
V=2=>L=0:W=90△KL"BP"W"ANG"
N=0：Z[1]=0：Z[2]=0：Z[3]=0：Z[4]=0 進入坐標迭代計算
Lbl2
N=N+1：H=2(K-A)/M：R=NH/2+A：R=C+180/π*(I+(J-I)/2(B-A)*(R-A))*(R-A)
Int(N/2)=N/2=>Z[1]=Z[1]+cosR：Z[2]=Z[2]+sinR：≠=>Z[3]=Z[3]+cosR：Z[4]=Z[4]+ sinR△
N=M=>Goto3：≠=>Goto 2
Lbl3
X=F+H/6*(cosC+4Z[3]+2Z[1]-cosR)+Lcos(R+W)
Y=G+H/6*(sinC+4Z[4]+2Z[2]-sinR)+Lsin(R+W)
V=2=>Goto 6△
X"X="◢
Y"Y="◢
R"FW"=R-360Intg(R/360◢
Goto 0
Lbl 6 進入樁號求算
Pol(T"XF"-X,U"YF"-Y
O=Icos(J-R
AbsO≤1e-4=>O"DL"◢K=K+O◢
O"LP"=Isin(J-R◢
{TU}
Goto 6:≠=>K=K+O:L=0:Goto 4
三、數據文件：
線元要素數據文件每行為一個線元段，逐句執行賦值，直至不滿足、運行完成。
主線"M":
K>=5865=>A=5865:B=586.526:C=101`49`58`:D=0:E=0:F=3361488.152:G=449202.642
K>B=>A=B:B=6235.526:C=101`49`58`:D=0:E=-3250:F=3361488.044:G=449203.157
K>B=>A=B:B=7674.84:C=98`34`16.9`:D=-3250:E=-3250:F=3361419.068:G=449566.617
K>B=>A=B:B=8044.84:C=73`11`49.3`:D=-3250:E=0:F=3361521.529:G=450990.516
K>B=>A=B:B=8273.243:C=69`56`7.9`:D=0:E=0:F=3361641.834:G=451340.355
"8+273=8+100 S=1 2":{S}
S=2=>A=8100:B=9610:C=C:D=0:E=0:F=3361720.194:G=451554.895
注：主線中S=1、2為斷鏈1、2段(如K8+273.243=K8+100)
匝"AB"
K≥0:A= 0:B=716.396:C=85`46`15`:D=800:E=800:F=3361721.181:G=450482.643
K>B=>A=B:B=776.396:C=137`4`44`:D=0:E=-60:F=3361468.151:G=451127.474
K>B=>A=B:B=842.929:C=108`25`51.6`:D=E:E=-60:F=3361431.988:G=451174.516
K>B=>A=B:B=902.929:C=44`53`50`:D=E:E=0:F=3361446.56:G=451235.989
K>B=>A=B:B=947.929:C=16`14`58`:D=E:E=180:F=3361499.992:G=451261.793
K>B=>A=B:B=1047.145:C=23`24`41.1`:D=E:E=180:F=3361542.602:G=451276.164
K>B=>A=B:B=1141.034:C=54`59`34.26`:D=E:E=0:F=3361618.517:G=451338.083
匝"C"
K≥842.929:A= 842.929:B=932.843:C=44`53`50`:D=-60:E=-60:F=3361446.56:G=451235.989
K>B=>A=B:B=1067.843:C=319`02`08.44`:D=E:E=0:F=3361528.246:G=451238.793
K>B=>A=B:B=1204.663:C=254`34`40`:D=E:E=0:F=3361541.195:G=451111.894
匝"D"
K≥0:A= 0:B=46.596:C=115`7`37`:D=791.747:E=60:F=3361640.231:G=450877.857
K>B=>A=B:B=101.697:C=139`3`39.2`:D=E:E=60:F=3361614.603:G=450916.347
K>B=>A=B:B=234.195:C=191`40`42.5`:D=E:E=3237.41:F=3361563.143:G=450929.781
匝"E"
K≥0:A= 0:B=53.571:C=115`30`28`:D=0:E=-420:F=3361652.781:G=450889.678
K>B=>A=B:B=303.928:C=111`51`13.6`:D=E:E=-420:F=3361630.748:G=450938.498
K>B=>A=B:B=385.416:C=77`42`1.6`:D=E:E=0:F=3361610.206:G=451184.308
K>B=>A=B:B=516.85:C=72`8`31.93`:D=E:E=0:F=3361632.665:G=451262.605
匝"FH"
K≥0:A= 0:B=687.778:C=265`46`15`:D=0:E=0:F=3361721.183:G=450482.658
K>B=>A=B:B=747.778:C=265`46`15`:D=E:E=-60:F=3361670.463:G=449796.752
K>B=>A=B:B=956.691:C=237`7`22.6`:D=E:E=-60:F=3361656.351:G=449739.118
K>B=>A=B:B=1038.357:C=37`37`32.61`:D=E:E=0:F=3361569.331:G=449819.209
K>B=>A=B:B=1242.528:C=358`37`59`:D=E:E=0:F=3361647.7:G=449835.267
匝"G"
K≥838.838:A= 838.838:B=898.838:C=150`10`2`:D=-60:E=0:F=3361576.113:G=449719.64
K>B=>A=B:B=949.463:C=121`31`9.6`:D=E:E=160:F=3361537.147:G=449764.388
K>B=>A=B:B=1033.468:C=130`35`1.34`:D=E:E=160:F=3361508.476:G=449806.043
K>B=>A=B:B=1121.958:C=160`39`56.8`:D=E:E=986.28538:F=3361439.935:G=449852.931
匝"I"
K≥0:A= 0:B=66.176:C=265`46`15`:D=0:E=85:F=3361694.914:G=450016.224
K>B=>A=B:B=105.294:C=288`4`27`:D=E:E=85:F=3361698.579:G=449950.594
K>B=>A=B:B=157.408:C=314`26`32.5`:D=E:E=400:F=3361718.698:G=449917.45
K>B=>A=B:B=256.758:C=335`44`20`:D=E:E=400:F=3361762.343:G=449889.539
K>B=>A=B:B=337.758:C=349`58`11`:D=E:E=0:F=3361857.034:G=449860.326
匝"J"
K≥0:A= 0:B=66.176:C=265`46`15`:D=0:E=-85:F=3361678.559:G=450017.433
K>B=>A=B:B=104.77:C=243`28`2.8`:D=E:E=-85:F=3361665.281:G=449953.056
K>B=>A=B:B=161.779:C=217`27`9`:D=E:E=-350:F=3361640.922:G=449923.549
K>B=>A=B:B=215.557:C=193`34`20.8`:D=E:E=-350:F=3361588.865:G=449901.365
K>B=>A=B:B=279.843:C=184`46`8`:D=E:E=0:F=3361535.825:G=449892.803
K>B=>A=B:B=387.002:C=179`30`25`:D=E:E=0:F=3361471.579:G=449891.389
匝"CR"
K≥3080:A= 3080:B= 3666.041:C=175`46`15`:D=0:E=0:F=3362052.282:G=449829.361
K>B=>A=B:B= 3786.041:C=175`46`15`:D=E:E=900:F=3361467.83675:G=449872.57911
K>B=>A=B:B=3920:C=179`35`25.95`:D=E:E=900:F=3361348.02009:G=449878.76613
其他說明：數字前的乘號不能省。如A5、π6不當，應改為A*5、π*6或5A、6π。
括號前的乘號可省,三角函數前的乘號可省。如H(A+B)、LcosR正確
以下形式為連續除法：1/AB=1/A/B,1/5(B+C)=1/5/(B+C),1/A^2=1/A/A。
句末尾的括號可省。如Pol(T"XF"-X,U"YF"-Y。
=>前不用冒號
≠=>前有冒號
△前后不用冒號，且可不換行。
李甲生 lijs0483@163.com
